Skip to content

Commit d935763

Browse files
dsekihatdsekihat
authored
PWGEM/PhotonMeson: update PCM cuts (#4623)
1 parent 88f3cfa commit d935763

File tree

2 files changed

+32
-6
lines changed

2 files changed

+32
-6
lines changed

‎PWGEM/PhotonMeson/Core/CutsLibrary.cxx‎

Lines changed: 31 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -26,6 +26,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
2626
cut->SetMinNCrossedRowsOverFindableClustersTPC(0.8);
2727
cut->SetChi2PerClusterTPC(0.0, 4.0);
2828
cut->SetTPCNsigmaElRange(-3, +3);
29+
cut->SetChi2PerClusterITS(0.0, 5.0);
30+
cut->SetNClustersITS(2, 4);
2931
cut->SetIsWithinBeamPipe(true);
3032
// for v0
3133
cut->SetV0PtRange(0.1f, 1e10f);
@@ -43,6 +45,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
4345
cut->SetMinNCrossedRowsOverFindableClustersTPC(0.8);
4446
cut->SetChi2PerClusterTPC(0.0, 4.0);
4547
cut->SetTPCNsigmaElRange(-3, +3);
48+
cut->SetChi2PerClusterITS(0.0, 5.0);
49+
cut->SetNClustersITS(2, 4);
4650
cut->SetIsWithinBeamPipe(true);
4751
// for v0
4852
cut->SetV0PtRange(0.1f, 1e10f);
@@ -62,12 +66,14 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
6266
cut->SetChi2PerClusterTPC(0.0, 4.0);
6367
cut->SetTPCNsigmaElRange(-3, +3);
6468
cut->SetIsWithinBeamPipe(true);
69+
cut->SetChi2PerClusterITS(0.0, 5.0);
70+
cut->SetNClustersITS(2, 4);
6571
cut->SetRequireITSTPC(true);
6672
// for v0
6773
cut->SetV0PtRange(0.1f, 1e10f);
6874
cut->SetV0EtaRange(-0.9, +0.9);
6975
cut->SetMinCosPA(0.99);
70-
cut->SetMaxPCA(0.5);
76+
cut->SetMaxPCA(1.0);
7177
cut->SetRxyRange(1, 90);
7278
cut->SetAPRange(0.95, 0.01);
7379
return cut;
@@ -78,11 +84,13 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
7884
// cut->SetTrackEtaRange(-0.9, +0.9);
7985
cut->SetIsWithinBeamPipe(true);
8086
cut->SetRequireITSonly(true);
87+
cut->SetChi2PerClusterITS(0.0, 5.0);
88+
cut->SetNClustersITS(4, 4);
8189
// for v0
8290
cut->SetV0PtRange(0.1f, 1e10f);
8391
cut->SetV0EtaRange(-0.9, +0.9);
8492
cut->SetMinCosPA(0.99);
85-
cut->SetMaxPCA(0.5);
93+
cut->SetMaxPCA(1.0);
8694
cut->SetRxyRange(1, 90);
8795
cut->SetAPRange(0.95, 0.01);
8896
return cut;
@@ -115,6 +123,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
115123
cut->SetChi2PerClusterTPC(0.0, 4.0);
116124
cut->SetTPCNsigmaElRange(-3, +3);
117125
cut->SetNClustersITS(2, 4);
126+
cut->SetChi2PerClusterITS(0.0, 5.0);
127+
cut->SetNClustersITS(2, 4);
118128
// cut->SetIsWithinBeamPipe(true);
119129
// for v0
120130
cut->SetV0PtRange(0.1f, 1e10f);
@@ -135,6 +145,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
135145
cut->SetChi2PerClusterTPC(0.0, 4.0);
136146
cut->SetTPCNsigmaElRange(-3, +3);
137147
cut->SetNClustersITS(2, 4);
148+
cut->SetChi2PerClusterITS(0.0, 5.0);
149+
cut->SetNClustersITS(2, 4);
138150
// cut->SetIsWithinBeamPipe(true);
139151
// for v0
140152
cut->SetV0PtRange(0.1f, 1e10f);
@@ -155,6 +167,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
155167
cut->SetChi2PerClusterTPC(0.0, 4.0);
156168
cut->SetTPCNsigmaElRange(-3, +3);
157169
cut->SetNClustersITS(0, 2);
170+
cut->SetChi2PerClusterITS(0.0, 5.0);
171+
cut->SetNClustersITS(2, 2);
158172
// cut->SetIsWithinBeamPipe(true);
159173
// for v0
160174
cut->SetV0PtRange(0.1f, 1e10f);
@@ -176,6 +190,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
176190
cut->SetMinNCrossedRowsOverFindableClustersTPC(0.8);
177191
cut->SetChi2PerClusterTPC(0.0, 4.0);
178192
cut->SetTPCNsigmaElRange(-3, +3);
193+
cut->SetChi2PerClusterITS(0.0, 5.0);
194+
cut->SetNClustersITS(2, 4);
179195
cut->SetIsWithinBeamPipe(true);
180196
// for v0
181197
cut->SetV0PtRange(0.05f, 1e10f);
@@ -194,13 +210,15 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
194210
cut->SetMinNCrossedRowsOverFindableClustersTPC(0.8);
195211
cut->SetChi2PerClusterTPC(0.0, 4.0);
196212
cut->SetTPCNsigmaElRange(-3, +3);
213+
cut->SetChi2PerClusterITS(0.0, 5.0);
214+
cut->SetNClustersITS(2, 4);
197215
cut->SetIsWithinBeamPipe(true);
198216
cut->SetRequireITSTPC(true);
199217
// for v0
200218
cut->SetV0PtRange(0.05f, 1e10f);
201219
cut->SetV0EtaRange(-0.9, +0.9);
202220
cut->SetMinCosPA(0.99);
203-
cut->SetMaxPCA(0.5);
221+
cut->SetMaxPCA(1.0);
204222
cut->SetRxyRange(1, 90);
205223
cut->SetAPRange(0.95, 0.01);
206224
return cut;
@@ -211,11 +229,13 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
211229
// cut->SetTrackEtaRange(-0.9, +0.9);
212230
cut->SetIsWithinBeamPipe(true);
213231
cut->SetRequireITSonly(true);
232+
cut->SetChi2PerClusterITS(0.0, 5.0);
233+
cut->SetNClustersITS(4, 4);
214234
// for v0
215235
cut->SetV0PtRange(0.05f, 1e10f);
216236
cut->SetV0EtaRange(-0.9, +0.9);
217237
cut->SetMinCosPA(0.99);
218-
cut->SetMaxPCA(0.5);
238+
cut->SetMaxPCA(1.0);
219239
cut->SetRxyRange(1, 90);
220240
cut->SetAPRange(0.95, 0.01);
221241
return cut;
@@ -247,6 +267,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
247267
cut->SetMinNCrossedRowsTPC(20);
248268
cut->SetMinNCrossedRowsOverFindableClustersTPC(0.8);
249269
cut->SetChi2PerClusterTPC(0.0, 4.0);
270+
cut->SetChi2PerClusterITS(0.0, 5.0);
271+
cut->SetNClustersITS(2, 4);
250272
cut->SetIsWithinBeamPipe(true);
251273
// for v0
252274
cut->SetV0PtRange(0.1f, 1e10f);
@@ -265,6 +287,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
265287
cut->SetMinNCrossedRowsOverFindableClustersTPC(0.8);
266288
cut->SetChi2PerClusterTPC(0.0, 4.0);
267289
cut->SetTPCNsigmaElRange(-3, +3);
290+
cut->SetChi2PerClusterITS(0.0, 5.0);
291+
cut->SetNClustersITS(2, 4);
268292
cut->SetIsWithinBeamPipe(false);
269293
// for v0
270294
cut->SetV0PtRange(0.1f, 1e10f);
@@ -283,13 +307,15 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
283307
cut->SetMinNCrossedRowsOverFindableClustersTPC(0.8);
284308
cut->SetChi2PerClusterTPC(0.0, 4.0);
285309
cut->SetTPCNsigmaElRange(-3, +3);
310+
cut->SetChi2PerClusterITS(0.0, 5.0);
311+
cut->SetNClustersITS(2, 4);
286312
cut->SetIsWithinBeamPipe(true);
287313
cut->SetRequireITSTPC(true);
288314
// for v0
289315
cut->SetV0PtRange(0.1f, 1e10f);
290316
cut->SetV0EtaRange(-0.9, +0.9);
291317
cut->SetMinCosPA(0.99);
292-
cut->SetMaxPCA(0.5);
318+
cut->SetMaxPCA(1.0);
293319
cut->SetRxyRange(1, 90);
294320
cut->SetAPRange(0.95, 0.01);
295321
return cut;

‎PWGEM/PhotonMeson/Core/HistogramsLibrary.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -49,7 +49,7 @@ void o2::aod::emphotonhistograms::DefineHistograms(THashList* list, const char*
4949
}
5050
if (TString(histClass) == "V0Leg") {
5151
list->Add(new TH1F("hPt", "pT;p_{T} (GeV/c)", 1000, 0.0f, 10));
52-
list->Add(new TH1F("hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", 400, -20, 20));
52+
list->Add(new TH1F("hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", 1000, -50, 50));
5353
list->Add(new TH2F("hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", 180, 0, 2 * M_PI, 40, -2.0f, 2.0f));
5454
list->Add(new TH2F("hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", 200, -50.0f, 50.0f, 200, -50.0f, 50.0f));
5555
list->Add(new TH1F("hNclsTPC", "number of TPC clusters", 161, -0.5, 160.5));

0 commit comments

Comments
 (0)