@@ -447,7 +447,7 @@ struct HfCorrelatorLcHadrons {
447447 }
448448 }
449449 if (cntLc == 0 ) {
450- entryHadron (track.phi (), track.eta (), track.pt ()* track.sign (), poolBin, gCollisionId , timeStamp);
450+ entryHadron (track.phi (), track.eta (), track.pt () * track.sign (), poolBin, gCollisionId , timeStamp);
451451 if (fillTrkPID) {
452452 entryTrkPID (track.tpcNSigmaPr (), track.tpcNSigmaKa (), track.tpcNSigmaPi (), track.tofNSigmaPr (), track.tofNSigmaKa (), track.tofNSigmaPi ());
453453 }
@@ -585,30 +585,30 @@ struct HfCorrelatorLcHadrons {
585585 }
586586 }
587587 registry.fill (HIST (" hLcBin" ), poolBin);
588-
589- if (calTrkEff && !isLcSignal && calEffLcEvent) continue ;
590588
591- if (calTrkEff && cntLc==1 ){
592- // genrated tracks
593- for (const auto & track : mcParticles) {
594- if (std::abs (track.eta ()) > etaTrackMax || track.pt () < ptTrackMin || track.pt () > ptTrackMax) {
595- continue ;
596- }
597- if ((std::abs (track.pdgCode ()) != kElectron ) && (std::abs (track.pdgCode ()) != kMuonMinus ) && (std::abs (track.pdgCode ()) != kPiPlus ) && (std::abs (track.pdgCode ()) != kKPlus ) && (std::abs (track.pdgCode ()) != kProton )) {
598- continue ;
599- }
589+ if (calTrkEff && !isLcSignal && calEffLcEvent)
590+ continue ;
600591
601- if (pidTrkApplied && (std::abs (track.pdgCode ()) != kProton ))
602- continue ; // proton PID
592+ if (calTrkEff && cntLc == 1 ) {
593+ // genrated tracks
594+ for (const auto & track : mcParticles) {
595+ if (std::abs (track.eta ()) > etaTrackMax || track.pt () < ptTrackMin || track.pt () > ptTrackMax) {
596+ continue ;
597+ }
598+ if ((std::abs (track.pdgCode ()) != kElectron ) && (std::abs (track.pdgCode ()) != kMuonMinus ) && (std::abs (track.pdgCode ()) != kPiPlus ) && (std::abs (track.pdgCode ()) != kKPlus ) && (std::abs (track.pdgCode ()) != kProton )) {
599+ continue ;
600+ }
603601
604- if (!track.isPhysicalPrimary ()) {
605- continue ;
606- }
602+ if (pidTrkApplied && (std::abs (track.pdgCode ()) != kProton ))
603+ continue ; // proton PID
607604
608- int8_t chargeTrack = pdg->GetParticle (track.pdgCode ())->Charge (); // Retrieve charge
609- registry.fill (HIST (" hPtTracksVsSignGen" ), track.pt (),chargeTrack);
610- }
611-
605+ if (!track.isPhysicalPrimary ()) {
606+ continue ;
607+ }
608+
609+ int8_t chargeTrack = pdg->GetParticle (track.pdgCode ())->Charge (); // Retrieve charge
610+ registry.fill (HIST (" hPtTracksVsSignGen" ), track.pt (), chargeTrack);
611+ }
612612 }
613613
614614 // Lc-Hadron correlation dedicated section
@@ -625,13 +625,15 @@ struct HfCorrelatorLcHadrons {
625625 continue ;
626626 }
627627
628- if (calTrkEff && cntLc==1 && track.has_mcParticle ()){
629- auto mcParticle = track.template mcParticle_as <aod::McParticles>();
630- if (!mcParticle.isPhysicalPrimary () && isRecTrkPhyPrimary) continue ;
628+ if (calTrkEff && cntLc == 1 && track.has_mcParticle ()) {
629+ auto mcParticle = track.template mcParticle_as <aod::McParticles>();
630+ if (!mcParticle.isPhysicalPrimary () && isRecTrkPhyPrimary)
631+ continue ;
631632
632- registry.fill (HIST (" hPtTracksVsSignRec" ), track.pt (),track.sign ());
633- if (std::abs (mcParticle.pdgCode ()) == kProton ) registry.fill (HIST (" hPtTracksVsSignRecTrue" ), track.pt (),track.sign ());
634- }
633+ registry.fill (HIST (" hPtTracksVsSignRec" ), track.pt (), track.sign ());
634+ if (std::abs (mcParticle.pdgCode ()) == kProton )
635+ registry.fill (HIST (" hPtTracksVsSignRecTrue" ), track.pt (), track.sign ());
636+ }
635637
636638 // Removing Lc daughters by checking track indices
637639 if ((candidate.prong0Id () == track.globalIndex ()) || (candidate.prong1Id () == track.globalIndex ()) || (candidate.prong2Id () == track.globalIndex ())) {
0 commit comments