Skip to content

Commit 4653ef9

Browse files
yashpatleyvictor-gonzalez
yashpatley
and
victor-gonzalez
authored
[PWGCF] Update lambdaR2Correlation.cxx (AliceO2Group#11260)
Co-authored-by: Victor Gonzalez <[email protected]>
1 parent 01dd039 commit 4653ef9

File tree

1 file changed

+83
-119
lines changed

1 file changed

+83
-119
lines changed

‎PWGCF/TwoParticleCorrelations/Tasks/lambdaR2Correlation.cxx‎

Lines changed: 83 additions & 119 deletions
Original file line numberDiff line numberDiff line change
@@ -148,7 +148,6 @@ enum TrackLabels {
148148
kPassV0DauTrackSel,
149149
kPassV0KinCuts,
150150
kPassV0TopoSel,
151-
kAllV0SelPassed,
152151
kAllSelPassed,
153152
kNotPrimaryLambda,
154153
kNotSecondaryLambda,
@@ -160,6 +159,10 @@ enum TrackLabels {
160159
kEffCorrPtRap,
161160
kEffCorrPtRapPhi,
162161
kNoEffCorr,
162+
kPFCorrPt,
163+
kPFCorrPtRap,
164+
kPFCorrPtRapPhi,
165+
kNoPFCorr,
163166
kGenTotAccLambda,
164167
kGenLambdaNoDau,
165168
kGenLambdaToPrPi
@@ -260,18 +263,6 @@ struct LambdaTableProducer {
260263
Configurable<bool> cV0TypeSelFlag{"cV0TypeSelFlag", false, "V0 Type Selection Flag"};
261264
Configurable<int> cV0TypeSelection{"cV0TypeSelection", 1, "V0 Type Selection"};
262265

263-
// Cascade V0
264-
Configurable<bool> cRemoveCascLambda{"cRemoveCascLambda", false, "Remove Cascade V0s"};
265-
Configurable<float> cMinCascDcaPosToPV{"cMinCascDcaPosToPV", 0.03, "DCA Casc Pos To PV"};
266-
Configurable<float> cMinCascDcaNegToPV{"cMinCascDcaNegToPV", 0.03, "DCA Casc Neg To PV"};
267-
Configurable<float> cMinCascDcaBachToPV{"cMinCascDcaBachToPV", 0.03, "DCA Casc Bach To PV"};
268-
Configurable<float> cMaxCascDcaDaughters{"cMaxCascDcaDaughters", 1.0, "DCA Casc Dau at Casc Decay Vtx"};
269-
Configurable<float> cMinCascRadius{"cMinCascRadius", 0.4, "Casc Decay Radius"};
270-
Configurable<float> cMinCascV0Radius{"cMinCascV0Radius", 0.9, "Casc V0 Decay Radius"};
271-
Configurable<float> cMinCascCosPA{"cMinCascCosPA", 0.99, "Casc CosThetaPA"};
272-
Configurable<float> cMinCascV0CosPA{"cMinCascV0CosPA", 0.99, "Casc V0 CosThetaPA"};
273-
Configurable<float> cCascMassWindow{"cCascMassWindow", 0.005, "Casc Mass Window"};
274-
275266
// V0s MC
276267
Configurable<bool> cHasMcFlag{"cHasMcFlag", true, "Has Mc Tag"};
277268
Configurable<bool> cSelectTrueLambda{"cSelectTrueLambda", true, "Select True Lambda"};
@@ -289,9 +280,11 @@ struct LambdaTableProducer {
289280
Configurable<bool> cDoTrackMcMatching{"cDoTrackMcMatching", false, "Do Track Mc Matching Flag"};
290281

291282
// Efficiency Correction
292-
Configurable<bool> cCorrectionFlag{"cCorrectionFlag", false, "Efficiency Correction Flag"};
293-
Configurable<int> cCorrFactHist{"cCorrFactHist", 0, "Correction Factor Histogram"};
294-
Configurable<bool> cDoEtaCorr{"cDoEtaCorr", false, "Do Eta Corr"};
283+
Configurable<bool> cCorrectionFlag{"cCorrectionFlag", false, "Correction Flag"};
284+
Configurable<bool> cGetEffFact{"cGetEffFact", false, "Get Efficiency Factor Flag"};
285+
Configurable<bool> cGetPrimFrac{"cGetPrimFrac", false, "Get Primary Fraction Flag"};
286+
Configurable<int> cCorrFactHist{"cCorrFactHist", 0, "Efficiency Factor Histogram"};
287+
Configurable<int> cPrimFracHist{"cPrimFracHist", 0, "Primary Fraction Histogram"};
295288

296289
// CCDB
297290
Configurable<std::string> cUrlCCDB{"cUrlCCDB", "http://ccdb-test.cern.ch:8080", "url of ccdb"};
@@ -310,6 +303,13 @@ struct LambdaTableProducer {
310303
{"hEffVsPtYPhiLambda", "hEffVsPtYPhiAntiLambda"},
311304
{"hEffVsPtEtaPhiLambda", "hEffVsPtEtaPhiAntiLambda"}};
312305

306+
// initialize corr_factor objects
307+
std::vector<std::vector<std::string>> vPrimFracStrings = {{"hPrimFracVsPtLambda", "hPrimFracVsPtAntiLambda"},
308+
{"hPrimFracVsPtYLambda", "hPrimFracVsPtYAntiLambda"},
309+
{"hPrimFracVsPtEtaLambda", "hPrimFracVsPtEtaAntiLambda"},
310+
{"hPrimFracVsPtYPhiLambda", "hPrimFracVsPtYPhiAntiLambda"},
311+
{"hPrimFracVsPtEtaPhiLambda", "hPrimFracVsPtEtaPhiAntiLambda"}};
312+
313313
// Initialize Global Variables
314314
float cent = 0.;
315315

@@ -321,7 +321,7 @@ struct LambdaTableProducer {
321321

322322
// initialize axis specifications
323323
const AxisSpec axisCols(5, 0.5, 5.5, "");
324-
const AxisSpec axisTrks(25, 0.5, 25.5, "");
324+
const AxisSpec axisTrks(30, 0.5, 30.5, "");
325325
const AxisSpec axisCent(100, 0, 100, "FT0M (%)");
326326
const AxisSpec axisMult(10, 0, 10, "N_{#Lambda}");
327327
const AxisSpec axisVz(220, -11, 11, "V_{z} (cm)");
@@ -349,8 +349,6 @@ struct LambdaTableProducer {
349349
const AxisSpec axisNsigma(401, -10.025, 10.025, {"n#sigma"});
350350
const AxisSpec axisdEdx(360, 20, 200, "#frac{dE}{dx}");
351351

352-
const AxisSpec axisCascMass{100, 1.28, 1.38, "M_{#Xi}"};
353-
354352
// Create Histograms.
355353
// Event histograms
356354
histos.add("Events/h1f_collisions_info", "# of Collisions", kTH1F, {axisCols});
@@ -395,9 +393,6 @@ struct LambdaTableProducer {
395393
histos.add("QA/Lambda/h2f_pos_prong_tpc_nsigma_pi_vs_p", "TPC n#sigma Pos Prong", kTH2F, {axisMomPID, axisNsigma});
396394
histos.add("QA/Lambda/h2f_neg_prong_tpc_nsigma_pi_vs_p", "TPC n#sigma Neg Prong", kTH2F, {axisMomPID, axisNsigma});
397395

398-
histos.add("QA/Casc/hMassBeforeCuts", "Xi Mass", kTH1F, {axisCascMass});
399-
histos.add("QA/Casc/hMass", "Xi Mass", kTH1F, {axisCascMass});
400-
401396
// Kinematic Histograms
402397
histos.add("McRec/Lambda/hPt", "Transverse Momentum", kTH1F, {axisV0Pt});
403398
histos.add("McRec/Lambda/hEta", "Pseudorapidity", kTH1F, {axisV0Eta});
@@ -464,12 +459,15 @@ struct LambdaTableProducer {
464459
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kPassV0DauTrackSel, "kPassV0DauTrackSel");
465460
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kPassV0KinCuts, "kPassV0KinCuts");
466461
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kPassV0TopoSel, "kPassV0TopoSel");
467-
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kAllV0SelPassed, "kAllV0SelPassed");
468462
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kAllSelPassed, "kAllSelPassed");
469463
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kEffCorrPt, "kEffCorrPt");
470464
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kEffCorrPtRap, "kEffCorrPtRap");
471465
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kEffCorrPtRapPhi, "kEffCorrPtRapPhi");
472466
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kNoEffCorr, "kNoEffCorr");
467+
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kPFCorrPt, "kPFCorrPt");
468+
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kPFCorrPtRap, "kPFCorrPtRap");
469+
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kPFCorrPtRapPhi, "kPFCorrPtRapPhi");
470+
histos.get<TH1>(HIST("Tracks/h1f_tracks_info"))->GetXaxis()->SetBinLabel(TrackLabels::kNoPFCorr, "kNoPFCorr");
473471
}
474472

