// Set default cuts
fEtaRange.SetLimits(-0.8, 0.8);
+ fPtRange.SetLimits(0.15, 100.);
}
// Loop over all tracks (No cuts applied)
for(int itrk = 0; itrk < fInputEvent->GetNumberOfTracks(); ++itrk){
track = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(itrk));
- // first fill without pielup cut
- if(fEtaRange.IsInRange(track->Eta())) continue;
+ if(!fEtaRange.IsInRange(track->Eta())) continue;
+ if(!fPtRange.IsInRange(track->Pt())) continue;
if(triggers[0]) FillTrackHist("MinBias", track, zv, isPileupEvent, 0);
if(!triggerstrings.size()) // Non-EMCal-triggered
FillTrackHist("NoEMCal", track, zv, isPileupEvent, 0);
TIter trackIter(acceptedTracks.get());
while((track = dynamic_cast<AliESDtrack *>(trackIter()))){
if(!fEtaRange.IsInRange(track->Eta())) continue;
+ if(!fPtRange.IsInRange(track->Pt())) continue;
if(triggers[0]) FillTrackHist("MinBias", track, zv, isPileupEvent, icut + 1);
if(!triggerstrings.size()) // Non-EMCal-triggered
FillTrackHist("NoEMCal", track, zv, isPileupEvent, icut + 1);
void AddTrackCuts(AliESDtrackCuts *trackCuts) { fListTrackCuts->Add(trackCuts); }
void SetEtaRange(double etamin, double etamax) { fEtaRange.SetLimits(etamin, etamax); }
+ void SetPtRange(double ptmin, double ptmax) { fPtRange.SetLimits(ptmin, ptmax); }
private:
AliAnalysisTaskPtEMCalTrigger(const AliAnalysisTaskPtEMCalTrigger &);
// Cuts
AliCutValueRange<double> fEtaRange; // Eta Selection Range
+ AliCutValueRange<double> fPtRange; // Pt Selection Range
ClassDef(AliAnalysisTaskPtEMCalTrigger, 1); // Analysis of EMCal triggered events
};
EMCalTriggerPtAnalysis::AliAnalysisTaskPtEMCalTrigger *pttriggertask = new EMCalTriggerPtAnalysis::AliAnalysisTaskPtEMCalTrigger("ptemcaltriggertask");
pttriggertask->SelectCollisionCandidates(AliVEvent::kINT7 | AliVEvent::kEMC7); // Select both INT7 or EMC7 triggered events
mgr->AddTask(pttriggertask);
+ pttriggertask->SetPtRange(2., 100.);
// Create charged hadrons pPb standard track cuts
AliESDtrackCuts *standardTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(true, 1);