#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 20 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 7 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 91 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);
                 std::vector<float> qTemp;
                 qTemp.push_back(stripChargeVal);
                 m_v_stripTotalCharge.push_back(qTemp);
                 std::vector<float> tTemp;
                 tTemp.push_back(timeBin * m_timeResolution);
                 m_v_stripTimeThreshold.push_back(tTemp);
             }
         }
     }
 }

◆ calculateTimeSeries()

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

Definition at line 30 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 22 of file MM_StripResponse.cxx.

22 { return m_Electrons; }

◆ getMaxCharge()

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

Definition at line 134 of file MM_StripResponse.cxx.

134 { return m_stripMaxCharge; }

◆ getMaxChargeVec()

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

Definition at line 139 of file MM_StripResponse.cxx.

139 { return m_v_stripMaxCharge; }

◆ getNElectrons()

int MM_StripResponse::getNElectrons ( ) const

Definition at line 14 of file MM_StripResponse.cxx.

14 { return m_Electrons.size(); }

◆ getStripVec()

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

Definition at line 136 of file MM_StripResponse.cxx.

136 { return m_v_strip; }

◆ getTimeMaxCharge()

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

Definition at line 135 of file MM_StripResponse.cxx.

135 { return m_stripTimeMaxCharge; }

◆ getTimeMaxChargeVec()

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

Definition at line 140 of file MM_StripResponse.cxx.

140 { return m_v_stripTimeMaxCharge; }

◆ getTimeThreshold()

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

Definition at line 132 of file MM_StripResponse.cxx.

132 { return m_stripTimeThreshold; }

◆ getTimeThresholdVec()

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

Definition at line 137 of file MM_StripResponse.cxx.

137 { return m_v_stripTimeThreshold; }

◆ getTotalCharge()

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

Definition at line 133 of file MM_StripResponse.cxx.

133 { return m_stripTotalCharge; }

◆ getTotalChargeVec()

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

Definition at line 138 of file MM_StripResponse.cxx.

138 { return m_v_stripTotalCharge; }

◆ simulateCrossTalk()

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

Definition at line 50 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 24 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 16 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 52 of file MM_StripResponse.h.

◆ m_maxstripID

int MM_StripResponse::m_maxstripID {0}

private

Definition at line 50 of file MM_StripResponse.h.

◆ m_minstripID

int MM_StripResponse::m_minstripID {0}

private

Definition at line 49 of file MM_StripResponse.h.

◆ m_stripCharges

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

private

Definition at line 55 of file MM_StripResponse.h.

◆ m_stripID

int MM_StripResponse::m_stripID {0}

private

Definition at line 48 of file MM_StripResponse.h.

◆ m_stripMaxCharge

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

private

Definition at line 60 of file MM_StripResponse.h.

◆ m_stripPitch

float MM_StripResponse::m_stripPitch {0.f}

private

Definition at line 47 of file MM_StripResponse.h.

◆ m_stripTimeMaxCharge

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

private

Definition at line 61 of file MM_StripResponse.h.

◆ m_stripTimeThreshold

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

private

Definition at line 58 of file MM_StripResponse.h.

◆ m_stripTotalCharge

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

private

Definition at line 59 of file MM_StripResponse.h.

◆ m_timeResolution

float MM_StripResponse::m_timeResolution {0.f}

private

Definition at line 46 of file MM_StripResponse.h.

◆ m_v_strip

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

private

Definition at line 64 of file MM_StripResponse.h.

◆ m_v_stripMaxCharge

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

private

Definition at line 67 of file MM_StripResponse.h.

◆ m_v_stripTimeMaxCharge

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

private

Definition at line 68 of file MM_StripResponse.h.

◆ m_v_stripTimeThreshold

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

private

Definition at line 65 of file MM_StripResponse.h.

◆ m_v_stripTotalCharge

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

private

Definition at line 66 of file MM_StripResponse.h.

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