BesEmcWaveform Class Reference

#include <BesEmcWaveform.hh>

List of all members.

Public Member Functions

 BesEmcWaveform ()
 BesEmcWaveform (G4long, G4double, G4double)
virtual ~BesEmcWaveform ()
G4double & operator[] (G4long) const
virtual BesEmcWaveform & operator *= (const G4double &)
virtual BesEmcWaveform & operator/= (const G4double &)
virtual BesEmcWaveform & operator+= (const BesEmcWaveform &)
virtual BesEmcWaveform & operator= (const BesEmcWaveform &)
G4long length () const
G4double * GetWave () const
G4double GetTau () const
G4double GetSampleTime () const
G4double GetPeakTime () const
G4double GetTimeOffset () const
G4int GetBitNb () const
G4int GetGainFlag () const
G4double GetPhotonsPerMeV () const
G4double GetNonuniformity () const
G4double max (G4long &binOfMax) const
void updateWaveform (BesEmcHit *)
void updateWaveform (BesEmcDigi *)
void makeWaveform (G4double energy, G4double time)
void digitize ()
void addElecNoise (G4double, G4double)
void print ()

Protected Attributes

G4double * emcWave
G4long array_size

Private Attributes

G4double m_tau
G4double m_sampleTime
G4double m_peakTime
G4double m_timeOffset
G4int m_bitNb
G4int m_flag
G4double m_photonsPerMeV
G4double m_nonuniformity
G4double m_highRange
G4double m_highPrecision
G4double m_midRange
G4double m_midPrecision
G4double m_lowRange
G4double m_lowPrecision

Detailed Description

Definition at line 44 of file BesEmcWaveform.hh.

Constructor & Destructor Documentation

BesEmcWaveform::BesEmcWaveform (  ) 

Definition at line 21 of file BesEmcWaveform.cc.

References array_size, emcWave, BesEmcParameter::GetADCbit(), BesEmcParameter::GetArraySize(), BesEmcParameter::GetHighRange(), BesEmcParameter::GetInstance(), BesEmcParameter::GetLowRange(), BesEmcParameter::GetMidRange(), BesEmcParameter::GetNonuniformity(), BesEmcParameter::GetPeakTime(), BesEmcParameter::GetPhotonsPerMeV(), BesEmcParameter::GetSampleTime(), BesEmcParameter::GetTau(), BesEmcParameter::GetTimeOffset(), genRecEmupikp::i, m_bitNb, m_flag, m_highPrecision, m_highRange, m_lowPrecision, m_lowRange, m_midPrecision, m_midRange, m_nonuniformity, m_peakTime, m_photonsPerMeV, m_sampleTime, m_tau, and m_timeOffset.

00022 {
00023   BesEmcParameter& emcPara=BesEmcParameter::GetInstance();
00024   //BesEmcParameter emcPara;
00025   //emcPara.ReadData();
00026   array_size = emcPara.GetArraySize();
00027   m_tau = emcPara.GetTau();
00028   m_highRange = emcPara.GetHighRange();
00029   m_midRange = emcPara.GetMidRange();
00030   m_lowRange = emcPara.GetLowRange();
00031   m_sampleTime = emcPara.GetSampleTime();
00032   m_peakTime = emcPara.GetPeakTime();
00033   m_timeOffset = emcPara.GetTimeOffset();
00034   m_bitNb = emcPara.GetADCbit();
00035   m_photonsPerMeV = emcPara.GetPhotonsPerMeV();
00036   m_nonuniformity = emcPara.GetNonuniformity();
00037   m_flag = -1;
00038   m_highPrecision = m_highRange/((G4double)(1<<m_bitNb));
00039   m_midPrecision = m_midRange/((G4double)(1<<m_bitNb));
00040   m_lowPrecision = m_lowRange/((G4double)(1<<m_bitNb));
00041   emcWave=new G4double[array_size];
00042 
00043   for (G4long i=0; i<array_size;i++)
00044     emcWave[i]=0;
00045 }

BesEmcWaveform::BesEmcWaveform ( G4long  ,
G4double  ,
G4double   
)

Definition at line 47 of file BesEmcWaveform.cc.

References array_size, emcWave, and init.

