Skip to content

Commit 7d3147f

Browse files
zeptoterazeptotera
committed
changes of Configurables Adding existing table LCHadronPair
1 parent 1b079f1 commit 7d3147f

File tree

4 files changed

+96
-41
lines changed

4 files changed

+96
-41
lines changed

‎PWGHF/HFC/DataModel/CorrelationTables.h‎

Lines changed: 0 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -172,14 +172,6 @@ DECLARE_SOA_COLUMN(PiNsigmTOF, piNsigmTOF, float); //! A
172172
} // namespace hf_correlation_lc_hadron
173173

174174
DECLARE_SOA_TABLE(LcHadronPair, "AOD", "LCHPAIR", //! Lc-Hadrons pairs Informations
175-
aod::hf_correlation_lc_hadron::DeltaPhi,
176-
aod::hf_correlation_lc_hadron::DeltaEta,
177-
aod::hf_correlation_lc_hadron::PtLc,
178-
aod::hf_correlation_lc_hadron::PtHadron,
179-
aod::hf_correlation_lc_hadron::PoolBin,
180-
aod::hf_correlation_lc_hadron::IsAutoCorrelated);
181-
182-
DECLARE_SOA_TABLE(LcHadronPairC, "AOD", "LCHPAIRC", //! Lc-Hadrons pairs Informations
183175
aod::hf_correlation_lc_hadron::DeltaPhi,
184176
aod::hf_correlation_lc_hadron::DeltaEta,
185177
aod::hf_correlation_lc_hadron::PtLc,

‎PWGHF/HFC/TableProducer/correlatorLcHadrons.cxx‎

Lines changed: 24 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -201,7 +201,6 @@ struct HfCorrelatorLcHadronsSelection {
201201
// Lc-Hadron correlation pair builder - for real data and data-like analysis (i.e. reco-level w/o matching request via Mc truth)
202202
struct HfCorrelatorLcHadrons {
203203
Produces<aod::LcHadronPair> entryLcHadronPair;
204-
Produces<aod::LcHadronPairC> entryLcHadronPairC;
205204
Produces<aod::LcHadronPairY> entryLcHadronPairY;
206205
Produces<aod::LcHadronPairTrkPID> entryLcHadronPairTrkPID;
207206
Produces<aod::LcHadronRecoInfo> entryLcHadronRecoInfo;
@@ -409,7 +408,6 @@ struct HfCorrelatorLcHadrons {
409408

410409
float cent=evaluateCentralityColl(collision);
411410

412-
413411
int poolBin = corrBinning.getBin(std::make_tuple(collision.posZ(), collision.multFT0M()));
414412
int nTracks = 0;
415413
if (collision.numContrib() > 1) {
@@ -495,7 +493,7 @@ struct HfCorrelatorLcHadrons {
495493
}
496494
}
497495
if (candidate.isSelLcToPKPi() >= selectionFlagLc) {
498-
entryLcHadronPairC(getDeltaPhi(track.phi(), candidate.phi()),
496+
entryLcHadronPair(getDeltaPhi(track.phi(), candidate.phi()),
499497
track.eta() - candidate.eta(),
500498
candidate.pt() * chargeLc,
501499
track.pt() * track.sign(),
@@ -512,7 +510,7 @@ struct HfCorrelatorLcHadrons {
512510
}
513511
}
514512
if (candidate.isSelLcToPiKP() >= selectionFlagLc) {
515-
entryLcHadronPairC(getDeltaPhi(track.phi(), candidate.phi()),
513+
entryLcHadronPair(getDeltaPhi(track.phi(), candidate.phi()),
516514
track.eta() - candidate.eta(),
517515
candidate.pt() * chargeLc,
518516
track.pt() * track.sign(),
@@ -707,6 +705,7 @@ struct HfCorrelatorLcHadrons {
707705
}
708706
//}
709707
}
708+
double_t cent =100.0; //will be updated later
710709

711710
// Lc-Hadron correlation dedicated section
712711
// if the candidate is selected as Lc, search for Hadron ad evaluate correlations
@@ -756,7 +755,8 @@ struct HfCorrelatorLcHadrons {
756755
candidate.pt() * chargeLc,
757756
track.pt() * track.sign(),
758757
poolBin,
759-
correlationStatus);
758+
correlationStatus,
759+
cent);
760760
entryLcHadronPairY(track.rapidity(MassProton) - hfHelper.yLc(candidate)); // only for proton as of now
761761
entryLcHadronRecoInfo(hfHelper.invMassLcToPKPi(candidate), isLcSignal);
762762
if (fillTrkPID) {
@@ -786,7 +786,8 @@ struct HfCorrelatorLcHadrons {
786786
candidate.pt() * chargeLc,
787787
track.pt() * track.sign(),
788788
poolBin,
789-
correlationStatus);
789+
correlationStatus,
790+
cent);
790791
entryLcHadronPairY(track.rapidity(MassProton) - hfHelper.yLc(candidate)); // only for proton as of now
791792
entryLcHadronRecoInfo(hfHelper.invMassLcToPiKP(candidate), isLcSignal);
792793
if (fillTrkPID) {
@@ -918,7 +919,8 @@ struct HfCorrelatorLcHadrons {
918919
particle.pt() * chargeLc / std::abs(chargeLc),
919920
particleAssoc.pt() * chargeAssoc / std::abs(chargeAssoc),
920921
poolBin,
921-
correlationStatus);
922+
correlationStatus,
923+
mcCollision.multMCFT0A());
922924
entryLcHadronPairY(particleAssoc.y() - yL);
923925
entryLcHadronRecoInfo(MassLambdaCPlus, true);
924926
entryLcHadronGenInfo(isLcPrompt, particleAssoc.isPhysicalPrimary(), trackOrigin);
@@ -950,6 +952,8 @@ struct HfCorrelatorLcHadrons {
950952

951953
auto trackPos1 = trigLc.template prong0_as<TracksData>(); // positive daughter (negative for the antiparticles)
952954
int8_t chargeLc = trackPos1.sign(); // charge of 1st prong will be the charge of Lc candidate
955+
956+
double_t cent =100.0; //will be updated later
953957

954958
std::vector<float> outputMl = {-1., -1., -1.};
955959
// LcToPKPi and LcToPiKP division
@@ -959,7 +963,8 @@ struct HfCorrelatorLcHadrons {
959963
trigLc.pt() * chargeLc,
960964
assocParticle.pt() * assocParticle.sign(),
961965
poolBin,
962-
correlationStatus);
966+
correlationStatus,
967+
cent);
963968
entryLcHadronPairY(assocParticle.y() - hfHelper.yLc(trigLc));
964969
entryLcHadronRecoInfo(hfHelper.invMassLcToPKPi(trigLc), false);
965970
entryLcHadronGenInfo(false, false, 0);
@@ -978,7 +983,8 @@ struct HfCorrelatorLcHadrons {
978983
trigLc.pt() * chargeLc,
979984
assocParticle.pt() * assocParticle.sign(),
980985
poolBin,
981-
correlationStatus);
986+
correlationStatus,
987+
cent);
982988
entryLcHadronPairY(assocParticle.y() - hfHelper.yLc(trigLc));
983989
entryLcHadronRecoInfo(hfHelper.invMassLcToPiKP(trigLc), false);
984990
entryLcHadronGenInfo(false, false, 0);
@@ -1038,14 +1044,15 @@ struct HfCorrelatorLcHadrons {
10381044
}
10391045
auto tracksTuple = std::make_tuple(candidates, tracks);
10401046
Pair<SelCollisionsWithLc, CandidatesLcMcRec, TracksWithMc, BinningType> pairMcRec{corrBinning, numberEventsMixed, -1, collisions, tracksTuple, &cache};
1041-
1047+
10421048
for (const auto& [c1, tracks1, c2, tracks2] : pairMcRec) {
10431049
int poolBin = corrBinning.getBin(std::make_tuple(c2.posZ(), c2.multFT0M()));
10441050
int poolBinLc = corrBinning.getBin(std::make_tuple(c1.posZ(), c1.multFT0M()));
10451051
registry.fill(HIST("hMultFT0M"), c1.multFT0M());
10461052
registry.fill(HIST("hZvtx"), c1.posZ());
10471053
registry.fill(HIST("hTracksPoolBin"), poolBin);
10481054
registry.fill(HIST("hLcPoolBin"), poolBinLc);
1055+
double_t cent =100.0; //will be updated later
10491056
for (const auto& [candidate, pAssoc] : o2::soa::combinations(o2::soa::CombinationsFullIndexPolicy(tracks1, tracks2))) {
10501057
if (std::abs(hfHelper.yLc(candidate)) > yCandMax || candidate.pt() < ptCandMin || candidate.pt() > ptCandMax) {
10511058
continue;
@@ -1078,7 +1085,8 @@ struct HfCorrelatorLcHadrons {
10781085
candidate.pt() * chargeLc,
10791086
pAssoc.pt() * pAssoc.sign(),
10801087
poolBin,
1081-
correlationStatus);
1088+
correlationStatus,
1089+
cent);
10821090
entryLcHadronPairY(pAssoc.y() - hfHelper.yLc(candidate));
10831091
entryLcHadronRecoInfo(hfHelper.invMassLcToPKPi(candidate), isLcSignal);
10841092
entryLcHadronGenInfo(isLcPrompt, isPhysicalPrimary, trackOrigin);
@@ -1097,7 +1105,8 @@ struct HfCorrelatorLcHadrons {
10971105
candidate.pt() * chargeLc,
10981106
pAssoc.pt() * pAssoc.sign(),
10991107
poolBin,
1100-
correlationStatus);
1108+
correlationStatus,
1109+
cent);
11011110
entryLcHadronPairY(pAssoc.y() - hfHelper.yLc(candidate));
11021111
entryLcHadronRecoInfo(hfHelper.invMassLcToPiKP(candidate), isLcSignal);
11031112
entryLcHadronGenInfo(isLcPrompt, isPhysicalPrimary, trackOrigin);
@@ -1145,6 +1154,7 @@ struct HfCorrelatorLcHadrons {
11451154
}
11461155
int8_t chargeLc = pdg->GetParticle(candidate.pdgCode())->Charge(); // Retrieve charge
11471156
int8_t chargeAssoc = pdg->GetParticle(particleAssoc.pdgCode())->Charge(); // Retrieve charge
1157+
double_t cent =100.0; //will be updated later
11481158

11491159
int trackOrigin = RecoDecay::getCharmHadronOrigin(mcParticles, particleAssoc, true);
11501160
bool isLcPrompt = candidate.originMcGen() == RecoDecay::OriginType::Prompt;
@@ -1153,7 +1163,8 @@ struct HfCorrelatorLcHadrons {
11531163
candidate.pt() * chargeLc / std::abs(chargeLc),
11541164
particleAssoc.pt() * chargeAssoc / std::abs(chargeAssoc),
11551165
poolBin,
1156-
correlationStatus);
1166+
correlationStatus,
1167+
cent);
11571168
entryLcHadronPairY(particleAssoc.y() - yL);
11581169
entryLcHadronRecoInfo(MassLambdaCPlus, true);
11591170
entryLcHadronGenInfo(isLcPrompt, particleAssoc.isPhysicalPrimary(), trackOrigin);

‎PWGHF/HFC/TableProducer/correlatorLcScHadrons.cxx‎

Lines changed: 12 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -476,13 +476,15 @@ struct HfCorrelatorLcScHadrons {
476476
{
477477
bool isPhysicalPrimary = false;
478478
int trackOrigin = -1;
479+
double_t cent =100.0; //will be updated later
479480

480481
entryCandHadronPair(getDeltaPhi(track.phi(), candidate.phi()),
481482
track.eta() - candidate.eta(),
482483
candidate.pt(),
483484
track.pt() * track.sign(),
484485
binPool,
485-
correlStatus);
486+
correlStatus,
487+
cent);
486488
entryCandHadronPairY(track.rapidity(MassProton) - yCand);
487489
entryCandHadronMlInfo(outMl[0], outMl[1]);
488490
entryTrackRecoInfo(track.dcaXY(), track.dcaZ(), track.tpcNClsCrossedRows());
@@ -1016,6 +1018,10 @@ struct HfCorrelatorLcScHadrons {
10161018
}
10171019
}
10181020

1021+
int8_t chargeLc = pdg->GetParticle(particle.pdgCode())->Charge(); // Retrieve charge
1022+
int8_t chargeAssoc = pdg->GetParticle(particleAssoc.pdgCode())->Charge(); // Retrieve charge
1023+
double_t cent =100.0; //will be updated later
1024+
10191025
int trackOrigin = RecoDecay::getCharmHadronOrigin(mcParticles, particleAssoc, true);
10201026
int8_t chargeAssoc = pdg->GetParticle(particleAssoc.pdgCode())->Charge(); // Retrieve charge
10211027
chargeAssoc = chargeAssoc / std::abs(chargeAssoc);
@@ -1025,7 +1031,8 @@ struct HfCorrelatorLcScHadrons {
10251031
particle.pt(),
10261032
particleAssoc.pt() * chargeAssoc,
10271033
poolBin,
1028-
correlationStatus);
1034+
correlationStatus,
1035+
cent);
10291036
entryCandHadronPairY(particleAssoc.y() - yCand);
10301037
entryCandHadronRecoInfo(massCand, true);
10311038
entryCandHadronGenInfo(isPrompt, particleAssoc.isPhysicalPrimary(), trackOrigin);
@@ -1166,6 +1173,7 @@ struct HfCorrelatorLcScHadrons {
11661173
}
11671174
int8_t chargeLc = pdg->GetParticle(candidate.pdgCode())->Charge(); // Retrieve charge
11681175
int8_t chargeAssoc = pdg->GetParticle(particleAssoc.pdgCode())->Charge(); // Retrieve charge
1176+
double_t cent =100.0; //will be updated later
11691177

11701178
int trackOrigin = RecoDecay::getCharmHadronOrigin(mcParticles, particleAssoc, true);
11711179
bool isPrompt = candidate.originMcGen() == RecoDecay::OriginType::Prompt;
@@ -1174,7 +1182,8 @@ struct HfCorrelatorLcScHadrons {
11741182
candidate.pt() * chargeLc / std::abs(chargeLc),
11751183
particleAssoc.pt() * chargeAssoc / std::abs(chargeAssoc),
11761184
poolBin,
1177-
correlationStatus);
1185+
correlationStatus,
1186+
cent);
11781187
entryCandHadronPairY(particleAssoc.y() - yL);
11791188
entryCandHadronRecoInfo(MassLambdaCPlus, true);
11801189
entryCandHadronGenInfo(isPrompt, particleAssoc.isPhysicalPrimary(), trackOrigin);

0 commit comments

Comments
 (0)