Skip to content

Commit bbb7cc2

Browse files
gianniliverarogianniliveraro
committed
Using PresliceUnsorted instead of Preslice + more QA histos + clean up
1 parent 9bbff51 commit bbb7cc2

File tree

1 file changed

+63
-69
lines changed

1 file changed

+63
-69
lines changed

‎PWGLF/TableProducer/Strangeness/sigma0builder.cxx‎

Lines changed: 63 additions & 69 deletions
Original file line numberDiff line numberDiff line change
@@ -71,8 +71,8 @@ struct sigma0builder {
7171
Produces<aod::SigmaMCCores> sigma0mccores;
7272

7373
// For manual sliceBy
74-
Preslice<V0DerivedMCDatas> perCollisionMCDerived = o2::aod::v0data::straCollisionId;
75-
Preslice<V0StandardDerivedDatas> perCollisionSTDDerived = o2::aod::v0data::straCollisionId;
74+
PresliceUnsorted<V0DerivedMCDatas> perCollisionMCDerived = o2::aod::v0data::straCollisionId;
75+
PresliceUnsorted<V0StandardDerivedDatas> perCollisionSTDDerived = o2::aod::v0data::straCollisionId;
7676
PresliceUnsorted<soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraCollLabels>> perMcCollision = aod::v0data::straMCCollisionId;
7777

7878
// pack track quality but separte also afterburner
@@ -256,6 +256,7 @@ struct sigma0builder {
256256
histos.add("PhotonSel/hPhotonDCAPosToPV", "hPhotonDCAPosToPV", kTH1F, {axisDCAtoPV});
257257
histos.add("PhotonSel/hPhotonDCADau", "hPhotonDCADau", kTH1F, {axisDCAdau});
258258
histos.add("PhotonSel/hPhotonRadius", "hPhotonRadius", kTH1F, {axisRadius});
259+
histos.add("PhotonSel/h3dPhotonMass", "h3dPhotonMass", kTH3D, {axisCentrality, axisPt, axisPhotonMass});
259260

260261
histos.add("LambdaSel/hSelectionStatistics", "hSelectionStatistics", kTH1D, {axisCandSel});
261262
histos.get<TH1>(HIST("LambdaSel/hSelectionStatistics"))->GetXaxis()->SetBinLabel(1, "No Sel");
@@ -274,6 +275,8 @@ struct sigma0builder {
274275
histos.add("LambdaSel/hLambdaDCAPosToPV", "hLambdaDCAPosToPV", kTH1F, {axisDCAtoPV});
275276
histos.add("LambdaSel/hLambdaDCADau", "hLambdaDCADau", kTH1F, {axisDCAdau});
276277
histos.add("LambdaSel/hLambdaRadius", "hLambdaRadius", kTH1F, {axisRadius});
278+
histos.add("LambdaSel/h3dLambdaMass", "h3dLambdaMass", kTH3D, {axisCentrality, axisPt, axisLambdaMass});
279+
histos.add("LambdaSel/h3dALambdaMass", "h3dALambdaMass", kTH3D, {axisCentrality, axisPt, axisLambdaMass});
277280

278281
histos.add("SigmaSel/hSelectionStatistics", "hSelectionStatistics", kTH1D, {axisCandSel});
279282
histos.get<TH1>(HIST("SigmaSel/hSelectionStatistics"))->GetXaxis()->SetBinLabel(1, "No Sel");
@@ -286,7 +289,7 @@ struct sigma0builder {
286289
histos.add("SigmaSel/hSigmaMassWindow", "hSigmaMassWindow", kTH1F, {{200, -0.09f, 0.11f}});
287290
histos.add("SigmaSel/hSigmaY", "hSigmaY", kTH1F, {axisRapidity});
288291
histos.add("SigmaSel/hSigmaMassSelected", "hSigmaMassSelected", kTH1F, {axisSigmaMass});
289-
histos.add("SigmaSel/h3dMassSigma0AfterSel", "h3dMassSigma0AfterSel", kTH3F, {axisCentrality, axisPt, axisSigmaMass});
292+
histos.add("SigmaSel/h3dMassSigma0AfterSel", "h3dMassSigma0AfterSel", kTH3D, {axisCentrality, axisPt, axisSigmaMass});
290293

291294
if (fillQAhistos) {
292295
histos.add("GeneralQA/h2dMassGammaVsK0S", "h2dMassGammaVsK0S", kTH2D, {axisPhotonMass, axisK0SMass});
@@ -347,10 +350,10 @@ struct sigma0builder {
347350

348351
// For Pi0 QA
349352
if (doPi0QA) {
350-
histos.add("Pi0QA/h2dPtVsMassPi0BeforeSel_MCAssoc", "h2dPtVsMassPi0BeforeSel_MCAssoc", kTH2D, {axisPt, axisPi0Mass});
351-
histos.add("Pi0QA/h2dPtVsMassPi0AfterSel_MCAssoc", "h2dPtVsMassPi0AfterSel_MCAssoc", kTH2D, {axisPt, axisPi0Mass});
352-
histos.add("Pi0QA/h2dPtVsMassPi0BeforeSel_Candidates", "h2dPtVsMassPi0BeforeSel_Candidates", kTH2D, {axisPt, axisPi0Mass});
353-
histos.add("Pi0QA/h2dPtVsMassPi0AfterSel_Candidates", "h2dPtVsMassPi0AfterSel_Candidates", kTH2D, {axisPt, axisPi0Mass});
353+
histos.add("Pi0QA/h3dMassPi0BeforeSel_MCAssoc", "h3dMassPi0BeforeSel_MCAssoc", kTH3D, {axisCentrality, axisPt, axisPi0Mass});
354+
histos.add("Pi0QA/h3dMassPi0AfterSel_MCAssoc", "h3dMassPi0AfterSel_MCAssoc", kTH3D, {axisCentrality, axisPt, axisPi0Mass});
355+
histos.add("Pi0QA/h3dMassPi0BeforeSel_Candidates", "h3dMassPi0BeforeSel_Candidates", kTH3D, {axisCentrality, axisPt, axisPi0Mass});
356+
histos.add("Pi0QA/h3dMassPi0AfterSel_Candidates", "h3dMassPi0AfterSel_Candidates", kTH3D, {axisCentrality, axisPt, axisPi0Mass});
354357
}
355358

356359
if (doprocessGeneratedRun3) {
@@ -500,6 +503,7 @@ struct sigma0builder {
500503
if (fillHists)
501504
histos.fill(HIST("hEventSelection"), 19 /* Above max IR */);
502505

506+
histos.fill(HIST("hEventCentrality"), collision.centFT0C());
503507
return true;
504508
}
505509

@@ -731,8 +735,8 @@ struct sigma0builder {
731735
}
732736
}
733737

734-
template <typename TV0Object>
735-
void runPi0QA(TV0Object const& gamma1, TV0Object const& gamma2)
738+
template <typename TV0Object, typename TCollision>
739+
void runPi0QA(TV0Object const& gamma1, TV0Object const& gamma2, TCollision collision)
736740
{
737741
// Check if both V0s are made of the same tracks
738742
if (gamma1.posTrackExtraId() == gamma2.posTrackExtraId() ||
@@ -762,12 +766,12 @@ struct sigma0builder {
762766
gamma1MC.pdgCodeMother() == 111 && gamma2MC.pdgCodeMother() == 111 &&
763767
gamma1.motherMCPartId() == gamma2.motherMCPartId()) {
764768
fIsPi0 = true;
765-
histos.fill(HIST("Pi0QA/h2dPtVsMassPi0BeforeSel_MCAssoc"), pi0Pt, pi0Mass);
769+
histos.fill(HIST("Pi0QA/h3dMassPi0BeforeSel_MCAssoc"), collision.centFT0C(), pi0Pt, pi0Mass);
766770
}
767771
}
768772
}
769773

770-
histos.fill(HIST("Pi0QA/h2dPtVsMassPi0BeforeSel_Candidates"), pi0Pt, pi0Mass);
774+
histos.fill(HIST("Pi0QA/h3dMassPi0BeforeSel_Candidates"), collision.centFT0C(), pi0Pt, pi0Mass);
771775

772776
// Photon-specific selections
773777
auto posTrackGamma1 = gamma1.template posTrackExtra_as<dauTracks>();
@@ -776,8 +780,8 @@ struct sigma0builder {
776780
auto negTrackGamma2 = gamma2.template negTrackExtra_as<dauTracks>();
777781

778782
// Gamma1 Selection
779-
bool passedTPCGamma1 = (posTrackGamma1.tpcNSigmaEl() == -999.f || TMath::Abs(posTrackGamma1.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas) &&
780-
(negTrackGamma1.tpcNSigmaEl() == -999.f || TMath::Abs(negTrackGamma1.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas);
783+
bool passedTPCGamma1 = (TMath::Abs(posTrackGamma1.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas) ||
784+
(TMath::Abs(negTrackGamma1.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas);
781785

782786
if (TMath::Abs(gamma1.mGamma()) > Pi0PhotonMaxMass ||
783787
gamma1.qtarm() >= Pi0PhotonMaxQt ||
@@ -797,8 +801,8 @@ struct sigma0builder {
797801
}
798802

799803
// Gamma2 Selection
800-
bool passedTPCGamma2 = (posTrackGamma2.tpcNSigmaEl() == -999.f || TMath::Abs(posTrackGamma2.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas) &&
801-
(negTrackGamma2.tpcNSigmaEl() == -999.f || TMath::Abs(negTrackGamma2.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas);
804+
bool passedTPCGamma2 = (TMath::Abs(posTrackGamma2.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas) ||
805+
(TMath::Abs(negTrackGamma2.tpcNSigmaEl()) < Pi0PhotonMaxTPCNSigmas);
802806

803807
if (TMath::Abs(gamma2.mGamma()) > Pi0PhotonMaxMass ||
804808
gamma2.qtarm() >= Pi0PhotonMaxQt ||
@@ -823,14 +827,14 @@ struct sigma0builder {
823827
}
824828

825829
// Fill histograms
826-
histos.fill(HIST("Pi0QA/h2dPtVsMassPi0AfterSel_Candidates"), pi0Pt, pi0Mass);
830+
histos.fill(HIST("Pi0QA/h3dMassPi0AfterSel_Candidates"), collision.centFT0C(), pi0Pt, pi0Mass);
827831
if (fIsMC && fIsPi0)
828-
histos.fill(HIST("Pi0QA/h2dPtVsMassPi0AfterSel_MCAssoc"), pi0Pt, pi0Mass);
832+
histos.fill(HIST("Pi0QA/h3dMassPi0AfterSel_MCAssoc"), collision.centFT0C(), pi0Pt, pi0Mass);
829833
}
830834

831835
// Process photon candidate
832-
template <typename TV0Object>
833-
bool processPhotonCandidate(TV0Object const& gamma)
836+
template <typename TV0Object, typename TCollision>
837+
bool processPhotonCandidate(TV0Object const& gamma, TCollision collision)
834838
{
835839
if (gamma.v0Type() == 0)
836840
return false;
@@ -869,12 +873,13 @@ struct sigma0builder {
869873
return false;
870874
histos.fill(HIST("PhotonSel/hSelectionStatistics"), 6.);
871875
}
876+
histos.fill(HIST("PhotonSel/h3dPhotonMass"), collision.centFT0C(), gamma.pt(), gamma.mGamma());
872877
return true;
873878
}
874879

875880
// Process photon candidate
876-
template <typename TV0Object>
877-
bool processLambdaCandidate(TV0Object const& lambda)
881+
template <typename TV0Object, typename TCollision>
882+
bool processLambdaCandidate(TV0Object const& lambda, TCollision collision)
878883
{
879884
if (lambda.v0Type() != 1)
880885
return false;
@@ -912,12 +917,15 @@ struct sigma0builder {
912917
histos.fill(HIST("LambdaSel/hSelectionStatistics"), 6.);
913918
}
914919

920+
histos.fill(HIST("LambdaSel/h3dLambdaMass"), collision.centFT0C(), lambda.pt(), lambda.mLambda());
921+
histos.fill(HIST("LambdaSel/h3dALambdaMass"), collision.centFT0C(), lambda.pt(), lambda.mAntiLambda());
922+
915923
return true;
916924
}
917925
///////////
918926
// Process sigma candidate and store properties in object
919-
template <typename TV0Object>
920-
bool buildSigma0(TV0Object const& lambda, TV0Object const& gamma)
927+
template <typename TV0Object, typename TCollision>
928+
bool buildSigma0(TV0Object const& lambda, TV0Object const& gamma, TCollision collision)
921929
{
922930
// Checking if both V0s are made of the very same tracks
923931
if (gamma.posTrackExtraId() == lambda.posTrackExtraId() ||
@@ -929,31 +937,35 @@ struct sigma0builder {
929937
std::array<float, 3> pVecPhotons{gamma.px(), gamma.py(), gamma.pz()};
930938
std::array<float, 3> pVecLambda{lambda.px(), lambda.py(), lambda.pz()};
931939
auto arrMom = std::array{pVecPhotons, pVecLambda};
932-
float sigmamass = RecoDecay::m(arrMom, std::array{o2::constants::physics::MassPhoton, o2::constants::physics::MassLambda0});
933-
float sigmarap = RecoDecay::y(std::array{gamma.px() + lambda.px(), gamma.py() + lambda.py(), gamma.pz() + lambda.pz()}, o2::constants::physics::MassSigma0);
940+
float sigmaMass = RecoDecay::m(arrMom, std::array{o2::constants::physics::MassPhoton, o2::constants::physics::MassLambda0});
941+
float sigmaY = RecoDecay::y(std::array{gamma.px() + lambda.px(), gamma.py() + lambda.py(), gamma.pz() + lambda.pz()}, o2::constants::physics::MassSigma0);
942+
float SigmapT = RecoDecay::pt(array{gamma.px() + lambda.px(), gamma.py() + lambda.py()});
943+
944+
// Before any selection
945+
histos.fill(HIST("SigmaSel/h3dMassSigma0BeforeSel"), collision.centFT0C(), SigmapT, sigmaMass);
934946

935947
histos.fill(HIST("SigmaSel/hSelectionStatistics"), 1.);
936-
histos.fill(HIST("SigmaSel/hSigmaMass"), sigmamass);
937-
histos.fill(HIST("SigmaSel/hSigmaMassWindow"), sigmamass - 1.192642);
948+
histos.fill(HIST("SigmaSel/hSigmaMass"), sigmaMass);
949+
histos.fill(HIST("SigmaSel/hSigmaMassWindow"), sigmaMass - 1.192642);
938950

939951
if (fillQAhistos) {
940952
histos.fill(HIST("GeneralQA/h2dMassGammaVsK0S"), gamma.mGamma(), gamma.mK0Short());
941953
histos.fill(HIST("GeneralQA/h2dMassLambdaVsK0S"), lambda.mLambda(), lambda.mK0Short());
942954
histos.fill(HIST("GeneralQA/h2dMassGammaVsLambda"), gamma.mGamma(), gamma.mLambda());
943955
histos.fill(HIST("GeneralQA/h2dMassLambdaVsGamma"), lambda.mLambda(), lambda.mGamma());
944-
histos.fill(HIST("GeneralQA/h3dMassSigma0VsDaupTs"), gamma.pt(), lambda.pt(), sigmamass);
956+
histos.fill(HIST("GeneralQA/h3dMassSigma0VsDaupTs"), gamma.pt(), lambda.pt(), sigmaMass);
945957
}
946958

947-
if (TMath::Abs(sigmamass - 1.192642) > Sigma0Window)
959+
if (TMath::Abs(sigmaMass - 1.192642) > Sigma0Window)
948960
return false;
949961

950-
histos.fill(HIST("SigmaSel/hSigmaY"), sigmarap);
962+
histos.fill(HIST("SigmaSel/hSigmaY"), sigmaY);
951963
histos.fill(HIST("SigmaSel/hSelectionStatistics"), 2.);
952964

953-
if (TMath::Abs(sigmarap) > SigmaMaxRap)
965+
if (TMath::Abs(sigmaY) > SigmaMaxRap)
954966
return false;
955967

956-
histos.fill(HIST("SigmaSel/hSigmaMassSelected"), sigmamass);
968+
histos.fill(HIST("SigmaSel/hSigmaMassSelected"), sigmaMass);
957969
histos.fill(HIST("SigmaSel/hSelectionStatistics"), 3.);
958970

959971
if (fillQAhistos) {
@@ -963,6 +975,8 @@ struct sigma0builder {
963975
histos.fill(HIST("GeneralQA/h2dV0XY"), gamma.x(), gamma.y());
964976
}
965977

978+
histos.fill(HIST("SigmaSel/h3dMassSigma0AfterSel"), collision.centFT0C(), SigmapT, sigmaMass);
979+
966980
return true;
967981
}
968982

@@ -1127,9 +1141,7 @@ struct sigma0builder {
11271141
// Do analysis with collision-grouped V0s, retain full collision information
11281142
const uint64_t collIdx = coll.globalIndex();
11291143
auto V0s = fullV0s.sliceBy(perCollisionMCDerived, collIdx);
1130-
11311144
float centrality = coll.centFT0C();
1132-
histos.fill(HIST("hEventCentrality"), centrality);
11331145

11341146
bool fhasMCColl = false;
11351147
if (coll.has_straMCCollision())
@@ -1140,7 +1152,6 @@ struct sigma0builder {
11401152
float interactionRate = -1;
11411153
if (fGetIR) {
11421154
interactionRate = rateFetcher.fetch(ccdb.service, coll.timestamp(), coll.runNumber(), irSource, fIRCrashOnNull) * 1.e-3;
1143-
11441155
if (interactionRate < 0)
11451156
histos.get<TH1>(HIST("GeneralQA/hRunNumberNegativeIR"))->Fill(Form("%d", coll.runNumber()), 1);
11461157

@@ -1175,10 +1186,10 @@ struct sigma0builder {
11751186
histos.fill(HIST("MC/h2dPtVsCentralityBeforeSel_MCAssocALambda"), centrality, v0.pt());
11761187
}
11771188

1178-
if (processPhotonCandidate(v0)) // selecting photons
1189+
if (processPhotonCandidate(v0, coll)) // selecting photons
11791190
bestGammasArray.push_back(v0.globalIndex()); // Save indices of best gamma candidates
11801191

1181-
if (processLambdaCandidate(v0)) // selecting lambdas
1192+
if (processLambdaCandidate(v0, coll)) // selecting lambdas
11821193
bestLambdasArray.push_back(v0.globalIndex()); // Save indices of best lambda candidates
11831194
}
11841195

@@ -1189,7 +1200,7 @@ struct sigma0builder {
11891200
auto gamma1 = fullV0s.rawIteratorAt(bestGammasArray[i]);
11901201
for (size_t j = i + 1; j < bestGammasArray.size(); ++j) {
11911202
auto gamma2 = fullV0s.rawIteratorAt(bestGammasArray[j]);
1192-
runPi0QA(gamma1, gamma2);
1203+
runPi0QA(gamma1, gamma2, coll);
11931204
}
11941205
}
11951206
}
@@ -1268,16 +1279,11 @@ struct sigma0builder {
12681279
if (fIsAntiSigma)
12691280
histos.fill(HIST("MC/h2dPtVsCentralityBeforeSel_MCAssocASigma0"), centrality, SigmaMCpT);
12701281
}
1271-
1272-
histos.fill(HIST("SigmaSel/h3dMassSigma0BeforeSel"), centrality, SigmapT, SigmaMass);
1273-
1282+
12741283
// Build sigma0 candidate, please
1275-
if (!buildSigma0(lambda, gamma))
1284+
if (!buildSigma0(lambda, gamma, coll))
12761285
continue;
12771286

1278-
// Filling histos and tables
1279-
histos.fill(HIST("SigmaSel/h3dMassSigma0AfterSel"), centrality, SigmapT, SigmaMass);
1280-
12811287
if (SigmaY < 0.5) {
12821288
if (fIsSigma)
12831289
histos.fill(HIST("MC/h2dPtVsCentralityAfterSel_MCAssocSigma0"), centrality, SigmaMCpT);
@@ -1293,7 +1299,8 @@ struct sigma0builder {
12931299
PhotonCandPDGCode, PhotonCandPDGCodeMother, fIsPhotonPrimary, PhotonMCpT, fIsPhotonCorrectlyAssign,
12941300
LambdaCandPDGCode, LambdaCandPDGCodeMother, fIsLambdaPrimary, LambdaMCpT, fIsLambdaCorrectlyAssign);
12951301

1296-
fillTables(lambda, gamma, coll); // filling tables with accepted candidates
1302+
// Filling tables with accepted candidates
1303+
fillTables(lambda, gamma, coll);
12971304

12981305
nSigmaCandidates++;
12991306
if (nSigmaCandidates % 10000 == 0)
@@ -1308,21 +1315,18 @@ struct sigma0builder {
13081315
for (const auto& coll : collisions) {
13091316

13101317
if (!IsEventAccepted(coll, true))
1311-
continue;
1318+
continue;
13121319

13131320
// Do analysis with collision-grouped V0s, retain full collision information
13141321
const uint64_t collIdx = coll.globalIndex();
13151322
auto V0s = fullV0s.sliceBy(perCollisionSTDDerived, collIdx);
1316-
13171323
float centrality = coll.centFT0C();
1318-
histos.fill(HIST("hEventCentrality"), centrality);
13191324

13201325
//_______________________________________________
13211326
// Retrieving IR info
13221327
float interactionRate = -1;
13231328
if (fGetIR) {
13241329
interactionRate = rateFetcher.fetch(ccdb.service, coll.timestamp(), coll.runNumber(), irSource, fIRCrashOnNull) * 1.e-3;
1325-
13261330
if (interactionRate < 0)
13271331
histos.get<TH1>(HIST("GeneralQA/hRunNumberNegativeIR"))->Fill(Form("%d", coll.runNumber()), 1);
13281332

@@ -1336,10 +1340,10 @@ struct sigma0builder {
13361340
//_______________________________________________
13371341
// V0s loop
13381342
for (auto& v0 : V0s) {
1339-
if (processPhotonCandidate(v0)) // selecting photons
1343+
if (processPhotonCandidate(v0, coll)) // selecting photons
13401344
bestGammasArray.push_back(v0.globalIndex()); // Save indices of best gamma candidates
13411345

1342-
if (processLambdaCandidate(v0)) // selecting lambdas
1346+
if (processLambdaCandidate(v0, coll)) // selecting lambdas
13431347
bestLambdasArray.push_back(v0.globalIndex()); // Save indices of best lambda candidates
13441348
}
13451349

@@ -1350,35 +1354,25 @@ struct sigma0builder {
13501354
auto gamma1 = fullV0s.rawIteratorAt(bestGammasArray[i]);
13511355
for (size_t j = i + 1; j < bestGammasArray.size(); ++j) {
13521356
auto gamma2 = fullV0s.rawIteratorAt(bestGammasArray[j]);
1353-
runPi0QA(gamma1, gamma2);
1357+
runPi0QA(gamma1, gamma2, coll);
13541358
}
13551359
}
13561360
}
13571361

13581362
//_______________________________________________
13591363
// Sigma0 nested loop
13601364
for (size_t i = 0; i < bestGammasArray.size(); ++i) {
1361-
auto gamma = fullV0s.rawIteratorAt(bestGammasArray[i]);
1365+
auto gamma = fullV0s.rawIteratorAt(bestGammasArray[i]);
13621366

13631367
for (size_t j = 0; j < bestLambdasArray.size(); ++j) {
13641368
auto lambda = fullV0s.rawIteratorAt(bestLambdasArray[j]);
1365-
1366-
// Sigma0 candidate properties
1367-
std::array<float, 3> pVecPhotons{gamma.px(), gamma.py(), gamma.pz()};
1368-
std::array<float, 3> pVecLambda{lambda.px(), lambda.py(), lambda.pz()};
1369-
auto arrMom = std::array{pVecPhotons, pVecLambda};
1370-
float SigmaMass = RecoDecay::m(arrMom, std::array{o2::constants::physics::MassPhoton, o2::constants::physics::MassLambda0});
1371-
float SigmapT = RecoDecay::pt(array{gamma.px() + lambda.px(), gamma.py() + lambda.py()});
1372-
1373-
histos.fill(HIST("SigmaSel/h3dMassSigma0BeforeSel"), centrality, SigmapT, SigmaMass);
1374-
1369+
13751370
// Building sigma0 candidate
1376-
if (!buildSigma0(lambda, gamma))
1377-
continue;
1378-
1379-
fillTables(lambda, gamma, coll); // filling tables with accepted candidates
1371+
if (!buildSigma0(lambda, gamma, coll))
1372+
continue;
13801373

1381-
histos.fill(HIST("SigmaSel/h3dMassSigma0AfterSel"), centrality, SigmapT, SigmaMass);
1374+
// Filling tables with accepted candidates
1375+
fillTables(lambda, gamma, coll);
13821376

13831377
nSigmaCandidates++;
13841378
if (nSigmaCandidates % 10000 == 0)

0 commit comments

Comments
 (0)