Skip to content

Commit 06e871b

Browse files
JStaaJStaa
committed
Changed the RD method and the V0 MC efficiency
1 parent 8ecfb76 commit 06e871b

File tree

1 file changed

+40
-99
lines changed

1 file changed

+40
-99
lines changed

‎PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx‎

Lines changed: 40 additions & 99 deletions
Original file line numberDiff line numberDiff line change
@@ -59,8 +59,8 @@ struct ThreeParticleCorrelations {
5959
float pionPtMid = 1.5, kaonPtMid1 = 1.5, kaonPtMid2 = 2.0, protonPtMid = 0.7;
6060

6161
// RD filter parameters
62-
float dEtaMax = 0.05, dEtaMin = 0.022;
63-
float dPhiStarMinOS = 0.075, dPhiStarMinSS = 0.12;
62+
float dEtaMax = 0.05, dEtaMin = 0.023;
63+
float dPhiStarMinOS = 0.09, dPhiStarMinSS = 0.095;
6464
float rMin = 0.8, rMax = 2.5;
6565

6666
// Lambda invariant mass fit
@@ -208,27 +208,19 @@ struct ThreeParticleCorrelations {
208208
rQARegistry.add("hInvMassAntiLambda_MC", "hInvMassAntiLambda_MC", {HistType::kTH3D, {{lambdaInvMassAxis}, {v0PtAxis}, {centralityAxis}}});
209209

210210
// PhiStar
211-
rPhiStarRegistry.add("hSEProtonPreCut_OS", "hSEProtonPreCut_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
212-
rPhiStarRegistry.add("hSEProtonPreCut_SS", "hSEProtonPreCut_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
213-
rPhiStarRegistry.add("hSEProtonPreCut_SSP", "hSEProtonPreCut_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
214-
rPhiStarRegistry.add("hSEProtonPreCut_SSN", "hSEProtonPreCut_SSN", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
215-
rPhiStarRegistry.add("hSEProtonPostCut_OS", "hSEProtonPostCut_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
216-
rPhiStarRegistry.add("hSEProtonPostCut_SS", "hSEProtonPostCut_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
217-
rPhiStarRegistry.add("hSEProtonPostCut_SSP", "hSEProtonPostCut_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
218-
rPhiStarRegistry.add("hSEProtonPostCut_SSN", "hSEProtonPostCut_SSN", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
211+
rPhiStarRegistry.add("hSEPhiStarIR_OS", "hSEPhiStarIR_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
212+
rPhiStarRegistry.add("hSEPhiStarIR_SS", "hSEPhiStarIR_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
213+
rPhiStarRegistry.add("hSEPhiStarIR_SSP", "hSEPhiStarIR_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
214+
rPhiStarRegistry.add("hSEPhiStarIR_SSN", "hSEPhiStarIR_SSN", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
219215
rPhiStarRegistry.add("hSEPhiStarMean_OS", "hSEPhiStarMean_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
220216
rPhiStarRegistry.add("hSEPhiStarMean_SS", "hSEPhiStarMean_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
221217
rPhiStarRegistry.add("hSEPhiStarMean_SSP", "hSEPhiStarMean_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
222218
rPhiStarRegistry.add("hSEPhiStarMean_SSN", "hSEPhiStarMean_SSN", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
223219

224-
rPhiStarRegistry.add("hMEProtonPreCut_OS", "hMEProtonPreCut_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
225-
rPhiStarRegistry.add("hMEProtonPreCut_SS", "hMEProtonPreCut_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
226-
rPhiStarRegistry.add("hMEProtonPreCut_SSP", "hMEProtonPreCut_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
227-
rPhiStarRegistry.add("hMEProtonPreCut_SSN", "hMEProtonPreCut_SSN", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
228-
rPhiStarRegistry.add("hMEProtonPostCut_OS", "hMEProtonPostCut_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
229-
rPhiStarRegistry.add("hMEProtonPostCut_SS", "hMEProtonPostCut_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
230-
rPhiStarRegistry.add("hMEProtonPostCut_SSP", "hMEProtonPostCut_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
231-
rPhiStarRegistry.add("hMEProtonPostCut_SSN", "hMEProtonPostCut_SSN", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
220+
rPhiStarRegistry.add("hMEPhiStarIR_OS", "hMEPhiStarIR_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
221+
rPhiStarRegistry.add("hMEPhiStarIR_SS", "hMEPhiStarIR_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
222+
rPhiStarRegistry.add("hMEPhiStarIR_SSP", "hMEPhiStarIR_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
223+
rPhiStarRegistry.add("hMEPhiStarIR_SSN", "hMEPhiStarIR_SSN", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
232224
rPhiStarRegistry.add("hMEPhiStarMean_OS", "hMEPhiStarMean_OS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
233225
rPhiStarRegistry.add("hMEPhiStarMean_SS", "hMEPhiStarMean_SS", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
234226
rPhiStarRegistry.add("hMEPhiStarMean_SSP", "hMEPhiStarMean_SSP", {HistType::kTH2D, {{121, -0.3025, 0.3025}, {101, -0.0505, 0.0505}}});
@@ -260,10 +252,6 @@ struct ThreeParticleCorrelations {
260252
rMCRegistry.add("hPIDKaonN", "hPIDKaonN", {HistType::kTH3D, {{trackPtAxis}, {trackEtaAxis}, {centralityAxis}}});
261253
rMCRegistry.add("hPIDProtonP", "hPIDProtonP", {HistType::kTH3D, {{trackPtAxis}, {trackEtaAxis}, {centralityAxis}}});
262254
rMCRegistry.add("hPIDProtonN", "hPIDProtonN", {HistType::kTH3D, {{trackPtAxis}, {trackEtaAxis}, {centralityAxis}}});
263-
rMCRegistry.add("hPIDLambdaP_SGNL", "hPIDLambdaP_SGNL", {HistType::kTH3D, {{v0PtAxis}, {v0EtaAxis}, {centralityAxis}}});
264-
rMCRegistry.add("hPIDLambdaP_SB", "hPIDLambdaP_SB", {HistType::kTH3D, {{v0PtAxis}, {v0EtaAxis}, {centralityAxis}}});
265-
rMCRegistry.add("hPIDLambdaN_SGNL", "hPIDLambdaN_SGNL", {HistType::kTH3D, {{v0PtAxis}, {v0EtaAxis}, {centralityAxis}}});
266-
rMCRegistry.add("hPIDLambdaN_SB", "hPIDLambdaN_SB", {HistType::kTH3D, {{v0PtAxis}, {v0EtaAxis}, {centralityAxis}}});
267255

268256
// Purity
269257
rMCRegistry.add("hSelectPionP", "hSelectPionP", {HistType::kTH1D, {trackPtAxis}});
@@ -768,42 +756,22 @@ struct ThreeParticleCorrelations {
768756

769757
for (const auto& v0 : v0s) {
770758

771-
if (!v0.has_mcParticle() || v0.pt() <= v0PtMin || v0.pt() >= v0PtMax || std::abs(v0.eta()) >= v0EtaMax) {
759+
if (!v0.has_mcParticle()) {
772760
continue;
773761
}
774-
auto particle = v0.mcParticle();
762+
763+
if (v0Filters(collision, v0, tracks)) {
775764

776-
if (particle.isPhysicalPrimary()) {
777-
778-
// V0 efficiency - Reconstructed
779-
if (particle.pdgCode() == kLambda0) { // Lambdas
765+
// V0 efficiency - Reconstructed
766+
if (v0Sign(v0) == 1) { // Lambdas
767+
candMass = v0.mLambda();
768+
rQARegistry.fill(HIST("hInvMassLambda_MC"), v0.mLambda(), v0.pt(), collision.centFT0C());
780769
rMCRegistry.fill(HIST("hRecLambdaP"), v0.pt(), v0.eta(), collision.centFT0C());
781-
} else if (particle.pdgCode() == kLambda0Bar) { // AntiLambdas
770+
} else if (v0Sign(v0) == -1) { // AntiLambdas
771+
candMass = v0.mAntiLambda();
772+
rQARegistry.fill(HIST("hInvMassAntiLambda_MC"), v0.mAntiLambda(), v0.pt(), collision.centFT0C());
782773
rMCRegistry.fill(HIST("hRecLambdaN"), v0.pt(), v0.eta(), collision.centFT0C());
783774
}
784-
785-
if (v0Filters(collision, v0, tracks)) {
786-
787-
// V0 efficiency - Reconstructed
788-
if (v0Sign(v0) == 1) { // Lambdas
789-
candMass = v0.mLambda();
790-
rQARegistry.fill(HIST("hInvMassLambda_MC"), v0.mLambda(), v0.pt(), collision.centFT0C());
791-
if (candMass >= MassLambda0 - 4 * dGaussSigma && candMass <= MassLambda0 + 4 * dGaussSigma) {
792-
rMCRegistry.fill(HIST("hPIDLambdaP_SGNL"), v0.pt(), v0.eta(), collision.centFT0C());
793-
} else if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
794-
rMCRegistry.fill(HIST("hPIDLambdaP_SB"), v0.pt(), v0.eta(), collision.centFT0C());
795-
}
796-
797-
} else if (v0Sign(v0) == -1) { // AntiLambdas
798-
candMass = v0.mAntiLambda();
799-
rQARegistry.fill(HIST("hInvMassAntiLambda_MC"), v0.mAntiLambda(), v0.pt(), collision.centFT0C());
800-
if (candMass >= MassLambda0 - 4 * dGaussSigma && candMass <= MassLambda0 + 4 * dGaussSigma) {
801-
rMCRegistry.fill(HIST("hPIDLambdaN_SGNL"), v0.pt(), v0.eta(), collision.centFT0C());
802-
} else if (candMass >= MassLambda0 - 8 * dGaussSigma && candMass <= MassLambda0 + 8 * dGaussSigma) {
803-
rMCRegistry.fill(HIST("hPIDLambdaN_SB"), v0.pt(), v0.eta(), collision.centFT0C());
804-
}
805-
}
806-
}
807775
}
808776
}
809777
// End of the Monte-Carlo reconstructed QA
@@ -1132,68 +1100,35 @@ struct ThreeParticleCorrelations {
11321100
if (r == rMin) {
11331101
if (!Mix) { // Same-event
11341102
if (proton.sign() * track.sign() == -1) { // OS (Electric charge)
1135-
rPhiStarRegistry.fill(HIST("hSEProtonPreCut_OS"), dPhiStar, dEta);
1103+
rPhiStarRegistry.fill(HIST("hSEPhiStarIR_OS"), dPhiStar, dEta);
11361104
} else if (proton.sign() * track.sign() == 1) { // SS (Electric charge)
1137-
rPhiStarRegistry.fill(HIST("hSEProtonPreCut_SS"), dPhiStar, dEta);
1105+
rPhiStarRegistry.fill(HIST("hSEPhiStarIR_SS"), dPhiStar, dEta);
11381106
if (proton.sign() == 1) { // Positive
1139-
rPhiStarRegistry.fill(HIST("hSEProtonPreCut_SSP"), dPhiStar, dEta);
1107+
rPhiStarRegistry.fill(HIST("hSEPhiStarIR_SSP"), dPhiStar, dEta);
11401108
} else if (proton.sign() == -1) { // Negative
1141-
rPhiStarRegistry.fill(HIST("hSEProtonPreCut_SSN"), dPhiStar, dEta);
1109+
rPhiStarRegistry.fill(HIST("hSEPhiStarIR_SSN"), dPhiStar, dEta);
11421110
}
11431111
}
11441112

11451113
} else { // Mixed-event
11461114
if (proton.sign() * track.sign() == -1) { // OS (Electric charge)
1147-
rPhiStarRegistry.fill(HIST("hMEProtonPreCut_OS"), dPhiStar, dEta);
1115+
rPhiStarRegistry.fill(HIST("hMEPhiStarIR_OS"), dPhiStar, dEta);
11481116
} else if (proton.sign() * track.sign() == 1) { // SS (Electric charge)
1149-
rPhiStarRegistry.fill(HIST("hMEProtonPreCut_SS"), dPhiStar, dEta);
1117+
rPhiStarRegistry.fill(HIST("hMEPhiStarIR_SS"), dPhiStar, dEta);
11501118
if (proton.sign() == 1) { // Positive
1151-
rPhiStarRegistry.fill(HIST("hMEProtonPreCut_SSP"), dPhiStar, dEta);
1119+
rPhiStarRegistry.fill(HIST("hMEPhiStarIR_SSP"), dPhiStar, dEta);
11521120
} else if (proton.sign() == -1) { // Negative
1153-
rPhiStarRegistry.fill(HIST("hMEProtonPreCut_SSN"), dPhiStar, dEta);
1121+
rPhiStarRegistry.fill(HIST("hMEPhiStarIR_SSN"), dPhiStar, dEta);
11541122
}
11551123
}
11561124
}
1157-
}
11581125

1159-
if (std::abs(dEta) <= dEtaMin) {
1160-
if (proton.sign() * track.sign() == -1) { // OS (Electric charge)
1161-
if (std::abs(dPhiStar) <= dPhiStarMinOS) {
1162-
pass = false;
1163-
}
1164-
} else if (proton.sign() * track.sign() == 1) { // SS (Electric charge)
1165-
if (std::abs(dPhiStar) <= dPhiStarMinSS) {
1166-
pass = false;
1167-
}
1168-
}
1169-
}
1170-
1171-
if (r == rMin && pass) {
1172-
if (!Mix) { // Same-event
1173-
if (proton.sign() * track.sign() == -1) { // OS (Electric charge)
1174-
rPhiStarRegistry.fill(HIST("hSEProtonPostCut_OS"), dPhiStar, dEta);
1175-
} else if (proton.sign() * track.sign() == 1) { // SS (Electric charge)
1176-
rPhiStarRegistry.fill(HIST("hSEProtonPostCut_SS"), dPhiStar, dEta);
1177-
if (proton.sign() == 1) { // Positive
1178-
rPhiStarRegistry.fill(HIST("hSEProtonPostCut_SSP"), dPhiStar, dEta);
1179-
} else if (proton.sign() == -1) { // Negative
1180-
rPhiStarRegistry.fill(HIST("hSEProtonPostCut_SSN"), dPhiStar, dEta);
1181-
}
1182-
}
1183-
1184-
} else { // Mixed-event
1185-
if (proton.sign() * track.sign() == -1) { // OS (Electric charge)
1186-
rPhiStarRegistry.fill(HIST("hMEProtonPostCut_OS"), dPhiStar, dEta);
1187-
} else if (proton.sign() * track.sign() == 1) { // SS (Electric charge)
1188-
rPhiStarRegistry.fill(HIST("hMEProtonPostCut_SS"), dPhiStar, dEta);
1189-
if (proton.sign() == 1) { // Positive
1190-
rPhiStarRegistry.fill(HIST("hMEProtonPostCut_SSP"), dPhiStar, dEta);
1191-
} else if (proton.sign() == -1) { // Negative
1192-
rPhiStarRegistry.fill(HIST("hMEProtonPostCut_SSN"), dPhiStar, dEta);
1193-
}
1194-
}
1195-
}
1196-
}
1126+
if (proton.sign() * track.sign() == -1) { // OS (Electric charge)
1127+
if (std::abs(dEta) < dEtaMin && std::abs(dPhiStar) < dPhiStarMinOS) {
1128+
pass = false;
1129+
}
1130+
}
1131+
}
11971132

11981133
dPhiStarMean += (dPhiStar / 170);
11991134
}
@@ -1223,6 +1158,12 @@ struct ThreeParticleCorrelations {
12231158
}
12241159
}
12251160
}
1161+
1162+
if (proton.sign() * track.sign() == 1) { // SS (Electric charge)
1163+
if (std::abs(dEta) < dEtaMin && std::abs(dPhiStarMean) < dPhiStarMinSS) {
1164+
pass = false;
1165+
}
1166+
}
12261167
}
12271168

12281169
return pass;

0 commit comments

Comments
 (0)