@@ -589,38 +589,38 @@ struct HfCorrelatorDplusHadrons {
589589 }
590590 counterDaughters += 1 ;
591591 prongsId[counterDaughters - 1 ] = daughI.globalIndex ();
592+ }
593+ if (!isDaughtersOk)
594+ continue ; // Skip this D+ candidate if any daughter fails eta cut
595+ counterDplusHadron++;
596+ // Dplus Hadron correlation dedicated section
597+ // if it's a Dplus particle, search for Hadron and evaluate correlations
598+ registry.fill (HIST (" hcountDplustriggersMCGen" ), 0 , particle1.pt ()); // to count trigger Dplus for normalisation)
599+ for (const auto & particleAssoc : mcParticles) {
600+ if (std::abs (particleAssoc.eta ()) > etaTrackMax || particleAssoc.pt () < ptTrackMin || particleAssoc.pt () > ptTrackMax) {
601+ continue ;
592602 }
593- if (!isDaughtersOk)
594- continue ; // Skip this D+ candidate if any daughter fails eta cut
595- counterDplusHadron++;
596- // Dplus Hadron correlation dedicated section
597- // if it's a Dplus particle, search for Hadron and evaluate correlations
598- registry.fill (HIST (" hcountDplustriggersMCGen" ), 0 , particle1.pt ()); // to count trigger Dplus for normalisation)
599- for (const auto & particleAssoc : mcParticles) {
600- if (std::abs (particleAssoc.eta ()) > etaTrackMax || particleAssoc.pt () < ptTrackMin || particleAssoc.pt () > ptTrackMax) {
601- continue ;
602- }
603- if (removeDaughters) {
604- if (particleAssoc.globalIndex () == prongsId[0 ] || particleAssoc.globalIndex () == prongsId[1 ] || particleAssoc.globalIndex () == prongsId[2 ]) {
605- continue ;
606- }
607- }
608- if ((std::abs (particleAssoc.pdgCode ()) != kElectron ) && (std::abs (particleAssoc.pdgCode ()) != kMuonMinus ) && (std::abs (particleAssoc.pdgCode ()) != kPiPlus ) && (std::abs (particleAssoc.pdgCode ()) != kKPlus ) && (std::abs (particleAssoc.pdgCode ()) != kProton )) {
609- continue ;
610- }
611- if (!particleAssoc.isPhysicalPrimary ()) {
603+ if (removeDaughters) {
604+ if (particleAssoc.globalIndex () == prongsId[0 ] || particleAssoc.globalIndex () == prongsId[1 ] || particleAssoc.globalIndex () == prongsId[2 ]) {
612605 continue ;
613606 }
607+ }
608+ if ((std::abs (particleAssoc.pdgCode ()) != kElectron ) && (std::abs (particleAssoc.pdgCode ()) != kMuonMinus ) && (std::abs (particleAssoc.pdgCode ()) != kPiPlus ) && (std::abs (particleAssoc.pdgCode ()) != kKPlus ) && (std::abs (particleAssoc.pdgCode ()) != kProton )) {
609+ continue ;
610+ }
611+ if (!particleAssoc.isPhysicalPrimary ()) {
612+ continue ;
613+ }
614614
615- int trackOrigin = RecoDecay::getCharmHadronOrigin (mcParticles, particleAssoc, true );
616- registry.fill (HIST (" hPtParticleAssocMcGen" ), particleAssoc.pt ());
617- entryDplusHadronPair (getDeltaPhi (particleAssoc.phi (), particle1.phi ()),
618- particleAssoc.eta () - particle1.eta (),
619- particle1.pt (),
620- particleAssoc.pt (),
621- poolBin);
622- entryDplusHadronRecoInfo (MassDPlus, true );
623- entryDplusHadronGenInfo (isDplusPrompt, particleAssoc.isPhysicalPrimary (), trackOrigin);
615+ int trackOrigin = RecoDecay::getCharmHadronOrigin (mcParticles, particleAssoc, true );
616+ registry.fill (HIST (" hPtParticleAssocMcGen" ), particleAssoc.pt ());
617+ entryDplusHadronPair (getDeltaPhi (particleAssoc.phi (), particle1.phi ()),
618+ particleAssoc.eta () - particle1.eta (),
619+ particle1.pt (),
620+ particleAssoc.pt (),
621+ poolBin);
622+ entryDplusHadronRecoInfo (MassDPlus, true );
623+ entryDplusHadronGenInfo (isDplusPrompt, particleAssoc.isPhysicalPrimary (), trackOrigin);
624624 } // end associated loop
625625 } // end trigger
626626 registry.fill (HIST (" hcountDplusHadronPerEvent" ), counterDplusHadron);
0 commit comments