diff --git a/PWGCF/Flow/TableProducer/zdcQVectors.cxx b/PWGCF/Flow/TableProducer/zdcQVectors.cxx index 07befdd4d31..2cc9b01a38c 100644 --- a/PWGCF/Flow/TableProducer/zdcQVectors.cxx +++ b/PWGCF/Flow/TableProducer/zdcQVectors.cxx @@ -20,8 +20,8 @@ #include "Common/CCDB/RCTSelectionFlags.h" #include "Common/DataModel/Centrality.h" #include "Common/DataModel/EventSelection.h" -#include "Common/DataModel/TrackSelectionTables.h" #include "Common/DataModel/Multiplicity.h" +#include "Common/DataModel/TrackSelectionTables.h" #include #include @@ -170,10 +170,10 @@ struct ZdcQVectors { // define my..... // Filter collisionFilter = nabs(aod::collision::posZ) <; - + using UsedCollisions = soa::Join; using BCsRun3 = soa::Join; - Filter trackFilter = nabs(aod::track::eta) < cfgTrackSelsEta && aod::track::pt > cfgTrackSelsPtmin && aod::track::pt < cfgTrackSelsPtmax && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t)true)) && nabs(aod::track::dcaXY) < cfgTrackSelsDCAxy && nabs(aod::track::dcaZ) < cfgTrackSelsDCAz; + Filter trackFilter = nabs(aod::track::eta) < cfgTrackSelsEta && aod::track::pt > cfgTrackSelsPtmin&& aod::track::pt < cfgTrackSelsPtmax && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t)true)) && nabs(aod::track::dcaXY) < cfgTrackSelsDCAxy&& nabs(aod::track::dcaZ) < cfgTrackSelsDCAz; using UnfilteredTracks = soa::Join; using UsedTracks = soa::Filtered; @@ -411,7 +411,7 @@ struct ZdcQVectors { fMultPVCutHigh = std::make_unique("fMultPVCutHigh", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x", 0, 100); fMultCutLow = std::make_unique("fMultCutLow", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x", 0, 100); fMultCutHigh = std::make_unique("fMultCutHigh", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x", 0, 100); - + std::vector paramsMultPVCut = cfgEvSelsMultPv; std::vector paramsMultCut = cfgEvSelsMult; @@ -426,8 +426,8 @@ struct ZdcQVectors { fMultPVCutHigh->SetParameters(paramsMultPVCut[5], paramsMultPVCut[6], paramsMultPVCut[7], paramsMultPVCut[8], paramsMultPVCut[9]); fMultCutLow->SetParameters(paramsMultCut[0], paramsMultCut[1], paramsMultCut[2], paramsMultCut[3], paramsMultCut[4]); fMultCutHigh->SetParameters(paramsMultCut[5], paramsMultCut[6], paramsMultCut[7], paramsMultCut[8], paramsMultCut[9]); + } } -} double rescaleTimestamp(uint64_t timestamp, int runnumber) { @@ -504,7 +504,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_Zvtx); fillCutAnalysis(collision, bunchCrossing, evSel_Zvtx); - } else if(cfgSelVec.value[evSel_Zvtx]){ + } else if (cfgSelVec.value[evSel_Zvtx]) { isEventSelected = false; } @@ -512,7 +512,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_sel8); fillCutAnalysis(collision, bunchCrossing, evSel_sel8); - } else if(cfgSelVec.value[evSel_sel8]){ + } else if (cfgSelVec.value[evSel_sel8]) { isEventSelected = false; } @@ -521,7 +521,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_occupancy); fillCutAnalysis(collision, bunchCrossing, evSel_occupancy); - } else if(cfgSelVec.value[evSel_occupancy]){ + } else if (cfgSelVec.value[evSel_occupancy]) { isEventSelected = false; } @@ -529,7 +529,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_kNoSameBunchPileup); fillCutAnalysis(collision, bunchCrossing, evSel_kNoSameBunchPileup); - } else if(cfgSelVec.value[evSel_kNoSameBunchPileup]){ + } else if (cfgSelVec.value[evSel_kNoSameBunchPileup]) { isEventSelected = false; } @@ -537,7 +537,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_kIsGoodZvtxFT0vsPV); fillCutAnalysis(collision, bunchCrossing, evSel_kIsGoodZvtxFT0vsPV); - } else if(cfgSelVec.value[evSel_kIsGoodZvtxFT0vsPV]){ + } else if (cfgSelVec.value[evSel_kIsGoodZvtxFT0vsPV]) { isEventSelected = false; } @@ -545,7 +545,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_kNoCollInTimeRangeStandard); fillCutAnalysis(collision, bunchCrossing, evSel_kNoCollInTimeRangeStandard); - } else if(cfgSelVec.value[evSel_kNoCollInTimeRangeStandard]){ + } else if (cfgSelVec.value[evSel_kNoCollInTimeRangeStandard]) { isEventSelected = false; } @@ -553,7 +553,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_kNoCollInTimeRangeNarrow); fillCutAnalysis(collision, bunchCrossing, evSel_kNoCollInTimeRangeNarrow); - } else if(cfgSelVec.value[evSel_kNoCollInTimeRangeNarrow]){ + } else if (cfgSelVec.value[evSel_kNoCollInTimeRangeNarrow]) { isEventSelected = false; } @@ -561,7 +561,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_kIsVertexITSTPC); fillCutAnalysis(collision, bunchCrossing, evSel_kIsVertexITSTPC); - } else if(cfgSelVec.value[evSel_kIsVertexITSTPC]){ + } else if (cfgSelVec.value[evSel_kIsVertexITSTPC]) { isEventSelected = false; } @@ -569,7 +569,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_kIsGoodITSLayersAll); fillCutAnalysis(collision, bunchCrossing, evSel_kIsGoodITSLayersAll); - } else if(cfgSelVec.value[evSel_kIsGoodITSLayersAll]){ + } else if (cfgSelVec.value[evSel_kIsGoodITSLayersAll]) { isEventSelected = false; } @@ -577,7 +577,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_kIsGoodITSLayer0123); fillCutAnalysis(collision, bunchCrossing, evSel_kIsGoodITSLayer0123); - } else if(cfgSelVec.value[evSel_kIsGoodITSLayer0123]){ + } else if (cfgSelVec.value[evSel_kIsGoodITSLayer0123]) { isEventSelected = false; } @@ -585,7 +585,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_RCTFlagsZDC); fillCutAnalysis(collision, bunchCrossing, evSel_RCTFlagsZDC); - } else if(cfgSelVec.value[evSel_RCTFlagsZDC]){ + } else if (cfgSelVec.value[evSel_RCTFlagsZDC]) { isEventSelected = false; } @@ -600,8 +600,8 @@ struct ZdcQVectors { } auto multNTracksPV = collision.multNTracksPV(); - selected = true; - + selected = true; + if (vtxz > cfgVtxZ || vtxz < -cfgVtxZ) selected = false; if (multNTracksPV < fMultPVCutLow->Eval(collision.centFT0C())) @@ -616,7 +616,7 @@ struct ZdcQVectors { if (selected) { selectionBits |= static_cast(0x1u << evSel_MultCut); fillCutAnalysis(collision, bunchCrossing, evSel_MultCut); - } else if(cfgSelVec.value[evSel_MultCut]){ + } else if (cfgSelVec.value[evSel_MultCut]) { isEventSelected = false; } @@ -909,7 +909,7 @@ struct ZdcQVectors { // Use this bool to check for given set of event selections of Q-vectors would be selected // Enable plotting only if event would be selected - bool isEventSelected = true; + bool isEventSelected = true; uint16_t eventSelectionFlags = eventSelected(collision, foundBC.zdc(), isEventSelected, tracks.size()); @@ -996,7 +996,7 @@ struct ZdcQVectors { calibtower++; } - if (cfgFillHistRegistry && !cfgFillNothing && isEventSelected) { + if (cfgFillHistRegistry && !cfgFillNothing && isEventSelected) { for (int i = 0; i < nTowersPerSide; i++) { float bincenter = i + .5; registry.fill(HIST("QA/ZNA_Energy"), bincenter, eZN[i]); @@ -1055,7 +1055,7 @@ struct ZdcQVectors { } } - if (cfgFillHistRegistry && !cfgFillNothing && isEventSelected) { + if (cfgFillHistRegistry && !cfgFillNothing && isEventSelected) { registry.get(HIST("QA/before/ZNA_Qx"))->Fill(Form("%d", runnumber), q[0]); registry.get(HIST("QA/before/ZNA_Qy"))->Fill(Form("%d", runnumber), q[1]); registry.get(HIST("QA/before/ZNC_Qx"))->Fill(Form("%d", runnumber), q[2]); @@ -1087,7 +1087,7 @@ struct ZdcQVectors { lastRunNumber = runnumber; return; } else { - if (cfgFillHistRegistry && isEventSelected) + if (cfgFillHistRegistry && isEventSelected) fillCommonRegistry(q[0], q[1], q[2], q[3], v, centrality, rsTimestamp); // vector of 4