Skip to content

Commit 7bbd9da

Browse files
creetz16creetz16
authored
[PWGLF] fix histogram registry in decay3bodybuilder (#9416)
1 parent 4339d46 commit 7bbd9da

File tree

1 file changed

+32
-30
lines changed

1 file changed

+32
-30
lines changed

‎PWGLF/TableProducer/Nuspex/decay3bodybuilder.cxx‎

Lines changed: 32 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -952,7 +952,7 @@ struct decay3bodyBuilder {
952952
registry.fill(HIST("QA/Event/hVtxCovYZ"), collision.covYZ());
953953
}
954954

955-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxAll);
955+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxAll);
956956

957957
auto trackParCovPos = getTrackParCov(trackPos);
958958
auto trackParCovNeg = getTrackParCov(trackNeg);
@@ -996,13 +996,13 @@ struct decay3bodyBuilder {
996996
if (trackPos.sign() != +1 || trackNeg.sign() != -1) {
997997
return;
998998
}
999-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxCharge);
999+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxCharge);
10001000

10011001
// track eta
10021002
if (abs(trackPos.eta()) > kfparticleConfigurations.maxEta || abs(trackNeg.eta()) > kfparticleConfigurations.maxEta || abs(trackBach.eta()) > kfparticleConfigurations.maxEta) {
10031003
return;
10041004
}
1005-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxEta);
1005+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxEta);
10061006

10071007
// number of TPC clusters
10081008
if (trackBach.tpcNClsFound() <= kfparticleConfigurations.mintpcNClsBach) {
@@ -1013,7 +1013,7 @@ struct decay3bodyBuilder {
10131013
} else if (!isMatter && ((kfparticleConfigurations.useTPCforPion && trackPos.tpcNClsFound() <= kfparticleConfigurations.mintpcNClsPion) || trackNeg.tpcNClsFound() <= kfparticleConfigurations.mintpcNClsProton)) {
10141014
return;
10151015
}
1016-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxTPCNcls);
1016+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxTPCNcls);
10171017

10181018
// number of TPC crossed rows
10191019
if (trackBach.tpcNClsCrossedRows() <= kfparticleConfigurations.mintpcCrossedRows) {
@@ -1024,7 +1024,7 @@ struct decay3bodyBuilder {
10241024
} else if (!isMatter && ((kfparticleConfigurations.useTPCforPion && trackPos.tpcNClsCrossedRows() <= kfparticleConfigurations.mintpcCrossedRowsPion) || trackNeg.tpcNClsCrossedRows() <= kfparticleConfigurations.mintpcCrossedRows)) {
10251025
return;
10261026
}
1027-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxTPCRows);
1027+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxTPCRows);
10281028

10291029
// TPC PID
10301030
float tpcNsigmaProton;
@@ -1051,7 +1051,7 @@ struct decay3bodyBuilder {
10511051
return;
10521052
}
10531053
}
1054-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxTPCPID);
1054+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxTPCPID);
10551055
LOG(debug) << "Basic track selections done.";
10561056

10571057
// TOF PID of deuteron (set motherhyp correctly)
@@ -1105,13 +1105,13 @@ struct decay3bodyBuilder {
11051105
} else if (!isMatter && (fabs(TrackPosDcaXY) <= kfparticleConfigurations.mindcaXYPionPV || fabs(TrackNegDcaXY) <= kfparticleConfigurations.mindcaXYProtonPV)) {
11061106
return;
11071107
}
1108-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxDCAxyPV);
1108+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxDCAxyPV);
11091109
if (isMatter && (fabs(TrackNegDcaZ) <= kfparticleConfigurations.mindcaZPionPV || fabs(TrackPosDcaZ) <= kfparticleConfigurations.mindcaZProtonPV)) {
11101110
return;
11111111
} else if (!isMatter && (fabs(TrackPosDcaZ) <= kfparticleConfigurations.mindcaZPionPV || fabs(TrackNegDcaZ) <= kfparticleConfigurations.mindcaZProtonPV)) {
11121112
return;
11131113
}
1114-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxDCAzPV);
1114+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxDCAzPV);
11151115
// SVertexer selection bachelor track for event mixing
11161116
if (isEventMixing && kfparticleConfigurations.applySVertexerCuts && TrackBachDca < 0.05) {
11171117
return;
@@ -1181,7 +1181,7 @@ struct decay3bodyBuilder {
11811181
LOG(debug) << "V0 mass constraint applied.";
11821182
KFV0 = KFV0Mass;
11831183
}
1184-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxV0MassConst);
1184+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxV0MassConst);
11851185

