Skip to content

Commit 0527ec6

Browse files
skundu692skundu692
authored
[PWGLF] remove centrality information from phi meson table producer and add mass information of daughter phi meson in analysis task (#9148)
1 parent 8fa7453 commit 0527ec6

File tree

2 files changed

+39
-17
lines changed

2 files changed

+39
-17
lines changed

‎PWGLF/TableProducer/Resonances/doublephitable.cxx‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -58,8 +58,8 @@ struct doublephitable {
5858

5959
// events
6060
Configurable<float> cfgCutVertex{"cfgCutVertex", 10.0f, "Accepted z-vertex range"};
61-
Configurable<float> cfgCutCentralityMax{"cfgCutCentralityMax", 0.0f, "Accepted maximum Centrality"};
62-
Configurable<float> cfgCutCentralityMin{"cfgCutCentralityMin", 100.0f, "Accepted minimum Centrality"};
61+
// Configurable<float> cfgCutCentralityMax{"cfgCutCentralityMax", 0.0f, "Accepted maximum Centrality"};
62+
// Configurable<float> cfgCutCentralityMin{"cfgCutCentralityMin", 100.0f, "Accepted minimum Centrality"};
6363
// track
6464
Configurable<bool> useGlobalTrack{"useGlobalTrack", true, "use Global track"};
6565
Configurable<float> cfgCutTOFBeta{"cfgCutTOFBeta", 0.0, "cut TOF beta"};
@@ -80,12 +80,12 @@ struct doublephitable {
8080
Configurable<float> maxPhiMass{"maxPhiMass", 1.03, "Maximum phi mass"};
8181

8282
Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex;
83-
Filter centralityFilter = (nabs(aod::cent::centFT0C) < cfgCutCentralityMax && nabs(aod::cent::centFT0C) > cfgCutCentralityMin);
83+
// Filter centralityFilter = (nabs(aod::cent::centFT0C) < cfgCutCentralityMax && nabs(aod::cent::centFT0C) > cfgCutCentralityMin);
8484
Filter acceptanceFilter = (nabs(aod::track::eta) < cfgCutEta && nabs(aod::track::pt) > cfgCutPT);
8585
Filter DCAcutFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz);
8686
Filter PIDcutFilter = nabs(aod::pidtpc::tpcNSigmaKa) < nsigmaCutTPC;
8787

88-
using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0Cs>>;
88+
using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::Mults>>;
8989
using TrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTOFbeta, aod::pidTPCFullKa, aod::pidTOFFullKa>>;
9090

9191
SliceCache cache;

‎PWGLF/Tasks/Resonances/doublephimeson.cxx‎

Lines changed: 35 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -41,6 +41,7 @@ struct doublephimeson {
4141
HistogramRegistry histos{"histos", {}, OutputObjHandlingPolicy::AnalysisObject};
4242

4343
Configurable<bool> fillRotation{"fillRotation", 1, "Fill rotation"};
44+
Configurable<bool> fillPhiMass{"fillPhiMass", 1, "Fill phi mass"};
4445
Configurable<int> strategyPID1{"strategyPID1", 0, "PID strategy 1"};
4546
Configurable<int> strategyPID2{"strategyPID2", 0, "PID strategy 2"};
4647
Configurable<float> minPhiMass{"minPhiMass", 1.01, "Minimum phi mass"};
@@ -78,11 +79,17 @@ struct doublephimeson {
7879
// const AxisSpec thnAxisKstar{configThnAxisKstar, "#it{k}^{*} (GeV/#it{c})"};
7980
const AxisSpec thnAxisDeltaR{configThnAxisDeltaR, "#Delta R)"};
8081
const AxisSpec thnAxisPhiMult{configThnAxisPhiMult, "#Phi Multiplicity)"};
81-
82-
histos.add("SEMassUnlike", "SEMassUnlike", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisDaughterPt, thnAxisDaughterPt, thnAxisDeltaR, thnAxisPhiMult});
83-
histos.add("SEMassRot", "SEMassRot", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisDaughterPt, thnAxisDaughterPt, thnAxisDeltaR, thnAxisPhiMult});
84-
histos.add("MEMassUnlike", "MEMassUnlike", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisDaughterPt, thnAxisDaughterPt, thnAxisDeltaR, thnAxisPhiMult});
8582
histos.add("SEMass", "SEMass", HistType::kTHnSparseF, {thnAxisInvMassPhi, thnAxisInvMassPhi, thnAxisDaughterPt, thnAxisDaughterPt});
83+
if (!fillPhiMass) {
84+
histos.add("SEMassUnlike", "SEMassUnlike", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisDaughterPt, thnAxisDaughterPt, thnAxisDeltaR, thnAxisPhiMult});
85+
histos.add("SEMassRot", "SEMassRot", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisDaughterPt, thnAxisDaughterPt, thnAxisDeltaR, thnAxisPhiMult});
86+
histos.add("MEMassUnlike", "MEMassUnlike", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisDaughterPt, thnAxisDaughterPt, thnAxisDeltaR, thnAxisPhiMult});
87+
}
88+
if (fillPhiMass) {
89+
histos.add("SEMassUnlike", "SEMassUnlike", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisInvMassPhi, thnAxisInvMassPhi, thnAxisDeltaR, thnAxisPhiMult});
90+
histos.add("SEMassRot", "SEMassRot", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisInvMassPhi, thnAxisInvMassPhi, thnAxisDeltaR, thnAxisPhiMult});
91+
histos.add("MEMassUnlike", "MEMassUnlike", HistType::kTHnSparseF, {thnAxisInvMass, thnAxisPt, thnAxisInvMassPhi, thnAxisInvMassPhi, thnAxisDeltaR, thnAxisPhiMult});
92+
}
8693
}
8794

