Skip to content

Commit 1462322

Browse files
lorber98lorber98
committed
Add random cone histograms (phi randomized)
1 parent 4133d27 commit 1462322

1 file changed

Lines changed: 171 additions & 3 deletions

File tree

‎PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx‎

Lines changed: 171 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -154,6 +154,7 @@ struct StrangenessInJetsIons {
154154
Configurable<bool> useV0inJetRec{"useV0inJetRec", true, "Include V0s in jet reconstruction"};
155155
Configurable<bool> doThermalToyModel{"doThermalToyModel", false, "Use the thermal toy model to embed background particles to the jet finder input list"};
156156
Configurable<bool> doPlotRho{"doPlotRho", false, "Plot rho distributions computed with perpendicular cone and area median method."};
157+
Configurable<bool> doRandomConeSys{"doRandomConeSys", false, "Study particle contribution outside jets using random cones with same pseudorapidity and phi generated in [pi/3, 2*pi/3]."};
157158
Configurable<bool> saveChargedParticleMB{"saveChargedParticleMB", false, "Store charged particle information to build inclusive spectra."};
158159

159160
// Event selection
@@ -336,9 +337,10 @@ struct StrangenessInJetsIons {
336337
registryData.add("h2_centrality_deltaPt_RandomCone", "h2_centrality_deltaPt_RandomCone", HistType::kTH2F, {multAxis, deltaPtAxis});
337338
registryData.add("h2_centrality_rhoPerp", "h2_centrality_rhoPerp", HistType::kTH2F, {multAxis, rhoAxis});
338339

339-
registryData.add("rho_perp", "rho_perp", HistType::kTH2F, {multAxis, rhoAxis});
340-
registryData.add("rho_median", "rho_median", HistType::kTH2F, {multAxis, rhoAxis});
341-
340+
if (doPlotRho) {
341+
registryData.add("rho_perp", "rho_perp", HistType::kTH2F, {multAxis, rhoAxis});
342+
registryData.add("rho_median", "rho_median", HistType::kTH2F, {multAxis, rhoAxis});
343+
}
342344
// Armenteros-Podolanski plot
343345
// registryQC.add("ArmenterosPreSel_DATA", "ArmenterosPreSel_DATA", HistType::kTH2F, {alphaArmAxis, qtarmAxis});
344346

@@ -350,6 +352,11 @@ struct StrangenessInJetsIons {
350352
registryData.add("AntiLambda_in_ue", "AntiLambda_in_ue", HistType::kTH3F, {multAxis, ptAxis, invMassLambdaAxis});
351353
registryData.add("K0s_in_jet", "K0s_in_jet", HistType::kTH3F, {multAxis, ptAxis, invMassK0sAxis});
352354
registryData.add("K0s_in_ue", "K0s_in_ue", HistType::kTH3F, {multAxis, ptAxis, invMassK0sAxis});
355+
if (doRandomConeSys) {
356+
registryData.add("Lambda_in_rc", "Lambda_in_rc", HistType::kTH3F, {multAxis, ptAxis, invMassLambdaAxis});
357+
registryData.add("AntiLambda_in_rc", "AntiLambda_in_rc", HistType::kTH3F, {multAxis, ptAxis, invMassLambdaAxis});
358+
registryData.add("K0s_in_rc", "K0s_in_rc", HistType::kTH3F, {multAxis, ptAxis, invMassK0sAxis});
359+
}
353360
}
354361
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
355362
registryData.add("XiPos_in_jet", "XiPos_in_jet", HistType::kTH3F, {multAxis, ptAxis, invMassXiAxis});
@@ -414,6 +421,11 @@ struct StrangenessInJetsIons {
414421
registryMC.add("Lambda_generated_ue", "Lambda_generated_ue", HistType::kTH2F, {multAxis, ptAxis});
415422
registryMC.add("AntiLambda_generated_jet", "AntiLambda_generated_jet", HistType::kTH2F, {multAxis, ptAxis});
416423
registryMC.add("AntiLambda_generated_ue", "AntiLambda_generated_ue", HistType::kTH2F, {multAxis, ptAxis});
424+
if (doRandomConeSys) {
425+
registryMC.add("Lambda_generated_rc", "Lambda_generated_rc", HistType::kTH2F, {multAxis, ptAxis});
426+
registryMC.add("AntiLambda_generated_rc", "AntiLambda_generated_rc", HistType::kTH2F, {multAxis, ptAxis});
427+
registryMC.add("K0s_generated_rc", "K0s_generated_rc", HistType::kTH2F, {multAxis, ptAxis});
428+
}
417429

418430
// --- Histograms for the full event (without jets)
419431
registryMC.add("K0s_generated_MB", "K0s_generated_MB", HistType::kTH2F, {multAxis, ptAxis});
@@ -521,6 +533,11 @@ struct StrangenessInJetsIons {
521533
registryMC.add("Lambda_reconstructed_ue", "Lambda_reconstructed_ue", HistType::kTH2F, {multAxis, ptAxis});
522534
registryMC.add("AntiLambda_reconstructed_jet", "AntiLambda_reconstructed_jet", HistType::kTH2F, {multAxis, ptAxis});
523535
registryMC.add("AntiLambda_reconstructed_ue", "AntiLambda_reconstructed_ue", HistType::kTH2F, {multAxis, ptAxis});
536+
if (doRandomConeSys) {
537+
registryMC.add("Lambda_reconstructed_rc", "Lambda_reconstructed_rc", HistType::kTH2F, {multAxis, ptAxis});
538+
registryMC.add("AntiLambda_reconstructed_rc", "AntiLambda_reconstructed_rc", HistType::kTH2F, {multAxis, ptAxis});
539+
registryMC.add("K0s_reconstructed_rc", "K0s_reconstructed_rc", HistType::kTH2F, {multAxis, ptAxis});
540+
}
524541

525542
// Histograms for secondary hadrons
526543
registryMC.add("K0s_reconstructed_jet_incl", "K0s_reconstructed_jet_incl", HistType::kTH2F, {multAxis, ptAxis});
@@ -529,6 +546,11 @@ struct StrangenessInJetsIons {
529546
registryMC.add("Lambda_reconstructed_ue_incl", "Lambda_reconstructed_ue_incl", HistType::kTH2F, {multAxis, ptAxis});
530547
registryMC.add("AntiLambda_reconstructed_jet_incl", "AntiLambda_reconstructed_jet_incl", HistType::kTH2F, {multAxis, ptAxis});
531548
registryMC.add("AntiLambda_reconstructed_ue_incl", "AntiLambda_reconstructed_ue_incl", HistType::kTH2F, {multAxis, ptAxis});
549+
if (doRandomConeSys) {
550+
registryMC.add("Lambda_reconstructed_rc_incl", "Lambda_reconstructed_rc_incl", HistType::kTH2F, {multAxis, ptAxis});
551+
registryMC.add("AntiLambda_reconstructed_rc_incl", "AntiLambda_reconstructed_rc_incl", HistType::kTH2F, {multAxis, ptAxis});
552+
registryMC.add("K0s_reconstructed_rc_incl", "K0s_reconstructed_rc_incl", HistType::kTH2F, {multAxis, ptAxis});
553+
}
532554

533555
// Histograms for generated particles in reconstructed events
534556
registryMC.add("K0s_gen_recoEvent_jet", "K0s_gen_recoEvent_jet", HistType::kTH2F, {multAxis, ptAxis});
@@ -537,6 +559,11 @@ struct StrangenessInJetsIons {
537559
registryMC.add("Lambda_gen_recoEvent_ue", "Lambda_gen_recoEvent_ue", HistType::kTH2F, {multAxis, ptAxis});
538560
registryMC.add("AntiLambda_gen_recoEvent_jet", "AntiLambda_gen_recoEvent_jet", HistType::kTH2F, {multAxis, ptAxis});
539561
registryMC.add("AntiLambda_gen_recoEvent_ue", "AntiLambda_gen_recoEvent_ue", HistType::kTH2F, {multAxis, ptAxis});
562+
if (doRandomConeSys) {
563+
registryMC.add("Lambda_gen_recoEvent_rc", "Lambda_gen_recoEvent_rc", HistType::kTH2F, {multAxis, ptAxis});
564+
registryMC.add("AntiLambda_gen_recoEvent_rc", "AntiLambda_gen_recoEvent_rc", HistType::kTH2F, {multAxis, ptAxis});
565+
registryMC.add("K0s_gen_recoEvent_rc", "K0s_gen_recoEvent_rc", HistType::kTH2F, {multAxis, ptAxis});
566+
}
540567

541568
// --- Histograms for the full event (without jets)
542569
registryMC.add("K0s_reconstructed_MB", "K0s_reconstructed_MB", HistType::kTH2F, {multAxis, ptAxis});
@@ -755,6 +782,23 @@ struct StrangenessInJetsIons {
755782
u2.SetXYZ(u2x, u2y, pz);
756783
}
757784

785+
// Compute two random directions in phi to vector p
786+
void getRandomConeDirections(const fastjet::PseudoJet& jet, TVector3& u1, TVector3& u2, double minDeltaPhi, double maxDeltaPhi)
787+
{
788+
// Generate random deltaPhi in specified range
789+
double randomDeltaPhi = fRng.Uniform(minDeltaPhi, maxDeltaPhi);
790+
791+
u1.SetPtEtaPhi(jet.pt(), jet.eta(), jet.phi() + randomDeltaPhi);
792+
u2.SetPtEtaPhi(jet.pt(), jet.eta(), jet.phi() - randomDeltaPhi);
793+
}
794+
795+
float getDeltaR(double eta1, double phi1, double etaAxis, double phiAxis)
796+
{
797+
float deltaEta = eta1 - etaAxis;
798+
float deltaPhi = getDeltaPhi(phi1, phiAxis);
799+
return std::sqrt(deltaEta * deltaEta + deltaPhi * deltaPhi);
800+
}
801+
758802
void computeRandomConeDeltaPt(const std::vector<fastjet::PseudoJet>& fjParticles,
759803
const std::vector<fastjet::PseudoJet>& jets,
760804
float multiplicity, double rhoPerp)
@@ -2222,6 +2266,8 @@ struct StrangenessInJetsIons {
22222266
std::vector<TVector3> selectedJet;
22232267
std::vector<TVector3> ue1;
22242268
std::vector<TVector3> ue2;
2269+
std::vector<TVector3> rcSys1; // randomn cone 1
2270+
std::vector<TVector3> rcSys2; // randomn cone 2
22252271
std::vector<double> jetPt;
22262272

22272273
// Loop over reconstructed jets
@@ -2246,6 +2292,13 @@ struct StrangenessInJetsIons {
22462292
continue;
22472293
}
22482294

2295+
TVector3 rcSysAxis1(0, 0, 0), rcSysAxis2(0, 0, 0);
2296+
if (doRandomConeSys) {
2297+
getRandomConeDirections(jet, rcSysAxis1, rcSysAxis2, TMath::Pi() / 3.0, 2.0 * TMath::Pi() / 3.0);
2298+
rcSys1.emplace_back(rcSysAxis1);
2299+
rcSys2.emplace_back(rcSysAxis2);
2300+
}
2301+
22492302
// Store jet and UE axes
22502303
selectedJet.emplace_back(jetAxis);
22512304
ue1.emplace_back(ueAxis1);
@@ -2287,6 +2340,13 @@ struct StrangenessInJetsIons {
22872340
const float deltaPhiUe2 = getDeltaPhi(v0dir.Phi(), ue2[i].Phi());
22882341
const float deltaRue2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
22892342

2343+
bool isInRC = false;
2344+
if (doRandomConeSys) {
2345+
const float deltaRrc1 = getDeltaR(v0dir.Eta(), v0dir.Phi(), rcSys1[i].Eta(), rcSys1[i].Phi());
2346+
const float deltaRrc2 = getDeltaR(v0dir.Eta(), v0dir.Phi(), rcSys2[i].Eta(), rcSys2[i].Phi());
2347+
isInRC = (deltaRrc1 < rJet) || (deltaRrc2 < rJet);
2348+
}
2349+
22902350
// Fill Armenteros-Podolanski TH2
22912351
// registryQC.fill(HIST("ArmenterosPreSel_DATA"), v0.alpha(), v0.qtarm());
22922352

@@ -2298,6 +2358,9 @@ struct StrangenessInJetsIons {
22982358
if (deltaRue1 < rJet || deltaRue2 < rJet) {
22992359
registryData.fill(HIST("K0s_in_ue"), multiplicity, v0.pt(), v0.mK0Short());
23002360
}
2361+
if (doRandomConeSys && isInRC) {
2362+
registryData.fill(HIST("K0s_in_rc"), multiplicity, v0.pt(), v0.mK0Short());
2363+
}
23012364
}
23022365
// Lambda
23032366
if (passedLambdaSelection(v0, pos, neg, vtxPos)) {
@@ -2307,6 +2370,9 @@ struct StrangenessInJetsIons {
23072370
if (deltaRue1 < rJet || deltaRue2 < rJet) {
23082371
registryData.fill(HIST("Lambda_in_ue"), multiplicity, v0.pt(), v0.mLambda());
23092372
}
2373+
if (doRandomConeSys && isInRC) {
2374+
registryData.fill(HIST("Lambda_in_rc"), multiplicity, v0.pt(), v0.mLambda());
2375+
}
23102376
}
23112377
// AntiLambda
23122378
if (passedAntiLambdaSelection(v0, pos, neg, vtxPos)) {
@@ -2316,6 +2382,9 @@ struct StrangenessInJetsIons {
23162382
if (deltaRue1 < rJet || deltaRue2 < rJet) {
23172383
registryData.fill(HIST("AntiLambda_in_ue"), multiplicity, v0.pt(), v0.mAntiLambda());
23182384
}
2385+
if (doRandomConeSys && isInRC) {
2386+
registryData.fill(HIST("AntiLambda_in_rc"), multiplicity, v0.pt(), v0.mAntiLambda());
2387+
}
23192388
}
23202389
}
23212390
}
@@ -2600,6 +2669,11 @@ struct StrangenessInJetsIons {
26002669
continue;
26012670
}
26022671

2672+
TVector3 rcSysAxis1(0, 0, 0), rcSysAxis2(0, 0, 0);
2673+
if (doRandomConeSys) {
2674+
getRandomConeDirections(jet, rcSysAxis1, rcSysAxis2, TMath::Pi() / 3.0, 2.0 * TMath::Pi() / 3.0);
2675+
}
2676+
26032677
// Loop over strange hadrons
26042678
int index = -1;
26052679
for (const auto& hadron : strHadronMomentum) {
@@ -2617,6 +2691,13 @@ struct StrangenessInJetsIons {
26172691
double deltaPhiUe2 = getDeltaPhi(hadron.Phi(), ueAxis2.Phi());
26182692
double deltaRUe2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
26192693

2694+
bool isInRC = false;
2695+
if (doRandomConeSys) {
2696+
const float deltaRrc1 = getDeltaR(hadron.Eta(), hadron.Phi(), rcSysAxis1.Eta(), rcSysAxis1.Phi());
2697+
const float deltaRrc2 = getDeltaR(hadron.Eta(), hadron.Phi(), rcSysAxis2.Eta(), rcSysAxis2.Phi());
2698+
isInRC = (deltaRrc1 < rJet) || (deltaRrc2 < rJet);
2699+
}
2700+
26202701
// Select particles inside jet
26212702
if (deltaRJet < coneRadius) {
26222703
switch (pdg[index]) {
@@ -2762,6 +2843,29 @@ struct StrangenessInJetsIons {
27622843
break;
27632844
}
27642845
}
2846+
2847+
// Select particles inside random cones (RC)
2848+
if (doRandomConeSys && isInRC) {
2849+
switch (pdg[index]) {
2850+
case kK0Short:
2851+
if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) {
2852+
registryMC.fill(HIST("K0s_generated_rc"), genMultiplicity, hadron.Pt());
2853+
}
2854+
break;
2855+
case kLambda0:
2856+
if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) {
2857+
registryMC.fill(HIST("Lambda_generated_rc"), genMultiplicity, hadron.Pt());
2858+
}
2859+
break;
2860+
case kLambda0Bar:
2861+
if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) {
2862+
registryMC.fill(HIST("AntiLambda_generated_rc"), genMultiplicity, hadron.Pt());
2863+
}
2864+
break;
2865+
default:
2866+
break;
2867+
}
2868+
}
27652869
}
27662870
}
27672871
// Fill jet counter
@@ -2782,6 +2886,8 @@ struct StrangenessInJetsIons {
27822886
std::vector<TVector3> selectedJet;
27832887
std::vector<TVector3> ue1;
27842888
std::vector<TVector3> ue2;
2889+
std::vector<TVector3> rcSys1; // randomn cone 1
2890+
std::vector<TVector3> rcSys2; // randomn cone 2
27852891

27862892
// Jet and area definitions
27872893
fastjet::JetDefinition jetDef(fastjet::antikt_algorithm, rJet);
@@ -2804,6 +2910,8 @@ struct StrangenessInJetsIons {
28042910
selectedJet.clear();
28052911
ue1.clear();
28062912
ue2.clear();
2913+
rcSys1.clear();
2914+
rcSys2.clear();
28072915

28082916
// Fill event counter before any selection
28092917
registryMC.fill(HIST("number_of_events_mc_rec"), 0.5);
@@ -2907,6 +3015,13 @@ struct StrangenessInJetsIons {
29073015
continue;
29083016
}
29093017

3018+
TVector3 rcSysAxis1(0, 0, 0), rcSysAxis2(0, 0, 0);
3019+
if (doRandomConeSys) {
3020+
getRandomConeDirections(jet, rcSysAxis1, rcSysAxis2, TMath::Pi() / 3.0, 2.0 * TMath::Pi() / 3.0);
3021+
rcSys1.emplace_back(rcSysAxis1);
3022+
rcSys2.emplace_back(rcSysAxis2);
3023+
}
3024+
29103025
// Store selected jet and UE cone axes
29113026
selectedJet.emplace_back(jetAxis);
29123027
ue1.emplace_back(ueAxis1);
@@ -2950,6 +3065,13 @@ struct StrangenessInJetsIons {
29503065
const double deltaPhiUe2 = getDeltaPhi(momVec.Phi(), ue2[i].Phi());
29513066
const double deltaRUe2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
29523067

3068+
bool isInRC = false;
3069+
if (doRandomConeSys) {
3070+
const float deltaRrc1 = getDeltaR(momVec.Eta(), momVec.Phi(), rcSys1[i].Eta(), rcSys1[i].Phi());
3071+
const float deltaRrc2 = getDeltaR(momVec.Eta(), momVec.Phi(), rcSys2[i].Eta(), rcSys2[i].Phi());
3072+
isInRC = (deltaRrc1 < rJet) || (deltaRrc2 < rJet);
3073+
}
3074+
29533075
// Select particles inside jet
29543076
if (deltaRJet < rJet) {
29553077
switch (particle.pdgCode()) {
@@ -2993,6 +3115,27 @@ struct StrangenessInJetsIons {
29933115
break;
29943116
}
29953117
}
3118+
if (doRandomConeSys && isInRC) {
3119+
switch (particle.pdgCode()) {
3120+
case kK0Short:
3121+
if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) {
3122+
registryMC.fill(HIST("K0s_gen_recoEvent_rc"), multiplicity, momVec.Pt());
3123+
}
3124+
break;
3125+
case kLambda0:
3126+
if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) {
3127+
registryMC.fill(HIST("Lambda_gen_recoEvent_rc"), multiplicity, momVec.Pt());
3128+
}
3129+
break;
3130+
case kLambda0Bar:
3131+
if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) {
3132+
registryMC.fill(HIST("AntiLambda_gen_recoEvent_rc"), multiplicity, momVec.Pt());
3133+
}
3134+
break;
3135+
default:
3136+
break;
3137+
}
3138+
}
29963139
}
29973140
// ----------------------------------------
29983141

@@ -3047,6 +3190,13 @@ struct StrangenessInJetsIons {
30473190
double deltaPhiUe2 = getDeltaPhi(v0dir.Phi(), ue2[i].Phi());
30483191
double deltaRue2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
30493192

3193+
bool isInRC = false;
3194+
if (doRandomConeSys) {
3195+
const float deltaRrc1 = getDeltaR(v0dir.Eta(), v0dir.Phi(), rcSys1[i].Eta(), rcSys1[i].Phi());
3196+
const float deltaRrc2 = getDeltaR(v0dir.Eta(), v0dir.Phi(), rcSys2[i].Eta(), rcSys2[i].Phi());
3197+
isInRC = (deltaRrc1 < rJet) || (deltaRrc2 < rJet);
3198+
}
3199+
30503200
// Fill Armenteros-Podolanski TH2
30513201
// registryQC.fill(HIST("ArmenterosPreSel_REC"), v0.alpha(), v0.qtarm());
30523202

@@ -3058,6 +3208,9 @@ struct StrangenessInJetsIons {
30583208
if (deltaRue1 < rJet || deltaRue2 < rJet) {
30593209
registryMC.fill(HIST("K0s_reconstructed_ue"), multiplicity, v0.pt());
30603210
}
3211+
if (doRandomConeSys && isInRC) {
3212+
registryMC.fill(HIST("K0s_reconstructed_rc"), multiplicity, v0.pt());
3213+
}
30613214
}
30623215
// Lambda
30633216
if (passedLambdaSelection(v0, pos, neg, vtxPos) && pdgParent == kLambda0 && isPhysPrim) {
@@ -3067,6 +3220,9 @@ struct StrangenessInJetsIons {
30673220
if (deltaRue1 < rJet || deltaRue2 < rJet) {
30683221
registryMC.fill(HIST("Lambda_reconstructed_ue"), multiplicity, v0.pt());
30693222
}
3223+
if (doRandomConeSys && isInRC) {
3224+
registryMC.fill(HIST("Lambda_reconstructed_rc"), multiplicity, v0.pt());
3225+
}
30703226
}
30713227
// AntiLambda
30723228
if (passedAntiLambdaSelection(v0, pos, neg, vtxPos) && pdgParent == kLambda0Bar && isPhysPrim) {
@@ -3076,6 +3232,9 @@ struct StrangenessInJetsIons {
30763232
if (deltaRue1 < rJet || deltaRue2 < rJet) {
30773233
registryMC.fill(HIST("AntiLambda_reconstructed_ue"), multiplicity, v0.pt());
30783234
}
3235+
if (doRandomConeSys && isInRC) {
3236+
registryMC.fill(HIST("AntiLambda_reconstructed_rc"), multiplicity, v0.pt());
3237+
}
30793238
}
30803239

30813240
// Fill inclusive spectra
@@ -3087,6 +3246,9 @@ struct StrangenessInJetsIons {
30873246
if (deltaRue1 < rJet || deltaRue2 < rJet) {
30883247
registryMC.fill(HIST("K0s_reconstructed_ue_incl"), multiplicity, v0.pt());
30893248
}
3249+
if (doRandomConeSys && isInRC) {
3250+
registryMC.fill(HIST("K0s_reconstructed_rc_incl"), multiplicity, v0.pt());
3251+
}
30903252
}
30913253
// Lambda
30923254
if (passedLambdaSelection(v0, pos, neg, vtxPos) && pdgParent == kLambda0) {
@@ -3096,6 +3258,9 @@ struct StrangenessInJetsIons {
30963258
if (deltaRue1 < rJet || deltaRue2 < rJet) {
30973259
registryMC.fill(HIST("Lambda_reconstructed_ue_incl"), multiplicity, v0.pt());
30983260
}
3261+
if (doRandomConeSys && isInRC) {
3262+
registryMC.fill(HIST("Lambda_reconstructed_rc_incl"), multiplicity, v0.pt());
3263+
}
30993264
}
31003265
// AntiLambda
31013266
if (passedAntiLambdaSelection(v0, pos, neg, vtxPos) && pdgParent == kLambda0Bar) {
@@ -3105,6 +3270,9 @@ struct StrangenessInJetsIons {
31053270
if (deltaRue1 < rJet || deltaRue2 < rJet) {
31063271
registryMC.fill(HIST("AntiLambda_reconstructed_ue_incl"), multiplicity, v0.pt());
31073272
}
3273+
if (doRandomConeSys && isInRC) {
3274+
registryMC.fill(HIST("AntiLambda_reconstructed_rc_incl"), multiplicity, v0.pt());
3275+
}
31083276
}
31093277
}
31103278
}

0 commit comments

Comments
 (0)