Skip to content

Commit dec23a2

Browse files
author
Prottay Das
committed
added few QA histograms
1 parent 0688089 commit dec23a2

File tree

1 file changed

+96
-0
lines changed

1 file changed

+96
-0
lines changed

‎PWGLF/Tasks/Strangeness/lambdapolsp.cxx‎

Lines changed: 96 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -190,6 +190,16 @@ struct lambdapolsp {
190190
ConfigurableAxis configbinAxis{"configbinAxis", {VARIABLE_WIDTH, -0.8, -0.4, -0.2, 0, 0.2, 0.4, 0.8}, "BA"};
191191
*/
192192
// ConfigurableAxis configphiAxis{"configphiAxis", {VARIABLE_WIDTH, 0.0, 0.2, 0.4, 0.8, 1.0, 2.0, 2.5, 3.0, 4.0, 5.0, 5.5, 6.28}, "PhiAxis"};
193+
194+
struct : ConfigurableGroup {
195+
Configurable<bool> isQA{"isQA", true, "Flag to fill QA"};
196+
ConfigurableAxis centfineAxis{"centfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "Cent V0M"};
197+
ConfigurableAxis vxfineAxis{"vxfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "vx fine axis"};
198+
ConfigurableAxis vyfineAxis{"vyfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "vy fine axis"};
199+
ConfigurableAxis vzfineAxis{"vzfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "vz fine axis"};
200+
ConfigurableAxis qxZDCAxis{"qxZDCAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "qx axis"};
201+
ConfigurableAxis psiAxis{"psiAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "psi axis"};
202+
} QAgrp;
193203
struct : ConfigurableGroup {
194204
Configurable<bool> requireRCTFlagChecker{"requireRCTFlagChecker", true, "Check event quality in run condition table"};
195205
Configurable<std::string> cfgEvtRCTFlagCheckerLabel{"cfgEvtRCTFlagCheckerLabel", "CBT_hadronPID", "Evt sel: RCT flag checker label"};
@@ -456,6 +466,32 @@ struct lambdapolsp {
456466
histos.add("hSparseRecLambda", "hSparseRecLambda", HistType::kTHnSparseF, runaxes2, true);
457467
histos.add("hSparseRecAntiLambda", "hSparseRecAntiLambda", HistType::kTHnSparseF, runaxes2, true);
458468

469+
if (QAgrp.isQA) {
470+
histos.add("hCentQxZDCA", "hCentQxZDCA", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});
471+
histos.add("hCentQyZDCA", "hCentQyZDCA", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});
472+
histos.add("hCentQxZDCC", "hCentQxZDCC", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});
473+
histos.add("hCentQyZDCC", "hCentQyZDCC", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});
474+
475+
histos.add("hvxQxZDCA", "hvxQxZDCA", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});
476+
histos.add("hvxQyZDCA", "hvxQyZDCA", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});
477+
histos.add("hvxQxZDCC", "hvxQxZDCC", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});
478+
histos.add("hvxQyZDCC", "hvxQyZDCC", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});
479+
480+
histos.add("hvyQxZDCA", "hvyQxZDCA", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});
481+
histos.add("hvyQyZDCA", "hvyQyZDCA", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});
482+
histos.add("hvyQxZDCC", "hvyQxZDCC", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});
483+
histos.add("hvyQyZDCC", "hvyQyZDCC", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});
484+
485+
histos.add("hvzQxZDCA", "hvzQxZDCA", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});
486+
histos.add("hvzQyZDCA", "hvzQyZDCA", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});
487+
histos.add("hvzQxZDCC", "hvzQxZDCC", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});
488+
histos.add("hvzQyZDCC", "hvzQyZDCC", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});
489+
490+
histos.add("PsiZDCC", "PsiZDCC", kTH2F, {QAgrp.centfineAxis, QAgrp.psiAxis});
491+
histos.add("PsiZDCA", "PsiZDCA", kTH2F, {QAgrp.centfineAxis, QAgrp.psiAxis});
492+
histos.add("PsiZDC", "PsiZDC", kTH2F, {QAgrp.centfineAxis, QAgrp.psiAxis});
493+
}
494+
459495
ccdb->setURL(cfgCcdbParam.cfgURL);
460496
ccdbApi.init("http://alice-ccdb.cern.ch");
461497
ccdb->setCaching(true);
@@ -939,6 +975,10 @@ struct lambdapolsp {
939975
// currentRunNumber = collision.foundBC_as<BCsRun3>().runNumber();
940976
auto bc = collision.foundBC_as<BCsRun3>();
941977

978+
auto vz = collision.vz();
979+
auto vx = collision.vx();
980+
auto vy = collision.vy();
981+
942982
auto qxZDCA = collision.qxZDCA();
943983
auto qxZDCC = collision.qxZDCC();
944984
auto qyZDCA = collision.qyZDCA();
@@ -963,6 +1003,32 @@ struct lambdapolsp {
9631003
psiZDC = psiZDCC - psiZDCA;
9641004
}*/
9651005

1006+
if (QAgrp.isQA) {
1007+
histos.fill(HIST("hCentQxZDCA"), centrality, modqxZDCA);
1008+
histos.fill(HIST("hCentQyZDCA"), centrality, modqyZDCA);
1009+
histos.fill(HIST("hCentQxZDCC"), centrality, modqxZDCC);
1010+
histos.fill(HIST("hCentQyZDCC"), centrality, modqyZDCC);
1011+
1012+
histos.fill(HIST("hvxQxZDCA"), vx, modqxZDCA);
1013+
histos.fill(HIST("hvxQyZDCA"), vx, modqyZDCA);
1014+
histos.fill(HIST("hvxQxZDCC"), vx, modqxZDCC);
1015+
histos.fill(HIST("hvxQyZDCC"), vx, modqyZDCC);
1016+
1017+
histos.fill(HIST("hvyQxZDCA"), vy, modqxZDCA);
1018+
histos.fill(HIST("hvyQyZDCA"), vy, modqyZDCA);
1019+
histos.fill(HIST("hvyQxZDCC"), vy, modqxZDCC);
1020+
histos.fill(HIST("hvyQyZDCC"), vy, modqyZDCC);
1021+
1022+
histos.fill(HIST("hvzQxZDCA"), vz, modqxZDCA);
1023+
histos.fill(HIST("hvzQyZDCA"), vz, modqyZDCA);
1024+
histos.fill(HIST("hvzQxZDCC"), vz, modqxZDCC);
1025+
histos.fill(HIST("hvzQyZDCC"), vz, modqyZDCC);
1026+
1027+
histos.fill(HIST("PsiZDCA"), centrality, psiZDCA);
1028+
histos.fill(HIST("PsiZDCC"), centrality, psiZDCC);
1029+
histos.fill(HIST("PsiZDC"), centrality, psiZDC);
1030+
}
1031+
9661032
histos.fill(HIST("hCentrality"), centrality);
9671033
if (!checkwithpub) {
9681034
// histos.fill(HIST("hVtxZ"), collision.posZ());
@@ -1260,6 +1326,10 @@ struct lambdapolsp {
12601326
else if (centestim == 3)
12611327
centrality = collision.centFV0A();
12621328

1329+
auto vz = collision.posZ();
1330+
auto vx = collision.posX();
1331+
auto vy = collision.posY();
1332+
12631333
auto runnumber = collision.runNumber();
12641334
if (!collision.triggereventsp()) { // provided by StraZDCSP
12651335
return;
@@ -1319,6 +1389,32 @@ struct lambdapolsp {
13191389

13201390
auto psiZDC = TMath::ATan2((modqyZDCC - modqyZDCA), (modqxZDCC - modqxZDCA)); // full event plane
13211391

1392+
if (QAgrp.isQA) {
1393+
histos.fill(HIST("hCentQxZDCA"), centrality, modqxZDCA);
1394+
histos.fill(HIST("hCentQyZDCA"), centrality, modqyZDCA);
1395+
histos.fill(HIST("hCentQxZDCC"), centrality, modqxZDCC);
1396+
histos.fill(HIST("hCentQyZDCC"), centrality, modqyZDCC);
1397+
1398+
histos.fill(HIST("hvxQxZDCA"), vx, modqxZDCA);
1399+
histos.fill(HIST("hvxQyZDCA"), vx, modqyZDCA);
1400+
histos.fill(HIST("hvxQxZDCC"), vx, modqxZDCC);
1401+
histos.fill(HIST("hvxQyZDCC"), vx, modqyZDCC);
1402+
1403+
histos.fill(HIST("hvyQxZDCA"), vy, modqxZDCA);
1404+
histos.fill(HIST("hvyQyZDCA"), vy, modqyZDCA);
1405+
histos.fill(HIST("hvyQxZDCC"), vy, modqxZDCC);
1406+
histos.fill(HIST("hvyQyZDCC"), vy, modqyZDCC);
1407+
1408+
histos.fill(HIST("hvzQxZDCA"), vz, modqxZDCA);
1409+
histos.fill(HIST("hvzQyZDCA"), vz, modqyZDCA);
1410+
histos.fill(HIST("hvzQxZDCC"), vz, modqxZDCC);
1411+
histos.fill(HIST("hvzQyZDCC"), vz, modqyZDCC);
1412+
1413+
histos.fill(HIST("PsiZDCA"), centrality, psiZDCA);
1414+
histos.fill(HIST("PsiZDCC"), centrality, psiZDCC);
1415+
histos.fill(HIST("PsiZDC"), centrality, psiZDC);
1416+
}
1417+
13221418
// fill histograms
13231419
histos.fill(HIST("hCentrality"), centrality);
13241420
if (!checkwithpub) {

0 commit comments

Comments
 (0)