@@ -61,6 +61,12 @@ struct femtoUniversePairTaskTrackPhi {
6161 SliceCache cache;
6262 Preslice<FemtoFullParticles> perCol = aod::femtouniverseparticle::fdCollisionId;
6363
64+ Configurable<bool > ConfIsCPR{" ConfIsCPR" , true , " Close Pair Rejection" };
65+ Configurable<bool > ConfCPRPlotPerRadii{" ConfCPRPlotPerRadii" , false , " Plot CPR per radii" };
66+ Configurable<float > ConfCPRdeltaPhiCut{" ConfCPRdeltaPhiCut" , 0.01 , " Delta Phi cut for Close Pair Rejection" };
67+ Configurable<float > ConfCPRdeltaEtaCut{" ConfCPRdeltaEtaCut" , 0.004 , " Delta Eta cut for Close Pair Rejection" };
68+ Configurable<float > ConfCPRChosenRadii{" ConfCPRChosenRadii" , 0.80 , " Delta Eta cut for Close Pair Rejection" };
69+
6470 // / Table for both particles
6571 struct : o2::framework::ConfigurableGroup {
6672 Configurable<float > ConfNsigmaCombinedProton{" ConfNsigmaCombinedProton" , 3.0 , " TPC and TOF Proton Sigma (combined) for momentum > 0.5" };
@@ -81,7 +87,6 @@ struct femtoUniversePairTaskTrackPhi {
8187 struct : o2::framework::ConfigurableGroup {
8288 Configurable<bool > ConfIsSame{" ConfIsSame" , false , " Pairs of the same particle" };
8389 Configurable<int > ConfPDGCodeTrack{" ConfPDGCodeTrack" , 2212 , " Particle 2 - PDG code" };
84- // Configurable<uint32_t> ConfCutTrack{"ConfCutTrack", 5542474, "Particle 2 - Selection bit"};
8590 Configurable<int > ConfPIDTrack{" ConfPIDTrack" , 2 , " Particle 2 - Read from cutCulator" }; // we also need the possibility to specify whether the bit is true/false ->std>>vector<std::pair<int, int>>
8691 Configurable<int > ConfTrackSign{" ConfTrackSign" , 1 , " Track sign" };
8792 Configurable<bool > ConfIsTrackIdentified{" ConfIsTrackIdentified" , true , " Enable PID for the track" };
@@ -90,10 +95,7 @@ struct femtoUniversePairTaskTrackPhi {
9095 // / Particle 2 --- PHI
9196 struct : o2::framework::ConfigurableGroup {
9297 Configurable<int > ConfPDGCodePhi{" ConfPDGCodePhi" , 333 , " Phi meson - PDG code" };
93- // Configurable<uint32_t> ConfCutPhi{"ConfCutPhi", 5542474, "Phi meson - Selection bit from cutCulator"};
9498 Configurable<int > ConfPIDPhi{" ConfPIDPhi" , 2 , " Phi meson - Read from cutCulator" }; // we also need the possibility to specify whether the bit is true/false ->std>>vector<std::pair<int, int>>int>>
95- // Configurable<float> ConfPtLowPhi{"ConfPtLowPhi", 0.5, "Lower limit for Pt for the Phi meson"}; // change according to wrzesa cuts
96- // Configurable<float> ConfPtHighPhi{"ConfPtHighPhi", 1.5, "Higher limit for Pt for the Phi meson"};
9799 } ConfPhi;
98100
99101 // / Partitions for particle 1
@@ -128,7 +130,6 @@ struct femtoUniversePairTaskTrackPhi {
128130
129131 // / Correlation part
130132 ConfigurableAxis ConfMultBins{" ConfMultBins" , {VARIABLE_WIDTH, 0 .0f , 4 .0f , 8 .0f , 12 .0f , 16 .0f , 20 .0f , 24 .0f , 28 .0f , 32 .0f , 36 .0f , 40 .0f , 44 .0f , 48 .0f , 52 .0f , 56 .0f , 60 .0f , 64 .0f , 68 .0f , 72 .0f , 76 .0f , 80 .0f , 84 .0f , 88 .0f , 92 .0f , 96 .0f , 100 .0f , 200 .0f , 99999 .f }, " Mixing bins - multiplicity" }; // \todo to be obtained from the hash task
131- // ConfigurableAxis ConfMultBins{"CfgMultBins", {VARIABLE_WIDTH, 0.0f, 20.0f, 40.0f, 60.0f, 80.0f, 100.0f, 200.0f, 99999.f}, "Mixing bins - multiplicity"};
132133 ConfigurableAxis ConfVtxBins{" ConfVtxBins" , {VARIABLE_WIDTH, -10 .0f , -8 .f , -6 .f , -4 .f , -2 .f , 0 .f , 2 .f , 4 .f , 6 .f , 8 .f , 10 .f }, " Mixing bins - z-vertex" };
133134 ConfigurableAxis ConfmTBins3D{" ConfmTBins3D" , {VARIABLE_WIDTH, 1 .02f , 1 .14f , 1 .20f , 1 .26f , 1 .38f , 1 .56f , 1 .86f , 4 .50f }, " mT Binning for the 3Dimensional plot: k* vs multiplicity vs mT (set <<ConfBothTracks.ConfUse3D>> to true in order to use)" };
134135 ConfigurableAxis ConfmultBins3D{" ConfmultBins3D" , {VARIABLE_WIDTH, 0 .0f , 20 .0f , 30 .0f , 40 .0f , 99999 .0f }, " multiplicity Binning for the 3Dimensional plot: k* vs multiplicity vs mT (set <<ConfBothTracks.ConfUse3D>> to true in order to use)" };
@@ -138,11 +139,6 @@ struct femtoUniversePairTaskTrackPhi {
138139 ConfigurableAxis ConfkstarBins{" ConfkstarBins" , {1500 , 0 ., 6 .}, " binning kstar" };
139140 ConfigurableAxis ConfkTBins{" ConfkTBins" , {150 , 0 ., 9 .}, " binning kT" };
140141 ConfigurableAxis ConfmTBins{" ConfmTBins" , {225 , 0 ., 7.5 }, " binning mT" };
141- Configurable<bool > ConfIsCPR{" ConfIsCPR" , true , " Close Pair Rejection" };
142- Configurable<bool > ConfCPRPlotPerRadii{" ConfCPRPlotPerRadii" , false , " Plot CPR per radii" };
143- Configurable<float > ConfCPRdeltaPhiCut{" ConfCPRdeltaPhiCut" , 0.01 , " Delta Phi cut for Close Pair Rejection" };
144- Configurable<float > ConfCPRdeltaEtaCut{" ConfCPRdeltaEtaCut" , 0.004 , " Delta Eta cut for Close Pair Rejection" };
145- Configurable<float > ConfCPRChosenRadii{" ConfCPRChosenRadii" , 0.80 , " Delta Eta cut for Close Pair Rejection" };
146142
147143 FemtoUniverseFemtoContainer<femtoUniverseFemtoContainer::EventType::same, femtoUniverseFemtoContainer::Observable::kstar> sameEventFemtoCont;
148144 FemtoUniverseFemtoContainer<femtoUniverseFemtoContainer::EventType::mixed, femtoUniverseFemtoContainer::Observable::kstar> mixedEventFemtoCont;
@@ -160,13 +156,6 @@ struct femtoUniversePairTaskTrackPhi {
160156 // PID for protons
161157 bool IsProtonNSigma (float mom, float nsigmaTPCPr, float nsigmaTOFPr) // previous version from: https://github.com/alisw/AliPhysics/blob/master/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoMJTrackCut.cxx
162158 {
163- // |nsigma_TPC| < 3 for p < 0.5 GeV/c
164- // |nsigma_combined| < 3 for p > 0.5
165-
166- // using configurables:
167- // ConfNsigmaTPCProton -> TPC Kaon Sigma for momentum < 0.5
168- // ConfNsigmaCombinedProton -> TPC and TOF Kaon Sigma (combined) for momentum > 0.5
169-
170159 if (mom < 0.5 ) {
171160 if (TMath::Abs (nsigmaTPCPr) < ConfBothTracks.ConfNsigmaTPCProton ) {
172161 return true ;
@@ -175,7 +164,6 @@ struct femtoUniversePairTaskTrackPhi {
175164 }
176165 } else if (mom > 0.4 ) {
177166 if (TMath::Hypot (nsigmaTOFPr, nsigmaTPCPr) < ConfBothTracks.ConfNsigmaCombinedProton ) {
178- // if (TMath::Abs(nsigmaTPCPr) < ConfBothTracks.ConfNsigmaCombinedProton) {
179167 return true ;
180168 } else {
181169 return false ;
@@ -225,12 +213,6 @@ struct femtoUniversePairTaskTrackPhi {
225213
226214 bool IsPionNSigma (float mom, float nsigmaTPCPi, float nsigmaTOFPi)
227215 {
228- // |nsigma_TPC| < 3 for p < 0.5 GeV/c
229- // |nsigma_combined| < 3 for p > 0.5
230-
231- // using configurables:
232- // ConfNsigmaTPCPion -> TPC Kaon Sigma for momentum < 0.5
233- // ConfNsigmaCombinedPion -> TPC and TOF Pion Sigma (combined) for momentum > 0.5
234216 if (true ) {
235217 if (mom < 0.5 ) {
236218 if (TMath::Abs (nsigmaTPCPi) < ConfBothTracks.ConfNsigmaTPCPion ) {
@@ -240,7 +222,6 @@ struct femtoUniversePairTaskTrackPhi {
240222 }
241223 } else if (mom > 0.5 ) {
242224 if (TMath::Hypot (nsigmaTOFPi, nsigmaTPCPi) < ConfBothTracks.ConfNsigmaCombinedPion ) {
243- // if (TMath::Abs(nsigmaTPCPi) < ConfBothTracks.ConfNsigmaCombinedPion) {
244225 return true ;
245226 } else {
246227 return false ;
0 commit comments