Skip to content

Commit c4fbe95

Browse files
committed
conditionally store number of tracks in the event
1 parent 180568f commit c4fbe95

File tree

1 file changed

+8
-3
lines changed

1 file changed

+8
-3
lines changed

PWGLF/TableProducer/Nuspex/ebyeMaker.cxx

Lines changed: 8 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -220,6 +220,7 @@ struct ebyeMaker {
220220
int mRunNumber;
221221
float d_bz;
222222
uint8_t nTrackletsColl;
223+
uint8_t nTracksColl;
223224
// o2::base::MatLayerCylSet* lut = nullptr;
224225

225226
Configurable<int> cfgMaterialCorrection{"cfgMaterialCorrection", static_cast<int>(o2::base::Propagator::MatCorrType::USEMatCorrNONE), "Type of material correction"};
@@ -244,7 +245,7 @@ struct ebyeMaker {
244245
Configurable<float> etaMaxV0dau{"etaMaxV0dau", 0.8f, "maximum eta V0 daughters"};
245246
Configurable<float> outerPIDMin{"outerPIDMin", -4.f, "minimum outer PID"};
246247

247-
Configurable<bool> fillOnlySignal{"fillOnlySignal", false, "fill histograms only for true signal candidates (MC)"};
248+
Configurable<bool> storeTracksNum{"storeTracksNum", false, "store the number of tracks instead of tracklets"};
248249
Configurable<std::string> genName{"genname", "", "Genearator name: HIJING, PYTHIA8, ... Default: \"\""};
249250

250251
Configurable<uint8_t> triggerCut{"triggerCut", 0x0, "trigger cut to select"};
@@ -663,6 +664,7 @@ struct ebyeMaker {
663664

664665
gpu::gpustd::array<float, 2> dcaInfo;
665666
uint8_t nTracklets[2]{0, 0};
667+
uint8_t nTracks{0};
666668
for (const auto& track : tracks) {
667669

668670
if (track.trackType() == 255 && std::abs(track.eta()) < 1.2) { // tracklet
@@ -692,6 +694,8 @@ struct ebyeMaker {
692694
}
693695
histos.fill(HIST("QA/tpcSignal"), track.tpcInnerParam(), track.tpcSignal());
694696

697+
nTracks++;
698+
695699
for (int iP{0}; iP < kNpart; ++iP) {
696700
if (trackPt < ptMin[iP] || trackPt > ptMax[iP]) {
697701
continue;
@@ -752,6 +756,7 @@ struct ebyeMaker {
752756
if (doprocessRun2 || doprocessMcRun2 || doprocessMiniRun2 || doprocessMiniMcRun2) {
753757
histos.fill(HIST("QA/nTrklCorrelation"), nTracklets[0], nTracklets[1]);
754758
nTrackletsColl = nTracklets[1];
759+
nTracksColl = nTracks;
755760
}
756761

757762
if (lambdaPtMax > lambdaPtMin) {
@@ -1234,7 +1239,7 @@ struct ebyeMaker {
12341239
if (triggerCut != 0x0 && (trigger & triggerCut) != triggerCut) {
12351240
continue;
12361241
}
1237-
miniCollTable(static_cast<int8_t>(collision.posZ() * 10), trigger, nTrackletsColl, cV0M);
1242+
miniCollTable(static_cast<int8_t>(collision.posZ() * 10), trigger, storeTracksNum ? nTracksColl : nTrackletsColl, cV0M);
12381243

12391244
for (auto& candidateTrack : candidateTracks[0]) { // protons
12401245
auto tk = tracks.rawIteratorAt(candidateTrack.globalIndex);
@@ -1426,7 +1431,7 @@ struct ebyeMaker {
14261431
fillMcEvent(collision, tracks, V0Table_thisCollision, cV0M, mcParticles, mcLab);
14271432
fillMcGen(mcParticles, mcLab, collision.mcCollisionId());
14281433

1429-
miniCollTable(static_cast<int8_t>(collision.posZ() * 10), 0x0, nTrackletsColl, cV0M);
1434+
miniCollTable(static_cast<int8_t>(collision.posZ() * 10), 0x0, storeTracksNum ? nTracksColl : nTrackletsColl, cV0M);
14301435

14311436
for (auto& candidateTrack : candidateTracks[0]) { // protons
14321437
int selMask = -1;

0 commit comments

Comments
 (0)