|
20 | 20 | #include <TH2F.h> |
21 | 21 | #include <TLorentzVector.h> |
22 | 22 | #include <TPDGCode.h> |
| 23 | +#include <string> |
| 24 | +#include <vector> |
23 | 25 | #include <TDatabasePDG.h> |
24 | 26 | #include <cmath> |
25 | 27 | #include <array> |
@@ -83,6 +85,7 @@ struct kstarpbpb { |
83 | 85 | Configurable<float> cfgCutDCAxy{"cfgCutDCAxy", 2.0f, "DCAxy range for tracks"}; |
84 | 86 | Configurable<float> cfgCutDCAz{"cfgCutDCAz", 2.0f, "DCAz range for tracks"}; |
85 | 87 | Configurable<bool> useGlobalTrack{"useGlobalTrack", true, "use Global track"}; |
| 88 | + Configurable<float> nsigmaCutTOF{"nsigmacutTOF", 3.0, "Value of the TOF Nsigma cut"}; |
86 | 89 | Configurable<float> nsigmaCutTPC{"nsigmacutTPC", 3.0, "Value of the TPC Nsigma cut"}; |
87 | 90 | Configurable<float> nsigmaCutCombined{"nsigmaCutCombined", 3.0, "Value of the TOF Nsigma cut"}; |
88 | 91 | Configurable<int> cfgNoMixedEvents{"cfgNoMixedEvents", 1, "Number of mixed events per event"}; |
@@ -269,17 +272,11 @@ struct kstarpbpb { |
269 | 272 | bool selectionPID(const T& candidate, int PID) |
270 | 273 | { |
271 | 274 | if (PID == 0) { |
272 | | - if (!candidate.hasTOF() && TMath::Abs(candidate.tpcNSigmaKa()) < nsigmaCutTPC) { |
273 | | - return true; |
274 | | - } |
275 | | - if (candidate.hasTOF() && ((candidate.tofNSigmaKa() * candidate.tofNSigmaKa()) + (candidate.tpcNSigmaKa() * candidate.tpcNSigmaKa())) < (nsigmaCutCombined * nsigmaCutCombined)) { |
| 275 | + if (candidate.hasTOF() && TMath::Abs(candidate.tofNSigmaKa()) < nsigmaCutTOF) { |
276 | 276 | return true; |
277 | 277 | } |
278 | 278 | } else if (PID == 1) { |
279 | | - if (!candidate.hasTOF() && TMath::Abs(candidate.tpcNSigmaPi()) < nsigmaCutTPC) { |
280 | | - return true; |
281 | | - } |
282 | | - if (candidate.hasTOF() && ((candidate.tofNSigmaPi() * candidate.tofNSigmaPi()) + (candidate.tpcNSigmaPi() * candidate.tpcNSigmaPi())) < (nsigmaCutCombined * nsigmaCutCombined)) { |
| 279 | + if (candidate.hasTOF() && TMath::Abs(candidate.tofNSigmaPi()) < nsigmaCutTOF) { |
283 | 280 | return true; |
284 | 281 | } |
285 | 282 | } |
|
0 commit comments