@@ -49,15 +49,13 @@ struct MultandptFluctuations {
4949
5050
5151 // ------ Event Cuts
52-
5352 Configurable<float > vtxZcut{" vertexZcut" , 10 .0f , " Vertex Z" };
5453 Configurable<bool > cfgNoSameBunchPileup{" cfgNoSameBunchPileup" , true , " kNoSameBunchPileup" };
5554 Configurable<bool > cfgEvSelUseGoodZvtxFT0vsPV{" cfgEvSelUseGoodZvtxFT0vsPV" ,true , " Good Zvtx FT0 vs PV" };
5655 Configurable<bool > cfgUseGoodITSLayerAllCut{" cfgUseGoodITSLayerAllCut" , true , " Good ITS Layers" };
5756
58-
59- // ----- Centrality estimator
6057
58+ // ----- Centrality estimator
6159 Configurable<bool > cFT0M{" cFT0M" , false , " FT0M centrality" }; // for pp
6260 Configurable<bool > cFT0C{" cFT0C" , true , " Use FT0C centrality" };
6361
@@ -77,16 +75,13 @@ struct MultandptFluctuations {
7775 Configurable<bool > requireITS{" requireITS" ,true , " Require ITS hit" };
7876 Configurable<bool > requireTPC{" requireTPC" , true , " Require TPC hit" };
7977 Configurable<bool > requireInnerITS{" requireInnerITS" , true ," At least one hit in ITS layers 0,1,2" };
80-
8178
8279
8380 // ----- Histogram Registry
84-
8581 HistogramRegistry histos{
8682 " histos" , {}, OutputObjHandlingPolicy::AnalysisObject
8783 };
8884
89-
9085 void init (InitContext const &)
9186 {
9287 const AxisSpec axisEta{30 , -1.5 , +1.5 , " #eta" };
@@ -98,34 +93,25 @@ struct MultandptFluctuations {
9893 const AxisSpec axisNch{500 , 0 , 500 , " Nch" };
9994 const AxisSpec axisCent{ nCentBins, 0.0 , 100.0 , " Centrality (%)" }; // ------------- centrality bins
10095
101-
96+
10297 histos.add (" QA/BeforeCut/VtxZ" , " Vertex Z" , kTH1F , {axisVtxz});
10398 histos.add (" QA/AfterCut/VtxZ" , " Vertex Z" , kTH1F , {axisVtxz});
104-
10599 histos.add (" QA/BeforeCut/Cent" , " Centrality" , kTH1F , {axisCent});
106100 histos.add (" QA/AfterCut/Cent" , " Centrality" , kTH1F , {axisCent});
107-
108101 histos.add (" QA/BeforeCut/Eta" , " Eta" , kTH1F , {axisEta});
109102 histos.add (" QA/AfterCut/Eta" , " Eta" , kTH1F , {axisEta});
110-
111103 histos.add (" QA/BeforeCut/Pt" , " Pt" , kTH1F , {axisPt});
112104 histos.add (" QA/AfterCut/Pt" , " Pt" , kTH1F , {axisPt});
113-
114105 histos.add (" QA/BeforeCut/Phi" , " Phi" , kTH1F , {axisPhi});
115106 histos.add (" QA/AfterCut/Phi" , " Phi" , kTH1F , {axisPhi});
116-
117107 histos.add (" QA/BeforeCut/DcaXY" , " DCAxy" , kTH1F , {axisdcaXY});
118108 histos.add (" QA/AfterCut/DcaXY" , " DCAxy" , kTH1F , {axisdcaXY});
119-
120109 histos.add (" QA/BeforeCut/DcaZ" , " DCAz" , kTH1F , {axisdcaZ});
121- histos.add (" QA/AfterCut/DcaZ" , " DCAz" , kTH1F , {axisdcaZ});
122-
110+ histos.add (" QA/AfterCut/DcaZ" , " DCAz" , kTH1F , {axisdcaZ});
123111 histos.add (" hNch" , " Nch" , kTH1F , {axisNch});
124112 histos.add (" h2_DcaZ" , " DCA_{Z}" , kTH2D , {{axisPt}, {axisdcaZ}});
125113 histos.add (" h2_DcaXY" , " DCA_{XY}" , kTH2D , {{axisPt}, {axisdcaXY}});
126-
127114 histos.add (" hEventCounter" , " Number of events vs centrality" , kTH1F , {axisCent});
128-
129115 histos.add (" p_a" , " <A> vs centrality" , kTProfile , {axisCent});
130116 histos.add (" p_b" , " <B> vs centrality" , kTProfile , {axisCent});
131117 histos.add (" p_a2" , " <A^{2}> vs centrality" , kTProfile , {axisCent});
@@ -135,12 +121,10 @@ struct MultandptFluctuations {
135121
136122 }
137123
138- // ---------- Filters
139-
124+ // ---------- Filters
140125 Filter ptFilter =(aod::track::pt > ptMinCut) && (aod::track::pt < ptMaxCut);
141126 Filter etaFilter =nabs(aod::track::eta) < etaCut;
142127 Filter posZFilter = nabs(aod::collision::posZ) < vtxZcut;
143-
144128
145129 using myColsData = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms,aod::CentFT0Cs>;
146130 using myTracksData = soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA>;
@@ -152,8 +136,7 @@ struct MultandptFluctuations {
152136
153137 float cent = -1 .0f ;
154138
155- if (cFT0M) {
156- cent = col.centFT0M ();
139+ if (cFT0M) { cent = col.centFT0M ();
157140 }
158141
159142 if (cFT0C) {
@@ -166,13 +149,13 @@ struct MultandptFluctuations {
166149 }
167150
168151 histos.fill (HIST (" QA/BeforeCut/Cent" ), cent);
169-
152+
170153 if (!col.sel8 ()) return ;
171-
154+
172155 if (std::abs (col.posZ ()) > vtxZcut) return ;
173156
174157 if (cfgNoSameBunchPileup && !col.selection_bit (o2::aod::evsel::kNoSameBunchPileup )) return ;
175-
158+
176159 if (cfgEvSelUseGoodZvtxFT0vsPV && !col.selection_bit (o2::aod::evsel::kIsGoodZvtxFT0vsPV )) return ;
177160
178161 if (cfgUseGoodITSLayerAllCut && !col.selection_bit (o2::aod::evsel::kIsGoodITSLayersAll )) return ;
@@ -181,10 +164,10 @@ struct MultandptFluctuations {
181164
182165 histos.fill (HIST (" QA/AfterCut/VtxZ" ), col.posZ ());
183166 histos.fill (HIST (" QA/AfterCut/Cent" ), cent);
184-
167+
185168 double A = 0.0 ;
186- double B = 0.0 ;
187-
169+ double B = 0.0 ;
170+
188171 // if (!track.isGlobalTrack()) continue;
189172 for (auto & track : tracks) {
190173
@@ -195,7 +178,7 @@ struct MultandptFluctuations {
195178 histos.fill (HIST (" QA/BeforeCut/DcaZ" ), track.dcaZ ());
196179
197180 // Track quality cuts
198-
181+
199182 if (requireITS && !track.hasITS ()) continue ;
200183 if (requireTPC && !track.hasTPC ()) continue ;
201184
@@ -222,7 +205,7 @@ struct MultandptFluctuations {
222205 continue ;
223206 }
224207 }
225-
208+
226209
227210 /* if (eta >= 0.6 && eta < 0.8)
228211 nf++;
@@ -239,10 +222,10 @@ struct MultandptFluctuations {
239222
240223 histos.fill (HIST (" h2_DcaZ" ), track.pt (),track.dcaZ ());
241224 histos.fill (HIST (" h2_DcaXY" ), track.pt (),track.dcaXY ());
242-
225+
243226 A++; // A-nch, B =pt
244227 B += track.pt ();
245-
228+
246229 }
247230 // LOG(info) << "Nch = " << nch;
248231 /* histos.fill(HIST("hNch"),nch);
0 commit comments