Skip to content

Commit f2594ed

Browse files
author
Chiara De Martin
committed
add resolution histogrmas
1 parent d2a36a2 commit f2594ed

File tree

1 file changed

+73
-31
lines changed

1 file changed

+73
-31
lines changed

‎PWGLF/TableProducer/Strangeness/cascadeflow.cxx‎

Lines changed: 73 additions & 31 deletions
Original file line numberDiff line numberDiff line change
@@ -476,15 +476,41 @@ struct cascadeFlow {
476476
double ApplyShiftCorrection(TCollision coll, double psiT0C){
477477
int nmode = 2;
478478
auto deltapsiFT0C = 0.0;
479+
479480
for (int ishift = 1; ishift <= 10; ishift++) {
480481
auto coeffshiftxFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 0.5, ishift - 0.5));
481482
auto coeffshiftyFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 1.5, ishift - 0.5));
482-
483+
483484
deltapsiFT0C += ((1 / (1.0 * ishift)) * (-coeffshiftxFT0C * TMath::Cos(ishift * static_cast<float>(nmode) * psiT0C) + coeffshiftyFT0C * TMath::Sin(ishift * static_cast<float>(nmode) * psiT0C)));
484485
}
485486
return psiT0C + deltapsiFT0C;
486487
}
487-
488+
489+
template <typename TCollision>
490+
double ComputeEPResolutionwShifts(TCollision coll, double psiT0C, double psiTPCA, double psiTPCC){
491+
int nmode = 2;
492+
auto deltapsiFT0C = 0.0;
493+
auto deltapsiTPCA = 0.0;
494+
auto deltapsiTPCC = 0.0;
495+
for (int ishift = 1; ishift <= 10; ishift++) {
496+
auto coeffshiftxFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 0.5, ishift - 0.5));
497+
auto coeffshiftyFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 1.5, ishift - 0.5));
498+
auto coeffshiftxTPCA = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 2.5, ishift - 0.5));
499+
auto coeffshiftyTPCA = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 3.5, ishift - 0.5));
500+
auto coeffshiftxTPCC = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 4.5, ishift - 0.5));
501+
auto coeffshiftyTPCC = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 5.5, ishift - 0.5));
502+
deltapsiFT0C += ((1 / (1.0 * ishift)) * (-coeffshiftxFT0C * TMath::Cos(ishift * static_cast<float>(nmode) * psiT0C) + coeffshiftyFT0C * TMath::Sin(ishift * static_cast<float>(nmode) * psiT0C)));
503+
deltapsiTPCA += ((1 / (1.0 * ishift)) * (-coeffshiftxTPCA * TMath::Cos(ishift * static_cast<float>(nmode) * psiTPCA) + coeffshiftyTPCA * TMath::Sin(ishift * static_cast<float>(nmode) * psiTPCA)));
504+
deltapsiTPCC += ((1 / (1.0 * ishift)) * (-coeffshiftxTPCC * TMath::Cos(ishift * static_cast<float>(nmode) * psiTPCC) + coeffshiftyTPCC * TMath::Sin(ishift * static_cast<float>(nmode) * psiTPCC)));
505+
}
506+
//histos.fill(HIST("psi2/QA/EP_FT0C_shifted"), coll.centFT0C(), psiT0C + deltapsiFT0C);
507+
//histos.fill(HIST("psi2/QA/EP_TPCA_shifted"), coll.centFT0C(), psiTPCA + deltapsiTPCA);
508+
//histos.fill(HIST("psi2/QA/EP_TPCC_shifted"), coll.centFT0C(), psiTPCC + deltapsiTPCC);
509+
resolution.fill(HIST("QVectorsT0CTPCA_Shifted"), coll.centFT0C(), TMath::Cos(static_cast<float>(nmode) * (psiT0C + deltapsiFT0C - psiTPCA - deltapsiTPCA)));
510+
resolution.fill(HIST("QVectorsT0CTPCC_Shifted"), coll.centFT0C(), TMath::Cos(static_cast<float>(nmode) * (psiT0C + deltapsiFT0C - psiTPCC - deltapsiTPCC)));
511+
resolution.fill(HIST("QVectorsTPCAC_Shifted"), coll.centFT0C(), TMath::Cos(static_cast<float>(nmode) * (psiTPCA + deltapsiTPCA - psiTPCC - deltapsiTPCC)));
512+
}
513+
488514
// objects to use for acceptance correction
489515
TH2F* hAcceptanceXi;
490516
TH2F* hAcceptanceOmega;
@@ -684,6 +710,9 @@ struct cascadeFlow {
684710
resolution.add("QVectorsNormT0CTPCC", "QVectorsNormT0CTPCC", HistType::kTH2F, {axisQVsNorm, CentAxis});
685711
resolution.add("QVectorsNormTPCAC", "QVectorsNormTPCCB", HistType::kTH2F, {axisQVsNorm, CentAxis});
686712
resolution.add("QVectorsSpecPlane", "QVectorsSpecPlane", HistType::kTH2F, {axisQVsNorm, CentAxis});
713+
resolution.add("QVectorsT0CTPCA_Shifted", "QVectorsT0CTPCA_Shifted", HistType::kTH2F, {axisQVs, CentAxis});
714+
resolution.add("QVectorsT0CTPCC_Shifted", "QVectorsT0CTPCC_Shifted", HistType::kTH2F, {axisQVs, CentAxis});
715+
resolution.add("QVectorsTPCAC_Shifted", "QVectorsTPCAC_Shifted", HistType::kTH2F, {axisQVs, CentAxis});
687716

688717
histos.add("hNEvents", "hNEvents", {HistType::kTH1F, {{10, 0.f, 10.f}}});
689718
for (Int_t n = 1; n <= histos.get<TH1>(HIST("hNEvents"))->GetNbinsX(); n++) {
@@ -995,8 +1024,11 @@ struct cascadeFlow {
9951024
ROOT::Math::XYZVector spectatorplaneVecZDCA{std::cos(coll.psiZDCA()), std::sin(coll.psiZDCA()), 0}; // eta positive = projectile
9961025
ROOT::Math::XYZVector spectatorplaneVecZDCC{std::cos(coll.psiZDCC()), std::sin(coll.psiZDCC()), 0}; // eta negative = target
9971026

998-
const float PsiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
999-
float PsiT0CCorr = PsiT0C;
1027+
const float psiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1028+
const float psiTPCA = std::atan2(coll.qvecBPosIm(), coll.qvecBPosRe()) * 0.5f;
1029+
const float psiTPCC = std::atan2(coll.qvecBNegIm(), coll.qvecBNegRe()) * 0.5f;
1030+
float psiT0CCorr = psiT0C;
1031+
10001032
if (cfgShiftCorr) {
10011033
currentRunNumber = coll.runNumber();
10021034
if (currentRunNumber != lastRunNumber) {
@@ -1012,12 +1044,13 @@ struct cascadeFlow {
10121044
}
10131045
}
10141046
if (cfgShiftCorr){
1015-
PsiT0CCorr = ApplyShiftCorrection(coll, PsiT0C);
1047+
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
1048+
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
10161049
}
10171050

1018-
histos.fill(HIST("hPsiT0C"), PsiT0CCorr);
1051+
histos.fill(HIST("hPsiT0C"), psiT0CCorr);
10191052
histos.fill(HIST("hPsiZDCA_vs_ZDCC"), coll.psiZDCC(), coll.psiZDCA());
1020-
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0CCorr);
1053+
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), psiT0CCorr);
10211054

10221055
resolution.fill(HIST("QVectorsT0CTPCA"), eventplaneVecT0C.Dot(eventplaneVecTPCA), coll.centFT0C());
10231056
resolution.fill(HIST("QVectorsT0CTPCC"), eventplaneVecT0C.Dot(eventplaneVecTPCC), coll.centFT0C());
@@ -1093,7 +1126,7 @@ struct cascadeFlow {
10931126

10941127
ROOT::Math::XYZVector cascQvec{std::cos(2 * casc.phi()), std::sin(2 * casc.phi()), 0};
10951128
auto v2CSP = cascQvec.Dot(eventplaneVecT0C); // not normalised by amplitude
1096-
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0CCorr);
1129+
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - psiT0CCorr);
10971130
auto v2CEP = std::cos(2.0 * cascminuspsiT0C);
10981131
ROOT::Math::XYZVector cascUvec{std::cos(casc.phi()), std::sin(casc.phi()), 0};
10991132
auto v1SP_ZDCA = cascUvec.Dot(spectatorplaneVecZDCA);
@@ -1149,11 +1182,11 @@ struct cascadeFlow {
11491182
int chargeIndex = 0;
11501183
if (casc.sign() > 0)
11511184
chargeIndex = 1;
1152-
double pzs2Xi = cosThetaStarLambda[0] * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaXi[chargeIndex] / meanCos2ThetaLambdaFromXi;
1153-
double pzs2Omega = cosThetaStarLambda[1] * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaOmega[chargeIndex] / meanCos2ThetaLambdaFromOmega;
1185+
double pzs2Xi = cosThetaStarLambda[0] * std::sin(2 * (casc.phi() - psiT0CCorr)) / cascadev2::AlphaXi[chargeIndex] / meanCos2ThetaLambdaFromXi;
1186+
double pzs2Omega = cosThetaStarLambda[1] * std::sin(2 * (casc.phi() - psiT0CCorr)) / cascadev2::AlphaOmega[chargeIndex] / meanCos2ThetaLambdaFromOmega;
11541187
double cos2ThetaXi = cosThetaStarLambda[0] * cosThetaStarLambda[0];
11551188
double cos2ThetaOmega = cosThetaStarLambda[1] * cosThetaStarLambda[1];
1156-
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
1189+
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - psiT0CCorr)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
11571190
double cos2ThetaLambda = cosThetaStarProton * cosThetaStarProton;
11581191

