Skip to content

Commit ad70bf8

Browse files
authored
[PWGCF] fill trigger separately, fix FT0A-FT0C container bug (#13237)
1 parent 5e6c4cf commit ad70bf8

File tree

1 file changed

+9
-4
lines changed

1 file changed

+9
-4
lines changed

PWGCF/TwoParticleCorrelations/Tasks/longRangeDihadronCor.cxx

Lines changed: 9 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -279,7 +279,8 @@ struct LongRangeDihadronCor {
279279
registry.add("Centrality_used", hCentTitle.c_str(), {HistType::kTH1D, {{100, 0, 100}}}); // histogram to see how many events are in the same and mixed event
280280
registry.add("zVtx", "zVtx", {HistType::kTH1D, {axisVertex}});
281281
registry.add("zVtx_used", "zVtx_used", {HistType::kTH1D, {axisVertex}});
282-
registry.add("Trig_hist", "", {HistType::kTHnSparseF, {{axisSample, axisVertex, axisPtTrigger}}});
282+
registry.add("Trig_hist_TPC_FT0A", "", {HistType::kTHnSparseF, {{axisSample, axisVertex, axisPtTrigger}}});
283+
registry.add("Trig_hist_TPC_FT0C", "", {HistType::kTHnSparseF, {{axisSample, axisVertex, axisPtTrigger}}});
283284
registry.add("FT0Amp", "", {HistType::kTH2F, {axisChID, axisFit}});
284285
registry.add("FT0AmpCorrect", "", {HistType::kTH2F, {axisChID, axisFit}});
285286
registry.add("FT0Cmp", "", {HistType::kTH2F, {axisChID, axisFit}});
@@ -558,7 +559,11 @@ struct LongRangeDihadronCor {
558559
if (!getEfficiencyCorrection(triggerWeight, track1.eta(), track1.pt(), posZ))
559560
continue;
560561
if (system == SameEvent) {
561-
registry.fill(HIST("Trig_hist"), fSampleIndex, posZ, track1.pt(), eventWeight * triggerWeight);
562+
if (corType == kFT0C) {
563+
registry.fill(HIST("Trig_hist_TPC_FT0C"), fSampleIndex, posZ, track1.pt(), eventWeight * triggerWeight);
564+
} else if (corType == kFT0A) {
565+
registry.fill(HIST("Trig_hist_TPC_FT0A"), fSampleIndex, posZ, track1.pt(), eventWeight * triggerWeight);
566+
}
562567
}
563568

564569
std::size_t channelSize = 0;
@@ -630,10 +635,10 @@ struct LongRangeDihadronCor {
630635
float deltaEta = etaA - etaC;
631636
// fill the right sparse and histograms
632637
if (system == SameEvent) {
633-
sameFt0aFt0c->getPairHist()->Fill(step, fSampleIndex, posZ, 0., 0., deltaPhi, deltaEta, amplA * amplC * eventWeight * triggerWeight);
638+
sameFt0aFt0c->getPairHist()->Fill(step, fSampleIndex, posZ, 0.5, 0.5, deltaPhi, deltaEta, amplA * amplC * eventWeight * triggerWeight);
634639
registry.fill(HIST("deltaEta_deltaPhi_same_FT0A_FT0C"), deltaPhi, deltaEta, amplA * amplC * eventWeight * triggerWeight);
635640
} else if (system == MixedEvent) {
636-
mixedFt0aFt0c->getPairHist()->Fill(step, fSampleIndex, posZ, 0., 0., deltaPhi, deltaEta, amplA * amplC * eventWeight * triggerWeight);
641+
mixedFt0aFt0c->getPairHist()->Fill(step, fSampleIndex, posZ, 0.5, 0.5, deltaPhi, deltaEta, amplA * amplC * eventWeight * triggerWeight);
637642
registry.fill(HIST("deltaEta_deltaPhi_mixed_FT0A_FT0C"), deltaPhi, deltaEta, amplA * amplC * eventWeight * triggerWeight);
638643
}
639644
}

0 commit comments

Comments
 (0)