@@ -57,7 +57,6 @@ struct MultandPtFluctuations {
5757 Configurable<bool > cFT0C{" cFT0C" , true , " Use FT0C centrality" };
5858
5959 // ------ Track cuts
60-
6160 Configurable<float > ptMinCut{" ptMinCut" , 0 .2f , " Minimum pT" };
6261 Configurable<float > ptMaxCut{" ptMaxCut" , 5 .0f , " Maximum pT" };
6362 Configurable<float > etaCut{" etaCut" , 0 .8f , " Maximum |eta|" };
@@ -73,11 +72,7 @@ struct MultandPtFluctuations {
7372 Configurable<bool > requireInnerITS{" requireInnerITS" , true , " At least one hit in ITS layers 0,1,2" };
7473
7574 // ----- Histogram Registry
76- HistogramRegistry histos{
77- " histos" ,
78- {},
79- OutputObjHandlingPolicy::AnalysisObject};
80-
75+ HistogramRegistry histos{" histos" , {}, OutputObjHandlingPolicy::AnalysisObject};
8176 void init (InitContext const &)
8277 {
8378 const AxisSpec axisEta{30 , -1.5 , +1.5 , " #eta" };
@@ -88,7 +83,6 @@ struct MultandPtFluctuations {
8883 const AxisSpec axisdcaZ{ndcaZ, -3.0 , 3.0 , " DCAz" };
8984 const AxisSpec axisNch{500 , 0 , 500 , " Nch" };
9085 const AxisSpec axisCent{nCentBins, 0.0 , 100.0 , " Centrality (%)" }; // ------------- centrality bins
91-
9286 histos.add (" QA/BeforeCut/VtxZ" , " Vertex Z" , kTH1F , {axisVtxz});
9387 histos.add (" QA/AfterCut/VtxZ" , " Vertex Z" , kTH1F , {axisVtxz});
9488 histos.add (" QA/BeforeCut/Cent" , " Centrality" , kTH1F , {axisCent});
@@ -127,41 +121,31 @@ struct MultandPtFluctuations {
127121 void process (myColsData::iterator const & col, myFilteredTracksData const & tracks)
128122 {
129123 histos.fill (HIST (" QA/BeforeCut/VtxZ" ), col.posZ ());
130-
131124 float cent = -1 .0f ;
132-
133125 if (cFT0M) {
134126 cent = col.centFT0M ();
135127 }
136-
137128 if (cFT0C) {
138129 cent = col.centFT0C ();
139130 // if (!track.isGlobalTrack()) continue;
140131 }
141-
142132 if (cent < 0 ) {
143133 return ;
144134 }
145135
146136 histos.fill (HIST (" QA/BeforeCut/Cent" ), cent);
147-
148137 if (!col.sel8 ())
149138 return ;
150-
151139 if (std::abs (col.posZ ()) > vtxZcut)
152140 return ;
153-
154141 if (cfgNoSameBunchPileup && !col.selection_bit (o2::aod::evsel::kNoSameBunchPileup ))
155142 return ;
156-
157143 if (cfgEvSelUseGoodZvtxFT0vsPV && !col.selection_bit (o2::aod::evsel::kIsGoodZvtxFT0vsPV ))
158144 return ;
159-
160145 if (cfgUseGoodITSLayerAllCut && !col.selection_bit (o2::aod::evsel::kIsGoodITSLayersAll ))
161146 return ;
162147
163148 histos.fill (HIST (" hEventCounter" ), cent);
164-
165149 histos.fill (HIST (" QA/AfterCut/VtxZ" ), col.posZ ());
166150 histos.fill (HIST (" QA/AfterCut/Cent" ), cent);
167151
@@ -183,32 +167,24 @@ struct MultandPtFluctuations {
183167 continue ;
184168 if (requireTPC && !track.hasTPC ())
185169 continue ;
186-
187170 if (track.itsNCls () < itsNClsCut)
188171 continue ;
189172 if (track.tpcNClsCrossedRows () < tpcCrossCut)
190173 continue ;
191-
192174 if (track.tpcCrossedRowsOverFindableCls () < crossedRowsOverFindableCut)
193175 continue ;
194-
195176 if (std::abs (track.dcaZ ()) > dcaZCut)
196177 continue ;
197178 if (std::abs (track.dcaXY ()) > dcaXYCut)
198179 continue ;
199-
200180 if (track.tpcChi2NCl () > tpcChiCut)
201181 continue ;
202182 if (track.itsChi2NCl () > itsChiCut)
203183 continue ;
204-
205184 if (std::abs (track.eta ()) > etaCut)
206185 continue ;
207-
208186 if (requireInnerITS) {
209-
210187 auto itsMap = track.itsClusterMap ();
211-
212188 if (!(itsMap & (1 << 0 )) &&
213189 !(itsMap & (1 << 1 )) &&
214190 !(itsMap & (1 << 2 ))) {
@@ -222,10 +198,8 @@ struct MultandPtFluctuations {
222198 histos.fill (HIST (" QA/AfterCut/Phi" ), track.phi ());
223199 histos.fill (HIST (" QA/AfterCut/DcaXY" ), track.dcaXY ());
224200 histos.fill (HIST (" QA/AfterCut/DcaZ" ), track.dcaZ ());
225-
226201 histos.fill (HIST (" h2_DcaZ" ), track.pt (), track.dcaZ ());
227202 histos.fill (HIST (" h2_DcaXY" ), track.pt (), track.dcaXY ());
228-
229203 a++; // A-nch, B =pt
230204 b += track.pt ();
231205 }
0 commit comments