Main Page | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages

BesEmcWaveform Class Reference

#include <BesEmcWaveform.hh>

List of all members.

Public Member Functions

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

Protected Attributes

G4long array_size
G4double * emcWave
G4double * emcWave

Private Attributes

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

Constructor & Destructor Documentation

BesEmcWaveform::BesEmcWaveform (   ) 
  

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 
) 
  

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]
  

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

BesEmcWaveform::BesEmcWaveform (   ) 
 

BesEmcWaveform::BesEmcWaveform (  G4long  ,
G4double  ,
G4double 
) 
 

virtual BesEmcWaveform::~BesEmcWaveform (   )  [virtual]
 

Member Function Documentation

void BesEmcWaveform::addElecNoise (  G4double  ,
G4double 
) 
 

void BesEmcWaveform::addElecNoise (  G4double  ,
G4double 
) 
  

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 (   ) 
 

void BesEmcWaveform::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]
  

00069 { return m_bitNb; }

G4int BesEmcWaveform::GetBitNb (   )  const [inline]
  

00069 { return m_bitNb; }

G4int BesEmcWaveform::GetGainFlag (   )  const [inline]
  

00070 { return m_flag; }

G4int BesEmcWaveform::GetGainFlag (   )  const [inline]
  

00070 { return m_flag; }

G4double BesEmcWaveform::GetNonuniformity (   )  const [inline]
  

00072 { return m_nonuniformity; }

G4double BesEmcWaveform::GetNonuniformity (   )  const [inline]
  

00072 { return m_nonuniformity; }

G4double BesEmcWaveform::GetPeakTime (   )  const [inline]
  

00067 { return m_peakTime; }

G4double BesEmcWaveform::GetPeakTime (   )  const [inline]
  

00067 { return m_peakTime; }

G4double BesEmcWaveform::GetPhotonsPerMeV (   )  const [inline]
  

00071 { return m_photonsPerMeV; }

G4double BesEmcWaveform::GetPhotonsPerMeV (   )  const [inline]
  

00071 { return m_photonsPerMeV; }

G4double BesEmcWaveform::GetSampleTime (   )  const [inline]
  

00066 { return m_sampleTime; }

G4double BesEmcWaveform::GetSampleTime (   )  const [inline]
  

00066 { return m_sampleTime; }

G4double BesEmcWaveform::GetTau (   )  const [inline]
  

00065 { return m_tau; }

G4double BesEmcWaveform::GetTau (   )  const [inline]
  

00065 { return m_tau; }

G4double BesEmcWaveform::GetTimeOffset (   )  const [inline]
  

00068 { return m_timeOffset; }

G4double BesEmcWaveform::GetTimeOffset (   )  const [inline]
  

00068 { return m_timeOffset; }

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

00064 { return emcWave; }

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

00064 { return emcWave; }

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

00063 {  return array_size;}

G4long BesEmcWaveform::length (   )  const [inline]
  

00063 {  return array_size;}

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

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

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
 

G4double BesEmcWaveform::max (  G4long &  binOfMax  )  const
  

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 }

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

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

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

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

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

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

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

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

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

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

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

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]
 

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

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 (   ) 
 

void BesEmcWaveform::print (   ) 
  

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 *   ) 
 

void BesEmcWaveform::updateWaveform (  BesEmcHit *   ) 
 

void BesEmcWaveform::updateWaveform (  BesEmcDigi *   ) 
  

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 *   ) 
  

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]
 

G4double* BesEmcWaveform::emcWave [protected]
 

G4double* BesEmcWaveform::emcWave [protected]
 

G4int BesEmcWaveform::m_bitNb [private]
 

G4int BesEmcWaveform::m_flag [private]
 

G4double BesEmcWaveform::m_highPrecision [private]
 

G4double BesEmcWaveform::m_highRange [private]
 

G4double BesEmcWaveform::m_lowPrecision [private]
 

G4double BesEmcWaveform::m_lowRange [private]
 

G4double BesEmcWaveform::m_midPrecision [private]
 

G4double BesEmcWaveform::m_midRange [private]
 

G4double BesEmcWaveform::m_nonuniformity [private]
 

G4double BesEmcWaveform::m_peakTime [private]
 

G4double BesEmcWaveform::m_photonsPerMeV [private]
 

G4double BesEmcWaveform::m_sampleTime [private]
 

G4double BesEmcWaveform::m_tau [private]
 

G4double BesEmcWaveform::m_timeOffset [private]
 

The documentation for this class was generated from the following files:
Generated on Wed Feb 2 15:52:24 2011 for BOSS6.5.5 by  doxygen 1.3.9.1