@@ -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