add cut to reject tracks close to TPC sector edge

[u/mrichter/AliRoot.git] / PWG / EMCAL / AliParticleContainer.h

diff --git a/PWG/EMCAL/AliParticleContainer.h b/PWG/EMCAL/AliParticleContainer.h
index 8a4a3c17df99e2c2148bc2038a7b93fb679f1cc4..683e47dec0f4dc4ac43fa664716d81ec962855f7 100644 (file)
--- a/PWG/EMCAL/AliParticleContainer.h
+++ b/PWG/EMCAL/AliParticleContainer.h
@@ -6,6 +6,8 @@
 class AliVEvent;
 class AliVParticle;
 
+#include "AliAODMCParticle.h"
+
 #include "AliEmcalContainer.h"
 
 class AliParticleContainer : public AliEmcalContainer {
@@ -14,13 +16,7 @@ class AliParticleContainer : public AliEmcalContainer {
   AliParticleContainer(const char *name); 
   virtual ~AliParticleContainer(){;}
 
-  void  SetParticlePtCut(Double_t cut)                    { fParticlePtCut = cut ; }
-  void  SetParticleEtaLimits(Double_t min, Double_t max)  { fParticleMaxEta = max ; fParticleMinEta = min ; }
-  void  SetParticlePhiLimits(Double_t min, Double_t max)  { fParticleMaxPhi = max ; fParticleMinPhi = min ; }
-  void  SetTrackBitMap(UInt_t m)                          { fTrackBitMap     = m ; }
-  void  SetMCTrackBitMap(UInt_t m)                        { fMCTrackBitMap   = m ; }
-  void  SetMinMCLabel(Int_t s)                            { fMinMCLabel      = s ; }
-
+  Bool_t                      AcceptParticle(AliVParticle         *vp)       ;
   Double_t                    GetParticlePtCut()                        const   { return fParticlePtCut; }
   Double_t                    GetParticleEtaMin()                       const   { return fParticleMinEta; }
   Double_t                    GetParticleEtaMax()                       const   { return fParticleMaxEta; }
@@ -28,16 +24,28 @@ class AliParticleContainer : public AliEmcalContainer {
   Double_t                    GetParticlePhiMax()                       const   { return fParticleMaxPhi; }
   AliVParticle               *GetLeadingParticle(const char* opt="")         ;
   AliVParticle               *GetParticle(Int_t i)                      const;
-  AliVParticle               *GetAcceptParticle(Int_t i)                const;
+  AliVParticle               *GetAcceptParticle(Int_t i)                     ;
   AliVParticle               *GetParticleWithLabel(Int_t lab)           const;
-  AliVParticle               *GetAcceptParticleWithLabel(Int_t lab)     const;
+  AliVParticle               *GetAcceptParticleWithLabel(Int_t lab)          ;
   AliVParticle               *GetNextAcceptParticle(Int_t i=-1)              ;
   AliVParticle               *GetNextParticle(Int_t i=-1)                    ;
   void                        GetMomentum(TLorentzVector &mom, Int_t i) const;
-  Bool_t                      AcceptParticle(AliVParticle         *vp)  const;
   Int_t                       GetNParticles()                           const   {return GetNEntries();}
   Int_t                       GetNAcceptedParticles()                   ;
   void                        SetClassName(const char *clname);
+  void                        SetMCTrackBitMap(UInt_t m)                        { fMCTrackBitMap   = m ; }
+  void                        SetMinMCLabel(Int_t s)                            { fMinMCLabel      = s ; }
+  void                        SetParticlePtCut(Double_t cut)                    { fParticlePtCut = cut ; }
+  void                        SetParticleEtaLimits(Double_t min, Double_t max)  { fParticleMaxEta = max ; fParticleMinEta = min ; }
+  void                        SetParticlePhiLimits(Double_t min, Double_t max)  { fParticleMaxPhi = max ; fParticleMinPhi = min ; }
+  void                        SetMinDistanceTPCSectorEdge(Double_t min)         { fMinDistanceTPCSectorEdge = min; }
+  void                        SetTrackBitMap(UInt_t m)                          { fTrackBitMap     = m ; }
+  void                        SetMCFlag(UInt_t m)                               { fMCFlag          = m ; }
+  void                        SelectHIJING(Bool_t s)                            { if (s) fGeneratorIndex = 0; else fGeneratorIndex = -1; }
+  void                        SetGeneratorIndex(Short_t i)                      { fGeneratorIndex = i  ; }
+  void                        SelectPhysicalPrimaries(Bool_t s)                 { if (s) fMCFlag |=  AliAODMCParticle::kPhysicalPrim ; 
+                                                                                  else   fMCFlag &= ~AliAODMCParticle::kPhysicalPrim ; }
+  void                        SetCharge(Short_t c)                              { fCharge = c         ; }
 
  protected:
   Double_t                    fParticlePtCut;                 // cut on particle pt
@@ -45,15 +53,19 @@ class AliParticleContainer : public AliEmcalContainer {
   Double_t                    fParticleMaxEta;                // cut on particle eta
   Double_t                    fParticleMinPhi;                // cut on particle phi
   Double_t                    fParticleMaxPhi;                // cut on particle phi
+  Double_t                    fMinDistanceTPCSectorEdge;      // require minimum distance to edge of TPC sector edge
   UInt_t                      fTrackBitMap;                   // bit map of accepted tracks (non MC)
   UInt_t                      fMCTrackBitMap;                 // bit map of accepted MC tracks
   Int_t                       fMinMCLabel;                    // minimum MC label value for the tracks/clusters being considered MC particles
+  UInt_t                      fMCFlag;                        // select MC particles with flags
+  Short_t                     fGeneratorIndex;                // select MC particles with generator index (default = -1 = switch off selection)
+  Short_t                     fCharge;                        // select particles with charge=fCharge
 
  private:
   AliParticleContainer(const AliParticleContainer& obj); // copy constructor
   AliParticleContainer& operator=(const AliParticleContainer& other); // assignment
 
-  ClassDef(AliParticleContainer,1);
+  ClassDef(AliParticleContainer,4);
 
 };