Skip to content

Commit a980e15

Browse files
fuchuncuifuchuncui
authored
[PWGCF] Fix two bugs at flowGfwOmegaXi.cxx (#13337)
1 parent 60eb300 commit a980e15

File tree

1 file changed

+15
-1
lines changed

1 file changed

+15
-1
lines changed

‎PWGCF/Flow/Tasks/flowGfwOmegaXi.cxx‎

Lines changed: 15 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -423,6 +423,8 @@ struct FlowGfwOmegaXi {
423423
registry.add("QAhisto/V0/hqadcanegtoPVbefore", "", {HistType::kTH1D, {{1000, -10, 10}}});
424424
registry.add("QAhisto/V0/hqadcanegtoPVafter", "", {HistType::kTH1D, {{1000, -10, 10}}});
425425
// Cascade QA
426+
registry.add("QAhisto/Xi/hqaCascRadiusbefore", "", {HistType::kTH1D, {{200, -10, 10}}});
427+
registry.add("QAhisto/Xi/hqaCascRadiusafter", "", {HistType::kTH1D, {{200, -10, 10}}});
426428
registry.add("QAhisto/Xi/hqaCasccosPAbefore", "", {HistType::kTH1D, {{1000, 0.95, 1}}});
427429
registry.add("QAhisto/Xi/hqaCasccosPAafter", "", {HistType::kTH1D, {{1000, 0.95, 1}}});
428430
registry.add("QAhisto/Xi/hqaCascV0cosPAbefore", "", {HistType::kTH1D, {{1000, 0.95, 1}}});
@@ -436,6 +438,8 @@ struct FlowGfwOmegaXi {
436438
registry.add("QAhisto/Xi/hqadcaCascV0daubefore", "", {HistType::kTH1D, {{100, 0, 1}}});
437439
registry.add("QAhisto/Xi/hqadcaCascV0dauafter", "", {HistType::kTH1D, {{100, 0, 1}}});
438440

441+
registry.add("QAhisto/Xi/hqaCascRadiusbefore", "", {HistType::kTH1D, {{200, -10, 10}}});
442+
registry.add("QAhisto/Xi/hqaCascRadiusafter", "", {HistType::kTH1D, {{200, -10, 10}}});
439443
registry.add("QAhisto/Omega/hqaCasccosPAbefore", "", {HistType::kTH1D, {{1000, 0.95, 1}}});
440444
registry.add("QAhisto/Omega/hqaCasccosPAafter", "", {HistType::kTH1D, {{1000, 0.95, 1}}});
441445
registry.add("QAhisto/Omega/hqaCascV0cosPAbefore", "", {HistType::kTH1D, {{1000, 0.95, 1}}});
@@ -996,7 +1000,7 @@ struct FlowGfwOmegaXi {
9961000
if (!collision.sel8())
9971001
return;
9981002
registry.fill(HIST("hEventCount"), 1.5);
999-
if (eventSelected(collision, cent, interactionRate))
1003+
if (!eventSelected(collision, cent, interactionRate))
10001004
return;
10011005
TH1D* hLocalDensity = new TH1D("hphi", "hphi", 400, -constants::math::TwoPI, constants::math::TwoPI);
10021006
loadCorrections(bc.timestamp());
@@ -1254,6 +1258,7 @@ struct FlowGfwOmegaXi {
12541258
// fill QA
12551259
if (cfgOutputQA) {
12561260
if (isXi) {
1261+
registry.fill(HIST("QAhisto/Xi/hqaCascRadiusbefore"), casc.cascradius());
12571262
registry.fill(HIST("QAhisto/Xi/hqaCasccosPAbefore"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
12581263
registry.fill(HIST("QAhisto/Xi/hqaCascV0cosPAbefore"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
12591264
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0toPVbefore"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1262,6 +1267,7 @@ struct FlowGfwOmegaXi {
12621267
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0daubefore"), casc.dcaV0daughters());
12631268
}
12641269
if (isOmega) {
1270+
registry.fill(HIST("QAhisto/Omega/hqaCascRadiusbefore"), casc.cascradius());
12651271
registry.fill(HIST("QAhisto/Omega/hqaCasccosPAbefore"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
12661272
registry.fill(HIST("QAhisto/Omega/hqaCascV0cosPAbefore"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
12671273
registry.fill(HIST("QAhisto/Omega/hqadcaCascV0toPVbefore"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1326,6 +1332,7 @@ struct FlowGfwOmegaXi {
13261332
// fill QA
13271333
if (cfgOutputQA) {
13281334
if (isXi) {
1335+
registry.fill(HIST("QAhisto/Xi/hqaCascRadiusafter"), casc.cascradius());
13291336
registry.fill(HIST("QAhisto/Xi/hqaCasccosPAafter"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
13301337
registry.fill(HIST("QAhisto/Xi/hqaCascV0cosPAafter"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
13311338
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0toPVafter"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1334,6 +1341,7 @@ struct FlowGfwOmegaXi {
13341341
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0dauafter"), casc.dcaV0daughters());
13351342
}
13361343
if (isOmega) {
1344+
registry.fill(HIST("QAhisto/Omega/hqaCascRadiusafter"), casc.cascradius());
13371345
registry.fill(HIST("QAhisto/Omega/hqaCasccosPAafter"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
13381346
registry.fill(HIST("QAhisto/Omega/hqaCascV0cosPAafter"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
13391347
registry.fill(HIST("QAhisto/Omega/hqadcaCascV0toPVafter"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1688,6 +1696,7 @@ struct FlowGfwOmegaXi {
16881696
// fill QA
16891697
if (cfgOutputQA) {
16901698
if (std::abs(pdgCode) == kXiMinus) {
1699+
registry.fill(HIST("QAhisto/Xi/hqaCascRadiusbefore"), casc.dcabachtopv());
16911700
registry.fill(HIST("QAhisto/Xi/hqaCasccosPAbefore"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
16921701
registry.fill(HIST("QAhisto/Xi/hqaCascV0cosPAbefore"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
16931702
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0toPVbefore"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1696,6 +1705,7 @@ struct FlowGfwOmegaXi {
16961705
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0daubefore"), casc.dcaV0daughters());
16971706
}
16981707
if (std::abs(pdgCode) == kOmegaMinus) {
1708+
registry.fill(HIST("QAhisto/Omega/hqaCascRadiusbefore"), casc.dcabachtopv());
16991709
registry.fill(HIST("QAhisto/Omega/hqaCasccosPAbefore"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
17001710
registry.fill(HIST("QAhisto/Omega/hqaCascV0cosPAbefore"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
17011711
registry.fill(HIST("QAhisto/Omega/hqadcaCascV0toPVbefore"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1767,6 +1777,7 @@ struct FlowGfwOmegaXi {
17671777
}
17681778
// fill QA
17691779
if (cfgOutputQA) {
1780+
registry.fill(HIST("QAhisto/Omega/hqaCascRadiusafter"), casc.cascradius());
17701781
registry.fill(HIST("QAhisto/Omega/hqaCasccosPAafter"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
17711782
registry.fill(HIST("QAhisto/Omega/hqaCascV0cosPAafter"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
17721783
registry.fill(HIST("QAhisto/Omega/hqadcaCascV0toPVafter"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1789,6 +1800,7 @@ struct FlowGfwOmegaXi {
17891800
}
17901801
// fill QA
17911802
if (cfgOutputQA) {
1803+
registry.fill(HIST("QAhisto/Omega/hqaCascRadiusafter"), casc.cascradius());
17921804
registry.fill(HIST("QAhisto/Omega/hqaCasccosPAafter"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
17931805
registry.fill(HIST("QAhisto/Omega/hqaCascV0cosPAafter"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
17941806
registry.fill(HIST("QAhisto/Omega/hqadcaCascV0toPVafter"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1815,6 +1827,7 @@ struct FlowGfwOmegaXi {
18151827
}
18161828
// fill QA
18171829
if (cfgOutputQA) {
1830+
registry.fill(HIST("QAhisto/Xi/hqaCascRadiusafter"), casc.cascradius());
18181831
registry.fill(HIST("QAhisto/Xi/hqaCasccosPAafter"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
18191832
registry.fill(HIST("QAhisto/Xi/hqaCascV0cosPAafter"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
18201833
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0toPVafter"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));
@@ -1836,6 +1849,7 @@ struct FlowGfwOmegaXi {
18361849
registry.fill(HIST("MC/densityMCRecXi"), cascPt, nch, density, casc.mXi());
18371850
}
18381851
if (cfgOutputQA) {
1852+
registry.fill(HIST("QAhisto/Xi/hqaCascRadiusafter"), casc.cascradius());
18391853
registry.fill(HIST("QAhisto/Xi/hqaCasccosPAafter"), casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()));
18401854
registry.fill(HIST("QAhisto/Xi/hqaCascV0cosPAafter"), casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()));
18411855
registry.fill(HIST("QAhisto/Xi/hqadcaCascV0toPVafter"), casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()));

0 commit comments

Comments
 (0)