Skip to content

Commit 4bb7594

Browse files
ariedel-cernariedel-cern
authored
[PWGCF] Update femto framework (#14373)
1 parent e0c87b8 commit 4bb7594

32 files changed

+1997
-1267
lines changed

‎PWGCF/Femto/Core/cascadeBuilder.h‎

Lines changed: 11 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -90,17 +90,17 @@ struct ConfOmegaBits : o2::framework::ConfigurableGroup {
9090

9191
#undef CASCADE_DEFAULT_BITS
9292

93-
#define CASCADE_DEFAULT_SELECTION(defaultMassMin, defaultMassMax, defaultPdgCode) \
94-
o2::framework::Configurable<int> pdgCode{"pdgCode", defaultPdgCode, "Track PDG code"}; \
95-
o2::framework::Configurable<int> sign{"sign", 1, "Sign of the charge of the Cascade "}; \
96-
o2::framework::Configurable<float> ptMin{"ptMin", 0.f, "Minimum pT"}; \
97-
o2::framework::Configurable<float> ptMax{"ptMax", 999.f, "Maximum pT"}; \
98-
o2::framework::Configurable<float> etaMin{"etaMin", -10.f, "Minimum eta"}; \
99-
o2::framework::Configurable<float> etaMax{"etaMax", 10.f, "Maximum eta"}; \
100-
o2::framework::Configurable<float> phiMin{"phiMin", 0.f, "Minimum eta"}; \
101-
o2::framework::Configurable<float> phiMax{"phiMax", 1.f * o2::constants::math::TwoPI, "Maximum phi"}; \
102-
o2::framework::Configurable<float> massMin{"massMin", defaultMassMin, "Minimum invariant mass for Cascade"}; \
103-
o2::framework::Configurable<float> massMax{"massMax", defaultMassMax, "Maximum invariant mass for Cascade"}; \
93+
#define CASCADE_DEFAULT_SELECTION(defaultMassMin, defaultMassMax, defaultPdgCode) \
94+
o2::framework::Configurable<int> pdgCodeAbs{"pdgCodeAbs", defaultPdgCode, "Cascade PDG code. Set sign to +1 to select antiparticle"}; \
95+
o2::framework::Configurable<int> sign{"sign", -1, "Sign of the charge of the Cascade"}; \
96+
o2::framework::Configurable<float> ptMin{"ptMin", 0.f, "Minimum pT"}; \
97+
o2::framework::Configurable<float> ptMax{"ptMax", 999.f, "Maximum pT"}; \
98+
o2::framework::Configurable<float> etaMin{"etaMin", -10.f, "Minimum eta"}; \
99+
o2::framework::Configurable<float> etaMax{"etaMax", 10.f, "Maximum eta"}; \
100+
o2::framework::Configurable<float> phiMin{"phiMin", 0.f, "Minimum eta"}; \
101+
o2::framework::Configurable<float> phiMax{"phiMax", 1.f * o2::constants::math::TwoPI, "Maximum phi"}; \
102+
o2::framework::Configurable<float> massMin{"massMin", defaultMassMin, "Minimum invariant mass for Cascade"}; \
103+
o2::framework::Configurable<float> massMax{"massMax", defaultMassMax, "Maximum invariant mass for Cascade"}; \
104104
o2::framework::Configurable<o2::aod::femtodatatypes::CascadeMaskType> mask{"mask", 0x0, "Bitmask for cascade selection"};
105105

106106
struct ConfXiSelection : o2::framework::ConfigurableGroup {

‎PWGCF/Femto/Core/cascadeHistManager.h‎

Lines changed: 98 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -173,41 +173,122 @@ class CascadeHistManager
173173
CascadeHistManager() = default;
174174
~CascadeHistManager() = default;
175175

176-
template <modes::Mode mode>
176+
template <modes::Mode mode, typename T>
177177
void init(o2::framework::HistogramRegistry* registry,
178178
std::map<CascadeHist, std::vector<o2::framework::AxisSpec>> const& cascadeSpecs,
179+
T const& ConfCascadeSelection,
179180
std::map<trackhistmanager::TrackHist, std::vector<o2::framework::AxisSpec>> const& BachelorSpecs,
180181
std::map<trackhistmanager::TrackHist, std::vector<o2::framework::AxisSpec>> const& PosDauSpecs,
181182
std::map<trackhistmanager::TrackHist, std::vector<o2::framework::AxisSpec>> const& NegDauSpecs)
182183
{
183184
mHistogramRegistry = registry;
184-
mBachelorManager.template init<mode>(registry, BachelorSpecs, AbsChargeDaughters);
185-
mPosDauManager.template init<mode>(registry, PosDauSpecs, AbsChargeDaughters);
186-
mNegDauManager.template init<mode>(registry, NegDauSpecs, AbsChargeDaughters);
185+
mPdgCode = std::abs(ConfCascadeSelection.pdgCodeAbs.value);
186+
187+
int bachelorPdgCodeAbs = 0;
188+
int posDauPdgCodeAbs = 0;
189+
int negDauPdgCodeAbs = 0;
190+
const int absCharge = 1;
191+
int signBachelor = 0;
192+
const int signPlus = 1;
193+
const int signMinus = -1;
194+
195+
if (mPdgCode == PDG_t::kXiMinus) {
196+
if (ConfCascadeSelection.sign.value < 0) {
197+
bachelorPdgCodeAbs = std::abs(PDG_t::kPiMinus);
198+
signBachelor = -1;
199+
posDauPdgCodeAbs = std::abs(PDG_t::kProton);
200+
negDauPdgCodeAbs = std::abs(PDG_t::kPiMinus);
201+
} else {
202+
mPdgCode = -1 * mPdgCode; // Xi+ has negative pdg code
203+
bachelorPdgCodeAbs = std::abs(PDG_t::kPiPlus);
204+
signBachelor = 1;
205+
posDauPdgCodeAbs = std::abs(PDG_t::kPiPlus);
206+
negDauPdgCodeAbs = std::abs(PDG_t::kProtonBar);
207+
}
208+
} else if (mPdgCode == PDG_t::kOmegaMinus) {
209+
if (ConfCascadeSelection.sign.value < 0) {
210+
bachelorPdgCodeAbs = std::abs(PDG_t::kKMinus);
211+
signBachelor = -1;
212+
posDauPdgCodeAbs = std::abs(PDG_t::kProton);
213+
negDauPdgCodeAbs = std::abs(PDG_t::kPiMinus);
214+
} else {
215+
mPdgCode = -1 * mPdgCode; // Omega+ has negative pdg code
216+
bachelorPdgCodeAbs = std::abs(PDG_t::kKPlus);
217+
signBachelor = 1;
218+
posDauPdgCodeAbs = std::abs(PDG_t::kPiPlus);
219+
negDauPdgCodeAbs = std::abs(PDG_t::kProtonBar);
220+
}
221+
} else {
222+
LOG(fatal) << "PDG code for Cascade has to be either Xi or Omega";
223+
}
224+
225+
mBachelorManager.template init<mode>(registry, BachelorSpecs, absCharge, signBachelor, bachelorPdgCodeAbs);
226+
mPosDauManager.template init<mode>(registry, PosDauSpecs, absCharge, signPlus, posDauPdgCodeAbs);
227+
mNegDauManager.template init<mode>(registry, NegDauSpecs, absCharge, signMinus, negDauPdgCodeAbs);
228+
187229
if constexpr (modes::isFlagSet(mode, modes::Mode::kAnalysis)) {
188230
initAnalysis(cascadeSpecs);
189231
}
190-
if constexpr (modes::isFlagSet(mode, modes::Mode::kQa)) {
191-
initQa(cascadeSpecs);
192-
}
193232
}
194233

195-
template <modes::Mode mode, typename T1, typename T2, typename T3, typename T4>
234+
template <modes::Mode mode, typename T1, typename T2, typename T3, typename T4, typename T5>
196235
void init(o2::framework::HistogramRegistry* registry,
197236
std::map<CascadeHist, std::vector<o2::framework::AxisSpec>> const& cascadeSpecs,
198-
T1 const& CascadeConfBinningQa,
237+
T1 const& ConfCascadeSelection,
238+
T2 const& ConfCascadeBinningQa,
199239
std::map<trackhistmanager::TrackHist, std::vector<o2::framework::AxisSpec>> const& BachelorSpecs,
200-
T2 const& BachelorConfBinningQa,
240+
T3 const& ConfBachelorQaBinning,
201241
std::map<trackhistmanager::TrackHist, std::vector<o2::framework::AxisSpec>> const& PosDauSpecs,
202-
T3 const& PosDauConfBinningQa,
242+
T4& ConfPosDauQaBinning,
203243
std::map<trackhistmanager::TrackHist, std::vector<o2::framework::AxisSpec>> const& NegDauSpecs,
204-
T4 const& NegDauConfBinningQa)
244+
T5& ConfNegDauQaBinning)
205245
{
206246
mHistogramRegistry = registry;
207-
mBachelorManager.template init<mode>(registry, BachelorSpecs, BachelorConfBinningQa, AbsChargeDaughters);
208-
mPosDauManager.template init<mode>(registry, PosDauSpecs, PosDauConfBinningQa, AbsChargeDaughters);
209-
mNegDauManager.template init<mode>(registry, NegDauSpecs, NegDauConfBinningQa, AbsChargeDaughters);
210-
this->enableOptionalHistograms(CascadeConfBinningQa);
247+
mPdgCode = std::abs(ConfCascadeSelection.pdgCodeAbs.value);
248+
this->enableOptionalHistograms(ConfCascadeBinningQa);
249+
250+
int bachelorPdgCodeAbs = 0;
251+
int posDauPdgCodeAbs = 0;
252+
int negDauPdgCodeAbs = 0;
253+
const int absCharge = 1;
254+
int signBachelor = 0;
255+
const int signPlus = 1;
256+
const int signMinus = -1;
257+
258+
if (mPdgCode == PDG_t::kXiMinus) {
259+
if (ConfCascadeSelection.sign.value < 0) {
260+
bachelorPdgCodeAbs = std::abs(PDG_t::kPiMinus);
261+
signBachelor = -1;
262+
posDauPdgCodeAbs = std::abs(PDG_t::kProton);
263+
negDauPdgCodeAbs = std::abs(PDG_t::kPiMinus);
264+
} else {
265+
mPdgCode = -1 * mPdgCode; // Xi+ has negative pdg code
266+
bachelorPdgCodeAbs = std::abs(PDG_t::kPiPlus);
267+
signBachelor = 1;
268+
posDauPdgCodeAbs = std::abs(PDG_t::kPiPlus);
269+
negDauPdgCodeAbs = std::abs(PDG_t::kProtonBar);
270+
}
271+
} else if (mPdgCode == PDG_t::kOmegaMinus) {
272+
if (ConfCascadeSelection.sign.value < 0) {
273+
bachelorPdgCodeAbs = std::abs(PDG_t::kKMinus);
274+
signBachelor = -1;
275+
posDauPdgCodeAbs = std::abs(PDG_t::kProton);
276+
negDauPdgCodeAbs = std::abs(PDG_t::kPiMinus);
277+
} else {
278+
mPdgCode = -1 * mPdgCode; // Omega+ has negative pdg code
279+
bachelorPdgCodeAbs = std::abs(PDG_t::kKPlus);
280+
signBachelor = 1;
281+
posDauPdgCodeAbs = std::abs(PDG_t::kPiPlus);
282+
negDauPdgCodeAbs = std::abs(PDG_t::kProtonBar);
283+
}
284+
} else {
285+
LOG(fatal) << "PDG code for Cascade has to be either Xi or Omega";
286+
}
287+
288+
mBachelorManager.template init<mode>(registry, BachelorSpecs, absCharge, signBachelor, bachelorPdgCodeAbs, ConfBachelorQaBinning);
289+
mPosDauManager.template init<mode>(registry, PosDauSpecs, absCharge, signPlus, posDauPdgCodeAbs, ConfPosDauQaBinning);
290+
mNegDauManager.template init<mode>(registry, NegDauSpecs, absCharge, signMinus, negDauPdgCodeAbs, ConfNegDauQaBinning);
291+
211292
if constexpr (modes::isFlagSet(mode, modes::Mode::kAnalysis)) {
212293
initAnalysis(cascadeSpecs);
213294
}
@@ -321,6 +402,7 @@ class CascadeHistManager
321402
trackhistmanager::TrackHistManager<bachelorPrefix> mBachelorManager;
322403
trackhistmanager::TrackHistManager<posDauPrefix> mPosDauManager;
323404
trackhistmanager::TrackHistManager<negDauPrefix> mNegDauManager;
405+
int mPdgCode = 0;
324406
};
325407
}; // namespace cascadehistmanager
326408
}; // namespace o2::analysis::femto

