@@ -817,14 +817,13 @@ struct kinkAnalysis {
817817
818818 if ((chargeM == -1 ) && (chargeD == -1 )) {
819819 if (cfgIsMC) {
820- histos.fill (HIST (" hcodes" ), motherPdg, daughterPdg);
821820 if ((motherPdg == particlePdgCode || motherPdg == -3222 ) && (daughterPdg == -211 )) {
822821 histos.fill (HIST (" hPtMinusRecMcTrth" ), mass, sigmaPt);
823822 histos.fill (HIST (" hptMDtrue" ), sigmaPt, PionTr.getPt ());
824823 } else if ((motherPdg == particlePdgCode || motherPdg == -3222 ) && (daughterPdg != -211 )) {
825824 histos.fill (HIST (" hptMtrue" ), sigmaPt, PionTr.getPt ());
826825 histos.fill (HIST (" hPtMinusRecMcTrthM" ), mass, sigmaPt);
827- } else if (( motherPdg == particlePdgCode || motherPdg == -1010010030 ) && ( daughterPdg == -1000010030 || daughterPdg == 1000010030 ) ) {
826+ } else if (motherPdg == -1010010030 && daughterPdg == -1000010030 ) {
828827 histos.fill (HIST (" hpRes" ), sigmaPt, (mcMotherPt - sigmaPt) / mcMotherPt);
829828 } else {
830829 histos.fill (HIST (" hptMDelse" ), sigmaPt, PionTr.getPt ());
@@ -837,6 +836,8 @@ struct kinkAnalysis {
837836 if (cfgIsMC) {
838837 if (motherPdg == particlePdgCode)
839838 histos.fill (HIST (" hPtPlusRecMcTrth" ), mass, sigmaPt);
839+ if (motherPdg == 1010010030 && daughterPdg == 1000010030 )
840+ histos.fill (HIST (" hpRes" ), sigmaPt, (mcMotherPt - sigmaPt) / mcMotherPt);
840841 }
841842 histos.fill (HIST (" hMassPlusPt" ), mass, sigmaPt);
842843 }
0 commit comments