Skip to content

Commit ee9b4a6

Browse files
rolavickrolavick
authored
Fixing few more bugs
1 parent dbb4f76 commit ee9b4a6

File tree

1 file changed

+30
-32
lines changed

1 file changed

+30
-32
lines changed

‎PWGUD/Tasks/sgExclusiveJpsiMidrapidity.cxx‎

Lines changed: 30 additions & 32 deletions
Original file line numberDiff line numberDiff line change
@@ -41,14 +41,14 @@ struct SgExclusiveJpsiMidrapidity {
4141
SGSelector sgSelector;
4242

4343
// Number of Selection cuts
44-
const int numSelectionCuts = 16;
44+
static constexpr int numSelectionCuts = 16;
4545

4646
// Gapside rejection
47-
const int gapSideLow = 0;
48-
const int gapSideHigh = 2;
47+
static constexpr int gapSideLow = 0;
48+
static constexpr int gapSideHigh = 2;
4949

5050
// Numbers for selections
51-
const int two = 2;
51+
static constexpr int two = 2;
5252

5353
Configurable<float> fv0Cut{"fv0Cut", 100., "fv0aThreshold"};
5454
Configurable<float> ft0aCut{"ft0aCut", 100., "ft0aThreshold"};
@@ -174,7 +174,6 @@ struct SgExclusiveJpsiMidrapidity {
174174
//____________________________________________________________________________________
175175

176176
// Create LorentzVector to store all tracks, Pion tracks and TPC Pion PID
177-
std::vector<LorentzVector> allTracks;
178177
std::vector<LorentzVector> onlyPionTracks;
179178
std::vector<float> onlyPionSigma;
180179
std::vector<decltype(tracks.begin())> rawPionTracks;
@@ -186,7 +185,7 @@ struct SgExclusiveJpsiMidrapidity {
186185

187186
for (const auto& t : tracks) {
188187
// Apply good track selection criteria
189-
if (!trackselector(t, parameters))
188+
if (!sgSelector.trackselector(t, parameters))
190189
continue;
191190

192191
double dEdx = t.tpcSignal();
@@ -195,7 +194,6 @@ struct SgExclusiveJpsiMidrapidity {
195194

196195
// Create Lorentz vector for this track (use constructor, portable)
197196
LorentzVector a(t.px(), t.py(), t.pz(), o2::constants::physics::MassPionCharged);
198-
allTracks.push_back(a);
199197

200198
// Apply TPC pion sigma
201199
auto nSigmaPi = t.tpcNSigmaPi();
@@ -251,46 +249,46 @@ struct SgExclusiveJpsiMidrapidity {
251249
if (rawPionTracks[0].sign() != rawPionTracks[1].sign()) {
252250

253251
registry.fill(HIST("hSelectionCounter"), 7);
254-
registry.fill(HIST("Pion/hMassUnlike"), p.M());
252+
registry.fill(HIST("TwoPion/hMassUnlike"), p.M());
255253

256254
// Flexible mass limits, can be selected in the configurable
257255
if ((p.M() > massMin) && (p.M() < massMax)) {
258256

259257
registry.fill(HIST("hSelectionCounter"), 8);
260258

261-
registry.fill(HIST("Pion/hPt"), p.Pt());
262-
registry.fill(HIST("Pion/hEta"), p.Eta());
263-
registry.fill(HIST("Pion/hRap"), p.Rapidity());
264-
registry.fill(HIST("Pion/hPhiSystem"), p.Phi());
265-
registry.fill(HIST("Pion/hMPt"), p.M(), p.Pt());
259+
registry.fill(HIST("TwoPion/hPt"), p.Pt());
260+
registry.fill(HIST("TwoPion/hEta"), p.Eta());
261+
registry.fill(HIST("TwoPion/hRap"), p.Rapidity());
262+
registry.fill(HIST("TwoPion/hPhiSystem"), p.Phi());
263+
registry.fill(HIST("TwoPion/hMPt"), p.M(), p.Pt());
266264

267265
// flexible pt limit for selecting coherent Rho(0)
268266
if (p.Pt() < ptCoherent) {
269267

270268
registry.fill(HIST("hSelectionCounter"), 9);
271269

272270
// Quality Control plots after coherent Rho(0) selection
273-
registry.fill(HIST("Pion/hEta_t1"), onlyPionTracks[0].Eta());
274-
registry.fill(HIST("Pion/hEta_t2"), onlyPionTracks[1].Eta());
275-
registry.fill(HIST("Pion/hPtsingle_track1"), onlyPionTracks[0].Pt());
276-
registry.fill(HIST("Pion/hPtsingle_track2"), onlyPionTracks[1].Pt());
277-
278-
registry.fill(HIST("Pion/hNsigMuvsPt1"), onlyPionTracks[0].Pt(), rawPionTracks[0].tpcNSigmaPi());
279-
registry.fill(HIST("Pion/hNsigMuvsPt2"), onlyPionTracks[1].Pt(), rawPionTracks[1].tpcNSigmaPi());
280-
registry.fill(HIST("Pion/hNsigElvsPt1"), onlyPionTracks[0].Pt(), rawPionTracks[0].tpcNSigmaEl());
281-
registry.fill(HIST("Pion/hNsigElvsPt2"), onlyPionTracks[1].Pt(), rawPionTracks[1].tpcNSigmaEl());
282-
registry.fill(HIST("Pion/hNsigEl1vsEl2"), rawPionTracks[0].tpcNSigmaPi(), rawPionTracks[1].tpcNSigmaPi());
283-
284-
registry.fill(HIST("Pion/hP1"), onlyPionTracks[0].P(), rawPionTracks[0].tpcSignal());
285-
registry.fill(HIST("Pion/hP2"), onlyPionTracks[1].P(), rawPionTracks[1].tpcSignal());
286-
registry.fill(HIST("Pion/hTPCsig1"), rawPionTracks[0].tpcSignal(), rawPionTracks[1].tpcSignal());
287-
288-
registry.fill(HIST("Pion/Coherent/hMassUnlikeCoherent"), p.M());
271+
registry.fill(HIST("TwoPion/hEta_t1"), onlyPionTracks[0].Eta());
272+
registry.fill(HIST("TwoPion/hEta_t2"), onlyPionTracks[1].Eta());
273+
registry.fill(HIST("TwoPion/hPtsingle_track1"), onlyPionTracks[0].Pt());
274+
registry.fill(HIST("TwoPion/hPtsingle_track2"), onlyPionTracks[1].Pt());
275+
276+
registry.fill(HIST("TwoPion/hNsigMuvsPt1"), onlyPionTracks[0].Pt(), rawPionTracks[0].tpcNSigmaPi());
277+
registry.fill(HIST("TwoPion/hNsigMuvsPt2"), onlyPionTracks[1].Pt(), rawPionTracks[1].tpcNSigmaPi());
278+
registry.fill(HIST("TwoPion/hNsigElvsPt1"), onlyPionTracks[0].Pt(), rawPionTracks[0].tpcNSigmaEl());
279+
registry.fill(HIST("TwoPion/hNsigElvsPt2"), onlyPionTracks[1].Pt(), rawPionTracks[1].tpcNSigmaEl());
280+
registry.fill(HIST("TwoPion/hNsigEl1vsEl2"), rawPionTracks[0].tpcNSigmaPi(), rawPionTracks[1].tpcNSigmaPi());
281+
282+
registry.fill(HIST("TwoPion/hP1"), onlyPionTracks[0].P(), rawPionTracks[0].tpcSignal());
283+
registry.fill(HIST("TwoPion/hP2"), onlyPionTracks[1].P(), rawPionTracks[1].tpcSignal());
284+
registry.fill(HIST("TwoPion/hTPCsig1"), rawPionTracks[0].tpcSignal(), rawPionTracks[1].tpcSignal());
285+
286+
registry.fill(HIST("TwoPion/Coherent/hMassUnlikeCoherent"), p.M());
289287
}
290288
// Incoherent Rho(0) selection
291289
if (p.Pt() > ptCoherent) {
292290
registry.fill(HIST("hSelectionCounter"), 10);
293-
registry.fill(HIST("Pion/Incoherent/hMassUnlikeInCoherent"), p.M());
291+
registry.fill(HIST("TwoPion/Incoherent/hMassUnlikeInCoherent"), p.M());
294292
}
295293
}
296294
}
@@ -299,13 +297,13 @@ struct SgExclusiveJpsiMidrapidity {
299297
if (rawPionTracks[0].sign() == rawPionTracks[1].sign()) {
300298

301299
registry.fill(HIST("hSelectionCounter"), 11);
302-
registry.fill(HIST("Pion/hMassLike"), p.M());
300+
registry.fill(HIST("TwoPion/hMassLike"), p.M());
303301

304302
// Mass limit
305303
if ((p.M() > massMin) && (p.M() < massMax)) {
306304

307305
registry.fill(HIST("hSelectionCounter"), 12);
308-
registry.fill(HIST("Pion/hPtLike"), p.Pt());
306+
registry.fill(HIST("TwoPion/hPtLike"), p.Pt());
309307

310308
// Coherent Rho(0) selection
311309
if (p.Pt() < ptCoherent) {

0 commit comments

Comments
 (0)