2 #ifndef ALIANALYSELEADINGTRACKUE_H
3 #define ALIANALYSELEADINGTRACKUE_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 ////////////////////////////////////////////////
8 //---------------------------------------------
9 // Class for transverse regions analysis w.r.t leading track
10 //---------------------------------------------
11 ////////////////////////////////////////////////
13 // --- ROOT system ---
17 class AliAODInputHandler;
19 class AliESDtrackCuts;
22 class AliGenPythiaEventHeader;
23 class AliInputEventHandler;
36 class AliAnalyseLeadingTrackUE : public TObject {
40 AliAnalyseLeadingTrackUE(); //constructor
41 AliAnalyseLeadingTrackUE(const AliAnalyseLeadingTrackUE & g); //copy constructor
42 AliAnalyseLeadingTrackUE & operator = (const AliAnalyseLeadingTrackUE & g); //assignment operator
43 virtual ~AliAnalyseLeadingTrackUE(); //virtual destructor
46 void SetParticleSelectionCriteria(Int_t filterbit, Bool_t onlyhadrons, Double_t etacut, Double_t etacutMin=-1., Double_t ptMin = 0) { fFilterBit = filterbit;
47 fOnlyHadrons = onlyhadrons;
48 fTrackEtaCut = etacut;
49 fTrackEtaCutMin = etacutMin;
52 void SetEventSelection(UInt_t bits) { fEventSelection = bits; }
53 void SetDebug(Int_t debug) { fDebug = debug; }
54 Bool_t ApplyCuts(TObject* track); // Reproduces the cuts of the corresponding bit in the ESD->AOD filtering
55 void DefineESDCuts(Int_t filterbit); // Emulate filterbit
56 TObjArray* FindLeadingObjects(TObject* obj); // Looks for leading track or MC particle
57 TObjArray* GetMinMaxRegion(TList* transv1, TList* transv2); // Sorts the transverse regions in MIN and MAX
58 Int_t NParticles(TObject *obj); // Counts tracks or MC particles
59 AliVParticle* ParticleWithCuts(TObject* obj, Int_t ipart, Bool_t onlyprimaries = kTRUE, Int_t particleSpecies = -1); // Returns track or MC particle at position "ipart" if passes selection criteria
60 void QSortTracks(TObjArray &a, Int_t first, Int_t last); // Sort by pT an array of AliVParticles
61 TObjArray* SortRegions(const AliVParticle* leading, TObject* obj, TObject* arrayMC, Bool_t onlyprimaries = kTRUE); // Assign particles to towards, away or transverse regions
62 TObjArray* GetAcceptedParticles(TObject* obj, TObject* arrayMC, Bool_t onlyprimaries = kTRUE, Int_t particleSpecies = -1, Bool_t useEtaPtCuts = kFALSE, Bool_t speciesOnTracks = kTRUE, Double_t evtPlane = -999.);
63 TObjArray* GetFakeParticles(TObject* obj, TObject* arrayMC, Bool_t onlyprimaries, Int_t particleSpecies, Bool_t useEtaPtCuts);
64 Bool_t TriggerSelection(const TObject* obj); // Select good triggers with AliPhysicsSelection class
65 Bool_t VertexSelection(const TObject* obj, Int_t ntracks, Double_t zed); // Vertex selection: see implementation
66 void RemoveInjectedSignals(TObjArray* tracks, TObject* arrayMC, Int_t maxLabel);
67 void RemoveWeakDecays(TObjArray* tracks, TObject* mcObj);
68 AliHelperPID* GetHelperPID() { return fHelperPID; }
69 void SetHelperPID(AliHelperPID* pid) { fHelperPID = pid; }
70 void SetTrackStatus(UInt_t status) { fTrackStatus = status; }
71 UInt_t GetTrackStatus() { return fTrackStatus; }
72 void SetCheckMotherPDG(Bool_t checkpdg) { fCheckMotherPDG = checkpdg; }
73 Bool_t GetCheckMotherPDG() { return fCheckMotherPDG; }
74 void NextEvent() { fEventCounter++; }
75 UInt_t GetEventCounter() { return fEventCounter; }
76 void SetDCAXYCut(TFormula* value) { fDCAXYCut = value; }
77 void SetSharedClusterCut(Double_t value) { fSharedClusterCut = value; }
78 void SetCrossedRowsCut(Int_t value) { fCrossedRowsCut = value; }
79 void SetFoundFractionCut(Double_t value) { fFoundFractionCut = value; }
80 void SetParticlePhiCutEventPlane(Double_t min, Double_t max) { fTrackPhiCutEvPlMin = min; fTrackPhiCutEvPlMax = max; }
83 Bool_t CheckTrack(AliVParticle * part);
86 Int_t fDebug; // debug flag
87 Int_t fFilterBit; // track selection cuts
88 UInt_t fTrackStatus; // if non-0, the bits set in this variable are required for each track
89 Bool_t fOnlyHadrons; // consider only charged Pions, Protons and Kaons
90 Bool_t fCheckMotherPDG; // Check the PDG code of mother for secondaries
91 Double_t fTrackEtaCut; // pseudo-rapidity limit of transverse regions
92 Double_t fTrackEtaCutMin; // minimum of the pseudo-rapidity limit of transverse regions
93 Double_t fTrackPhiCutEvPlMin; // Minimum Phi cut on particles with respect to the Event Plane (values between 0 and Pi/2)
94 Double_t fTrackPhiCutEvPlMax; // Maximum Phi cut on particles with respect to the Event Plane (values between 0 and Pi/2)
95 Double_t fTrackPtMin; // pt limit for selecting particles
96 UInt_t fEventSelection; // bit for physics selection
97 TFormula* fDCAXYCut; // additional pt dependent cut on DCA XY (only for AOD)
98 Double_t fSharedClusterCut; // cut on shared clusters (only for AOD)
99 Int_t fCrossedRowsCut; // cut on crossed rows (only for AOD)
100 Double_t fFoundFractionCut; // cut on crossed rows/findable clusters (only for AOD)
102 AliESDtrackCuts *fEsdTrackCuts; // set of cuts when reading ESD
103 AliESDtrackCuts *fEsdTrackCutsExtra1; // set of cuts when reading ESD
104 AliESDtrackCuts *fEsdTrackCutsExtra2; // set of cuts when reading ESD
105 AliHelperPID *fHelperPID; // PID Helper object
106 UInt_t fEventCounter; //! processed events (counter)
108 ClassDef(AliAnalyseLeadingTrackUE,0)