Skip to content

Commit d2a36a2

Browse files
author
Chiara De Martin
committed
add shift correction for qvector
1 parent 9fd8fa2 commit d2a36a2

File tree

1 file changed

+125
-29
lines changed

1 file changed

+125
-29
lines changed

‎PWGLF/TableProducer/Strangeness/cascadeflow.cxx‎

Lines changed: 125 additions & 29 deletions
Original file line numberDiff line numberDiff line change
@@ -42,10 +42,10 @@ using namespace o2::framework::expressions;
4242
using std::array;
4343

4444
using DauTracks = soa::Join<aod::DauTrackExtras, aod::DauTrackTPCPIDs>;
45-
using CollEventPlane = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraFT0CQVsEv, aod::StraTPCQVs>::iterator;
46-
using CollEventPlaneCentralFW = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraTPCQVs>::iterator;
47-
using CollEventAndSpecPlane = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraFT0CQVsEv, aod::StraTPCQVs, aod::StraZDCSP>::iterator;
48-
using CollEventAndSpecPlaneCentralFW = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraTPCQVs, aod::StraZDCSP>::iterator;
45+
using CollEventPlane = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraFT0CQVsEv, aod::StraTPCQVs, aod::StraStamps>::iterator;
46+
using CollEventPlaneCentralFW = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraTPCQVs, aod::StraStamps>::iterator;
47+
using CollEventAndSpecPlane = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraFT0CQVsEv, aod::StraTPCQVs, aod::StraZDCSP, aod::StraStamps>::iterator;
48+
using CollEventAndSpecPlaneCentralFW = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraFT0CQVs, aod::StraTPCQVs, aod::StraZDCSP, aod::StraStamps>::iterator;
4949
using MCCollisionsStra = soa::Join<aod::StraMCCollisions, aod::StraMCCollMults>;
5050
using V0Candidates = soa::Join<aod::V0CollRefs, aod::V0Cores, aod::V0Extras>;
5151
using CascCandidates = soa::Join<aod::CascCollRefs, aod::CascCores, aod::CascExtras, aod::CascBBs>;
@@ -170,6 +170,11 @@ struct cascadeFlow {
170170
ConfigurableAxis axisQVs{"axisQVs", {500, -10.f, 10.f}, "axisQVs"};
171171
ConfigurableAxis axisQVsNorm{"axisQVsNorm", {200, -1.f, 1.f}, "axisQVsNorm"};
172172

173+
// Configurable for shift correction
174+
Configurable<bool> cfgShiftCorr{"cfgShiftCorr", 0, ""};
175+
Configurable<std::string> cfgShiftPath{"cfgShiftPath", "Users/j/junlee/Qvector/QvecCalib/Shift", "Path for Shift"};
176+
Configurable<int> cfgnMods{"cfgnMods", 1, "The number of modulations of interest starting from 2"};
177+
173178
// THN axes
174179
ConfigurableAxis thnConfigAxisFT0C{"thnConfigAxisFT0C", {8, 0, 80}, "FT0C centrality (%)"};
175180
ConfigurableAxis thnConfigAxisEta{"thnConfigAxisEta", {8, -0.8, 0.8}, "pseudorapidity"};
@@ -462,6 +467,24 @@ struct cascadeFlow {
462467
return phi;
463468
}
464469

470+
int currentRunNumber = -999;
471+
int lastRunNumber = -999;
472+
std::vector<TProfile3D*> shiftprofile{};
473+
std::string fullCCDBShiftCorrPath;
474+
475+
template <typename TCollision>
476+
double ApplyShiftCorrection(TCollision coll, double psiT0C){
477+
int nmode = 2;
478+
auto deltapsiFT0C = 0.0;
479+
for (int ishift = 1; ishift <= 10; ishift++) {
480+
auto coeffshiftxFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 0.5, ishift - 0.5));
481+
auto coeffshiftyFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 1.5, ishift - 0.5));
482+
483+
deltapsiFT0C += ((1 / (1.0 * ishift)) * (-coeffshiftxFT0C * TMath::Cos(ishift * static_cast<float>(nmode) * psiT0C) + coeffshiftyFT0C * TMath::Sin(ishift * static_cast<float>(nmode) * psiT0C)));
484+
}
485+
return psiT0C + deltapsiFT0C;
486+
}
487+
465488
// objects to use for acceptance correction
466489
TH2F* hAcceptanceXi;
467490
TH2F* hAcceptanceOmega;
@@ -973,9 +996,28 @@ struct cascadeFlow {
973996
ROOT::Math::XYZVector spectatorplaneVecZDCC{std::cos(coll.psiZDCC()), std::sin(coll.psiZDCC()), 0}; // eta negative = target
974997

975998
const float PsiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
976-
histos.fill(HIST("hPsiT0C"), PsiT0C);
999+
float PsiT0CCorr = PsiT0C;
1000+
if (cfgShiftCorr) {
1001+
currentRunNumber = coll.runNumber();
1002+
if (currentRunNumber != lastRunNumber) {
1003+
shiftprofile.clear();
1004+
for (int i = 2; i < cfgnMods + 2; i++) {
1005+
fullCCDBShiftCorrPath = cfgShiftPath;
1006+
fullCCDBShiftCorrPath += "/v";
1007+
fullCCDBShiftCorrPath += std::to_string(i);
1008+
auto objshift = ccdb->getForTimeStamp<TProfile3D>(fullCCDBShiftCorrPath, coll.timestamp());
1009+
shiftprofile.push_back(objshift);
1010+
}
1011+
lastRunNumber = currentRunNumber;
1012+
}
1013+
}
1014+
if (cfgShiftCorr){
1015+
PsiT0CCorr = ApplyShiftCorrection(coll, PsiT0C);
1016+
}
1017+
1018+
histos.fill(HIST("hPsiT0C"), PsiT0CCorr);
9771019
histos.fill(HIST("hPsiZDCA_vs_ZDCC"), coll.psiZDCC(), coll.psiZDCA());
978-
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0C);
1020+
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0CCorr);
9791021

