#include <EmcBhabha.h>
Public Member Functions | |
const double & | calibEnergy () const |
const double & | calibEnergy () const |
EmcBhabha (const EmcBhabha &otherBhabha) | |
EmcBhabha () | |
EmcBhabha (const EmcBhabha &otherBhabha) | |
EmcBhabha () | |
const double & | errorOnCalibEnergy () const |
const double & | errorOnCalibEnergy () const |
double | errorOnShowerEnergy () const |
double | errorOnShowerEnergy () const |
const bool & | found () |
const bool & | found () |
const double & | phi () const |
const double & | phi () const |
const unsigned int & | phiIndex () const |
const unsigned int & | phiIndex () const |
void | print () |
void | print () |
void | setCalibEnergy (double energy) |
void | setCalibEnergy (double energy) |
void | setErrorOnCalibEnergy (double error) |
void | setErrorOnCalibEnergy (double error) |
void | setFound (bool what) |
void | setFound (bool what) |
void | setPhi (double phi) |
void | setPhi (double phi) |
void | setPhiIndex (unsigned int phiIndex) |
void | setPhiIndex (unsigned int phiIndex) |
void | setShower (EmcShower aShower) |
void | setShower (EmcShower aShower) |
void | setTheta (double theta) |
void | setTheta (double theta) |
void | setThetaIndex (unsigned int thetaIndex) |
void | setThetaIndex (unsigned int thetaIndex) |
EmcShower | shower () const |
EmcShower | shower () const |
double | sigma2 () const |
double | sigma2 () const |
const double & | theta () const |
const double & | theta () const |
const unsigned int & | thetaIndex () const |
const unsigned int & | thetaIndex () const |
~EmcBhabha () | |
~EmcBhabha () | |
Private Attributes | |
double | m_calibEnergy |
double | m_errCalibEnergy |
bool | m_found |
double | m_phi |
unsigned int | m_phiIndex |
EmcShower | m_shower |
double | m_theta |
unsigned int | m_thetaIndex |
|
00031 :m_calibEnergy(0), 00032 m_errCalibEnergy(0), 00033 m_shower(EmcShower::EmcShower()), 00034 m_theta(-1.), 00035 m_found(false), 00036 m_thetaIndex(9999), 00037 m_phiIndex(9999) 00038 { 00039 00040 }
|
|
00043 { 00044 m_calibEnergy = otherBhabha.m_calibEnergy; 00045 m_errCalibEnergy = otherBhabha.m_errCalibEnergy; 00046 m_shower = otherBhabha.m_shower; 00047 m_theta = otherBhabha.m_theta; 00048 m_found = otherBhabha.m_found; 00049 m_thetaIndex = otherBhabha.m_thetaIndex; 00050 m_phiIndex = otherBhabha.m_phiIndex; 00051 }
|
|
00056 { 00057 00058 }
|
|
|
|
|
|
|
|
00046 {return m_calibEnergy;}
|
|
00046 {return m_calibEnergy;}
|
|
00049 {return m_errCalibEnergy;}
|
|
00049 {return m_errCalibEnergy;}
|
|
|
|
00061 { 00062 00063 double sig=0,sig1=0,sig2=0; 00064 00065 //energy measurement error of calorimeter 00066 /* 00067 sig1 = 0.0034/shower().energy(); 00068 sig1 *= sig1; 00069 sig2 = (0.0070/pow(shower().energy(),0.25)); 00070 sig2 *= sig2; 00071 sig = sig1+sig2; 00072 sig += (0.0134*0.0134); 00073 sig = sqrt(sig) * shower().energy(); 00074 */ 00075 return sig; 00076 }
|
|
00043 {return m_found;}
|
|
00043 {return m_found;}
|
|
00069 {
00070 return m_phi;
00071 }
|
|
00069 {
00070 return m_phi;
00071 }
|
|
00079 {
00080 return m_phiIndex;
00081 }
|
|
00079 {
00080 return m_phiIndex;
00081 }
|
|
|
|
00095 { 00096 00097 cout << "EmcBhabha: " << endl 00098 << "Energies: calib: " << calibEnergy() 00099 << " shower: " << m_shower.energy() 00100 << " Theta: " << theta() 00101 << " nr Xtals: " << m_shower.digiList().size() 00102 << endl 00103 << "Errors: calibEne: " << errorOnCalibEnergy() 00104 << " shower: " << errorOnShowerEnergy() 00105 << " sum: " << sigma2() 00106 << endl; 00107 00108 }
|
|
00090 {m_calibEnergy=energy;}
|
|
00090 {m_calibEnergy=energy;}
|
|
00091 { m_errCalibEnergy = error;}
|
|
00091 { m_errCalibEnergy = error;}
|
|
00089 {m_found=what;}
|
|
00089 {m_found=what;}
|
|
00094 {m_phi=phi;}
|
|
00094 {m_phi=phi;}
|
|
00096 {m_phiIndex=phiIndex;}
|
|
00096 {m_phiIndex=phiIndex;}
|
|
00092 { m_shower=aShower;}
|
|
00092 { m_shower=aShower;}
|
|
00093 {m_theta=theta;}
|
|
00093 {m_theta=theta;}
|
|
00095 {m_thetaIndex=thetaIndex;}
|
|
00095 {m_thetaIndex=thetaIndex;}
|
|
00055 {return m_shower;}
|
|
00055 {return m_shower;}
|
|
|
|
00079 { 00080 00081 double sigSh=0., sigDepo=0., sig2=0.; 00082 00083 //error on measured shower energy 00084 sigSh = errorOnShowerEnergy(); 00085 00086 //error of energy to calibrate on 00087 sigDepo = errorOnCalibEnergy(); 00088 00089 sig2 = (sigSh * sigSh) + (sigDepo * sigDepo); 00090 return sig2; 00091 }
|
|
00060 { 00061 if ( m_theta < 0 ) 00062 { 00063 cout<<"warning: " << " EmcBhabha::theta not yet set !" << endl; 00064 } 00065 return m_theta; 00066 }
|
|
00060 { 00061 if ( m_theta < 0 ) 00062 { 00063 cout<<"warning: " << " EmcBhabha::theta not yet set !" << endl; 00064 } 00065 return m_theta; 00066 }
|
|
00074 {
00075 return m_thetaIndex;
00076 }
|
|
00074 {
00075 return m_thetaIndex;
00076 }
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|