Skip to content

Commit 506729e

Browse files
Merge pull request #33 from alibuild/alibot-cleanup-12630
Please consider the following formatting changes to #12630
2 parents f2594ed + 026d204 commit 506729e

File tree

1 file changed

+14
-12
lines changed

1 file changed

+14
-12
lines changed

PWGLF/TableProducer/Strangeness/cascadeflow.cxx

Lines changed: 14 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -174,7 +174,7 @@ struct cascadeFlow {
174174
Configurable<bool> cfgShiftCorr{"cfgShiftCorr", 0, ""};
175175
Configurable<std::string> cfgShiftPath{"cfgShiftPath", "Users/j/junlee/Qvector/QvecCalib/Shift", "Path for Shift"};
176176
Configurable<int> cfgnMods{"cfgnMods", 1, "The number of modulations of interest starting from 2"};
177-
177+
178178
// THN axes
179179
ConfigurableAxis thnConfigAxisFT0C{"thnConfigAxisFT0C", {8, 0, 80}, "FT0C centrality (%)"};
180180
ConfigurableAxis thnConfigAxisEta{"thnConfigAxisEta", {8, -0.8, 0.8}, "pseudorapidity"};
@@ -473,21 +473,23 @@ struct cascadeFlow {
473473
std::string fullCCDBShiftCorrPath;
474474

475475
template <typename TCollision>
476-
double ApplyShiftCorrection(TCollision coll, double psiT0C){
476+
double ApplyShiftCorrection(TCollision coll, double psiT0C)
477+
{
477478
int nmode = 2;
478479
auto deltapsiFT0C = 0.0;
479480

480481
for (int ishift = 1; ishift <= 10; ishift++) {
481482
auto coeffshiftxFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 0.5, ishift - 0.5));
482483
auto coeffshiftyFT0C = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 1.5, ishift - 0.5));
483-
484+
484485
deltapsiFT0C += ((1 / (1.0 * ishift)) * (-coeffshiftxFT0C * TMath::Cos(ishift * static_cast<float>(nmode) * psiT0C) + coeffshiftyFT0C * TMath::Sin(ishift * static_cast<float>(nmode) * psiT0C)));
485486
}
486487
return psiT0C + deltapsiFT0C;
487488
}
488489

489490
template <typename TCollision>
490-
double ComputeEPResolutionwShifts(TCollision coll, double psiT0C, double psiTPCA, double psiTPCC){
491+
double ComputeEPResolutionwShifts(TCollision coll, double psiT0C, double psiTPCA, double psiTPCC)
492+
{
491493
int nmode = 2;
492494
auto deltapsiFT0C = 0.0;
493495
auto deltapsiTPCA = 0.0;
@@ -498,14 +500,14 @@ struct cascadeFlow {
498500
auto coeffshiftxTPCA = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 2.5, ishift - 0.5));
499501
auto coeffshiftyTPCA = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 3.5, ishift - 0.5));
500502
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));
503+
auto coeffshiftyTPCC = shiftprofile.at(nmode - 2)->GetBinContent(shiftprofile.at(nmode - 2)->FindBin(coll.centFT0C(), 5.5, ishift - 0.5));
502504
deltapsiFT0C += ((1 / (1.0 * ishift)) * (-coeffshiftxFT0C * TMath::Cos(ishift * static_cast<float>(nmode) * psiT0C) + coeffshiftyFT0C * TMath::Sin(ishift * static_cast<float>(nmode) * psiT0C)));
503505
deltapsiTPCA += ((1 / (1.0 * ishift)) * (-coeffshiftxTPCA * TMath::Cos(ishift * static_cast<float>(nmode) * psiTPCA) + coeffshiftyTPCA * TMath::Sin(ishift * static_cast<float>(nmode) * psiTPCA)));
504506
deltapsiTPCC += ((1 / (1.0 * ishift)) * (-coeffshiftxTPCC * TMath::Cos(ishift * static_cast<float>(nmode) * psiTPCC) + coeffshiftyTPCC * TMath::Sin(ishift * static_cast<float>(nmode) * psiTPCC)));
505507
}
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);
508+
// histos.fill(HIST("psi2/QA/EP_FT0C_shifted"), coll.centFT0C(), psiT0C + deltapsiFT0C);
509+
// histos.fill(HIST("psi2/QA/EP_TPCA_shifted"), coll.centFT0C(), psiTPCA + deltapsiTPCA);
510+
// histos.fill(HIST("psi2/QA/EP_TPCC_shifted"), coll.centFT0C(), psiTPCC + deltapsiTPCC);
509511
resolution.fill(HIST("QVectorsT0CTPCA_Shifted"), coll.centFT0C(), TMath::Cos(static_cast<float>(nmode) * (psiT0C + deltapsiFT0C - psiTPCA - deltapsiTPCA)));
510512
resolution.fill(HIST("QVectorsT0CTPCC_Shifted"), coll.centFT0C(), TMath::Cos(static_cast<float>(nmode) * (psiT0C + deltapsiFT0C - psiTPCC - deltapsiTPCC)));
511513
resolution.fill(HIST("QVectorsTPCAC_Shifted"), coll.centFT0C(), TMath::Cos(static_cast<float>(nmode) * (psiTPCA + deltapsiTPCA - psiTPCC - deltapsiTPCC)));
@@ -1043,7 +1045,7 @@ struct cascadeFlow {
10431045
lastRunNumber = currentRunNumber;
10441046
}
10451047
}
1046-
if (cfgShiftCorr){
1048+
if (cfgShiftCorr) {
10471049
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
10481050
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
10491051
}
@@ -1346,7 +1348,7 @@ struct cascadeFlow {
13461348
lastRunNumber = currentRunNumber;
13471349
}
13481350
}
1349-
if (cfgShiftCorr){
1351+
if (cfgShiftCorr) {
13501352
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
13511353
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
13521354
}
@@ -1636,7 +1638,7 @@ struct cascadeFlow {
16361638
lastRunNumber = currentRunNumber;
16371639
}
16381640
}
1639-
if (cfgShiftCorr){
1641+
if (cfgShiftCorr) {
16401642
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
16411643
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
16421644
}
@@ -1841,7 +1843,7 @@ struct cascadeFlow {
18411843
lastRunNumber = currentRunNumber;
18421844
}
18431845
}
1844-
if (cfgShiftCorr){
1846+
if (cfgShiftCorr) {
18451847
psiT0CCorr = ApplyShiftCorrection(coll, psiT0C);
18461848
ComputeEPResolutionwShifts(coll, psiT0C, psiTPCA, psiTPCC);
18471849
}

0 commit comments

Comments
 (0)