1 #ifndef STANDALONE_ECALCORR 14 if (runNumber <= 203755 && runNumber > 197770)
16 else if (runNumber > 203755)
28 if (runRange ==
RunAB) {
30 result = (iso-0.2827)/1.0949;
32 result = (iso-0.0931)/1.0738;
33 }
else if (runRange ==
RunC) {
35 result = (iso-0.5690)/0.9217;
37 result = (iso-0.1824)/0.9279;
38 }
else if (runRange ==
RunD) {
40 result = (iso-0.9997)/0.8781;
42 result = (iso-0.0944)/0.8140;
45 std::cout <<
"Warning: you should correct MC to data" << std::endl;
49 if (runRange ==
RunAB) {
51 result = (iso+0.1174)/1.0012;
53 result = (iso+0.2736)/0.9948;
54 }
else if (runRange ==
RunC) {
56 result = (iso+0.2271)/0.9684;
58 result = (iso+0.5962)/0.9568;
59 }
else if (runRange ==
RunD) {
61 result = (iso+0.2907)/1.0005;
63 result = (iso+0.9098)/0.9395;
66 std::cout <<
"Warning: you should correct MC to data" << std::endl;
73 #ifndef STANDALONE_ECALCORR 78 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
89 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
94 if (runName ==
"RunAB")
96 else if (runName ==
"RunC")
98 else if (runName ==
"RunD")
101 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
111 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
125 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
136 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
141 if (runName ==
"RunAB")
143 else if (runName ==
"RunC")
145 else if (runName ==
"RunD")
148 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
158 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
172 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
183 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
188 if (runName ==
"RunAB")
190 else if (runName ==
"RunC")
192 else if (runName ==
"RunD")
195 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
205 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
219 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
230 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
235 if (runName ==
"RunAB")
237 else if (runName ==
"RunC")
239 else if (runName ==
"RunD")
242 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
252 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
273 if (runName ==
"RunAB")
275 else if (runName ==
"RunC")
277 else if (runName ==
"RunD")
280 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
302 if (runRange ==
RunAB) {
304 result = iso*0.8499-0.6510;
306 result = iso*0.8504-0.5658;
307 }
else if (runRange ==
RunC) {
309 result = iso*0.9346-0.9987;
311 result = iso*0.8529-0.6816;
312 }
else if (runRange ==
RunD) {
314 result = iso*0.8318-0.9999;
316 result = iso*0.8853-0.8783;
319 if (runRange ==
RunAB) {
321 result = iso*0.9849-0.6871;
323 result = iso*0.8542-0.3558;
324 }
else if (runRange ==
RunC) {
326 result = iso*0.9996-0.8485;
328 result = iso*0.9994-0.5085;
329 }
else if (runRange ==
RunD) {
331 result = iso*0.9467-0.9998;
333 result = iso*0.8574-0.4862;
340 #ifndef STANDALONE_ECALCORR 345 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
359 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
363 if (runName ==
"RunAB")
365 else if (runName ==
"RunC")
367 else if (runName ==
"RunD")
370 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
383 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
400 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
414 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
418 if (runName ==
"RunAB")
420 else if (runName ==
"RunC")
422 else if (runName ==
"RunD")
425 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
438 std::cerr <<
"Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
455 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
469 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
474 if (runName ==
"RunAB")
476 else if (runName ==
"RunC")
478 else if (runName ==
"RunD")
481 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
494 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
512 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
526 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
531 if (runName ==
"RunAB")
533 else if (runName ==
"RunC")
535 else if (runName ==
"RunD")
538 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
551 std::cerr <<
"Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
578 if (runName ==
"RunAB")
580 else if (runName ==
"RunC")
582 else if (runName ==
"RunD")
585 std::cerr <<
"Error: Unknown run range " << runName << std::endl;
Analysis-level Photon class.
bool isBarrel(GeomDetEnumerators::SubDetector m)
float correctForNoise(reco::GsfElectron e, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
RunRange checkRunRange(int runNumber)
float correctForHLTDefinition(reco::GsfElectron e, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
Analysis-level electron class.
float dr03EcalRecHitSumEt() const
float ecalRecHitSumEtConeDR03() const
Isolation variables in cone dR=0.3.
EcalIsolationCorrector(bool forElectrons)