@@ -47,14 +47,20 @@ struct f0980analysis {
4747
4848 void init (o2::framework::InitContext&)
4949 {
50+ std::vector<double > ptBinning = {
51+ 0.0 , 0.2 , 0.4 , 0.6 , 0.8 ,
52+ 1.0 , 1.5 , 2.0 , 2.5 , 3.0 ,
53+ 3.5 , 4.0 , 4.5 , 5.0 , 6.0 ,
54+ 7.0 , 8.0 , 10.0 , 13.0 , 20.0 };
55+
5056 AxisSpec centAxis = {20 , 0 , 100 };
51- AxisSpec ptAxis = {100 , 0 , 20 };
57+ AxisSpec ptAxis = {ptBinning };
5258 AxisSpec massAxis = {400 , 0.2 , 2.2 };
5359 AxisSpec epAxis = {20 , -constants::math::PI , constants::math::PI };
5460
55- histos.add (" hInvMass_f0980_US" , " unlike invariant mass" , {HistType::kTHnF , {massAxis, ptAxis, centAxis, epAxis }});
56- histos.add (" hInvMass_f0980_LSpp" , " ++ invariant mass" , {HistType::kTHnF , {massAxis, ptAxis, centAxis, epAxis }});
57- histos.add (" hInvMass_f0980_LSmm" , " -- invariant mass" , {HistType::kTHnF , {massAxis, ptAxis, centAxis, epAxis }});
61+ histos.add (" hInvMass_f0980_US" , " unlike invariant mass" , {HistType::kTH3F , {massAxis, ptAxis, centAxis}});
62+ histos.add (" hInvMass_f0980_LSpp" , " ++ invariant mass" , {HistType::kTH3F , {massAxis, ptAxis, centAxis}});
63+ histos.add (" hInvMass_f0980_LSmm" , " -- invariant mass" , {HistType::kTH3F , {massAxis, ptAxis, centAxis}});
5864
5965 histos.print ();
6066 }
@@ -66,7 +72,7 @@ struct f0980analysis {
6672 {
6773 if (track.pt () < cfgMinPt)
6874 return false ;
69- if (track.pt () < cfgMaxPt)
75+ if (track.pt () > cfgMaxPt)
7076 return false ;
7177 if (std::fabs (track.eta ()) > cfgMaxEta)
7278 return false ;
@@ -90,7 +96,6 @@ struct f0980analysis {
9096 return false ;
9197 }
9298 }
93-
9499 return true ;
95100 }
96101
@@ -111,11 +116,11 @@ struct f0980analysis {
111116 if (Reco.Rapidity () > cfgMaxRap || Reco.Rapidity () < cfgMinRap)
112117 continue ;
113118 if (trk1.sign () * trk2.sign () < 0 ) {
114- histos.fill (HIST (" hInvMass_f0980_US" ), Reco.M (), Reco.Pt (), collision.multV0M (), 0 );
119+ histos.fill (HIST (" hInvMass_f0980_US" ), Reco.M (), Reco.Pt (), collision.multV0M ());
115120 } else if (trk1.sign () > 0 && trk2.sign () > 0 ) {
116- histos.fill (HIST (" hInvMass_f0980_LSpp" ), Reco.M (), Reco.Pt (), collision.multV0M (), 0 );
121+ histos.fill (HIST (" hInvMass_f0980_LSpp" ), Reco.M (), Reco.Pt (), collision.multV0M ());
117122 } else if (trk1.sign () < 0 && trk2.sign () < 0 ) {
118- histos.fill (HIST (" hInvMass_f0980_LSmm" ), Reco.M (), Reco.Pt (), collision.multV0M (), 0 );
123+ histos.fill (HIST (" hInvMass_f0980_LSmm" ), Reco.M (), Reco.Pt (), collision.multV0M ());
119124 }
120125 }
121126 }
@@ -125,8 +130,8 @@ struct f0980analysis {
125130 {
126131 LOGF (debug, " [DATA] Processing %d collisions" , collisions.size ());
127132 for (auto & collision : collisions) {
133+ Partition<aod::ResoTracks> selectedTracks = o2::aod::track::pt > static_cast <float_t >(cfgMinPt) && (nabs (o2::aod::track::dcaZ) > static_cast <float_t >(cfgMinDCAzToPVcut)) && (nabs (o2::aod::track::dcaZ) > static_cast <float_t >(cfgMaxDCAzToPVcut)) && (nabs (o2::aod::track::dcaXY) < static_cast <float_t >(cfgMaxDCArToPVcut));
128134
129- Partition<aod::ResoTracks> selectedTracks = requireTOFPIDPionCutInFilter ();
130135 selectedTracks.bindTable (resotracks);
131136 auto colTracks = selectedTracks->sliceByCached (aod::resodaughter::resoCollisionId, collision.globalIndex (), cache);
132137 fillHistograms<false >(collision, colTracks);
0 commit comments