diff --git a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx index b42396d0457..f1a866e6191 100644 --- a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx +++ b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx @@ -191,7 +191,6 @@ struct lambdapolsp { Configurable ConfNUA{"ConfNUA", "Users/p/prottay/My/Object/NUAwgtschk", "Path to NUA"}; Configurable fillnominal{"fillnominal", true, "flag to fill default histograms"}; - Configurable filldefwithsyst{"filldefwithsyst", true, "flag to fill default histograms in syst"}; } QAgrp; struct : ConfigurableGroup { @@ -500,7 +499,18 @@ struct lambdapolsp { histos.add("hpQxytpvscent", "hpQxytpvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true); histos.add("hpQxtQypvscent", "hpQxtQypvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true); histos.add("hpQxpQytvscent", "hpQxpQytvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true); + + histos.add("hSparseLambdaPolSP_xAwgt", "hSparseLambdaPolSP_xAwgt", HistType::kTHnSparseF, runaxes, true); + histos.add("hSparseLambdaPolSP_yAwgt", "hSparseLambdaPolSP_yAwgt", HistType::kTHnSparseF, runaxes, true); + histos.add("hSparseLambdaPolSP_xCwgt", "hSparseLambdaPolSP_xCwgt", HistType::kTHnSparseF, runaxes, true); + histos.add("hSparseLambdaPolSP_yCwgt", "hSparseLambdaPolSP_yCwgt", HistType::kTHnSparseF, runaxes, true); + + histos.add("hSparseAntiLambdaPolSP_xAwgt", "hSparseAntiLambdaPolSP_xAwgt", HistType::kTHnSparseF, runaxes, true); + histos.add("hSparseAntiLambdaPolSP_yAwgt", "hSparseAntiLambdaPolSP_yAwgt", HistType::kTHnSparseF, runaxes, true); + histos.add("hSparseAntiLambdaPolSP_xCwgt", "hSparseAntiLambdaPolSP_xCwgt", HistType::kTHnSparseF, runaxes, true); + histos.add("hSparseAntiLambdaPolSP_yCwgt", "hSparseAntiLambdaPolSP_yCwgt", HistType::kTHnSparseF, runaxes, true); } + if (usesubdet) { histos.add("hSparseLambdaCosPsiA", "hSparseLambdaCosPsiA", HistType::kTHnSparseF, runaxes, true); histos.add("hSparseLambdaSinPsiA", "hSparseLambdaSinPsiA", HistType::kTHnSparseF, runaxes, true); @@ -545,12 +555,12 @@ struct lambdapolsp { // histos.add("hSparseLambda_corr2b", "hSparseLambda_corr2b", HistType::kTHnSparseF, runaxes, true); histos.add("hSparseAntiLambda_corr2a", "hSparseAntiLambda_corr2a", HistType::kTHnSparseF, runaxes, true); // histos.add("hSparseAntiLambda_corr2b", "hSparseAntiLambda_corr2b", HistType::kTHnSparseF, runaxes, true); - if (randGrp.useSP) { + /*if (randGrp.useSP) { histos.add("hSparseAntiLambda_avgux", "hSparseAntiLambda_avgux", HistType::kTHnSparseF, {thnAxisInvMass, axisGrp.configthnAxispT, axisGrp.configthnAxisPol, axisGrp.configcentAxis}, true); histos.add("hSparseAntiLambda_avguy", "hSparseAntiLambda_avguy", HistType::kTHnSparseF, {thnAxisInvMass, axisGrp.configthnAxispT, axisGrp.configthnAxisPol, axisGrp.configcentAxis}, true); histos.add("hSparseLambda_avgux", "hSparseLambda_avgux", HistType::kTHnSparseF, {thnAxisInvMass, axisGrp.configthnAxispT, axisGrp.configthnAxisPol, axisGrp.configcentAxis}, true); histos.add("hSparseLambda_avguy", "hSparseLambda_avguy", HistType::kTHnSparseF, {thnAxisInvMass, axisGrp.configthnAxispT, axisGrp.configthnAxisPol, axisGrp.configcentAxis}, true); - } + }*/ } if (distGrp.filldist) { @@ -720,6 +730,14 @@ struct lambdapolsp { HistType::kTHnSparseF, runaxesSyst, true); + + histos.add("hSparseAntiLambda_corr1aSyst", "hSparseAntiLambda_corr1aSyst", HistType::kTHnSparseF, runaxesSyst, true); + histos.add("hSparseAntiLambda_corr1bSyst", "hSparseAntiLambda_corr1bSyst", HistType::kTHnSparseF, runaxesSyst, true); + histos.add("hSparseAntiLambda_corr2aSyst", "hSparseAntiLambda_corr2aSyst", HistType::kTHnSparseF, runaxesSyst, true); + + histos.add("hSparseLambda_corr1aSyst", "hSparseLambda_corr1aSyst", HistType::kTHnSparseF, runaxesSyst, true); + histos.add("hSparseLambda_corr1bSyst", "hSparseLambda_corr1bSyst", HistType::kTHnSparseF, runaxesSyst, true); + histos.add("hSparseLambda_corr2aSyst", "hSparseLambda_corr2aSyst", HistType::kTHnSparseF, runaxesSyst, true); } } @@ -1080,6 +1098,20 @@ struct lambdapolsp { auto PolSP_A = uy * modqxZDCA - ux * modqyZDCA; // u_y QxA - u_x QyA auto PolSP_C = uy * modqxZDCC - ux * modqyZDCC; // u_y QxC - u_x QyC + // SP numerator components separately + auto PolSP_xA = -ux * modqyZDCA; + auto PolSP_yA = uy * modqxZDCA; + + auto PolSP_xC = -ux * modqyZDCC; + auto PolSP_yC = uy * modqxZDCC; + + // acceptance-corrected versions + auto PolSP_xAwgt = PolSP_xA / acvalue; + auto PolSP_yAwgt = PolSP_yA / acvalue; + + auto PolSP_xCwgt = PolSP_xC / acvalue; + auto PolSP_yCwgt = PolSP_yC / acvalue; + if (randGrp.useSP) { Pol = PolSP; PolA = PolSP_A; @@ -1088,6 +1120,7 @@ struct lambdapolsp { PolAwgt = PolSP_A / acvalue; PolCwgt = PolSP_C / acvalue; } + ////////////////////////////// if (QAgrp.fillnominal) { @@ -1116,9 +1149,28 @@ struct lambdapolsp { histos.fill(HIST("hSparseAntiLambda_corr2a"), candmass, candpt, sinThetaStar, centrality, desbinvalue, wgtfactor); // histos.fill(HIST("hSparseAntiLambda_corr2b"), candmass, candpt, sinThetaStarcosphiphiStar, centrality, desbinvalue, wgtfactor); if (randGrp.useSP) { - histos.fill(HIST("hSparseAntiLambda_avgux"), candmass, candpt, ux, centrality); - histos.fill(HIST("hSparseAntiLambda_avguy"), candmass, candpt, uy, centrality); + // histos.fill(HIST("hSparseAntiLambda_avgux"), candmass, candpt, ux, centrality); + // histos.fill(HIST("hSparseAntiLambda_avguy"), candmass, candpt, uy, centrality); + + histos.fill(HIST("hSparseAntiLambdaPolwgt"), candmass, candpt, Polwgt, centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_xAwgt"), + candmass, candpt, PolSP_xAwgt, + centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_yAwgt"), + candmass, candpt, PolSP_yAwgt, + centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_xCwgt"), + candmass, candpt, PolSP_xCwgt, + centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_yCwgt"), + candmass, candpt, PolSP_yCwgt, + centrality, desbinvalue, wgtfactor); } + } else { if (usesubdet) { histos.fill(HIST("hSparseAntiLambdaCosPsiA"), candmass, candpt, (TMath::Cos(GetPhiInRange(psiZDCA))), centrality, wgtfactor); @@ -1142,8 +1194,26 @@ struct lambdapolsp { histos.fill(HIST("hSparseAntiLambda_corr2a"), candmass, candpt, sinThetaStar, centrality, wgtfactor); // histos.fill(HIST("hSparseAntiLambda_corr2b"), candmass, candpt, sinThetaStarcosphiphiStar, centrality, wgtfactor); if (randGrp.useSP) { - histos.fill(HIST("hSparseAntiLambda_avgux"), candmass, candpt, ux, centrality); - histos.fill(HIST("hSparseAntiLambda_avguy"), candmass, candpt, uy, centrality); + // histos.fill(HIST("hSparseAntiLambda_avgux"), candmass, candpt, ux, centrality); + // histos.fill(HIST("hSparseAntiLambda_avguy"), candmass, candpt, uy, centrality); + + histos.fill(HIST("hSparseAntiLambdaPolwgt"), candmass, candpt, Polwgt, centrality, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_xAwgt"), + candmass, candpt, PolSP_xAwgt, + centrality, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_yAwgt"), + candmass, candpt, PolSP_yAwgt, + centrality, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_xCwgt"), + candmass, candpt, PolSP_xCwgt, + centrality, wgtfactor); + + histos.fill(HIST("hSparseAntiLambdaPolSP_yCwgt"), + candmass, candpt, PolSP_yCwgt, + centrality, wgtfactor); } } } @@ -1171,8 +1241,26 @@ struct lambdapolsp { histos.fill(HIST("hSparseLambda_corr2a"), candmass, candpt, sinThetaStar, centrality, desbinvalue, wgtfactor); // histos.fill(HIST("hSparseLambda_corr2b"), candmass, candpt, sinThetaStarcosphiphiStar, centrality, desbinvalue, wgtfactor); if (randGrp.useSP) { - histos.fill(HIST("hSparseLambda_avgux"), candmass, candpt, ux, centrality); - histos.fill(HIST("hSparseLambda_avguy"), candmass, candpt, uy, centrality); + // histos.fill(HIST("hSparseLambda_avgux"), candmass, candpt, ux, centrality); + // histos.fill(HIST("hSparseLambda_avguy"), candmass, candpt, uy, centrality); + + histos.fill(HIST("hSparseLambdaPolwgt"), candmass, candpt, Polwgt, centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_xAwgt"), + candmass, candpt, PolSP_xAwgt, + centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_yAwgt"), + candmass, candpt, PolSP_yAwgt, + centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_xCwgt"), + candmass, candpt, PolSP_xCwgt, + centrality, desbinvalue, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_yCwgt"), + candmass, candpt, PolSP_yCwgt, + centrality, desbinvalue, wgtfactor); } } else { if (usesubdet) { @@ -1197,8 +1285,26 @@ struct lambdapolsp { histos.fill(HIST("hSparseLambda_corr2a"), candmass, candpt, sinThetaStar, centrality, wgtfactor); // histos.fill(HIST("hSparseLambda_corr2b"), candmass, candpt, sinThetaStarcosphiphiStar, centrality, wgtfactor); if (randGrp.useSP) { - histos.fill(HIST("hSparseLambda_avgux"), candmass, candpt, ux, centrality); - histos.fill(HIST("hSparseLambda_avguy"), candmass, candpt, uy, centrality); + // histos.fill(HIST("hSparseLambda_avgux"), candmass, candpt, ux, centrality); + // histos.fill(HIST("hSparseLambda_avguy"), candmass, candpt, uy, centrality); + + histos.fill(HIST("hSparseLambdaPolwgt"), candmass, candpt, Polwgt, centrality, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_xAwgt"), + candmass, candpt, PolSP_xAwgt, + centrality, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_yAwgt"), + candmass, candpt, PolSP_yAwgt, + centrality, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_xCwgt"), + candmass, candpt, PolSP_xCwgt, + centrality, wgtfactor); + + histos.fill(HIST("hSparseLambdaPolSP_yCwgt"), + candmass, candpt, PolSP_yCwgt, + centrality, wgtfactor); } } } @@ -1223,18 +1329,9 @@ struct lambdapolsp { static_cast(isyst), wgtfactor); - histos.fill(HIST("hSparseLambdaPolwgtSyst"), - candmass, - candpt, - Polwgt, - centrality, - desbinvalue, - static_cast(isyst), - wgtfactor); - if (QAgrp.filldefwithsyst) { - histos.fill(HIST("hSparseLambdaPol"), candmass, candpt, Pol, centrality, desbinvalue, wgtfactor); - histos.fill(HIST("hSparseLambdaPolwgt"), candmass, candpt, Polwgt, centrality, desbinvalue, wgtfactor); - } + histos.fill(HIST("hSparseLambda_corr1aSyst"), candmass, candpt, sinPhiStar, centrality, desbinvalue, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseLambda_corr1bSyst"), candmass, candpt, cosPhiStar, centrality, desbinvalue, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseLambda_corr2aSyst"), candmass, candpt, sinThetaStar, centrality, desbinvalue, static_cast(isyst), wgtfactor); } else { @@ -1246,18 +1343,9 @@ struct lambdapolsp { static_cast(isyst), wgtfactor); - histos.fill(HIST("hSparseLambdaPolwgtSyst"), - candmass, - candpt, - Polwgt, - centrality, - static_cast(isyst), - wgtfactor); - - if (QAgrp.filldefwithsyst) { - histos.fill(HIST("hSparseLambdaPol"), candmass, candpt, Pol, centrality, wgtfactor); - histos.fill(HIST("hSparseLambdaPolwgt"), candmass, candpt, Polwgt, centrality, wgtfactor); - } + histos.fill(HIST("hSparseLambda_corr1aSyst"), candmass, candpt, sinPhiStar, centrality, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseLambda_corr1bSyst"), candmass, candpt, cosPhiStar, centrality, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseLambda_corr2aSyst"), candmass, candpt, sinThetaStar, centrality, static_cast(isyst), wgtfactor); } } @@ -1275,18 +1363,9 @@ struct lambdapolsp { static_cast(isyst), wgtfactor); - histos.fill(HIST("hSparseAntiLambdaPolwgtSyst"), - candmass, - candpt, - Polwgt, - centrality, - desbinvalue, - static_cast(isyst), - wgtfactor); - if (QAgrp.filldefwithsyst) { - histos.fill(HIST("hSparseAntiLambdaPol"), candmass, candpt, Pol, centrality, desbinvalue, wgtfactor); - histos.fill(HIST("hSparseAntiLambdaPolwgt"), candmass, candpt, Polwgt, centrality, desbinvalue, wgtfactor); - } + histos.fill(HIST("hSparseAntiLambda_corr1aSyst"), candmass, candpt, sinPhiStar, centrality, desbinvalue, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseAntiLambda_corr1bSyst"), candmass, candpt, cosPhiStar, centrality, desbinvalue, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseAntiLambda_corr2aSyst"), candmass, candpt, sinThetaStar, centrality, desbinvalue, static_cast(isyst), wgtfactor); } else { @@ -1298,18 +1377,9 @@ struct lambdapolsp { static_cast(isyst), wgtfactor); - histos.fill(HIST("hSparseAntiLambdaPolwgtSyst"), - candmass, - candpt, - Polwgt, - centrality, - static_cast(isyst), - wgtfactor); - - if (QAgrp.filldefwithsyst) { - histos.fill(HIST("hSparseAntiLambdaPol"), candmass, candpt, Pol, centrality, wgtfactor); - histos.fill(HIST("hSparseAntiLambdaPolwgt"), candmass, candpt, Polwgt, centrality, wgtfactor); - } + histos.fill(HIST("hSparseAntiLambda_corr1aSyst"), candmass, candpt, sinPhiStar, centrality, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseAntiLambda_corr1bSyst"), candmass, candpt, cosPhiStar, centrality, static_cast(isyst), wgtfactor); + histos.fill(HIST("hSparseAntiLambda_corr2aSyst"), candmass, candpt, sinThetaStar, centrality, static_cast(isyst), wgtfactor); } } }