]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - HBTAN/AliHBTAnalysis.h
changes for proper protection against failed retrieval of CDB Reco object (moved...
[u/mrichter/AliRoot.git] / HBTAN / AliHBTAnalysis.h
index 175165b97fef3d3a6d2325e67676180df173494c..7f0782738e58ca2890cb019a69cd440c6d0b7747 100644 (file)
 // ALICE Off-Line framework:
 //
 // Piotr.Skowronski@cern.ch
-// more info: http://alisoft.cern.ch/people/skowron/analyzer/index.html
+// more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
 //
 ////////////////////////////////////////////////////////////////////////////
 //_________________________________________________________
 
 #include <AliAnalysis.h>
-#include "AliAODPairCut.h"
-#include "AliAODParticleCut.h"
 
 class AliHBTCut;
 //class AliHBTPair;
@@ -31,6 +29,7 @@ class AliHBTRun;
 class AliAOD;
 class AliReader;
 class AliEventBuffer;
+class AliEventCut;
 
 class AliHBTOnePairFctn;      
 class AliHBTTwoPairFctn;
@@ -49,52 +48,45 @@ class AliHBTAnalysis: public AliAnalysis
      virtual ~AliHBTAnalysis();
 
      Int_t  Init();
-     Int_t ProcessEvent(AliAOD* aodrec, AliAOD* aodsim = 0x0);
-     Int_t Finish();
+     Int_t  ProcessEvent(AliAOD* aodrec, AliAOD* aodsim = 0x0);
+     Int_t  Finish();
      
-     enum EProcessOption{kReconstructed,kSimulated,kSimulatedAndReconstructed};
-     void  SetProcessOption(EProcessOption option){fProcessOption = option;}//Init must be called to make effect
+     enum   EProcessOption{kReconstructed,kSimulated,kSimulatedAndReconstructed};
+     void   SetProcessOption(EProcessOption option){fProcessOption = option;}//Init must be called to make effect
      
-     virtual void Process(Option_t* option = "TracksAndParticles");//Stand alone HBT analysis 
+     void   Process(Option_t* option = "TracksAndParticles");//Stand alone HBT analysis 
      
-     void SetGlobalPairCut(AliAODPairCut* cut);
+     void   SetGlobalPairCut(AliAODPairCut* cut);
      
-     void AddTrackFunction(AliHBTOnePairFctn* f);
-     void AddParticleFunction(AliHBTOnePairFctn* f);
-     void AddParticleAndTrackFunction(AliHBTTwoPairFctn* f);
+     void   AddTrackFunction(AliHBTOnePairFctn* f);
+     void   AddParticleFunction(AliHBTOnePairFctn* f);
+     void   AddParticleAndTrackFunction(AliHBTTwoPairFctn* f);
      
-     void AddParticleMonitorFunction(AliHBTMonOneParticleFctn* f);    //z.ch.
-     void AddTrackMonitorFunction(AliHBTMonOneParticleFctn* f);    //z.ch.
-     void AddParticleAndTrackMonitorFunction(AliHBTMonTwoParticleFctn* f);//z.ch.
+     void   AddParticleMonitorFunction(AliHBTMonOneParticleFctn* f);    //z.ch.
+     void   AddTrackMonitorFunction(AliHBTMonOneParticleFctn* f);    //z.ch.
+     void   AddParticleAndTrackMonitorFunction(AliHBTMonTwoParticleFctn* f);//z.ch.
 
-     void AddResolutionFunction(AliHBTTwoPairFctn* f){AddParticleAndTrackFunction(f);}
+     void   AddResolutionFunction(AliHBTTwoPairFctn* f){AddParticleAndTrackFunction(f);}
      
-     void SetReader(AliReader* r){fReader = r;}
+     void   SetReader(AliReader* r){fReader = r;}
      
-     void WriteFunctions();
-     
-     void SetBufferSize(Int_t buffsize){fBufferSize=buffsize;}
-     void SetOwner(Bool_t owner=kTRUE){fIsOwner=owner;}
+     void   WriteFunctions();
+     void   SetOutputFileName(const char* fname);
+          
+     void   SetBufferSize(Int_t buffsize){fBufferSize=buffsize;}
+     void   SetOwner(Bool_t owner=kTRUE){fIsOwner=owner;}
      Bool_t IsOwner() const {return fIsOwner;}
      Bool_t IsNonIdentAnalysis();
      void   ResetFunctions();
      void   SetDisplayInfo(Int_t howoften){fDisplayMixingInfo = howoften;}//defines every each line info about mixing is displayed
      
-     void   SetCutsOnParticles(); // -- aplies only to Process Tracks And Particles
-     void   SetCutsOnTracks();// -- aplies only to Process Tracks And Particles
-     void   SetCutsOnTracksAndParticles();// Default // -- aplies only to Process Tracks And Particles
+     void   SetApparentVertex(Double_t x, Double_t y, Double_t z);//Sets apparent vertex
      
      static void PressAnyKey();//small utility function that helps to make comfortable macros
-   protected:
-     
-     Bool_t RunCoherencyCheck();
+
      
-     void FilterOut(AliAOD* outpart1, AliAOD* outpart2, AliAOD* inpart,
-                    AliAOD* outtrack1, AliAOD* outtrack2, AliAOD* intrack)const;
-     void FilterOut(AliAOD* out1, AliAOD* out2, AliAOD* in)const;
-     void DeleteFunctions();
+   protected:
      
-
      /**********************************************/
      /*      E V E N T     P R O C E S S I N G     */
      /**********************************************/
@@ -112,14 +104,21 @@ class AliHBTAnalysis: public AliAnalysis
 
      // OLD legacy schema
 
-     virtual void ProcessTracks();
-     virtual void ProcessParticles();
-     virtual void ProcessTracksAndParticles();
+     void   ProcessTracks();
+     void   ProcessParticles();
+     void   ProcessTracksAndParticles();
      
-     virtual void ProcessTracksAndParticlesNonIdentAnal();
-     virtual void ProcessParticlesNonIdentAnal();
-     virtual void ProcessTracksNonIdentAnal();
+     void   ProcessTracksAndParticlesNonIdentAnal();
+     void   ProcessParticlesNonIdentAnal();
+     void   ProcessTracksNonIdentAnal();
 
+     Bool_t RunCoherencyCheck();
+     
+     void   FilterOut(AliAOD* outpart1, AliAOD* outpart2, AliAOD* inpart,
+                      AliAOD* outtrack1, AliAOD* outtrack2, AliAOD* intrack)const;
+     void   FilterOut(AliAOD* out1, AliAOD* out2, AliAOD* in)const;
+     void   DeleteFunctions();
+     
 
      AliReader* fReader;//! Pointer to reader
      
@@ -140,48 +139,28 @@ class AliHBTAnalysis: public AliAnalysis
      AliHBTMonOneParticleFctn**  fTrackMonitorFunctions; //! which are used for single particle analysis,
      AliHBTMonTwoParticleFctn**  fParticleAndTrackMonitorFunctions;  //! cut monitoring, etc.
 
+     AliEventCut*     fBkgEventCut;// We can narrow class of events used in 
+
+     AliEventBuffer*  fPartBuffer;//Sim Particles event buffer
+     AliEventBuffer*  fTrackBuffer;//Rec Tracks event buffer
 
      /**********************************************/
-     /* Control parameters  */
+     /* Control parameters                         */
      /**********************************************/
-
-     AliAODPairCut*   fPairCut;//! Pair cut applied for all mixed particles
       
-     Int_t  fBufferSize; //!defines the size of buffer for mixed events; -1==MIX All
-     Int_t  fDisplayMixingInfo;//!defines every which particle mixing info is displayed
-     Bool_t fIsOwner;//!defines of all functions are supposed to be deleted while by the way of analysis defaulr false
-    
-     AliEventBuffer* fPartBuffer;//Sim Particles event buffer
-     AliEventBuffer* fTrackBuffer;//Rec Tracks event buffer
+     Int_t            fBufferSize; //defines the size of buffer for mixed events; -1==MIX All
+     Int_t            fDisplayMixingInfo;//!defines every which particle mixing info is displayed
+     Bool_t           fIsOwner;//!defines of all functions are supposed to be deleted while by the way of analysis defaulr false
      
-     EProcessOption  fProcessOption;//Option that says waht analysis to do (Rec, Sim or SimAndRec)
-     Bool_t          fNoCorrfctns;//Internal flag indicating that no cfs are set by the user (only monitor ones)
+     EProcessOption   fProcessOption;//Option that says waht analysis to do (Rec, Sim or SimAndRec)
+     Bool_t           fNoCorrfctns;//Internal flag indicating that no cfs are set by the user (only monitor ones)
+     TString*         fOutputFileName;//Fiele name where to dump results, if not specified reults are written to gDirectory
+     
+     Double_t         fVertexX;//X position of apparent vertex
+     Double_t         fVertexY;//Y position of apparent vertex
+     Double_t         fVertexZ;//Z position of apparent vertex
      
    private:
-   
-     /**********************************************/
-     /*                C U T S                     */
-     /**********************************************/
-   
-     Bool_t (AliHBTAnalysis::*fkPass)(AliAODPair* partpair, AliAODPair* trackpair) const;//Pointer to function that performes pair cut
-     Bool_t (AliHBTAnalysis::*fkPass1)(AliVAODParticle* partpair, AliVAODParticle* trackpair) const;//Pointer to function that performes cut on first particle
-     Bool_t (AliHBTAnalysis::*fkPass2)(AliVAODParticle* partpair, AliVAODParticle* trackpair) const;//Pointer to function that performes cut on second particle
-     Bool_t (AliHBTAnalysis::*fkPassPairProp)(AliAODPair* partpair, AliAODPair* trackpair) const;//Pointer to function that performes pair cut
-     
-     Bool_t PassPartAndTrack (AliAODPair* partpair, AliAODPair* trackpair) const {return (fPairCut->Pass((AliAODPair*)partpair))?kTRUE:fPairCut->Pass((AliAODPair*)trackpair);}
-     Bool_t PassPartAndTrack1(AliVAODParticle* part, AliVAODParticle* track) const;
-     Bool_t PassPartAndTrack2(AliVAODParticle* part, AliVAODParticle* track) const;
-     Bool_t PassPairPropPartAndTrack (AliAODPair* partpair, AliAODPair* trackpair) const {return (fPairCut->PassPairProp((AliAODPair*)partpair))?kTRUE:fPairCut->PassPairProp((AliAODPair*)trackpair);}
-     
-     Bool_t PassPart (AliAODPair* partpair, AliAODPair* /*trackpair*/) const {return fPairCut->Pass((AliAODPair*)partpair);}
-     Bool_t PassPart1(AliVAODParticle* part, AliVAODParticle* /*track*/) const {return fPairCut->GetFirstPartCut()->Pass(part);}
-     Bool_t PassPart2(AliVAODParticle* part, AliVAODParticle* /*track*/) const {return fPairCut->GetSecondPartCut()->Pass(part);}
-     Bool_t PassPairPropPart (AliAODPair* partpair, AliAODPair* /*trackpair*/) const {return fPairCut->PassPairProp((AliAODPair*)partpair);}
-     
-     Bool_t PassTrack (AliAODPair* /*partpair*/, AliAODPair* trackpair) const {return fPairCut->Pass((AliAODPair*)trackpair);}
-     Bool_t PassTrack1(AliVAODParticle* /*part*/, AliVAODParticle* track) const {return fPairCut->GetFirstPartCut()->Pass(track);}
-     Bool_t PassTrack2(AliVAODParticle* /*part*/, AliVAODParticle* track) const {return fPairCut->GetSecondPartCut()->Pass(track);}
-     Bool_t PassPairPropTrack (AliAODPair* /*partpair*/, AliAODPair* trackpair) const {return fPairCut->PassPairProp((AliAODPair*)trackpair);}
 
      static const UInt_t fgkFctnArraySize;//!
      static const UInt_t fgkDefaultMixingInfo;//!
@@ -189,20 +168,4 @@ class AliHBTAnalysis: public AliAnalysis
 
      ClassDef(AliHBTAnalysis,0)
  };
-inline Bool_t AliHBTAnalysis::PassPartAndTrack1(AliVAODParticle* part,AliVAODParticle* track) const
-{
-//Checks first particle from both, particle and track pairs
-  AliAODParticleCut* pc = fPairCut->GetFirstPartCut();
-  return (pc->Pass(part))?kTRUE:pc->Pass(track);
-}
-/*************************************************************************************/ 
-inline Bool_t AliHBTAnalysis::PassPartAndTrack2(AliVAODParticle* part,AliVAODParticle* track) const
-{
-//Checks second particle from both, particle and track pairs
-  AliAODParticleCut* pc = fPairCut->GetSecondPartCut();
-  return (pc->Pass(part))?kTRUE:pc->Pass(track);
-}
-/*************************************************************************************/ 
 #endif