diff --git a/PWGLF/Tasks/Strangeness/v0topologicalcuts.cxx b/PWGLF/Tasks/Strangeness/v0topologicalcuts.cxx index 4a2c8b18d91..be49007ddb1 100644 --- a/PWGLF/Tasks/Strangeness/v0topologicalcuts.cxx +++ b/PWGLF/Tasks/Strangeness/v0topologicalcuts.cxx @@ -30,41 +30,56 @@ This analysis includes two processes, one for Real Data and one for MC Data swit // namespaces to be used for the plot names and topological cuts that will be given by a configurable string namespace cuthistoskzerosh { -std::shared_ptr cospaCut[20]; +std::shared_ptr cospaCut_MC[20]; +std::shared_ptr cospaCut_Data[20]; static std::vector cospacuts; -std::shared_ptr dcaCut[20]; +std::shared_ptr dcaCut_MC[20]; +std::shared_ptr dcaCut_Data[20]; static std::vector dcacuts; -std::shared_ptr v0radiusCut[20]; +std::shared_ptr v0radiusCut_MC[20]; +std::shared_ptr v0radiusCut_Data[20]; static std::vector v0radiuscuts; -std::shared_ptr dcapostopCut[20]; +std::shared_ptr dcapostopCut_MC[20]; +std::shared_ptr dcapostopCut_Data[20]; static std::vector dcapostopvcuts; -std::shared_ptr dcanegtopCut[20]; +std::shared_ptr dcanegtopCut_MC[20]; +std::shared_ptr dcanegtopCut_Data[20]; static std::vector dcanegtopvcuts; } // namespace cuthistoskzerosh namespace cuthistoslambda { -std::shared_ptr cospaCut[20]; +std::shared_ptr cospaCut_MC[20]; +std::shared_ptr cospaCut_Data[20]; static std::vector cospacuts; -std::shared_ptr dcaCut[20]; +std::shared_ptr dcaCut_MC[20]; +std::shared_ptr dcaCut_Data[20]; static std::vector dcacuts; -std::shared_ptr v0radiusCut[20]; +std::shared_ptr v0radiusCut_MC[20]; +std::shared_ptr v0radiusCut_Data[20]; static std::vector v0radiuscuts; -std::shared_ptr dcapostopCut[20]; +std::shared_ptr dcapostopCut_MC[20]; +std::shared_ptr dcapostopCut_Data[20]; static std::vector dcapostopvcuts; -std::shared_ptr dcanegtopCut[20]; +std::shared_ptr dcanegtopCut_MC[20]; +std::shared_ptr dcanegtopCut_Data[20]; static std::vector dcanegtopvcuts; } // namespace cuthistoslambda namespace cuthistosantilambda { -std::shared_ptr cospaCut[20]; +std::shared_ptr cospaCut_MC[20]; +std::shared_ptr cospaCut_Data[20]; static std::vector cospacuts; -std::shared_ptr dcaCut[20]; +std::shared_ptr dcaCut_MC[20]; +std::shared_ptr dcaCut_Data[20]; static std::vector dcacuts; -std::shared_ptr v0radiusCut[20]; +std::shared_ptr v0radiusCut_MC[20]; +std::shared_ptr v0radiusCut_Data[20]; static std::vector v0radiuscuts; -std::shared_ptr dcapostopCut[20]; +std::shared_ptr dcapostopCut_MC[20]; +std::shared_ptr dcapostopCut_Data[20]; static std::vector dcapostopvcuts; -std::shared_ptr dcanegtopCut[20]; +std::shared_ptr dcanegtopCut_MC[20]; +std::shared_ptr dcanegtopCut_Data[20]; static std::vector dcanegtopvcuts; } // namespace cuthistosantilambda using namespace o2; @@ -80,25 +95,46 @@ struct v0topologicalcuts { HistogramRegistry rV0Parameters_Data{"rV0Parameters_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; // kzero cut Histogram Registry with MC-matching, each will include 20 histograms for 20 different cuts - HistogramRegistry rKzeroShort_cospaCut{"KzeroShort_cospaCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rKzeroShort_dcaCut{"KzeroShort_dcaCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rKzeroShort_v0radiusCut{"KzeroShort_v0radiusCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rKzeroShort_dcapostopCut{"KzeroShort_dcapostopvCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rKzeroShort_dcanegtopCut{"KzeroShort_dcanegtopvCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_cospaCut_MC{"KzeroShort_cospaCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_dcaCut_MC{"KzeroShort_dcaCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_v0radiusCut_MC{"KzeroShort_v0radiusCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_dcapostopCut_MC{"KzeroShort_dcapostopvCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_dcanegtopCut_MC{"KzeroShort_dcanegtopvCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; // lambdas cut histograms with MC-matching (same as in Kzeros above) - HistogramRegistry rLambda_cospaCut{"Lambda_cospaCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rLambda_dcaCut{"Lambda_dcaCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rLambda_v0radiusCut{"Lambda_v0radiusCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rLambda_dcapostopCut{"Lambda_dcapostopvCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rLambda_dcanegtopCut{"Lambda_dcanegtopvCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_cospaCut_MC{"Lambda_cospaCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_dcaCut_MC{"Lambda_dcaCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_v0radiusCut_MC{"Lambda_v0radiusCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_dcapostopCut_MC{"Lambda_dcapostopvCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_dcanegtopCut_MC{"Lambda_dcanegtopvCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; // antilambdas cut histograms with MC-matching (same as in Lambdas an Kzeros above) - HistogramRegistry rAntiLambda_cospaCut{"AntiLambda_cospaCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rAntiLambda_dcaCut{"AntiLambda_dcaCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rAntiLambda_v0radiusCut{"AntiLambda_v0radiusCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rAntiLambda_dcapostopCut{"AntiLambda_dcapostopvCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; - HistogramRegistry rAntiLambda_dcanegtopCut{"AntiLambda_dcanegtopvCuts", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_cospaCut_MC{"AntiLambda_cospaCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_dcaCut_MC{"AntiLambda_dcaCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_v0radiusCut_MC{"AntiLambda_v0radiusCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_dcapostopCut_MC{"AntiLambda_dcapostopvCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_dcanegtopCut_MC{"AntiLambda_dcanegtopvCuts_MC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + + // kzero cut Histogram Registry for data, each will include 20 histograms for 20 different cuts + HistogramRegistry rKzeroShort_cospaCut_Data{"KzeroShort_cospaCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_dcaCut_Data{"KzeroShort_dcaCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_v0radiusCut_Data{"KzeroShort_v0radiusCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_dcapostopCut_Data{"KzeroShort_dcapostopvCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rKzeroShort_dcanegtopCut_Data{"KzeroShort_dcanegtopvCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + + // lambdas cut histograms for data (same as in Kzeros above) + HistogramRegistry rLambda_cospaCut_Data{"Lambda_cospaCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_dcaCut_Data{"Lambda_dcaCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_v0radiusCut_Data{"Lambda_v0radiusCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_dcapostopCut_Data{"Lambda_dcapostopvCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rLambda_dcanegtopCut_Data{"Lambda_dcanegtopvCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + + // antilambdas cut histograms for data (same as in Lambdas an Kzeros above) + HistogramRegistry rAntiLambda_cospaCut_Data{"AntiLambda_cospaCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_dcaCut_Data{"AntiLambda_dcaCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_v0radiusCut_Data{"AntiLambda_v0radiusCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_dcapostopCut_Data{"AntiLambda_dcapostopvCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; + HistogramRegistry rAntiLambda_dcanegtopCut_Data{"AntiLambda_dcanegtopvCuts_Data", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; // Configurable for histograms Configurable nBins{"nBins", 100, "N bins in all histos"}; @@ -107,22 +143,22 @@ struct v0topologicalcuts { Configurable kzeroshsetting_cospacuts_string{"kzerosetting_cospacuts", {"0_98,0_981,0_982,0_983,0_984,0_985,0_986,0_987,0_988,0_989,0_99,0_991,0_992,0_993,0_994,0_995,0_996,0_997,0_998,0_999"}, "Kzero cosPA Cut Values"}; Configurable kzeroshsetting_dcacuts_string{"kzerosetting_dcacuts", {"0_3,0_285,0_27,0_255,0_24,0_225,0_21,0_195,0_18,0_165,0_15,0_135,0_12,0_105,0_09,0_075,0_06,0_045,0_03,0_015"}, "Kzero DCA Cut Values"}; Configurable kzeroshsetting_v0radius_string{"kzerosetting_v0radiuscuts", {"0_5,0_51,0_52,0_53,0_54,0_55,0_56,0_57,0_58,0_59,0_6,0_61,0_62,0_63,0_64,0_65,0_66,0_67,0_68,0_69"}, "Kzero V0Radius Cut Values"}; - Configurable kzeroshsetting_dcapostopv_string{"kzerosetting_dcapostopvcuts", {"0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Kzero DCA Pos to PV Cut Values"}; - Configurable kzeroshsetting_dcanegtopv_string{"kzerosetting_dcanegtopvcuts", {"0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "KzeroDCA Neg to PV Cut Values"}; + Configurable kzeroshsetting_dcapostopv_string{"kzerosetting_dcapostopvcuts", {"0_0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Kzero DCA Pos to PV Cut Values"}; + Configurable kzeroshsetting_dcanegtopv_string{"kzerosetting_dcanegtopvcuts", {"0_0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "KzeroDCA Neg to PV Cut Values"}; // Configurable strings for Lambdacuts Configurable lambdasetting_cospacuts_string{"lambdasetting_cospacuts", {"0_98,0_981,0_982,0_983,0_984,0_985,0_986,0_987,0_988,0_989,0_99,0_991,0_992,0_993,0_994,0_995,0_996,0_997,0_998,0_999"}, "Lambda cosPA Cut Values"}; Configurable lambdasetting_dcacuts_string{"lambdasetting_dcacuts", {"0_3,0_285,0_27,0_255,0_24,0_225,0_21,0_195,0_18,0_165,0_15,0_135,0_12,0_105,0_09,0_075,0_06,0_045,0_03,0_015"}, "Lambda DCA Cut Values"}; Configurable lambdasetting_v0radius_string{"lambdasetting_v0radiuscuts", {"0_5,0_51,0_52,0_53,0_54,0_55,0_56,0_57,0_58,0_59,0_6,0_61,0_62,0_63,0_64,0_65,0_66,0_67,0_68,0_69"}, "Lambda V0Radius Cut Values"}; - Configurable lambdasetting_dcapostopv_string{"lambdasetting_dcapostopvcuts", {"0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Lambda DCA Pos to PV Cut Values"}; - Configurable lambdasetting_dcanegtopv_string{"lambdasetting_dcanegtopvcuts", {"0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Lambda DCA Neg to PV Cut Values"}; + Configurable lambdasetting_dcapostopv_string{"lambdasetting_dcapostopvcuts", {"0_0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Lambda DCA Pos to PV Cut Values"}; + Configurable lambdasetting_dcanegtopv_string{"lambdasetting_dcanegtopvcuts", {"0_0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Lambda DCA Neg to PV Cut Values"}; // Configurable strings for AntiLambdacuts Configurable antilambdasetting_cospacuts_string{"antilambdasetting_cospacuts", {"0_98,0_981,0_982,0_983,0_984,0_985,0_986,0_987,0_988,0_989,0_99,0_991,0_992,0_993,0_994,0_995,0_996,0_997,0_998,0_999"}, "Antilambda cosPA Cut Values"}; Configurable antilambdasetting_dcacuts_string{"antilambdasetting_dcacuts", {"0_3,0_285,0_27,0_255,0_24,0_225,0_21,0_195,0_18,0_165,0_15,0_135,0_12,0_105,0_09,0_075,0_06,0_045,0_03,0_015"}, "Antilambda DCA Cut Values"}; Configurable antilambdasetting_v0radius_string{"antilambdasetting_v0radiuscuts", {"0_5,0_51,0_52,0_53,0_54,0_55,0_56,0_57,0_58,0_59,0_6,0_61,0_62,0_63,0_64,0_65,0_66,0_67,0_68,0_69"}, "Antilambda V0Radius Cut Values"}; - Configurable antilambdasetting_dcapostopv_string{"antilambdasetting_dcapostopvcuts", {"0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Antilambda DCA Pos to PV Cut Values"}; - Configurable antilambdasetting_dcanegtopv_string{"antilambdasetting_dcanegtopvcuts", {"0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Antilambda DCA Neg to PV Cut Values"}; + Configurable antilambdasetting_dcapostopv_string{"antilambdasetting_dcapostopvcuts", {"0_0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Antilambda DCA Pos to PV Cut Values"}; + Configurable antilambdasetting_dcanegtopv_string{"antilambdasetting_dcanegtopvcuts", {"0_0,0_01,0_02,0_03,0_04,0_05,0_06,0_07,0_08,0_09,0_1,0_11,0_12,0_13,0_14,0_15,0_16,0_17,0_18,0_19"}, "Antilambda DCA Neg to PV Cut Values"}; void init(InitContext const&) { @@ -245,36 +281,61 @@ struct v0topologicalcuts { } // Axes for the three invariant mass plots - AxisSpec K0ShortMassAxis = {nBins, 0.45f, 0.55f, "#it{M}_{inv} [GeV/#it{c}^{2}]"}; - AxisSpec LambdaMassAxis = {nBins, 1.085f, 1.145f, "#it{M}_{inv} [GeV/#it{c}^{2}]"}; - AxisSpec AntiLambdaMassAxis = {nBins, 1.085f, 1.145f, "#it{M}_{inv} [GeV/#it{c}^{2}]"}; + AxisSpec K0ShortMassAxis = {nBins, 0.45f, 0.55f, "#it{M} #pi^{+}#pi^{-} [GeV/#it{c}^{2}]"}; + AxisSpec LambdaMassAxis = {nBins, 1.085f, 1.145f, "#it{M} p^{+}#pi^{-} [GeV/#it{c}^{2}]"}; + AxisSpec AntiLambdaMassAxis = {nBins, 1.085f, 1.145f, "#it{M} p^{-}#pi^{+} [GeV/#it{c}^{2}]"}; - // adding the invariant mass histograms to their Registries using the namespace for kzeros, lambdas and antilambdas + // adding the MC invariant mass histograms to their Registries using the namespace for kzeros, lambdas and antilambdas for (int i = 0; i < 20; i++) { - cuthistoskzerosh::cospaCut[i] = rKzeroShort_cospaCut.add(fmt::format("hKzerocospaCut_{}", cuthistoskzerosh::cospacuts[i]).data(), fmt::format("hKzerocospaCut_{}", cuthistoskzerosh::cospacuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); - cuthistoskzerosh::dcaCut[i] = rKzeroShort_dcaCut.add(fmt::format("hKzerodcaCut_{}", cuthistoskzerosh::dcacuts[i]).data(), fmt::format("hKzerodcaCut_{}", cuthistoskzerosh::dcacuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); - cuthistoskzerosh::v0radiusCut[i] = rKzeroShort_v0radiusCut.add(fmt::format("hKzerov0radiusCut_{}", cuthistoskzerosh::v0radiuscuts[i]).data(), fmt::format("hKzerov0radiusCuts_{}", cuthistoskzerosh::v0radiuscuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); - cuthistoskzerosh::dcapostopCut[i] = rKzeroShort_dcapostopCut.add(fmt::format("hKzerodcapostopCut_{}", cuthistoskzerosh::dcapostopvcuts[i]).data(), fmt::format("hKzerodcapostopCut_{}", cuthistoskzerosh::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); - cuthistoskzerosh::dcanegtopCut[i] = rKzeroShort_dcanegtopCut.add(fmt::format("hKzerodcanegtopCut_{}", cuthistoskzerosh::dcanegtopvcuts[i]).data(), fmt::format("hKzerodcanegtopCut_{}", cuthistoskzerosh::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::cospaCut_MC[i] = rKzeroShort_cospaCut_MC.add(fmt::format("hKzerocospaCut_{}", cuthistoskzerosh::cospacuts[i]).data(), fmt::format("hKzerocospaCut_{}", cuthistoskzerosh::cospacuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::dcaCut_MC[i] = rKzeroShort_dcaCut_MC.add(fmt::format("hKzerodcaCut_{}", cuthistoskzerosh::dcacuts[i]).data(), fmt::format("hKzerodcaCut_{}", cuthistoskzerosh::dcacuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::v0radiusCut_MC[i] = rKzeroShort_v0radiusCut_MC.add(fmt::format("hKzerov0radiusCut_{}", cuthistoskzerosh::v0radiuscuts[i]).data(), fmt::format("hKzerov0radiusCut_{}", cuthistoskzerosh::v0radiuscuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::dcapostopCut_MC[i] = rKzeroShort_dcapostopCut_MC.add(fmt::format("hKzerodcapostopCut_{}", cuthistoskzerosh::dcapostopvcuts[i]).data(), fmt::format("hKzerodcapostopCut_{}", cuthistoskzerosh::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::dcanegtopCut_MC[i] = rKzeroShort_dcanegtopCut_MC.add(fmt::format("hKzerodcanegtopCut_{}", cuthistoskzerosh::dcanegtopvcuts[i]).data(), fmt::format("hKzerodcanegtopCut_{}", cuthistoskzerosh::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); } for (int i = 0; i < 20; i++) { - cuthistoslambda::cospaCut[i] = rLambda_cospaCut.add(fmt::format("hLambdacospaCut_{}", cuthistoslambda::cospacuts[i]).data(), fmt::format("hLambdacospaCut_{}", cuthistoslambda::cospacuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); - cuthistoslambda::dcaCut[i] = rLambda_dcaCut.add(fmt::format("hLambdadcaCut_{}", cuthistoslambda::dcacuts[i]).data(), fmt::format("hLambdadcaCut_{}", cuthistoslambda::dcacuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); - cuthistoslambda::v0radiusCut[i] = rLambda_v0radiusCut.add(fmt::format("hLambdav0radiusCut_{}", cuthistoslambda::v0radiuscuts[i]).data(), fmt::format("hLambdav0radiusCuts_{}", cuthistoslambda::v0radiuscuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); - cuthistoslambda::dcapostopCut[i] = rLambda_dcapostopCut.add(fmt::format("hLambdadcapostopCut_{}", cuthistoslambda::dcapostopvcuts[i]).data(), fmt::format("hLambdadcapostopCut_{}", cuthistoslambda::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); - cuthistoslambda::dcanegtopCut[i] = rLambda_dcanegtopCut.add(fmt::format("hLambdadcanegtopCut_{}", cuthistoslambda::dcanegtopvcuts[i]).data(), fmt::format("hLambdadcanegtopCut_{}", cuthistoslambda::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::cospaCut_MC[i] = rLambda_cospaCut_MC.add(fmt::format("hLambdacospaCut_{}", cuthistoslambda::cospacuts[i]).data(), fmt::format("hLambdacospaCut_{}", cuthistoslambda::cospacuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::dcaCut_MC[i] = rLambda_dcaCut_MC.add(fmt::format("hLambdadcaCut_{}", cuthistoslambda::dcacuts[i]).data(), fmt::format("hLambdadcaCut_{}", cuthistoslambda::dcacuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::v0radiusCut_MC[i] = rLambda_v0radiusCut_MC.add(fmt::format("hLambdav0radiusCut_{}", cuthistoslambda::v0radiuscuts[i]).data(), fmt::format("hLambdav0radiusCut_{}", cuthistoslambda::v0radiuscuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::dcapostopCut_MC[i] = rLambda_dcapostopCut_MC.add(fmt::format("hLambdadcapostopCut_{}", cuthistoslambda::dcapostopvcuts[i]).data(), fmt::format("hLambdadcapostopCut_{}", cuthistoslambda::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::dcanegtopCut_MC[i] = rLambda_dcanegtopCut_MC.add(fmt::format("hLambdadcanegtopCut_{}", cuthistoslambda::dcanegtopvcuts[i]).data(), fmt::format("hLambdadcanegtopCut_{}", cuthistoslambda::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); } for (int i = 0; i < 20; i++) { - cuthistosantilambda::cospaCut[i] = rAntiLambda_cospaCut.add(fmt::format("hAntiLambdacospaCut_{}", cuthistosantilambda::cospacuts[i]).data(), fmt::format("hAntiLambdacospaCut_{}", cuthistosantilambda::cospacuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); - cuthistosantilambda::dcaCut[i] = rAntiLambda_dcaCut.add(fmt::format("hAntiLambdadcaCut_{}", cuthistosantilambda::dcacuts[i]).data(), fmt::format("hAntiLambdadcaCut_{}", cuthistosantilambda::dcacuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); - cuthistosantilambda::v0radiusCut[i] = rAntiLambda_v0radiusCut.add(fmt::format("hAntiLambdav0radiusCut_{}", cuthistosantilambda::v0radiuscuts[i]).data(), fmt::format("hAntiLambdav0radiusCuts_{}", cuthistosantilambda::v0radiuscuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); - cuthistosantilambda::dcapostopCut[i] = rAntiLambda_dcapostopCut.add(fmt::format("hAntiLambdadcapostopCut_{}", cuthistosantilambda::dcapostopvcuts[i]).data(), fmt::format("hAntiLambdadcapostopCut_{}", cuthistosantilambda::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); - cuthistosantilambda::dcanegtopCut[i] = rAntiLambda_dcanegtopCut.add(fmt::format("hAntiLambdadcanegtopCut_{}", cuthistosantilambda::dcanegtopvcuts[i]).data(), fmt::format("hAntiLambdadcanegtopCut_{}", cuthistosantilambda::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::cospaCut_MC[i] = rAntiLambda_cospaCut_MC.add(fmt::format("hAntiLambdacospaCut_{}", cuthistosantilambda::cospacuts[i]).data(), fmt::format("hAntiLambdacospaCut_{}", cuthistosantilambda::cospacuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::dcaCut_MC[i] = rAntiLambda_dcaCut_MC.add(fmt::format("hAntiLambdadcaCut_{}", cuthistosantilambda::dcacuts[i]).data(), fmt::format("hAntiLambdadcaCut_{}", cuthistosantilambda::dcacuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::v0radiusCut_MC[i] = rAntiLambda_v0radiusCut_MC.add(fmt::format("hAntiLambdav0radiusCut_{}", cuthistosantilambda::v0radiuscuts[i]).data(), fmt::format("hAntiLambdav0radiusCut_{}", cuthistosantilambda::v0radiuscuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::dcapostopCut_MC[i] = rAntiLambda_dcapostopCut_MC.add(fmt::format("hAntiLambdadcapostopCut_{}", cuthistosantilambda::dcapostopvcuts[i]).data(), fmt::format("hAntiLambdadcapostopCut_{}", cuthistosantilambda::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::dcanegtopCut_MC[i] = rAntiLambda_dcanegtopCut_MC.add(fmt::format("hAntiLambdadcanegtopCut_{}", cuthistosantilambda::dcanegtopvcuts[i]).data(), fmt::format("hAntiLambdadcanegtopCut_{}", cuthistosantilambda::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); } - // K0s topological/PID cut histograms added and MC-matched + // adding the Data invariant mass histograms to their Registries using the namespace for kzeros, lambdas and antilambdas + for (int i = 0; i < 20; i++) { + cuthistoskzerosh::cospaCut_Data[i] = rKzeroShort_cospaCut_Data.add(fmt::format("hKzerocospaCut_{}", cuthistoskzerosh::cospacuts[i]).data(), fmt::format("hKzerocospaCut_{}", cuthistoskzerosh::cospacuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::dcaCut_Data[i] = rKzeroShort_dcaCut_Data.add(fmt::format("hKzerodcaCut_{}", cuthistoskzerosh::dcacuts[i]).data(), fmt::format("hKzerodcaCut_{}", cuthistoskzerosh::dcacuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::v0radiusCut_Data[i] = rKzeroShort_v0radiusCut_Data.add(fmt::format("hKzerov0radiusCut_{}", cuthistoskzerosh::v0radiuscuts[i]).data(), fmt::format("hKzerov0radiusCut_{}", cuthistoskzerosh::v0radiuscuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::dcapostopCut_Data[i] = rKzeroShort_dcapostopCut_Data.add(fmt::format("hKzerodcapostopCut_{}", cuthistoskzerosh::dcapostopvcuts[i]).data(), fmt::format("hKzerodcapostopCut_{}", cuthistoskzerosh::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + cuthistoskzerosh::dcanegtopCut_Data[i] = rKzeroShort_dcanegtopCut_Data.add(fmt::format("hKzerodcanegtopCut_{}", cuthistoskzerosh::dcanegtopvcuts[i]).data(), fmt::format("hKzerodcanegtopCut_{}", cuthistoskzerosh::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{K0ShortMassAxis}}}); + } + + for (int i = 0; i < 20; i++) { + cuthistoslambda::cospaCut_Data[i] = rLambda_cospaCut_Data.add(fmt::format("hLambdacospaCut_{}", cuthistoslambda::cospacuts[i]).data(), fmt::format("hLambdacospaCut_{}", cuthistoslambda::cospacuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::dcaCut_Data[i] = rLambda_dcaCut_Data.add(fmt::format("hLambdadcaCut_{}", cuthistoslambda::dcacuts[i]).data(), fmt::format("hLambdadcaCut_{}", cuthistoslambda::dcacuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::v0radiusCut_Data[i] = rLambda_v0radiusCut_Data.add(fmt::format("hLambdav0radiusCut_{}", cuthistoslambda::v0radiuscuts[i]).data(), fmt::format("hLambdav0radiusCut_{}", cuthistoslambda::v0radiuscuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::dcapostopCut_Data[i] = rLambda_dcapostopCut_Data.add(fmt::format("hLambdadcapostopCut_{}", cuthistoslambda::dcapostopvcuts[i]).data(), fmt::format("hLambdadcapostopCut_{}", cuthistoslambda::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + cuthistoslambda::dcanegtopCut_Data[i] = rLambda_dcanegtopCut_Data.add(fmt::format("hLambdadcanegtopCut_{}", cuthistoslambda::dcanegtopvcuts[i]).data(), fmt::format("hLambdadcanegtopCut_{}", cuthistoslambda::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{LambdaMassAxis}}}); + } + + for (int i = 0; i < 20; i++) { + cuthistosantilambda::cospaCut_Data[i] = rAntiLambda_cospaCut_Data.add(fmt::format("hAntiLambdacospaCut_{}", cuthistosantilambda::cospacuts[i]).data(), fmt::format("hAntiLambdacospaCut_{}", cuthistosantilambda::cospacuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::dcaCut_Data[i] = rAntiLambda_dcaCut_Data.add(fmt::format("hAntiLambdadcaCut_{}", cuthistosantilambda::dcacuts[i]).data(), fmt::format("hAntiLambdadcaCut_{}", cuthistosantilambda::dcacuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::v0radiusCut_Data[i] = rAntiLambda_v0radiusCut_Data.add(fmt::format("hAntiLambdav0radiusCut_{}", cuthistosantilambda::v0radiuscuts[i]).data(), fmt::format("hAntiLambdav0radiusCut_{}", cuthistosantilambda::v0radiuscuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::dcapostopCut_Data[i] = rAntiLambda_dcapostopCut_Data.add(fmt::format("hAntiLambdadcapostopCut_{}", cuthistosantilambda::dcapostopvcuts[i]).data(), fmt::format("hAntiLambdadcapostopCut_{}", cuthistosantilambda::dcapostopvcuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + cuthistosantilambda::dcanegtopCut_Data[i] = rAntiLambda_dcanegtopCut_Data.add(fmt::format("hAntiLambdadcanegtopCut_{}", cuthistosantilambda::dcanegtopvcuts[i]).data(), fmt::format("hAntiLambdadcanegtopCut_{}", cuthistosantilambda::dcanegtopvcuts[i]).data(), {HistType::kTH1D, {{AntiLambdaMassAxis}}}); + } + + // K0s topological cut histograms added and MC-matched rV0Parameters_MC_V0match.add("hDCAV0Daughters_V0_Match", "hDCAV0Daughters_No_Match", {HistType::kTH1F, {{nBins, 0.0f, 1.2f}}}); rV0Parameters_MC_V0match.add("hV0CosPA_V0_Match", "hV0CosPA_No_Match", {HistType::kTH1F, {{nBins, 0.95f, 1.f}}}); rV0Parameters_MC_V0match.add("hV0Radius_V0_Match", "hV0Radius_No_Match", {HistType::kTH1F, {{nBins, 0.0f, 5.0f}}}); @@ -307,6 +368,7 @@ struct v0topologicalcuts { rV0Parameters_Data.add("hDCAV0Daughters_V0_Data", "hDCAV0Daughters_V0_Data", {HistType::kTH1F, {{nBins, 0.0f, 1.2f}}}); rV0Parameters_Data.add("hV0CosPA_V0_Data", "hV0CosPA_V0_Data", {HistType::kTH1F, {{nBins, 0.95f, 1.f}}}); rV0Parameters_Data.add("hV0Radius_V0_Data", "hV0Radius_V0_Data", {HistType::kTH1F, {{nBins, 0.2f, 5.0f}}}); + rV0Parameters_Data.add("hV0Radius_Full_V0_Data", "hV0Radius_Full_V0_Data", {HistType::kTH1F, {{nBins, 0.2f, 5.0f}}}); rV0Parameters_Data.add("hDCAPostoPV_V0_Data", "hDCAPostoPV_V0_Data", {HistType::kTH1F, {{nBins, 0.0f, 5.0f}}}); rV0Parameters_Data.add("hDCANegtoPV_V0_Data", "hDCANegtoPV_V0_Data", {HistType::kTH1F, {{nBins, 0.0f, 5.0f}}}); rV0Parameters_Data.add("hMassK0ShortNoCuts_V0_Data", "hMassK0ShortNoCuts_V0_Data", {HistType::kTH1F, {{K0ShortMassAxis}}}); @@ -319,7 +381,7 @@ struct v0topologicalcuts { Filter eventFilter = (o2::aod::evsel::sel8 == true); // Defining the type of the daughter tracks - using DaughterTracks = soa::Join; + using DaughterTracks = soa::Join; // This is the Process for the MC reconstructed Data void RecMCprocess(soa::Filtered>::iterator const&, @@ -349,12 +411,12 @@ struct v0topologicalcuts { rV0Parameters_MC_K0Smatch.fill(HIST("hDCANegtoPV_KzeroMC_Match"), v0.dcanegtopv()); for (int j = 0; j < 20; j++) { - std::string cospacut = cuthistoskzerosh::cospacuts[j]; // Get the current cut value from the namespace - size_t pos = cospacut.find("_"); // find the "_" which needs to change to a "." for it to be a number - cospacut[pos] = '.'; // change the "_" into an "." - const float cospacutvalue = std::stod(cospacut); // make the string into a float value - if (v0.v0cosPA() > cospacutvalue) { // enforce the cut value - cuthistoskzerosh::cospaCut[j]->Fill(v0.mK0Short()); // fill the corresponding histo from the namespace with the invariant mass (of a Kzero here) + std::string cospacut = cuthistoskzerosh::cospacuts[j]; // Get the current cut value from the namespace + size_t pos = cospacut.find("_"); // find the "_" which needs to change to a "." for it to be a number + cospacut[pos] = '.'; // change the "_" into an "." + const float cospacutvalue = std::stod(cospacut); // make the string into a float value + if (v0.v0cosPA() > cospacutvalue) { // enforce the cut value + cuthistoskzerosh::cospaCut_MC[j]->Fill(v0.mK0Short()); // fill the corresponding histo from the namespace with the invariant mass (of a Kzero here) } } for (int j = 0; j < 20; j++) { @@ -363,7 +425,7 @@ struct v0topologicalcuts { dcacut[pos] = '.'; const float dcacutvalue = std::stod(dcacut); if (v0.dcaV0daughters() < dcacutvalue) { - cuthistoskzerosh::dcaCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::dcaCut_MC[j]->Fill(v0.mK0Short()); } } for (int j = 0; j < 20; j++) { @@ -372,7 +434,7 @@ struct v0topologicalcuts { v0radiuscut[pos] = '.'; const float v0radiuscutvalue = std::stod(v0radiuscut); if (v0.v0radius() > v0radiuscutvalue) { - cuthistoskzerosh::v0radiusCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::v0radiusCut_MC[j]->Fill(v0.mK0Short()); } } for (int j = 0; j < 20; j++) { @@ -381,7 +443,7 @@ struct v0topologicalcuts { dcapostopcut[pos] = '.'; const float dcapostopcutvalue = std::stod(dcapostopcut); if (v0.dcapostopv() > dcapostopcutvalue) { - cuthistoskzerosh::dcapostopCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::dcapostopCut_MC[j]->Fill(v0.mK0Short()); } } for (int j = 0; j < 20; j++) { @@ -390,7 +452,7 @@ struct v0topologicalcuts { dcanegtopcut[pos] = '.'; const float dcanegtopcutvalue = std::stod(dcanegtopcut); if (v0.dcanegtopv() > dcanegtopcutvalue) { - cuthistoskzerosh::dcanegtopCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::dcanegtopCut_MC[j]->Fill(v0.mK0Short()); } } } @@ -408,7 +470,7 @@ struct v0topologicalcuts { cospacutlambda[pos] = '.'; const float cospacutlambdavalue = std::stod(cospacutlambda); if (v0.v0cosPA() > cospacutlambdavalue) { - cuthistoslambda::cospaCut[j]->Fill(v0.mLambda()); + cuthistoslambda::cospaCut_MC[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -417,7 +479,7 @@ struct v0topologicalcuts { dcacutlambda[pos] = '.'; const float dcacutlambdavalue = std::stod(dcacutlambda); if (v0.dcaV0daughters() < dcacutlambdavalue) { - cuthistoslambda::dcaCut[j]->Fill(v0.mLambda()); + cuthistoslambda::dcaCut_MC[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -426,7 +488,7 @@ struct v0topologicalcuts { v0radiuscutlambda[pos] = '.'; const float v0radiuscutlambdavalue = std::stod(v0radiuscutlambda); if (v0.v0radius() > v0radiuscutlambdavalue) { - cuthistoslambda::v0radiusCut[j]->Fill(v0.mLambda()); + cuthistoslambda::v0radiusCut_MC[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -435,7 +497,7 @@ struct v0topologicalcuts { dcapostopcutlambda[pos] = '.'; const float dcapostopcutlambdavalue = std::stod(dcapostopcutlambda); if (v0.dcapostopv() > dcapostopcutlambdavalue) { - cuthistoslambda::dcapostopCut[j]->Fill(v0.mLambda()); + cuthistoslambda::dcapostopCut_MC[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -444,7 +506,7 @@ struct v0topologicalcuts { dcanegtopcutlambda[pos] = '.'; const float dcanegtopcutlambdavalue = std::stod(dcanegtopcutlambda); if (v0.dcanegtopv() > dcanegtopcutlambdavalue) { - cuthistoslambda::dcanegtopCut[j]->Fill(v0.mLambda()); + cuthistoslambda::dcanegtopCut_MC[j]->Fill(v0.mLambda()); } } } @@ -461,7 +523,7 @@ struct v0topologicalcuts { cospacutantilambda[pos] = '.'; const float cospacutantilambdavalue = std::stod(cospacutantilambda); if (v0.v0cosPA() > cospacutantilambdavalue) { - cuthistosantilambda::cospaCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::cospaCut_MC[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -470,7 +532,7 @@ struct v0topologicalcuts { dcacutantilambda[pos] = '.'; const float dcacutantilambdavalue = std::stod(dcacutantilambda); if (v0.dcaV0daughters() < dcacutantilambdavalue) { - cuthistosantilambda::dcaCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::dcaCut_MC[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -479,7 +541,7 @@ struct v0topologicalcuts { v0radiusantilambda[pos] = '.'; const float v0radiuscutantilambdavalue = std::stod(v0radiusantilambda); if (v0.v0radius() > v0radiuscutantilambdavalue) { - cuthistosantilambda::v0radiusCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::v0radiusCut_MC[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -488,7 +550,7 @@ struct v0topologicalcuts { dcapostopantilambda[pos] = '.'; const float dcapostopcutantilambdavalue = std::stod(dcapostopantilambda); if (v0.dcapostopv() > dcapostopcutantilambdavalue) { - cuthistosantilambda::dcapostopCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::dcapostopCut_MC[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -497,7 +559,7 @@ struct v0topologicalcuts { dcanegtopantilambda[pos] = '.'; const float dcanegtopcutantilambdavalue = std::stod(dcanegtopantilambda); if (v0.dcanegtopv() > dcanegtopcutantilambdavalue) { - cuthistosantilambda::dcanegtopCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::dcanegtopCut_MC[j]->Fill(v0.mAntiLambda()); } } } @@ -512,7 +574,7 @@ struct v0topologicalcuts { for (const auto& v0 : V0s) { rV0Parameters_Data.fill(HIST("hMassK0ShortNoCuts_V0_Data"), v0.mK0Short()); rV0Parameters_Data.fill(HIST("hMassLambdaNoCuts_V0_Data"), v0.mLambda()); - rV0Parameters_Data.fill(HIST("hMassAntiLambdaNoCuts_V0_Data"), v0.mAntiLambda()); + rV0Parameters_Data.fill(HIST("hMassAntilambdaNoCuts_V0_Data"), v0.mAntiLambda()); rV0Parameters_Data.fill(HIST("hDCAV0Daughters_V0_Data"), v0.dcaV0daughters()); rV0Parameters_Data.fill(HIST("hV0CosPA_V0_Data"), v0.v0cosPA()); rV0Parameters_Data.fill(HIST("hV0Radius_V0_Data"), v0.v0radius()); @@ -527,7 +589,7 @@ struct v0topologicalcuts { cospacut[pos] = '.'; const float cospacutvalue = std::stod(cospacut); if (v0.v0cosPA() > cospacutvalue) { - cuthistoskzerosh::cospaCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::cospaCut_Data[j]->Fill(v0.mK0Short()); } } for (int j = 0; j < 20; j++) { @@ -536,7 +598,7 @@ struct v0topologicalcuts { dcacut[pos] = '.'; const float dcacutvalue = std::stod(dcacut); if (v0.dcaV0daughters() < dcacutvalue) { - cuthistoskzerosh::dcaCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::dcaCut_Data[j]->Fill(v0.mK0Short()); } } for (int j = 0; j < 20; j++) { @@ -545,7 +607,7 @@ struct v0topologicalcuts { v0radiuscut[pos] = '.'; const float v0radiuscutvalue = std::stod(v0radiuscut); if (v0.v0radius() > v0radiuscutvalue) { - cuthistoskzerosh::v0radiusCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::v0radiusCut_Data[j]->Fill(v0.mK0Short()); } } for (int j = 0; j < 20; j++) { @@ -554,7 +616,7 @@ struct v0topologicalcuts { dcapostopcut[pos] = '.'; const float dcapostopcutvalue = std::stod(dcapostopcut); if (v0.dcapostopv() > dcapostopcutvalue) { - cuthistoskzerosh::dcapostopCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::dcapostopCut_Data[j]->Fill(v0.mK0Short()); } } for (int j = 0; j < 20; j++) { @@ -563,7 +625,7 @@ struct v0topologicalcuts { dcanegtopcut[pos] = '.'; const float dcanegtopcutvalue = std::stod(dcanegtopcut); if (v0.dcanegtopv() > dcanegtopcutvalue) { - cuthistoskzerosh::dcanegtopCut[j]->Fill(v0.mK0Short()); + cuthistoskzerosh::dcanegtopCut_Data[j]->Fill(v0.mK0Short()); } } // Filling the five Lambda invariant mass plots for different cuts (which are taken from namespace), same as with Kzeros above,for full explanation see the first kzero cut filling in the MC process @@ -573,7 +635,7 @@ struct v0topologicalcuts { cospacutlambda[pos] = '.'; const float cospacutlambdavalue = std::stod(cospacutlambda); if (v0.v0cosPA() > cospacutlambdavalue) { - cuthistoslambda::cospaCut[j]->Fill(v0.mLambda()); + cuthistoslambda::cospaCut_Data[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -582,7 +644,7 @@ struct v0topologicalcuts { dcacutlambda[pos] = '.'; const float dcacutlambdavalue = std::stod(dcacutlambda); if (v0.dcaV0daughters() < dcacutlambdavalue) { - cuthistoslambda::dcaCut[j]->Fill(v0.mLambda()); + cuthistoslambda::dcaCut_Data[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -591,7 +653,7 @@ struct v0topologicalcuts { v0radiuscutlambda[pos] = '.'; const float v0radiuscutlambdavalue = std::stod(v0radiuscutlambda); if (v0.v0radius() > v0radiuscutlambdavalue) { - cuthistoslambda::v0radiusCut[j]->Fill(v0.mLambda()); + cuthistoslambda::v0radiusCut_Data[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -600,7 +662,7 @@ struct v0topologicalcuts { dcapostopcutlambda[pos] = '.'; const float dcapostopcutlambdavalue = std::stod(dcapostopcutlambda); if (v0.dcapostopv() > dcapostopcutlambdavalue) { - cuthistoslambda::dcapostopCut[j]->Fill(v0.mLambda()); + cuthistoslambda::dcapostopCut_Data[j]->Fill(v0.mLambda()); } } for (int j = 0; j < 20; j++) { @@ -609,17 +671,17 @@ struct v0topologicalcuts { dcanegtopcutlambda[pos] = '.'; const float dcanegtopcutlambdavalue = std::stod(dcanegtopcutlambda); if (v0.dcanegtopv() > dcanegtopcutlambdavalue) { - cuthistoslambda::dcanegtopCut[j]->Fill(v0.mLambda()); + cuthistoslambda::dcanegtopCut_Data[j]->Fill(v0.mLambda()); } } - // Filling the five Lambda invariant mass plots for different cuts (which are taken from namespace), same as with Kzeros and Lambdas above,for full explanation see the first kzero cut filling in the MC process + // Filling the five Anti-Lambda invariant mass plots for different cuts (which are taken from namespace), same as with Kzeros and Lambdas above,for full explanation see the first kzero cut filling in the MC process for (int j = 0; j < 20; j++) { std::string cospacutantilambda = cuthistosantilambda::cospacuts[j]; size_t pos = cospacutantilambda.find("_"); cospacutantilambda[pos] = '.'; const float cospacutantilambdavalue = std::stod(cospacutantilambda); if (v0.v0cosPA() > cospacutantilambdavalue) { - cuthistosantilambda::cospaCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::cospaCut_Data[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -628,7 +690,7 @@ struct v0topologicalcuts { dcacutantilambda[pos] = '.'; const float dcacutantilambdavalue = std::stod(dcacutantilambda); if (v0.dcaV0daughters() < dcacutantilambdavalue) { - cuthistosantilambda::dcaCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::dcaCut_Data[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -637,7 +699,7 @@ struct v0topologicalcuts { v0radiusantilambda[pos] = '.'; const float v0radiuscutantilambdavalue = std::stod(v0radiusantilambda); if (v0.v0radius() > v0radiuscutantilambdavalue) { - cuthistosantilambda::v0radiusCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::v0radiusCut_Data[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -646,7 +708,7 @@ struct v0topologicalcuts { dcapostopantilambda[pos] = '.'; const float dcapostopcutantilambdavalue = std::stod(dcapostopantilambda); if (v0.dcapostopv() > dcapostopcutantilambdavalue) { - cuthistosantilambda::dcapostopCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::dcapostopCut_Data[j]->Fill(v0.mAntiLambda()); } } for (int j = 0; j < 20; j++) { @@ -655,12 +717,11 @@ struct v0topologicalcuts { dcanegtopantilambda[pos] = '.'; const float dcanegtopcutantilambdavalue = std::stod(dcanegtopantilambda); if (v0.dcanegtopv() > dcanegtopcutantilambdavalue) { - cuthistosantilambda::dcanegtopCut[j]->Fill(v0.mAntiLambda()); + cuthistosantilambda::dcanegtopCut_Data[j]->Fill(v0.mAntiLambda()); } } } } - PROCESS_SWITCH(v0topologicalcuts, RecMCprocess, "Process Run 3 MC:Reconstructed", true); PROCESS_SWITCH(v0topologicalcuts, Dataprocess, "Process Run 3 Data,", false); };