‎PWGCF/Femto/Core/collisionBuilder.h‎

Lines changed: 12 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -89,6 +89,7 @@ struct ConfCcdb : o2::framework::ConfigurableGroup {
8989
o2::framework::Configurable<std::string> ccdbUrl{"ccdbUrl", "http://alice-ccdb.cern.ch", "URL to ccdb"};
9090
o2::framework::Configurable<std::string> grpPath{"grpPath", "GLO/Config/GRPMagField", "Path to GRP object (Run3 -> GLO/Config/GRPMagField, Run2 -> GLO/GRP/GRP"};
9191
o2::framework::Configurable<std::string> triggerPath{"triggerPath", "EventFiltering/Zorro/", "CCDB path for trigger information"};
92+
o2::framework::Configurable<int> magFieldForced{"magFieldForced", 0, "Force value for magnetic field (kG). This will skip calls to the ccdb. Deactivate by setting value to 0"};
9293
};
9394

9495
struct ConfCollisionRctFlags : o2::framework::ConfigurableGroup {
@@ -390,6 +391,7 @@ class CollisionBuilder
390391
mUseRctFlags = true;
391392
mRctFlagsChecker.init(confRct.label.value, confRct.useZdc.value, confRct.treatLimitedAcceptanceAsBad.value);
392393
}
394+
mMagFieldForced = confCcdb.magFieldForced.value;
393395
mGrpPath = confCcdb.grpPath.value;
394396
mSubGeneratorId = confFilter.subGeneratorId.value;
395397

