Skip to content

Commit e18f2f8

Browse files
nidamalikknidamalikk
committed
use different polynomial for fitting multiplicity correlation plots
2 parents 4fda1aa + f2d8005 commit e18f2f8

File tree

58 files changed

+4475
-1117
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

58 files changed

+4475
-1117
lines changed

‎ALICE3/macros/a3geo.ini‎

Lines changed: 148 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,148 @@
1+
# Example of the ALICE3 geometry file in TEnv format
2+
bpipe0.r: 0.48
3+
bpipe0.z: 250
4+
bpipe0.x0: 0.001592
5+
bpipe0.xrho: 0
6+
bpipe0.resRPhi: 0.0
7+
bpipe0.resZ: 0.0
8+
bpipe0.eff: 0.0
9+
bpipe0.type: 0
10+
bpipe0.deadPhiRegions: /tmp/asd.root
11+
12+
B00.r: 0.5
13+
B00.z: 250
14+
B00.x0: 0.00076
15+
B00.xrho: 0
16+
B00.resRPhi: 0.00025
17+
B00.resZ: 0.00025
18+
B00.eff: 1.
19+
B00.type: 1
20+
21+
B01.r: 1.2
22+
B01.z: 250
23+
B01.x0: 0.00096
24+
B01.xrho: 0
25+
B01.resRPhi: 0.00025
26+
B01.resZ: 0.00025
27+
B01.eff: 1.
28+
B01.type: 1
29+
30+
B02.r: 2.5
31+
B02.z: 250
32+
B02.x0: 0.00096
33+
B02.xrho: 0
34+
B02.resRPhi: 0.00025
35+
B02.resZ: 0.00025
36+
B02.eff: 1.
37+
B02.type: 1
38+
39+
coldplate.r: 2.6
40+
coldplate.z: 250
41+
coldplate.x0: 0.02
42+
coldplate.xrho: 0
43+
coldplate.resRPhi: 0
44+
coldplate.resZ: 0
45+
coldplate.eff: 0
46+
coldplate.type: 0
47+
48+
petal.r: 3.7
49+
petal.z: 250
50+
petal.x0: 0.001592
51+
petal.xrho: 0
52+
petal.resRPhi: 0
53+
petal.resZ: 0
54+
petal.eff: 0
55+
petal.type: 0
56+
57+
bpipe1.r: 5.7
58+
bpipe1.z: 250
59+
bpipe1.x0: 0.0023
60+
bpipe1.xrho: 0
61+
bpipe1.resRPhi: 0.0
62+
bpipe1.resZ: 0.0
63+
bpipe1.eff: 0.0
64+
bpipe1.type: 0
65+
66+
B03.r: 7.
67+
B03.z: 250
68+
B03.x0: 0.01
69+
B03.xrho: 0
70+
B03.resRPhi: 0.001
71+
B03.resZ: 0.001
72+
B03.eff: 1.
73+
B03.type: 1
74+
75+
B04.r: 9.
76+
B04.z: 250
77+
B04.x0: 0.01
78+
B04.xrho: 0
79+
B04.resRPhi: 0.001
80+
B04.resZ: 0.001
81+
B04.eff: 1.
82+
B04.type: 1
83+
84+
B05.r: 12.
85+
B05.z: 250
86+
B05.x0: 0.01
87+
B05.xrho: 0
88+
B05.resRPhi: 0.001
89+
B05.resZ: 0.001
90+
B05.eff: 1.
91+
B05.type: 1
92+
93+
iTOF.r: 19.
94+
iTOF.z: 250
95+
iTOF.x0: 0.03
96+
iTOF.xrho: 0
97+
iTOF.resRPhi: 0.001
98+
iTOF.resZ: 0.001
99+
iTOF.eff: 1.
100+
iTOF.type: 0
101+
102+
B06.r: 20.
103+
B06.z: 250
104+
B06.x0: 0.01
105+
B06.xrho: 0
106+
B06.resRPhi: 0.001
107+
B06.resZ: 0.001
108+
B06.eff: 1.
109+
B06.type: 1
110+
111+
B07.r: 30.
112+
B07.z: 250
113+
B07.x0: 0.01
114+
B07.xrho: 0
115+
B07.resRPhi: 0.001
116+
B07.resZ: 0.001
117+
B07.eff: 1.
118+
B07.type: 1
119+
120+
B08.r: 45.
121+
B08.z: 250
122+
B08.x0: 0.01
123+
B08.xrho: 0
124+
B08.resRPhi: 0.001
125+
B08.resZ: 0.001
126+
B08.eff: 1.
127+
B08.type: 1
128+
129+
B09.r: 60.
130+
B09.z: 250
131+
B09.x0: 0.01
132+
B09.xrho: 0
133+
B09.resRPhi: 0.001
134+
B09.resZ: 0.001
135+
B09.eff: 1.
136+
B09.type: 1
137+
138+
B10.r: 80.
139+
B10.z: 250
140+
B10.x0: 0.01
141+
B10.xrho: 0
142+
B10.resRPhi: 0.001
143+
B10.resZ: 0.001
144+
B10.eff: 1.
145+
B10.type: 1
146+
147+
148+

