Skip to content
Closed
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
24 changes: 12 additions & 12 deletions PWGHF/HFC/TableProducer/correlatorD0Hadrons.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -137,15 +137,15 @@ struct HfCorrelatorD0HadronsSelection {
}
float cent = 0.;
if (useCentrality) {
cent = collision.centFT0M();
cent = collision.centFT0M();
}
if (useSel8) {
isSel8 = collision.sel8();
}
if (selNoSameBunchPileUpColl) {
isNosameBunchPileUp = collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup);
}
isCentInRange = (cent >= centMin && cent <=centMax);
isCentInRange = (cent >= centMin && cent <= centMax);
isSelColl = isD0Found && isSel8 && isNosameBunchPileUp && isCentInRange;
collisionsWithSelD0(isSelColl);
}
Expand Down Expand Up @@ -295,7 +295,7 @@ struct HfCorrelatorD0Hadrons {
AxisSpec axisBdtScoreBkg = {100, 0., 1., "Bdt score background"};
AxisSpec axisBdtScorePrompt = {100, 0., 1., "Bdt score prompt"};
AxisSpec axisOrigin = {10, 0., 10., "Candidate origin"};
AxisSpec axisCent = {binsCentFt0m,"Centrality"};
AxisSpec axisCent = {binsCentFt0m, "Centrality"};

// Histograms for Data
registry.add("hPtCand", "D0, D0bar candidates", {HistType::kTH1F, {axisPtD}});
Expand All @@ -305,9 +305,9 @@ struct HfCorrelatorD0Hadrons {
registry.add("hEta", "D0,D0bar candidates", {HistType::kTH1F, {axisEta}});
registry.add("hPhi", "D0,D0bar candidates", {HistType::kTH1F, {axisPhi}});
registry.add("hY", "D0,D0bar candidates", {HistType::kTH1F, {axisRapidity}});
registry.add("hCentFT0M","Centrality FT0M;centrality;entries", {HistType::kTH1F,{{100, 0., 100.}}});
registry.add("hCentFT0M", "Centrality FT0M;centrality;entries", {HistType::kTH1F, {{100, 0., 100.}}});
registry.add("hMultiplicityPreSelection", "multiplicity prior to selection;multiplicity;entries", {HistType::kTH1F, {axisMultiplicity}});
registry.add("hMultiplicity", "multiplicity;multiplicity;entries", {HistType::kTH1F, {axisMultiplicity}});
registry.add("hMultiplicity", "multiplicity;multiplicity;entries", {HistType::kTH1F, {axisMultiplicity}});
registry.add("hMass", "D0, D0bar candidates massVsPt", {HistType::kTH2F, {{axisMassD}, {axisPtD}}});
registry.add("hMass1D", "D0, D0bar candidates mass", {HistType::kTH1F, {axisMassD}});
registry.add("hMassD01D", "D0 candidates mass", {HistType::kTH1F, {axisMassD}});
Expand Down Expand Up @@ -370,7 +370,7 @@ struct HfCorrelatorD0Hadrons {
leadingIndex = findLeadingParticle(tracks, etaTrackMax.value);
}
float cent = 0.;
if (useCentrality){
if (useCentrality) {
cent = collision.centFT0M();
}

Expand Down Expand Up @@ -511,7 +511,7 @@ struct HfCorrelatorD0Hadrons {
entryD0HadronGenInfo(false, false, 0);
entryD0HadronMlInfo(outputMlD0[0], outputMlD0[1], outputMlD0[2], outputMlD0bar[0], outputMlD0bar[1], outputMlD0bar[2]);
entryTrackRecoInfo(track.dcaXY(), track.dcaZ(), track.tpcNClsCrossedRows());
registry.fill(HIST("hCentFT0M"),cent );
registry.fill(HIST("hCentFT0M"), cent);

} // end inner loop (tracks)

Expand Down Expand Up @@ -637,7 +637,7 @@ struct HfCorrelatorD0Hadrons {

flagD0 = candidate.flagMcMatchRec() == o2::hf_decay::hf_cand_2prong::DecayChannelMain::D0ToPiK; // flagD0Signal 'true' if candidate matched to D0 (particle)
flagD0bar = candidate.flagMcMatchRec() == -o2::hf_decay::hf_cand_2prong::DecayChannelMain::D0ToPiK; // flagD0Reflection 'true' if candidate, selected as D0 (particle), is matched to D0bar (antiparticle)
float cent = 100.0; //Centrality Placeholder: will be updated later
float cent = 100.0; // Centrality Placeholder: will be updated later
// ========== track loop starts here ========================

for (const auto& track : tracks) {
Expand Down Expand Up @@ -753,7 +753,7 @@ struct HfCorrelatorD0Hadrons {
bool isD0Prompt = false;
bool isD0NonPrompt = false;
int trackOrigin = -1;
float cent = 100.; //Centrality Placeholder: will be updated later
float cent = 100.; // Centrality Placeholder: will be updated later

for (const auto& particleTrigg : mcParticles) {
if (std::abs(particleTrigg.pdgCode()) != Pdg::kD0) {
Expand Down Expand Up @@ -879,7 +879,7 @@ struct HfCorrelatorD0Hadrons {
invMassDstar2 = std::sqrt((eKPi + ePion) * (eKPi + ePion) - pSum2);
std::vector<float> outputMlD0 = {-1., -1., -1.};
std::vector<float> outputMlD0bar = {-1., -1., -1.};
float cent = 100.; //Centrality Placeholder: will be updated later
float cent = 100.; // Centrality Placeholder: will be updated later

if (candidate.isSelD0() >= selectionFlagD0) {
if ((std::abs(invMassDstar1 - HfHelper::invMassD0ToPiK(candidate)) - softPiMass) < ptSoftPionMax) {
Expand Down Expand Up @@ -945,7 +945,7 @@ struct HfCorrelatorD0Hadrons {
registry.fill(HIST("hD0PoolBin"), poolBinD0);
registry.fill(HIST("hMultFT0M"), c1.multFT0M());
registry.fill(HIST("hZvtx"), c1.posZ());
float cent = 100.; //Centrality Placeholder: will be updated later
float cent = 100.; // Centrality Placeholder: will be updated later

for (const auto& [candidate, particleAssoc] : o2::soa::combinations(o2::soa::CombinationsFullIndexPolicy(tracks1, tracks2))) {

Expand Down Expand Up @@ -1092,7 +1092,7 @@ struct HfCorrelatorD0Hadrons {
if (std::abs(particleAssoc.pdgCode()) == kPiPlus && indexMotherPi >= 0 && indexMotherD0 >= 0 && indexMotherPi == indexMotherD0) {
continue;
}
float cent = 100.; //Centrality Placeholder: will be updated later
float cent = 100.; // Centrality Placeholder: will be updated later
bool correlationStatus = false;
int trackOrigin = RecoDecay::getCharmHadronOrigin(mcParticles, particleAssoc, true);
bool isD0Prompt = particleTrigg.originMcGen() == RecoDecay::OriginType::Prompt;
Expand Down
4 changes: 2 additions & 2 deletions PWGHF/HFC/Tasks/taskCorrelationD0Hadrons.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -172,7 +172,7 @@ struct HfTaskCorrelationD0Hadrons {
registry.add("hDeltaEtaPtIntSidebands", "D0-h deltaEta sidebands", {HistType::kTH1F, {axisDeltaEta}});
registry.add("hDeltaPhiPtIntSidebands", "D0-h deltaPhi sidebands", {HistType::kTH1F, {axisDeltaPhi}});
registry.add("hCorrel2DPtIntSidebands", "D0-h deltaPhi vs deltaEta sidebands", {HistType::kTH2F, {{axisDeltaPhi}, {axisDeltaEta}}});
registry.add("hCorrel2DVsPtSidebands", "D0-h correlations sidebands", {HistType::kTHnSparseD, {{axisDeltaPhi}, {axisDeltaEta}, {axisPtD}, {axisPtHadron}, {axisPoolBin},{axisCentFT0M}}});
registry.add("hCorrel2DVsPtSidebands", "D0-h correlations sidebands", {HistType::kTHnSparseD, {{axisDeltaPhi}, {axisDeltaEta}, {axisPtD}, {axisPtHadron}, {axisPoolBin}, {axisCentFT0M}}});
registry.add("hDeltaEtaPtIntSidebandsSoftPi", "D0-h deltaEta sidebands soft pi only", {HistType::kTH1F, {axisDeltaEta}});
registry.add("hDeltaPhiPtIntSidebandsSoftPi", "D0-h deltaPhi sidebands soft pi only", {HistType::kTH1F, {axisDeltaPhi}});
registry.add("hCorrel2DPtIntSidebandsSoftPi", "D0-h deltaPhi vs deltaEta sidebands soft pi only", {HistType::kTH2F, {{axisDeltaPhi}, {axisDeltaEta}}});
Expand Down Expand Up @@ -327,7 +327,7 @@ struct HfTaskCorrelationD0Hadrons {
float cent = 0.;
if (useCentrality) {
cent = pairEntry.cent();
}
}
// define variables for widely used quantities
double const deltaPhi = pairEntry.deltaPhi();
double const deltaEta = pairEntry.deltaEta();
Expand Down