11591192
double cosThetaXiWithAlpha = cosThetaStarLambda[0] / cascadev2::AlphaXi[chargeIndex];
@@ -1259,7 +1292,7 @@ struct cascadeFlow {
12591292

12601293
if (isSelectedCasc[0] || isSelectedCasc[1]) {
12611294
if (fillingConfigs.isFillTree)
1262-
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, PsiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
1295+
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, psiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
12631296
}
12641297
}
12651298
}
@@ -1295,8 +1328,10 @@ struct cascadeFlow {
12951328
ROOT::Math::XYZVector eventplaneVecTPCA{coll.qvecBPosRe(), coll.qvecBPosIm(), 0};
12961329
ROOT::Math::XYZVector eventplaneVecTPCC{coll.qvecBNegRe(), coll.qvecBNegIm(), 0};
12971330

1298-
const float PsiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1299-
float PsiT0CCorr = PsiT0C;
1331+
const float psiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1332+
const float psiTPCA = std::atan2(coll.qvecBPosIm(), coll.qvecBPosRe()) * 0.5f;
1333+
const float psiTPCC = std::atan2(coll.qvecBNegIm(), coll.qvecBNegRe()) * 0.5f;
1334+
float psiT0CCorr = psiT0C;
13001335
if (cfgShiftCorr) {
13011336
currentRunNumber = coll.runNumber();
13021337
if (currentRunNumber != lastRunNumber) {
@@ -1312,10 +1347,11 @@ struct cascadeFlow {
13121347
}
13131348
}
13141349
if (cfgShiftCorr){
1315-
PsiT0CCorr = ApplyShiftCorrection(coll, PsiT0C);
1350+
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
1351+
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
13161352
}
1317-
histos.fill(HIST("hPsiT0C"), PsiT0CCorr);
1318-
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0CCorr);
1353+
histos.fill(HIST("hPsiT0C"), psiT0CCorr);
1354+
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), psiT0CCorr);
13191355

