MDT_Response Class Reference

#include <MDT_Response.h>

Collaboration diagram for MDT_Response:

Public Types
typedef std::vector< double >	clusterVec

Public Member Functions
	MDT_Response ()
	MDT_Response (double timewindow, double binsize)
	~MDT_Response ()
void	SetSegment (double r, double x)
void	SetTubeRadius (double radius)
void	SetClusterDensity (double dens)
void	SetAttLength (double len)
void	SetRtMode (unsigned int mode)
void	SetDifRMS (double rms)
void	SetTriggerElectron (double el)
void	SetIntegrationWindow (double win)
void	SetRtParameters (int npar, double *par)
bool	GetSignal (CLHEP::HepRandomEngine *rndmEngine)
bool	GetSignal (double ParticleCharge, double ParticleGamma, CLHEP::HepRandomEngine *rndmEngine)
double	DriftTime () const
double	TThreshold () const
double	TFirst () const
double	T0 () const
double	AdcResponse () const
double	Charge () const
double	V_r (double r)
double	DoStep (CLHEP::HepRandomEngine *rndmEngine) const
double	DampingFactor (double x)
double	PropagationDelay (double x)
double	RtoT (double r)
double	Diffusion (double r, CLHEP::HepRandomEngine *rndmEngine) const
double	SigmaDiffusion (double r) const
int	GenerateQ (CLHEP::HepRandomEngine *rndmEngine) const
const double *	RtParameters () const

Private Member Functions
void	InitTubeParameters ()
void	InitClusters (double timewindow, double binsize)
void	InitRt ()
void	InitdEdxTable ()
void	DoStepping (CLHEP::HepRandomEngine *rndmEngine)
void	DoStepping (double ParticleCharge, double ParticleGamma, CLHEP::HepRandomEngine *rndmEngine)
void	Reset ()

Private Attributes
double	m_radius = 0.0
double	m_rhit
double	m_xhit
double	m_pathLength
double	m_clusterDensity = 0.0
std::vector< double >	m_gammaFactorVec
std::vector< double >	m_numberOfClustersPerMmVec
double	m_attLength = 0.0
double	m_signalSpeed = 0.0
int	m_rtMode = 0
unsigned int	m_rtNpar = 0U
double *	m_rtParameters
double	m_difSmearing = 0.0
double	m_triggerElectron = 0.0
double	m_integrationWindow = 0.0
double	m_binsize = 0.0
double	m_timeWindow = 0.0
int	m_offset = 0
int	m_bins = 0
Amplifier	m_amplifier
clusterVec	m_clusters
double	m_t0

Detailed Description

Definition at line 20 of file MDT_Response.h.

Member Typedef Documentation

clusterVec

typedef std::vector<double> MDT_Response::clusterVec

Definition at line 23 of file MDT_Response.h.

Constructor & Destructor Documentation

MDT_Response() [1/2]

MDT_Response::MDT_Response

(

)

Definition at line 8 of file MDT_Response.cxx.

                           : 
  m_rhit(0.0), 
  m_xhit(0.0), 
  m_pathLength(0.0), 
  m_rtParameters(nullptr),
  m_t0(0.0)
{
  InitTubeParameters();
  InitClusters(200.,0.2);
  InitdEdxTable(); 
}

MDT_Response() [2/2]

MDT_Response::MDT_Response	(	double	timewindow,
		double	binsize )

Definition at line 20 of file MDT_Response.cxx.

                                                            : 
  m_rhit(0.0), 
  m_xhit(0.0), 
  m_pathLength(0.0), 
  m_rtParameters(nullptr),
  m_t0(0.0)
{
  InitTubeParameters();
  InitClusters(timewindow,binsize);
  InitdEdxTable(); 
}

~MDT_Response()

MDT_Response::~MDT_Response

(

)

Definition at line 32 of file MDT_Response.cxx.

{
  delete[] m_rtParameters; m_rtParameters=nullptr;
}

Member Function Documentation

AdcResponse()

double MDT_Response::AdcResponse ( ) const

inline

Definition at line 191 of file MDT_Response.h.

{ return m_amplifier.AdcResponse(); } 

Charge()

double MDT_Response::Charge ( ) const

inline

Definition at line 194 of file MDT_Response.h.

{ return m_amplifier.Charge(); }