11861186
// -------- STEP 3: fit three body vertex --------
11871187
// Create KFParticle object from deuteron track
@@ -1204,6 +1204,7 @@ struct decay3bodyBuilder {
12041204
kfPionDaughters.push_back(kfpPion); // transported to HtVtx
12051205
kfDeuteronDaughters.push_back(kfpDeuteron); // transported to HtVtx
12061206
kf3bodyVtxCandidates.push_back(KFHt);
1207+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxhasSV);
12071208
}
12081209
} else if (kfparticleConfigurations.rotateDeuteron) { // Case 2
12091210
for (int i = 0; i < nRotations; i++) {
@@ -1215,6 +1216,7 @@ struct decay3bodyBuilder {
12151216
kfPionDaughters.push_back(kfpPion); // transported to HtVtx
12161217
kfDeuteronDaughters.push_back(kfpDeuteron); // rotated and transported to HtVtx
12171218
kf3bodyVtxCandidates.push_back(KFHt);
1219+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxhasSV);
12181220
}
12191221
} else { // Case 3
12201222
KFParticle KFHt;
@@ -1223,7 +1225,7 @@ struct decay3bodyBuilder {
12231225
kfPionDaughters.push_back(kfpPion); // transported to HtVtx
12241226
kfDeuteronDaughters.push_back(kfpDeuteron); // transported to HtVtx
12251227
kf3bodyVtxCandidates.push_back(KFHt);
1226-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxhasSV);
1228+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxhasSV);
12271229
}
12281230

12291231
for (int i = 0; i < nRotations; i++) {
@@ -1234,21 +1236,21 @@ struct decay3bodyBuilder {
12341236
continue;
12351237
}
12361238
float DCAvtxDaughters3D = kfProtonDaughters[i].GetDistanceFromParticle(kfPionDaughters[i]) + kfProtonDaughters[i].GetDistanceFromParticle(kfDeuteronDaughters[i]) + kfPionDaughters[i].GetDistanceFromParticle(kfDeuteronDaughters[i]);
1237-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxDcaDau);
1239+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxDcaDau);
12381240
LOG(debug) << "DCA selection after vertex fit applied.";
12391241

12401242
// daughter DCAs to vertex
12411243
if (kfProtonDaughters[i].GetDistanceFromVertexXY(kf3bodyVtxCandidates[i]) >= kfparticleConfigurations.maxDcaXYSVDau || kfPionDaughters[i].GetDistanceFromVertexXY(kf3bodyVtxCandidates[i]) >= kfparticleConfigurations.maxDcaXYSVDau || kfDeuteronDaughters[i].GetDistanceFromVertexXY(kf3bodyVtxCandidates[i]) >= kfparticleConfigurations.maxDcaXYSVDau) {
12421244
continue;
12431245
}
1244-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxDcaDauVtx);
1246+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxDcaDauVtx);
12451247
LOG(debug) << "DCA to vertex selection after vertex fit applied.";
12461248

12471249
// daughter pT
12481250
if (kfProtonDaughters[i].GetPt() < kfparticleConfigurations.minPtProton || kfProtonDaughters[i].GetPt() > kfparticleConfigurations.maxPtProton || kfPionDaughters[i].GetPt() < kfparticleConfigurations.minPtPion || kfPionDaughters[i].GetPt() > kfparticleConfigurations.maxPtPion || kfDeuteronDaughters[i].GetPt() < kfparticleConfigurations.minPtDeuteron || kfDeuteronDaughters[i].GetPt() > kfparticleConfigurations.maxPtDeuteron) {
12491251
continue;
12501252
}
1251-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxDauPt);
1253+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxDauPt);
12521254
LOG(debug) << "Daughter pT selection applied.";
12531255

12541256
// -------- STEP 5: candidate selection and constraint after geometrical vertex fit --------
@@ -1257,40 +1259,40 @@ struct decay3bodyBuilder {
12571259
if (std::abs(rapHt) > kfparticleConfigurations.maxRapidityHt) {
12581260
continue;
12591261
}
1260-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxRap);
1262+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxRap);
12611263

12621264
// Pt selection
12631265
if (kf3bodyVtxCandidates[i].GetPt() <= kfparticleConfigurations.minPtHt || kf3bodyVtxCandidates[i].GetPt() >= kfparticleConfigurations.maxPtHt) {
12641266
continue;
12651267
}
1266-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxPt);
1268+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxPt);
12671269

12681270
// Mass window
12691271
float massHt, sigmaMassHt;
12701272
kf3bodyVtxCandidates[i].GetMass(massHt, sigmaMassHt);
12711273
if (massHt <= kfparticleConfigurations.minMassHt || massHt >= kfparticleConfigurations.maxMassHt) {
12721274
continue;
12731275
}
1274-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxMass);
1276+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxMass);
12751277

12761278
// cos(PA) to PV
12771279
if (std::abs(cpaFromKF(kf3bodyVtxCandidates[i], kfpv)) <= kfparticleConfigurations.minCosPA) {
12781280
continue;
12791281
}
1280-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxCosPA);
1282+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxCosPA);
12811283

12821284
// cos(PA) xy to PV
12831285
if (std::abs(cpaXYFromKF(kf3bodyVtxCandidates[i], kfpv)) <= kfparticleConfigurations.minCosPAxy) {
12841286
continue;
12851287
}
1286-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxCosPAXY);
1288+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxCosPAXY);
12871289