8895
// get kstar
@@ -122,30 +129,30 @@ struct doublephimeson {
122129
}
123130
}
124131
if (PIDStrategy == 1) {
125-
if (ptcand < 0.6) {
132+
if (ptcand < 0.5) {
126133
if (TOFHit != 1 && TMath::Abs(nsigmaTPC) < cutNsigmaTPC) {
127134
return true;
128135
}
129136
if (TOFHit == 1 && TMath::Abs(nsigmaTOF) < cutNsigmaTOF) {
130137
return true;
131138
}
132139
}
133-
if (ptcand >= 0.6) {
140+
if (ptcand >= 0.5) {
134141
if (TOFHit == 1 && TMath::Abs(nsigmaTOF) < cutNsigmaTOF) {
135142
return true;
136143
}
137144
}
138145
}
139146
if (PIDStrategy == 2) {
140-
if (ptcand < 0.6) {
147+
if (ptcand < 0.5) {
141148
if (TOFHit != 1 && TMath::Abs(nsigmaTPC) < cutNsigmaTPC) {
142149
return true;
143150
}
144151
if (TOFHit == 1 && TMath::Abs(nsigmaTOF) < cutNsigmaTOF) {
145152
return true;
146153
}
147154
}
148-
if (ptcand >= 0.6 && ptcand < 1.2) {
155+
if (ptcand >= 0.5 && ptcand < 1.2) {
149156
if (TOFHit == 1 && TMath::Abs(nsigmaTOF) < cutNsigmaTOF) {
150157
return true;
151158
}
@@ -163,15 +170,15 @@ struct doublephimeson {
163170
}
164171
}
165172
if (PIDStrategy == 3) {
166-
if (ptcand < 0.6) {
173+
if (ptcand < 0.5) {
167174
if (TOFHit != 1 && TMath::Abs(nsigmaTPC) < cutNsigmaTPC) {
168175
return true;
169176
}
170177
if (TOFHit == 1 && TMath::Abs(nsigmaTOF) < cutNsigmaTOF) {
171178
return true;
172179
}
173180
}
174-
if (ptcand >= 0.6 && ptcand < 1.2) {
181+
if (ptcand >= 0.5 && ptcand < 1.2) {
175182
if (TOFHit == 1 && TMath::Abs(nsigmaTOF) < cutNsigmaTOF) {
176183
return true;
177184
}
@@ -243,7 +250,12 @@ struct doublephimeson {
243250
exotic = Phid1 + Phid2;
244251
// auto kstar = getkstar(Phid1, Phid2);
245252
auto deltaR = TMath::Sqrt(TMath::Power(Phid1.Phi() - Phid2.Phi(), 2.0) + TMath::Power(Phid1.Eta() - Phid2.Eta(), 2.0));
246-
histos.fill(HIST("SEMassUnlike"), exotic.M(), exotic.Pt(), Phid1.Pt(), Phid2.Pt(), deltaR, phimult);
253+
if (!fillPhiMass) {
254+
histos.fill(HIST("SEMassUnlike"), exotic.M(), exotic.Pt(), Phid1.Pt(), Phid2.Pt(), deltaR, phimult);
255+
}
256+
if (fillPhiMass) {
257+
histos.fill(HIST("SEMassUnlike"), exotic.M(), exotic.Pt(), Phid1.M(), Phid2.M(), deltaR, phimult);
258+
}
247259
histos.fill(HIST("SEMass"), Phid1.M(), Phid2.M(), Phid1.Pt(), Phid2.Pt());
248260
if (fillRotation) {
249261
for (int nrotbkg = 0; nrotbkg < 5; nrotbkg++) {
@@ -257,7 +269,12 @@ struct doublephimeson {
257269
exoticRot = Phid1Rot + Phid2;
258270
// auto kstar_rot = getkstar(Phid1Rot, Phid2);
259271
auto deltaR_rot = TMath::Sqrt(TMath::Power(Phid1Rot.Phi() - Phid2.Phi(), 2.0) + TMath::Power(Phid1Rot.Eta() - Phid2.Eta(), 2.0));
260-
histos.fill(HIST("SEMassRot"), exoticRot.M(), exoticRot.Pt(), Phid1Rot.Pt(), Phid2.Pt(), deltaR_rot, phimult);
272+
if (!fillPhiMass) {
273+
histos.fill(HIST("SEMassRot"), exoticRot.M(), exoticRot.Pt(), Phid1Rot.Pt(), Phid2.Pt(), deltaR_rot, phimult);
274+
}
275+
if (fillPhiMass) {
276+
histos.fill(HIST("SEMassRot"), exoticRot.M(), exoticRot.Pt(), Phid1Rot.M(), Phid2.M(), deltaR_rot, phimult);
277+
}
261278
}
262279
}
263280
}
@@ -311,7 +328,12 @@ struct doublephimeson {
311328
exotic = Phid1 + Phid2;
312329
// auto kstar = getkstar(Phid1, Phid2);
313330
auto deltaR = TMath::Sqrt(TMath::Power(Phid1.Phi() - Phid2.Phi(), 2.0) + TMath::Power(Phid1.Eta() - Phid2.Eta(), 2.0));
314-
histos.fill(HIST("MEMassUnlike"), exotic.M(), exotic.Pt(), Phid1.Pt(), Phid2.Pt(), deltaR, phimult);
331+
if (!fillPhiMass) {
332+
histos.fill(HIST("MEMassUnlike"), exotic.M(), exotic.Pt(), Phid1.Pt(), Phid2.Pt(), deltaR, phimult);
333+
}
334+
if (fillPhiMass) {
335+
histos.fill(HIST("MEMassUnlike"), exotic.M(), exotic.Pt(), Phid1.M(), Phid2.M(), deltaR, phimult);
336+
}
315337
}
316338
}
317339
}

0 commit comments

Comments
 (0)