DampingFactor()

double MDT_Response::DampingFactor ( double x )

inline

Definition at line 136 of file MDT_Response.h.

{ return exp(-1.*x/m_attLength); }

Diffusion()

double MDT_Response::Diffusion ( double r, CLHEP::HepRandomEngine * rndmEngine ) const

inline

Definition at line 161 of file MDT_Response.h.

{
  // Smearing with maximum of three sigma 
  double sigma = SigmaDiffusion(r);
  double t = CLHEP::RandGaussZiggurat::shoot(rndmEngine,0.,sigma); 
  if(fabs(t) > 3*sigma) {
    t = Diffusion(r, rndmEngine);
  }
  return t;
}

DoStep()

double MDT_Response::DoStep ( CLHEP::HepRandomEngine * rndmEngine ) const

inline

Definition at line 114 of file MDT_Response.h.

{
  return (-1./m_clusterDensity)*log( CLHEP::RandFlat::shoot(rndmEngine) );
}

DoStepping() [1/2]

void MDT_Response::DoStepping ( CLHEP::HepRandomEngine * rndmEngine )

private

Definition at line 151 of file MDT_Response.cxx.

{
  //double propDelay = PropagationDelay(m_xhit);
  double damping = DampingFactor(m_xhit);
  double r2(m_rhit*m_rhit);
  m_t0 = RtoT(m_rhit);
   
  for(int i=0;i<2;++i){        // splite track into to segments
    double cl(0.);             // current position along track
    while(cl < m_pathLength){   
      cl += DoStep(rndmEngine);    // Do step along track
      double r = sqrt(cl*cl + r2); // calculate corresponding r
      int q = GenerateQ(rndmEngine);// generate cluster size
      double t = RtoT(r);     
      // check if value lies within window
      if(t-m_t0 > m_bins*m_binsize + 3.*SigmaDiffusion(r) ) {
    break;
      }
      double tc = t + Diffusion(r, rndmEngine) - m_t0;
      int bin( (int)(tc/m_binsize+m_offset) );
    
      // fill value into corresponding bin
      if(bin < (int)m_clusters.size() && bin >=0 ){
          m_clusters[bin] += q*damping;
      }else if(bin < 0) {
          std::cout << "out of range " << tc << " bin " << bin << std::endl;
      }
    }
  }
}

DoStepping() [2/2]

void MDT_Response::DoStepping ( double ParticleCharge, double ParticleGamma, CLHEP::HepRandomEngine * rndmEngine )

private

Definition at line 183 of file MDT_Response.cxx.

{
// extract  Cluster Density for the given value of particle gamma from the vector gammaFactorVec()
    double correctedClusterDensity;
    int km=0;
    int kmm=0;
    if(ParticleGamma<0.90852E+06){
        while(ParticleGamma>m_gammaFactorVec.at(km)){
            ++km;
        }
        int mm=km-1;    
        if(km!=0){
            double deltagamma=fabs(m_gammaFactorVec.at(km)-m_gammaFactorVec.at(mm))/2.;
            kmm=km-1;
            if(ParticleGamma > (m_gammaFactorVec.at(mm)+deltagamma)){
                kmm=km;
            }
        }
    correctedClusterDensity = m_numberOfClustersPerMmVec.at(kmm);
    }else{
        correctedClusterDensity = m_clusterDensity; // response plateaus: default to value for 100 GeV muon
    }
    
  //double propDelay = PropagationDelay(m_xhit);
  double damping = DampingFactor(m_xhit);
  double r2(m_rhit*m_rhit);
  m_t0 = RtoT(m_rhit);
  
  for(int i=0;i<2;++i){        // splite track into to segments
    double cl(0.);             // current position along track
    while(cl < m_pathLength){ 
 
    if(fabs(ParticleCharge)!=1.){    
      cl += m_clusterDensity/(correctedClusterDensity*pow(ParticleCharge,2))*DoStep(rndmEngine);                // Do step along track
    }else{  
      cl += m_clusterDensity/(correctedClusterDensity)*DoStep(rndmEngine);                // Do step along track
    }  
      
      double r = sqrt(cl*cl + r2); // calculate corresponding r
      int q = GenerateQ(rndmEngine);         // generate cluster size
      double t = RtoT(r);     
      // check if value lies within window
      if(t-m_t0 > m_bins*m_binsize + 3.*SigmaDiffusion(r) ) {
    break;
      }
      double tc = t + Diffusion(r, rndmEngine) - m_t0;
      int bin( (int)(tc/m_binsize+m_offset) );
    
      // fill value into corresponding bin
      if(bin < (int)m_clusters.size() && bin >=0 ){
    m_clusters[bin] += q*damping;
      }else if(bin < 0) {
    std::cout << "out of range " << tc << " bin " << bin << std::endl;
      }
    }
  }
}

