Changeset 69 in svn for trunk/src/SmearUtil.cc
- Timestamp:
- Dec 2, 2008, 12:17:27 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SmearUtil.cc
r62 r69 52 52 ELG_Ncen = 0.25 ; // N term for central ECAL 53 53 ELG_Ccen = 0.0055 ; // C term for central ECAL 54 ELG_ Cfwd = 0.107 ;// S term for forward ECAL55 ELG_ Sfwd = 2.084 ; // C term for forwardECAL56 ELG_ Nfwd = 0.0 ; // N term for centralECAL54 ELG_Sfwd = 1.5; // S term for forward ECAL 55 ELG_Nfwd = 0.0; // N term for central ECAL 56 ELG_Cfwd = 0.06; // C term for forward ECAL 57 57 58 58 HAD_Shcal = 1.5 ; // S term for central HCAL // hadronic calorimeter … … 182 182 PROTOJET_PTMIN = 0.0; 183 183 184 DOTRIGGER=1; 184 185 185 186 } 186 187 187 188 void RESOLution::Logfile(string LogName) { 188 //void RESOLution::Logfile(string outputfilename) {189 189 190 190 ofstream f_out(LogName.c_str()); … … 446 446 float energy = electron.E(); // before smearing 447 447 float energyS = 0.0; // after smearing // \sigma/E = C + N/E + S/\sqrt{E} 448 449 if(fabs(electron.Eta()) < MAX_TRACKER) { // if the electron is inside the tracker 448 float eta=fabs(electron.Eta()); 449 450 if(eta < MAX_TRACKER) { // if the electron is inside the tracker 450 451 energyS = gRandom->Gaus(energy, sqrt( 451 452 pow(ELG_Ncen,2) + … … 453 454 pow(ELG_Scen*sqrt(energy),2) )); 454 455 } 455 if( fabs(electron.Eta()) > MAX_TRACKER && fabs(electron.Eta())< MAX_CALO_FWD){456 if(eta > MAX_TRACKER && eta < MAX_CALO_FWD){ 456 457 energyS = gRandom->Gaus(energy, sqrt( 457 458 pow(ELG_Nfwd,2) +
Note:
See TracChangeset
for help on using the changeset viewer.