Skip to content

[PWGLF] Chk892Flow - Add RCT flag Checker #10867

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
BongHwi merged 1 commit into from
Apr 15, 2025
Merged

[PWGLF] Chk892Flow - Add RCT flag Checker #10867

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
16 changes: 15 additions & 1 deletion PWGLF/Tasks/Resonances/chk892Flow.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -78,6 +78,7 @@
using namespace o2::framework::expressions;
using namespace o2::soa;
using namespace o2::constants::physics;
using namespace o2::aod::rctsel;

struct Chk892Flow {
enum BinType : unsigned int {
Expand Down Expand Up @@ -148,7 +149,12 @@
Configurable<bool> cfgEvtCollInTimeRangeStandard{"cfgEvtCollInTimeRangeStandard", true, "Evt sel: apply NoCollInTimeRangeStandard"};
Configurable<float> cfgEventCentralityMin{"cfgEventCentralityMin", 0.0f, "Event sel: minimum centrality"};
Configurable<float> cfgEventCentralityMax{"cfgEventCentralityMax", 80.0f, "Event sel: maximum centrality"};
Configurable<bool> cfgEvtUseRCTFlagChecker{"cfgEvtUseRCTFlagChecker", false, "Evt sel: use RCT flag checker"};
Configurable<std::string> cfgEvtRCTFlagCheckerLabel{"cfgEvtRCTFlagCheckerLabel", "CBT_hadronPID", "Evt sel: RCT flag checker label"};
Configurable<bool> cfgEvtRCTFlagCheckerZDCCheck{"cfgEvtRCTFlagCheckerZDCCheck", false, "Evt sel: RCT flag checker ZDC check"};
Configurable<bool> cfgEvtRCTFlagCheckerLimitAcceptAsBad{"cfgEvtRCTFlagCheckerLimitAcceptAsBad", false, "Evt sel: RCT flag checker treat Limited Acceptance As Bad"};
} EventCuts;
RCTFlagsChecker rctChecker;

/// PID Selections, pion
struct : ConfigurableGroup {
Expand Down Expand Up @@ -262,6 +268,8 @@
colCuts.setApplyCollInTimeRangeStandard(EventCuts.cfgEvtCollInTimeRangeStandard);
colCuts.printCuts();

rctChecker.init(EventCuts.cfgEvtRCTFlagCheckerLabel, EventCuts.cfgEvtRCTFlagCheckerZDCCheck, EventCuts.cfgEvtRCTFlagCheckerLimitAcceptAsBad);

AxisSpec centAxis = {AxisConfig.cfgBinsCent, "T0M (%)"};
AxisSpec vtxzAxis = {AxisConfig.cfgBinsVtxZ, "Z Vertex (cm)"};
AxisSpec epAxis = {100, -1.0 * constants::math::PI, constants::math::PI};
Expand Down Expand Up @@ -458,7 +466,7 @@
lRefAId = 4;
lRefBId = 5;
}
if (EventPlaneConfig.cfgNQvec < 2) {

Check warning on line 469 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOG(fatal) << "nMode must be larger than 1, current input (cfgNQvec): " << EventPlaneConfig.cfgNQvec;
}
LOGF(info, "lDetId: %d, lRefAId: %d, lRefBId: %d", lDetId, lRefAId, lRefBId);
Expand Down Expand Up @@ -534,7 +542,7 @@
{
if (AnalysisConfig.cfgCentEst == 1) {
return collision.centFT0C();
} else if (AnalysisConfig.cfgCentEst == 2) {

Check warning on line 545 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
return collision.centFT0M();
} else {
return -999;
Expand Down Expand Up @@ -758,16 +766,16 @@
if (std::abs(motherbTrack.pdgCode()) != kKstarPlus) // Are you charged Kstar's daughter?
return false; // Apply first since it's more restrictive

if (std::abs(motherkV0.pdgCode()) != 310) // Is it K0s?

Check warning on line 769 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.

Check warning on line 769 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/explicit-code] 

Avoid using hard-coded PDG codes. Use named values from PDG_t or o2::constants::physics::Pdg instead.
return false;
// Check if K0s's mother is K0 (311)
auto motherK0 = motherkV0.template mothers_as<aod::McParticles>();
if (std::abs(motherK0.pdgCode()) != 311)

Check warning on line 773 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.

Check warning on line 773 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/explicit-code] 

Avoid using hard-coded PDG codes. Use named values from PDG_t or o2::constants::physics::Pdg instead.
return false;

// Check if K0's mother is Kstar (323)
auto motherKstar = motherK0.template mothers_as<aod::McParticles>();
if (std::abs(motherKstar.pdgCode()) != 323)

Check warning on line 778 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.

Check warning on line 778 in PWGLF/Tasks/Resonances/chk892Flow.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/explicit-code] 

Avoid using hard-coded PDG codes. Use named values from PDG_t or o2::constants::physics::Pdg instead.
return false;

// Check if bTrack and K0 have the same mother (global index)
Expand Down Expand Up @@ -1066,12 +1074,15 @@
{
if (!colCuts.isSelected(collision)) // Default event selection
return;
if (EventCuts.cfgEvtUseRCTFlagChecker && !rctChecker(collision)) {
return;
}
if (AnalysisConfig.cfgQvecSel && (collision.qvecAmp()[lDetId] < 1e-4 || collision.qvecAmp()[lRefAId] < 1e-4 || collision.qvecAmp()[lRefBId] < 1e-4))
return; // If we don't have a Q-vector
colCuts.fillQA(collision);
lCentrality = getCentrality(collision);
if (lCentrality < EventCuts.cfgEventCentralityMin || lCentrality > EventCuts.cfgEventCentralityMax)
return;
colCuts.fillQA(collision);

fillHistograms<false, false>(collision, tracks, v0s, EventPlaneConfig.cfgnMods); // second order
}
Expand All @@ -1082,6 +1093,9 @@
MCTrackCandidates const& tracks,
MCV0Candidates const& v0s)
{
if (EventCuts.cfgEvtUseRCTFlagChecker && !rctChecker(collision)) {
return;
}
fillHistograms<true, false>(collision, tracks, v0s, EventPlaneConfig.cfgnMods);
}
PROCESS_SWITCH(Chk892Flow, processMC, "Process Event for MC", false);
Expand Down
Loading