Skip to content

Commit 7bb2b64

Browse files
lorber98lorber98
committed
Add cascades random cone histograms (MC)
1 parent c8f93c1 commit 7bb2b64

1 file changed

Lines changed: 78 additions & 0 deletions

File tree

‎PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx‎

Lines changed: 78 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -468,6 +468,12 @@ struct StrangenessInJetsIons {
468468
registryMC.add("OmegaPos_generated_ue", "OmegaPos_generated_ue", HistType::kTH2F, {multAxis, ptAxis});
469469
registryMC.add("OmegaNeg_generated_jet", "OmegaNeg_generated_jet", HistType::kTH2F, {multAxis, ptAxis});
470470
registryMC.add("OmegaNeg_generated_ue", "OmegaNeg_generated_ue", HistType::kTH2F, {multAxis, ptAxis});
471+
if (doRandomConeSys) {
472+
registryMC.add("XiPos_generated_rc", "XiPos_generated_rc", HistType::kTH2F, {multAxis, ptAxis});
473+
registryMC.add("XiNeg_generated_rc", "XiNeg_generated_rc", HistType::kTH2F, {multAxis, ptAxis});
474+
registryMC.add("OmegaPos_generated_rc", "OmegaPos_generated_rc", HistType::kTH2F, {multAxis, ptAxis});
475+
registryMC.add("OmegaNeg_generated_rc", "OmegaNeg_generated_rc", HistType::kTH2F, {multAxis, ptAxis});
476+
}
471477

472478
// Histograms for the full event (without jets)
473479
registryMC.add("OmegaPos_generated_MB", "OmegaPos_generated_MB", HistType::kTH2F, {multAxis, ptAxis});
@@ -607,6 +613,19 @@ struct StrangenessInJetsIons {
607613
registryMC.add("OmegaPos_reconstructed_ue", "OmegaPos_reconstructed_ue", HistType::kTH2F, {multAxis, ptAxis});
608614
registryMC.add("OmegaNeg_reconstructed_jet", "OmegaNeg_reconstructed_jet", HistType::kTH2F, {multAxis, ptAxis});
609615
registryMC.add("OmegaNeg_reconstructed_ue", "OmegaNeg_reconstructed_ue", HistType::kTH2F, {multAxis, ptAxis});
616+
if (doRandomConeSys) {
617+
registryMC.add("XiPos_reconstructed_rc", "XiPos_reconstructed_rc", HistType::kTH2F, {multAxis, ptAxis});
618+
registryMC.add("XiNeg_reconstructed_rc", "XiNeg_reconstructed_rc", HistType::kTH2F, {multAxis, ptAxis});
619+
registryMC.add("OmegaPos_reconstructed_rc", "OmegaPos_reconstructed_rc", HistType::kTH2F, {multAxis, ptAxis});
620+
registryMC.add("OmegaNeg_reconstructed_rc", "OmegaNeg_reconstructed_rc", HistType::kTH2F, {multAxis, ptAxis});
621+
}
622+
623+
if (doRandomConeSys) {
624+
registryMC.add("XiPos_gen_recoEvent_rc", "XiPos_gen_recoEvent_rc", HistType::kTH2F, {multAxis, ptAxis});
625+
registryMC.add("XiNeg_gen_recoEvent_rc", "XiNeg_gen_recoEvent_rc", HistType::kTH2F, {multAxis, ptAxis});
626+
registryMC.add("OmegaPos_gen_recoEvent_rc", "OmegaPos_gen_recoEvent_rc", HistType::kTH2F, {multAxis, ptAxis});
627+
registryMC.add("OmegaNeg_gen_recoEvent_rc", "OmegaNeg_gen_recoEvent_rc", HistType::kTH2F, {multAxis, ptAxis});
628+
}
610629

611630
// Histograms for the full event (without jets)
612631
registryMC.add("XiPos_reconstructed_MB", "XiPos_reconstructed_MB", HistType::kTH2F, {multAxis, ptAxis});
@@ -2971,6 +2990,26 @@ struct StrangenessInJetsIons {
29712990
registryMC.fill(HIST("AntiLambda_generated_rc"), genMultiplicity, hadron.Pt());
29722991
}
29732992
break;
2993+
case kXiMinus:
2994+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
2995+
registryMC.fill(HIST("XiNeg_generated_rc"), genMultiplicity, hadron.Pt());
2996+
}
2997+
break;
2998+
case kXiPlusBar:
2999+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
3000+
registryMC.fill(HIST("XiPos_generated_rc"), genMultiplicity, hadron.Pt());
3001+
}
3002+
break;
3003+
case kOmegaMinus:
3004+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
3005+
registryMC.fill(HIST("OmegaNeg_generated_rc"), genMultiplicity, hadron.Pt());
3006+
}
3007+
break;
3008+
case kOmegaPlusBar:
3009+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
3010+
registryMC.fill(HIST("OmegaPos_generated_rc"), genMultiplicity, hadron.Pt());
3011+
}
3012+
break;
29743013
default:
29753014
break;
29763015
}
@@ -3244,6 +3283,26 @@ struct StrangenessInJetsIons {
32443283
registryMC.fill(HIST("AntiLambda_gen_recoEvent_rc"), multiplicity, momVec.Pt());
32453284
}
32463285
break;
3286+
case kXiMinus:
3287+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
3288+
registryMC.fill(HIST("XiNeg_gen_recoEvent_rc"), multiplicity, momVec.Pt());
3289+
}
3290+
break;
3291+
case kXiPlusBar:
3292+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
3293+
registryMC.fill(HIST("XiPos_gen_recoEvent_rc"), multiplicity, momVec.Pt());
3294+
}
3295+
break;
3296+
case kOmegaMinus:
3297+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
3298+
registryMC.fill(HIST("OmegaNeg_gen_recoEvent_rc"), multiplicity, momVec.Pt());
3299+
}
3300+
break;
3301+
case kOmegaPlusBar:
3302+
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
3303+
registryMC.fill(HIST("OmegaPos_gen_recoEvent_rc"), multiplicity, momVec.Pt());
3304+
}
3305+
break;
32473306
default:
32483307
break;
32493308
}
@@ -3437,6 +3496,13 @@ struct StrangenessInJetsIons {
34373496
double deltaPhiUe2 = getDeltaPhi(cascadeDir.Phi(), ue2[i].Phi());
34383497
double deltaRue2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
34393498

3499+
bool isInRC = false;
3500+
if (doRandomConeSys) {
3501+
const float deltaRrc1 = getDeltaR(cascadeDir.Eta(), cascadeDir.Phi(), rcSys1[i].Eta(), rcSys1[i].Phi());
3502+
const float deltaRrc2 = getDeltaR(cascadeDir.Eta(), cascadeDir.Phi(), rcSys2[i].Eta(), rcSys2[i].Phi());
3503+
isInRC = (deltaRrc1 < rJet) || (deltaRrc2 < rJet);
3504+
}
3505+
34403506
// Xi+
34413507
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() > 0 && pdgParent == kXiPlusBar) {
34423508
if (deltaRjet < rJet) {
@@ -3445,6 +3511,9 @@ struct StrangenessInJetsIons {
34453511
if (deltaRue1 < rJet || deltaRue2 < rJet) {
34463512
registryMC.fill(HIST("XiPos_reconstructed_ue"), multiplicity, casc.pt());
34473513
}
3514+
if (doRandomConeSys && isInRC) {
3515+
registryMC.fill(HIST("XiPos_reconstructed_rc"), multiplicity, casc.pt());
3516+
}
34483517
}
34493518
// Xi-
34503519
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() < 0 && pdgParent == kXiMinus) {
@@ -3454,6 +3523,9 @@ struct StrangenessInJetsIons {
34543523
if (deltaRue1 < rJet || deltaRue2 < rJet) {
34553524
registryMC.fill(HIST("XiNeg_reconstructed_ue"), multiplicity, casc.pt());
34563525
}
3526+
if (doRandomConeSys && isInRC) {
3527+
registryMC.fill(HIST("XiNeg_reconstructed_rc"), multiplicity, casc.pt());
3528+
}
34573529
}
34583530
// Omega+
34593531
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() > 0 && pdgParent == kOmegaPlusBar) {
@@ -3463,6 +3535,9 @@ struct StrangenessInJetsIons {
34633535
if (deltaRue1 < rJet || deltaRue2 < rJet) {
34643536
registryMC.fill(HIST("OmegaPos_reconstructed_ue"), multiplicity, casc.pt());
34653537
}
3538+
if (doRandomConeSys && isInRC) {
3539+
registryMC.fill(HIST("OmegaPos_reconstructed_rc"), multiplicity, casc.pt());
3540+
}
34663541
}
34673542
// Omega-
34683543
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() < 0 && pdgParent == kOmegaMinus) {
@@ -3472,6 +3547,9 @@ struct StrangenessInJetsIons {
34723547
if (deltaRue1 < rJet || deltaRue2 < rJet) {
34733548
registryMC.fill(HIST("OmegaNeg_reconstructed_ue"), multiplicity, casc.pt());
34743549
}
3550+
if (doRandomConeSys && isInRC) {
3551+
registryMC.fill(HIST("OmegaNeg_reconstructed_rc"), multiplicity, casc.pt());
3552+
}
34753553
}
34763554
}
34773555
}

0 commit comments

Comments
 (0)