snd::analysis_tools::DSPlane Class Reference

#include <sndDSPlane.h>

Collaboration diagram for snd::analysis_tools::DSPlane:

Classes
struct	DSHit
struct	rl_pair

Public Member Functions
	DSPlane (std::vector< MuFilterHit * > snd_hits, const Configuration &configuration, MuFilter *muon_filter_geometry, int station)
const int	GetNHits () const
const int	GetStation () const
const double	GetTotQdc () const
const std::vector< DSHit >	GetHits () const
const int	GetBarNHits (int bar_to_compute) const
void	TimeFilter (double min_timestamp, double max_timestamp)
const int	GetNHitBars () const

Private Attributes
std::vector< DSHit >	hits_
Configuration	configuration_
int	station_

Detailed Description

Definition at line 12 of file sndDSPlane.h.

Constructor & Destructor Documentation

DSPlane()

snd::analysis_tools::DSPlane::DSPlane	(	std::vector< MuFilterHit * >	snd_hits,
		const Configuration &	configuration,
		MuFilter *	muon_filter_geometry,
		int	station
	)

Definition at line 13 of file sndDSPlane.cxx.

                                                                                                                                                   : configuration_(configuration), station_(station)
{
    for ( auto mu_hit : snd_hits)
    {
        TVector3 A, B;
        int detectorID = mu_hit->GetDetectorID();
        muon_filter_geometry->GetPosition(detectorID, A, B);
        const int n_sipms  = muon_filter_geometry->GetnSiPMs(detectorID);
        const int n_sides = muon_filter_geometry->GetnSides(detectorID);
        for (int i{0}; i < n_sipms * n_sides; ++i)
        {
            if (mu_hit->isMasked(i) || mu_hit->GetSignal(i) < -990.) continue;
            DSHit hit;
            hit.bar = static_cast<int>(detectorID % 1000);
            hit.channel_index = n_sipms * n_sides * hit.bar + i;
            hit.is_right = (i >= n_sipms) ? true : false;
            hit.timestamp = configuration_.is_mc ? mu_hit->GetTime(i) / ShipUnit::snd_TDC2ns : mu_hit->GetTime(i);
            hit.qdc = mu_hit->GetSignal(i);
            // use the left and right measurements to calculate the x coordinate along the bar
            if (!mu_hit->isVertical()) {
                hit.is_x = false;
                float tmp_x = mu_hit->GetImpactXpos(true, true, false, configuration_.is_mc);
                hit.x = (tmp_x < -990.) ? std::nan("") : A.X() - tmp_x; 
                hit.y = A.Y();
            }
            else {
                hit.is_x = true;
                hit.x = A.X();
                hit.y = std::nan("");
            }
            hit.z = A.Z();
            hits_.push_back(hit);
        }
    }
}

Member Function Documentation

GetBarNHits()

const int snd::analysis_tools::DSPlane::GetBarNHits

(

int

bar_to_compute

)

const

Definition at line 56 of file sndDSPlane.cxx.

{
    int bar_hit = std::count_if(hits_.begin(), hits_.end(),
                                [bar_to_compute](const auto &hit) {return hit.bar == bar_to_compute;});
    return bar_hit;
}

GetHits()

const std::vector< DSHit > snd::analysis_tools::DSPlane::GetHits ( ) const

inline

Definition at line 48 of file sndDSPlane.h.

{ return hits_; };

GetNHitBars()

const int snd::analysis_tools::DSPlane::GetNHitBars

(

)

const

Definition at line 71 of file sndDSPlane.cxx.

{
    int count{0};
    for (int bar{0}; bar < 2*configuration_.ds_bar_per_station; ++bar) {
        if (GetBarNHits(bar) > 0) count++;
    }
    return count;
}

GetNHits()

const int snd::analysis_tools::DSPlane::GetNHits ( ) const

inline

Definition at line 44 of file sndDSPlane.h.

{ return hits_.size(); };

GetStation()

const int snd::analysis_tools::DSPlane::GetStation ( ) const

inline

Definition at line 45 of file sndDSPlane.h.

{ return station_; }

GetTotQdc()

const double snd::analysis_tools::DSPlane::GetTotQdc

(

)

const

Definition at line 49 of file sndDSPlane.cxx.

{
    double tot_qdc = std::accumulate(hits_.begin(), hits_.end(), 0.0,
                               [](double sum, const auto &b) {return sum + b.qdc;});
    return tot_qdc;
}

TimeFilter()

void snd::analysis_tools::DSPlane::TimeFilter	(	double	min_timestamp,
		double	max_timestamp
	)

Definition at line 63 of file sndDSPlane.cxx.

{
    hits_.erase(std::remove_if(hits_.begin(), hits_.end(),
                               [&](auto &hit)
                               { return hit.timestamp < min_timestamp || hit.timestamp > max_timestamp; }),
                hits_.end());
}

Member Data Documentation

configuration_

Configuration snd::analysis_tools::DSPlane::configuration_

private

Definition at line 56 of file sndDSPlane.h.

hits_

std::vector<DSHit> snd::analysis_tools::DSPlane::hits_

private

Definition at line 55 of file sndDSPlane.h.

station_

int snd::analysis_tools::DSPlane::station_

private

Definition at line 57 of file sndDSPlane.h.

---

The documentation for this class was generated from the following files:

• analysis/tools/sndDSPlane.h
• analysis/tools/sndDSPlane.cxx