DriftTime()

double MDT_Response::DriftTime ( ) const

inline

Definition at line 179 of file MDT_Response.h.

{ return m_t0 + m_amplifier.TThreshold() - m_offset*m_binsize; }

GenerateQ()

int MDT_Response::GenerateQ ( CLHEP::HepRandomEngine * rndmEngine ) const

inline

Definition at line 120 of file MDT_Response.h.

{
  double v = CLHEP::RandFlat::shoot(rndmEngine);
  double p;
  if(v<0.08){
    p = CLHEP::RandPoisson::shoot(rndmEngine,13.);
  }else if(v<0.28){
    double v1 = CLHEP::RandFlat::shoot(rndmEngine);
    p = (1./(v1+1.e-7));
  }else{
    p = 1+CLHEP::RandPoisson::shoot(rndmEngine,0.05);
  }
  return (int)p;
}

GetSignal() [1/2]

bool MDT_Response::GetSignal

(

CLHEP::HepRandomEngine *

rndmEngine

)

Definition at line 253 of file MDT_Response.cxx.

{
  // generate and propagate clusters
  DoStepping(rndmEngine);
 
  // get amplifier response and see if it passed threshold
  bool hasSignal = m_amplifier.AddClusters(m_clusters);  
  return hasSignal;
}

GetSignal() [2/2]

bool MDT_Response::GetSignal	(	double	ParticleCharge,
		double	ParticleGamma,
		CLHEP::HepRandomEngine *	rndmEngine )

Definition at line 264 of file MDT_Response.cxx.

{
  // generate and propagate clusters
  DoStepping(ParticleCharge,ParticleGamma,rndmEngine); 
 
  // get amplifier response and see if it passed threshold
  bool hasSignal = m_amplifier.AddClusters(m_clusters);  
  return hasSignal;
}

InitClusters()

void MDT_Response::InitClusters ( double timewindow, double binsize )

private

Definition at line 37 of file MDT_Response.cxx.

{
  m_binsize = binsize;
  m_timeWindow = timewindow;
  m_offset = (int)(3*m_difSmearing/m_binsize);
  m_bins = (int)(m_timeWindow/m_binsize) + m_offset;
  m_clusters.resize(m_bins);
  m_amplifier.InitResponse(m_bins,binsize);
  //std::cout << "MDT_Response Initialize: new cluster vector with " << m_clusters.size() << " bins " << m_bins << std::endl;
}

InitdEdxTable()

void MDT_Response::InitdEdxTable ( )

private

Definition at line 64 of file MDT_Response.cxx.

