1 #ifndef ALIANALYSISETSELECTOR_H
2 #define ALIANALYSISETSELECTOR_H
3 //_________________________________________________________________________
4 // Utility Class for transverse energy studies
8 //*-- Authors: Oystein Djuvsland (Bergen)
9 //_________________________________________________________________________
11 #include "AliAnalysisEtCommon.h"
12 #include "AliESDEvent.h"
14 class AliESDCaloCluster;
17 class AliAnalysisEtCuts;
22 class AliAnalysisEtSelector : public AliAnalysisEtCommon
27 // Constructor takes cuts object
28 AliAnalysisEtSelector(AliAnalysisEtCuts *cuts);
29 AliAnalysisEtSelector();
32 virtual ~AliAnalysisEtSelector();
34 // Set the current event
35 virtual void SetEvent(const AliESDEvent *event);
38 virtual void Init() {}
41 virtual Int_t Init(const AliESDEvent *event) { fRunNumber = event->GetRunNumber(); return 0; }
43 // Return CaloClusters for the detector
44 virtual TRefArray* GetClusters() { return 0; }
46 virtual Float_t ShiftAngle(Float_t phi);//Always returns an angle in radians between -pi<phi<pi
48 // Return true if cluster has energy > cut
49 virtual Bool_t PassMinEnergyCut(const AliESDCaloCluster &/*cluster*/) const { return true; }
51 // Return true if cluster has energy > cut
52 virtual Bool_t PassMinEnergyCut(const TParticle &/*part*/) const { return true; }
54 // Cut on distance to bad channel
55 virtual Bool_t PassDistanceToBadChannelCut(const AliESDCaloCluster &/*cluster*/) const { return true; }
57 // Cut on track matching
58 virtual Bool_t PassTrackMatchingCut(const AliESDCaloCluster &/*cluster*/) const { return true; }
60 // Cut on neutral monte carlo particle
61 virtual Bool_t IsNeutralMcParticle(Int_t pIdx, AliStack& s, const TParticlePDG& pdg) const;
63 // Is it an EM E_T particle
64 virtual Bool_t IsEmEtParticle(const Int_t pdgCode) const;
66 // Does the particle come from an EM E_T primary ?
67 virtual Bool_t PrimaryIsEmEtParticle(const Int_t pIdx, AliStack &stack) const;
69 // Get the index of primary particle for the particle
70 Int_t GetPrimary(const Int_t partIdx, AliStack &stack) const;
72 // Cut on geometrical acceptance
73 virtual Bool_t CutGeometricalAcceptance(const TParticle &/*part*/) { return true; }
75 // Cut on geometrical acceptance
76 virtual Bool_t CutGeometricalAcceptance(const AliVTrack &/*part*/) { return true; }
78 // Cut on geometrical acceptance
79 virtual Bool_t CutGeometricalAcceptance(const AliESDCaloCluster &/*part*/) { return true; }
81 // From secondary vertex?
82 //virtual Bool_t FromSecondaryInteraction(const TParticle& part, AliStack& stack) const;
83 virtual Bool_t FromSecondaryInteraction(Int_t partID, AliStack& stack) const;
85 Int_t GetMother(Int_t partID, AliStack& stack) const;
86 Bool_t IsFromDetectorCover(Int_t partID, AliStack& stack) const;
87 Int_t GetFirstMotherNotFromDetectorCover(Int_t partID, AliStack& stack) const;
89 // Cluster is in correct detector
90 virtual Bool_t IsDetectorCluster(const AliESDCaloCluster &cluster) const = 0;
92 // Get correct cluster label - PHOS needs different method
93 virtual UInt_t GetLabel(const AliESDCaloCluster *cluster, AliStack *stack){if(!stack){return 0;}else{return TMath::Abs(cluster->GetLabel());}}
97 const AliVEvent *fEvent; // Pointer to current event
99 TRefArray *fClusterArray; // Array of clusters
101 AliAnalysisEtCuts *fCuts; // Pointer to the cuts object; DS: also in base class?
103 Bool_t SuspiciousDecayInChain(const UInt_t suspectMotherPdg, const UInt_t suspectDaughterPdg, const TParticle& part, AliStack& stack) const;
107 Bool_t fInitialized; // matrix initialized
111 //AliAnalysisEtSelector(); // Prohibited
112 AliAnalysisEtSelector(const AliAnalysisEtSelector& other);// Prohibited
113 AliAnalysisEtSelector& operator=(const AliAnalysisEtSelector& other);// Prohibited
114 bool operator==(const AliAnalysisEtSelector& other) const;// Prohibited
116 ClassDef(AliAnalysisEtSelector, 1);
119 #endif // ALIANALYSISETSELECTOR_H