Skip to content

[PWGHF/HFL] Z task: Add charge in histograms #12068

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
vkucera merged 15 commits into from
Jul 17, 2025
Merged

[PWGHF/HFL] Z task: Add charge in histograms #12068

Changes from all commits
Commits
Show all changes
15 commits
Select commit Hold shift + click to select a range
4fd72ee
implememted Z-h correlation
sashingo Jul 7, 2025
32c1341
fixed O2 linter error
sashingo Jul 7, 2025
187dfc9
fixed O2 linter errors
sashingo Jul 7, 2025
1101165
implemetned comments by Vit
sashingo Jul 7, 2025
e6d27cc
revert the changes
sashingo Jul 7, 2025
ae09c84
minor corrections
sashingo Jul 8, 2025
3033f8d
modified some name of parameters
sashingo Jul 8, 2025
a7ceea5
fixed O2 linter errors
sashingo Jul 8, 2025
8483257
implemented comment
sashingo Jul 9, 2025
2dddc15
adding missing braces
sashingo Jul 10, 2025
da8c3f8
Merge branch 'AliceO2Group:master' into master
sashingo Jul 14, 2025
c934352
Merge branch 'AliceO2Group:master' into master
sashingo Jul 15, 2025
8a9aab4
included charge info in histogrmas for Z measurement
sashingo Jul 15, 2025
c3b7601
changed name of axis for Z pt histogram
sashingo Jul 15, 2025
144dfe4
changed TH1F to TH1D for histograms of event counters
sashingo Jul 16, 2025
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
25 changes: 15 additions & 10 deletions PWGHF/HFL/Tasks/taskElectronWeakBoson.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -213,12 +213,13 @@ struct HfTaskElectronWeakBoson {
const AxisSpec axisTrigger{3, -0.5, 2.5, "Trigger status of zorro"};
const AxisSpec axisDPhiZh{64, -o2::constants::math::PIHalf, 3 * o2::constants::math::PIHalf, "#Delta#phi(Z-h)"};
const AxisSpec axisPtHadron{50, 0, 50, "p_{T,hadron} (GeV/c)"};
const AxisSpec axisZpt{150, 0, 150, "p_{T,Z} (GeV/c)"};
const AxisSpec axisPtZ{150, 0, 150, "p_{T,Z} (GeV/c)"};
const AxisSpec axisSign{2, -2, 2, "charge sign"};

// create registrygrams
registry.add("hZvtx", "Z vertex", kTH1F, {axisZvtx});
registry.add("hEventCounterInit", "hEventCounterInit", kTH1F, {axisCounter});
registry.add("hEventCounter", "hEventCounter", kTH1F, {axisCounter});
registry.add("hZvtx", "Z vertex", kTH1D, {axisZvtx});
registry.add("hEventCounterInit", "hEventCounterInit", kTH1D, {axisCounter});
registry.add("hEventCounter", "hEventCounter", kTH1D, {axisCounter});
registry.add("hITSchi2", "ITS #chi^{2}", kTH1F, {axisChi2});
registry.add("hTPCchi2", "TPC #chi^{2}", kTH1F, {axisChi2});
registry.add("hTPCnCls", "TPC NCls", kTH1F, {axisCluster});
Expand Down Expand Up @@ -248,11 +249,11 @@ struct HfTaskElectronWeakBoson {
registry.add("hTHnTrMatch", "Track EMC Match", HistType::kTHnSparseF, {axisPt, axisdPhi, axisdEta});

// Z-hadron correlation histograms
registry.add("hZHadronDphi", "Z-hadron #Delta#phi correlation", kTH2F, {{axisZpt}, {axisDPhiZh}});
registry.add("hZptSpectrum", "Z boson p_{T} spectrum", kTH1F, {axisZpt});
registry.add("hZHadronDphi", "Z-hadron #Delta#phi correlation", HistType::kTHnSparseF, {axisSign, axisPtZ, axisDPhiZh});
registry.add("hZptSpectrum", "Z boson p_{T} spectrum", kTH2F, {{axisSign}, {axisPtZ}});

// hisotgram for EMCal trigger
registry.add("hEMCalTrigger", "EMCal trigger", kTH1F, {axisTrigger});
registry.add("hEMCalTrigger", "EMCal trigger", kTH1D, {axisTrigger});
}

double getIsolatedCluster(const o2::aod::EMCALCluster& cluster,
Expand Down Expand Up @@ -351,6 +352,9 @@ struct HfTaskElectronWeakBoson {
if (zeeKF.GetNDF() < 1) {
continue;
}
if (zeeKF.GetChi2() < 0) {
continue;
}
if (chiSqNdf > chiSqNdfMax) {
continue;
}
Expand Down Expand Up @@ -605,9 +609,10 @@ struct HfTaskElectronWeakBoson {
if (reconstructedZ.size() > 0) {
for (const auto& zBoson : reconstructedZ) {
// Z boson selection
if (zBoson.mass < massZMin || zBoson.mass > massZMax)
if (zBoson.mass < massZMin || zBoson.mass > massZMax) {
continue;
registry.fill(HIST("hZptSpectrum"), zBoson.pt);
}
registry.fill(HIST("hZptSpectrum"), zBoson.charge, zBoson.pt);
for (const auto& trackAss : selectedElectronsAss) {
if (std::abs(trackAss.pt - zBoson.ptchild0) < ptMatch) {
continue;
Expand All @@ -617,7 +622,7 @@ struct HfTaskElectronWeakBoson {
}
// calculate Z-h correlation
double deltaPhi = RecoDecay::constrainAngle(trackAss.phi - zBoson.phi, -o2::constants::math::PIHalf);
registry.fill(HIST("hZHadronDphi"), zBoson.pt, deltaPhi);
registry.fill(HIST("hZHadronDphi"), zBoson.charge, zBoson.pt, deltaPhi);
}
}
} // end of Z-hadron correlation
Expand Down
Loading