00048   :m_tau(tau),m_sampleTime(sampleTime)
00049 {
00050   if(size>0){
00051     array_size=size;
00052     emcWave=new G4double[array_size];
00053     G4double *init=emcWave+array_size;
00054     while(init!=emcWave) *--init=0.0;
00055   }
00056   else{
00057     G4Exception("Invalid size");
00058   }
00059 }

BesEmcWaveform::~BesEmcWaveform (  )  [virtual]

Definition at line 62 of file BesEmcWaveform.cc.

References emcWave.

00062                                {
00063   delete []emcWave;
00064   emcWave=0;
00065 }

Member Function Documentation

void BesEmcWaveform::addElecNoise ( G4double  ,
G4double   
)

Definition at line 196 of file BesEmcWaveform.cc.

References array_size, emcWave, and genRecEmupikp::i.

Referenced by BesEmcDigitizer::AddNoise5x5(), BesEmcDigitizer::AddNoiseAll(), and BesEmcDigitizer::Digitize().

00197 {
00198   for(G4int i=0;i<array_size;i++)
00199   {
00200     emcWave[i] += RandGauss::shoot()*width;
00201     emcWave[i] += coherentNoise;
00202     emcWave[i] = emcWave[i]>0?emcWave[i]:0;
00203   }
00204 }

void BesEmcWaveform::digitize (  ) 

Definition at line 159 of file BesEmcWaveform.cc.

References array_size, emcWave, energy, genRecEmupikp::i, m_bitNb, m_flag, m_highPrecision, m_highRange, m_lowPrecision, m_lowRange, m_midPrecision, and m_midRange.

Referenced by BesEmcDigitizer::AddNoise5x5(), BesEmcDigitizer::AddNoiseAll(), and BesEmcDigitizer::Digitize().

00160 {
00161   G4double oneBitResolution;
00162   oneBitResolution = m_midRange*2/((G4double)(1<<m_bitNb));   //not used now
00163   G4double energy;
00164   
00165   for(G4long i=0;i<array_size;i++)
00166   {
00167     energy = emcWave[i];
00168     if(energy > m_highRange)
00169       G4cout<<"---In BesEmcWaveform: Over measurement!--- energy="<<energy<<G4endl;
00170     else if(energy > m_midRange)
00171     {
00172       m_flag = 2;
00173       emcWave[i] = (G4double)((G4long)(emcWave[i]/m_highPrecision))*m_highPrecision; 
00174     }
00175     else if(energy > m_lowRange)
00176     {
00177       m_flag = 1;
00178       emcWave[i] = (G4double)((G4long)(emcWave[i]/m_midPrecision))*m_midPrecision;
00179     }
00180     else
00181     {
00182       m_flag = 0;
00183       emcWave[i] = (G4double)((G4long)(emcWave[i]/m_lowPrecision))*m_lowPrecision;
00184     }
00185   }
00186 }

G4int BesEmcWaveform::GetBitNb (  )  const [inline]

Definition at line 69 of file BesEmcWaveform.hh.

References m_bitNb.

00069 { return m_bitNb; }

G4int BesEmcWaveform::GetGainFlag (  )  const [inline]

Definition at line 70 of file BesEmcWaveform.hh.

References m_flag.

00070 { return m_flag; }

G4double BesEmcWaveform::GetNonuniformity (  )  const [inline]

Definition at line 72 of file BesEmcWaveform.hh.

References m_nonuniformity.

00072 { return m_nonuniformity; }

G4double BesEmcWaveform::GetPeakTime (  )  const [inline]

Definition at line 67 of file BesEmcWaveform.hh.

References m_peakTime.

00067 { return m_peakTime; }

G4double BesEmcWaveform::GetPhotonsPerMeV (  )  const [inline]

Definition at line 71 of file BesEmcWaveform.hh.

References m_photonsPerMeV.

00071 { return m_photonsPerMeV; }

G4double BesEmcWaveform::GetSampleTime (  )  const [inline]

Definition at line 66 of file BesEmcWaveform.hh.

References m_sampleTime.

00066 { return m_sampleTime; }

G4double BesEmcWaveform::GetTau (  )  const [inline]

Definition at line 65 of file BesEmcWaveform.hh.

References m_tau.

00065 { return m_tau; }

G4double BesEmcWaveform::GetTimeOffset (  )  const [inline]

Definition at line 68 of file BesEmcWaveform.hh.

References m_timeOffset.

00068 { return m_timeOffset; }

