@@ -745,20 +745,24 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
745745 if (subGroupStr.Contains (" tpcpid_fine" )) {
746746 // fine binning for pIN: steps in 10 MeV/c from 0 to 1 GeV/c and 100 MeV/c up to 10 GeV/c
747747 std::array<double , 281 > pIN_bins;
748- for (int i = 0 ; i <= 200 ; i++)
748+ for (int i = 0 ; i <= 200 ; i++) {
749749 pIN_bins[i] = 0.01 * i;
750- for (int i = 1 ; i <= 80 ; i++)
750+ }
751+ for (int i = 1 ; i <= 80 ; i++) {
751752 pIN_bins[200 + i] = 2 . + 0.1 * i;
753+ }
752754 int nbins_pIN = static_cast <int >(pIN_bins.size ()) - 1 ;
753755
754756 std::array<double , 201 > TPCdEdx_bins;
755- for (int i = 0 ; i <= 200 ; i++)
757+ for (int i = 0 ; i <= 200 ; i++) {
756758 TPCdEdx_bins[i] = i;
759+ }
757760 int nbins_TPCdEdx = static_cast <int >(TPCdEdx_bins.size ()) - 1 ;
758761
759762 std::array<double , 101 > nSigma_bins;
760- for (int i = 0 ; i <= 100 ; i++)
763+ for (int i = 0 ; i <= 100 ; i++) {
761764 nSigma_bins[i] = -5 . + 0.1 * i;
765+ }
762766 int nbins_nSigma = static_cast <int >(nSigma_bins.size ()) - 1 ;
763767
764768 hm->AddHistogram (histClass, " TPCdedx_pIN" , " TPC dE/dx vs pIN" , false , nbins_pIN, pIN_bins.data (), VarManager::kPin , nbins_TPCdEdx, TPCdEdx_bins.data (), VarManager::kTPCsignal );
@@ -819,8 +823,9 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
819823 const int kNvarsPID = 4 ;
820824 const int kTPCnsigmaNbins = 70 ;
821825 std::array<double , kTPCnsigmaNbins + 1 > tpcNsigmaBinLims;
822- for (int i = 0 ; i <= kTPCnsigmaNbins ; ++i)
826+ for (int i = 0 ; i <= kTPCnsigmaNbins ; ++i) {
823827 tpcNsigmaBinLims[i] = -7.0 + 0.2 * i;
828+ }
824829
825830 const int kPinEleNbins = 20 ;
826831 std::array<double , kPinEleNbins + 1 > pinEleBinLims = {0.1 , 0.2 , 0.4 , 0.6 , 0.8 , 1.0 , 1.2 , 1.5 , 2.0 , 2.5 , 3.0 , 3.5 , 4.0 , 4.5 , 5.0 , 6.0 , 8.0 , 10.0 , 12.0 , 16.0 , 20.0 };
@@ -830,8 +835,9 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
830835
831836 const int kTPCnClusterbins = 16 ;
832837 std::array<double , kTPCnClusterbins + 1 > tpcNclusterBinLims;
833- for (int i = 0 ; i <= kTPCnClusterbins ; ++i)
838+ for (int i = 0 ; i <= kTPCnClusterbins ; ++i) {
834839 tpcNclusterBinLims[i] = 10 * i;
840+ }
835841
836842 std::array<TArrayD, kNvarsPID > nSigBinLimits;
837843 nSigBinLimits[0 ] = TArrayD (kTPCnsigmaNbins + 1 , tpcNsigmaBinLims.data ());
@@ -877,20 +883,24 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
877883 if (subGroupStr.Contains (" tofpid_fine" )) {
878884 // fine binning for pIN: steps in 10 MeV/c from 0 to 1 GeV/c and 100 MeV/c up to 10 GeV/c
879885 std::array<double , 281 > pIN_bins;
880- for (int i = 0 ; i <= 200 ; i++)
886+ for (int i = 0 ; i <= 200 ; i++) {
881887 pIN_bins[i] = 0.01 * i;
882- for (int i = 1 ; i <= 80 ; i++)
888+ }
889+ for (int i = 1 ; i <= 80 ; i++) {
883890 pIN_bins[200 + i] = 2 . + 0.1 * i;
891+ }
884892 int nbins_pIN = static_cast <int >(pIN_bins.size ()) - 1 ;
885893
886894 std::array<double , 241 > TOFbeta_bins;
887- for (int i = 0 ; i <= 240 ; i++)
895+ for (int i = 0 ; i <= 240 ; i++) {
888896 TOFbeta_bins[i] = 0.005 * i;
897+ }
889898 int nbins_TOFbeta = static_cast <int >(TOFbeta_bins.size ()) - 1 ;
890899
891900 std::array<double , 101 > nSigma_bins;
892- for (int i = 0 ; i <= 100 ; i++)
901+ for (int i = 0 ; i <= 100 ; i++) {
893902 nSigma_bins[i] = -5 . + 0.1 * i;
903+ }
894904 int nbins_nSigma = static_cast <int >(nSigma_bins.size ()) - 1 ;
895905
896906 hm->AddHistogram (histClass, " TOFbeta_pIN" , " TOF #beta vs pIN" , false , nbins_pIN, pIN_bins.data (), VarManager::kPin , nbins_TOFbeta, TOFbeta_bins.data (), VarManager::kTOFbeta );
@@ -910,19 +920,24 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
910920 const int kNvarsPID = 3 ;
911921 const int kNbins_pIN = 169 ;
912922 std::array<double , kNbins_pIN + 1 > pIN_bins;
913- for (int i = 0 ; i <= 140 ; i++)
923+ for (int i = 0 ; i <= 140 ; i++) {
914924 pIN_bins[i] = 0.01 * i + 0.1 ;
915- for (int i = 1 ; i <= 15 ; i++)
925+ }
926+ for (int i = 1 ; i <= 15 ; i++) {
916927 pIN_bins[140 + i] = 1.5 + 0.1 * i;
917- for (int i = 1 ; i <= 14 ; i++)
928+ }
929+ for (int i = 1 ; i <= 14 ; i++) {
918930 pIN_bins[155 + i] = 3 . + 0.5 * i;
931+ }
919932
920933 const int kNbins_pINmore = 135 ;
921934 std::array<double , kNbins_pINmore + 1 > pIN_binsmore;
922- for (int i = 0 ; i <= 120 ; i++)
935+ for (int i = 0 ; i <= 120 ; i++) {
923936 pIN_binsmore[i] = 0.01 * i + 0.3 ;
924- for (int i = 1 ; i <= 10 ; i++)
937+ }
938+ for (int i = 1 ; i <= 10 ; i++) {
925939 pIN_binsmore[120 + i] = 1.5 + 0.2 * i;
940+ }
926941 pIN_binsmore[131 ] = 4 .;
927942 pIN_binsmore[132 ] = 5 .;
928943 pIN_binsmore[133 ] = 6 .;
@@ -931,28 +946,33 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
931946
932947 const int kNbins_nSigma = 100 ;
933948 std::array<double , kNbins_nSigma + 1 > nSigma_bins;
934- for (int i = 0 ; i <= kNbins_nSigma ; i++)
949+ for (int i = 0 ; i <= kNbins_nSigma ; i++) {
935950 nSigma_bins[i] = -5 . + 0.1 * i;
951+ }
936952
937953 const int kNbins_nSigmamore = 50 ;
938954 std::array<double , kNbins_nSigmamore + 1 > nSigma_binsmore;
939- for (int i = 0 ; i <= kNbins_nSigmamore ; i++)
955+ for (int i = 0 ; i <= kNbins_nSigmamore ; i++) {
940956 nSigma_binsmore[i] = -5 . + 0.2 * i;
957+ }
941958
942959 const int kNbins_nSigmagrob = 24 ;
943960 std::array<double , kNbins_nSigmagrob + 1 > nSigma_binsgrob;
944- for (int i = 0 ; i <= kNbins_nSigmagrob ; i++)
961+ for (int i = 0 ; i <= kNbins_nSigmagrob ; i++) {
945962 nSigma_binsgrob[i] = -6 . + 0.5 * i;
963+ }
946964
947965 const int kNbins_TOFbeta = 120 ;
948966 std::array<double , kNbins_TOFbeta + 1 > TOFbeta_bins;
949- for (int i = 0 ; i <= kNbins_TOFbeta ; i++)
967+ for (int i = 0 ; i <= kNbins_TOFbeta ; i++) {
950968 TOFbeta_bins[i] = 0.01 * i;
969+ }
951970
952971 const int kNbins_TPCdEdx = 140 ;
953972 std::array<double , kNbins_TPCdEdx + 1 > TPCdEdx_bins;
954- for (int i = 0 ; i <= kNbins_TPCdEdx ; i++)
973+ for (int i = 0 ; i <= kNbins_TPCdEdx ; i++) {
955974 TPCdEdx_bins[i] = i + 20 ;
975+ }
956976
957977 std::array<TArrayD, kNvarsPID > nSigmaBinLimits;
958978 nSigmaBinLimits[0 ] = TArrayD (kNbins_pIN + 1 , pIN_bins.data ());
@@ -1544,8 +1564,9 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
15441564
15451565 const int kTauNBins = 500 ;
15461566 std::array<double , kTauNBins + 1 > TauBinLims;
1547- for (int i = 0 ; i <= kTauNBins ; ++i)
1567+ for (int i = 0 ; i <= kTauNBins ; ++i) {
15481568 TauBinLims[i] = -0.3 + (0.0015 * i);
1569+ }
15491570
15501571 std::array<TArrayD, kNvarsPair > nCutsBinLimits;
15511572 nCutsBinLimits[0 ] = TArrayD (kInvMassNbins + 1 , InvMassBinLims.data ());
@@ -2059,19 +2080,23 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
20592080 // every 100 MeV from 0.2 to 1. GeV/c2
20602081 // every 500 GeV from 1 to 5 GeV/c2
20612082 std::array<double , 37 > mee_bins;
2062- for (int i = 0 ; i <= 20 ; i++)
2083+ for (int i = 0 ; i <= 20 ; i++) {
20632084 mee_bins[i] = 0.01 * i;
2064- for (int i = 1 ; i <= 8 ; i++)
2085+ }
2086+ for (int i = 1 ; i <= 8 ; i++) {
20652087 mee_bins[20 + i] = 0.2 + 0.1 * i;
2066- for (int i = 1 ; i <= 8 ; i++)
2088+ }
2089+ for (int i = 1 ; i <= 8 ; i++) {
20672090 mee_bins[28 + i] = 1 . + 0.5 * i;
2091+ }
20682092 int nbins_mee = static_cast <int >(mee_bins.size ()) - 1 ;
20692093
20702094 // binning for ptee at large scales:
20712095 // every 0.2 GeV/c from 0 to 10 GeV/c
20722096 std::array<double , 51 > ptee_bins;
2073- for (int i = 0 ; i <= 50 ; i++)
2097+ for (int i = 0 ; i <= 50 ; i++) {
20742098 ptee_bins[i] = 0.2 * i;
2099+ }
20752100 int nbins_ptee = static_cast <int >(ptee_bins.size ()) - 1 ;
20762101
20772102 // binning for phiv:
@@ -2103,37 +2128,46 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
21032128 // every 10 MeV from 2.7 to 3.2 GeV/c2
21042129 // every 50 MeV from 3.2 to 12 GeV/c2
21052130 std::array<double , 369 > mee_bins;
2106- for (int i = 0 ; i <= 110 ; i++)
2131+ for (int i = 0 ; i <= 110 ; i++) {
21072132 mee_bins[i] = 0.01 * i;
2108- for (int i = 1 ; i <= 32 ; i++)
2133+ }
2134+ for (int i = 1 ; i <= 32 ; i++) {
21092135 mee_bins[110 + i] = 1.1 + 0.05 * i;
2110- for (int i = 1 ; i <= 50 ; i++)
2136+ }
2137+ for (int i = 1 ; i <= 50 ; i++) {
21112138 mee_bins[142 + i] = 2.7 + 0.01 * i;
2112- for (int i = 1 ; i <= 176 ; i++)
2139+ }
2140+ for (int i = 1 ; i <= 176 ; i++) {
21132141 mee_bins[192 + i] = 3.2 + 0.05 * i;
2142+ }
21142143 int nbins_mee = static_cast <int >(mee_bins.size ()) - 1 ;
21152144
21162145 // binning for ptee at large scales:
21172146 // every 0.1 GeV/c from 0 to 10 GeV/c
21182147 // every 0.5 GeV/c from 10 to 30 GeV/c
21192148 std::array<double , 201 > ptee_bins;
2120- for (int i = 0 ; i <= 100 ; i++)
2149+ for (int i = 0 ; i <= 100 ; i++) {
21212150 ptee_bins[i] = 0.1 * i;
2122- for (int i = 1 ; i <= 100 ; i++)
2151+ }
2152+ for (int i = 1 ; i <= 100 ; i++) {
21232153 ptee_bins[100 + i] = 10 + 0.2 * i;
2154+ }
21242155 int nbins_ptee = static_cast <int >(ptee_bins.size ()) - 1 ;
21252156
21262157 // binning for dca at large scales:
21272158 // every 0.1 sigma from 0 to 5 sigma
21282159 // every 0.5 sigma from 5 to 10 sigma
21292160 // every 1.0 sigma from 10 to 40 sigma
21302161 std::array<double , 91 > dca_bins;
2131- for (int i = 0 ; i <= 50 ; i++)
2162+ for (int i = 0 ; i <= 50 ; i++) {
21322163 dca_bins[i] = 0.1 * i;
2133- for (int i = 1 ; i <= 10 ; i++)
2164+ }
2165+ for (int i = 1 ; i <= 10 ; i++) {
21342166 dca_bins[50 + i] = 5 + 0.5 * i;
2135- for (int i = 1 ; i <= 30 ; i++)
2167+ }
2168+ for (int i = 1 ; i <= 30 ; i++) {
21362169 dca_bins[60 + i] = 10 + 1 * i;
2170+ }
21372171 int nbins_dca = static_cast <int >(dca_bins.size ()) - 1 ;
21382172
21392173 hm->AddHistogram (histClass, " Mass_QuadDCAsigXY" , " " , false , nbins_mee, mee_bins.data (), VarManager::kMass , nbins_dca, dca_bins.data (), VarManager::kQuadDCAsigXY );
@@ -2218,8 +2252,9 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
22182252 const int kNvarsTripletCuts = 4 ;
22192253 const int kInvMassNbins = 100 ;
22202254 std::array<double , kInvMassNbins + 1 > InvMassBinLims;
2221- for (int i = 0 ; i <= kInvMassNbins ; ++i)
2255+ for (int i = 0 ; i <= kInvMassNbins ; ++i) {
22222256 InvMassBinLims[i] = 4.0 + 0.02 * i;
2257+ }
22232258
22242259 const int kPtNbins = 6 ;
22252260 std::array<double , kPtNbins + 1 > PtBinLims = {0 ., 2 ., 4 ., 6 ., 8 ., 10 ., 20 .};
@@ -2245,13 +2280,15 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
22452280 if (subGroupStr.Contains (" dilepton-hadron-array-correlation" )) {
22462281 const int kInvMassBins = 500 ;
22472282 std::array<double , kInvMassBins + 1 > InvMassBinLims;
2248- for (int i = 0 ; i <= kInvMassBins ; i++)
2283+ for (int i = 0 ; i <= kInvMassBins ; i++) {
22492284 InvMassBinLims[i] = 0 + i * 0.01 ;
2285+ }
22502286
22512287 const int kDelEtaBins = 20 ;
22522288 std::array<double , kDelEtaBins + 1 > DelEtaBinLims;
2253- for (int i = 0 ; i <= kDelEtaBins ; i++)
2289+ for (int i = 0 ; i <= kDelEtaBins ; i++) {
22542290 DelEtaBinLims[i] = -2 + i * 0.2 ;
2291+ }
22552292
22562293 const int kDelPhiBins = 26 ;
22572294 std::array<double , kDelPhiBins + 1 > DelPhiBinLims = {-1.69647 , -1.44513 , -1.19381 , -0.94248 , -0.69115 , -0.43982 , -0.18850 , 0.06283 , 0.31416 , 0.56549 , 0.81681 , 1.06814 , 1.31947 , 1.57080 , 1.82212 , 2.07345 , 2.32478 , 2.57611 , 2.82743 , 3.07876 , 3.33009 , 3.58142 , 3.83274 , 4.08407 , 4.33540 , 4.58673 , 4.83806 };
0 commit comments