@@ -45,7 +45,7 @@ struct ThreeParticleCorrelations {
4545
4646 // Track PID parameters
4747 double pionID = 0.0 , kaonID = 1.0 , protonID = 2.0 ;
48- float nSigma0 = 0.0 , nSigma2 = 2.0 , nSigma4 = 4.0 , nSigma5 = 5.0 ;
48+ float nSigma0 = 0.0 , nSigma1 = 1.0 , nSigma2 = 2.0 , nSigma4 = 4.0 , nSigma5 = 5.0 ;
4949
5050 // V0 filter parameters
5151 float tpcNCrossedRowsMin = 70.0 ;
@@ -55,8 +55,8 @@ struct ThreeParticleCorrelations {
5555 int dcaV0DauMax = 1 ;
5656
5757 // Track filter parameters
58- float pionPtMin = 0.3 , pionPtMax = 2.3 , kaonPtMin = 0.5 , kaonPtMax = 2.5 , protonPtMin = 0.5 , protonPtMax = 2.5 ;
59- float pionPtMid = 1.5 , kaonPtMid1 = 1.5 , kaonPtMid2 = 2.0 , protonPtMid = 0.7 ;
58+ float pionPtMin = 0.3 , pionPtMax = 2.3 , kaonPtMin = 0.5 , kaonPtMax = 2.3 , protonPtMin = 0.6 ;
59+ float pionPtMid1 = 1.6 , pionPtMid2 = 2.0 , kaonPtMid1 = 1.5 , kaonPtMid2 = 2.0 , protonPtMid = 2.3 ;
6060
6161 // RD filter parameters
6262 float dEtaMax = 0.05 , dEtaMin = 0.023 ;
@@ -1002,11 +1002,15 @@ struct ThreeParticleCorrelations {
10021002 }
10031003 if (track.pt () < pionPtMin) {
10041004 return false ;
1005- } else if (track.pt () > pionPtMin && track.pt () < pionPtMid ) {
1005+ } else if (track.pt () > pionPtMin && track.pt () < pionPtMid1 ) {
10061006 if (std::abs (track.tofNSigmaPi ()) >= nSigma4) {
10071007 return false ;
10081008 }
1009- } else if (track.pt () > pionPtMid && track.pt () < pionPtMax) {
1009+ } else if (track.pt () > pionPtMid1 && track.pt () < pionPtMid2) {
1010+ if (track.tofNSigmaPi () <= -nSigma4 || track.tofNSigmaPi () >= -nSigma2) {
1011+ return false ;
1012+ }
1013+ } else if (track.pt () > pionPtMid2 && track.pt () < pionPtMax) {
10101014 if (track.tofNSigmaPi () <= -nSigma4 || track.tofNSigmaPi () >= nSigma0) {
10111015 return false ;
10121016 }
@@ -1025,7 +1029,7 @@ struct ThreeParticleCorrelations {
10251029 return false ;
10261030 }
10271031 } else if (track.pt () > kaonPtMid1 && track.pt () < kaonPtMid2) {
1028- if (track.tofNSigmaKa () <= -nSigma2 || track.tofNSigmaKa () >= nSigma4) {
1032+ if (track.tofNSigmaKa () <= -nSigma1 || track.tofNSigmaKa () >= nSigma4) {
10291033 return false ;
10301034 }
10311035 } else if (track.pt () > kaonPtMid2 && track.pt () < kaonPtMax) {
@@ -1043,14 +1047,10 @@ struct ThreeParticleCorrelations {
10431047 if (track.pt () < protonPtMin) {
10441048 return false ;
10451049 } else if (track.pt () > protonPtMin && track.pt () < protonPtMid) {
1046- if (track.tofNSigmaPr () <= -nSigma2 || track.tofNSigmaPr () >= nSigma4) {
1047- return false ;
1048- }
1049- } else if (track.pt () > protonPtMid && track.pt () < protonPtMax) {
10501050 if (std::abs (track.tofNSigmaPr ()) >= nSigma4) {
10511051 return false ;
10521052 }
1053- } else if (track.pt () > protonPtMax ) {
1053+ } else if (track.pt () > protonPtMid ) {
10541054 if (track.tofNSigmaPr () <= -nSigma2 || track.tofNSigmaPr () >= nSigma4) {
10551055 return false ;
10561056 }
0 commit comments