snd::analysis_tools::USPlane Class Reference

#include <sndUSPlane.h>

Collaboration diagram for snd::analysis_tools::USPlane:

Classes
struct	rl_pair
struct	sl_pair
struct	USHit

Public Member Functions
	USPlane (std::vector< MuFilterHit * > snd_hits, const Configuration &configuration, MuFilter *muon_filter_geometry, int station, bool use_small_sipms=false)
const sl_pair< int >	GetNHits () const
const int	GetStation () const
const sl_pair< double >	GetTotQdc () const
const sl_pair< double >	GetTotEnergy () const
const rl_pair< double >	GetSideQdc () const
const rl_pair< double >	GetBarQdc (int bar_to_compute) const
const sl_pair< int >	GetBarNHits (int bar_to_compute) const
const std::vector< USHit >	GetHits () const
double	HasShower () const
void	FindCentroid ()
ROOT::Math::XYZPoint	GetCentroid () const
ROOT::Math::XYZPoint	GetCentroidError () const
const int	GetNHitBars () const
void	TimeFilter (double min_timestamp, double max_timestamp)

Private Attributes
std::vector< USHit >	hits_
Configuration	configuration_
ROOT::Math::XYZPoint	centroid_
ROOT::Math::XYZPoint	centroid_error_
int	station_

Detailed Description

Definition at line 13 of file sndUSPlane.h.

Constructor & Destructor Documentation

USPlane()

snd::analysis_tools::USPlane::USPlane	(	std::vector< MuFilterHit * >	snd_hits,
		const Configuration &	configuration,
		MuFilter *	muon_filter_geometry,
		int	station,
		bool	use_small_sipms = false
	)

Definition at line 13 of file sndUSPlane.cxx.

                                                                                                                                                                         : configuration_(configuration), centroid_(std::nan(""), std::nan(""), std::nan("")), centroid_error_(std::nan(""), std::nan(""),std::nan("")), station_(station)
{
    for ( auto mu_hit : snd_hits)
    {
        TVector3 A, B;
        int detectorID = mu_hit->GetDetectorID();
        muon_filter_geometry->GetPosition(detectorID, A, B);
        for (int i{0}; i < 16; ++i)
        {
            if (mu_hit->isMasked(i) || mu_hit->GetSignal(i) < -990.) continue;
            USHit hit;
            hit.bar = static_cast<int>(detectorID % 1000);
            hit.channel_index = 16 * hit.bar + i;
            hit.is_large = !mu_hit->isShort(i);
            hit.is_right = i > 7 ? true : false;
 
            if (!hit.is_large && !use_small_sipms)
            {
                hit.timestamp = std::nan("");
                hit.qdc = std::nan("");
                hit.x = std::nan("");
                hit.y = std::nan("");
                hit.z = std::nan("");
            }
            else
            {
                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
                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();
                hit.z = A.Z();
            }
            hits_.push_back(hit);
        }
    }
}

Member Function Documentation

FindCentroid()

void snd::analysis_tools::USPlane::FindCentroid

(

)

Definition at line 52 of file sndUSPlane.cxx.

{
    // select large SiPM channels with positive qdc
    std::vector<USHit> cleaned_hits = hits_;
 
    cleaned_hits.erase(std::remove_if(cleaned_hits.begin(), cleaned_hits.end(),
                               [&](auto &hit)
                               { return (hit.qdc <= 0 || hit.is_large == false || std::isnan(hit.x)); }),
                       cleaned_hits.end());
 
    // min number of hit in the plane to attempt to find a centroid
    if (static_cast<int>(cleaned_hits.size()) < configuration_.us_min_n_hits_for_centroid)
    {
        return;
    }
 
    // weigthed sum calculated per plane
    double  weighted_sum_x{0.0}, weighted_sum_y{0.0}, weighted_sum_z{0.0};
    double total_qdc_positive{0.0}, sum_qdc2_positive{0.0};
    // loop over hits in the plane
    for (const auto &hit : cleaned_hits)
    {
        weighted_sum_x += hit.x * hit.qdc;
        weighted_sum_y += hit.y * hit.qdc;
        weighted_sum_z += hit.z * hit.qdc;
        total_qdc_positive += hit.qdc;
        sum_qdc2_positive += hit.qdc*hit.qdc;
    }
    weighted_sum_x /= total_qdc_positive;
    weighted_sum_y /= total_qdc_positive;
    weighted_sum_z /= total_qdc_positive;
    centroid_.SetXYZ(weighted_sum_x, weighted_sum_y, weighted_sum_z);
    auto qdc_error_scaler = sqrt(sum_qdc2_positive)/total_qdc_positive;
    centroid_error_.SetXYZ(configuration_.us_centroid_error_x*qdc_error_scaler,
                           configuration_.us_centroid_error_y*qdc_error_scaler,
                           configuration_.us_centroid_error_z*qdc_error_scaler);
}

GetBarNHits()

const snd::analysis_tools::USPlane::sl_pair< int > snd::analysis_tools::USPlane::GetBarNHits

(

int

bar_to_compute

)

const

Definition at line 152 of file sndUSPlane.cxx.

