Skip to content

Commit 838920d

Browse files
hkoleyhkoley
authored
PWGLF: Update TPCNClsFound variable (#5226)
1 parent 61caa22 commit 838920d

1 file changed

Lines changed: 24 additions & 14 deletions

File tree

‎PWGLF/Tasks/Resonances/lambda1520analysis.cxx‎

Lines changed: 24 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -34,19 +34,21 @@ struct lambda1520analysis {
3434

3535
// Configurables
3636
// switches
37-
Configurable<bool> isEtaAssym{"isEtaAssym", false, "Turn on/off EtaAssym calculation"};
37+
Configurable<bool> cEtaAssym{"cEtaAssym", false, "Turn on/off EtaAssym calculation"};
3838
// Configurable<bool> isFillQA{"isFillQA", false, "Turn on/off QA plots"};
39-
Configurable<bool> IsAddlTrackcut{"IsAddlTrackcut", false, "Switch to turn on/off Additional track cut"};
40-
Configurable<bool> IsOldPIDcut{"IsOldPIDcut", false, "Switch to turn on/off old PID cut to apply pt dependent cut"};
41-
Configurable<bool> IsDCAr7SigCut{"IsDCAr7SigCut", false, "Track DCAr 7 Sigma cut to PV Maximum"};
39+
Configurable<bool> cAddlTrackcut{"cAddlTrackcut", false, "Switch to turn on/off Additional track cut"};
40+
Configurable<bool> cOldPIDcut{"cOldPIDcut", false, "Switch to turn on/off old PID cut to apply pt dependent cut"};
41+
Configurable<bool> cDCAr7SigCut{"cDCAr7SigCut", false, "Track DCAr 7 Sigma cut to PV Maximum"};
4242
Configurable<bool> cKinCuts{"cKinCuts", false, "Kinematic Cuts for p-K pair opening angle"};
43+
Configurable<bool> cTPCNClsFound{"cTPCNClsFound", false, "Switch to turn on/off TPCNClsFound cut"};
4344

4445
// Pre-selection Track cuts
4546
Configurable<float> cMinPtcut{"cMinPtcut", 0.15f, "Minimal pT for tracks"};
4647
Configurable<int> cMinTPCncr{"cMinTPCncr", 70, "Minimum number of TPC X rows"};
4748
Configurable<float> cMinRtpccut{"cMinRtpccut", 0.8f, "minimum ratio of number of Xrows to findable clusters in TPC"};
4849
Configurable<float> cMaxChi2ITScut{"cMaxChi2ITScut", 36.0f, "Maximal pT for Chi2/cluster for ITS"};
4950
Configurable<float> cMaxChi2TPCcut{"cMaxChi2TPCcut", 4.0f, "Maximal pT for Chi2/cluster for TPC"};
51+
Configurable<float> cMinTPCNClsFound{"cMinTPCNClsFound", 120, "minimum TPCNClsFound value for good track"};
5052

5153
// DCA Selections
5254
// DCAr to PV
@@ -186,7 +188,7 @@ struct lambda1520analysis {
186188
histos.add("Result/Data/h3lambda1520invmassME", "Invariant mass of #Lambda(1520) mixed event K^{#pm}p^{#mp}", HistType::kTH3F, {axisMult, axisPt, axisMassLambda1520});
187189
}
188190

189-
if (isEtaAssym) {
191+
if (cEtaAssym) {
190192
histos.add("Result/Data/hlambda1520invmassUnlikeSignAside", "Invariant mass of #Lambda(1520) Unlike Sign A side", {HistType::kTH1F, {axisMassLambda1520}});
191193
histos.add("Result/Data/hlambda1520invmassLikeSignAside", "Invariant mass of #Lambda(1520) Like Sign A side", {HistType::kTH1F, {axisMassLambda1520}});
192194
histos.add("Result/Data/hlambda1520invmassUnlikeSignCside", "Invariant mass of #Lambda(1520) Unlike Sign C side", {HistType::kTH1F, {axisMassLambda1520}});
@@ -239,7 +241,7 @@ struct lambda1520analysis {
239241
// basic track cuts
240242
if (std::abs(track.pt()) < cMinPtcut)
241243
return false;
242-
if (IsDCAr7SigCut) {
244+
if (cDCAr7SigCut) {
243245
if (std::abs(track.dcaXY()) > (0.004f + 0.0130f / (track.pt()))) // 7 - Sigma cut
244246
return false;
245247
} else {
@@ -250,7 +252,7 @@ struct lambda1520analysis {
250252
return false;
251253
if (track.tpcNClsCrossedRows() < cMinTPCncr)
252254
return false;
253-
if (IsAddlTrackcut) {
255+
if (cAddlTrackcut) {
254256
if (!track.passedITSRefit() || !track.passedTPCRefit())
255257
return false;
256258
if (track.tpcCrossedRowsOverFindableCls() < cMinRtpccut)
@@ -260,6 +262,8 @@ struct lambda1520analysis {
260262
if (track.tpcChi2NCl() > cMaxChi2TPCcut)
261263
return false;
262264
}
265+
if (cTPCNClsFound && (track.tpcNClsFound() < cMinTPCNClsFound))
266+
return false;
263267
if (cfgPrimaryTrack && !track.isPrimaryTrack())
264268
return false;
265269
if (cfgGlobalWoDCATrack && !track.isGlobalTrackWoDCA())
@@ -469,7 +473,7 @@ struct lambda1520analysis {
469473
continue;
470474
if (cUseOnlyTOFTrackKa && !isTrk2hasTOF)
471475
continue;
472-
if (IsOldPIDcut) {
476+
if (cOldPIDcut) {
473477
if (!selectionoldPIDProton(trk1) || !selectionoldPIDKaon(trk2))
474478
continue;
475479
} else {
@@ -511,6 +515,12 @@ struct lambda1520analysis {
511515
histos.fill(HIST("QA/QAafter/Kaon/eta"), trk2.eta());
512516
}
513517

518+
// TPCncluster distributions
519+
histos.fill(HIST("TPCncluster/TPCnclusterpr"), trk1.tpcNClsFound());
520+
histos.fill(HIST("TPCncluster/TPCnclusterka"), trk2.tpcNClsFound());
521+
histos.fill(HIST("TPCncluster/TPCnclusterPhipr"), trk1.tpcNClsFound(), trk1.phi());
522+
histos.fill(HIST("TPCncluster/TPCnclusterPhika"), trk2.tpcNClsFound(), trk2.phi());
523+
514524
// Multiplicity correlation calibration plots
515525
int counterglo1 = 0;
516526
int counterglo2 = 0;
@@ -552,20 +562,20 @@ struct lambda1520analysis {
552562
if constexpr (!IsMix) {
553563
histos.fill(HIST("Result/Data/lambda1520invmass"), lResonance.M());
554564
histos.fill(HIST("Result/Data/h3lambda1520invmass"), collision.cent(), lResonance.Pt(), lResonance.M());
555-
if (isEtaAssym && trk1.eta() > 0.2 && trk1.eta() < 0.8 && trk2.eta() > 0.2 && trk2.eta() < 0.8) { // Eta-range will be updated
565+
if (cEtaAssym && trk1.eta() > 0.2 && trk1.eta() < 0.8 && trk2.eta() > 0.2 && trk2.eta() < 0.8) { // Eta-range will be updated
556566
histos.fill(HIST("Result/Data/hlambda1520invmassUnlikeSignAside"), lResonance.M());
557567
histos.fill(HIST("Result/Data/h3lambda1520invmassUnlikeSignAside"), collision.cent(), lResonance.Pt(), lResonance.M());
558-
} else if (isEtaAssym && trk1.eta() > -0.6 && trk1.eta() < 0.0 && trk2.eta() > -0.6 && trk2.eta() < 0.0) { // Eta-range will be updated
568+
} else if (cEtaAssym && trk1.eta() > -0.6 && trk1.eta() < 0.0 && trk2.eta() > -0.6 && trk2.eta() < 0.0) { // Eta-range will be updated
559569
histos.fill(HIST("Result/Data/hlambda1520invmassUnlikeSignCside"), lResonance.M());
560570
histos.fill(HIST("Result/Data/h3lambda1520invmassUnlikeSignCside"), collision.cent(), lResonance.Pt(), lResonance.M());
561571
}
562572
} else {
563573
histos.fill(HIST("Result/Data/lambda1520invmassME"), lResonance.M());
564574
histos.fill(HIST("Result/Data/h3lambda1520invmassME"), collision.cent(), lResonance.Pt(), lResonance.M());
565-
if (isEtaAssym && trk1.eta() > 0.2 && trk1.eta() < 0.8 && trk2.eta() > 0.2 && trk2.eta() < 0.8) { // Eta-range will be updated
575+
if (cEtaAssym && trk1.eta() > 0.2 && trk1.eta() < 0.8 && trk2.eta() > 0.2 && trk2.eta() < 0.8) { // Eta-range will be updated
566576
histos.fill(HIST("Result/Data/hlambda1520invmassMixedAside"), lResonance.M());
567577
histos.fill(HIST("Result/Data/h3lambda1520invmassMixedAside"), collision.cent(), lResonance.Pt(), lResonance.M());
568-
} else if (isEtaAssym && trk1.eta() > -0.6 && trk1.eta() < 0.0 && trk2.eta() > -0.6 && trk2.eta() < 0.0) { // Eta-range will be updated
578+
} else if (cEtaAssym && trk1.eta() > -0.6 && trk1.eta() < 0.0 && trk2.eta() > -0.6 && trk2.eta() < 0.0) { // Eta-range will be updated
569579
histos.fill(HIST("Result/Data/hlambda1520invmassMixedCside"), lResonance.M());
570580
histos.fill(HIST("Result/Data/h3lambda1520invmassMixedCside"), collision.cent(), lResonance.Pt(), lResonance.M());
571581
}
@@ -601,10 +611,10 @@ struct lambda1520analysis {
601611
}
602612
} else {
603613
if constexpr (!IsMix) {
604-
if (isEtaAssym && trk1.eta() > 0.2 && trk1.eta() < 0.8 && trk2.eta() > 0.2 && trk2.eta() < 0.8) { // Eta-range will be updated
614+
if (cEtaAssym && trk1.eta() > 0.2 && trk1.eta() < 0.8 && trk2.eta() > 0.2 && trk2.eta() < 0.8) { // Eta-range will be updated
605615
histos.fill(HIST("Result/Data/hlambda1520invmassLikeSignAside"), lResonance.M());
606616
histos.fill(HIST("Result/Data/h3lambda1520invmassLikeSignAside"), collision.cent(), lResonance.Pt(), lResonance.M());
607-
} else if (isEtaAssym && trk1.eta() > -0.6 && trk1.eta() < 0.0 && trk2.eta() > -0.6 && trk2.eta() < 0.0) { // Eta-range will be updated
617+
} else if (cEtaAssym && trk1.eta() > -0.6 && trk1.eta() < 0.0 && trk2.eta() > -0.6 && trk2.eta() < 0.0) { // Eta-range will be updated
608618
histos.fill(HIST("Result/Data/hlambda1520invmassLikeSignCside"), lResonance.M());
609619
histos.fill(HIST("Result/Data/h3lambda1520invmassLikeSignCside"), collision.cent(), lResonance.Pt(), lResonance.M());
610620
}

0 commit comments

Comments
 (0)