G4double* BesEmcWaveform::GetWave (  )  const [inline]

Definition at line 64 of file BesEmcWaveform.hh.

References emcWave.

00064 { return emcWave; }

G4long BesEmcWaveform::length ( void   )  const [inline]

Definition at line 63 of file BesEmcWaveform.hh.

References array_size.

00063 {  return array_size;}

void BesEmcWaveform::makeWaveform ( G4double  energy,
G4double  time 
)

Definition at line 128 of file BesEmcWaveform.cc.

References array_size, emcWave, exp(), genRecEmupikp::i, m_peakTime, m_photonsPerMeV, m_sampleTime, m_tau, and m_timeOffset.

Referenced by updateWaveform().

00129 {
00130   G4double amplitude;
00131   if(m_photonsPerMeV==0)
00132   {
00133     amplitude = energy;
00134   }
00135   else
00136   {
00137     G4double photons = energy*m_photonsPerMeV;
00138     if(photons>0) {
00139       G4double photonStatFactor = RandGauss::shoot(photons, sqrt(photons))/photons;
00140       amplitude = energy*photonStatFactor;
00141     } else {
00142       amplitude = 0;
00143     }
00144   }
00145 
00146   G4double tempTime;
00147   tempTime = 0;
00148 
00149   G4double peak = m_peakTime*m_peakTime*m_peakTime*m_peakTime*exp(-m_peakTime/m_tau)/24;
00150 
00151   for (G4long i=0;i<array_size;i++)
00152   {
00153     tempTime = i*m_sampleTime + m_timeOffset - time;
00154     if(tempTime>0)
00155       emcWave[i] += amplitude*tempTime*tempTime*tempTime*tempTime*exp(-tempTime/m_tau)/(24*peak);
00156   }
00157 }

G4double BesEmcWaveform::max ( G4long &  binOfMax  )  const

Definition at line 97 of file BesEmcWaveform.cc.

References array_size, emcWave, and genRecEmupikp::i.

Referenced by BesEmcDigitizer::AddNoise5x5(), BesEmcDigitizer::AddNoiseAll(), BesEmcDigitizer::Digitize(), and BesEmcDigi::operator+=().

00098 {
00099   G4double maxi=emcWave[0];
00100   binOfMax = 0;
00101   for (G4long i=1;i<array_size;i++)
00102   {
00103     if (emcWave[i]>maxi) 
00104     {
00105       maxi=emcWave[i];
00106       binOfMax = i;
00107     }
00108   }
00109   return maxi;
00110 }

BesEmcWaveform & BesEmcWaveform::operator *= ( const G4double &   )  [virtual]

Definition at line 68 of file BesEmcWaveform.cc.

References array_size, emcWave, and genRecEmupikp::i.

00069 {
00070   for (G4long i=0; i<array_size;i++) emcWave[i]*=scale;
00071   return *this;
00072 }

BesEmcWaveform & BesEmcWaveform::operator+= ( const BesEmcWaveform &   )  [virtual]

Definition at line 80 of file BesEmcWaveform.cc.

References array_size, emcWave, and genRecEmupikp::i.

00081 {
00082   for (G4long i=0; i<array_size;i++) emcWave[i]+=assign[i];
00083   return *this;
00084 }

BesEmcWaveform & BesEmcWaveform::operator/= ( const G4double &   )  [virtual]

Definition at line 74 of file BesEmcWaveform.cc.

References array_size, emcWave, and genRecEmupikp::i.

00075 {
00076   for (G4long i=0; i<array_size;i++) emcWave[i]/=scale;
00077   return *this;
00078 }

BesEmcWaveform & BesEmcWaveform::operator= ( const BesEmcWaveform &   )  [virtual]

Definition at line 86 of file BesEmcWaveform.cc.

References array_size, emcWave, and genRecEmupikp::i.

00087 {
00088   if (this != &assign) {
00089     if (emcWave!=0) delete []emcWave;
00090     emcWave=new G4double[assign.array_size];
00091     array_size=assign.array_size;
00092     for (G4long i=0;i<array_size;i++) emcWave[i]=assign[i];
00093   }
00094   return *this;
00095 }

G4double & BesEmcWaveform::operator[] ( G4long   )  const [inline]

Definition at line 112 of file BesEmcWaveform.hh.

References array_size, and emcWave.