{
    sl_pair<int> bar_hit{0, 0};
    for (const auto &hit : hits_)
    {
        if (hit.bar != bar_to_compute)
            continue;
        else
        {
            if (hit.is_large)
            
                    bar_hit.large++;
            else
                    bar_hit.small++;
            
        }
    }
    return bar_hit;
}

GetBarQdc()

const snd::analysis_tools::USPlane::rl_pair< double > snd::analysis_tools::USPlane::GetBarQdc

(

int

bar_to_compute

)

const

Definition at line 131 of file sndUSPlane.cxx.

{
    rl_pair<double> bar_qdc{0.0, 0.0};
    for (const auto &hit : hits_)
    {
        if (hit.bar != bar_to_compute)
            continue;
        else
        {
            if (hit.is_large)
            {
                if (hit.is_right)
                    bar_qdc.right += hit.qdc;
                else
                    bar_qdc.left += hit.qdc;
            }
        }
    }
    return bar_qdc;
}

GetCentroid()

ROOT::Math::XYZPoint snd::analysis_tools::USPlane::GetCentroid ( ) const

inline

Definition at line 64 of file sndUSPlane.h.

{ return centroid_; }

GetCentroidError()

ROOT::Math::XYZPoint snd::analysis_tools::USPlane::GetCentroidError ( ) const

inline

Definition at line 65 of file sndUSPlane.h.

{ return centroid_error_; }

GetHits()

const std::vector< USHit > snd::analysis_tools::USPlane::GetHits ( ) const

inline

Definition at line 60 of file sndUSPlane.h.

{ return hits_; };

GetNHitBars()

const int snd::analysis_tools::USPlane::GetNHitBars

(

)

const

Definition at line 190 of file sndUSPlane.cxx.

{
    int count{0};
    for (int bar{0}; bar < configuration_.us_bar_per_station; ++bar) {
        if (GetBarNHits(bar).large > configuration_.us_min_hit_on_bar) count++;
    }
    return count;
}

GetNHits()

const snd::analysis_tools::USPlane::sl_pair< int > snd::analysis_tools::USPlane::GetNHits

(

)

const

Definition at line 180 of file sndUSPlane.cxx.

{
    sl_pair<int> counts{0, 0};
    counts.large = std::count_if(hits_.begin(), hits_.end(), [](auto &hit)
                             { return hit.is_large; });
    counts.small = static_cast<int>(hits_.size()) - counts.large;
 
    return counts;
}

GetSideQdc()

const snd::analysis_tools::USPlane::rl_pair< double > snd::analysis_tools::USPlane::GetSideQdc

(

)

const

Definition at line 115 of file sndUSPlane.cxx.

{
    rl_pair<double> side_qdc{0.0, 0.0};
    for (const auto &hit : hits_)
    {
        if (hit.is_large)
        {
            if (hit.is_right)
                side_qdc.right += hit.qdc;
            else
                side_qdc.left += hit.qdc;
        }
    }
    return side_qdc;
}

GetStation()

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

inline

Definition at line 53 of file sndUSPlane.h.

{ return station_; }

GetTotEnergy()

const snd::analysis_tools::USPlane::sl_pair< double > snd::analysis_tools::USPlane::GetTotEnergy

(

)

const

Definition at line 103 of file sndUSPlane.cxx.

{
    sl_pair<double> tot_energy{0.0, 0.0};
    sl_pair<double> tot_qdc = GetTotQdc();
 
    tot_energy.large = tot_qdc.large *configuration_.us_qdc_to_gev; 
    tot_energy.small = tot_qdc.small *configuration_.us_qdc_to_gev;
 
    return tot_energy;
}

GetTotQdc()

const snd::analysis_tools::USPlane::sl_pair< double > snd::analysis_tools::USPlane::GetTotQdc

(

)

const

Definition at line 90 of file sndUSPlane.cxx.

{
    sl_pair<double> totQdc{0.0, 0.0};
    for (const auto &hit : hits_)
    {
        if (hit.is_large)
            totQdc.large += hit.qdc;
        else
            totQdc.small += hit.qdc;
    }
    return totQdc;
}

HasShower()

double snd::analysis_tools::USPlane::HasShower ( ) const

inline

Definition at line 61 of file sndUSPlane.h.

{ return GetNHits().large >= configuration_.us_min_n_hits_for_centroid; }

TimeFilter()

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

Definition at line 172 of file sndUSPlane.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

centroid_

ROOT::Math::XYZPoint snd::analysis_tools::USPlane::centroid_

private

Definition at line 74 of file sndUSPlane.h.

centroid_error_

ROOT::Math::XYZPoint snd::analysis_tools::USPlane::centroid_error_

private

Definition at line 75 of file sndUSPlane.h.

configuration_

Configuration snd::analysis_tools::USPlane::configuration_

private

Definition at line 73 of file sndUSPlane.h.

hits_

std::vector<USHit> snd::analysis_tools::USPlane::hits_

private

Definition at line 72 of file sndUSPlane.h.

station_

int snd::analysis_tools::USPlane::station_

private

Definition at line 76 of file sndUSPlane.h.

---

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

• analysis/tools/sndUSPlane.h
• analysis/tools/sndUSPlane.cxx