fMCMultRecTrackHist1(0),
// rec. track control histograms
- fRecTrackHist3(0)
+ fRecTrackHist3(0),
+
+ fTriggerAnalysis(0)
{
// default constructor
for(Int_t i=0; i<AlidNdPtHelper::kCutSteps; i++) {
fMCMultRecTrackHist1(0),
// rec. track control histograms
- fRecTrackHist3(0)
+ fRecTrackHist3(0),
+
+ fTriggerAnalysis(0)
{
//
// constructor
if(fRecTrackHist3) delete fRecTrackHist3; fRecTrackHist3=0;
//
if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;
+
+ if (fTriggerAnalysis) delete fTriggerAnalysis; fTriggerAnalysis = 0;
}
//_____________________________________________________________________________
// init folder
fAnalysisFolder = CreateFolder("folderdNdPt","Analysis dNdPt Folder");
+
+ // init trigger analysis (for zdc cut)
+ fTriggerAnalysis = new AliTriggerAnalysis;
}
//_____________________________________________________________________________
AliDebug(AliLog::kError, "esdEvent not available");
return;
}
+
+ // zdc cut
+ if (!fTriggerAnalysis->ZDCTimeTrigger(esdEvent)) {
+ return;
+ }
+ // track cuts from Jochen
+ const AliESDVertex* vtxESDTPC = esdEvent->GetPrimaryVertexTPC();
+ if( vtxESDTPC->GetNContributors() < 1 ) {
+ return;
+ }
+
+ // francesco prino cut
+ const AliMultiplicity* multESD = esdEvent->GetMultiplicity();
+ if( vtxESDTPC->GetNContributors() < (-10.+0.25*multESD->GetNumberOfITSClusters(0)) ) {
+ return;
+ }
// get selection cuts
AlidNdPtEventCuts *evtCuts = GetEventCuts();
class AliESDfriend;
class AliESDfriendTrack;
class AlidNdPtHelper;
+class AliTriggerAnalysis;
#include "AlidNdPt.h"
// track control histograms
THnSparseF *fRecTrackHist3; //-> nclust:chi2:Pt:Eta:Phi
+ AliTriggerAnalysis *fTriggerAnalysis; //! trigger analysis object;
+
AlidNdPtAnalysisPbPb(const AlidNdPtAnalysisPbPb&); // not implemented
AlidNdPtAnalysisPbPb& operator=(const AlidNdPtAnalysisPbPb&); // not implemented