#include <MM_StripResponse.h>

Collaboration diagram for MM_StripResponse:

Public Member Functions
	MM_StripResponse ()=default
	MM_StripResponse (std::vector< std::unique_ptr< MM_IonizationCluster > > &IonizationClusters, float timeResolution, float stripPitch, int stripID, int minstripID, int maxstripID)
void	timeOrderElectrons ()
void	calculateTimeSeries (float thetaD, int gasgap)
void	simulateCrossTalk (float crossTalk1, float crossTalk2)
void	calculateSummaries (float chargeThreshold)
const std::map< int, int > &	getTimeThreshold () const
const std::map< int, float > &	getTotalCharge () const
const std::map< int, float > &	getMaxCharge () const
const std::map< int, int > &	getTimeMaxCharge () const
const std::vector< int > &	getStripVec () const
const std::vector< std::vector< float > > &	getTimeThresholdVec () const
const std::vector< std::vector< float > > &	getTotalChargeVec () const
const std::vector< float > &	getMaxChargeVec () const
const std::vector< float > &	getTimeMaxChargeVec () const
int	getNElectrons () const
float	totalCharge () const
std::vector< std::unique_ptr< MM_Electron > > &	getElectrons ()

Private Attributes
float	m_timeResolution {0.f}
float	m_stripPitch {0.f}
int	m_stripID {0}
int	m_minstripID {0}
int	m_maxstripID {0}
std::vector< std::unique_ptr< MM_Electron > >	m_Electrons {}
std::map< int, std::map< int, float > >	m_stripCharges {}
std::map< int, int >	m_stripTimeThreshold {}
std::map< int, float >	m_stripTotalCharge {}
std::map< int, float >	m_stripMaxCharge {}
std::map< int, int >	m_stripTimeMaxCharge {}
std::vector< int >	m_v_strip {}
std::vector< std::vector< float > >	m_v_stripTimeThreshold {}
std::vector< std::vector< float > >	m_v_stripTotalCharge {}
std::vector< float >	m_v_stripMaxCharge {}
std::vector< float >	m_v_stripTimeMaxCharge {}

Detailed Description

Definition at line 21 of file MM_StripResponse.h.

Constructor & Destructor Documentation

◆ MM_StripResponse() [1/2]

MM_StripResponse::MM_StripResponse ( )

default

◆ MM_StripResponse() [2/2]

MM_StripResponse::MM_StripResponse	(	std::vector< std::unique_ptr< MM_IonizationCluster > > &	IonizationClusters,
		float	timeResolution,
		float	stripPitch,
		int	stripID,
		int	minstripID,
		int	maxstripID )

Definition at line 11 of file MM_StripResponse.cxx.

                                                                                                  :
    m_timeResolution(timeResolution), m_stripPitch(stripPitch), m_stripID(stripID), m_minstripID(minstripID), m_maxstripID(maxstripID) {
    for (auto& IonizationCluster : IonizationClusters)
        for (auto& Electron : IonizationCluster->getElectrons()) m_Electrons.push_back(std::move(Electron));
}

Member Function Documentation

◆ calculateSummaries()

void MM_StripResponse::calculateSummaries ( float chargeThreshold )

Definition at line 95 of file MM_StripResponse.cxx.

                                                               {
    std::map<int, std::map<int, float>> stripChargesCopy2;
    stripChargesCopy2.insert(m_stripCharges.begin(), m_stripCharges.end());
    for (auto& stripTimeSeries : stripChargesCopy2) {
        int timeBin = stripTimeSeries.first;
        for (auto& stripCharge : stripTimeSeries.second) {
            int stripVal = stripCharge.first;
            // remove dead (missing) strips
            // First active strip starts at m_minstripID
            // Last active strip numbrer is maxStripID-1
            if (stripVal < m_minstripID || stripVal > m_maxstripID - 1) continue;
            // remove PCB gap strips
            if (stripVal == 1024 || stripVal == 1025 || stripVal == 2048 || stripVal == 2049 || stripVal == 3072 || stripVal == 3073 ||
                stripVal == 4096 || stripVal == 4097)
                continue;
            float stripChargeVal = stripCharge.second;
            if (stripChargeVal < chargeThreshold) continue;
 
            bool found = false;
            for (size_t ii = 0; ii < m_v_strip.size(); ii++) {
                if (m_v_strip.at(ii) == stripVal) {
                    m_v_stripTimeThreshold.at(ii).push_back(timeBin * m_timeResolution);
                    m_v_stripTotalCharge.at(ii).push_back(stripChargeVal);
                    found = true;
                    break;
                }
            }
            if (!found) {  //   // strip not in vector, add new entry
                m_v_strip.push_back(stripVal);
                //construct vector of 1 value of stripChargeVal in-place
                m_v_stripTotalCharge.emplace_back(1,stripChargeVal);
                //construct vector of 1 value of timeBin * m_timeResolution in-place
                m_v_stripTimeThreshold.emplace_back(1, timeBin * m_timeResolution);
            }
        }
    }
}

