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. ;
293 ptcone20 = myptcone20;
296 Float_t myehad1 = 0.0;
297 Float_t myecore = 0.0;
305 Float_t myre233e237 = 0.0;
306 Float_t myre237e277 = 0.0;
315 if (e237!=0) myre233e237 = e233 / e237;
316 if (e277!=0) myre237e277 = e237 / e277;
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);
418 et = myet ;
eta = myeta ;
phi = myphi ; time = mytime ;
419 ptcone20 = myptcone20; topoetcone40 = mytopoetcone40 ;
421 ptcone20,topoetcone40,lbNProbeCandidates,is_TightID,is_Iso,is_IsoandTight);
423 return StatusCode::SUCCESS;