@@ -201,11 +201,8 @@ struct HfElectronSelectionWithTpcEmcal {
201201 registry.add (" hMcgenInElectron" , " Mc Gen Inclusive Electron" , {HistType::kTH1D , {{axisPt}}});
202202 registry.add (" hMcRecInElectron" , " Mc Rec Inclusive Electron" , {HistType::kTH1D , {{axisPt}}});
203203 registry.add (" hMcRecwithoutEMCalInElectron" , " Mc Rec Inclusive Electron without Emcal" , {HistType::kTH1D , {{axisPt}}});
204-
205204 registry.add (" hphiElectron" , " hphiElectron" , {HistType::kTH1D , {axisPhi}});
206205 registry.add (" hphiElectronPassEmcal" , " hphiElectron pass Emcal" , {HistType::kTH1D , {axisPhi}});
207-
208-
209206 registry.add (" hMcgenAllNonHfeElectron" , " Mc Gen All NonHf Electron" , {HistType::kTH1D , {{axisPt}}});
210207 registry.add (" hMcgenNonHfeElectron" , " Mc Gen NonHf Electron with mother" , {HistType::kTH1D , {{axisPt}}});
211208 registry.add (" hPi0eEmbTrkPt" , " Mc Gen Pi0 mother NonHf Electron" , {HistType::kTH1D , {{axisPt}}});
@@ -558,52 +555,6 @@ struct HfElectronSelectionWithTpcEmcal {
558555 bool isEmbEta = false ;
559556 bool isEmbPi0 = false ;
560557
561-
562- // Check first mother
563- if (mcParticle.has_mothers ()) {
564- auto const & mother = mcParticle.template mothers_first_as <aod::McParticles>();
565-
566- if (std::abs (mother.pdgCode ()) == kEtaLocal || std::abs (mother.pdgCode ()) == kPi0 || std::abs (mother.pdgCode ()) == kGamma ) {
567-
568- auto const & gmother = mother.template mothers_first_as <aod::McParticles>();
569- // cases to consider: eta->e, eta->pi0->e, eta->gamma->e, eta->pi0->gamma->e, pi0->e, pi0->gamma->e
570-
571- // ================= eta->e ======================================
572- if (std::abs (mother.pdgCode ()) == kEtaLocal ) {
573-
574- if (mother.isPhysicalPrimary ()) {
575- if ((std::abs (gmother.pdgCode ()) >= pdgCodeCharmMin && std::abs (gmother.pdgCode ()) < pdgCodeCharmMax) ||
576- (std::abs (gmother.pdgCode ()) >= pdgCodeBeautyMin && std::abs (gmother.pdgCode ()) < pdgCodeBeautyMax)) {
577- continue ;
578- }
579- isEmbEta = true ;
580- }
581- }
582-
583- // ================= eta->pi0->e ======================================
584-
585- if (std::abs (mother.pdgCode ()) == kPi0 ) {
586- if (mother.isPhysicalPrimary ()) {
587- if ((std::abs (gmother.pdgCode ()) >= pdgCodeCharmMin && std::abs (gmother.pdgCode ()) < pdgCodeCharmMax) ||
588- (std::abs (gmother.pdgCode ()) >= pdgCodeBeautyMin && std::abs (gmother.pdgCode ()) < pdgCodeBeautyMax)) {
589- continue ;
590- }
591- isEmbPi0 = true ; // pi0 -> e
592- }
593- if (std::abs (gmother.pdgCode ()) == kEtaLocal ) {
594- if (gmother.isPhysicalPrimary () || gmother.has_mothers ()) {
595- auto const & ggmother = gmother.template mothers_first_as <aod::McParticles>();
596- if ((std::abs (ggmother.pdgCode ()) >= pdgCodeCharmMin && std::abs (ggmother.pdgCode ()) < pdgCodeCharmMax) ||
597- (std::abs (ggmother.pdgCode ()) >= pdgCodeBeautyMin && std::abs (ggmother.pdgCode ()) < pdgCodeBeautyMax)) {
598- continue ;
599- }
600- isEmbEta = true ; // eta->pi0-> e
601- }
602- }
603- }
604-
605-
606-
607558 // Check first mother
608559 if (mcParticle.has_mothers ()) {
609560 auto const & mother = mcParticle.template mothers_first_as <aod::McParticles>();
@@ -647,7 +598,6 @@ struct HfElectronSelectionWithTpcEmcal {
647598 }
648599 }
649600
650- >
651601 // / ==================================== eta->gamma->e and eta->pi0->gamma->e============
652602 if (std::abs (mother.pdgCode ()) == kGamma ) {
653603
@@ -712,26 +662,12 @@ struct HfElectronSelectionWithTpcEmcal {
712662 if ((track.tpcNSigmaEl () < tpcNsigmaElectronMin || track.tpcNSigmaEl () > tpcNsigmaElectronMax)) {
713663 continue ;
714664 }
715-
716665 registry.fill (HIST (" hphiElectron" ), track.phi ());
717- =======
718-
719666 if constexpr (IsMc) {
720667 if (track.has_mcParticle ()) {
721668 auto mcParticle = track.template mcParticle_as <aod::McParticles>();
722669 if (std::abs (mcParticle.pdgCode ()) == kElectron ) {
723670
724-
725- registry.fill (HIST (" hMcRecwithoutEMCalInElectron" ), mcParticle.pt ());
726- bool isEmbEta = false ;
727- bool isEmbPi0 = false ;
728-
729- // Check first mother
730- if (mcParticle.has_mothers ()) {
731- auto const & mother = mcParticle.template mothers_first_as <aod::McParticles>();
732-
733-
734-
735671 registry.fill (HIST (" hMcRecwithoutEMCalInElectron" ), mcParticle.pt ());
736672 bool isEmbEta = false ;
737673 bool isEmbPi0 = false ;
0 commit comments