Skip to content
Merged
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
47 changes: 22 additions & 25 deletions PWGCF/Flow/Tasks/flowGfwTask.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -80,7 +80,7 @@ struct FlowGfwTask {
O2_DEFINE_CONFIGURABLE(cfgNoSameBunchPileup, bool, false, "kNoSameBunchPileup");
O2_DEFINE_CONFIGURABLE(cfgIsGoodZvtxFT0vsPV, bool, false, "kIsGoodZvtxFT0vsPV");
O2_DEFINE_CONFIGURABLE(cfgNoCollInTimeRangeStandard, bool, false, "kNoCollInTimeRangeStandard");
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodITSLayersAll, bool, false, "kIsGoodITSLayersAll")
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodITSLayersAll, bool, false, "kIsGoodITSLayersAll")
O2_DEFINE_CONFIGURABLE(cfgOccupancy, bool, false, "Bool for event selection on detector occupancy");
O2_DEFINE_CONFIGURABLE(cfgMultCut, bool, false, "Use additional event cut on mult correlations");
O2_DEFINE_CONFIGURABLE(FineBinning, bool, false, "Manually change to fine binning")
Expand Down Expand Up @@ -559,8 +559,8 @@ struct FlowGfwTask {
}
registry.fill(HIST("hEventCount"), kNOCOLLINTIMERANGESTANDART);
}
if (cfgEvSelkIsGoodITSLayersAll){
if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
if (cfgEvSelkIsGoodITSLayersAll) {
if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
// removes dead staves of ITS
return false;
}
Expand Down Expand Up @@ -644,7 +644,6 @@ struct FlowGfwTask {
}
}


// Apply process filters
Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex && (aod::cent::centFT0C > cfgMinCentFT0C) && (aod::cent::centFT0C < cfgMaxCentFT0C);
Filter trackFilter = (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtMin) && (aod::track::pt < cfgCutPtMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t) true)) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls) && (nabs(aod::track::dcaZ) < cfgCutDCAz) && (nabs(aod::track::dcaXY) < cfgCutDCAxy);
Expand Down Expand Up @@ -690,25 +689,24 @@ struct FlowGfwTask {
}

const auto& foundBC = collision.foundBC_as<BCsRun3>();
if (foundBC.has_zdc())
{
registry.fill(HIST("hEventCounterForZDC"), 1);
if (foundBC.has_zdc()) {
registry.fill(HIST("hEventCounterForZDC"), 1);

// FT0 amplitude to use in fine binning
double ft0aAmp = 0;
double ft0cAmp = 0;
double ft0mAmp = 0;

if (foundBC.has_ft0()) {
if (foundBC.has_ft0()) {
for (const auto& amplitude : foundBC.ft0().amplitudeA()) {
ft0aAmp += amplitude;
ft0aAmp += amplitude;
}
for (const auto& amplitude : foundBC.ft0().amplitudeC()) {
ft0cAmp += amplitude;
}
} else{
}
} else {
ft0aAmp = ft0cAmp = -999;
}
}

registry.fill(HIST("FT0AAmp"), ft0aAmp);
registry.fill(HIST("FT0CAmp"), ft0cAmp);
Expand All @@ -725,10 +723,10 @@ struct FlowGfwTask {
auto aZEM1 = zdcread.amplitudeZEM1();
auto aZEM2 = zdcread.amplitudeZEM2();

registry.fill(HIST("ZNAcoll"),aZNA);
registry.fill(HIST("ZNCcoll"),aZNC);
registry.fill(HIST("ZPAcoll"),aZPA);
registry.fill(HIST("ZPCcoll"),aZPC);
registry.fill(HIST("ZNAcoll"), aZNA);
registry.fill(HIST("ZNCcoll"), aZNC);
registry.fill(HIST("ZPAcoll"), aZPA);
registry.fill(HIST("ZPCcoll"), aZPC);

registry.fill(HIST("ZNvsFT0correl"), (ft0aAmp + ft0cAmp) / 100., aZNC + aZNA);

Expand All @@ -739,18 +737,18 @@ struct FlowGfwTask {

registry.fill(HIST("ZNvsZEMcoll"), aZEM1 + aZEM2, aZNA + aZNC);

if(centrality>=0 && centrality<=5){
if (centrality >= 0 && centrality <= 5) {
registry.fill(HIST("ZNvsZEMcoll05"), aZEM1 + aZEM2, aZNA + aZNC);
}else if (centrality>5 && centrality<=10){
} else if (centrality > 5 && centrality <= 10) {
registry.fill(HIST("ZNvsZEMcoll510"), aZEM1 + aZEM2, aZNA + aZNC);
}else if (centrality>10 && centrality<=20){
} else if (centrality > 10 && centrality <= 20) {
registry.fill(HIST("ZNvsZEMcoll1020"), aZEM1 + aZEM2, aZNA + aZNC);
}else if (centrality>20 && centrality<=30){
} else if (centrality > 20 && centrality <= 30) {
registry.fill(HIST("ZNvsZEMcoll2030"), aZEM1 + aZEM2, aZNA + aZNC);
}else{
} else {
registry.fill(HIST("ZNvsZEMcollrest"), aZEM1 + aZEM2, aZNA + aZNC);
}
}//End of ZDC
} // End of ZDC

float vtxz = collision.posZ();
float lRandom = fRndm->Rndm();
Expand Down Expand Up @@ -788,8 +786,7 @@ struct FlowGfwTask {
// track loop
int globalTracksNch = 0;

for (const auto& track : tracks)
{
for (const auto& track : tracks) {
if (!trackSelected(track))
continue;

Expand Down Expand Up @@ -891,7 +888,7 @@ struct FlowGfwTask {
fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32Nch05etagap], globalTracksNch);
fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34Nch05], globalTracksNch);

registry.fill(HIST("Nch"),globalTracksNch);
registry.fill(HIST("Nch"), globalTracksNch);
}

// Filling Flow Container
Expand Down
Loading