9801022
resolution.fill(HIST("QVectorsT0CTPCA"), eventplaneVecT0C.Dot(eventplaneVecTPCA), coll.centFT0C());
9811023
resolution.fill(HIST("QVectorsT0CTPCC"), eventplaneVecT0C.Dot(eventplaneVecTPCC), coll.centFT0C());
@@ -1051,7 +1093,7 @@ struct cascadeFlow {
10511093

10521094
ROOT::Math::XYZVector cascQvec{std::cos(2 * casc.phi()), std::sin(2 * casc.phi()), 0};
10531095
auto v2CSP = cascQvec.Dot(eventplaneVecT0C); // not normalised by amplitude
1054-
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0C);
1096+
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0CCorr);
10551097
auto v2CEP = std::cos(2.0 * cascminuspsiT0C);
10561098
ROOT::Math::XYZVector cascUvec{std::cos(casc.phi()), std::sin(casc.phi()), 0};
10571099
auto v1SP_ZDCA = cascUvec.Dot(spectatorplaneVecZDCA);
@@ -1107,11 +1149,11 @@ struct cascadeFlow {
11071149
int chargeIndex = 0;
11081150
if (casc.sign() > 0)
11091151
chargeIndex = 1;
1110-
double pzs2Xi = cosThetaStarLambda[0] * std::sin(2 * (casc.phi() - PsiT0C)) / cascadev2::AlphaXi[chargeIndex] / meanCos2ThetaLambdaFromXi;
1111-
double pzs2Omega = cosThetaStarLambda[1] * std::sin(2 * (casc.phi() - PsiT0C)) / cascadev2::AlphaOmega[chargeIndex] / meanCos2ThetaLambdaFromOmega;
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;
11121154
double cos2ThetaXi = cosThetaStarLambda[0] * cosThetaStarLambda[0];
11131155
double cos2ThetaOmega = cosThetaStarLambda[1] * cosThetaStarLambda[1];
1114-
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - PsiT0C)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
1156+
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
11151157
double cos2ThetaLambda = cosThetaStarProton * cosThetaStarProton;
11161158