{
 
    // initialize number of clusters per mm vs gamma table  
 
    double gammaFactor[] = {0.10110E+01,0.10133E+01,0.10162E+01,0.10196E+01,0.10237E+01,0.10287E+01,0.10348E+01,0.10422E+01,0.10511E+01,
    0.10619E+01,0.10750E+01,0.10909E+01,0.11101E+01,0.11334E+01,0.11616E+01,0.11957E+01,0.12371E+01,0.12873E+01,                       
    0.13481E+01,0.14217E+01,0.15109E+01,0.16190E+01,0.17499E+01,0.19085E+01,0.21007E+01,0.23335E+01,0.26156E+01,                       
    0.29573E+01,0.33714E+01,0.38730E+01,0.44807E+01,0.52170E+01,0.61090E+01,0.71897E+01,0.84989E+01,0.10085E+02,
    0.12007E+02,0.14335E+02,0.17156E+02,0.20573E+02,0.24714E+02,0.29730E+02,0.35807E+02,0.43170E+02,0.52090E+02,
    0.62897E+02,0.75989E+02,0.91852E+02,0.11107E+03,0.13435E+03,0.16256E+03,0.19673E+03,0.23814E+03,0.28830E+03,
    0.34907E+03,0.42270E+03,0.51190E+03,0.61997E+03,0.75089E+03,0.90952E+03,0.11017E+04,0.13345E+04,0.16166E+04,
    0.19583E+04,0.23724E+04,0.28740E+04,0.34817E+04,0.42180E+04,0.51100E+04,0.61907E+04,0.74999E+04,0.90862E+04,
    0.11008E+05,0.13336E+05,0.16157E+05,0.19574E+05,0.23715E+05,0.28731E+05,0.34808E+05,0.42171E+05,0.51091E+05,
    0.61898E+05,0.74990E+05,0.90853E+05,0.11007E+06,0.13335E+06,0.16156E+06,0.19574E+06,0.23714E+06,0.28730E+06,
    0.34807E+06,0.42170E+06,0.51090E+06,0.61897E+06,0.74990E+06,0.90852E+06};
 
    double numberOfClustersPerMm[] = {188.7333, 160.5230, 135.5958, 115.2599, 98.1318, 83.4168, 70.8906, 60.3481, 51.4957, 44.0192, 37.6769, 
  32.36505, 27.9281, 24.1556, 21.0378, 18.4179, 16.2422, 14.4153, 12.9175, 11.6626, 10.6435, 9.8294, 9.1667, 8.6551, 8.2926, 7.9775, 
  7.7862, 7.656, 7.5869, 7.5707, 7.6145, 7.6727, 7.7725, 7.8923, 8.0268, 8.2013, 8.3624, 8.5358, 8.7236, 8.9275, 9.1214, 9.2668, 9.3918, 9.5144, 
  9.6032, 9.6733, 9.7272, 9.7535, 9.7837, 9.8438, 9.8365, 9.8620, 9.8892, 9.9119, 9.9063, 9.9235, 9.9267, 9.9175, 9.9299, 9.9409, 9.9269, 9.9437, 
  9.9463, 9.9532, 9.9531, 9.9576, 9.9635, 9.9623, 9.9214, 9.9636, 9.9377, 9.9508, 9.9591, 9.9694, 9.937, 9.9606, 9.9611, 9.9481, 9.9492, 9.9469,
  9.9401, 9.9232, 9.9566, 9.9532, 9.9514, 9.9550, 9.9714, 9.9504, 9.9726, 9.9381, 9.9259, 9.9633, 9.9455, 9.9560, 9.9693, 9.9516};  
    
    for( int ik = 0; ik < 96; ++ik ) {
        m_gammaFactorVec.push_back( gammaFactor[ik] );
        m_numberOfClustersPerMmVec.push_back( numberOfClustersPerMm[ik] );
        }
}

InitRt()

void MDT_Response::InitRt ( )

private

Definition at line 95 of file MDT_Response.cxx.

                         {
  //std::cout << "Setting rt-relation" << std::endl;
  if(m_rtParameters) delete m_rtParameters;
  switch(m_rtMode){
  case 0:    /* linear rt */
    m_rtNpar = 1;
    m_rtParameters = new double[m_rtNpar];
    m_rtParameters[0] = 700./m_radius;
    break;
  case 1:    /* used for work presented in muonsw nov 2003 */
    m_rtNpar = 4;
    m_rtParameters = new double[m_rtNpar];
    m_rtParameters[0] = 0.;
    m_rtParameters[1] = 16.4;
    m_rtParameters[2] = 0.43;
    m_rtParameters[3] = 0.102;
    break;
  case 2: // garfield rt-relation
    m_rtNpar = 10;
    m_rtParameters = new double[m_rtNpar];
    m_rtParameters[0] = -1.79390e-01;
    m_rtParameters[1] = 2.04417e+01;
    m_rtParameters[2] = 4.96576e+00;
    m_rtParameters[3] = -2.59002e+00;
    m_rtParameters[4] = 5.45522e-01;
    m_rtParameters[5] = -4.66626e-02;
    m_rtParameters[6] = 5.66214e-04;
    m_rtParameters[7] = 1.89740e-04;
    m_rtParameters[8] = -1.30073e-05;
    m_rtParameters[9] = 2.68995e-07;
    break;
  case 3: // output rt-relation
    m_rtNpar = 10;
    m_rtParameters = new double[m_rtNpar];
    m_rtParameters[0] =  2.31708e+01;
    m_rtParameters[1] =  4.55828;
    m_rtParameters[2] =  1.41878e+01; 
    m_rtParameters[3] = -5.99392; 
    m_rtParameters[4] =  1.34392; 
    m_rtParameters[5] = -1.65248e-01; 
    m_rtParameters[6] =  1.15960e-02; 
    m_rtParameters[7] = -4.29922e-04; 
    m_rtParameters[8] =  6.14201e-06; 
    m_rtParameters[9] = 2.00429e-08; 
    break;
  default:
    std::cout << "MDT_Response:  Wrong rt-mode" << std::endl;
    m_rtNpar = 0;
    m_rtParameters = nullptr;
    break;
  }
//   std::cout << "MDT_Response: RT-relation -> Polynomial" << std::endl;
//   for(unsigned int i=0;i<m_rtNpar;++i)
//     std::cout << "    Par[" << i << "] = " << m_rtParameters[i] << std::endl;
}