13201356
resolution.fill(HIST("QVectorsT0CTPCA"), eventplaneVecT0C.Dot(eventplaneVecTPCA), coll.centFT0C());
13211357
resolution.fill(HIST("QVectorsT0CTPCC"), eventplaneVecT0C.Dot(eventplaneVecTPCC), coll.centFT0C());
@@ -1390,7 +1426,7 @@ struct cascadeFlow {
13901426

13911427
ROOT::Math::XYZVector cascQvec{std::cos(2 * casc.phi()), std::sin(2 * casc.phi()), 0};
13921428
auto v2CSP = cascQvec.Dot(eventplaneVecT0C); // not normalised by amplitude
1393-
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0CCorr);
1429+
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - psiT0CCorr);
13941430
auto v2CEP = std::cos(2.0 * cascminuspsiT0C);
13951431
ROOT::Math::XYZVector cascUvec{std::cos(casc.phi()), std::sin(casc.phi()), 0};
13961432

@@ -1440,11 +1476,11 @@ struct cascadeFlow {
14401476
int chargeIndex = 0;
14411477
if (casc.sign() > 0)
14421478
chargeIndex = 1;
1443-
double pzs2Xi = cosThetaStarLambda[0] * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaXi[chargeIndex] / meanCos2ThetaLambdaFromXi;
1444-
double pzs2Omega = cosThetaStarLambda[1] * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaOmega[chargeIndex] / meanCos2ThetaLambdaFromOmega;
1479+
double pzs2Xi = cosThetaStarLambda[0] * std::sin(2 * (casc.phi() - psiT0CCorr)) / cascadev2::AlphaXi[chargeIndex] / meanCos2ThetaLambdaFromXi;
1480+
double pzs2Omega = cosThetaStarLambda[1] * std::sin(2 * (casc.phi() - psiT0CCorr)) / cascadev2::AlphaOmega[chargeIndex] / meanCos2ThetaLambdaFromOmega;
14451481
double cos2ThetaXi = cosThetaStarLambda[0] * cosThetaStarLambda[0];
14461482
double cos2ThetaOmega = cosThetaStarLambda[1] * cosThetaStarLambda[1];
1447-
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
1483+
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - psiT0CCorr)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
14481484
double cos2ThetaLambda = cosThetaStarProton * cosThetaStarProton;
14491485