@@ -403,12 +405,16 @@ class CollisionBuilder
403405
{
404406
if (mRunNumber != bc.runNumber()) {
405407
mRunNumber = bc.runNumber();
406-
static o2::parameters::GRPMagField* grpo = nullptr;
407-
grpo = ccdb->template getForRun<o2::parameters::GRPMagField>(mGrpPath, mRunNumber);
408-
if (grpo == nullptr) {
409-
LOG(fatal) << "GRP object not found for Run " << mRunNumber;
408+
if (mMagFieldForced == 0) {
409+
static o2::parameters::GRPMagField* grpo = nullptr;
410+
grpo = ccdb->template getForRun<o2::parameters::GRPMagField>(mGrpPath, mRunNumber);
411+
if (grpo == nullptr) {
412+
LOG(fatal) << "GRP object not found for Run " << mRunNumber;
413+
}
414+
mMagField = static_cast<int>(grpo->getNominalL3Field()); // get magnetic field in kG
415+
} else {
416+
mMagField = mMagFieldForced;
410417
}
411-
mMagField = static_cast<int>(grpo->getNominalL3Field()); // get magnetic field in kG
412418

413419
if (mUseTrigger) {
414420
mZorro.initCCDB(ccdb.service, mRunNumber, bc.timestamp(), mTriggerNames);
@@ -534,6 +540,7 @@ class CollisionBuilder
534540
bool mUseTrigger = false;
535541
int mRunNumber = -1;
536542
std::string mGrpPath = std::string("");
543+
int mMagFieldForced = 0;
537544
int mMagField = 0;
538545
aod::rctsel::RCTFlagsChecker mRctFlagsChecker;
539546
bool mUseRctFlags = false;

‎PWGCF/Femto/Core/collisionHistManager.h‎

Lines changed: 25 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -52,8 +52,8 @@ enum ColHist {
5252
kCentVsSphericity,
5353
kMultVsSphericity,
5454
// mc
55-
kMcCent,
56-
kMcMult,
55+
kTrueCentVsCent,
56+
kTrueMultVsMult,
5757
kColHistLast
5858
};
5959

@@ -80,8 +80,8 @@ constexpr std::array<histmanager::HistInfo<ColHist>, kColHistLast> HistTable = {
8080
{kMultVsSphericity, o2::framework::kTH2F, "hMultVsSphericity", "Multiplicity vs Sphericity; Multiplicity; Sphericity"},
8181
{kCentVsSphericity, o2::framework::kTH2F, "hCentVsSphericity", "Centrality vs Sphericity; Centrality (%); Sphericity"},
8282
// mc
83-
{kMcCent, o2::framework::kTH1F, "hMcCent", "Monte Carlo Centrality; Centrality (%); Entries"},
84-
{kMcMult, o2::framework::kTH1F, "hMcMult", "Monte Carlo Multiplicity; Multiplicity; Entries"},
83+
{kTrueCentVsCent, o2::framework::kTH2F, "hTrueCentVsCent", "True centrality vs centrality; Centrality_{True} (%); Centrality (%)"},
84+
{kTrueMultVsMult, o2::framework::kTH2F, "hTrueMultVsMult", "True multiplicity vs multiplicity; Multiplicity_{True}; Multiplicity"},
8585
}};
8686

8787
#define COL_HIST_ANALYSIS_MAP(conf) \
@@ -102,9 +102,9 @@ constexpr std::array<histmanager::HistInfo<ColHist>, kColHistLast> HistTable = {
102102
{kMultVsSphericity, {confAnalysis.mult, confQa.sphericity}}, \
103103
{kCentVsSphericity, {confBinningAnalysis.cent, confQa.sphericity}},
104104

105-
#define COL_HIST_MC_QA_MAP(conf) \
106-
{kMcMult, {conf.mult}}, \
107-
{kMcCent, {conf.cent}},
105+
#define COL_HIST_MC_MAP(conf) \
106+
{kTrueMultVsMult, {conf.mult, conf.mult}}, \
107+
{kTrueCentVsCent, {conf.cent, conf.cent}},
108108

109109
template <typename T>
110110
auto makeColHistSpecMap(const T& confBinningAnalysis)
@@ -113,6 +113,14 @@ auto makeColHistSpecMap(const T& confBinningAnalysis)
113113
COL_HIST_ANALYSIS_MAP(confBinningAnalysis)};
114114
}
115115

116+
template <typename T>
117+
auto makeColMcHistSpecMap(const T& confBinningAnalysis)
118+
{
119+
return std::map<ColHist, std::vector<framework::AxisSpec>>{
120+
COL_HIST_ANALYSIS_MAP(confBinningAnalysis)
121+
COL_HIST_MC_MAP(confBinningAnalysis)};
122+
}
123+
116124
template <typename T1, typename T2>
117125
auto makeColQaHistSpecMap(const T1& confBinningAnalysis, const T2& confBinningQa)
118126
{
@@ -127,12 +135,12 @@ auto makeColMcQaHistSpecMap(const T1& confBinningAnalysis, const T2& confBinning
127135
return std::map<ColHist, std::vector<framework::AxisSpec>>{
128136
COL_HIST_ANALYSIS_MAP(confBinningAnalysis)
129137
COL_HIST_QA_MAP(confBinningAnalysis, confBinningQa)
130-
COL_HIST_MC_QA_MAP(confBinningAnalysis)};
138+
COL_HIST_MC_MAP(confBinningAnalysis)};
131139
}
132140

133141
#undef COL_HIST_ANALYSIS_MAP
134142
#undef COL_HIST_QA_MAP
135-
#undef COL_HIST_MC_QA_MAP
143+
#undef COL_HIST_MC_MAP
136144

137145
struct ConfCollisionBinning : o2::framework::ConfigurableGroup {
138146
std::string prefix = std::string("CollisionBinning");
@@ -240,8 +248,8 @@ class CollisionHistManager
240248
void initMc(std::map<ColHist, std::vector<o2::framework::AxisSpec>> const& Specs)
241249
{
242250
std::string mcDir = std::string(ColMcDir);
243-
mHistogramRegistry->add(mcDir + getHistNameV2(kMcMult, HistTable), getHistDesc(kMcMult, HistTable), getHistType(kMcMult, HistTable), {Specs.at(kMcMult)});
244-
mHistogramRegistry->add(mcDir + getHistNameV2(kMcCent, HistTable), getHistDesc(kMcCent, HistTable), getHistType(kMcCent, HistTable), {Specs.at(kMcCent)});
251+
mHistogramRegistry->add(mcDir + getHistNameV2(kTrueMultVsMult, HistTable), getHistDesc(kTrueMultVsMult, HistTable), getHistType(kTrueMultVsMult, HistTable), {Specs.at(kTrueMultVsMult)});
252+
mHistogramRegistry->add(mcDir + getHistNameV2(kTrueCentVsCent, HistTable), getHistDesc(kTrueCentVsCent, HistTable), getHistType(kTrueCentVsCent, HistTable), {Specs.at(kTrueCentVsCent)});
245253
}
246254

247255
template <typename T>
@@ -273,9 +281,12 @@ class CollisionHistManager
273281
template <typename T1, typename T2>
274282
void fillMc(T1 const& col, T2 const& /*mcCols*/)
275283
{
276-
auto genCol = col.template fMcCol_as<T2>();
277-
mHistogramRegistry->fill(HIST(ColMcDir) + HIST(getHistName(kMcMult, HistTable)), genCol.multMc());
278-
mHistogramRegistry->fill(HIST(ColMcDir) + HIST(getHistName(kMcCent, HistTable)), genCol.centMc());
284+
if (!col.has_fMcCol()) {
285+
return;
286+
}
287+
auto mcCol = col.template fMcCol_as<T2>();
288+
mHistogramRegistry->fill(HIST(ColMcDir) + HIST(getHistName(kTrueMultVsMult, HistTable)), mcCol.mult(), col.mult());
289+
mHistogramRegistry->fill(HIST(ColMcDir) + HIST(getHistName(kTrueCentVsCent, HistTable)), mcCol.cent(), col.cent());
279290
}
280291

281292
o2::framework::HistogramRegistry* mHistogramRegistry = nullptr;

‎PWGCF/Femto/Core/kinkBuilder.h‎

Lines changed: 12 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -94,16 +94,16 @@ struct ConfSigmaPlusBits : o2::framework::ConfigurableGroup {
9494
#undef KINK_DEFAULT_BITS
9595

9696
// base selection for analysis task for kinks
97-
#define KINK_DEFAULT_SELECTIONS(defaultMassMin, defaultMassMax, defaultPdgCode) \
98-
o2::framework::Configurable<int> pdgCode{"pdgCode", defaultPdgCode, "Kink PDG code"}; \
99-
o2::framework::Configurable<float> ptMin{"ptMin", 0.f, "Minimum pT"}; \
100-
o2::framework::Configurable<float> ptMax{"ptMax", 999.f, "Maximum pT"}; \
101-
o2::framework::Configurable<float> etaMin{"etaMin", -10.f, "Minimum eta"}; \
102-
o2::framework::Configurable<float> etaMax{"etaMax", 10.f, "Maximum eta"}; \
103-
o2::framework::Configurable<float> phiMin{"phiMin", 0.f, "Minimum phi"}; \
104-
o2::framework::Configurable<float> phiMax{"phiMax", 1.f * o2::constants::math::TwoPI, "Maximum phi"}; \
105-
o2::framework::Configurable<float> massMin{"massMin", defaultMassMin, "Minimum invariant mass for Sigma"}; \
106-
o2::framework::Configurable<float> massMax{"massMax", defaultMassMax, "Maximum invariant mass for Sigma"}; \
97+
#define KINK_DEFAULT_SELECTIONS(defaultMassMin, defaultMassMax, defaultPdgCode) \
98+
o2::framework::Configurable<int> pdgCodeAbs{"pdgCodeAbs", defaultPdgCode, "PDG code. Select antipartilce via sign"}; \
99+
o2::framework::Configurable<float> ptMin{"ptMin", 0.f, "Minimum pT"}; \
100+
o2::framework::Configurable<float> ptMax{"ptMax", 999.f, "Maximum pT"}; \
101+
o2::framework::Configurable<float> etaMin{"etaMin", -10.f, "Minimum eta"}; \
102+
o2::framework::Configurable<float> etaMax{"etaMax", 10.f, "Maximum eta"}; \
103+
o2::framework::Configurable<float> phiMin{"phiMin", 0.f, "Minimum phi"}; \
104+
o2::framework::Configurable<float> phiMax{"phiMax", 1.f * o2::constants::math::TwoPI, "Maximum phi"}; \
105+
o2::framework::Configurable<float> massMin{"massMin", defaultMassMin, "Minimum invariant mass for Sigma"}; \
106+
o2::framework::Configurable<float> massMax{"massMax", defaultMassMax, "Maximum invariant mass for Sigma"}; \
107107
o2::framework::Configurable<o2::aod::femtodatatypes::KinkMaskType> mask{"mask", 0x0, "Bitmask for kink selection"};
108108

109109
// base selection for analysis task for sigmas
@@ -498,11 +498,11 @@ class KinkBuilder
498498

499499
if constexpr (modes::isEqual(kinkType, modes::Kink::kSigma)) {
500500
fillSigma(collisionProducts, kinkProducts, kink, daughterIndex);
501-
mcBuilder.template fillMcSigmaWithLabel<system>(col, mcCols, kink, daughter, mcParticles, mcProducts);
501+
mcBuilder.template fillMcSigmaWithLabel<system>(col, mcCols, daughter, mcParticles, mcProducts);
502502
}
503503
if constexpr (modes::isEqual(kinkType, modes::Kink::kSigmaPlus)) {
504504
fillSigmaPlus(collisionProducts, kinkProducts, kink, daughterIndex);
505-
mcBuilder.template fillMcSigmaPlusWithLabel<system>(col, mcCols, kink, daughter, mcParticles, mcProducts);
505+
mcBuilder.template fillMcSigmaPlusWithLabel<system>(col, mcCols, daughter, mcParticles, mcProducts);
506506
}
507507
}
508508
}

0 commit comments

Comments
 (0)