◆ calculateTimeSeries()

void MM_StripResponse::calculateTimeSeries	(	float	thetaD,
		int	gasgap )

Definition at line 34 of file MM_StripResponse.cxx.

                                                         {
    for (auto& Electron : m_Electrons) {
        int timeBin = (int)(Electron->getTime() / m_timeResolution);
        // m_stripID defines the initial strip where the muon entered the gas gap
 
        int stripVal = 0;
        if (std::abs(Electron->getX()) > m_stripPitch / 2) {
            if (Electron->getX() > 0.0)
                stripVal = m_stripID + int((Electron->getX() - m_stripPitch / 2) / m_stripPitch) + 1;
            else
                stripVal = m_stripID + int((Electron->getX() + m_stripPitch / 2) / m_stripPitch) - 1;
        } else
            stripVal = m_stripID;
 
        // Only add the strips that are either read out, or can cross talk to the read out strips
        if (stripVal < m_minstripID - 2 || stripVal > m_maxstripID + 1) stripVal = -1;
        if (stripVal > 0) m_stripCharges[timeBin][stripVal] += Electron->getCharge();
    }
}

◆ getElectrons()

std::vector< std::unique_ptr< MM_Electron > > & MM_StripResponse::getElectrons ( )

Definition at line 26 of file MM_StripResponse.cxx.

26{ return m_Electrons; }

◆ getMaxCharge()

const std::map< int, float > & MM_StripResponse::getMaxCharge ( ) const

Definition at line 136 of file MM_StripResponse.cxx.

136{ return m_stripMaxCharge; }

MM_StripResponse::m_stripMaxCharge

std::map< int, float > m_stripMaxCharge

Definition MM_StripResponse.h:61

◆ getMaxChargeVec()

const std::vector< float > & MM_StripResponse::getMaxChargeVec ( ) const

Definition at line 141 of file MM_StripResponse.cxx.

141{ return m_v_stripMaxCharge; }

MM_StripResponse::m_v_stripMaxCharge

std::vector< float > m_v_stripMaxCharge

Definition MM_StripResponse.h:68

◆ getNElectrons()

int MM_StripResponse::getNElectrons ( ) const

Definition at line 18 of file MM_StripResponse.cxx.

18{ return m_Electrons.size(); }

◆ getStripVec()

const std::vector< int > & MM_StripResponse::getStripVec ( ) const

Definition at line 138 of file MM_StripResponse.cxx.

138{ return m_v_strip; }

◆ getTimeMaxCharge()

const std::map< int, int > & MM_StripResponse::getTimeMaxCharge ( ) const

Definition at line 137 of file MM_StripResponse.cxx.

137{ return m_stripTimeMaxCharge; }

MM_StripResponse::m_stripTimeMaxCharge

std::map< int, int > m_stripTimeMaxCharge

Definition MM_StripResponse.h:62

◆ getTimeMaxChargeVec()

const std::vector< float > & MM_StripResponse::getTimeMaxChargeVec ( ) const

Definition at line 142 of file MM_StripResponse.cxx.

142{ return m_v_stripTimeMaxCharge; }

MM_StripResponse::m_v_stripTimeMaxCharge

std::vector< float > m_v_stripTimeMaxCharge

Definition MM_StripResponse.h:69

