2121#include " Common/DataModel/EventSelection.h"
2222#include " Common/DataModel/PIDResponseTOF.h"
2323#include " Common/DataModel/TrackSelectionTables.h"
24- # include " MathUtils/BetheBlochAleph.h "
24+
2525#include " Framework/AnalysisDataModel.h"
2626#include " Framework/AnalysisTask.h"
2727#include " Framework/HistogramRegistry.h"
2828#include " Framework/runDataProcessing.h"
29+ #include " MathUtils/BetheBlochAleph.h"
30+
2931#include < TMCProcess.h>
3032#include < TParameter.h>
33+
3134#include < cmath>
3235#include < string>
3336#include < vector>
@@ -52,7 +55,7 @@ static const int kMinTpcCrossedRows = 70;
5255static const double kBetheBlochDefault [6 ]{-1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 };
5356static const std::vector<std::string> kParamNamesBB {" p0" , " p1" , " p2" , " p3" , " p4" , " resolution" };
5457static const std::vector<std::string> kParticleNamesBB {" He3" };
55- }
58+ } // namespace
5659
5760struct AntiNucleiTask {
5861 HistogramRegistry histo{" histos" , {}, OutputObjHandlingPolicy::AnalysisObject};
@@ -98,13 +101,20 @@ struct AntiNucleiTask {
98101 template <typename T>
99102 bool passTrackSelection (const T& track)
100103 {
101- if (std::abs (track.eta ()) > kMaxEta ) return false ;
102- if (track.tpcNClsFound () < mMinClsTPC ) return false ;
103- if (track.tpcNClsCrossedRows () < kMinTpcCrossedRows ) return false ;
104- if (track.itsNCls () < mMinClsITS ) return false ;
105- if (track.tpcChi2NCl () > mMaxChi2TPC || track.tpcChi2NCl () < mMinChi2TPC ) return false ;
106- if (track.itsChi2NCl () > mMaxChi2ITS ) return false ;
107- if (std::abs (track.dcaXY ()) > mMaxDCAXY || std::abs (track.dcaZ ()) > mMaxDCAZ ) return false ;
104+ if (std::abs (track.eta ()) > kMaxEta )
105+ return false ;
106+ if (track.tpcNClsFound () < mMinClsTPC )
107+ return false ;
108+ if (track.tpcNClsCrossedRows () < kMinTpcCrossedRows )
109+ return false ;
110+ if (track.itsNCls () < mMinClsITS )
111+ return false ;
112+ if (track.tpcChi2NCl () > mMaxChi2TPC || track.tpcChi2NCl () < mMinChi2TPC )
113+ return false ;
114+ if (track.itsChi2NCl () > mMaxChi2ITS )
115+ return false ;
116+ if (std::abs (track.dcaXY ()) > mMaxDCAXY || std::abs (track.dcaZ ()) > mMaxDCAZ )
117+ return false ;
108118 return true ;
109119 }
110120
@@ -120,10 +130,12 @@ struct AntiNucleiTask {
120130
121131 void process (CollisionWithEvSel::iterator const & collision, TotalTracks const & tracks)
122132 {
123- if (std::abs (collision.posZ ()) > kMaxVertexZ ) return ;
133+ if (std::abs (collision.posZ ()) > kMaxVertexZ )
134+ return ;
124135 histo.fill (HIST (" ZVtx_Raw" ), collision.posZ ());
125136
126- if (!collision.sel8 ()) return ;
137+ if (!collision.sel8 ())
138+ return ;
127139 histo.fill (HIST (" ZVtx" ), collision.posZ ());
128140
129141 for (const auto & track : tracks) {
@@ -159,7 +171,8 @@ struct AntiNucleiTask {
159171 void processMC (CollisionMC const & collisions, aod::McCollisions const &, TracksMC const & tracks, aod::McParticles const & particlesMC)
160172 {
161173 for (const auto & collision : collisions) {
162- if (std::abs (collision.posZ ()) > kMaxVertexZ || !collision.sel8 ()) continue ;
174+ if (std::abs (collision.posZ ()) > kMaxVertexZ || !collision.sel8 ())
175+ continue ;
163176
164177 for (const auto & particle : particlesMC) {
165178 if (particle.mcCollisionId () == collision.mcCollisionId () && isPrimaryHe3 (particle)) {
@@ -170,7 +183,8 @@ struct AntiNucleiTask {
170183 }
171184
172185 for (const auto & track : tracks) {
173- if (track.collisionId () != collision.index () || !passTrackSelection (track) || track.mcParticleId () == -1 ) continue ;
186+ if (track.collisionId () != collision.index () || !passTrackSelection (track) || track.mcParticleId () == -1 )
187+ continue ;
174188
175189 auto particle = particlesMC.iteratorAt (track.mcParticleId ());
176190 if (isPrimaryHe3 (particle)) {
0 commit comments