19 return StatusCode::SUCCESS;
30 return StatusCode::SUCCESS;
50 for (
const auto *
const e_iter : *
electrons) {
53 if (! e_iter->passSelection(isGood,
m_RecoName)) {
57 if(isGood) {mynp++;}
else continue;
60 float etloc = e_iter->pt();
66 }
else if(etloc>subl_et) {
73 if (mynp<2)
return StatusCode::SUCCESS;
82 float lead_eta = lead_el->
eta();
83 float lead_phi = lead_el->
phi();
84 float subl_eta = subl_el->
eta();
85 float subl_phi = subl_el->
phi();
86 ATH_MSG_DEBUG(
"Leading electron (eta,phi,et,q): (" << lead_eta <<
", " << lead_phi <<
", " << lead_et <<
", " << lead_el->
charge() <<
")");
87 ATH_MSG_DEBUG(
"Subleading electron (eta,phi,et,q): (" << subl_eta <<
", " << subl_phi <<
", " << subl_et <<
", " << subl_el->
charge() <<
")");
90 if(fabs(lead_eta)>2.47 || fabs(subl_eta)>2.47)
return StatusCode::SUCCESS;
92 if (lead_el->
charge()*subl_el->
charge()>=0)
return StatusCode::SUCCESS;
95 Float_t
mass = (lead_el->
p4()+subl_el->
p4()).M();
104 fill(
"MonitorTnP",TnPMass);
109 TnPMass_barrel =
mass;
110 fill(
"MonitorTnP", TnPMass_barrel);
113 TnPMass_endcap =
mass;
114 fill(
"MonitorTnP", TnPMass_endcap);
117 TnPMass_crack =
mass;
118 fill(
"MonitorTnP", TnPMass_crack);
126 if(mass<m_MassLowerCut || mass>
m_MassUpperCut)
return StatusCode::SUCCESS;
133 fill(
"MonitorTnP",lbNZ);
135 bool lead_isLHTight =
false;
136 if ( !lead_el->
passSelection(lead_isLHTight,
"LHTight") )
return StatusCode::FAILURE;
138 bool subl_isLHTight =
false;
139 if ( !subl_el->
passSelection(subl_isLHTight,
"LHTight") )
return StatusCode::FAILURE;
150 return StatusCode::SUCCESS;
166 fill(
"MonitorTnP",EtaZ,PhiZ);
171 bool lead_isLHTight =
false;
172 if ( !
el->passSelection(lead_isLHTight,
"LHTight") )
return StatusCode::FAILURE;
173 is_TightID = lead_isLHTight;
178 Float_t myet =
el->pt();
179 Float_t myptcone20 = -999.;
180 Float_t mytopoetcone40 = -999.;
182 bool isIsolated =
false;
187 if ( myptcone20/myet < 0.15 ) isIsolated =
true;
191 is_IsoandTight = isIsolated && lead_isLHTight;
275 Float_t myeta =
el->eta();
276 Float_t myphi =
el->phi();
279 is_pt_gt_4gev = myet > 4000. ;
282 is_pt_gt_4gev_barrel = myet > 4000. ;
285 is_pt_gt_4gev_endcap = myet > 4000. ;
288 is_pt_gt_4gev_crack = myet > 4000. ;
296 Float_t myehad1 = 0.0;
297 Float_t myecore = 0.0;
305 Float_t myre233e237 = 0.0;
306 Float_t myre237e277 = 0.0;
319 Float_t mydeltaeta1 = -999.0;
321 Float_t mydeltaphi2 = -999.0;
324 Float_t myeoverp = -999.0;
325 Float_t mytrackd0 = -999.0;
330 u_int8_t mynofblayerhits=-1;
331 u_int8_t mynofpixelhits=-1;
332 u_int8_t mynofscthits=-1;
333 u_int8_t mynoftrthits=-1;
334 u_int8_t mynoftrthighthresholdhits=-1;
337 trackp =
t->pt()*cosh(
t->eta());
347 ATH_MSG_DEBUG(
"Error attempting to retrieve associated track");
354 mytime = aCluster->
time();
357 double ec = aCluster->
et()*cosh(aCluster->
eta());
358 if (ec!=0) myf0 = aCluster->
energyBE(0)/ec;
359 if (ec!=0) myf1 = aCluster->
energyBE(1)/ec;
360 if (ec!=0) myf2 = aCluster->
energyBE(2)/ec;
361 if (ec!=0) myf3 = aCluster->
energyBE(3)/ec;
363 if(trackp !=0) myeoverp = aCluster->
e()/trackp;
373 nofblayerhits_barrel = mynofblayerhits ; nofpixelhits_barrel = mynofpixelhits ; nofscthits_barrel = mynofscthits ;
374 noftrthits_barrel = mynoftrthits ; noftrthighthresholdhits_barrel = mynoftrthighthresholdhits;
375 deltaeta1_barrel = mydeltaeta1; deltaphi2_barrel = mydeltaphi2; trackd0_barrel = mytrackd0;
376 eoverp_barrel = myeoverp;
377 et_barrel = myet ; eta_barrel = myeta ; phi_barrel = myphi ;
378 time_barrel = mytime; ehad1_barrel = myehad1; coreem_barrel = myecore;
379 f0_barrel = myf0; f1_barrel = myf1; f2_barrel = myf2; f3_barrel = myf3; re233e237_barrel = myre233e237; re237e277_barrel = myre237e277;
380 fill(
"MonitorTnP",nofblayerhits_barrel, nofpixelhits_barrel, nofscthits_barrel, noftrthits_barrel, noftrthighthresholdhits_barrel,
381 deltaeta1_barrel, deltaphi2_barrel,trackd0_barrel,eoverp_barrel,
382 et_barrel,eta_barrel,phi_barrel, time_barrel, ehad1_barrel,coreem_barrel, f0_barrel,f1_barrel,f2_barrel, f3_barrel,
383 re233e237_barrel,re237e277_barrel,is_pt_gt_4gev_barrel);
388 nofblayerhits_endcap = mynofblayerhits; nofpixelhits_endcap = mynofpixelhits ; nofscthits_endcap = mynofscthits ;
389 noftrthits_endcap = mynoftrthits ; noftrthighthresholdhits_endcap = mynoftrthighthresholdhits;
390 deltaeta1_endcap = mydeltaeta1; deltaphi2_endcap = mydeltaphi2; trackd0_endcap = mytrackd0;
391 eoverp_endcap = myeoverp;
392 et_endcap = myet ; eta_endcap = myeta ; phi_endcap = myphi ;
393 time_endcap = mytime; ehad1_endcap = myehad1; coreem_endcap = myecore;
394 f0_endcap = myf0; f1_endcap = myf1; f2_endcap = myf2; f3_endcap = myf3; re233e237_endcap = myre233e237; re237e277_endcap = myre237e277;
395 fill(
"MonitorTnP",nofblayerhits_endcap,nofpixelhits_endcap,nofscthits_endcap,noftrthits_endcap,noftrthighthresholdhits_endcap,
396 deltaeta1_endcap, deltaphi2_endcap,trackd0_endcap,eoverp_endcap,et_endcap,eta_endcap,phi_endcap,
397 time_endcap, ehad1_endcap,coreem_endcap,
398 f0_endcap,f1_endcap,f2_endcap,f3_endcap,re233e237_endcap,re237e277_endcap,is_pt_gt_4gev_endcap);
402 nofblayerhits_crack = mynofblayerhits; nofpixelhits_crack = mynofpixelhits ; nofscthits_crack = mynofscthits ;
403 noftrthits_crack = mynoftrthits ; noftrthighthresholdhits_crack = mynoftrthighthresholdhits;
404 deltaeta1_crack = mydeltaeta1; deltaphi2_crack = mydeltaphi2; trackd0_crack = mytrackd0;
405 eoverp_crack = myeoverp;
406 et_crack = myet ; eta_crack = myeta ; phi_crack = myphi ;
407 time_crack = mytime; ehad1_crack = myehad1; coreem_crack = myecore;
408 f0_crack = myf0; f1_crack = myf1; f2_crack = myf2; f3_crack = myf3; re233e237_crack = myre233e237; re237e277_crack = myre237e277;
409 fill(
"MonitorTnP",nofblayerhits_crack, nofpixelhits_crack,nofscthits_crack,noftrthits_crack,noftrthighthresholdhits_crack,
410 deltaeta1_crack, deltaphi2_crack,trackd0_crack,eoverp_crack,et_crack,eta_crack,phi_crack,
411 time_crack, ehad1_crack,coreem_crack,f0_crack,f1_crack,f2_crack,f3_crack,re233e237_crack,re237e277_crack,is_pt_gt_4gev_crack);
423 return StatusCode::SUCCESS;