12881290
// chi2 geometrical
12891291
float chi2geoNDF = kf3bodyVtxCandidates[i].GetChi2() / kf3bodyVtxCandidates[i].GetNDF();
12901292
if (chi2geoNDF >= kfparticleConfigurations.maxChi2geo) {
12911293
continue;
12921294
}
1293-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxChi2geo);
1295+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxChi2geo);
12941296
LOG(debug) << "Basic selections after vertex fit done.";
12951297

12961298
// ctau before topo constraint
@@ -1306,14 +1308,14 @@ struct decay3bodyBuilder {
13061308
LOG(error) << "Exception caught KFParticle process call: Topological constraint failed";
13071309
continue;
13081310
}
1309-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxTopoConstr); // to check if topo constraint fails
1311+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxTopoConstr); // to check if topo constraint fails
13101312
// get topological chi2
13111313
float chi2topoNDF = KFHtPV.GetChi2() / KFHtPV.GetNDF();
13121314
KFHtPV.TransportToDecayVertex();
13131315
if (kfparticleConfigurations.applyTopoSel && chi2topoNDF >= kfparticleConfigurations.maxChi2topo) {
13141316
continue;
13151317
}
1316-
registry.fill(HIST("Counter/hVtx3BodyCounterKFParticle"), kKfVtxChi2topo);
1318+
registry.fill(HIST("Counters/hVtx3BodyCounterKFParticle"), kKfVtxChi2topo);
13171319

13181320
// additional selections from SVertexer for event mixing
13191321
float radius3body = sqrt(kf3bodyVtxCandidates[i].GetX() * kf3bodyVtxCandidates[i].GetX() + kf3bodyVtxCandidates[i].GetY() * kf3bodyVtxCandidates[i].GetY());
@@ -1432,7 +1434,7 @@ struct decay3bodyBuilder {
14321434
} // end loop over rotations
14331435

14341436
// fill event counter hist (has selected candidate) --> only filled once per vertex even in rotation case
1435-
registry.fill(HIST("Counter/hEventCounterKFParticle"), 3.5);
1437+
registry.fill(HIST("Counters/hEventCounterKFParticle"), 3.5);
14361438
} // end buildVtx3BodyDataTableKFParticle
14371439

14381440
//------------------------------------------------------------------
@@ -1588,23 +1590,23 @@ struct decay3bodyBuilder {
15881590
if (kfparticleConfigurations.cfgSkimmedProcessing) {
15891591
isZorroSelected = zorro.isSelected(collision.template bc_as<aod::BCsWithTimestamps>().globalBC());
15901592
if (isZorroSelected) {
1591-
registry.fill(HIST("Counter/hEventCounterZorro"), 0.);
1593+
registry.fill(HIST("Counters/hEventCounterZorro"), 0.);
15921594
}
15931595
}
15941596

15951597
// event selection
1596-
registry.fill(HIST("Counter/hEventCounterKFParticle"), 0.5);
1598+
registry.fill(HIST("Counters/hEventCounterKFParticle"), 0.5);
15971599
if (kfparticleConfigurations.doSel8selection && !collision.sel8()) {
15981600
continue;
15991601
}
1600-
registry.fill(HIST("Counter/hEventCounterKFParticle"), 1.5);
1602+
registry.fill(HIST("Counters/hEventCounterKFParticle"), 1.5);
16011603
if (kfparticleConfigurations.doPosZselection && abs(collision.posZ()) > 10.f) {
16021604
continue;
16031605
}
1604-
registry.fill(HIST("Counter/hEventCounterKFParticle"), 2.5);
1606+
registry.fill(HIST("Counters/hEventCounterKFParticle"), 2.5);
16051607

16061608
if (isZorroSelected) {
1607-
registry.fill(HIST("hCounter/EventCounterZorro"), 1.);
1609+
registry.fill(HIST("Counters/hEventCounterZorro"), 1.);
16081610
}
16091611

16101612
// slice Decay3Body table by collision
@@ -1687,15 +1689,15 @@ struct decay3bodyBuilder {
16871689
// same process as processRun3withKFParticle but with nRotations set according to configurable
16881690
for (const auto& collision : collisions) {
16891691
// event selection
1690-
registry.fill(HIST("Counter/hEventCounterKFParticle"), 0.5);
1692+
registry.fill(HIST("Counters/hEventCounterKFParticle"), 0.5);
16911693
if (kfparticleConfigurations.doSel8selection && !collision.sel8()) {
16921694
continue;
16931695
}
1694-
registry.fill(HIST("Counter/hEventCounterKFParticle"), 1.5);
1696+
registry.fill(HIST("Counters/hEventCounterKFParticle"), 1.5);
16951697
if (kfparticleConfigurations.doPosZselection && abs(collision.posZ()) > 10.f) {
16961698
continue;
16971699
}
1698-
registry.fill(HIST("Counter/hEventCounterKFParticle"), 2.5);
1700+
registry.fill(HIST("Counters/hEventCounterKFParticle"), 2.5);
16991701

17001702
auto bc = collision.bc_as<aod::BCsWithTimestamps>();
17011703
initCCDB(bc);

0 commit comments

Comments
 (0)