- Timestamp:
- Nov 17, 2008, 8:49:25 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SmearUtil.cc
r31 r32 122 122 else if(strstr(temp_string.c_str(),"ELG_Cfwd")){curstring >> varname >> value; ELG_Cfwd = value;} 123 123 else if(strstr(temp_string.c_str(),"ELG_Nfwd")){curstring >> varname >> value; ELG_Nfwd = value;} 124 else if(strstr(temp_string.c_str(),"HAD_Secal")){curstring >> varname >> value; HAD_Secal = value;}125 else if(strstr(temp_string.c_str(),"HAD_Necal")){curstring >> varname >> value; HAD_Necal = value;}126 else if(strstr(temp_string.c_str(),"HAD_Cecal")){curstring >> varname >> value; HAD_Cecal = value;}127 124 else if(strstr(temp_string.c_str(),"HAD_Shcal")){curstring >> varname >> value; HAD_Shcal = value;} 128 125 else if(strstr(temp_string.c_str(),"HAD_Nhcal")){curstring >> varname >> value; HAD_Nhcal = value;} … … 227 224 228 225 energyS2 = gRandom->Gaus(energy_ecal, sqrt( 229 pow( HAD_Necal,2) +230 pow( HAD_Cecal*energy_ecal,2) +231 pow( HAD_Secal*sqrt(energy_ecal),2) ) );226 pow(ELG_Ncen,2) + 227 pow(ELG_Ccen*energy_ecal,2) + 228 pow(ELG_Scen*sqrt(energy_ecal),2) ) ); 232 229 233 230 energyS = ((energyS1>0)?energyS1:0) + ((energyS2>0)?energyS2:0);
Note:
See TracChangeset
for help on using the changeset viewer.