◆ getTimeThreshold()

const std::map< int, int > & MM_StripResponse::getTimeThreshold ( ) const

Definition at line 134 of file MM_StripResponse.cxx.

134{ return m_stripTimeThreshold; }

MM_StripResponse::m_stripTimeThreshold

std::map< int, int > m_stripTimeThreshold

Definition MM_StripResponse.h:59

◆ getTimeThresholdVec()

const std::vector< std::vector< float > > & MM_StripResponse::getTimeThresholdVec ( ) const

Definition at line 139 of file MM_StripResponse.cxx.

139{ return m_v_stripTimeThreshold; }

◆ getTotalCharge()

const std::map< int, float > & MM_StripResponse::getTotalCharge ( ) const

Definition at line 135 of file MM_StripResponse.cxx.

135{ return m_stripTotalCharge; }

MM_StripResponse::m_stripTotalCharge

std::map< int, float > m_stripTotalCharge

Definition MM_StripResponse.h:60

◆ getTotalChargeVec()

const std::vector< std::vector< float > > & MM_StripResponse::getTotalChargeVec ( ) const

Definition at line 140 of file MM_StripResponse.cxx.

140{ return m_v_stripTotalCharge; }

◆ simulateCrossTalk()

void MM_StripResponse::simulateCrossTalk	(	float	crossTalk1,
		float	crossTalk2 )

Definition at line 54 of file MM_StripResponse.cxx.

                                                                           {
    // if no cross talk is simulate just skip everything and keep m_stripCharges as it was
    if (crossTalk1 > 0.) {
        // Unfortunately get stuck in the loop if you edit the map in the loop
        //     So make a copy!
 
        std::map<int, std::map<int, float>> stripChargesCopy1;
        stripChargesCopy1.insert(m_stripCharges.begin(), m_stripCharges.end());
 
        // clear strip charge map since charge on the main strip needs to be scaled
        m_stripCharges.clear();
 
        for (auto& stripTimeSeries : stripChargesCopy1) {
            int timeBin = stripTimeSeries.first;
            for (auto& stripCharge : stripTimeSeries.second) {
                int stripVal = stripCharge.first;
                float stripChargeVal = stripCharge.second;
 
                if (stripChargeVal == 0.) continue;
 
                // scale factcor for the charge on the main strip to account for the cross talk charge
                float chargeScaleFactor = 1.0 / (1. + ((stripVal - 1 > 0) + (stripVal + 1 < m_maxstripID)) * crossTalk1 +
                                                 ((stripVal - 2 > 0) + (stripVal + 2 < m_maxstripID)) * crossTalk2);
                stripChargeVal *= chargeScaleFactor;
 
                m_stripCharges[timeBin][stripVal] += stripChargeVal;
 
                // Allow crosstalk between strips that exist.
                // Will check for read out strips in calculateSummaries function
                if (stripVal - 1 > 0) m_stripCharges[timeBin][stripVal - 1] += stripChargeVal * crossTalk1;
                if (stripVal + 1 < m_maxstripID) m_stripCharges[timeBin][stripVal + 1] += stripChargeVal * crossTalk1;
 
                if (crossTalk2 > 0.) {
                    if (stripVal - 2 > 0) m_stripCharges[timeBin][stripVal - 2] += stripChargeVal * crossTalk2;
                    if (stripVal + 2 < m_maxstripID) m_stripCharges[timeBin][stripVal + 2] += stripChargeVal * crossTalk2;
                }
            }
        }
    }
}

◆ timeOrderElectrons()

void MM_StripResponse::timeOrderElectrons ( )

Definition at line 28 of file MM_StripResponse.cxx.

                                          {
    std::sort(
        m_Electrons.begin(), m_Electrons.end(),
        [](const std::unique_ptr<MM_Electron>& a, const std::unique_ptr<MM_Electron>& b) -> bool { return a->getTime() < b->getTime(); });
}

◆ totalCharge()

float MM_StripResponse::totalCharge ( ) const

Definition at line 20 of file MM_StripResponse.cxx.

                                          {
    float qtot = 0;
    for (const std::unique_ptr<MM_Electron>& electron : m_Electrons) { qtot += electron->getCharge(); }
    return qtot;
}

