Skip to content

[PWGUD] Modify FwdMuonsUPC.cxx #10691

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

Closed
ariffero wants to merge 11 commits into from
Closed

[PWGUD] Modify FwdMuonsUPC.cxx #10691

Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
e6f0540
Add finer pt binning histo
Feb 7, 2025
3630bee
Merge branch 'AliceO2Group:master' into master
ariffero Mar 6, 2025
34611b4
Modify names of saved gen varaibles.
Mar 6, 2025
96b879f
Please consider the following formatting changes
alibuild Mar 28, 2025
e6c28b0
Merge pull request #11 from alibuild/alibot-cleanup-10691
ariffero Mar 28, 2025
43be741
Merge branch 'AliceO2Group:master' into master
ariffero Apr 3, 2025
64ae9ef
Merge branch 'AliceO2Group:master' into master
ariffero Apr 4, 2025
67f9f95
Merge branch 'AliceO2Group:master' into master
ariffero Apr 7, 2025
af35137
Merge branch 'AliceO2Group:master' into master
ariffero May 15, 2025
4509294
Merge branch 'AliceO2Group:master' into master
ariffero May 16, 2025
b246d90
Merge branch 'AliceO2Group:master' into master
ariffero May 18, 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
46 changes: 30 additions & 16 deletions PWGUD/Tasks/FwdMuonsUPC.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -86,28 +86,28 @@ DECLARE_SOA_TABLE(DiMu, "AOD", "DIMU",
namespace gendimu
{
// dimuon
DECLARE_SOA_COLUMN(M, m, float);
DECLARE_SOA_COLUMN(Pt, pt, float);
DECLARE_SOA_COLUMN(Rap, rap, float);
DECLARE_SOA_COLUMN(Phi, phi, float);
DECLARE_SOA_COLUMN(PhiAv, phiAv, float);
DECLARE_SOA_COLUMN(PhiCh, phiCh, float);
DECLARE_SOA_COLUMN(GenM, genM, float);
DECLARE_SOA_COLUMN(GenPt, genPt, float);
DECLARE_SOA_COLUMN(GenRap, genRap, float);
DECLARE_SOA_COLUMN(GenPhi, genPhi, float);
DECLARE_SOA_COLUMN(GenPhiAv, genPhiAv, float);
DECLARE_SOA_COLUMN(GenPhiCh, genPhiCh, float);
// tracks positive (p) and negative (n)
DECLARE_SOA_COLUMN(Ptp, ptp, float);
DECLARE_SOA_COLUMN(Etap, etap, float);
DECLARE_SOA_COLUMN(Phip, phip, float);
DECLARE_SOA_COLUMN(Ptn, ptn, float);
DECLARE_SOA_COLUMN(Etan, etan, float);
DECLARE_SOA_COLUMN(Phin, phin, float);
DECLARE_SOA_COLUMN(GenPtp, genPtp, float);
DECLARE_SOA_COLUMN(GenEtap, genEtap, float);
DECLARE_SOA_COLUMN(GenPhip, genPhip, float);
DECLARE_SOA_COLUMN(GenPtn, genPtn, float);
DECLARE_SOA_COLUMN(GenEtan, genEtan, float);
DECLARE_SOA_COLUMN(GenPhin, genPhin, float);
} // namespace gendimu

namespace o2::aod
{
DECLARE_SOA_TABLE(GenDimu, "AOD", "GENDIMU",
gendimu::M, gendimu::Pt, gendimu::Rap, gendimu::Phi,
gendimu::PhiAv, gendimu::PhiCh,
gendimu::Ptp, gendimu::Etap, gendimu::Phip,
gendimu::Ptn, gendimu::Etan, gendimu::Phin);
gendimu::GenM, gendimu::GenPt, gendimu::GenRap, gendimu::GenPhi,
gendimu::GenPhiAv, gendimu::GenPhiCh,
gendimu::GenPtp, gendimu::GenEtap, gendimu::GenPhip,
gendimu::GenPtn, gendimu::GenEtan, gendimu::GenPhin);
} // namespace o2::aod

// for saving tree with info on reco MC
Expand Down Expand Up @@ -237,9 +237,19 @@ struct FwdMuonsUPC {
0.60, 0.70, 0.80, 0.90, 1.00, 1.20, 1.40, 1.60, 1.80, 2.00, 2.50,
3.00, 3.50};

