Skip to content
Closed
Changes from all commits
Commits
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
26 changes: 18 additions & 8 deletions PWGUD/Tasks/flowCumulantsUpc.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -70,9 +70,11 @@ struct FlowCumulantsUpc {

O2_DEFINE_CONFIGURABLE(cfgCutVertex, float, 10.0f, "Accepted z-vertex range")
O2_DEFINE_CONFIGURABLE(cfgIfVertex, bool, false, "choose vertex or not")
O2_DEFINE_CONFIGURABLE(cfgCutPtPOIMin, float, 0.1f, "Minimal pT for poi tracks")
O2_DEFINE_CONFIGURABLE(cfgPtCutMin, float, 0.2f, "minimum accepted track pT")
O2_DEFINE_CONFIGURABLE(cfgPtCutMax, float, 10.0f, "maximum accepted track pT")
O2_DEFINE_CONFIGURABLE(cfgCutPtPOIMin, float, 0.2f, "Minimal pT for poi tracks")
O2_DEFINE_CONFIGURABLE(cfgCutPtPOIMax, float, 10.0f, "Maximal pT for poi tracks")
O2_DEFINE_CONFIGURABLE(cfgCutPtRefMin, float, 0.1f, "Minimal pT for ref tracks")
O2_DEFINE_CONFIGURABLE(cfgCutPtRefMin, float, 0.2f, "Minimal pT for ref tracks")
O2_DEFINE_CONFIGURABLE(cfgCutPtRefMax, float, 3.0f, "Maximal pT for ref tracks")
O2_DEFINE_CONFIGURABLE(cfgCutEta, float, 0.8f, "max Eta range for tracks")
O2_DEFINE_CONFIGURABLE(cfgCutTPCCrossedRows, float, 70.0f, "minimum number of crossed TPC Rows")
Expand Down Expand Up @@ -101,6 +103,11 @@ struct FlowCumulantsUpc {
O2_DEFINE_CONFIGURABLE(cfgZdcTime, bool, false, "choose zdc time cut")
O2_DEFINE_CONFIGURABLE(cfgZdcTimeCut, float, 2.0, "zdc time cut")
O2_DEFINE_CONFIGURABLE(cfgIRMaxCut, double, 50, "maximum interaction rate for UPC events")
O2_DEFINE_CONFIGURABLE(cfgMaxDebugEvents, int, -1, "Maximum number of events for debug (-1: all events)")
O2_DEFINE_CONFIGURABLE(cfgvtxITSTPC, bool, true, "require vertex from ITS+TPC")
O2_DEFINE_CONFIGURABLE(cfgsbp, bool, true, "require sbp")
O2_DEFINE_CONFIGURABLE(cfgitsROFb, bool, true, "require itsROFb")
O2_DEFINE_CONFIGURABLE(cfgtfb, bool, true, "require tfb")

Configurable<std::vector<std::string>> cfgUserDefineGFWCorr{"cfgUserDefineGFWCorr", std::vector<std::string>{"refN02 {2} refP02 {-2}", "refN12 {2} refP12 {-2}"}, "User defined GFW CorrelatorConfig"};
Configurable<std::vector<std::string>> cfgUserDefineGFWName{"cfgUserDefineGFWName", std::vector<std::string>{"Ch02Gap22", "Ch12Gap22"}, "User defined GFW Name"};
Expand Down Expand Up @@ -173,8 +180,8 @@ struct FlowCumulantsUpc {
// Event QA
registry.add("hEventCount", "Number of Event;; Count", {HistType::kTH1D, {{13, 0, 13}}});
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(1, "Filtered event");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(2, "after ZDC timing");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(3, "after gapside");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(2, "after gapside");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(3, "after ZDC timing");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(4, "after Vz");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(5, "after occupancy");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(6, "after load corr");
Expand Down Expand Up @@ -597,6 +604,9 @@ struct FlowCumulantsUpc {
bool trackSelected(TTrack track)
{
registry.fill(HIST("hTrackCount"), 0.5);
if (track.pt() < cfgPtCutMin || track.pt() > cfgPtCutMax) {
return false;
}
// UPC selection
if (!track.isPVContributor()) {
return false;
Expand Down Expand Up @@ -748,19 +758,19 @@ struct FlowCumulantsUpc {
}
}
registry.fill(HIST("hEventCount"), 5.5);
if (!collision.vtxITSTPC()) {
if (cfgvtxITSTPC && !collision.vtxITSTPC()) {
return;
}
registry.fill(HIST("hEventCount"), 6.5);
if (!collision.sbp()) {
if (cfgsbp && !collision.sbp()) {
return;
}
registry.fill(HIST("hEventCount"), 7.5);
if (!collision.itsROFb()) {
if (cfgitsROFb && !collision.itsROFb()) {
return;
}
registry.fill(HIST("hEventCount"), 8.5);
if (!collision.tfb()) {
if (cfgtfb && !collision.tfb()) {
return;
}
registry.fill(HIST("hEventCount"), 9.5);
Expand Down
Loading