1313// / \brief Task reading AOD/ZDCEXTRA table
1414// / \author Uliana Dmitrieva <uliana.dmitrieva@cern.ch>, INFN Torino
1515
16- #include " Common/CCDB/EventSelectionParams.h"
1716#include " Common/Core/RecoDecay.h"
18- #include " Common/DataModel/Centrality.h"
19- #include " Common/DataModel/EventSelection.h"
2017#include " Common/DataModel/ZDCExtra.h"
2118
2219#include < CCDB/BasicCCDBManager.h>
23- #include < CCDB/CcdbApi.h>
2420#include < CommonConstants/MathConstants.h>
2521#include < Framework/AnalysisDataModel.h>
2622#include < Framework/AnalysisHelpers.h>
2925#include < Framework/HistogramRegistry.h>
3026#include < Framework/HistogramSpec.h>
3127#include < Framework/InitContext.h>
32- #include < Framework/OutputObjHeader.h>
3328#include < Framework/runDataProcessing.h>
3429
3530#include < TH1.h>
3934#include < TProfile3D.h>
4035
4136#include < Rtypes.h>
37+ #include < TAxis.h>
38+ #include < TString.h>
4239
40+ #include < array>
41+ #include < cmath>
42+ #include < cstddef>
43+ #include < cstdint>
4344#include < memory>
4445#include < string>
4546#include < unordered_map>
@@ -185,8 +186,8 @@ double getMeanQFromMap(THn* h, double cent, double vx, double vy, double vz)
185186 int binVy = axVy->FindFixBin (vy);
186187 int binVz = axVz->FindFixBin (vz);
187188
188- int idx[ 4 ] = {binCent, binVx, binVy, binVz};
189- return h->GetBinContent (idx);
189+ std::array< int , 4 > idx = {binCent, binVx, binVy, binVz};
190+ return h->GetBinContent (idx. data () );
190191}
191192
192193// Helper for 1D recentering maps: returns mean Q for coordinate x
@@ -258,7 +259,7 @@ struct ZdcExtraTableReader {
258259 Configurable<std::string> qRecenteringCcdb{" qRecenteringCcdb" , " Users/u/udmitrie/ZDC/LHC24ar_apass2" , " Recentering maps containing step folder" };
259260
260261 // CCDB
261- Service<o2::ccdb::BasicCCDBManager> ccdb;
262+ Service<o2::ccdb::BasicCCDBManager> ccdb{} ;
262263
263264 // Struct to hold calibration data for a single step
264265 struct CalibStepData {
@@ -351,8 +352,9 @@ struct ZdcExtraTableReader {
351352 template <typename T>
352353 T* safeClone (TObject* obj)
353354 {
354- if (!obj)
355+ if (!obj) {
355356 return nullptr ;
357+ }
356358 T* cloned = dynamic_cast <T*>(obj->Clone ());
357359 if (cloned) {
358360
@@ -740,17 +742,19 @@ struct ZdcExtraTableReader {
740742 constexpr float QvectorMaxValue = 990.0 ;
741743
742744 if (isZNAhit) {
743- int activeTowersZNA = (zdc.znaTow1 () > 0 .) + (zdc.znaTow2 () > 0 .) + (zdc.znaTow3 () > 0 .) + (zdc.znaTow4 () > 0 .);
745+ int activeTowersZNA = static_cast < int > (zdc.znaTow1 () > 0 .) + static_cast < int > (zdc.znaTow2 () > 0 .) + static_cast < int > (zdc.znaTow3 () > 0 .) + static_cast < int > (zdc.znaTow4 () > 0 .);
744746 float znaSum = zdc.znaTow1 () + zdc.znaTow2 () + zdc.znaTow3 () + zdc.znaTow4 ();
745- if (activeTowersZNA >= minNTowersFired && znaSum > 0 && zdc.znaQx () < QvectorMaxValue)
747+ if (activeTowersZNA >= minNTowersFired && znaSum > 0 && zdc.znaQx () < QvectorMaxValue) {
746748 isZNASpDeterminable = true ;
749+ }
747750 }
748751
749752 if (isZNChit) {
750- int activeTowersZNC = (zdc.zncTow1 () > 0 .) + (zdc.zncTow2 () > 0 .) + (zdc.zncTow3 () > 0 .) + (zdc.zncTow4 () > 0 .);
753+ int activeTowersZNC = static_cast < int > (zdc.zncTow1 () > 0 .) + static_cast < int > (zdc.zncTow2 () > 0 .) + static_cast < int > (zdc.zncTow3 () > 0 .) + static_cast < int > (zdc.zncTow4 () > 0 .);
751754 float zncSum = zdc.zncTow1 () + zdc.zncTow2 () + zdc.zncTow3 () + zdc.zncTow4 ();
752- if (activeTowersZNC >= minNTowersFired && zncSum > 0 && zdc.zncQx () < QvectorMaxValue)
755+ if (activeTowersZNC >= minNTowersFired && zncSum > 0 && zdc.zncQx () < QvectorMaxValue) {
753756 isZNCSpDeterminable = true ;
757+ }
754758 }
755759
756760 if (plotPMs) {
@@ -805,8 +809,9 @@ struct ZdcExtraTableReader {
805809
806810 int cacheIdx = step - 1 ;
807811 // Check if index is valid within cached vector
808- if (cacheIdx >= static_cast <int >(calibCache.size ()))
812+ if (cacheIdx >= static_cast <int >(calibCache.size ())) {
809813 continue ;
814+ }
810815
811816 const auto & calib = calibCache[cacheIdx];
812817
@@ -832,8 +837,8 @@ struct ZdcExtraTableReader {
832837 }
833838 }
834839
835- double valuesQxZNA[ 5 ] = {cent, vx, vy, vz, qxZNArec};
836- double valuesQyZNA[ 5 ] = {cent, vx, vy, vz, qyZNArec};
840+ std::array< double , 5 > valuesQxZNA = {cent, vx, vy, vz, qxZNArec};
841+ std::array< double , 5 > valuesQyZNA = {cent, vx, vy, vz, qyZNArec};
837842
838843 gCurrentCentroidZNA ->Fill (qxZNArec, qyZNArec);
839844
@@ -851,8 +856,8 @@ struct ZdcExtraTableReader {
851856 gCurrentQyVsTimeZNA ->Fill (timeInMinutes, qyZNArec);
852857
853858 if (plot5D) {
854- gCurrentQxZNA ->Fill (valuesQxZNA);
855- gCurrentQyZNA ->Fill (valuesQyZNA);
859+ gCurrentQxZNA ->Fill (valuesQxZNA. data () );
860+ gCurrentQyZNA ->Fill (valuesQyZNA. data () );
856861 }
857862
858863 // Calculate raw/recentered angle
@@ -955,8 +960,8 @@ struct ZdcExtraTableReader {
955960 }
956961 }
957962
958- double valuesQxZNC[ 5 ] = {cent, vx, vy, vz, qxZNCrec};
959- double valuesQyZNC[ 5 ] = {cent, vx, vy, vz, qyZNCrec};
963+ std::array< double , 5 > valuesQxZNC = {cent, vx, vy, vz, qxZNCrec};
964+ std::array< double , 5 > valuesQyZNC = {cent, vx, vy, vz, qyZNCrec};
960965
961966 gCurrentCentroidZNC ->Fill (qxZNCrec, qyZNCrec);
962967
@@ -974,8 +979,8 @@ struct ZdcExtraTableReader {
974979 gCurrentQyVsTimeZNC ->Fill (timeInMinutes, qyZNCrec);
975980
976981 if (plot5D) {
977- gCurrentQxZNC ->Fill (valuesQxZNC);
978- gCurrentQyZNC ->Fill (valuesQyZNC);
982+ gCurrentQxZNC ->Fill (valuesQxZNC. data () );
983+ gCurrentQyZNC ->Fill (valuesQyZNC. data () );
979984 }
980985
981986 // Calculate raw/recentered angle
0 commit comments