‎CODEOWNERS‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -61,7 +61,7 @@
6161
/PWGMM/Lumi @alibuild @sustripathy @skundu692 @aalkin @jgcn @gbencedi @abmodak
6262
/PWGMM/UE @alibuild @sustripathy @skundu692 @aalkin @aortizve @jgcn
6363

64-
/PWGUD @alibuild @pbuehler @nystrand @rolavick
64+
/PWGUD @alibuild @amatyja @rolavick
6565
/PWGJE @alibuild @lhavener @maoyx @nzardosh @fjonasALICE @mfasDa @mhemmer-cern
6666
/Tools/PIDML @alibuild @saganatt
6767
/Tools/ML @alibuild @fcatalan92 @fmazzasc
@@ -72,4 +72,4 @@
7272
/Tutorials/PWGJE @alibuild @lhavener @maoyx @nzardosh @mfasDa @fjonasALICE
7373
/Tutorials/PWGLF @alibuild @alcaliva @lbariogl @chiarapinto @BongHwi @lbarnby @ercolessi @iravasen @njacazio @romainschotter @skundu692
7474
/Tutorials/PWGMM @alibuild @aalkin @ddobrigk
75-
/Tutorials/PWGUD @alibuild @pbuehler
75+
/Tutorials/PWGUD @alibuild @pbuehler @amatyja

‎Common/CCDB/RCTSelectionFlags.h‎

Lines changed: 36 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -75,7 +75,15 @@ enum RCTSelectionFlags {
7575
kTPCLimAccMCRepr,
7676
kTRDBad,
7777
kZDCBad,
78-
kNRCTSelectionFlags
78+
kNRCTSelectionFlags,
79+
kDummy24,
80+
kDummy25,
81+
kDummy26,
82+
kDummy27,
83+
kDummy28,
84+
kDummy29,
85+
kDummy30,
86+
kCcdbObjectLoaded
7987
};
8088

8189
template <typename T>
@@ -104,9 +112,14 @@ class RCTFlagsChecker : public o2::utils::EnumFlags<RCTSelectionFlags>
104112
// The checkZDC boolean flag controls whether to iclude the ZDC quality in all the pre-defined selections (for Pb-Pb data)
105113
// The treatLimitedAcceptanceAsBad boolean flag controls whether "LimitedAcceptanceMCReproducible" flags should be
106114
// treated as Bad and the corresponding events excluded
107-
explicit RCTFlagsChecker(const std::string& label, bool checkZDC = false, bool treatLimitedAcceptanceAsBad = false)
115+
// The checkTableValidity boolean flag controls whether events without a corresponding valid RCT CCDB object should be
116+
// treated as Bad and excluded
117+
explicit RCTFlagsChecker(const std::string& label,
118+
bool checkZDC = false,
119+
bool treatLimitedAcceptanceAsBad = false,
120+
bool checkTableValidity = false)
108121
{
109-
init(label, checkZDC, treatLimitedAcceptanceAsBad);
122+
init(label, checkZDC, treatLimitedAcceptanceAsBad, checkTableValidity);
110123
}
111124

