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 (TClonesArray *snd_hits, Configuration configuration, MuFilter *muon_filter_geometry, int index_begin, int index_end, int station)
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
void	FindCentroid ()
const TVector3	GetCentroid () const
const TVector3	GetCentroidError () const
void	TimeFilter (double min_timestamp, double max_timestamp)

Private Attributes
std::vector< USHit >	hits_
Configuration	configuration_
TVector3	centroid_
TVector3	centroid_error_
int	station_

Detailed Description

Definition at line 13 of file sndUSPlane.h.

Constructor & Destructor Documentation

USPlane()

snd::analysis_tools::USPlane::USPlane	(	TClonesArray *	snd_hits,
		snd::Configuration	configuration,
		MuFilter *	muon_filter_geometry,
		int	index_begin,
		int	index_end,
		int	station
	)

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)
{
    if (index_begin > index_end)
    {
        throw std::runtime_error{"Begin index > end index"};
    }
    for (int j{index_begin}; j < index_end; ++j)
    {
        auto mu_hit = static_cast<MuFilterHit *>(snd_hits->At(j));
        
        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>(mu_hit->GetDetectorID() % 1000);
            hit.channel_index = 16 * hit.bar + i;
            hit.timestamp = mu_hit->GetTime(i);
            hit.qdc = mu_hit->GetSignal(i);
            hit.is_large = !mu_hit->isShort(i);
        
            TVector3 A, B;
            int detectorID = mu_hit->GetDetectorID();
            muon_filter_geometry->GetPosition(detectorID, A, B);
            hit.is_right = i > 7 ? true : false;
            hit.x = hit.is_right ? B.X() : A.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 45 of file sndUSPlane.cxx.

{
    // min number of hit in the plane to attempt to find a centroid
    if (hits_.size() < configuration_.us_min_n_hits_for_centroid)
    {
        // std::cout<<"Not enough hits in US plane " << station_ <<" to find centroid\n";
        return;
    }
    
    double total_qdc = GetTotQdc().large;
    if (total_qdc > 0.0)
    {
        for (const auto &hit : hits_)
        {
            // weigthed sum
            double  weighted_sum_x{0.0}, weighted_sum_y{0.0}, weighted_sum_z{0.0};
            double total_qdc_positive{0.0};
 
            if (hit.qdc > 0.0)
            {
                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;
            }
 
            double x = weighted_sum_x / total_qdc_positive;
            double y = weighted_sum_y / total_qdc_positive;
            double z = weighted_sum_z / total_qdc_positive;
            centroid_.SetXYZ(x, y, z);
        }
    }
    centroid_error_.SetXYZ(configuration_.us_centroid_error_x, configuration_.us_centroid_error_y, configuration_.us_centroid_error_z);
}

GetBarNHits()

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

(

int

bar_to_compute

)

const

Definition at line 142 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 121 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()

const TVector3 snd::analysis_tools::USPlane::GetCentroid ( ) const

inline

Definition at line 61 of file sndUSPlane.h.

{ return centroid_; };

GetCentroidError()

const TVector3 snd::analysis_tools::USPlane::GetCentroidError ( ) const

inline

Definition at line 62 of file sndUSPlane.h.

{ return centroid_error_; };

GetHits()

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

inline

Definition at line 58 of file sndUSPlane.h.

{ return hits_; };

GetNHits()

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

(

)

const

Definition at line 170 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 = hits_.size() - counts.large;
 
    return counts;
}

GetSideQdc()

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

(

)

const

Definition at line 105 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 51 of file sndUSPlane.h.

{ return station_; };

GetTotEnergy()

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

(

)

const

Definition at line 93 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 80 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;
}

TimeFilter()

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

Definition at line 162 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_

TVector3 snd::analysis_tools::USPlane::centroid_

private

Definition at line 69 of file sndUSPlane.h.

centroid_error_

TVector3 snd::analysis_tools::USPlane::centroid_error_

private

Definition at line 70 of file sndUSPlane.h.

configuration_

Configuration snd::analysis_tools::USPlane::configuration_

private

Definition at line 68 of file sndUSPlane.h.

hits_

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

private

Definition at line 67 of file sndUSPlane.h.

station_

int snd::analysis_tools::USPlane::station_

private

Definition at line 71 of file sndUSPlane.h.

---

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

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