Skip to content

Commit 41a08aa

Browse files
author
Sawan Sawan
committed
added particle f2(1270) in MC. Added bin labels for event histogram
1 parent ed009f1 commit 41a08aa

File tree

1 file changed

+39
-27
lines changed

1 file changed

+39
-27
lines changed

‎PWGLF/Tasks/Resonances/higherMassResonances.cxx‎

Lines changed: 39 additions & 27 deletions
Original file line numberDiff line numberDiff line change
@@ -151,8 +151,8 @@ struct HigherMassResonances {
151151
Configurable<bool> iscTVXEvsel{"iscTVXEvsel", true, "Triggger selection"};
152152
Configurable<bool> isavoidsplitrackMC{"isavoidsplitrackMC", false, "avoid split track in MC"};
153153
Configurable<bool> isapplyRapidityMC{"isapplyRapidityMC", true, "Apply rapidity cut on generated and reconstructed particles"};
154-
Configurable<int> selectMCparticles{"selectMCparticles", 1, "0: f0(1710), 1: f2(1525), 2: a2(1320), 3: f0(1370), 4: f0(1500)"};
155-
std::vector<int> pdgCodes = {10331, 335, 115, 10221, 9030221};
154+
Configurable<int> selectMCparticles{"selectMCparticles", 1, "0: f0(1710), 1: f2(1525), 2: a2(1320), 3: f0(1370), 4: f0(1500), 5: f2(1270)"};
155+
std::vector<int> pdgCodes = {10331, 335, 115, 10221, 9030221, 225};
156156

157157
// output THnSparses
158158
Configurable<bool> activateHelicityFrame{"activateHelicityFrame", false, "Activate the THnSparse with cosThStar w.r.t. helicity axis"};
@@ -245,7 +245,7 @@ struct HigherMassResonances {
245245
rEventSelection.add("htrackscheck_v0", "htrackscheck_v0", kTH1I, {{15, 0, 15}});
246246
rEventSelection.add("htrackscheck_v0_daughters", "htrackscheck_v0_daughters", kTH1I, {{15, 0, 15}});
247247
hMChists.add("events_check", "No. of events in the generated MC", kTH1I, {{20, 0, 20}});
248-
hMChists.add("events_checkrec", "No. of events in the reconstructed MC", kTH1I, {{25, 0, 25}});
248+
hMChists.add("events_checkrec", "No. of events in the reconstructed MC", kTH1I, {{20, 0, 20}});
249249

250250
rEventSelection.add("hEventCut", "No. of event after cuts", kTH1I, {{20, 0, 20}});
251251
std::shared_ptr<TH1> hCutFlow = rEventSelection.get<TH1>(HIST("hEventCut"));
@@ -292,12 +292,23 @@ struct HigherMassResonances {
292292
hv0labelmcrec->GetXaxis()->SetBinLabel(2, "V0Daughter Sel.");
293293
hv0labelmcrec->GetXaxis()->SetBinLabel(3, "V0 Sel.");
294294
hv0labelmcrec->GetXaxis()->SetBinLabel(4, "V0 PDG");
295-
hv0labelmcrec->GetXaxis()->SetBinLabel(5, "Mother PDG");
296-
hv0labelmcrec->GetXaxis()->SetBinLabel(6, "Same Mother");
297-
hv0labelmcrec->GetXaxis()->SetBinLabel(7, "Split Track");
298-
hv0labelmcrec->GetXaxis()->SetBinLabel(8, "Global Index");
299-
hv0labelmcrec->GetXaxis()->SetBinLabel(9, "Generator");
300-
hv0labelmcrec->GetXaxis()->SetBinLabel(10, "Rapidity");
295+
hv0labelmcrec->GetXaxis()->SetBinLabel(5, "All Mothers");
296+
hv0labelmcrec->GetXaxis()->SetBinLabel(6, "Mother PDG");
297+
hv0labelmcrec->GetXaxis()->SetBinLabel(7, "Same Mother");
298+
hv0labelmcrec->GetXaxis()->SetBinLabel(8, "Split Track");
299+
hv0labelmcrec->GetXaxis()->SetBinLabel(9, "Global Index");
300+
hv0labelmcrec->GetXaxis()->SetBinLabel(10, "Generator");
301+
hv0labelmcrec->GetXaxis()->SetBinLabel(11, "Rapidity");
302+
303+
std::shared_ptr<TH1> hv0labelmcgen = hMChists.get<TH1>(HIST("events_check"));
304+
hv0labelmcgen->GetXaxis()->SetBinLabel(1, "All Events");
305+
hv0labelmcgen->GetXaxis()->SetBinLabel(2, "Event Sel.");
306+
hv0labelmcgen->GetXaxis()->SetBinLabel(3, "Event reconstructed");
307+
hv0labelmcgen->GetXaxis()->SetBinLabel(4, "PDG check");
308+
hv0labelmcgen->GetXaxis()->SetBinLabel(5, "Rapidity");
309+
hv0labelmcgen->GetXaxis()->SetBinLabel(6, "Daughters2");
310+
hv0labelmcgen->GetXaxis()->SetBinLabel(7, "PhysicalPrimary");
311+
hv0labelmcgen->GetXaxis()->SetBinLabel(8, "Daughters K0s");
301312
}
302313

303314
hglue.add("h3glueInvMassDS", "h3glueInvMassDS", kTHnSparseF, {multiplicityAxis, ptAxis, glueballMassAxis, thnAxisPOL, thnAxisPhi}, true);
@@ -367,8 +378,8 @@ struct HigherMassResonances {
367378
if (fillHist)
368379
rEventSelection.fill(HIST("hEventCut"), 1);
369380

370-
if (!collision.sel8())
371-
return false;
381+
// if (!collision.sel8())
382+
// return false;
372383
if (fillHist)
373384
rEventSelection.fill(HIST("hEventCut"), 2);
374385

@@ -405,8 +416,8 @@ struct HigherMassResonances {
405416
if (fillHist)
406417
rEventSelection.fill(HIST("hEventCut"), 8);
407418

408-
if (config.isTriggerTVX && !collision.selection_bit(aod::evsel::kIsTriggerTVX))
409-
return false;
419+
// if (config.isTriggerTVX && !collision.selection_bit(aod::evsel::kIsTriggerTVX))
420+
// return false;
410421
if (fillHist)
411422
rEventSelection.fill(HIST("hEventCut"), 9);
412423

@@ -1147,42 +1158,42 @@ struct HigherMassResonances {
11471158
selectedEvents[nevts++] = collision.mcCollision_as<aod::McCollisions>().globalIndex();
11481159
}
11491160
selectedEvents.resize(nevts);
1150-
hMChists.fill(HIST("events_check"), 3.5);
1161+
hMChists.fill(HIST("events_check"), 1.5);
11511162
const auto evtReconstructedAndSelected = std::find(selectedEvents.begin(), selectedEvents.end(), mcCollision.globalIndex()) != selectedEvents.end();
11521163

11531164
if (!config.isallGenCollisions && !evtReconstructedAndSelected) { // Check that the event is reconstructed and that the reconstructed events pass the selection
11541165
return;
11551166
}
1156-
hMChists.fill(HIST("events_check"), 4.5);
1167+
hMChists.fill(HIST("events_check"), 2.5);
11571168
for (const auto& mcParticle : mcParticles) {
11581169

11591170
if (std::abs(mcParticle.pdgCode()) != config.pdgCodes[config.selectMCparticles]) // f2(1525), f0(1710)
11601171
{
11611172
continue;
11621173
}
1163-
hMChists.fill(HIST("events_check"), 5.5);
1174+
hMChists.fill(HIST("events_check"), 3.5);
11641175

11651176
if (config.isapplyRapidityMC && std::abs(mcParticle.y()) >= config.rapidityMotherData) {
11661177
continue;
11671178
}
1168-
hMChists.fill(HIST("events_check"), 6.5);
1179+
hMChists.fill(HIST("events_check"), 4.5);
11691180

11701181
auto kDaughters = mcParticle.daughters_as<aod::McParticles>();
11711182
if (kDaughters.size() != config.noOfDaughters) {
11721183
continue;
11731184
}
1174-
hMChists.fill(HIST("events_check"), 7.5);
1185+
hMChists.fill(HIST("events_check"), 5.5);
11751186

11761187
for (const auto& kCurrentDaughter : kDaughters) {
11771188
// int daupdg = std::abs(kCurrentDaughter.pdgCode());
11781189

11791190
if (!kCurrentDaughter.isPhysicalPrimary()) {
11801191
continue;
11811192
}
1182-
hMChists.fill(HIST("events_check"), 8.5);
1193+
hMChists.fill(HIST("events_check"), 6.5);
11831194
if (std::abs(kCurrentDaughter.pdgCode()) == PDG_t::kK0Short) {
11841195
passKs.push_back(true);
1185-
hMChists.fill(HIST("events_check"), 9.5);
1196+
hMChists.fill(HIST("events_check"), 7.5);
11861197
if (passKs.size() == 1) {
11871198
daughter1 = ROOT::Math::PxPyPzMVector(kCurrentDaughter.px(), kCurrentDaughter.py(), kCurrentDaughter.pz(), o2::constants::physics::MassK0Short);
11881199
} else if (static_cast<int>(passKs.size()) == config.noOfDaughters) {
@@ -1323,45 +1334,46 @@ struct HigherMassResonances {
13231334
}
13241335

13251336
for (const auto& mothertrack2 : mctrackv02.mothers_as<aod::McParticles>()) {
1337+
hMChists.fill(HIST("events_checkrec"), 4.5);
13261338

13271339
if (mothertrack1.pdgCode() != config.pdgCodes[config.selectMCparticles]) {
13281340
continue;
13291341
}
1330-
hMChists.fill(HIST("events_checkrec"), 4.5);
1342+
hMChists.fill(HIST("events_checkrec"), 5.5);
13311343

13321344
if (mothertrack1.pdgCode() != mothertrack2.pdgCode()) {
13331345
continue;
13341346
}
1335-
hMChists.fill(HIST("events_checkrec"), 5.5);
1347+
hMChists.fill(HIST("events_checkrec"), 6.5);
13361348

13371349
gindex2.push_back(mothertrack2.globalIndex());
13381350
if (gindex2.size() > 1) {
13391351
if (std::find(gindex2.begin(), gindex2.end(), mothertrack2.globalIndex()) != gindex2.end()) {
13401352
continue;
13411353
}
13421354
}
1343-
hMChists.fill(HIST("events_checkrec"), 6.5);
1355+
hMChists.fill(HIST("events_checkrec"), 7.5);
13441356

13451357
if (mothertrack1.globalIndex() != mothertrack2.globalIndex()) {
13461358
continue;
13471359
}
1348-
hMChists.fill(HIST("events_checkrec"), 7.5);
1360+
hMChists.fill(HIST("events_checkrec"), 8.5);
13491361

13501362
if (!mothertrack1.producedByGenerator()) {
13511363
continue;
13521364
}
1353-
hMChists.fill(HIST("events_checkrec"), 8.5);
1365+
hMChists.fill(HIST("events_checkrec"), 9.5);
13541366

13551367
if (config.isapplyRapidityMC && std::abs(mothertrack1.y()) >= config.rapidityMotherData) {
13561368
continue;
13571369
}
1358-
hMChists.fill(HIST("events_checkrec"), 9.5);
1370+
hMChists.fill(HIST("events_checkrec"), 10.5);
13591371

13601372
// if (config.isavoidsplitrackMC && oldindex == mothertrack1.globalIndex()) {
13611373
// hMChists.fill(HIST("h1Recsplit"), mothertrack1.pt());
13621374
// continue;
13631375
// }
1364-
// hMChists.fill(HIST("events_checkrec"), 20.5);
1376+
// hMChists.fill(HIST("events_checkrec"), 11.5);
13651377
// oldindex = mothertrack1.globalIndex(); // split tracks is already handled using gindex1 and gindex2
13661378

13671379
daughter1 = ROOT::Math::PxPyPzMVector(v01.px(), v01.py(), v01.pz(), o2::constants::physics::MassK0Short);

0 commit comments

Comments
 (0)