InitTubeParameters()

void MDT_Response::InitTubeParameters ( )

private

Definition at line 48 of file MDT_Response.cxx.

{
  m_radius = 14.6275;        // tube radius in atlas
  m_clusterDensity = 10.;           // clusters per mm for 100 GeV muon
  m_attLength = 30000.;   // mm
  m_signalSpeed = 300.;   // mm/ns
  m_rtMode = 2;
  m_difSmearing = 10.;   
  m_integrationWindow = 18.5; // ns
  m_triggerElectron = 20.;
  m_amplifier.SetTriggerElectron(m_triggerElectron);
  m_amplifier.SetIntegrationWindow(m_integrationWindow);
  InitRt();
    
}

PropagationDelay()

double MDT_Response::PropagationDelay ( double x )

inline

Definition at line 139 of file MDT_Response.h.

{ return x/m_signalSpeed; }

Reset()

void MDT_Response::Reset ( )

private

Definition at line 245 of file MDT_Response.cxx.

{
  clusterVec::iterator cit = m_clusters.begin();
  for( ; cit != m_clusters.end() ; ++cit)
    *cit = 0;
  m_amplifier.Reset();
} 

RtoT()

double MDT_Response::RtoT ( double r )

inline

Definition at line 142 of file MDT_Response.h.

{
  int j = m_rtNpar-1;
  double time = m_rtParameters[j];
  while(j>0) { --j; time = time*r + m_rtParameters[j]; }
  return time;
} 

RtParameters()

const double * MDT_Response::RtParameters ( ) const

inline

Definition at line 227 of file MDT_Response.h.

{
  return m_rtParameters;
}

SetAttLength()

void MDT_Response::SetAttLength ( double len )

inline

Definition at line 204 of file MDT_Response.h.

{ m_attLength = len ; }

SetClusterDensity()

void MDT_Response::SetClusterDensity ( double dens )

inline

Definition at line 201 of file MDT_Response.h.

{ m_clusterDensity = dens ; }    

SetDifRMS()

void MDT_Response::SetDifRMS ( double rms )

inline

Definition at line 210 of file MDT_Response.h.

{ m_difSmearing = rms; InitClusters(m_timeWindow,m_binsize); }

SetIntegrationWindow()

void MDT_Response::SetIntegrationWindow ( double win )

inline

Definition at line 216 of file MDT_Response.h.

{ m_integrationWindow = win; m_amplifier.SetIntegrationWindow(win); }

SetRtMode()

void MDT_Response::SetRtMode ( unsigned int mode )

inline

Definition at line 207 of file MDT_Response.h.

{ m_rtMode = mode;   InitRt(); }

SetRtParameters()

void MDT_Response::SetRtParameters ( int npar, double * par )

inline

Definition at line 219 of file MDT_Response.h.

{
  m_rtNpar = npar;
  if(m_rtParameters) delete m_rtParameters;
  m_rtParameters = new double[m_rtNpar];
  for(unsigned int i=0;i<m_rtNpar;++i)
    m_rtParameters[i] = par[i]; 
}

SetSegment()

void MDT_Response::SetSegment ( double r, double x )

inline

Definition at line 105 of file MDT_Response.h.

{
  m_rhit = r;
  m_xhit = x;
  m_pathLength = sqrt(m_radius*m_radius - r*r);
  Reset();
}