11171159
double cosThetaXiWithAlpha = cosThetaStarLambda[0] / cascadev2::AlphaXi[chargeIndex];
@@ -1217,7 +1259,7 @@ struct cascadeFlow {
12171259

12181260
if (isSelectedCasc[0] || isSelectedCasc[1]) {
12191261
if (fillingConfigs.isFillTree)
1220-
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, PsiT0C, BDTresponse[0], BDTresponse[1], 0);
1262+
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, PsiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
12211263
}
12221264
}
12231265
}
@@ -1254,8 +1296,26 @@ struct cascadeFlow {
12541296
ROOT::Math::XYZVector eventplaneVecTPCC{coll.qvecBNegRe(), coll.qvecBNegIm(), 0};
12551297

12561298
const float PsiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1257-
histos.fill(HIST("hPsiT0C"), PsiT0C);
1258-
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0C);
1299+
float PsiT0CCorr = PsiT0C;
1300+
if (cfgShiftCorr) {
1301+
currentRunNumber = coll.runNumber();
1302+
if (currentRunNumber != lastRunNumber) {
1303+
shiftprofile.clear();
1304+
for (int i = 2; i < cfgnMods + 2; i++) {
1305+
fullCCDBShiftCorrPath = cfgShiftPath;
1306+
fullCCDBShiftCorrPath += "/v";
1307+
fullCCDBShiftCorrPath += std::to_string(i);
1308+
auto objshift = ccdb->getForTimeStamp<TProfile3D>(fullCCDBShiftCorrPath, coll.timestamp());
1309+
shiftprofile.push_back(objshift);
1310+
}
1311+
lastRunNumber = currentRunNumber;
1312+
}
1313+
}
1314+
if (cfgShiftCorr){
1315+
PsiT0CCorr = ApplyShiftCorrection(coll, PsiT0C);
1316+
}
1317+
histos.fill(HIST("hPsiT0C"), PsiT0CCorr);
1318+
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0CCorr);
12591319

12601320
resolution.fill(HIST("QVectorsT0CTPCA"), eventplaneVecT0C.Dot(eventplaneVecTPCA), coll.centFT0C());
12611321
resolution.fill(HIST("QVectorsT0CTPCC"), eventplaneVecT0C.Dot(eventplaneVecTPCC), coll.centFT0C());
@@ -1330,7 +1390,7 @@ struct cascadeFlow {
13301390

13311391
ROOT::Math::XYZVector cascQvec{std::cos(2 * casc.phi()), std::sin(2 * casc.phi()), 0};
13321392
auto v2CSP = cascQvec.Dot(eventplaneVecT0C); // not normalised by amplitude
1333-
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0C);
1393+
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0CCorr);
13341394
auto v2CEP = std::cos(2.0 * cascminuspsiT0C);
13351395
ROOT::Math::XYZVector cascUvec{std::cos(casc.phi()), std::sin(casc.phi()), 0};
13361396

@@ -1380,11 +1440,11 @@ struct cascadeFlow {
13801440
int chargeIndex = 0;
13811441
if (casc.sign() > 0)
13821442
chargeIndex = 1;
1383-
double pzs2Xi = cosThetaStarLambda[0] * std::sin(2 * (casc.phi() - PsiT0C)) / cascadev2::AlphaXi[chargeIndex] / meanCos2ThetaLambdaFromXi;
1384-
double pzs2Omega = cosThetaStarLambda[1] * std::sin(2 * (casc.phi() - PsiT0C)) / cascadev2::AlphaOmega[chargeIndex] / meanCos2ThetaLambdaFromOmega;
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;
13851445
double cos2ThetaXi = cosThetaStarLambda[0] * cosThetaStarLambda[0];
13861446
double cos2ThetaOmega = cosThetaStarLambda[1] * cosThetaStarLambda[1];
1387-
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - PsiT0C)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
1447+
double pzs2LambdaFromCasc = cosThetaStarProton * std::sin(2 * (casc.phi() - PsiT0CCorr)) / cascadev2::AlphaLambda[chargeIndex] / meanCos2ThetaProtonFromLambda;
13881448
double cos2ThetaLambda = cosThetaStarProton * cosThetaStarProton;
13891449

