Skip to content

Commit bac70b9

Browse files
Maxim VirtaMaxim Virta
authored andcommitted
Fixed eta cut for charged pt spectra
1 parent 72a7f0f commit bac70b9

File tree

1 file changed

+3
-5
lines changed

1 file changed

+3
-5
lines changed

PWGCF/GenericFramework/Tasks/flowGfwV02.cxx

Lines changed: 3 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -180,7 +180,7 @@ struct FlowGfwV02 {
180180
// // The analysis assumes the data has been subjected to a QA of its selection,
181181
// // and thus only the final distributions of the data for analysis are saved.
182182
o2::framework::expressions::Filter collFilter = (nabs(aod::collision::posZ) < cfgEventCuts.cfgZvtxMax);
183-
o2::framework::expressions::Filter trackFilter = (aod::track::pt > cfgTrackCuts.cfgPtMin) && (aod::track::pt < cfgTrackCuts.cfgPtMax) && (nabs(aod::track::eta) < cfgTrackCuts.cfgEtaMax && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t)true)) && (aod::track::itsChi2NCl < cfgTrackCuts.cfgChi2PrITSCls) && (aod::track::tpcChi2NCl < cfgTrackCuts.cfgChi2PrTPCCls) && nabs(aod::track::dcaZ) < cfgTrackCuts.cfgDCAz);
183+
o2::framework::expressions::Filter trackFilter = (aod::track::pt > cfgTrackCuts.cfgPtMin) && (aod::track::pt < cfgTrackCuts.cfgPtMax) && nabs(aod::track::eta) < cfgTrackCuts.cfgEtaMax && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t)true)) && (aod::track::itsChi2NCl < cfgTrackCuts.cfgChi2PrITSCls) && (aod::track::tpcChi2NCl < cfgTrackCuts.cfgChi2PrTPCCls) && nabs(aod::track::dcaZ) < cfgTrackCuts.cfgDCAz;
184184

185185
// Connect to ccdb
186186
Service<ccdb::BasicCCDBManager> ccdb;
@@ -878,7 +878,6 @@ struct FlowGfwV02 {
878878
}
879879
}
880880
// Fill the profiles for each pT bin
881-
// printf("Config name: %s\n", corrconfigs.at(0).Head.c_str());
882881
auto dnx = fGFW->Calculate(corrconfigs.at(0), 0, kTRUE).real();
883882
if (dnx == 0)
884883
return;
@@ -889,7 +888,6 @@ struct FlowGfwV02 {
889888
ptFraction = pidStates.hPtMid[PidCharged]->GetBinContent(i) / pidStates.hPtMid[PidCharged]->Integral();
890889
if (std::abs(val) < 1)
891890
registry.fill(HIST("v02pt"), fSecondAxis->GetBinCenter(i), centmult, val * ptFraction, (cfgUseMultiplicityFlowWeights) ? dnx : 1.0);
892-
// printf("bincenter hPtMid: %f, fsecondaxis: %f\n", hPtMid->GetBinCenter(i), fSecondAxis->GetBinCenter(i));
893891
registry.fill(HIST("nchMid"), fSecondAxis->GetBinCenter(i), centmult, ptFraction);
894892
}
895893
}
@@ -931,7 +929,8 @@ struct FlowGfwV02 {
931929
AcceptedTracks acceptedTracks{0, 0, 0, 0};
932930
for (const auto& track : tracks) {
933931
processTrack(track, vtxz, xaxis.multiplicity, run, acceptedTracks);
934-
pidStates.hPtMid[PidCharged]->Fill(track.pt(), getEfficiency(track));
932+
if (track.eta() > -0.4 && track.eta() < 0.4)
933+
pidStates.hPtMid[PidCharged]->Fill(track.pt(), getEfficiency(track));
935934
// If PID is identified, fill pt spectrum for the corresponding particle
936935
int pidInd = getNsigmaPID(track);
937936
if (pidInd != -1 && track.eta() > -0.4 && track.eta() < 0.4) {
@@ -1096,7 +1095,6 @@ struct FlowGfwV02 {
10961095
lastRun = run;
10971096
LOGF(info, "run = %d", run);
10981097
}
1099-
11001098
loadCorrections(bc);
11011099

11021100
registry.fill(HIST("eventQA/eventSel"), kFilteredEvent);

0 commit comments

Comments
 (0)