diff --git a/PWGLF/TableProducer/cascadebuilder.cxx b/PWGLF/TableProducer/cascadebuilder.cxx index 26bbb758875..5b3382b3596 100644 --- a/PWGLF/TableProducer/cascadebuilder.cxx +++ b/PWGLF/TableProducer/cascadebuilder.cxx @@ -1584,6 +1584,9 @@ struct cascadeBuilder { struct cascadePreselector { Produces casctags; // MC tags + // for bookkeeping + HistogramRegistry histos{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject}; + Configurable dIfMCgenerateXiMinus{"dIfMCgenerateXiMinus", true, "if MC, generate MC true XiMinus (yes/no)"}; Configurable dIfMCgenerateXiPlus{"dIfMCgenerateXiPlus", true, "if MC, generate MC true XiPlus (yes/no)"}; Configurable dIfMCgenerateOmegaMinus{"dIfMCgenerateOmegaMinus", true, "if MC, generate MC true OmegaMinus (yes/no)"}; @@ -1619,7 +1622,9 @@ struct cascadePreselector { bitdEdxOmegaPlus, bitUsedInTrackedCascade }; - void init(InitContext const&) {} + void init(InitContext const&) { + histos.add("hPreselectorStatistics", "hPreselectorStatistics", kTH1F, {{5, -0.5f, 4.5f}}); + } //*+-+*+-+*+-+*+-+*+-+*+-+*+-+*+-+*+-+*+-+* /// function to check track quality @@ -1762,19 +1767,24 @@ struct cascadePreselector { { // parse + publish tag table now for (int ii = 0; ii < selectionMask.size(); ii++) { + histos.fill(HIST("hPreselectorStatistics"), 0.0f); // All cascades bool validCascade = bitcheck(selectionMask[ii], bitTrackQuality); + histos.fill(HIST("hPreselectorStatistics"), 1.0f); // pass MC assoc (if requested) if (doprocessBuildMCAssociated || doprocessBuildValiddEdxMCAssociated) validCascade = validCascade && ((bitcheck(selectionMask[ii], bitTrueXiMinus) && dIfMCgenerateXiMinus) || (bitcheck(selectionMask[ii], bitTrueXiPlus) && dIfMCgenerateXiPlus) || (bitcheck(selectionMask[ii], bitTrueOmegaMinus) && dIfMCgenerateOmegaMinus) || (bitcheck(selectionMask[ii], bitTrueOmegaPlus) && dIfMCgenerateOmegaPlus)); + histos.fill(HIST("hPreselectorStatistics"), 2.0f); // pass MC assoc (if requested) if (doprocessBuildValiddEdx || doprocessBuildValiddEdxMCAssociated) validCascade = validCascade && ((bitcheck(selectionMask[ii], bitdEdxXiMinus) && ddEdxPreSelectXiMinus) || (bitcheck(selectionMask[ii], bitdEdxXiPlus) && ddEdxPreSelectXiPlus) || (bitcheck(selectionMask[ii], bitdEdxOmegaMinus) && ddEdxPreSelectOmegaMinus) || (bitcheck(selectionMask[ii], bitdEdxOmegaPlus) && ddEdxPreSelectOmegaPlus)); + histos.fill(HIST("hPreselectorStatistics"), 3.0f); // pass dEdx (if requested) if (doprocessSkipCascadesNotUsedInTrackedCascades) validCascade = validCascade && bitcheck(selectionMask[ii], bitUsedInTrackedCascade); + histos.fill(HIST("hPreselectorStatistics"), 4.0f); // All cascades casctags(validCascade, bitcheck(selectionMask[ii], bitTrueXiMinus), bitcheck(selectionMask[ii], bitTrueXiPlus), bitcheck(selectionMask[ii], bitTrueOmegaMinus), bitcheck(selectionMask[ii], bitTrueOmegaPlus), diff --git a/PWGLF/TableProducer/lambdakzerobuilder.cxx b/PWGLF/TableProducer/lambdakzerobuilder.cxx index c48e8a72ef9..ea190957891 100644 --- a/PWGLF/TableProducer/lambdakzerobuilder.cxx +++ b/PWGLF/TableProducer/lambdakzerobuilder.cxx @@ -251,7 +251,7 @@ struct lambdakzeroBuilder { const AxisSpec axisVsPtCoarse{static_cast(dQANBinsPtCoarse), 0, dQAMaxPt, "#it{p}_{T} (GeV/c)"}; const AxisSpec axisGammaMass{static_cast(dQANBinsMass), 0.000f, 0.400f, "Inv. Mass (GeV/c^{2})"}; const AxisSpec axisK0ShortMass{static_cast(dQANBinsMass), 0.400f, 0.600f, "Inv. Mass (GeV/c^{2})"}; - const AxisSpec axisLambdaMass{static_cast(dQANBinsMass), 1.01f, 1.21f, "Inv. Mass (GeV/c^{2})"}; + const AxisSpec axisLambdaMass{static_cast(5*dQANBinsMass), 1.01f, 2.01f, "Inv. Mass (GeV/c^{2})"}; const AxisSpec axisHypertritonMass{static_cast(dQANBinsMass), 2.900f, 3.300f, "Inv. Mass (GeV/c^{2})"}; registry.add("h2dGammaMass", "h2dGammaMass", kTH2F, {axisVsPtCoarse, axisGammaMass});