00113 {
00114   if(index<0 || index>array_size-1){
00115     G4cout << "Array bounds exceeded. Index " << index << G4endl;
00116     ::abort();
00117   }
00118   return emcWave[index];
00119 }

void BesEmcWaveform::print (  ) 

Definition at line 188 of file BesEmcWaveform.cc.

References array_size, emcWave, and genRecEmupikp::i.

00189 {
00190   G4cout<<"New Wave!"<<G4endl;
00191   for(G4long i=0;i<array_size;i++)
00192     G4cout<<emcWave[i]<<"\t";
00193   G4cout<<G4endl;
00194 }

void BesEmcWaveform::updateWaveform ( BesEmcDigi *   ) 

Definition at line 120 of file BesEmcWaveform.cc.

References energy, BesEmcDigi::GetEnergy(), BesEmcDigi::GetTime(), m_peakTime, m_sampleTime, m_timeOffset, and makeWaveform().

00121 {
00122   G4double energy = digi->GetEnergy();
00123   G4double time = digi->GetTime()*m_sampleTime+m_timeOffset-m_peakTime;
00124 
00125   makeWaveform(energy, time);
00126 }

void BesEmcWaveform::updateWaveform ( BesEmcHit *   ) 

Definition at line 112 of file BesEmcWaveform.cc.

References energy, BesEmcHit::GetEdepCrystal(), BesEmcHit::GetTimeCrystal(), and makeWaveform().

Referenced by BesEmcDigitizer::AddNoise5x5(), BesEmcDigitizer::AddNoiseAll(), BesEmcDigitizer::Digitize(), BesEmcDigi::MakeWaveform(), and BesEmcDigi::operator+=().

00113 {
00114   G4double energy = hit->GetEdepCrystal();
00115   G4double time = hit->GetTimeCrystal();
00116 
00117   makeWaveform(energy, time);
00118 }

Member Data Documentation

G4long BesEmcWaveform::array_size [protected]

Definition at line 88 of file BesEmcWaveform.hh.

Referenced by addElecNoise(), BesEmcWaveform(), digitize(), length(), makeWaveform(), max(), operator *=(), operator+=(), operator/=(), operator=(), operator[](), and print().

G4double* BesEmcWaveform::emcWave [protected]

Definition at line 87 of file BesEmcWaveform.hh.

Referenced by addElecNoise(), BesEmcWaveform(), digitize(), GetWave(), makeWaveform(), max(), operator *=(), operator+=(), operator/=(), operator=(), operator[](), print(), and ~BesEmcWaveform().

G4int BesEmcWaveform::m_bitNb [private]

Definition at line 96 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), digitize(), and GetBitNb().

G4int BesEmcWaveform::m_flag [private]

Definition at line 97 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), digitize(), and GetGainFlag().

G4double BesEmcWaveform::m_highPrecision [private]

Definition at line 104 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), and digitize().

G4double BesEmcWaveform::m_highRange [private]

Definition at line 103 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), and digitize().

G4double BesEmcWaveform::m_lowPrecision [private]

Definition at line 108 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), and digitize().

G4double BesEmcWaveform::m_lowRange [private]

Definition at line 107 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), and digitize().

G4double BesEmcWaveform::m_midPrecision [private]

Definition at line 106 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), and digitize().

G4double BesEmcWaveform::m_midRange [private]

Definition at line 105 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), and digitize().

G4double BesEmcWaveform::m_nonuniformity [private]

Definition at line 101 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), and GetNonuniformity().

G4double BesEmcWaveform::m_peakTime [private]

Definition at line 94 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), GetPeakTime(), makeWaveform(), and updateWaveform().

G4double BesEmcWaveform::m_photonsPerMeV [private]

Definition at line 100 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), GetPhotonsPerMeV(), and makeWaveform().

G4double BesEmcWaveform::m_sampleTime [private]

Definition at line 93 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), GetSampleTime(), makeWaveform(), and updateWaveform().

G4double BesEmcWaveform::m_tau [private]

Definition at line 92 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), GetTau(), and makeWaveform().

G4double BesEmcWaveform::m_timeOffset [private]

Definition at line 95 of file BesEmcWaveform.hh.

Referenced by BesEmcWaveform(), GetTimeOffset(), makeWaveform(), and updateWaveform().

Generated on Tue Nov 29 23:17:30 2016 for BOSS_7.0.2 by  doxygen 1.4.7