Skip to content

Commit 60c568a

Browse files
JStaaJStaa
committed
Changed a number of cut parameters to be configurable
1 parent 47c97e8 commit 60c568a

File tree

1 file changed

+52
-43
lines changed

1 file changed

+52
-43
lines changed

‎PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx‎

Lines changed: 52 additions & 43 deletions
Original file line numberDiff line numberDiff line change
@@ -53,11 +53,16 @@ struct ThreeParticleCorrelations {
5353
float nSigma0 = 0.0, nSigma1 = 1.0, nSigma2 = 2.0, nSigma4 = 4.0, nSigma5 = 5.0;
5454

5555
// V0 filter parameters
56-
float tpcNCrossedRowsMin = 70.0;
57-
float decayRMin = 1.2, ctauMax = 30.0;
58-
float cosPAMin = 0.995;
59-
float dcaProtonMin = 0.05, dcaPionMin = 0.2;
60-
int dcaV0DauMax = 1;
56+
struct : ConfigurableGroup {
57+
std::string prefix = "V0Selection";
58+
Configurable<float> tpcNCrossedRows{"tpcNCrossedRows", 70.0, "Minimum number of TPC crossed rows"};
59+
Configurable<float> decayR{"decayR", 1.2, "Minimum V0 decay radius (cm)"};
60+
Configurable<float> ctau{"ctau", 30.0, "Maximum V0 proper lifetime (cm)"};
61+
Configurable<float> cosPA{"cosPA", 0.995, "Minimum V0 cosine of pointing angle"};
62+
Configurable<float> dcaProton{"dcaProton", 0.05, "Minimum DCA of proton daughter (cm)"};
63+
Configurable<float> dcaPion{"dcaPion", 0.2, "Minimum DCA of pion daughter (cm)"};
64+
Configurable<float> dcaV0Dau{"dcaV0Dau", 1.0, "Maximum DCA between V0 daughters"};
65+
} v0SelGroup;
6166

6267
// Track filter parameters
6368
float pionPtMin = 0.3, pionPtMax = 2.3, kaonPtMin = 0.5, kaonPtMax = 2.3, protonPtMin = 0.6;
@@ -69,7 +74,8 @@ struct ThreeParticleCorrelations {
6974
float rMin = 0.8, rMax = 2.5;
7075

7176
// Lambda invariant mass fit
72-
double dGaussSigma = 0.0021;
77+
Configurable<float> invMassNSigma{"invMassNSigma", 4.0, "Number of standard deviations from the mean of the Lambda invariant mass peak"};
78+
double dGaussSigma = 0.0019;
7379

7480
// Histogram registry
7581
HistogramRegistry rMECorrRegistry{"MECorrRegistry", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};
@@ -137,10 +143,13 @@ struct ThreeParticleCorrelations {
137143
SameKindPair<MyFilteredMCGenCollisions, MyFilteredMCParticles, BinningTypeMC> pairMC{collBinningMC, 5, -1, &cache};
138144

139145
// Process configurables
140-
Configurable<int> confBfieldSwitch{"confBfieldSwitch", 0, "Switch for the detector magnetic field (1 if Pos, -1 if Neg, 0 if both)"};
141-
Configurable<bool> confRatioCorrectionSwitch{"confRatioCorrectionSwitch", false, "Switch for correcting the negative spectra back to the positive spectra"};
142-
Configurable<bool> confFakeV0Switch{"confFakeV0Switch", false, "Switch for the fakeV0Filter function"};
143-
Configurable<bool> confRDSwitch{"confRDSwitch", true, "Switch for the radialDistanceFilter function"};
146+
struct : ConfigurableGroup {
147+
std::string prefix = "processSwitchBoard";
148+
Configurable<int> confBfieldSwitch{"confBfieldSwitch", 0, "Switch for the detector magnetic field (1 if Pos, -1 if Neg, 0 if both)"};
149+
Configurable<bool> confRatioCorrectionSwitch{"confRatioCorrectionSwitch", false, "Switch for correcting the negative spectra back to the positive spectra"};
150+
Configurable<bool> confFakeV0Switch{"confFakeV0Switch", false, "Switch for the fakeV0Filter function"};
151+
Configurable<bool> confRDSwitch{"confRDSwitch", true, "Switch for the radialDistanceFilter function"};
152+
} switchGroup;
144153

145154
// Efficiency histograms
146155
TH3D** hEffPions = new TH3D*[2];
@@ -187,8 +196,8 @@ struct ThreeParticleCorrelations {
187196
rQARegistry.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(2, "kIsGoodZvtxFT0vsPV");
188197
rQARegistry.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(3, "kNoSameBunchPileup");
189198

190-
rQARegistry.add("hEventCentrality", "hEventCentrality", {HistType::kTH1D, {{centralityAxis}}});
191-
rQARegistry.add("hEventCentrality_MC", "hEventCentrality_MC", {HistType::kTH1D, {{centralityAxis}}});
199+
rQARegistry.add("hEventCentrality", "hEventCentrality", {HistType::kTH1D, {{fineCentralityAxis}}});
200+
rQARegistry.add("hEventCentrality_MC", "hEventCentrality_MC", {HistType::kTH1D, {{fineCentralityAxis}}});
192201
rQARegistry.add("hEventZvtx", "hEventZvtx", {HistType::kTH1D, {{zvtxAxis}}});
193202
rQARegistry.add("hEventBfield", "hEventBfield", {HistType::kTH1D, {{2, -1, 1}}});
194203
rQARegistry.add("hTrackPt", "hTrackPt", {HistType::kTH1D, {{100, 0, 4}}});
@@ -359,8 +368,8 @@ struct ThreeParticleCorrelations {
359368

360369
auto bc = collision.bc_as<aod::BCsWithTimestamps>();
361370
auto bField = getMagneticField(bc.timestamp());
362-
if (confBfieldSwitch != 0) {
363-
if (std::signbit(static_cast<double>(confBfieldSwitch)) != std::signbit(bField)) {
371+
if (switchGroup.confBfieldSwitch != 0) {
372+
if (std::signbit(static_cast<double>(switchGroup.confBfieldSwitch)) != std::signbit(bField)) {
364373
return;
365374
}
366375
}
@@ -433,7 +442,7 @@ struct ThreeParticleCorrelations {
433442
deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf);
434443
deltaEta = trigger.eta() - associate.eta();
435444

436-
if (candMass >= MassLambda0 - 4 * dGaussSigma && candMass <= MassLambda0 + 4 * dGaussSigma) {
445+
if (candMass >= MassLambda0 - invMassNSigma*dGaussSigma && candMass <= MassLambda0 + invMassNSigma*dGaussSigma) {
437446
if (assocPID[0] == pionID) { // Pions
438447
rSECorrRegistry.fill(HIST("hSameLambdaPion_SGNL"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, collision.centFT0C()) / (trackEff(hEffPions, associate, collision.centFT0C()) * v0Efficiency));
439448
} else if (assocPID[0] == kaonID) { // Kaons
@@ -442,28 +451,28 @@ struct ThreeParticleCorrelations {
442451
rSECorrRegistry.fill(HIST("hSameLambdaProton_SGNL"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, collision.centFT0C()) / (trackEff(hEffProtons, associate, collision.centFT0C()) * v0Efficiency));
443452
}
444453

445-
} else if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
454+
} else if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
446455
if (assocPID[0] == pionID) { // Pions
447456
rSECorrRegistry.fill(HIST("hSameLambdaPion_SB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, collision.centFT0C()) / (trackEff(hEffPions, associate, collision.centFT0C()) * v0Efficiency));
448-
if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass < MassLambda0 - 4 * dGaussSigma) {
457+
if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass < MassLambda0 - invMassNSigma*dGaussSigma) {
449458
rSECorrRegistry.fill(HIST("hSameLambdaPion_leftSB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, collision.centFT0C()) / (trackEff(hEffPions, associate, collision.centFT0C()) * v0Efficiency));
450-
} else if (candMass > MassLambda0 + 4 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
459+
} else if (candMass > MassLambda0 + invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
451460
rSECorrRegistry.fill(HIST("hSameLambdaPion_rightSB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, collision.centFT0C()) / (trackEff(hEffPions, associate, collision.centFT0C()) * v0Efficiency));
452461
}
453462

454463
} else if (assocPID[0] == kaonID) { // Kaons
455464
rSECorrRegistry.fill(HIST("hSameLambdaKaon_SB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionKaons, associate, collision.centFT0C()) / (trackEff(hEffKaons, associate, collision.centFT0C()) * v0Efficiency));
456-
if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass < MassLambda0 - 4 * dGaussSigma) {
465+
if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass < MassLambda0 - invMassNSigma*dGaussSigma) {
457466
rSECorrRegistry.fill(HIST("hSameLambdaKaon_leftSB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionKaons, associate, collision.centFT0C()) / (trackEff(hEffKaons, associate, collision.centFT0C()) * v0Efficiency));
458-
} else if (candMass > MassLambda0 + 4 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
467+
} else if (candMass > MassLambda0 + invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
459468
rSECorrRegistry.fill(HIST("hSameLambdaKaon_rightSB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionKaons, associate, collision.centFT0C()) / (trackEff(hEffKaons, associate, collision.centFT0C()) * v0Efficiency));
460469
}
461470

462471
} else if (assocPID[0] == protonID) { // Protons
463472
rSECorrRegistry.fill(HIST("hSameLambdaProton_SB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, collision.centFT0C()) / (trackEff(hEffProtons, associate, collision.centFT0C()) * v0Efficiency));
464-
if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass < MassLambda0 - 4 * dGaussSigma) {
473+
if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass < MassLambda0 - invMassNSigma*dGaussSigma) {
465474
rSECorrRegistry.fill(HIST("hSameLambdaProton_leftSB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, collision.centFT0C()) / (trackEff(hEffProtons, associate, collision.centFT0C()) * v0Efficiency));
466-
} else if (candMass > MassLambda0 + 4 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
475+
} else if (candMass > MassLambda0 + invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
467476
rSECorrRegistry.fill(HIST("hSameLambdaProton_rightSB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, collision.centFT0C()) / (trackEff(hEffProtons, associate, collision.centFT0C()) * v0Efficiency));
468477
}
469478
}
@@ -487,8 +496,8 @@ struct ThreeParticleCorrelations {
487496

488497
auto bc = coll_1.bc_as<aod::BCsWithTimestamps>();
489498
auto bField = getMagneticField(bc.timestamp());
490-
if (confBfieldSwitch != 0) {
491-
if (std::signbit(static_cast<double>(confBfieldSwitch)) != std::signbit(bField)) {
499+
if (switchGroup.confBfieldSwitch != 0) {
500+
if (std::signbit(static_cast<double>(switchGroup.confBfieldSwitch)) != std::signbit(bField)) {
492501
return;
493502
}
494503
}
@@ -510,7 +519,7 @@ struct ThreeParticleCorrelations {
510519
deltaPhi = RecoDecay::constrainAngle(trigger.phi() - associate.phi(), -constants::math::PIHalf);
511520
deltaEta = trigger.eta() - associate.eta();
512521

513-
if (candMass >= MassLambda0 - 4 * dGaussSigma && candMass <= MassLambda0 + 4 * dGaussSigma) {
522+
if (candMass >= MassLambda0 - invMassNSigma*dGaussSigma && candMass <= MassLambda0 + invMassNSigma*dGaussSigma) {
514523
if (assocPID[0] == pionID) { // Pions
515524
rMECorrRegistry.fill(HIST("hMixLambdaPion_SGNL"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, coll_1.centFT0C()) / (trackEff(hEffPions, associate, coll_1.centFT0C()) * v0Efficiency));
516525
} else if (assocPID[0] == kaonID) { // Kaons
@@ -519,28 +528,28 @@ struct ThreeParticleCorrelations {
519528
rMECorrRegistry.fill(HIST("hMixLambdaProton_SGNL"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, coll_1.centFT0C()) / (trackEff(hEffProtons, associate, coll_1.centFT0C()) * v0Efficiency));
520529
}
521530

522-
} else if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
531+
} else if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
523532
if (assocPID[0] == pionID) { // Pions
524533
rMECorrRegistry.fill(HIST("hMixLambdaPion_SB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, coll_1.centFT0C()) / (trackEff(hEffPions, associate, coll_1.centFT0C()) * v0Efficiency));
525-
if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass < MassLambda0 - 4 * dGaussSigma) {
534+
if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass < MassLambda0 - invMassNSigma*dGaussSigma) {
526535
rMECorrRegistry.fill(HIST("hMixLambdaPion_leftSB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, coll_1.centFT0C()) / (trackEff(hEffPions, associate, coll_1.centFT0C()) * v0Efficiency));
527-
} else if (candMass > MassLambda0 + 4 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
536+
} else if (candMass > MassLambda0 + invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
528537
rMECorrRegistry.fill(HIST("hMixLambdaPion_rightSB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionPions, associate, coll_1.centFT0C()) / (trackEff(hEffPions, associate, coll_1.centFT0C()) * v0Efficiency));
529538
}
530539

531540
} else if (assocPID[0] == kaonID) { // Kaons
532541
rMECorrRegistry.fill(HIST("hMixLambdaKaon_SB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionKaons, associate, coll_1.centFT0C()) / (trackEff(hEffKaons, associate, coll_1.centFT0C()) * v0Efficiency));
533-
if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass < MassLambda0 - 4 * dGaussSigma) {
542+
if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass < MassLambda0 - invMassNSigma*dGaussSigma) {
534543
rMECorrRegistry.fill(HIST("hMixLambdaKaon_leftSB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionKaons, associate, coll_1.centFT0C()) / (trackEff(hEffKaons, associate, coll_1.centFT0C()) * v0Efficiency));
535-
} else if (candMass > MassLambda0 + 4 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
544+
} else if (candMass > MassLambda0 + invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
536545
rMECorrRegistry.fill(HIST("hMixLambdaKaon_rightSB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionKaons, associate, coll_1.centFT0C()) / (trackEff(hEffKaons, associate, coll_1.centFT0C()) * v0Efficiency));
537546
}
538547

539548
} else if (assocPID[0] == protonID) { // Protons
540549
rMECorrRegistry.fill(HIST("hMixLambdaProton_SB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, coll_1.centFT0C()) / (trackEff(hEffProtons, associate, coll_1.centFT0C()) * v0Efficiency));
541-
if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass < MassLambda0 - 4 * dGaussSigma) {
550+
if (candMass >= MassLambda0 - 2*invMassNSigma*dGaussSigma && candMass < MassLambda0 - invMassNSigma*dGaussSigma) {
542551
rMECorrRegistry.fill(HIST("hMixLambdaProton_leftSB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, coll_1.centFT0C()) / (trackEff(hEffProtons, associate, coll_1.centFT0C()) * v0Efficiency));
543-
} else if (candMass > MassLambda0 + 4 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
552+
} else if (candMass > MassLambda0 + invMassNSigma*dGaussSigma && candMass <= MassLambda0 + 2*invMassNSigma*dGaussSigma) {
544553
rMECorrRegistry.fill(HIST("hMixLambdaProton_rightSB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), ratioCorrection(hCorrectionProtons, associate, coll_1.centFT0C()) / (trackEff(hEffProtons, associate, coll_1.centFT0C()) * v0Efficiency));
545554
}
546555
}
@@ -909,7 +918,7 @@ struct ThreeParticleCorrelations {
909918
{
910919

911920
double ratioCorrection = 1.0;
912-
if (confRatioCorrectionSwitch) {
921+
if (switchGroup.confRatioCorrectionSwitch) {
913922

914923
int index = -999;
915924
if (track.sign() > 0) {
@@ -1009,37 +1018,37 @@ struct ThreeParticleCorrelations {
10091018
if (std::abs(posDaughter.eta()) >= trackEtaMax || std::abs(negDaughter.eta()) >= trackEtaMax) {
10101019
return false;
10111020
}
1012-
if (posDaughter.tpcNClsCrossedRows() <= tpcNCrossedRowsMin || negDaughter.tpcNClsCrossedRows() <= tpcNCrossedRowsMin) {
1021+
if (posDaughter.tpcNClsCrossedRows() <= v0SelGroup.tpcNCrossedRows || negDaughter.tpcNClsCrossedRows() <= v0SelGroup.tpcNCrossedRows) {
10131022
return false;
10141023
}
10151024
if (v0Sign(v0) == 1) {
10161025
if (std::abs(posDaughter.tpcNSigmaPr()) >= nSigma5 || std::abs(negDaughter.tpcNSigmaPi()) >= nSigma5) {
10171026
return false;
10181027
}
1019-
if (std::abs(v0.dcapostopv()) <= dcaProtonMin || std::abs(v0.dcanegtopv()) <= dcaPionMin) {
1028+
if (std::abs(v0.dcapostopv()) <= v0SelGroup.dcaProton || std::abs(v0.dcanegtopv()) <= v0SelGroup.dcaPion) {
10201029
return false;
10211030
}
10221031
} else if (v0Sign(v0) == -1) {
10231032
if (std::abs(posDaughter.tpcNSigmaPi()) >= nSigma5 || std::abs(negDaughter.tpcNSigmaPr()) >= nSigma5) {
10241033
return false;
10251034
}
1026-
if (std::abs(v0.dcapostopv()) <= dcaPionMin || std::abs(v0.dcanegtopv()) <= dcaProtonMin) {
1035+
if (std::abs(v0.dcapostopv()) <= v0SelGroup.dcaPion || std::abs(v0.dcanegtopv()) <= v0SelGroup.dcaProton) {
10271036
return false;
10281037
}
10291038
}
10301039

10311040
// Topological cuts
10321041
float ctau = v0.distovertotmom(col.posX(), col.posY(), col.posZ()) * MassLambda0;
1033-
if (v0.v0radius() <= decayRMin) {
1042+
if (v0.v0radius() <= v0SelGroup.decayR) {
10341043
return false;
10351044
}
1036-
if (ctau >= ctauMax) {
1045+
if (ctau >= v0SelGroup.ctau) {
10371046
return false;
10381047
}
1039-
if (v0.v0cosPA() <= cosPAMin) {
1048+
if (v0.v0cosPA() <= v0SelGroup.cosPA) {
10401049
return false;
10411050
}
1042-
if (v0.dcaV0daughters() >= dcaV0DauMax) {
1051+
if (v0.dcaV0daughters() >= v0SelGroup.dcaV0Dau) {
10431052
return false;
10441053
}
10451054

@@ -1133,7 +1142,7 @@ struct ThreeParticleCorrelations {
11331142
bool fakeV0Filter(const V0Cand& v0, const TrackCand& track)
11341143
{
11351144

1136-
if (confFakeV0Switch) {
1145+
if (switchGroup.confFakeV0Switch) {
11371146

11381147
if (trackPID(track)[0] == kaonID) { // Kaons
11391148
return true;
@@ -1165,7 +1174,7 @@ struct ThreeParticleCorrelations {
11651174
}
11661175

11671176
double invMass = RecoDecay::m(std::array{dMomArray, aMomArray}, massArray);
1168-
if (invMass >= MassLambda0 - 4 * dGaussSigma && invMass <= MassLambda0 + 4 * dGaussSigma) {
1177+
if (invMass >= MassLambda0 - invMassNSigma*dGaussSigma && invMass <= MassLambda0 + invMassNSigma*dGaussSigma) {
11691178
return false;
11701179
}
11711180
}
@@ -1178,7 +1187,7 @@ struct ThreeParticleCorrelations {
11781187
{
11791188

11801189
bool pass = true;
1181-
if (confRDSwitch) {
1190+
if (switchGroup.confRDSwitch) {
11821191

11831192
auto proton = v0.template posTrack_as<MyFilteredTracks>();
11841193
if (v0Sign(v0) == -1) {

0 commit comments

Comments
 (0)