Skip to content

[PWGHF] Z: Add eta and dE/dx axes. Add electron momentum correction factor. #12303

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 21 commits into from
Aug 1, 2025
Merged

[PWGHF] Z: Add eta and dE/dx axes. Add electron momentum correction factor. #12303

Changes from all commits
Commits
Show all changes
21 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
7ca039a
Merge branch 'AliceO2Group:master' into master
sashingo Jul 20, 2025
953b290
updated to check eID
sashingo Jul 29, 2025
964a2ed
Merge branch 'AliceO2Group:master' into master
sashingo Jul 29, 2025
37da52b
changed ptZeeCorrfac to ptEleCorrfac
sashingo Jul 29, 2025
5136521
Update taskElectronWeakBoson.cxx
sashingo Jul 30, 2025
f80040b
Update taskElectronWeakBoson.cxx
sashingo Jul 30, 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
10 changes: 6 additions & 4 deletions PWGHF/HFL/Tasks/taskElectronWeakBoson.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -103,6 +103,7 @@ struct HfTaskElectronWeakBoson {

Configurable<float> massZMin{"massZMin", 60.0, "Minimum Z mass (GeV/c^2)"};
Configurable<float> massZMax{"massZMax", 120.0, "Maximum Z mass (GeV/c^2)"};
Configurable<float> correctionPtElectron{"correctionPtElectron", 1.0, "momentum correction factor for decay electrons from Z boson"};

// flag for THn
Configurable<bool> isTHnElectron{"isTHnElectron", true, "Enables THn for electrons"};
Expand Down Expand Up @@ -195,6 +196,7 @@ struct HfTaskElectronWeakBoson {
const AxisSpec axisEta{20, -1.0, 1.0, "#eta"};
const AxisSpec axisPt{nBinsPt, 0, binPtmax, "p_{T}"};
const AxisSpec axisNsigma{100, -5, 5, "N#sigma"};
const AxisSpec axisDedx{150, 0, 150, "dEdx"};
const AxisSpec axisE{nBinsE, 0, binEmax, "Energy"};
const AxisSpec axisM02{100, 0, 1, "M02"};
const AxisSpec axisdPhi{100, -0.5, 0.5, "dPhi"};
Expand Down Expand Up @@ -245,7 +247,7 @@ struct HfTaskElectronWeakBoson {
registry.add("hInvMassZeeUls", "invariant mass for Z ULS pair", kTH2F, {{axisPt}, {axisInvMassZ}});
registry.add("hKfInvMassZeeLs", "invariant mass for Z LS pair KFp", kTH2F, {{axisPt}, {axisInvMassZ}});
registry.add("hKfInvMassZeeUls", "invariant mass for Z ULS pair KFp", kTH2F, {{axisPt}, {axisInvMassZ}});
registry.add("hTHnElectrons", "electron info", HistType::kTHnSparseF, {axisPt, axisNsigma, axisM02, axisEop, axisIsoEnergy, axisIsoTrack});
registry.add("hTHnElectrons", "electron info", HistType::kTHnSparseF, {axisPt, axisNsigma, axisM02, axisEop, axisIsoEnergy, axisIsoTrack, axisEta, axisDedx});
registry.add("hTHnTrMatch", "Track EMC Match", HistType::kTHnSparseF, {axisPt, axisdPhi, axisdEta});

// Z-hadron correlation histograms
Expand Down Expand Up @@ -338,8 +340,8 @@ struct HfTaskElectronWeakBoson {
KFPTrack kfpTrackAssEle = createKFPTrackFromTrack(track);
KFParticle kfpAssEle(kfpTrackAssEle, pdgAss);
// reco by RecoDecay
auto child1 = RecoDecayPtEtaPhi::pVector(kfpIsoEle.GetPt(), kfpIsoEle.GetEta(), kfpIsoEle.GetPhi());
auto child2 = RecoDecayPtEtaPhi::pVector(kfpAssEle.GetPt(), kfpAssEle.GetEta(), kfpAssEle.GetPhi());
auto child1 = RecoDecayPtEtaPhi::pVector(kfpIsoEle.GetPt() * correctionPtElectron, kfpIsoEle.GetEta(), kfpIsoEle.GetPhi());
auto child2 = RecoDecayPtEtaPhi::pVector(kfpAssEle.GetPt() * correctionPtElectron, kfpAssEle.GetEta(), kfpAssEle.GetPhi());
double invMassEE = RecoDecay::m(std::array{child1, child2}, std::array{o2::constants::physics::MassElectron, o2::constants::physics::MassElectron});

// reco by KFparticle
Expand Down Expand Up @@ -556,7 +558,7 @@ struct HfTaskElectronWeakBoson {
int trackCount = getIsolatedTrack(track.eta(), track.phi(), track.pt(), tracks) - 1;

if (match.track_as<TrackEle>().pt() > ptTHnThresh && isTHnElectron) {
registry.fill(HIST("hTHnElectrons"), match.track_as<TrackEle>().pt(), match.track_as<TrackEle>().tpcNSigmaEl(), m02Emc, eop, isoEnergy, trackCount);
registry.fill(HIST("hTHnElectrons"), match.track_as<TrackEle>().pt(), match.track_as<TrackEle>().tpcNSigmaEl(), m02Emc, eop, isoEnergy, trackCount, track.eta(), track.tpcSignal());
}
// LOG(info) << "E/p" << eop;
registry.fill(HIST("hEopNsigTPC"), match.track_as<TrackEle>().tpcNSigmaEl(), eop);
Expand Down
Loading