13901450
double cosThetaXiWithAlpha = cosThetaStarLambda[0] / cascadev2::AlphaXi[chargeIndex];
@@ -1488,7 +1548,7 @@ struct cascadeFlow {
14881548

14891549
if (isSelectedCasc[0] || isSelectedCasc[1]) {
14901550
if (fillingConfigs.isFillTree)
1491-
fillAnalysedTable(coll, hasEventPlane, 0, casc, v2CSP, v2CEP, 0, 0, PsiT0C, BDTresponse[0], BDTresponse[1], 0);
1551+
fillAnalysedTable(coll, hasEventPlane, 0, casc, v2CSP, v2CEP, 0, 0, PsiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
14921552
}
14931553
}
14941554
}
@@ -1523,8 +1583,26 @@ struct cascadeFlow {
15231583
ROOT::Math::XYZVector eventplaneVecTPCC{coll.qvecBNegRe(), coll.qvecBNegIm(), 0};
15241584

15251585
const float psiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1526-
histos.fill(HIST("hPsiT0C"), psiT0C);
1527-
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), psiT0C);
1586+
float psiT0CCorr = psiT0C;
1587+
if (cfgShiftCorr) {
1588+
currentRunNumber = coll.runNumber();
1589+
if (currentRunNumber != lastRunNumber) {
1590+
shiftprofile.clear();
1591+
for (int i = 2; i < cfgnMods + 2; i++) {
1592+
fullCCDBShiftCorrPath = cfgShiftPath;
1593+
fullCCDBShiftCorrPath += "/v";
1594+
fullCCDBShiftCorrPath += std::to_string(i);
1595+
auto objshift = ccdb->getForTimeStamp<TProfile3D>(fullCCDBShiftCorrPath, coll.timestamp());
1596+
shiftprofile.push_back(objshift);
1597+
}
1598+
lastRunNumber = currentRunNumber;
1599+
}
1600+
}
1601+
if (cfgShiftCorr){
1602+
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
1603+
}
1604+
histos.fill(HIST("hPsiT0C"), psiT0CCorr);
1605+
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), psiT0CCorr);
15281606

15291607
resolution.fill(HIST("QVectorsT0CTPCA"), eventplaneVecT0C.Dot(eventplaneVecTPCA), coll.centFT0C());
15301608
resolution.fill(HIST("QVectorsT0CTPCC"), eventplaneVecT0C.Dot(eventplaneVecTPCC), coll.centFT0C());
@@ -1595,7 +1673,7 @@ struct cascadeFlow {
15951673

15961674
ROOT::Math::XYZVector lambdaQvec{std::cos(2 * v0.phi()), std::sin(2 * v0.phi()), 0};
15971675
auto v2CSP = lambdaQvec.Dot(eventplaneVecT0C); // not normalised by amplitude
1598-
auto lambdaminuspsiT0C = GetPhiInRange(v0.phi() - psiT0C);
1676+
auto lambdaminuspsiT0C = GetPhiInRange(v0.phi() - psiT0CCorr);
15991677
auto v2CEP = std::cos(2.0 * lambdaminuspsiT0C);
16001678
ROOT::Math::XYZVector lambdaUvec{std::cos(v0.phi()), std::sin(v0.phi()), 0};
16011679

@@ -1625,15 +1703,15 @@ struct cascadeFlow {
16251703
double cos2ThetaLambda = 0;
16261704
double cosThetaLambda = 0;
16271705
if (chargeIndex == 0) {
1628-
pzs2Lambda = cosThetaStarProton[0] * std::sin(2 * (v0.phi() - psiT0C)) / lambdav2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
1706+
pzs2Lambda = cosThetaStarProton[0] * std::sin(2 * (v0.phi() - psiT0CCorr)) / lambdav2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
16291707
cos2ThetaLambda = cosThetaStarProton[0] * cosThetaStarProton[0];
16301708
cosThetaLambda = cosThetaStarProton[0] / cascadev2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
16311709
} else if (chargeIndex == 1) {
1632-
pzs2Lambda = cosThetaStarProton[1] * std::sin(2 * (v0.phi() - psiT0C)) / lambdav2::AlphaLambda[1] / meanCos2ThetaProtonFromLambda;
1710+
pzs2Lambda = cosThetaStarProton[1] * std::sin(2 * (v0.phi() - psiT0CCorr)) / lambdav2::AlphaLambda[1] / meanCos2ThetaProtonFromLambda;
16331711
cos2ThetaLambda = cosThetaStarProton[1] * cosThetaStarProton[1];
16341712
cosThetaLambda = cosThetaStarProton[1] / cascadev2::AlphaLambda[1] / meanCos2ThetaProtonFromLambda;
16351713
} else { // I treat these bkg candidates as Lambdas for the purpose of calculating Pz
1636-
pzs2Lambda = cosThetaStarProton[0] * std::sin(2 * (v0.phi() - psiT0C)) / lambdav2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
1714+
pzs2Lambda = cosThetaStarProton[0] * std::sin(2 * (v0.phi() - psiT0CCorr)) / lambdav2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
16371715
cos2ThetaLambda = cosThetaStarProton[0] * cosThetaStarProton[0];
16381716
cosThetaLambda = cosThetaStarProton[0] / cascadev2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
16391717
}
@@ -1660,7 +1738,7 @@ struct cascadeFlow {
16601738
}
16611739

16621740
if (fillingConfigs.isFillTree)
1663-
fillAnalysedLambdaTable(coll, hasEventPlane, hasSpectatorPlane, chargeIndex, v0, v2CEP, psiT0C, pzs2Lambda, cos2ThetaLambda, cosThetaLambda);
1741+
fillAnalysedLambdaTable(coll, hasEventPlane, hasSpectatorPlane, chargeIndex, v0, v2CEP, psiT0CCorr, pzs2Lambda, cos2ThetaLambda, cosThetaLambda);
16641742
}
16651743
}
16661744

