@@ -203,10 +203,10 @@ struct LfTreeCreatorClusterStudies {
203203 Configurable<float > cascsetting_cascCosPA{" cascsetting_cascCosPA" , 0 .99f , " Minimum cCosine of the pointing angle of the cascade" };
204204 Configurable<float > cascsetting_v0cosPA{" cascsetting_v0cosPA" , 0 .97f , " Minimum cCosine of the pointing angle of the v0" };
205205 Configurable<float > cascsetting_dcaMaxBachelorToV0{" cascsetting_dcaMaxBachelorToV0" , 0 .8f , " DCA of the bachelor to V0" };
206- // Configurable<float> cascsetting_dcaMinBachelorToProton{"cascsetting_dcaMinBachelorToProton", 0.015f, "DCA of the bachelor to proton"};
206+ // Configurable<float> cascsetting_dcaMinBachelorToProton{"cascsetting_dcaMinBachelorToProton", 0.015f, "DCA of the bachelor to proton"};
207207 Configurable<float > cascsetting_radiusV0{" cascsetting_radiusV0" , 1 .2f , " Minimum radius of the V0 accepted" };
208208 Configurable<float > cascsetting_radiusCasc{" cascsetting_radiusCasc" , 0 .5f , " Minimum radius of the cascade accepted" };
209-
209+
210210 Configurable<float > cascsetting_massWindowOmega{" cascsetting_massWindowOmega" , 0 .01f , " Mass window for the Omega" };
211211 Configurable<float > cascsetting_massWindowXi{" cascsetting_massWindowXi" , 0 .01f , " Mass window for the Xi" };
212212 Configurable<float > cascsetting_nsigmatpc{" cascsetting_nsigmaTPC" , 3 .f , " Number of sigmas for the TPC PID" };
@@ -243,34 +243,35 @@ struct LfTreeCreatorClusterStudies {
243243
244244 HistogramRegistry m_hAnalysis{
245245 " LFTreeCreator" ,
246- {{" collision_selections" , " Collision selection; selection; counts" , {HistType::kTH1F , {{Selections::kAll , -0.5 , static_cast <double >(Selections::kAll ) - 0.5 }}}},
247- {" v0_selections" , " V0 selection; selection; counts" , {HistType::kTH1F , {{V0Selections::kV0All , -0.5 , static_cast <double >(V0Selections::kV0All ) - 0.5 }}}},
248- {" casc_selections" , " Cascade selection; selection; counts" , {HistType::kTH1F , {{CascSelections::kCascAll , -0.5 , static_cast <double >(CascSelections::kCascAll ) - 0.5 }}}},
249- {" e_selections" , " e^{#pm} selection; selection; counts" , {HistType::kTH1F , {{ESelections::kEAll , -0.5 , static_cast <double >(ESelections::kEAll ) - 0.5 }}}},
250- {" v0_type" , " Selected V0; particle; counts" , {HistType::kTH1F , {{V0Type::V0TypeAll, -0.5 , static_cast <double >(V0Type::V0TypeAll) - 0.5 }}}},
251- {" radiusV0" , " Decay radius (xy) V0; radius (cm); counts" , {HistType::kTH1F , {{100 , 0 ., 100 .}}}},
252- {" massLambda" , " #Lambda invariant mass; signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Lambda} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {200 , 1 .08f , 1 .18f }}}},
253- {" massLambdaMc" , " #Lambda invariant mass (MC); signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Lambda} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {200 , 1 .08f , 1 .18f }}}},
254- {" Lambda_vs_K0s" , " Mass #Lambda vs K^{0}_s; #it{m}_{K^{0}_{s}} (GeV/#it{c}^{2}); #it{m}_{#Lambda} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{50 , 0 .f , 1 .f }, {70 , 0 .6f , 2 .f }}}},
255- {" armenteros_plot_before_selections" , " Armenteros-Podolanski plot; #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
256- {" armenteros_plot" , " Armenteros-Podolanski plot; #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
257- {" armenteros_plot_lambda" , " Armenteros-Podolanski plot (#Lambda only); #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
258- {" armenteros_plot_gamma" , " Armenteros-Podolanski plot (#gamma only); #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
259- {" photon_radiusV0" , " Photon conversion radius (xy) V0; #it{r} (cm); counts" , {HistType::kTH1F , {{100 , 0 ., 100 .}}}},
260- {" photon_conversion_position" , " Photon conversion position; #it{x} (cm); #it{y} (cm)" , {HistType::kTH2F , {{250 , -5 .f , 5 .f }, {250 , -5 .f , 5 .f }}}},
261- {" photon_conversion_position_layer" , " Photon conversion position (ITS layers); #it{x} (cm); #it{y} (cm)" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {100 , -5 .f , 5 .f }}}},
262- {" casc_dca_daughter_pairs" , " DCA (xy) for cascade daughter pairs; DCA_{#it{xy}} (cm); counts" , {HistType::kTH1F , {{100 , -0.1 , 0.1 }}}},
263- {" Xi_vs_Omega" , " Mass Xi vs Omega; mass Omega (GeV/#it{c}^{2}); #it{m}_#Xi (GeV/#it{c}^{2})" , {HistType::kTH2F , {{50 , 1 .f , 2 .f }, {50 , 1 .f , 2 .f }}}},
264- {" massOmega" , " Mass #Omega; signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Omega} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {400 , 1 .62f , 1 .72f }}}},
265- {" massOmegaMc" , " Mass #Omega (MC); signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Omega} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {400 , 1 .62f , 1 .72f }}}},
266- {" massPi0" , " Mass #pi^{0}; #it{m}_{#pi^{0}} (GeV/#it{c}^{2})" , {HistType::kTH1F , {{100 , 0 .0f , 0 .200f }}}},
267- {" massPi0Mc" , " Mass #pi^{0} (MC); #it{m}_{#pi^{0}} (GeV/#it{c}^{2})" , {HistType::kTH1F , {{100 , 0 .0f , 0 .200f }}}},
268- {" massPi0WithBkg" , " Mass #pi^{0} with Background; #it{m}_{#pi^{0}} (GeV/#it{c}^{2}); counts" , {HistType::kTH1F , {{100 , 0 .0f , 0 .200f }}}},
269- {" zVtx" , " Binning for the vertex z in cm; #it{z}_{vertex} (cm)" , {HistType::kTH1F , {{100 , -20 .f , 20 .f }}}},
270- {" isPositive" , " is the candidate positive?; isPositive; counts" , {HistType::kTH1F , {{2 , -0 .5f , 1 .5f }}}},
271-
272- {" electron/DCAxyBeforeSelection" , " DCA (xy) for cascade daughter pairs; DCA_{#it{xy}} (cm); counts" , {HistType::kTH1F , {{100 , -0.1 , 0.1 }}}},
273- {" electron/DCAzBeforeSelection" , " DCA (z) for cascade daughter pairs; DCA_{#it{z}} (cm); counts" , {HistType::kTH1F , {{200 , -0.2 , 0.2 }}}},
246+ {
247+ {" collision_selections" , " Collision selection; selection; counts" , {HistType::kTH1F , {{Selections::kAll , -0.5 , static_cast <double >(Selections::kAll ) - 0.5 }}}},
248+ {" v0_selections" , " V0 selection; selection; counts" , {HistType::kTH1F , {{V0Selections::kV0All , -0.5 , static_cast <double >(V0Selections::kV0All ) - 0.5 }}}},
249+ {" casc_selections" , " Cascade selection; selection; counts" , {HistType::kTH1F , {{CascSelections::kCascAll , -0.5 , static_cast <double >(CascSelections::kCascAll ) - 0.5 }}}},
250+ {" e_selections" , " e^{#pm} selection; selection; counts" , {HistType::kTH1F , {{ESelections::kEAll , -0.5 , static_cast <double >(ESelections::kEAll ) - 0.5 }}}},
251+ {" v0_type" , " Selected V0; particle; counts" , {HistType::kTH1F , {{V0Type::V0TypeAll, -0.5 , static_cast <double >(V0Type::V0TypeAll) - 0.5 }}}},
252+ {" radiusV0" , " Decay radius (xy) V0; radius (cm); counts" , {HistType::kTH1F , {{100 , 0 ., 100 .}}}},
253+ {" massLambda" , " #Lambda invariant mass; signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Lambda} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {200 , 1 .08f , 1 .18f }}}},
254+ {" massLambdaMc" , " #Lambda invariant mass (MC); signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Lambda} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {200 , 1 .08f , 1 .18f }}}},
255+ {" Lambda_vs_K0s" , " Mass #Lambda vs K^{0}_s; #it{m}_{K^{0}_{s}} (GeV/#it{c}^{2}); #it{m}_{#Lambda} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{50 , 0 .f , 1 .f }, {70 , 0 .6f , 2 .f }}}},
256+ {" armenteros_plot_before_selections" , " Armenteros-Podolanski plot; #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
257+ {" armenteros_plot" , " Armenteros-Podolanski plot; #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
258+ {" armenteros_plot_lambda" , " Armenteros-Podolanski plot (#Lambda only); #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
259+ {" armenteros_plot_gamma" , " Armenteros-Podolanski plot (#gamma only); #alpha; #it{q}_{T} (GeV/#it{c})" , {HistType::kTH2F , {{100 , -1 .f , 1 .f }, {100 , 0 .f , 0 .3f }}}},
260+ {" photon_radiusV0" , " Photon conversion radius (xy) V0; #it{r} (cm); counts" , {HistType::kTH1F , {{100 , 0 ., 100 .}}}},
261+ {" photon_conversion_position" , " Photon conversion position; #it{x} (cm); #it{y} (cm)" , {HistType::kTH2F , {{250 , -5 .f , 5 .f }, {250 , -5 .f , 5 .f }}}},
262+ {" photon_conversion_position_layer" , " Photon conversion position (ITS layers); #it{x} (cm); #it{y} (cm)" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {100 , -5 .f , 5 .f }}}},
263+ {" casc_dca_daughter_pairs" , " DCA (xy) for cascade daughter pairs; DCA_{#it{xy}} (cm); counts" , {HistType::kTH1F , {{100 , -0.1 , 0.1 }}}},
264+ {" Xi_vs_Omega" , " Mass Xi vs Omega; mass Omega (GeV/#it{c}^{2}); #it{m}_#Xi (GeV/#it{c}^{2})" , {HistType::kTH2F , {{50 , 1 .f , 2 .f }, {50 , 1 .f , 2 .f }}}},
265+ {" massOmega" , " Mass #Omega; signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Omega} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {400 , 1 .62f , 1 .72f }}}},
266+ {" massOmegaMc" , " Mass #Omega (MC); signed #it{p}_{T} (GeV/#it{c}); #it{m}_{#Omega} (GeV/#it{c}^{2})" , {HistType::kTH2F , {{100 , -5 .f , 5 .f }, {400 , 1 .62f , 1 .72f }}}},
267+ {" massPi0" , " Mass #pi^{0}; #it{m}_{#pi^{0}} (GeV/#it{c}^{2})" , {HistType::kTH1F , {{100 , 0 .0f , 0 .200f }}}},
268+ {" massPi0Mc" , " Mass #pi^{0} (MC); #it{m}_{#pi^{0}} (GeV/#it{c}^{2})" , {HistType::kTH1F , {{100 , 0 .0f , 0 .200f }}}},
269+ {" massPi0WithBkg" , " Mass #pi^{0} with Background; #it{m}_{#pi^{0}} (GeV/#it{c}^{2}); counts" , {HistType::kTH1F , {{100 , 0 .0f , 0 .200f }}}},
270+ {" zVtx" , " Binning for the vertex z in cm; #it{z}_{vertex} (cm)" , {HistType::kTH1F , {{100 , -20 .f , 20 .f }}}},
271+ {" isPositive" , " is the candidate positive?; isPositive; counts" , {HistType::kTH1F , {{2 , -0 .5f , 1 .5f }}}},
272+
273+ {" electron/DCAxyBeforeSelection" , " DCA (xy) for cascade daughter pairs; DCA_{#it{xy}} (cm); counts" , {HistType::kTH1F , {{100 , -0.1 , 0.1 }}}},
274+ {" electron/DCAzBeforeSelection" , " DCA (z) for cascade daughter pairs; DCA_{#it{z}} (cm); counts" , {HistType::kTH1F , {{200 , -0.2 , 0.2 }}}},
274275 },
275276 OutputObjHandlingPolicy::AnalysisObject,
276277 false };
@@ -381,7 +382,7 @@ struct LfTreeCreatorClusterStudies {
381382 }
382383
383384 template <typename Track>
384- bool selectPidV0Daughters (Candidate& candidatePos, Candidate& candidateNeg,
385+ bool selectPidV0Daughters (Candidate& candidatePos, Candidate& candidateNeg,
385386 aod::V0Datas::iterator const & v0, const Track& posTrack, const Track& negTrack, uint8_t v0Bitmask)
386387 {
387388 if (TESTBIT (v0Bitmask, Lambda)) {
@@ -785,7 +786,7 @@ struct LfTreeCreatorClusterStudies {
785786 m_hAnalysis.fill (HIST (" v0_selections" ), V0Selections::kV0Topology );
786787
787788 std::array<float , 3 > momPos{v0.pxpos (), v0.pypos (), v0.pzpos ()},
788- momNeg{v0.pxneg (), v0.pyneg (), v0.pzneg ()};
789+ momNeg{v0.pxneg (), v0.pyneg (), v0.pzneg ()};
789790
790791 m_hAnalysis.fill (HIST (" armenteros_plot_before_selections" ), v0.alpha (), v0.qtarm ());
791792 m_hAnalysis.fill (HIST (" Lambda_vs_K0s" ), v0.mK0Short (), v0.mAntiLambda ());
@@ -850,7 +851,7 @@ struct LfTreeCreatorClusterStudies {
850851 const float & massXi = cascade.mXi ();
851852 const float & massOmega = cascade.mOmega ();
852853 m_hAnalysis.fill (HIST (" Xi_vs_Omega" ), massOmega, massXi);
853-
854+
854855 if (std::abs (massXi - o2::constants::physics::MassXiMinus) < cascsetting_massWindowXi) {
855856 return ;
856857 }
@@ -1094,7 +1095,7 @@ struct LfTreeCreatorClusterStudies {
10941095
10951096 // =========================================================================================================
10961097
1097- void processDataV0Casc (CollisionsCustom::iterator const & collision/* s*/ , TracksFullIU const & tracks, aod::V0Datas const & v0s, aod::CascDatas const & cascades, aod::BCsWithTimestamps const &)
1098+ void processDataV0Casc (CollisionsCustom::iterator const & collision /* s*/ , TracksFullIU const & tracks, aod::V0Datas const & v0s, aod::CascDatas const & cascades, aod::BCsWithTimestamps const &)
10981099 {
10991100 if (!collisionSelection (collision)) {
11001101 return ;
@@ -1184,7 +1185,7 @@ struct LfTreeCreatorClusterStudies {
11841185 }
11851186 PROCESS_SWITCH (LfTreeCreatorClusterStudies, processDataPKPi, " process Data p, K, pi" , false );
11861187
1187- void processMcV0Casc (CollisionsCustomMc::iterator const & collision/* s*/ , TracksFullIUMc const & tracks, aod::V0Datas const & v0s, aod::CascDatas const & cascades, aod::BCsWithTimestamps const &, aod::McParticles const &, aod::McCollisions const &)
1188+ void processMcV0Casc (CollisionsCustomMc::iterator const & collision /* s*/ , TracksFullIUMc const & tracks, aod::V0Datas const & v0s, aod::CascDatas const & cascades, aod::BCsWithTimestamps const &, aod::McParticles const &, aod::McCollisions const &)
11881189 {
11891190 if (!collisionSelection (collision)) {
11901191 return ;
0 commit comments