@@ -107,7 +107,6 @@ struct FlowPbpbPikp {
107107 O2_DEFINE_CONFIGURABLE (cfgUseStrictPID, bool , true , " Use strict PID cuts for TPC" )
108108 O2_DEFINE_CONFIGURABLE (cfgV0AT0Acut, int , 5 , " V0AT0A cut" )
109109 O2_DEFINE_CONFIGURABLE (cfgUseAsymmetricPID, bool , false , " Use asymmetric PID cuts" );
110-
111110
112111 Configurable<std::vector<double >> cfgTrackDensityP0{" cfgTrackDensityP0" , std::vector<double >{0.7217476707 , 0.7384792571 , 0.7542625668 , 0.7640680200 , 0.7701951667 , 0.7755299053 , 0.7805901710 , 0.7849446786 , 0.7957356586 , 0.8113039262 , 0.8211968966 , 0.8280558878 , 0.8329342135 }, " parameter 0 for track density efficiency correction" };
113112 Configurable<std::vector<double >> cfgTrackDensityP1{" cfgTrackDensityP1" , std::vector<double >{-2.169488e-05 , -2.191913e-05 , -2.295484e-05 , -2.556538e-05 , -2.754463e-05 , -2.816832e-05 , -2.846502e-05 , -2.843857e-05 , -2.705974e-05 , -2.477018e-05 , -2.321730e-05 , -2.203315e-05 , -2.109474e-05 }, " parameter 1 for track density efficiency correction" };
@@ -135,7 +134,7 @@ struct FlowPbpbPikp {
135134 std::vector<double > tpcNsigmaCut = cfgTofNsigmaCut;
136135
137136 Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex;
138- Filter trackFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz) && (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtPOIMin) && (aod::track::pt < cfgCutPtPOIMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t ) true )) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls);
137+ Filter trackFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz) && (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtPOIMin) && (aod::track::pt < cfgCutPtPOIMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t )true )) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls);
139138
140139 using AodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::TPCMults, aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0Cs, aod::CentFT0As, aod::Mults>>;
141140 using AodTracksWithoutBayes = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA, aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr, aod::pidTOFbeta, aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>>;
@@ -440,7 +439,6 @@ struct FlowPbpbPikp {
440439 return 0 ;
441440 }
442441 histos.fill (HIST (" hTrackCount" ), 6.5 ); // After TPC crossed rows selection
443-
444442 return 1 ;
445443 }
446444
@@ -514,21 +512,21 @@ struct FlowPbpbPikp {
514512 }
515513
516514 template <typename TTrack>
517- int getNsigmaPIDTpcTof_assymmetric (TTrack track)
515+ int getNsigmaPIDTpcTof_Assymmetric (TTrack track)
518516 {
519517 // Computing Nsigma arrays for pion, kaon, and protons
520518 std::array<float , 3 > nSigmaTPC = {track.tpcNSigmaPi (), track.tpcNSigmaKa (), track.tpcNSigmaPr ()};
521519 std::array<float , 3 > nSigmaTOF = {track.tofNSigmaPi (), track.tofNSigmaKa (), track.tofNSigmaPr ()};
522520 int pid = -1 ;
523521
524522 bool isPion, isKaon, isProton;
525- bool isTpcPion = nSigmaTPC[0 ] < tpcNsigmaCut[0 ] && nSigmaTPC[0 ] > tpcNsigmaCut[0 + 3 ];
526- bool isTpcKaon = nSigmaTPC[1 ] < tpcNsigmaCut[1 ] && nSigmaTPC[1 ] > tpcNsigmaCut[1 + 3 ];
527- bool isTpcProton = nSigmaTPC[2 ] < tpcNsigmaCut[2 ] && nSigmaTPC[2 ] > tpcNsigmaCut[2 + 3 ];
523+ bool isTpcPion = nSigmaTPC[0 ] < tpcNsigmaCut[0 ] && nSigmaTPC[0 ] > tpcNsigmaCut[0 + 3 ];
524+ bool isTpcKaon = nSigmaTPC[1 ] < tpcNsigmaCut[1 ] && nSigmaTPC[1 ] > tpcNsigmaCut[1 + 3 ];
525+ bool isTpcProton = nSigmaTPC[2 ] < tpcNsigmaCut[2 ] && nSigmaTPC[2 ] > tpcNsigmaCut[2 + 3 ];
528526
529- bool isTofPion = nSigmaTOF[0 ] < tofNsigmaCut[0 ] && nSigmaTOF[0 ] > tofNsigmaCut[0 + 3 ];
530- bool isTofKaon = nSigmaTOF[1 ] < tofNsigmaCut[1 ] && nSigmaTOF[1 ] > tofNsigmaCut[1 + 3 ];
531- bool isTofProton = nSigmaTOF[2 ] < tofNsigmaCut[2 ] && nSigmaTOF[2 ] > tofNsigmaCut[2 + 3 ];
527+ bool isTofPion = nSigmaTOF[0 ] < tofNsigmaCut[0 ] && nSigmaTOF[0 ] > tofNsigmaCut[0 + 3 ];
528+ bool isTofKaon = nSigmaTOF[1 ] < tofNsigmaCut[1 ] && nSigmaTOF[1 ] > tofNsigmaCut[1 + 3 ];
529+ bool isTofProton = nSigmaTOF[2 ] < tofNsigmaCut[2 ] && nSigmaTOF[2 ] > tofNsigmaCut[2 + 3 ];
532530
533531 if (track.pt () > cfgTofPtCut && !track.hasTOF ()) {
534532 return 0 ;
@@ -922,7 +920,7 @@ struct FlowPbpbPikp {
922920 bool withinPtPOI = (cfgCutPtPOIMin < pt) && (pt < cfgCutPtPOIMax); // within POI pT range
923921 bool withinPtRef = (cfgCutPtMin < pt) && (pt < cfgCutPtMax); // within RF pT range
924922
925- pidIndex = cfgUseAsymmetricPID ? getNsigmaPIDTpcTof_assymmetric (track) : getNsigmaPIDTpcTof (track);
923+ pidIndex = cfgUseAsymmetricPID ? getNsigmaPIDTpcTof_Assymmetric (track) : getNsigmaPIDTpcTof (track);
926924
927925 weff = 1 ; // Initializing weff for each track
928926 // NUA weights
0 commit comments