SetTriggerElectron()

void MDT_Response::SetTriggerElectron ( double el )

inline

Definition at line 213 of file MDT_Response.h.

{  m_triggerElectron = el; m_amplifier.SetTriggerElectron(el); }

SetTubeRadius()

void MDT_Response::SetTubeRadius ( double radius )

inline

Definition at line 197 of file MDT_Response.h.

{ m_radius = radius ; }    

SigmaDiffusion()

double MDT_Response::SigmaDiffusion ( double r ) const

inline

Definition at line 173 of file MDT_Response.h.

{
  return m_difSmearing*r/m_radius;
}

T0()

double MDT_Response::T0 ( ) const

inline

Definition at line 188 of file MDT_Response.h.

{ return m_t0; }

TFirst()

double MDT_Response::TFirst ( ) const

inline

Definition at line 185 of file MDT_Response.h.

{ return m_amplifier.TFirst(); }

TThreshold()

double MDT_Response::TThreshold ( ) const

inline

Definition at line 182 of file MDT_Response.h.

{ return m_amplifier.TThreshold() - m_offset*m_binsize; }

V_r()

double MDT_Response::V_r ( double r )

inline

Definition at line 151 of file MDT_Response.h.

{
  int j = m_rtNpar-1;
  double p = m_rtParameters[j];
  double dp = 0.;
  while(j>0) {dp = dp*r + p;--j;p=p*r+m_rtParameters[j]; }
  return 1./dp;
} 

Member Data Documentation

m_amplifier

Amplifier MDT_Response::m_amplifier

private

Definition at line 98 of file MDT_Response.h.

m_attLength

double MDT_Response::m_attLength = 0.0

private

Definition at line 82 of file MDT_Response.h.

m_bins

int MDT_Response::m_bins = 0

private

Definition at line 96 of file MDT_Response.h.

m_binsize

double MDT_Response::m_binsize = 0.0

private

Definition at line 93 of file MDT_Response.h.

m_clusterDensity

double MDT_Response::m_clusterDensity = 0.0

private

Definition at line 78 of file MDT_Response.h.

m_clusters

clusterVec MDT_Response::m_clusters

private

Definition at line 99 of file MDT_Response.h.

m_difSmearing

double MDT_Response::m_difSmearing = 0.0

private

Definition at line 89 of file MDT_Response.h.

m_gammaFactorVec

std::vector<double> MDT_Response::m_gammaFactorVec

private

Definition at line 79 of file MDT_Response.h.

m_integrationWindow

double MDT_Response::m_integrationWindow = 0.0

private

Definition at line 92 of file MDT_Response.h.

m_numberOfClustersPerMmVec

std::vector<double> MDT_Response::m_numberOfClustersPerMmVec

private

Definition at line 80 of file MDT_Response.h.

m_offset

int MDT_Response::m_offset = 0

private

Definition at line 95 of file MDT_Response.h.

m_pathLength

double MDT_Response::m_pathLength

private

Definition at line 76 of file MDT_Response.h.

m_radius

double MDT_Response::m_radius = 0.0

private

Definition at line 72 of file MDT_Response.h.

m_rhit

double MDT_Response::m_rhit

private

Definition at line 74 of file MDT_Response.h.

m_rtMode

int MDT_Response::m_rtMode = 0

private

Definition at line 85 of file MDT_Response.h.

m_rtNpar

unsigned int MDT_Response::m_rtNpar = 0U

private

Definition at line 86 of file MDT_Response.h.

m_rtParameters

double* MDT_Response::m_rtParameters

private

Definition at line 87 of file MDT_Response.h.

m_signalSpeed

double MDT_Response::m_signalSpeed = 0.0

private

Definition at line 83 of file MDT_Response.h.

m_t0

double MDT_Response::m_t0

private

Definition at line 101 of file MDT_Response.h.

m_timeWindow

double MDT_Response::m_timeWindow = 0.0

private

Definition at line 94 of file MDT_Response.h.

m_triggerElectron

double MDT_Response::m_triggerElectron = 0.0

private

Definition at line 91 of file MDT_Response.h.

m_xhit

double MDT_Response::m_xhit

private

Definition at line 75 of file MDT_Response.h.

---

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

• MDT_Response.h
• MDT_Response.cxx