std::vector<double> ptFitBinningHalfWidth = {
0.00, 0.005, 0.01, 0.015, 0.02, 0.025, 0.03, 0.035, 0.04, 0.045, 0.05,
0.055, 0.06, 0.065, 0.07, 0.075, 0.08, 0.085, 0.09, 0.095, 0.10,
0.105, 0.11, 0.115, 0.12, 0.125, 0.13, 0.135, 0.14, 0.145, 0.15,
0.1625, 0.175, 0.1875, 0.20, 0.225, 0.25, 0.275, 0.30, 0.35, 0.40,
0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00,
1.10, 1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 2.25,
2.50, 2.75, 3.00, 3.25, 3.50};

// axis
const AxisSpec axisPt{nBinsPt, lowPt, highPt, "#it{p}_{T} GeV/#it{c}"};
const AxisSpec axisPtFit = {ptFitBinning, "#it{p}_{T} (GeV/c)"};
const AxisSpec axisPtFit2 = {ptFitBinningHalfWidth, "#it{p}_{T} (GeV/c)"};
const AxisSpec axisMass{nBinsMass, lowMass, highMass, "m_{#mu#mu} GeV/#it{c}^{2}"};
const AxisSpec axisEta{nBinsEta, lowEta, highEta, "#eta"};
const AxisSpec axisRapidity{nBinsRapidity, lowRapidity, highRapidity, "Rapidity"};
Expand All @@ -256,6 +266,7 @@ struct FwdMuonsUPC {
registry.add("hMass", "Invariant mass of muon pairs;;#counts", kTH1D, {axisMass});
registry.add("hPt", "Transverse momentum of muon pairs;;#counts", kTH1D, {axisPt});
registry.add("hPtFit", "Transverse momentum of muon pairs;;#counts", kTH1D, {axisPtFit});
registry.add("hPtFit2", "Transverse momentum of muon pairs;;#counts", kTH1D, {axisPtFit2});
registry.add("hEta", "Pseudorapidty of muon pairs;;#counts", kTH1D, {axisEta});
registry.add("hRapidity", "Rapidty of muon pairs;;#counts", kTH1D, {axisRapidity});
registry.add("hPhi", "#varphi of muon pairs;;#counts", kTH1D, {axisPhi});
Expand Down Expand Up @@ -314,6 +325,7 @@ struct FwdMuonsUPC {
mcRecoRegistry.add("hMass", "Invariant mass of muon pairs;;#counts", kTH1D, {axisMass});
mcRecoRegistry.add("hPt", "Transverse momentum of muon pairs;;#counts", kTH1D, {axisPt});
mcRecoRegistry.add("hPtFit", "Transverse momentum of muon pairs;;#counts", kTH1D, {axisPtFit});
mcRecoRegistry.add("hPtFit2", "Transverse momentum of muon pairs;;#counts", kTH1D, {axisPtFit2});
mcRecoRegistry.add("hEta", "Pseudorapidty of muon pairs;;#counts", kTH1D, {axisEta});
mcRecoRegistry.add("hRapidity", "Rapidty of muon pairs;;#counts", kTH1D, {axisRapidity});
mcRecoRegistry.add("hPhi", "#varphi of muon pairs;;#counts", kTH1D, {axisPhi});
Expand Down Expand Up @@ -615,6 +627,7 @@ struct FwdMuonsUPC {
registry.fill(HIST("hMass"), p.M());
registry.fill(HIST("hPt"), p.Pt());
registry.fill(HIST("hPtFit"), p.Pt());
registry.fill(HIST("hPtFit2"), p.Pt());
registry.fill(HIST("hEta"), p.Eta());
registry.fill(HIST("hRapidity"), p.Rapidity());
registry.fill(HIST("hPhi"), p.Phi());
Expand Down Expand Up @@ -812,6 +825,7 @@ struct FwdMuonsUPC {
mcRecoRegistry.fill(HIST("hMass"), p.M());
mcRecoRegistry.fill(HIST("hPt"), p.Pt());
mcRecoRegistry.fill(HIST("hPtFit"), p.Pt());
mcRecoRegistry.fill(HIST("hPtFit2"), p.Pt());
mcRecoRegistry.fill(HIST("hEta"), p.Eta());
mcRecoRegistry.fill(HIST("hRapidity"), p.Rapidity());
mcRecoRegistry.fill(HIST("hPhi"), p.Phi());
Expand Down