4343#include < TProfile3D.h>
4444#include < TString.h>
4545
46- #include < memory>
4746#include < array>
4847#include < chrono>
4948#include < cmath>
5049#include < cstddef>
5150#include < cstdint>
51+ #include < memory>
5252#include < numeric>
5353#include < string>
5454#include < string_view>
@@ -157,7 +157,6 @@ struct ZdcQVectors {
157157 Configurable<std::vector<double >> cfgEvSelsMultPv{" cfgEvSelsMultPv" , std::vector<double >{2223.49 , -75.1444 , 0.963572 , -0.00570399 , 1.34877e-05 , 3790.99 , -137.064 , 2.13044 , -0.017122 , 5.82834e-05 }, " Multiplicity cuts (PV) first 5 parameters cutLOW last 5 cutHIGH (Default is +-2sigma pass5) " };
158158 Configurable<std::vector<double >> cfgEvSelsMult{" cfgEvSelsMult" , std::vector<double >{1301.56 , -41.4615 , 0.478224 , -0.00239449 , 4.46966e-06 , 2967.6 , -102.927 , 1.47488 , -0.0106534 , 3.28622e-05 }, " Multiplicity cuts (Global) first 5 parameters cutLOW last 5 cutHIGH (Default is +-2sigma pass5) " };
159159
160-
161160 // define my.....
162161 // Filter collisionFilter = nabs(aod::collision::posZ) <;
163162
@@ -167,7 +166,6 @@ struct ZdcQVectors {
167166 using UnfilteredTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TrackSelection, aod::TracksDCA>;
168167 using UsedTracks = soa::Filtered<UnfilteredTracks>;
169168
170-
171169 enum SelectionCriteria {
172170 evSel_FilteredEvent,
173171 evSel_BCHasZDC,
@@ -408,7 +406,6 @@ struct ZdcQVectors {
408406 fMultCutLow ->SetParameters (paramsMultCut[0 ], paramsMultCut[1 ], paramsMultCut[2 ], paramsMultCut[3 ], paramsMultCut[4 ]);
409407 fMultCutHigh ->SetParameters (paramsMultCut[5 ], paramsMultCut[6 ], paramsMultCut[7 ], paramsMultCut[8 ], paramsMultCut[9 ]);
410408 }
411-
412409 }
413410
414411 double rescaleTimestamp (uint64_t timestamp, int runnumber)
@@ -471,7 +468,7 @@ struct ZdcQVectors {
471468 }
472469 }
473470
474- template <typename TCollision, typename TBunchCrossing>
471+ template <typename TCollision, typename TBunchCrossing>
475472 uint16_t eventSelected (TCollision collision, TBunchCrossing bunchCrossing, bool & isEventSelected, const int & multTrk)
476473 {
477474 uint16_t selectionBits = 0 ;
@@ -943,13 +940,13 @@ struct ZdcQVectors {
943940 registry.get <TProfile>(HIST (" vmean/hvertex_vz" ))->Fill (Form (" %d" , runnumber), v[2 ]);
944941
945942 // Fill to get mean energy per tower in 1% centrality bins
946-
943+
947944 registry.get <TProfile2D>(HIST (" Energy/hZNA_mean_t0_cent" ))->Fill (Form (" %d" , runnumber), cent, zdcCol.energyCommonZNA (), 1 );
948945 registry.get <TProfile2D>(HIST (" Energy/hZNA_mean_t1_cent" ))->Fill (Form (" %d" , runnumber), cent, eZN[0 ], 1 );
949946 registry.get <TProfile2D>(HIST (" Energy/hZNA_mean_t2_cent" ))->Fill (Form (" %d" , runnumber), cent, eZN[1 ], 1 );
950947 registry.get <TProfile2D>(HIST (" Energy/hZNA_mean_t3_cent" ))->Fill (Form (" %d" , runnumber), cent, eZN[2 ], 1 );
951948 registry.get <TProfile2D>(HIST (" Energy/hZNA_mean_t4_cent" ))->Fill (Form (" %d" , runnumber), cent, eZN[3 ], 1 );
952-
949+
953950 registry.get <TProfile2D>(HIST (" Energy/hZNC_mean_t0_cent" ))->Fill (Form (" %d" , runnumber), cent, zdcCol.energyCommonZNC (), 1 );
954951 registry.get <TProfile2D>(HIST (" Energy/hZNC_mean_t1_cent" ))->Fill (Form (" %d" , runnumber), cent, eZN[4 ], 1 );
955952 registry.get <TProfile2D>(HIST (" Energy/hZNC_mean_t2_cent" ))->Fill (Form (" %d" , runnumber), cent, eZN[5 ], 1 );
0 commit comments