@@ -1707,8 +1785,26 @@ struct cascadeFlow {
17071785
float NormQvTPCC = std::sqrt(eventplaneVecTPCC.Dot(eventplaneVecTPCC));
17081786

17091787
const float PsiT0C = std::atan2(coll.qvecFT0CIm(), coll.qvecFT0CRe()) * 0.5f;
1710-
histos.fill(HIST("hPsiT0C"), PsiT0C);
1711-
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0C);
1788+
float PsiT0CCorr = PsiT0C;
1789+
if (cfgShiftCorr) {
1790+
currentRunNumber = coll.runNumber();
1791+
if (currentRunNumber != lastRunNumber) {
1792+
shiftprofile.clear();
1793+
for (int i = 2; i < cfgnMods + 2; i++) {
1794+
fullCCDBShiftCorrPath = cfgShiftPath;
1795+
fullCCDBShiftCorrPath += "/v";
1796+
fullCCDBShiftCorrPath += std::to_string(i);
1797+
auto objshift = ccdb->getForTimeStamp<TProfile3D>(fullCCDBShiftCorrPath, coll.timestamp());
1798+
shiftprofile.push_back(objshift);
1799+
}
1800+
lastRunNumber = currentRunNumber;
1801+
}
1802+
}
1803+
if (cfgShiftCorr){
1804+
PsiT0CCorr = ApplyShiftCorrection(coll, PsiT0C);
1805+
}
1806+
histos.fill(HIST("hPsiT0C"), PsiT0CCorr);
1807+
histos.fill(HIST("hPsiT0CvsCentFT0C"), coll.centFT0C(), PsiT0CCorr);
17121808

17131809
resolution.fill(HIST("QVectorsT0CTPCA"), eventplaneVecT0C.Dot(eventplaneVecTPCA), coll.centFT0C());
17141810
resolution.fill(HIST("QVectorsT0CTPCC"), eventplaneVecT0C.Dot(eventplaneVecTPCC), coll.centFT0C());
@@ -1784,7 +1880,7 @@ struct cascadeFlow {
17841880

17851881
ROOT::Math::XYZVector cascQvec{std::cos(2 * casc.phi()), std::sin(2 * casc.phi()), 0};
17861882
auto v2CSP = cascQvec.Dot(eventplaneVecT0C);
1787-
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0C);
1883+
auto cascminuspsiT0C = GetPhiInRange(casc.phi() - PsiT0CCorr);
17881884
auto v2CEP = std::cos(2.0 * cascminuspsiT0C);
17891885
ROOT::Math::XYZVector cascUvec{std::cos(casc.phi()), std::sin(casc.phi()), 0};
17901886
auto v1SP_ZDCA = cascUvec.Dot(spectatorplaneVecZDCA);
@@ -1816,7 +1912,7 @@ struct cascadeFlow {
18161912
}
18171913
if (isSelectedCasc[0] || isSelectedCasc[1])
18181914
if (fillingConfigs.isFillTree)
1819-
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, PsiT0C, BDTresponse[0], BDTresponse[1], 0);
1915+
fillAnalysedTable(coll, hasEventPlane, hasSpectatorPlane, casc, v2CSP, v2CEP, v1SP_ZDCA, v1SP_ZDCC, PsiT0CCorr, BDTresponse[0], BDTresponse[1], 0);
18201916
}
18211917
}
18221918

0 commit comments

Comments
 (0)