diff --git a/PWGCF/Flow/Tasks/flowGfwTask.cxx b/PWGCF/Flow/Tasks/flowGfwTask.cxx index f9e35f7ccf8..a7cd60e2653 100644 --- a/PWGCF/Flow/Tasks/flowGfwTask.cxx +++ b/PWGCF/Flow/Tasks/flowGfwTask.cxx @@ -615,50 +615,50 @@ struct FlowGfwTask { return false; } registry.fill(HIST("hEventCount"), kISVERTEXITSTPC); - if (cfgNoCollInTimeRangeStandard) { - if (!collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) { - // no collisions in specified time range - return false; + if (cfgNoCollInTimeRangeStandard) { + if (!collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) { + // no collisions in specified time range + return false; + } + registry.fill(HIST("hEventCount"), kNOCOLLINTIMERANGESTANDART); } - registry.fill(HIST("hEventCount"), kNOCOLLINTIMERANGESTANDART); - } - if (cfgEvSelkIsGoodITSLayersAll) { - if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) { - // removes dead staves of ITS - return false; + if (cfgEvSelkIsGoodITSLayersAll) { + if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) { + // removes dead staves of ITS + return false; + } + registry.fill(HIST("hEventCount"), kISGOODITSLAYERSALL); } - registry.fill(HIST("hEventCount"), kISGOODITSLAYERSALL); - } - - float vtxz = -999; - if (collision.numContrib() > 1) { - vtxz = collision.posZ(); - float zRes = std::sqrt(collision.covZZ()); - if (zRes > 0.25 && collision.numContrib() < 20) - vtxz = -999; - } - auto multNTracksPV = collision.multNTracksPV(); + float vtxz = -999; + if (collision.numContrib() > 1) { + vtxz = collision.posZ(); + float zRes = std::sqrt(collision.covZZ()); + if (zRes > 0.25 && collision.numContrib() < 20) + vtxz = -999; + } - if (std::abs(vtxz) > cfgCutVertex) - return false; + auto multNTracksPV = collision.multNTracksPV(); - if (cfgMultCut) { - if (multNTracksPV < fMultPVCutLow->Eval(centrality)) - return false; - if (multNTracksPV > fMultPVCutHigh->Eval(centrality)) - return false; - if (multTrk < fMultCutLow->Eval(centrality)) + if (std::abs(vtxz) > cfgCutVertex) return false; - if (multTrk > fMultCutHigh->Eval(centrality)) - return false; - registry.fill(HIST("hEventCount"), kAFTERMULTCUTS); - } - // V0A T0A 5 sigma cut - if(cfgV0AT0A5Sigma) - if (std::abs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > 5 * fT0AV0ASigma->Eval(collision.multFT0A())) - return false; + if (cfgMultCut) { + if (multNTracksPV < fMultPVCutLow->Eval(centrality)) + return false; + if (multNTracksPV > fMultPVCutHigh->Eval(centrality)) + return false; + if (multTrk < fMultCutLow->Eval(centrality)) + return false; + if (multTrk > fMultCutHigh->Eval(centrality)) + return false; + registry.fill(HIST("hEventCount"), kAFTERMULTCUTS); + } + + // V0A T0A 5 sigma cut + if (cfgV0AT0A5Sigma) + if (std::abs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > 5 * fT0AV0ASigma->Eval(collision.multFT0A())) + return false; } return true; @@ -706,13 +706,13 @@ struct FlowGfwTask { if (cfgTrackSel) { return myTrackSel.IsSelected(track); - }else if (cfgGlobalplusITS) { + } else if (cfgGlobalplusITS) { return ((track.tpcNClsFound() >= cfgCutTPCclu) && (track.itsNCls() >= cfgCutITSclu)); - }else if (cfgGlobalonly) { + } else if (cfgGlobalonly) { return ((track.tpcNClsFound() >= cfgCutTPCclu)); - }else if (cfgITSonly) { + } else if (cfgITSonly) { return ((track.itsNCls() >= cfgCutITSclu)); - }else{ + } else { return false; } } @@ -938,7 +938,7 @@ struct FlowGfwTask { } } } else { - if (cfgITSonly) { + if (cfgITSonly) { if (withinPtRef) { globalTracksNch++; registry.fill(HIST("ITSonly"), centrality, globalTracksNch);