475473
template <RunType run, typename C>
@@ -809,57 +807,6 @@ struct LambdaTableProducer {
809807
return true;
810808
}
811809

812-
// Remove Lambda from Xi ---> Lambda Pi
813-
template <typename C, typename V, typename X, typename T>
814-
bool removeCascLambda(C const& collision, V const& v0, X const& cascs, T const&)
815-
{
816-
// Loop over cascades
817-
for (auto const& casc : cascs) {
818-
// cascade daughters
819-
auto cascposdau = casc.template posTrack_as<T>();
820-
auto cascnegdau = casc.template negTrack_as<T>();
821-
auto bachelor = casc.template bachelor_as<T>();
822-
823-
// V0 daughters
824-
auto v0posdau = v0.template posTrack_as<T>();
825-
auto v0negdau = v0.template negTrack_as<T>();
826-
827-
// topological selection
828-
if (std::abs(casc.dcapostopv()) <= cMinCascDcaPosToPV ||
829-
std::abs(casc.dcanegtopv()) <= cMinCascDcaNegToPV ||
830-
std::abs(casc.dcabachtopv()) <= cMinCascDcaBachToPV ||
831-
casc.dcaV0daughters() >= cMaxV0DcaDaughters ||
832-
casc.dcacascdaughters() >= cMaxCascDcaDaughters ||
833-
casc.cascradius() <= cMinCascRadius ||
834-
casc.v0radius() <= cMinCascV0Radius ||
835-
casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()) <= cMinCascCosPA ||
836-
casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()) <= cMinCascV0CosPA ||
837-
std::abs(cascposdau.eta()) >= cTrackEtaCut ||
838-
std::abs(cascnegdau.eta()) >= cTrackEtaCut ||
839-
std::abs(bachelor.eta()) >= cTrackEtaCut) {
840-
continue;
841-
}
842-
843-
histos.fill(HIST("QA/Casc/hMassBeforeCuts"), casc.mXi());
844-
845-
// apply mass window selection to Xi
846-
if (std::abs(casc.mXi() - MassXiMinus) >= cCascMassWindow) {
847-
continue;
848-
}
849-
850-
histos.fill(HIST("QA/Casc/hMass"), casc.mXi());
851-
852-
// Check if the daughters of cascades match with v0 daughters
853-
// The "OR" Logic is Temperory, we can further see the effect of "AND" Logic
854-
if (v0posdau.index() == cascposdau.index() || v0negdau.index() == cascnegdau.index()) {
855-
return true;
856-
}
857-
}
858-
859-
// No Lambda from Cascade
860-
return false;
861-
}
862-
863810
template <typename T>
864811
bool getMcMatch(T const& vrec, T const& vgen)
865812
{
@@ -899,7 +846,7 @@ struct LambdaTableProducer {
899846
template <ParticleType part, typename V>
900847
float getCorrectionFactors(V const& v0)
901848
{
902-
// Check for efficiency correction flag and Rec/Gen Data
849+
// Check for efficiency correction flag
903850
if (!cCorrectionFlag) {
904851
return 1.;
905852
}
@@ -913,29 +860,56 @@ struct LambdaTableProducer {
913860
return 1.;
914861
}
915862

916-
// get ccdb object
917-
TObject* obj = reinterpret_cast<TObject*>(ccdbObj->FindObject(Form("%s", vCorrFactStrings[cCorrFactHist][part].c_str())));
918-
TH1F* hist = reinterpret_cast<TH1F*>(obj->Clone());
919-
float retVal = 0.;
920-
float rap = (cDoEtaCorr) ? v0.eta() : v0.yLambda();
921-
922-
if (hist->GetDimension() == OneDimCorr) {
923-
histos.fill(HIST("Tracks/h1f_tracks_info"), kEffCorrPt);
924-
retVal = hist->GetBinContent(hist->FindBin(v0.pt()));
925-
} else if (hist->GetDimension() == TwoDimCorr) {
926-
histos.fill(HIST("Tracks/h1f_tracks_info"), kEffCorrPtRap);
927-
retVal = hist->GetBinContent(hist->FindBin(v0.pt(), rap));
928-
} else if (hist->GetDimension() == ThreeDimCorr) {
929-
histos.fill(HIST("Tracks/h1f_tracks_info"), kEffCorrPtRapPhi);
930-
retVal = hist->GetBinContent(hist->FindBin(v0.pt(), rap, v0.phi()));
931-
} else {
932-
histos.fill(HIST("Tracks/h1f_tracks_info"), kNoEffCorr);
933-
LOGF(warning, "CCDB OBJECT IS NOT A HISTOGRAM !!!");
934-
retVal = 1.;
863+
// initialize efficiency factor and primary fraction values
864+
float effFact = 1., primFrac = 1.;
865+
float rap = (cDoEtaAnalysis) ? v0.eta() : v0.yLambda();
866+
867+
// Get Efficiency Factor
868+
if (cGetEffFact) {
869+
TObject* objEff = reinterpret_cast<TObject*>(ccdbObj->FindObject(Form("%s", vCorrFactStrings[cCorrFactHist][part].c_str())));
870+
TH1F* histEff = reinterpret_cast<TH1F*>(objEff->Clone());
871+
if (histEff->GetDimension() == OneDimCorr) {
872+
histos.fill(HIST("Tracks/h1f_tracks_info"), kEffCorrPt);
873+
effFact = histEff->GetBinContent(histEff->FindBin(v0.pt()));
874+
} else if (histEff->GetDimension() == TwoDimCorr) {
875+
histos.fill(HIST("Tracks/h1f_tracks_info"), kEffCorrPtRap);
876+
effFact = histEff->GetBinContent(histEff->FindBin(v0.pt(), rap));
877+
} else if (histEff->GetDimension() == ThreeDimCorr) {
878+
histos.fill(HIST("Tracks/h1f_tracks_info"), kEffCorrPtRapPhi);
879+
effFact = histEff->GetBinContent(histEff->FindBin(v0.pt(), rap, v0.phi()));
880+
} else {
881+
histos.fill(HIST("Tracks/h1f_tracks_info"), kNoEffCorr);
882+
LOGF(warning, "CCDB OBJECT IS NOT A HISTOGRAM !!!");
883+
effFact = 1.;
884+
}
885+
delete objEff;
886+
delete histEff;
887+
}
888+
889+
// Get Primary Fraction
890+
// (The dimension of this could be different than efficiency because of large errors !!!)
891+
if (cGetPrimFrac) {
892+
TObject* objPrm = reinterpret_cast<TObject*>(ccdbObj->FindObject(Form("%s", vPrimFracStrings[cPrimFracHist][part].c_str())));
893+
TH1F* histPrm = reinterpret_cast<TH1F*>(objPrm->Clone());
894+
if (histPrm->GetDimension() == OneDimCorr) {
895+
histos.fill(HIST("Tracks/h1f_tracks_info"), kPFCorrPt);
896+
primFrac = histPrm->GetBinContent(histPrm->FindBin(v0.pt()));
897+
} else if (histPrm->GetDimension() == TwoDimCorr) {
898+
histos.fill(HIST("Tracks/h1f_tracks_info"), kPFCorrPtRap);
899+
primFrac = histPrm->GetBinContent(histPrm->FindBin(v0.pt(), rap));
900+
} else if (histPrm->GetDimension() == ThreeDimCorr) {
901+
histos.fill(HIST("Tracks/h1f_tracks_info"), kPFCorrPtRapPhi);
902+
primFrac = histPrm->GetBinContent(histPrm->FindBin(v0.pt(), rap, v0.phi()));
903+
} else {
904+
histos.fill(HIST("Tracks/h1f_tracks_info"), kNoPFCorr);
905+
LOGF(warning, "CCDB OBJECT IS NOT A HISTOGRAM !!!");
906+
primFrac = 1.;
907+
}
908+
delete objPrm;
909+
delete histPrm;
935910
}
936911

937-
delete hist;
938-
return retVal;
912+
return primFrac / effFact;
939913
}
940914

941915
template <ParticleType part, typename V>
@@ -1013,8 +987,8 @@ struct LambdaTableProducer {
1013987
}
1014988

1015989
// Reconstructed Level Tables
1016-
template <RunType run, DMCType dmc, typename C, typename V, typename T, typename X>
1017-
void fillLambdaRecoTables(C const& collision, V const& v0tracks, T const& tracks, X const& cascs)
990+
template <RunType run, DMCType dmc, typename C, typename V, typename T>
991+
void fillLambdaRecoTables(C const& collision, V const& v0tracks, T const& tracks)
1018992
{
1019993
// Total Collisions
1020994
histos.fill(HIST("Events/h1f_collisions_info"), kTotCol);
@@ -1036,7 +1010,6 @@ struct LambdaTableProducer {
10361010
ParticleType v0Type = kLambda;
10371011
float mass = 0., corr_fact = 1.;
10381012

1039-
// Loop Over V0s
10401013
for (auto const& v0 : v0tracks) {
10411014
// check for corresponding MCGen Particle
10421015
if constexpr (dmc == kMC) {
@@ -1059,13 +1032,6 @@ struct LambdaTableProducer {
10591032
continue;
10601033
}
10611034

1062-
histos.fill(HIST("Tracks/h1f_tracks_info"), kAllV0SelPassed);
1063-
1064-
// Remove Lambda V0 coming from Xi -> Lambda Pi
1065-
if (cRemoveCascLambda && removeCascLambda(collision, v0, cascs, tracks)) {
1066-
continue;
1067-
}
1068-
10691035
histos.fill(HIST("Tracks/h1f_tracks_info"), kAllSelPassed);
10701036

10711037
// we have v0 as lambda
@@ -1223,8 +1189,8 @@ struct LambdaTableProducer {
12231189
}
12241190
}
12251191

1226-
template <RunType run, DMCType dmc, typename M, typename C, typename V, typename T, typename X, typename P>
1227-
void analyzeMcRecoGen(M const& mcCollision, C const& collisions, V const& V0s, T const& tracks, X const& cascs, P const& mcParticles)
1192+
template <RunType run, DMCType dmc, typename M, typename C, typename V, typename T, typename P>
1193+
void analyzeMcRecoGen(M const& mcCollision, C const& collisions, V const& V0s, T const& tracks, P const& mcParticles)
12281194
{
12291195
// Number of Rec Collisions Associated to the McGen Collision
12301196
int nRecCols = collisions.size();
@@ -1251,7 +1217,7 @@ struct LambdaTableProducer {
12511217
histos.fill(HIST("McGen/h2f_collision_posZ"), mcCollision.posZ(), collisions.begin().posZ());
12521218
histos.fill(HIST("McGen/h2f_collision_cent"), mcCollision.centFT0M(), cent);
12531219
auto v0Tracks = V0s.sliceBy(perCollision, collisions.begin().globalIndex());
1254-
fillLambdaRecoTables<run, dmc>(collisions.begin(), v0Tracks, tracks, cascs);
1220+
fillLambdaRecoTables<run, dmc>(collisions.begin(), v0Tracks, tracks);
12551221
fillLambdaMcGenTables<run>(mcCollision, mcParticles);
12561222
}
12571223

@@ -1260,42 +1226,40 @@ struct LambdaTableProducer {
12601226

12611227
using CollisionsRun3 = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms>;
12621228
using CollisionsRun2 = soa::Join<aod::Collisions, aod::EvSels, aod::CentRun2V0Ms>;
1263-
using CascsTracks = aod::CascDataExt;
12641229
using Tracks = soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA, aod::pidTPCPi, aod::pidTPCPr, aod::pidTOFPi, aod::pidTOFPr>;
12651230
using McV0Tracks = soa::Join<aod::V0Datas, aod::McV0Labels>;
1266-
using McCascTracks = soa::Join<aod::CascDataExt, aod::McCascLabels>;
12671231
using TracksMC = soa::Join<Tracks, aod::McTrackLabels>;
12681232

1269-
void processDataRun3(CollisionsRun3::iterator const& collision, aod::V0Datas const& V0s, Tracks const& tracks, CascsTracks const& cascs)
1233+
void processDataRun3(CollisionsRun3::iterator const& collision, aod::V0Datas const& V0s, Tracks const& tracks)
12701234
{
1271-
fillLambdaRecoTables<kRun3, kData>(collision, V0s, tracks, cascs);
1235+
fillLambdaRecoTables<kRun3, kData>(collision, V0s, tracks);
12721236
}
12731237

12741238
PROCESS_SWITCH(LambdaTableProducer, processDataRun3, "Process for Run3 DATA", true);
12751239

1276-
void processDataRun2(CollisionsRun2::iterator const& collision, aod::V0Datas const& V0s, Tracks const& tracks, CascsTracks const& cascs)
1240+
void processDataRun2(CollisionsRun2::iterator const& collision, aod::V0Datas const& V0s, Tracks const& tracks)
12771241
{
1278-
fillLambdaRecoTables<kRun2, kData>(collision, V0s, tracks, cascs);
1242+
fillLambdaRecoTables<kRun2, kData>(collision, V0s, tracks);
12791243
}
12801244

12811245
PROCESS_SWITCH(LambdaTableProducer, processDataRun2, "Process for Run2 DATA", false);
12821246

12831247
void processMCRun3(soa::Join<aod::McCollisions, aod::McCentFT0Ms>::iterator const& mcCollision,
12841248
soa::SmallGroups<soa::Join<CollisionsRun3, aod::McCollisionLabels>> const& collisions,
1285-
McV0Tracks const& V0s, TracksMC const& tracks, McCascTracks const& cascs,
1249+
McV0Tracks const& V0s, TracksMC const& tracks,
12861250
aod::McParticles const& mcParticles)
12871251
{
1288-
analyzeMcRecoGen<kRun3, kMC>(mcCollision, collisions, V0s, tracks, cascs, mcParticles);
1252+
analyzeMcRecoGen<kRun3, kMC>(mcCollision, collisions, V0s, tracks, mcParticles);
12891253
}
12901254

12911255
PROCESS_SWITCH(LambdaTableProducer, processMCRun3, "Process for Run3 MC RecoGen", false);
12921256

12931257
void processMCRun2(soa::Join<aod::McCollisions, aod::McCentFT0Ms>::iterator const& mcCollision,
12941258
soa::SmallGroups<soa::Join<CollisionsRun2, aod::McCollisionLabels>> const& collisions,
1295-
McV0Tracks const& V0s, TracksMC const& tracks, McCascTracks const& cascs,
1259+
McV0Tracks const& V0s, TracksMC const& tracks,
12961260
aod::McParticles const& mcParticles)
12971261
{
1298-
analyzeMcRecoGen<kRun2, kMC>(mcCollision, collisions, V0s, tracks, cascs, mcParticles);
1262+
analyzeMcRecoGen<kRun2, kMC>(mcCollision, collisions, V0s, tracks, mcParticles);
12991263
}
13001264

13011265
PROCESS_SWITCH(LambdaTableProducer, processMCRun2, "Process for Run2 MC RecoGen", false);

0 commit comments

Comments
 (0)