112125
// Initialize the object from an initializer list of RCTSelectionFlags values
@@ -121,13 +134,18 @@ class RCTFlagsChecker : public o2::utils::EnumFlags<RCTSelectionFlags>
121134
// - "CBT"
122135
// - "CBT_hadronPID"
123136
// - "CBT_electronPID"
124-
// - "CBT_calo"
137+
// - "CCBT_calo"
125138
// - "CBT_muon"
126139
// - "CBT_muon_glo"
127140
// The checkZDC boolean flag controls whether to iclude the ZDC quality in all the pre-defined selections (for Pb-Pb data)
128141
// The treatLimitedAcceptanceAsBad boolean flag controls whether "LimitedAcceptanceMCReproducible" flags should be
129142
// treated as Bad and the corresponding events excluded
130-
void init(const std::string& label, bool checkZDC = false, bool treatLimitedAcceptanceAsBad = false)
143+
// The checkTableValidity boolean flag controls whether events without a corresponding valid RCT CCDB object should be
144+
// treated as Bad and excluded
145+
void init(const std::string& label,
146+
bool checkZDC = false,
147+
bool treatLimitedAcceptanceAsBad = false,
148+
bool checkTableValidity = false)
131149
{
132150
auto setFlags = [this](std::initializer_list<RCTSelectionFlags> flags) {
133151
std::for_each(flags.begin(),
@@ -182,14 +200,19 @@ class RCTFlagsChecker : public o2::utils::EnumFlags<RCTSelectionFlags>
182200
if (checkZDC) {
183201
set(kZDCBad);
184202
}
203+
204+
if (checkTableValidity) {
205+
set(kCcdbObjectLoaded);
206+
}
185207
}
186208

187209
// Check the RCT column of a given event selection table.
188210
// The function returns true if none of the checked flags is set in the RCT column.
189211
bool checkTable(const HasRCTFlags auto& table)
190212
{
213+
// throw an exception if none of the bits in the checker mask is set
191214
if (!any()) {
192-
throw std::out_of_range("RCTFlagsCheckerAlt with empty RCTSelectionFlags bits mask");
215+
throw std::out_of_range("RCTFlagsChecker has empty RCTSelectionFlags bits mask");
193216
}
194217

195218
// bitmask of the current table
@@ -201,6 +224,13 @@ class RCTFlagsChecker : public o2::utils::EnumFlags<RCTSelectionFlags>
201224
return ((tableBits & flagsBits) == 0);
202225
}
203226

227+
// Check the validity of the RCT column of a given event selection table.
228+
// The function returns true if the RCT CCDB object was correctly fetched (bit31 == 0).
229+
bool isTableValid(const HasRCTFlags auto& table)
230+
{
231+
return (table.rct_bit(kCcdbObjectLoaded) == 0);
232+
}
233+
204234
bool operator()(const HasRCTFlags auto& table)
205235
{
206236
return checkTable(table);

‎DPG/Tasks/AOTEvent/eventSelectionQa.cxx‎

Lines changed: 16 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -999,20 +999,24 @@ struct EventSelectionQaTask {
999999
for (const auto& track : tracks) {
10001000
auto mapAmbTrIdsIt = mapAmbTrIds.find(track.globalIndex());
10011001
int ambTrId = mapAmbTrIdsIt == mapAmbTrIds.end() ? -1 : mapAmbTrIdsIt->second;
1002+
1003+
// special check to avoid crashes (in particular, on some MC Pb-Pb datasets)
1004+
// (related to shifts in ambiguous tracks association to bc slices (off by 1) - see https://mattermost.web.cern.ch/alice/pl/g9yaaf3tn3g4pgn7c1yex9copy
1005+
if (ambTrId >= 0 && (ambTracks.iteratorAt(ambTrId).bcIds()[0] >= bcs.size()))
1006+
continue;
1007+
10021008
int indexBc = ambTrId < 0 ? track.collision_as<ColEvSels>().bc_as<BCsRun3>().globalIndex() : ambTracks.iteratorAt(ambTrId).bc_as<BCsRun3>().begin().globalIndex();
1003-
if (ambTrId < 0) { // temprorary limitation, to avoid crashes, in particular, on MC Pb-Pb datasets
1004-
auto bc = bcs.iteratorAt(indexBc);
1005-
int64_t globalBC = bc.globalBC() + floor(track.trackTime() / o2::constants::lhc::LHCBunchSpacingNS);
1009+
auto bc = bcs.iteratorAt(indexBc);
1010+
int64_t globalBC = bc.globalBC() + floor(track.trackTime() / o2::constants::lhc::LHCBunchSpacingNS);
10061011

1007-
int32_t indexClosestTVX = findClosest(globalBC, mapGlobalBcWithTVX);
1008-
int bcDiff = static_cast<int>(globalBC - vGlobalBCs[indexClosestTVX]);
1009-
if (track.hasTOF() || track.hasTRD() || !track.hasITS() || !track.hasTPC() || track.pt() < 1)
1010-
continue;
1011-
histos.fill(HIST("hTrackBcDiffVsEtaAll"), track.eta(), bcDiff);
1012-
if (track.eta() < -0.2 || track.eta() > 0.2)
1013-
continue;
1014-
histos.fill(HIST("hSecondsTVXvsBcDifAll"), bc.timestamp() / 1000., bcDiff);
1015-
}
1012+
int32_t indexClosestTVX = findClosest(globalBC, mapGlobalBcWithTVX);
1013+
int bcDiff = static_cast<int>(globalBC - vGlobalBCs[indexClosestTVX]);
1014+
if (track.hasTOF() || track.hasTRD() || !track.hasITS() || !track.hasTPC() || track.pt() < 1)
1015+
continue;
1016+
histos.fill(HIST("hTrackBcDiffVsEtaAll"), track.eta(), bcDiff);
1017+
if (track.eta() < -0.2 || track.eta() > 0.2)
1018+
continue;
1019+
histos.fill(HIST("hSecondsTVXvsBcDifAll"), bc.timestamp() / 1000., bcDiff);
10161020
}
10171021

10181022
// collision-based event selection qa

‎DPG/Tasks/AOTEvent/timeDependentQa.cxx‎

Lines changed: 42 additions & 39 deletions
Original file line numberDiff line numberDiff line change
@@ -305,39 +305,40 @@ struct TimeDependentQaTask {
305305

306306
// ### QA RCT flags
307307
int nRctFlagsTotal = enNumRctFlagsTotal;
308-
histos.add("hSecondsRCTflags", "", kTH2F, {axisSecondsWideBins, {nRctFlagsTotal + 1, -0.5, nRctFlagsTotal + 1 - 0.5, "Monitoring of RCT flags"}});
308+
histos.add("hSecondsRCTflags", "", kTH2F, {axisSecondsWideBins, {nRctFlagsTotal + 2, -0.5, nRctFlagsTotal + 2 - 0.5, "Monitoring of RCT flags"}});
309309
axRctFlags = reinterpret_cast<TAxis*>(histos.get<TH2>(HIST("hSecondsRCTflags"))->GetYaxis());
310310
axRctFlags->SetBinLabel(1, "NcollisionsSel8");
311-
axRctFlags->SetBinLabel(2 + kCPVBad, "CPVBad");
312-
axRctFlags->SetBinLabel(2 + kEMCBad, "EMCBad");
313-
axRctFlags->SetBinLabel(2 + kEMCLimAccMCRepr, "EMCLimAccMCRepr");
314-
axRctFlags->SetBinLabel(2 + kFDDBad, "FDDBad");
315-
axRctFlags->SetBinLabel(2 + kFT0Bad, "FT0Bad");
316-
axRctFlags->SetBinLabel(2 + kFV0Bad, "FV0Bad");
317-
axRctFlags->SetBinLabel(2 + kHMPBad, "HMPBad");
318-
axRctFlags->SetBinLabel(2 + kITSBad, "ITSBad");
319-
axRctFlags->SetBinLabel(2 + kITSLimAccMCRepr, "ITSLimAccMCRepr");
320-
axRctFlags->SetBinLabel(2 + kMCHBad, "MCHBad");
321-
axRctFlags->SetBinLabel(2 + kMCHLimAccMCRepr, "MCHLimAccMCRepr");
322-
axRctFlags->SetBinLabel(2 + kMFTBad, "MFTBad");
323-
axRctFlags->SetBinLabel(2 + kMFTLimAccMCRepr, "MFTLimAccMCRepr");
324-
axRctFlags->SetBinLabel(2 + kMIDBad, "MIDBad");
325-
axRctFlags->SetBinLabel(2 + kMIDLimAccMCRepr, "MIDLimAccMCRepr");
326-
axRctFlags->SetBinLabel(2 + kPHSBad, "PHSBad");
327-
axRctFlags->SetBinLabel(2 + kTOFBad, "TOFBad");
328-
axRctFlags->SetBinLabel(2 + kTOFLimAccMCRepr, "TOFLimAccMCRepr");
329-
axRctFlags->SetBinLabel(2 + kTPCBadTracking, "TPCBadTracking");
330-
axRctFlags->SetBinLabel(2 + kTPCBadPID, "TPCBadPID");
331-
axRctFlags->SetBinLabel(2 + kTPCLimAccMCRepr, "TPCLimAccMCRepr");
332-
axRctFlags->SetBinLabel(2 + kTRDBad, "TRDBad");
333-
axRctFlags->SetBinLabel(2 + kZDCBad, "ZDCBad");
311+
axRctFlags->SetBinLabel(2, "CcdbNotFound");
312+
axRctFlags->SetBinLabel(3 + kCPVBad, "CPVBad");
313+
axRctFlags->SetBinLabel(3 + kEMCBad, "EMCBad");
314+
axRctFlags->SetBinLabel(3 + kEMCLimAccMCRepr, "EMCLimAccMCRepr");
315+
axRctFlags->SetBinLabel(3 + kFDDBad, "FDDBad");
316+
axRctFlags->SetBinLabel(3 + kFT0Bad, "FT0Bad");
317+
axRctFlags->SetBinLabel(3 + kFV0Bad, "FV0Bad");
318+
axRctFlags->SetBinLabel(3 + kHMPBad, "HMPBad");
319+
axRctFlags->SetBinLabel(3 + kITSBad, "ITSBad");
320+
axRctFlags->SetBinLabel(3 + kITSLimAccMCRepr, "ITSLimAccMCRepr");
321+
axRctFlags->SetBinLabel(3 + kMCHBad, "MCHBad");
322+
axRctFlags->SetBinLabel(3 + kMCHLimAccMCRepr, "MCHLimAccMCRepr");
323+
axRctFlags->SetBinLabel(3 + kMFTBad, "MFTBad");
324+
axRctFlags->SetBinLabel(3 + kMFTLimAccMCRepr, "MFTLimAccMCRepr");
325+
axRctFlags->SetBinLabel(3 + kMIDBad, "MIDBad");
326+
axRctFlags->SetBinLabel(3 + kMIDLimAccMCRepr, "MIDLimAccMCRepr");
327+
axRctFlags->SetBinLabel(3 + kPHSBad, "PHSBad");
328+
axRctFlags->SetBinLabel(3 + kTOFBad, "TOFBad");
329+
axRctFlags->SetBinLabel(3 + kTOFLimAccMCRepr, "TOFLimAccMCRepr");
330+
axRctFlags->SetBinLabel(3 + kTPCBadTracking, "TPCBadTracking");
331+
axRctFlags->SetBinLabel(3 + kTPCBadPID, "TPCBadPID");
332+
axRctFlags->SetBinLabel(3 + kTPCLimAccMCRepr, "TPCLimAccMCRepr");
333+
axRctFlags->SetBinLabel(3 + kTRDBad, "TRDBad");
334+
axRctFlags->SetBinLabel(3 + kZDCBad, "ZDCBad");
334335
// combined flags
335-
axRctFlags->SetBinLabel(2 + enCBT, "CBT");
336-
axRctFlags->SetBinLabel(2 + enCBT_hadronPID, "CBT_hadronPID");
337-
axRctFlags->SetBinLabel(2 + enCBT_electronPID, "CBT_electronPID");
338-
axRctFlags->SetBinLabel(2 + enCBT_calo, "CBT_calo");
339-
axRctFlags->SetBinLabel(2 + enCBT_muon, "CBT_muon");
340-
axRctFlags->SetBinLabel(2 + enCBT_muon_glo, "CBT_muon_glo");
336+
axRctFlags->SetBinLabel(3 + enCBT, "CBT");
337+
axRctFlags->SetBinLabel(3 + enCBT_hadronPID, "CBT_hadronPID");
338+
axRctFlags->SetBinLabel(3 + enCBT_electronPID, "CBT_electronPID");
339+
axRctFlags->SetBinLabel(3 + enCBT_calo, "CBT_calo");
340+
axRctFlags->SetBinLabel(3 + enCBT_muon, "CBT_muon");
341+
axRctFlags->SetBinLabel(3 + enCBT_muon_glo, "CBT_muon_glo");
341342

342343
// QA for all tracks
343344
// const AxisSpec axisChi2ITS{40, 0., 20., "chi2/ndof"};
@@ -568,18 +569,20 @@ struct TimeDependentQaTask {
568569

569570
// check RCT flags
570571
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 0); // n collisions sel8
572+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1, col.rct_bit(kCcdbObjectLoaded)); // CCDB object not loaded
573+
LOGP(debug, "i = 1, bitValue = {}, binLabel={}, binCenter={}", col.rct_bit(kCcdbObjectLoaded), axRctFlags->GetBinLabel(2), axRctFlags->GetBinCenter(2));
571574
for (int iFlag = 0; iFlag < kNRCTSelectionFlags; iFlag++) {
572-
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1 + iFlag, col.rct_bit(iFlag));
573-
LOGP(debug, "i = {}, bitValue = {}, binLabel={}, binCenter={}", iFlag, col.rct_bit(iFlag), axRctFlags->GetBinLabel(2 + iFlag), axRctFlags->GetBinCenter(2 + iFlag));
575+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 2 + iFlag, col.rct_bit(iFlag));
576+
LOGP(debug, "i = {}, bitValue = {}, binLabel={}, binCenter={}", iFlag, col.rct_bit(iFlag), axRctFlags->GetBinLabel(3 + iFlag), axRctFlags->GetBinCenter(3 + iFlag));
574577
}
575578
LOGP(debug, "CBT_hadronPID = {}, kFT0Bad = {}, kITSBad = {}, kTPCBadTracking = {}, kTPCBadPID = {}, kTOFBad = {}, 1 + enCBT_hadronPID = {}, binLabel={}, binCenter={}", rctCheckerCBT_hadronPID(col),
576-
col.rct_bit(kFT0Bad), col.rct_bit(kITSBad), col.rct_bit(kTPCBadTracking), col.rct_bit(kTPCBadPID), col.rct_bit(kTOFBad), 1 + enCBT_hadronPID, axRctFlags->GetBinLabel(2 + enCBT_hadronPID), axRctFlags->GetBinCenter(2 + enCBT_hadronPID));
577-
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1 + enCBT, rctCheckerCBT(col));
578-
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1 + enCBT_hadronPID, rctCheckerCBT_hadronPID(col));
579-
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1 + enCBT_electronPID, rctCheckerCBT_electronPID(col));
580-
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1 + enCBT_calo, rctCheckerCBT_calo(col));
581-
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1 + enCBT_muon, rctCheckerCBT_muon(col));
582-
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 1 + enCBT_muon_glo, rctCheckerCBT_muon_glo(col));
579+
col.rct_bit(kFT0Bad), col.rct_bit(kITSBad), col.rct_bit(kTPCBadTracking), col.rct_bit(kTPCBadPID), col.rct_bit(kTOFBad), 2 + enCBT_hadronPID, axRctFlags->GetBinLabel(3 + enCBT_hadronPID), axRctFlags->GetBinCenter(3 + enCBT_hadronPID));
580+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 2 + enCBT, rctCheckerCBT(col));
581+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 2 + enCBT_hadronPID, rctCheckerCBT_hadronPID(col));
582+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 2 + enCBT_electronPID, rctCheckerCBT_electronPID(col));
583+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 2 + enCBT_calo, rctCheckerCBT_calo(col));
584+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 2 + enCBT_muon, rctCheckerCBT_muon(col));
585+
histos.fill(HIST("hSecondsRCTflags"), secFromSOR, 2 + enCBT_muon_glo, rctCheckerCBT_muon_glo(col));
583586

584587
// check hadronic rate
585588
double hadronicRate = mRateFetcher.fetch(ccdb.service, ts, runNumber, "ZNC hadronic") * 1.e-3; // kHz

0 commit comments

Comments
 (0)