14501486
double cosThetaXiWithAlpha = cosThetaStarLambda[0] / cascadev2::AlphaXi[chargeIndex];
@@ -1548,7 +1584,7 @@ struct cascadeFlow {
15481584

15491585
if (isSelectedCasc[0] || isSelectedCasc[1]) {
15501586
if (fillingConfigs.isFillTree)
1551-
fillAnalysedTable(coll, hasEventPlane, 0, casc, v2CSP, v2CEP, 0, 0, PsiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
1587+
fillAnalysedTable(coll, hasEventPlane, 0, casc, v2CSP, v2CEP, 0, 0, psiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
15521588
}
15531589
}
15541590
}
@@ -1583,6 +1619,8 @@ struct cascadeFlow {
15831619
ROOT::Math::XYZVector eventplaneVecTPCC{coll.qvecBNegRe(), coll.qvecBNegIm(), 0};
15841620

15851621
const float psiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1622+
const float psiTPCA = std::atan2(coll.qvecBPosIm(), coll.qvecBPosRe()) * 0.5f;
1623+
const float psiTPCC = std::atan2(coll.qvecBNegIm(), coll.qvecBNegRe()) * 0.5f;
15861624
float psiT0CCorr = psiT0C;
15871625
if (cfgShiftCorr) {
15881626
currentRunNumber = coll.runNumber();
@@ -1600,6 +1638,7 @@ struct cascadeFlow {
16001638
}
16011639
if (cfgShiftCorr){
16021640
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
1641+
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
16031642
}
16041643
histos.fill(HIST("hPsiT0C"), psiT0CCorr);
16051644
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), psiT0CCorr);
@@ -1784,8 +1823,10 @@ struct cascadeFlow {
17841823
float NormQvTPCA = std::sqrt(eventplaneVecTPCA.Dot(eventplaneVecTPCA));
17851824
float NormQvTPCC = std::sqrt(eventplaneVecTPCC.Dot(eventplaneVecTPCC));
17861825

1787-
const float PsiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1788-
float PsiT0CCorr = PsiT0C;
1826+
const float psiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1827+
const float psiTPCA = std::atan2(coll.qvecBPosIm(), coll.qvecBPosRe()) * 0.5f;
1828+
const float psiTPCC = std::atan2(coll.qvecBNegIm(), coll.qvecBNegRe()) * 0.5f;
1829+
float psiT0CCorr = psiT0C;
17891830
if (cfgShiftCorr) {
17901831
currentRunNumber = coll.runNumber();
17911832
if (currentRunNumber != lastRunNumber) {
@@ -1801,10 +1842,11 @@ struct cascadeFlow {
18011842
}
18021843
}
18031844
if (cfgShiftCorr){
1804-
PsiT0CCorr = ApplyShiftCorrection(coll, PsiT0C);
1845+
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
1846+
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
18051847
}
1806-
histos.fill(HIST("hPsiT0C"), PsiT0CCorr);
1807-
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0CCorr);
1848+
histos.fill(HIST("hpsiT0C"), psiT0CCorr);
1849+
histos.fill(HIST("hpsiT0CvsCentFT0C"), coll.centFT0C(), psiT0CCorr);
18081850

18091851
resolution.fill(HIST("QVectorsT0CTPCA"), eventplaneVecT0C.Dot(eventplaneVecTPCA), coll.centFT0C());
18101852
resolution.fill(HIST("QVectorsT0CTPCC"), eventplaneVecT0C.Dot(eventplaneVecTPCC), coll.centFT0C());
@@ -1880,7 +1922,7 @@ struct cascadeFlow {
18801922

18811923
ROOT::Math::XYZVector cascQvec{std::cos(2 * casc.phi()), std::sin(2 * casc.phi()), 0};
18821924
auto v2CSP = cascQvec.Dot(eventplaneVecT0C);
1883-
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0CCorr);
1925+
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - psiT0CCorr);
18841926
auto v2CEP = std::cos(2.0 * cascminuspsiT0C);
18851927
ROOT::Math::XYZVector cascUvec{std::cos(casc.phi()), std::sin(casc.phi()), 0};
18861928
auto v1SP_ZDCA = cascUvec.Dot(spectatorplaneVecZDCA);
@@ -1912,7 +1954,7 @@ struct cascadeFlow {
19121954
}
19131955
if (isSelectedCasc[0] || isSelectedCasc[1])
19141956
if (fillingConfigs.isFillTree)
1915-
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, PsiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
1957+
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, psiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
19161958
}
19171959
}
19181960

@@ -2030,7 +2072,7 @@ struct cascadeFlow {
20302072
histos.fill(HIST("hCascadePhi"), casc.phi());
20312073

20322074
float BDTresponse[nParticles]{0.f, 0.f};
2033-
const float PsiT0C = 0; // not defined in MC for now
2075+
const float psiT0C = 0; // not defined in MC for now
20342076
auto v2CSP = 0; // not defined in MC for now
20352077
auto v2CEP = 0; // not defined in MC for now
20362078
auto v1SP_ZDCA = 0; // not defined in MC for now
@@ -2045,7 +2087,7 @@ struct cascadeFlow {
20452087
continue;
20462088
}
20472089
if (isSelectedCasc[0] || isSelectedCasc[1])
2048-
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, PsiT0C, BDTresponse[0], BDTresponse[1], pdgCode);
2090+
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, psiT0C, BDTresponse[0], BDTresponse[1], pdgCode);
20492091
}
20502092
}
20512093

0 commit comments

Comments
 (0)