@@ -104,27 +104,27 @@ struct HadronNucleiCorrelation {
104104 Configurable<double > dcaPar0{" dcaPar0" , 0.004 , " par 0" };
105105 Configurable<double > dcaPar1{" dcaPar1" , 0.013 , " par 1" };
106106 Configurable<bool > doDCAZ{" doDCAZ" , true , " do DCA z cut" };
107- Configurable<int16_t > min_TPC_nClusters{ " min_TPC_nClusters " , 80 , " minimum number of found TPC clusters" };
108- Configurable<float > min_TPC_nCrossedRowsOverFindableCls{ " min_TPC_nCrossedRowsOverFindableCls " , 0.8 , " n TPC Crossed Rows Over Findable Cls" };
109- Configurable<float > max_chi2_TPC{ " max_chi2_TPC " , 4 .0f , " maximum TPC chi^2/Ncls" };
110- Configurable<float > max_chi2_ITS{ " max_chi2_ITS " , 36 .0f , " maximum ITS chi^2/Ncls" };
107+ Configurable<int16_t > minTPCnClusters{ " minTPCnClusters " , 80 , " minimum number of found TPC clusters" };
108+ Configurable<float > minTPCnCrossedRowsOverFindableCls{ " minTPCnCrossedRowsOverFindableCls " , 0.8 , " n TPC Crossed Rows Over Findable Cls" };
109+ Configurable<float > maxchi2TPC{ " maxchi2TPC " , 4 .0f , " maximum TPC chi^2/Ncls" };
110+ Configurable<float > maxchi2ITS{ " maxchi2ITS " , 36 .0f , " maximum ITS chi^2/Ncls" };
111111 Configurable<float > etaCut{" etaCut" , 0 .8f , " eta cut" };
112- Configurable<float > max_DCAxy{ " max_DCAxy " , 0 .14f , " Maximum DCAxy" };
113- Configurable<float > max_DCAz{ " max_DCAz " , 0 .1f , " Maximum DCAz" };
112+ Configurable<float > maxDCAxy{ " maxDCAxy " , 0 .14f , " Maximum DCAxy" };
113+ Configurable<float > maxDCAz{ " maxDCAz " , 0 .1f , " Maximum DCAz" };
114114 Configurable<float > nsigmaTPC{" nsigmaTPC" , 3 .0f , " cut nsigma TPC" };
115115 Configurable<float > nsigmaElPr{" nsigmaElPr" , 1 .0f , " cut nsigma TPC El for protons" };
116116 Configurable<float > nsigmaElDe{" nsigmaElDe" , 3 .0f , " cut nsigma TPC El for protons" };
117117 Configurable<float > nsigmaTOF{" nsigmaTOF" , 3 .5f , " cut nsigma TOF" };
118118 Configurable<float > nsigmaITSPr{" nsigmaITSPr" , -2 .0f , " cut nsigma ITS Pr" };
119119 Configurable<float > nsigmaITSDe{" nsigmaITSDe" , -2 .0f , " cut nsigma ITS De" };
120120 Configurable<bool > doITSPID{" doITSPID" , true , " do ITS PID" };
121- Configurable<float > pTthrpr_TOF{ " pTthrpr_TOF " , 0 .8f , " threshold pT proton to use TOF" };
122- Configurable<float > pTthrpr_TPCEl{ " pTthrpr_TPCEl " , 1 .0f , " threshold pT proton to use TPC El rejection" };
123- Configurable<float > pTthrde_TOF{ " pTthrde_TOF " , 1 .0f , " threshold pT deuteron to use TOF" };
124- Configurable<float > pTthrde_TPCEl{ " pTthrde_TPCEl " , 1 .0f , " threshold pT deuteron to use TPC El rejection" };
121+ Configurable<float > pTthrprTOF{ " pTthrprTOF " , 0 .8f , " threshold pT proton to use TOF" };
122+ Configurable<float > pTthrprTPCEl{ " pTthrprTPCEl " , 1 .0f , " threshold pT proton to use TPC El rejection" };
123+ Configurable<float > pTthrdeTOF{ " pTthrdeTOF " , 1 .0f , " threshold pT deuteron to use TOF" };
124+ Configurable<float > pTthrdeTPCEl{ " pTthrdeTPCEl " , 1 .0f , " threshold pT deuteron to use TPC El rejection" };
125125 Configurable<bool > rejectionEl{" rejectionEl" , true , " use TPC El rejection" };
126- Configurable<float > max_tpcSharedCls{ " max_tpcSharedCls " , 0.4 , " maximum fraction of TPC shared clasters" };
127- Configurable<int > min_itsNCls{ " min_itsNCls " , 0 , " minimum allowed number of ITS clasters" };
126+ Configurable<float > maxtpcSharedCls{ " maxtpcSharedCls " , 0.4 , " maximum fraction of TPC shared clasters" };
127+ Configurable<int > minitsNCls{ " minitsNCls " , 0 , " minimum allowed number of ITS clasters" };
128128 Configurable<int > maxmixcollsGen{" maxmixcollsGen" , 100 , " maxmixcollsGen" };
129129 Configurable<float > radiusTPC{" radiusTPC" , 1.2 , " TPC radius to calculate phi_star for" };
130130 Configurable<float > dEta{" dEta" , 0.01 , " minimum allowed difference in eta between two tracks in a pair" };
@@ -421,12 +421,12 @@ struct HadronNucleiCorrelation {
421421
422422 // Filters
423423 Filter vertexFilter = nabs(o2::aod::singletrackselector::posZ) <= cutzVertex;
424- Filter trackFilter = o2::aod::singletrackselector::tpcNClsFound >= min_TPC_nClusters &&
425- o2::aod::singletrackselector::unPack<singletrackselector::binning::chi2>(o2::aod::singletrackselector::storedTpcChi2NCl) <= max_chi2_TPC &&
426- o2::aod::singletrackselector::unPack<singletrackselector::binning::rowsOverFindable>(o2::aod::singletrackselector::storedTpcCrossedRowsOverFindableCls) >= min_TPC_nCrossedRowsOverFindableCls &&
427- o2::aod::singletrackselector::unPack<singletrackselector::binning::chi2>(o2::aod::singletrackselector::storedItsChi2NCl) <= max_chi2_ITS &&
428- nabs (o2::aod::singletrackselector::unPack<singletrackselector::binning::dca>(o2::aod::singletrackselector::storedDcaXY)) <= max_DCAxy &&
429- nabs (o2::aod::singletrackselector::unPack<singletrackselector::binning::dca>(o2::aod::singletrackselector::storedDcaXY)) <= max_DCAz &&
424+ Filter trackFilter = o2::aod::singletrackselector::tpcNClsFound >= minTPCnClusters &&
425+ o2::aod::singletrackselector::unPack<singletrackselector::binning::chi2>(o2::aod::singletrackselector::storedTpcChi2NCl) <= maxchi2TPC &&
426+ o2::aod::singletrackselector::unPack<singletrackselector::binning::rowsOverFindable>(o2::aod::singletrackselector::storedTpcCrossedRowsOverFindableCls) >= minTPCnCrossedRowsOverFindableCls &&
427+ o2::aod::singletrackselector::unPack<singletrackselector::binning::chi2>(o2::aod::singletrackselector::storedItsChi2NCl) <= maxchi2ITS &&
428+ nabs (o2::aod::singletrackselector::unPack<singletrackselector::binning::dca>(o2::aod::singletrackselector::storedDcaXY)) <= maxDCAxy &&
429+ nabs (o2::aod::singletrackselector::unPack<singletrackselector::binning::dca>(o2::aod::singletrackselector::storedDcaXY)) <= maxDCAz &&
430430 nabs (o2::aod::singletrackselector::eta) <= etaCut;
431431
432432 Filter simvertexFilter = nabs(o2::aod::mccollision::posZ) <= cutzVertex;
@@ -437,11 +437,11 @@ struct HadronNucleiCorrelation {
437437 bool isProton = false ;
438438 bool isTPCPID = std::abs (track.tpcNSigmaPr ()) < nsigmaTPC;
439439 bool isTOFPID = std::abs (track.tofNSigmaPr ()) < nsigmaTOF;
440- bool isTPCElRejection = rejectionEl && track.beta () < betahasTOFthr && track.pt () < pTthrpr_TPCEl && track.tpcNSigmaEl () >= nsigmaElPr;
440+ bool isTPCElRejection = rejectionEl && track.beta () < betahasTOFthr && track.pt () < pTthrprTPCEl && track.tpcNSigmaEl () >= nsigmaElPr;
441441 bool isITSPID = track.itsNSigmaPr () > nsigmaITSPr;
442442
443443 if (isTPCPID) {
444- if (track.pt () < pTthrpr_TOF ) {
444+ if (track.pt () < pTthrprTOF ) {
445445 if (!doITSPID || isITSPID) {
446446 if (sign > 0 ) {
447447 if (track.sign () > 0 ) {
@@ -496,11 +496,11 @@ struct HadronNucleiCorrelation {
496496 bool isDeuteron = false ;
497497 bool isTPCPID = std::abs (track.tpcNSigmaDe ()) < nsigmaTPC;
498498 bool isTOFPID = std::abs (track.tofNSigmaDe ()) < nsigmaTOF;
499- bool isTPCElRejection = rejectionEl && track.beta () < betahasTOFthr && track.pt () < pTthrde_TPCEl && track.tpcNSigmaEl () >= nsigmaElDe;
499+ bool isTPCElRejection = rejectionEl && track.beta () < betahasTOFthr && track.pt () < pTthrdeTPCEl && track.tpcNSigmaEl () >= nsigmaElDe;
500500 bool isITSPID = track.itsNSigmaDe () > nsigmaITSDe;
501501
502502 if (isTPCPID) {
503- if (track.pt () < pTthrde_TOF ) {
503+ if (track.pt () < pTthrdeTOF ) {
504504 if (!doITSPID || isITSPID) {
505505 if (sign > 0 ) {
506506 if (track.sign () > 0 ) {
@@ -733,9 +733,9 @@ struct HadronNucleiCorrelation {
733733 if (removeSameBunchPileup && !track.template singleCollSel_as <soa::Filtered<FilteredCollisions>>().isNoSameBunchPileup ())
734734 continue ;
735735
736- if (track.tpcFractionSharedCls () > max_tpcSharedCls )
736+ if (track.tpcFractionSharedCls () > maxtpcSharedCls )
737737 continue ;
738- if (track.itsNCls () < min_itsNCls )
738+ if (track.itsNCls () < minitsNCls )
739739 continue ;
740740
741741 if (IsProton (track, +1 ))
@@ -809,13 +809,13 @@ struct HadronNucleiCorrelation {
809809 if (removeSameBunchPileup && !part0.template singleCollSel_as <soa::Filtered<FilteredCollisions>>().isNoSameBunchPileup ())
810810 continue ;
811811
812- if (part0.tpcFractionSharedCls () > max_tpcSharedCls )
812+ if (part0.tpcFractionSharedCls () > maxtpcSharedCls )
813813 continue ;
814- if (part0.itsNCls () < min_itsNCls )
814+ if (part0.itsNCls () < minitsNCls )
815815 continue ;
816- if (part1.tpcFractionSharedCls () > max_tpcSharedCls )
816+ if (part1.tpcFractionSharedCls () > maxtpcSharedCls )
817817 continue ;
818- if (part1.itsNCls () < min_itsNCls )
818+ if (part1.itsNCls () < minitsNCls )
819819 continue ;
820820
821821 if (!applyDCAcut (part0))
@@ -854,13 +854,13 @@ struct HadronNucleiCorrelation {
854854 if (removeSameBunchPileup && !part0.template singleCollSel_as <soa::Filtered<FilteredCollisions>>().isNoSameBunchPileup ())
855855 continue ;
856856
857- if (part0.tpcFractionSharedCls () > max_tpcSharedCls )
857+ if (part0.tpcFractionSharedCls () > maxtpcSharedCls )
858858 continue ;
859- if (part0.itsNCls () < min_itsNCls )
859+ if (part0.itsNCls () < minitsNCls )
860860 continue ;
861- if (part1.tpcFractionSharedCls () > max_tpcSharedCls )
861+ if (part1.tpcFractionSharedCls () > maxtpcSharedCls )
862862 continue ;
863- if (part1.itsNCls () < min_itsNCls )
863+ if (part1.itsNCls () < minitsNCls )
864864 continue ;
865865
866866 if (!applyDCAcut (part0))
@@ -945,13 +945,13 @@ struct HadronNucleiCorrelation {
945945 if (removeSameBunchPileup && !part1.template singleCollSel_as <soa::Filtered<FilteredCollisions>>().isNoSameBunchPileup ())
946946 continue ;
947947
948- if (part0.tpcFractionSharedCls () > max_tpcSharedCls )
948+ if (part0.tpcFractionSharedCls () > maxtpcSharedCls )
949949 continue ;
950- if (part0.itsNCls () < min_itsNCls )
950+ if (part0.itsNCls () < minitsNCls )
951951 continue ;
952- if (part1.tpcFractionSharedCls () > max_tpcSharedCls )
952+ if (part1.tpcFractionSharedCls () > maxtpcSharedCls )
953953 continue ;
954- if (part1.itsNCls () < min_itsNCls )
954+ if (part1.itsNCls () < minitsNCls )
955955 continue ;
956956
957957 if (!applyDCAcut (part0))
@@ -1031,9 +1031,9 @@ struct HadronNucleiCorrelation {
10311031 if (std::abs (track.template singleCollSel_as <FilteredCollisions>().posZ ()) > cutzVertex)
10321032 continue ;
10331033
1034- if (track.tpcFractionSharedCls () > max_tpcSharedCls )
1034+ if (track.tpcFractionSharedCls () > maxtpcSharedCls )
10351035 continue ;
1036- if (track.itsNCls () < min_itsNCls )
1036+ if (track.itsNCls () < minitsNCls )
10371037 continue ;
10381038
10391039 if (IsProton (track, +1 ) && track.pdgCode () == PDG_t::kProton ) {
0 commit comments