Member Data Documentation

◆ m_Electrons

std::vector<std::unique_ptr<MM_Electron> > MM_StripResponse::m_Electrons {}

private

Definition at line 53 of file MM_StripResponse.h.

53{};

◆ m_maxstripID

int MM_StripResponse::m_maxstripID {0}

private

Definition at line 51 of file MM_StripResponse.h.

51{0};

◆ m_minstripID

int MM_StripResponse::m_minstripID {0}

private

Definition at line 50 of file MM_StripResponse.h.

50{0};

◆ m_stripCharges

std::map<int, std::map<int, float> > MM_StripResponse::m_stripCharges {}

private

Definition at line 56 of file MM_StripResponse.h.

56{};

◆ m_stripID

int MM_StripResponse::m_stripID {0}

private

Definition at line 49 of file MM_StripResponse.h.

49{0};

◆ m_stripMaxCharge

std::map<int, float> MM_StripResponse::m_stripMaxCharge {}

private

Definition at line 61 of file MM_StripResponse.h.

61{};

◆ m_stripPitch

float MM_StripResponse::m_stripPitch {0.f}

private

Definition at line 48 of file MM_StripResponse.h.

48{0.f};

◆ m_stripTimeMaxCharge

std::map<int, int> MM_StripResponse::m_stripTimeMaxCharge {}

private

Definition at line 62 of file MM_StripResponse.h.

62{};

◆ m_stripTimeThreshold

std::map<int, int> MM_StripResponse::m_stripTimeThreshold {}

private

Definition at line 59 of file MM_StripResponse.h.

59{};

◆ m_stripTotalCharge

std::map<int, float> MM_StripResponse::m_stripTotalCharge {}

private

Definition at line 60 of file MM_StripResponse.h.

60{};

◆ m_timeResolution

float MM_StripResponse::m_timeResolution {0.f}

private

Definition at line 47 of file MM_StripResponse.h.

47{0.f};

◆ m_v_strip

std::vector<int> MM_StripResponse::m_v_strip {}

private

Definition at line 65 of file MM_StripResponse.h.

65{};

◆ m_v_stripMaxCharge

std::vector<float> MM_StripResponse::m_v_stripMaxCharge {}

private

Definition at line 68 of file MM_StripResponse.h.

68{};

◆ m_v_stripTimeMaxCharge

std::vector<float> MM_StripResponse::m_v_stripTimeMaxCharge {}

private

Definition at line 69 of file MM_StripResponse.h.

69{};

◆ m_v_stripTimeThreshold

std::vector<std::vector<float> > MM_StripResponse::m_v_stripTimeThreshold {}

private

Definition at line 66 of file MM_StripResponse.h.

66{};

◆ m_v_stripTotalCharge

std::vector<std::vector<float> > MM_StripResponse::m_v_stripTotalCharge {}

private

Definition at line 67 of file MM_StripResponse.h.

67{};

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MM_StripResponse() [1/2]

◆ MM_StripResponse() [2/2]

Member Function Documentation

◆ calculateSummaries()

◆ calculateTimeSeries()

◆ getElectrons()

◆ getMaxCharge()

◆ getMaxChargeVec()

◆ getNElectrons()

◆ getStripVec()

◆ getTimeMaxCharge()

◆ getTimeMaxChargeVec()

◆ getTimeThreshold()

◆ getTimeThresholdVec()

◆ getTotalCharge()

◆ getTotalChargeVec()

◆ simulateCrossTalk()

◆ timeOrderElectrons()

◆ totalCharge()

Member Data Documentation

◆ m_Electrons

◆ m_maxstripID

◆ m_minstripID

◆ m_stripCharges

◆ m_stripID

◆ m_stripMaxCharge

◆ m_stripPitch

◆ m_stripTimeMaxCharge

◆ m_stripTimeThreshold

◆ m_stripTotalCharge

◆ m_timeResolution

◆ m_v_strip

◆ m_v_stripMaxCharge

◆ m_v_stripTimeMaxCharge

◆ m_v_stripTimeThreshold

◆ m_v_stripTotalCharge