Skip to content

Commit 207398b

Browse files
ananthapadmanabhan18ananthapadmanabhan18
authored
[PWGUD] Applied pT and y cuts to invariant mass distributions (AliceO2Group#9013)
1 parent e92e3ed commit 207398b

File tree

1 file changed

+63
-25
lines changed

1 file changed

+63
-25
lines changed

‎PWGUD/Tasks/exclusiveRhoTo4Pi.cxx‎

Lines changed: 63 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -82,27 +82,35 @@ struct UPCAnalysis {
8282
histos.add("tofNSigmaPi_WTS_PID_Pi", "TOF nSigma Pion with track selection and PID Selection of Pi; Entries", kTH1F, {{100, -15, 15}});
8383

8484
// Track Transverse Momentum
85-
histos.add("pT_track_WOTS", "pT without track selection; pT [GeV/c]; Events", kTH1F, {{100, 0, 2}});
86-
histos.add("pT_track_WTS", "pT with track selection; pT [GeV/c]; Events", kTH1F, {{100, 0, 2}});
87-
histos.add("pT_track_WTS_PID_Pi", "pT with track selection and PID selection of Pi; pT [GeV/c]; Events", kTH1F, {{100, 0, 2}});
85+
histos.add("pT_track_WOTS", "pT without track selection; pT [GeV/c]; Events", kTH1F, {{1000, 0, 2}});
86+
histos.add("pT_track_WTS", "pT with track selection; pT [GeV/c]; Events", kTH1F, {{1000, 0, 2}});
87+
histos.add("pT_track_WTS_PID_Pi", "pT with track selection and PID selection of Pi; pT [GeV/c]; Events", kTH1F, {{1000, 0, 2}});
8888

8989
// Zero charge Event Transverse Momentum
90-
histos.add("pT_event_0charge_WTS_PID_Pi", "Event pT in 0 Charge Events With Track Selection and PID Selection of Pi; pT [GeV/c]; Counts", kTH1F, {{100, 0, 2}});
90+
histos.add("pT_event_0charge_WTS_PID_Pi", "Event pT in 0 Charge Events With Track Selection and PID Selection of Pi; pT [GeV/c]; Counts", kTH1F, {{1000, 0, 2}});
9191

9292
// Non Zero charge Event Transverse Momentum
9393
histos.add("pT_event_non0charge_WTS_PID_Pi", "Event pT in Non 0 Charge Events With Track Selection and PID Selection of Pi; pT [GeV/c]; Counts", kTH1F, {{100, 0, 2}});
9494

9595
// Rapidity of 0 charge Events
96-
histos.add("rapidity_event_0charge_WTS_PID_Pi", "Rapidity of Events With Track Selection and PID Selection of Pi; y; Counts", kTH1F, {{100, -2.5, 2.5}});
96+
histos.add("rapidity_event_0charge_WTS_PID_Pi_domainA", "Rapidity of Events With Track Selection and PID Selection of Pi for p_{T} < 0.15 GeV/c; y; Counts", kTH1F, {{1000, -2.5, 2.5}});
97+
histos.add("rapidity_event_0charge_WTS_PID_Pi_domainB", "Rapidity of Events With Track Selection and PID Selection of Pi for 0.15< p_{T} < 0.80 GeV/c; y; Counts", kTH1F, {{1000, -2.5, 2.5}});
98+
histos.add("rapidity_event_0charge_WTS_PID_Pi_domainC", "Rapidity of Events With Track Selection and PID Selection of Pi for p_{T} > 0.80 GeV/c; y; Counts", kTH1F, {{1000, -2.5, 2.5}});
9799

98100
// Rapidity of non 0 charge Events
99-
histos.add("rapidity_event_non0charge_WTS_PID_Pi", "Rapidity of Events With Track Selection and PID Selection of Pi; y; Counts", kTH1F, {{100, -2.5, 2.5}});
101+
histos.add("rapidity_event_non0charge_WTS_PID_Pi_domainA", "Rapidity of Events With Track Selection and PID Selection of Pi for p_{T} < 0.15 GeV/c; y; Counts", kTH1F, {{1000, -2.5, 2.5}});
102+
histos.add("rapidity_event_non0charge_WTS_PID_Pi_domainB", "Rapidity of Events With Track Selection and PID Selection of Pi for 0.15< p_{T} < 0.80 GeV/c$; y; Counts", kTH1F, {{1000, -2.5, 2.5}});
103+
histos.add("rapidity_event_non0charge_WTS_PID_Pi_domainC", "Rapidity of Events With Track Selection and PID Selection of Pi for p_{T} > 0.80 GeV/c; y; Counts", kTH1F, {{1000, -2.5, 2.5}});
100104

101105
// Invariant Mass of 0 charge events
102-
histos.add("invMass_event_0charge_WTS_PID_Pi", "Invariant Mass Distribution of 0 charge Events with PID Selection of Pi; m(#pi^{+}#pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}});
106+
histos.add("invMass_event_0charge_WTS_PID_Pi_domainA", "Invariant Mass Distribution of 0 charge Events with PID Selection of Pi for p_{T} < 0.15 GeV/c; m(#pi^{+}#pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}}); // pT < 0.15GeV
107+
histos.add("invMass_event_0charge_WTS_PID_Pi_domainB", "Invariant Mass Distribution of 0 charge Events with PID Selection of Pi for 0.15< p_{T} < 0.80 GeV/c; m(#pi^{+}#pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}}); // 0.15GeV < pT < 0.8GeV
108+
histos.add("invMass_event_0charge_WTS_PID_Pi_domainC", "Invariant Mass Distribution of 0 charge Events with PID Selection of Pi for p_{T} > 0.80 GeV/c; m(#pi^{+}#pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}}); // 0.8GeV < pT
103109

104110
// Invariant mass of non 0 charge events
105-
histos.add("invMass_event_non0charge_WTS_PID_Pi", "Invariant Mass Distribution of non 0 charge Events with PID Selection of Pi; m(#pi^{+}#Pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}});
111+
histos.add("invMass_event_non0charge_WTS_PID_Pi_domainA", "Invariant Mass Distribution of non 0 charge Events with PID Selection of Pi for p_{T} < 0.15 GeV/c; m(#pi^{+}#pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}}); // pT < 0.15GeV
112+
histos.add("invMass_event_non0charge_WTS_PID_Pi_domainB", "Invariant Mass Distribution of non 0 charge Events with PID Selection of Pi for 0.15< p_{T} < 0.80 GeV/c; m(#pi^{+}#pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}}); // 0.15GeV < pT < 0.8GeV
113+
histos.add("invMass_event_non0charge_WTS_PID_Pi_domainC", "Invariant Mass Distribution of non 0 charge Events with PID Selection of Pi for p_{T} > 0.80 GeV/c; m(#pi^{+}#pi^{-}#pi^{+}#pi^{-}) [GeV/c]", kTH1F, {{1000, 0.8, 2.5}}); // 0.8GeV < pT
106114

107115
// tpc signal
108116
histos.add("tpcSignal", "TPC dEdx vs p; p [GeV/c]; dEdx [a.u.]", kTH2F, {{500, 0, 10}, {5000, 0.0, 5000.0}});
@@ -120,8 +128,10 @@ struct UPCAnalysis {
120128
histos.add("V0A", "V0A amplitude", kTH1F, {{1000, 0.0, 100}});
121129

122130
// Collin Soper Theta and Phi
123-
histos.add("CS_phi", "#phi Distribution; #phi; Entries", kTH1F, {{200, -3.2, 3.2}});
124-
histos.add("CS_costheta", "#theta Distribution;cos(#theta); Entries", kTH1F, {{200, -1, 1}});
131+
histos.add("CS_phi_pair_1", "#phi Distribution; #phi; Entries", kTH1F, {{200, -3.2, 3.2}});
132+
histos.add("CS_phi_pair_2", "#phi Distribution; #phi; Entries", kTH1F, {{200, -3.2, 3.2}});
133+
histos.add("CS_costheta_pair_1", "#theta Distribution;cos(#theta); Entries", kTH1F, {{200, -1, 1}});
134+
histos.add("CS_costheta_pair_2", "#theta Distribution;cos(#theta); Entries", kTH1F, {{200, -1, 1}});
125135

126136
} // End of init function
127137
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------
@@ -284,6 +294,7 @@ struct UPCAnalysis {
284294
return;
285295
}
286296

297+
// Selecting Events with net charge = 0
287298
if (len_Pi_minus == 2 && len_Pi_plus == 2) {
288299

289300
TLorentzVector p1, p2, p3, p4, p1234;
@@ -307,21 +318,35 @@ struct UPCAnalysis {
307318
k23 = k2 + k3;
308319
k24 = k2 + k4;
309320

310-
auto phi_pair_1 = PhiCollinsSoperFrame(k13, k24, k1234);
311-
auto phi_pair_2 = PhiCollinsSoperFrame(k14, k23, k1234);
312-
auto cos_theta_1 = CosThetaCollinsSoperFrame(k13, k24, k1234);
313-
auto cos_theta_2 = CosThetaCollinsSoperFrame(k14, k23, k1234);
314-
315-
histos.fill(HIST("pT_event_0charge_WTS_PID_Pi"), p1234.Pt());
316-
histos.fill(HIST("rapidity_event_0charge_WTS_PID_Pi"), p1234.Rapidity());
317-
histos.fill(HIST("invMass_event_0charge_WTS_PID_Pi"), p1234.M());
318-
histos.fill(HIST("CS_phi"), phi_pair_1);
319-
histos.fill(HIST("CS_phi"), phi_pair_2);
320-
histos.fill(HIST("CS_costheta"), cos_theta_1);
321-
histos.fill(HIST("CS_costheta"), cos_theta_2);
321+
if (fabs(p1234.Rapidity()) < 0.5) {
322+
histos.fill(HIST("pT_event_0charge_WTS_PID_Pi"), p1234.Pt());
323+
if (p1234.Pt() < 0.15) {
324+
histos.fill(HIST("rapidity_event_0charge_WTS_PID_Pi_domainA"), p1234.Rapidity());
325+
histos.fill(HIST("invMass_event_0charge_WTS_PID_Pi_domainA"), p1234.M());
326+
327+
auto phi_pair_1 = PhiCollinsSoperFrame(k13, k24, k1234);
328+
auto phi_pair_2 = PhiCollinsSoperFrame(k14, k23, k1234);
329+
auto cos_theta_1 = CosThetaCollinsSoperFrame(k13, k24, k1234);
330+
auto cos_theta_2 = CosThetaCollinsSoperFrame(k14, k23, k1234);
331+
332+
histos.fill(HIST("CS_phi_pair_1"), phi_pair_1);
333+
histos.fill(HIST("CS_phi_pair_2"), phi_pair_2);
334+
histos.fill(HIST("CS_costheta_pair_1"), cos_theta_1);
335+
histos.fill(HIST("CS_costheta_pair_2"), cos_theta_2);
336+
}
337+
if (p1234.Pt() > 0.15 && p1234.Pt() < 0.80) {
338+
histos.fill(HIST("rapidity_event_0charge_WTS_PID_Pi_domainB"), p1234.Rapidity());
339+
histos.fill(HIST("invMass_event_0charge_WTS_PID_Pi_domainB"), p1234.M());
340+
}
341+
if (p1234.Pt() > 0.80) {
342+
histos.fill(HIST("rapidity_event_0charge_WTS_PID_Pi_domainC"), p1234.Rapidity());
343+
histos.fill(HIST("invMass_event_0charge_WTS_PID_Pi_domainC"), p1234.M());
344+
}
345+
} // End of Rapidity range selection
322346

323347
} // End of Analysis for 0 charge events
324348

349+
// Selecting Events with net charge != 0 for estimation of background
325350
if (len_Pi_minus != 2 && len_Pi_plus != 2) {
326351

327352
TLorentzVector p1, p2, p3, p4, p1234;
@@ -332,9 +357,22 @@ struct UPCAnalysis {
332357

333358
p1234 = p1 + p2 + p3 + p4;
334359

335-
histos.fill(HIST("pT_event_non0charge_WTS_PID_Pi"), p1234.Pt());
336-
histos.fill(HIST("rapidity_event_non0charge_WTS_PID_Pi"), p1234.Rapidity());
337-
histos.fill(HIST("invMass_event_non0charge_WTS_PID_Pi"), p1234.M());
360+
if (fabs(p1234.Rapidity()) < 0.5) {
361+
histos.fill(HIST("pT_event_non0charge_WTS_PID_Pi"), p1234.Pt());
362+
363+
if (p1234.Pt() < 0.15) {
364+
histos.fill(HIST("rapidity_event_non0charge_WTS_PID_Pi_domainA"), p1234.Rapidity());
365+
histos.fill(HIST("invMass_event_non0charge_WTS_PID_Pi_domainA"), p1234.M());
366+
}
367+
if (p1234.Pt() > 0.15 && p1234.Pt() < 0.80) {
368+
histos.fill(HIST("rapidity_event_non0charge_WTS_PID_Pi_domainB"), p1234.Rapidity());
369+
histos.fill(HIST("invMass_event_non0charge_WTS_PID_Pi_domainB"), p1234.M());
370+
}
371+
if (p1234.Pt() > 0.80) {
372+
histos.fill(HIST("rapidity_event_non0charge_WTS_PID_Pi_domainC"), p1234.Rapidity());
373+
histos.fill(HIST("invMass_event_non0charge_WTS_PID_Pi_domainC"), p1234.M());
374+
}
375+
} // End of Rapidity range selection
338376

339377
} // End of Analysis for non 0 charge events
340378

0 commit comments

Comments
 (0)