Skip to content

Commit 1e6dd70

Browse files
mfagginMattia Faggin
authored and
smaff92
committed
PWGHF: fix filling of some MC histos and re-order inv. mass histos for signal. (AliceO2Group#3889)
Co-authored-by: Mattia Faggin <mfaggin@cern.ch>
1 parent 3472f35 commit 1e6dd70

File tree

1 file changed

+41
-47
lines changed

1 file changed

+41
-47
lines changed

‎PWGHF/D2H/Tasks/taskSigmac.cxx‎

Lines changed: 41 additions & 47 deletions
Original file line numberDiff line numberDiff line change
@@ -605,7 +605,7 @@ struct HfTaskSigmac {
605605
registry.fill(HIST("MC/reconstructed/hPtGenSc0Sig"), ptGenSc, origin, channel);
606606
registry.fill(HIST("MC/reconstructed/hEtaSc0Sig"), etaSc, origin, channel);
607607
registry.fill(HIST("MC/reconstructed/hPhiSc0Sig"), phiSc, origin, channel);
608-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0Sig"), deltaMass, ptSc);
608+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0Sig"), deltaMass, ptSc, channel);
609609
////////////////////////////////////////////////////////////////////////////////////////////////////// Σc0 signal
610610
registry.fill(HIST("MC/reconstructed/hPtSoftPiSc0Sig"), ptSoftPi, origin, channel);
611611
registry.fill(HIST("MC/reconstructed/hPtGenSoftPiSc0Sig"), ptGenSoftPi, origin, channel);
@@ -615,15 +615,15 @@ struct HfTaskSigmac {
615615
registry.fill(HIST("MC/reconstructed/hPtGenLcFromSc0Sig"), ptGenLc, origin, channel);
616616
registry.fill(HIST("MC/reconstructed/hEtaLcFromSc0Sig"), etaLc, origin, channel);
617617
registry.fill(HIST("MC/reconstructed/hPhiLcFromSc0Sig"), phiLc, origin, channel);
618-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0Sig"), deltaMass, ptLc);
618+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0Sig"), deltaMass, ptLc, channel);
619619
////////////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc0 signal
620620

621621
/// Fill the histograms for reconstructed Σc0,++ signal
622622
registry.fill(HIST("MC/reconstructed/hPtSc0PlusPlusSig"), ptSc, origin, channel);
623623
registry.fill(HIST("MC/reconstructed/hPtGenSc0PlusPlusSig"), ptGenSc, origin, channel);
624624
registry.fill(HIST("MC/reconstructed/hEtaSc0PlusPlusSig"), etaSc, origin, channel);
625625
registry.fill(HIST("MC/reconstructed/hPhiSc0PlusPlusSig"), phiSc, origin, channel);
626-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSig"), deltaMass, ptSc);
626+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSig"), deltaMass, ptSc, channel);
627627
////////////////////////////////////////////////////////////////////////////////////////////////////// Σc0,++ signal
628628
registry.fill(HIST("MC/reconstructed/hPtSoftPiSc0PlusPlusSig"), ptSoftPi, origin, channel);
629629
registry.fill(HIST("MC/reconstructed/hPtGenSoftPiSc0PlusPlusSig"), ptGenSoftPi, origin, channel);
@@ -638,15 +638,15 @@ struct HfTaskSigmac {
638638

639639
if (origin == RecoDecay::OriginType::Prompt) {
640640
/// prompt signal
641-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0SigPrompt"), deltaMass, ptSc); // Σc0 signal
642-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0SigPrompt"), deltaMass, ptLc); // Λc+ ← Σc0 signal
643-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigPrompt"), deltaMass, ptSc); // Σc0,++ signal
641+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0SigPrompt"), deltaMass, ptSc, channel); // Σc0 signal
642+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0SigPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0 signal
643+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigPrompt"), deltaMass, ptSc, channel); // Σc0,++ signal
644644
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0,++ signal
645645
} else if (origin == RecoDecay::OriginType::NonPrompt) {
646646
/// non-prompt signal
647-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0SigNonPrompt"), deltaMass, ptSc); // Σc0 signal
648-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0SigNonPrompt"), deltaMass, ptLc); // Λc+ ← Σc0 signal
649-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigNonPrompt"), deltaMass, ptSc); // Σc0,++ signal
647+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0SigNonPrompt"), deltaMass, ptSc, channel); // Σc0 signal
648+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0SigNonPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0 signal
649+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigNonPrompt"), deltaMass, ptSc, channel); // Σc0,++ signal
650650
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigNonPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0,++ signal
651651
}
652652

@@ -757,37 +757,25 @@ struct HfTaskSigmac {
757757
registry.fill(HIST("MC/reconstructed/hEtaScPlusPlusSig"), etaSc, origin, channel);
758758
registry.fill(HIST("MC/reconstructed/hPhiScPlusPlusSig"), phiSc, origin, channel);
759759
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSig"), deltaMass, ptSc, channel);
760-
if (origin == RecoDecay::OriginType::Prompt) {
761-
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigPrompt"), deltaMass, ptSc, channel);
762-
} else if (origin == RecoDecay::OriginType::NonPrompt) {
763-
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigNonPrompt"), deltaMass, ptSc, channel);
764-
} //////////////////////////////////////////////////////////////////////////////////////////////// Σc0 signal
760+
//////////////////////////////////////////////////////////////////////////////////////////////// Σc++ signal
765761
registry.fill(HIST("MC/reconstructed/hPtSoftPiScPlusPlusSig"), ptSoftPi, origin, channel);
766762
registry.fill(HIST("MC/reconstructed/hPtGenSoftPiScPlusPlusSig"), ptGenSoftPi, origin, channel);
767763
registry.fill(HIST("MC/reconstructed/hEtaSoftPiScPlusPlusSig"), etaSoftPi, origin, channel);
768-
registry.fill(HIST("MC/reconstructed/hPhiSoftPiScPlusPlusSig"), phiSoftPi, origin, channel); // π ← Σc0 signal
764+
registry.fill(HIST("MC/reconstructed/hPhiSoftPiScPlusPlusSig"), phiSoftPi, origin, channel); // π ← Σc++ signal
769765
registry.fill(HIST("MC/reconstructed/hPtLcFromScPlusPlusSig"), ptLc, origin, channel);
770766
registry.fill(HIST("MC/reconstructed/hPtGenLcFromScPlusPlusSig"), ptGenLc, origin, channel);
771767
registry.fill(HIST("MC/reconstructed/hEtaLcFromScPlusPlusSig"), etaLc, origin, channel);
772768
registry.fill(HIST("MC/reconstructed/hPhiLcFromScPlusPlusSig"), phiLc, origin, channel);
773769
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSig"), deltaMass, ptLc, channel);
774-
if (origin == RecoDecay::OriginType::Prompt) {
775-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigPrompt"), deltaMass, ptLc, channel);
776-
} else if (origin == RecoDecay::OriginType::NonPrompt) {
777-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigNonPrompt"), deltaMass, ptLc, channel);
778-
} //////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc0 signal
770+
//////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc++ signal
779771

780772
/// Fill the histograms for reconstructed Σc0,++ signal
781773
registry.fill(HIST("MC/reconstructed/hPtSc0PlusPlusSig"), ptSc, origin, channel);
782774
registry.fill(HIST("MC/reconstructed/hPtGenSc0PlusPlusSig"), ptGenSc, origin, channel);
783775
registry.fill(HIST("MC/reconstructed/hEtaSc0PlusPlusSig"), etaSc, origin, channel);
784776
registry.fill(HIST("MC/reconstructed/hPhiSc0PlusPlusSig"), phiSc, origin, channel);
785777
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSig"), deltaMass, ptSc, channel);
786-
if (origin == RecoDecay::OriginType::Prompt) {
787-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigPrompt"), deltaMass, ptSc, channel);
788-
} else if (origin == RecoDecay::OriginType::NonPrompt) {
789-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigNonPrompt"), deltaMass, ptSc, channel);
790-
} /////////////////////////////////////////////////////////////////////////////////////////////// Σc0,++ signal
778+
/////////////////////////////////////////////////////////////////////////////////////////////// Σc0,++ signal
791779
registry.fill(HIST("MC/reconstructed/hPtSoftPiSc0PlusPlusSig"), ptSoftPi, origin, channel);
792780
registry.fill(HIST("MC/reconstructed/hPtGenSoftPiSc0PlusPlusSig"), ptGenSoftPi, origin, channel);
793781
registry.fill(HIST("MC/reconstructed/hEtaSoftPiSc0PlusPlusSig"), etaSoftPi, origin, channel);
@@ -797,11 +785,21 @@ struct HfTaskSigmac {
797785
registry.fill(HIST("MC/reconstructed/hEtaLcFromSc0PlusPlusSig"), etaLc, origin, channel);
798786
registry.fill(HIST("MC/reconstructed/hPhiLcFromSc0PlusPlusSig"), phiLc, origin, channel);
799787
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSig"), deltaMass, ptLc, channel);
788+
////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc0,++ signal
789+
800790
if (origin == RecoDecay::OriginType::Prompt) {
801-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigPrompt"), deltaMass, ptLc, channel);
791+
/// prompt signal
792+
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigPrompt"), deltaMass, ptSc, channel); // Σc++ signal
793+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc++ signal
794+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigPrompt"), deltaMass, ptSc, channel); // Σc0,++ signal
795+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0,++ signal
802796
} else if (origin == RecoDecay::OriginType::NonPrompt) {
803-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigNonPrompt"), deltaMass, ptLc, channel);
804-
} ////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc0,++ signal
797+
/// non-prompt signal
798+
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigNonPrompt"), deltaMass, ptSc, channel); // Σc++ signal
799+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigNonPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc++ signal
800+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigNonPrompt"), deltaMass, ptSc, channel); // Σc0,++ signal
801+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigNonPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0,++ signal
802+
}
805803

806804
/// THn for candidate Σc0,++ cut variation
807805
if (enableTHn) {
@@ -821,37 +819,25 @@ struct HfTaskSigmac {
821819
registry.fill(HIST("MC/reconstructed/hEtaScPlusPlusSig"), etaSc, origin, channel);
822820
registry.fill(HIST("MC/reconstructed/hPhiScPlusPlusSig"), phiSc, origin, channel);
823821
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSig"), deltaMass, ptSc, channel);
824-
if (origin == RecoDecay::OriginType::Prompt) {
825-
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigPrompt"), deltaMass, ptSc, channel);
826-
} else if (origin == RecoDecay::OriginType::NonPrompt) {
827-
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigNonPrompt"), deltaMass, ptSc, channel);
828-
} ////////////////////////////////////////////////////////////////////////////////////////////////////// Σc0 signal
822+
////////////////////////////////////////////////////////////////////////////////////////////////////// Σc++ signal
829823
registry.fill(HIST("MC/reconstructed/hPtSoftPiScPlusPlusSig"), ptSoftPi, origin, channel);
830824
registry.fill(HIST("MC/reconstructed/hPtGenSoftPiScPlusPlusSig"), ptGenSoftPi, origin, channel);
831825
registry.fill(HIST("MC/reconstructed/hEtaSoftPiScPlusPlusSig"), etaSoftPi, origin, channel);
832-
registry.fill(HIST("MC/reconstructed/hPhiSoftPiScPlusPlusSig"), phiSoftPi, origin, channel); // π ← Σc0 signal
826+
registry.fill(HIST("MC/reconstructed/hPhiSoftPiScPlusPlusSig"), phiSoftPi, origin, channel); // π ← Σc++ signal
833827
registry.fill(HIST("MC/reconstructed/hPtLcFromScPlusPlusSig"), ptLc, origin, channel);
834828
registry.fill(HIST("MC/reconstructed/hPtGenLcFromScPlusPlusSig"), ptGenLc, origin, channel);
835829
registry.fill(HIST("MC/reconstructed/hEtaLcFromScPlusPlusSig"), etaLc, origin, channel);
836830
registry.fill(HIST("MC/reconstructed/hPhiLcFromScPlusPlusSig"), phiLc, origin, channel);
837831
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSig"), deltaMass, ptLc, channel);
838-
if (origin == RecoDecay::OriginType::Prompt) {
839-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigPrompt"), deltaMass, ptLc, channel);
840-
} else if (origin == RecoDecay::OriginType::NonPrompt) {
841-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigNonPrompt"), deltaMass, ptLc, channel);
842-
} ///////////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc0 signal
832+
///////////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc++ signal
843833

844834
/// Fill the histograms for reconstructed Σc0,++ signal
845835
registry.fill(HIST("MC/reconstructed/hPtSc0PlusPlusSig"), ptSc, origin, channel);
846836
registry.fill(HIST("MC/reconstructed/hPtGenSc0PlusPlusSig"), ptGenSc, origin, channel);
847837
registry.fill(HIST("MC/reconstructed/hEtaSc0PlusPlusSig"), etaSc, origin, channel);
848838
registry.fill(HIST("MC/reconstructed/hPhiSc0PlusPlusSig"), phiSc, origin, channel);
849839
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSig"), deltaMass, ptSc, channel);
850-
if (origin == RecoDecay::OriginType::Prompt) {
851-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigPrompt"), deltaMass, ptSc, channel);
852-
} else if (origin == RecoDecay::OriginType::NonPrompt) {
853-
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigNonPrompt"), deltaMass, ptSc, channel);
854-
} ///////////////////////////////////////////////////////////////////////////////////////////////////// Σc0,++ signal
840+
///////////////////////////////////////////////////////////////////////////////////////////////////// Σc0,++ signal
855841
registry.fill(HIST("MC/reconstructed/hPtSoftPiSc0PlusPlusSig"), ptSoftPi, origin, channel);
856842
registry.fill(HIST("MC/reconstructed/hPtGenSoftPiSc0PlusPlusSig"), ptGenSoftPi, origin, channel);
857843
registry.fill(HIST("MC/reconstructed/hEtaSoftPiSc0PlusPlusSig"), etaSoftPi, origin, channel);
@@ -861,11 +847,19 @@ struct HfTaskSigmac {
861847
registry.fill(HIST("MC/reconstructed/hEtaLcFromSc0PlusPlusSig"), etaLc, origin, channel);
862848
registry.fill(HIST("MC/reconstructed/hPhiLcFromSc0PlusPlusSig"), phiLc, origin, channel);
863849
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSig"), deltaMass, ptLc, channel);
850+
//////////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc0,++ signal
851+
864852
if (origin == RecoDecay::OriginType::Prompt) {
865-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigPrompt"), deltaMass, ptLc, channel);
853+
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigPrompt"), deltaMass, ptSc, channel); // Σc++ signal
854+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc++ signal
855+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigPrompt"), deltaMass, ptSc, channel); // Σc0,++ signal
856+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0,++ signal
866857
} else if (origin == RecoDecay::OriginType::NonPrompt) {
867-
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigNonPrompt"), deltaMass, ptLc, channel);
868-
} //////////////////////////////////////////////////////////////////////////////////////////////////// Λc+ ← Σc0,++ signal
858+
registry.fill(HIST("MC/reconstructed/hDeltaMassScPlusPlusSigNonPrompt"), deltaMass, ptSc, channel); // Σc++ signal
859+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromScPlusPlusSigNonPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc++ signal
860+
registry.fill(HIST("MC/reconstructed/hDeltaMassSc0PlusPlusSigNonPrompt"), deltaMass, ptSc, channel); // Σc0,++ signal
861+
registry.fill(HIST("MC/reconstructed/hDeltaMassLcFromSc0PlusPlusSigNonPrompt"), deltaMass, ptLc, channel); // Λc+ ← Σc0,++ signal
862+
}
869863

870864
/// THn for candidate Σc0,++ cut variation
871865
if (enableTHn) {

0 commit comments

Comments
 (0)