CORRFW/AliCFEventRecCuts.h

   1 #ifndef ALICFEVENTRECCUTS_H
   2 #define ALICFEVENTRECCUTS_H
   3 /**************************************************************************
   4  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  *                                                                        *
   6  * Author: The ALICE Off-line Project.                                    *
   7  * Contributors are mentioned in the code where appropriate.              *
   8  *                                                                        *
   9  * Permission to use, copy, modify and distribute this software and its   *
  10  * documentation strictly for non-commercial purposes is hereby granted   *
  11  * without fee, provided that the above copyright notice appears in all   *
  12  * copies and that both the copyright notice and this permission notice   *
  13  * appear in the supporting documentation. The authors make no claims     *
  14  * about the suitability of this software for any purpose. It is          *
  15  * provided "as is" without express or implied warranty.                  *
  16  **************************************************************************/
  17 // Cut on the Event at reconstructed level: for the moment
  18 // just the requirements on the number of charged tracks and on
  19 // the vertex position and resolution are implemented
  20 // The argument of IsSelected member function (passed object) is cast into
  21 // an AliESDEvent. In the future may be modified to use AliVEvent interface
  22 // and include more cut variables.
  23 // The class derives from AliCFCutBase
  24 // Author:S.Arcelli Silvia.Arcelli@cern.ch
  25
  26
  27 #include "AliCFCutBase.h"
  28 class TH1F;
  29 class TBits;
  30 //_____________________________________________________________________________
  31 class AliCFEventRecCuts: public AliCFCutBase
  32 {
  33  public :
  34   AliCFEventRecCuts() ;
  35   AliCFEventRecCuts(Char_t* name, Char_t* title) ;
  36   AliCFEventRecCuts(const AliCFEventRecCuts& c) ;
  37   AliCFEventRecCuts& operator=(const AliCFEventRecCuts& c) ;
  38   ~AliCFEventRecCuts();
  39   void GetBitMap(TObject *obj, TBits*bitmap);
  40   Bool_t IsSelected(TObject* obj);
  41   void Init();
  42   void AddQAHistograms(TList *list) const;
  43
  44   void SetNTracksCut(Int_t xMin=-1, Int_t xMax=1000000) {fNTracksMin=xMin; fNTracksMax=xMax;} // cut values setter
  45
  46   void SetRequireVtxCuts(Bool_t vtx=kFALSE) {fRequireVtxCuts=vtx;} // cut values setter
  47   void SetVertexXCut(Double_t xMin=-1.e99, Double_t xMax=1.e99) { fVtxXMin=xMin; fVtxXMax=xMax;} // cut values setter
  48   void SetVertexYCut(Double_t yMin=-1.e99, Double_t yMax=1.e99) { fVtxYMin=yMin; fVtxYMax=yMax;} // cut values setter
  49   void SetVertexZCut(Double_t zMin=-1.e99, Double_t zMax=1.e99) { fVtxZMin=zMin; fVtxZMax=zMax;} // cut values setter
  50
  51   void SetVertexXResCut(Double_t xMax=1.e99) {fVtxXResMax=xMax;} // cut values setter
  52   void SetVertexYResCut(Double_t yMax=1.e99){fVtxYResMax=yMax;} // cut values setter
  53   void SetVertexZResCut(Double_t zMax=1.e99){fVtxZResMax=zMax;} // cut values setter
  54
  55   Int_t    GetNTracksMin() const {return fNTracksMin;} // cut values getter
  56   Int_t    GetNTracksMax() const {return fNTracksMax;} // cut values getter
  57   Bool_t   GetRequireVtxCuts() const {return fRequireVtxCuts;} // cut value getter
  58   Double_t GetVertexXMax() const {return fVtxXMax;} // cut values getter
  59   Double_t GetVertexYMax() const {return fVtxYMax;} // cut values getter
  60   Double_t GetVertexZMax() const {return fVtxZMax;} // cut values getter
  61   Double_t GetVertexXMin() const {return fVtxXMin;} // cut values getter
  62   Double_t GetVertexYMin() const {return fVtxYMin;} // cut values getter
  63   Double_t GetVertexZMin() const {return fVtxZMin;} // cut values getter
  64   Double_t GetVertexXResMax() const {return fVtxXResMax;} // cut values getter
  65   Double_t GetVertexYResMax() const {return fVtxYResMax;} // cut values getter
  66   Double_t GetVertexZResMax() const {return fVtxZResMax;} // cut values getter
  67
  68   // QA histograms
  69   void FillHistogramsBeforeCuts(TObject* obj) {return FillHistograms(obj,kFALSE);}
  70   void FillHistogramsAfterCuts(TObject* obj)  {return FillHistograms(obj,kTRUE);}
  71   // QA histogram setter
  72   // please use indices from the enumeration below
  73   void SetHistogramBins(Int_t index, Int_t nbins, Double_t *bins);
  74   void SetHistogramBins(Int_t index, Int_t nbins, Double_t xmin, Double_t xmax);
  75   enum{kNTracks=0,
  76          kVtxPosX,
  77          kVtxPosY,
  78          kVtxPosZ,
  79          kVtxResX,
  80          kVtxResY,
  81          kVtxResZ,
  82          kNCuts=7,
  83          kNStepQA=2
  84          };
  85
  86  private:
  87   TBits *SelectionBitMap(TObject* obj);
  88   void DefineHistograms();              // books histograms
  89   void Initialise();                    // sets everything to 0
  90   void FillHistograms(TObject* obj, Bool_t b);
  91   Int_t fNTracksMin; //minimum number of esd tracks
  92   Int_t fNTracksMax; //maximum number of esd tracks
  93   Bool_t fRequireVtxCuts ; //The type of trigger to be checked
  94   Double_t fVtxXMax ; //X vertex position, maximum value
  95   Double_t fVtxYMax ; //Y vertex position, maximum value
  96   Double_t fVtxZMax ; //Z vertex position, maximum value
  97   Double_t fVtxXMin ; //X vertex position, minimum value
  98   Double_t fVtxYMin ; //Y vertex position, minimum value
  99   Double_t fVtxZMin ; //Z vertex position, minimum value
 100   Double_t fVtxXResMax ;//Maximum value of sigma_vtx in X
 101   Double_t fVtxYResMax ;//Maximum value of sigma_vtx in X
 102   Double_t fVtxZResMax ;//Maximum value of sigma_vtx in X
 103
 104
 105   TBits *fBitMap ; //cut mask
 106
 107
 108   TH1F* fhQA[kNCuts][kNStepQA];         // QA Histograms
 109   //QA Histogram parameters
 110   Int_t fhNBinsNTracks;//size of array of bin limits, N Tracks ESD
 111   Double_t *fhBinLimNTracks;//[fhNBinsNTracks] bin limits, N Tracks ESD
 112   Int_t fhNBinsVtxPosX;//size of array of bin limits, Vtx Pos X
 113   Double_t *fhBinLimVtxPosX;//[fhNBinsVtxPosX] bin limits, Vtx Pos X
 114   Int_t fhNBinsVtxPosY;//size of array of bin limits, Vtx Pos Y
 115   Double_t *fhBinLimVtxPosY;//[fhNBinsVtxPosY] bin limits, Vtx Pos Y
 116   Int_t fhNBinsVtxPosZ;//size of array of bin limits, Vtx Pos Z
 117   Double_t *fhBinLimVtxPosZ;//[fhNBinsVtxPosZ] bin limits, Vtx Pos Z
 118   Int_t fhNBinsVtxResX;//size of array of bin limits, Vtx Res X
 119   Double_t *fhBinLimVtxResX;//[fhNBinsVtxResX] bin limits, Vtx Res X
 120   Int_t fhNBinsVtxResY;//size of array of bin limits, Vtx Res Y
 121   Double_t *fhBinLimVtxResY;//[fhNBinsVtxResY] bin limits, Vtx Res Y
 122   Int_t fhNBinsVtxResZ;//size of array of bin limits, Vtx Res Z
 123   Double_t *fhBinLimVtxResZ;//[fhNBinsVtxResZ] bin limits, Vtx Res Z
 124
 125   ClassDef(AliCFEventRecCuts,2);
 126 };
 127
 128 #endif