PWGLF/SPECTRA/MultEvShape/AliAnalysisCentralCutEvtMC.h

   1 #ifndef ALIANALYSISCENTRALCUTEVTMC_H
   2 #define ALIANALYSISCENTRALCUTEVTMC_H
   3
   4
   5 /*
   6  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
   7  * See cxx source for full Copyright notice
   8  * $Id$
   9  */
  10
  11 // ---------------------------------------------------
  12 //  MC event level cuts for azimuthal isotropic
  13 //  expansion in highly central collisions analysis
  14 //  author: Cristian Andrei
  15 //          acristian@niham.nipne.ro
  16 // ----------------------------------------------------
  17
  18
  19
  20 #include "AliAnalysisCuts.h"
  21
  22
  23 class TObject;
  24 class AliMCEvent;
  25 class AliStack;
  26
  27
  28 class AliAnalysisCentralCutEvtMC: public AliAnalysisCuts {
  29  public:
  30         AliAnalysisCentralCutEvtMC(const char *name="AliAnalysisCentralCutEvtMC", const char *title="MC_cuts");
  31     virtual ~AliAnalysisCentralCutEvtMC();
  32
  33     Bool_t  IsSelected(TObject* obj);
  34     Bool_t  IsSelected(TList* /*list*/) {return kTRUE;}
  35
  36     void SetMultiplicityRange(Int_t r1=0, Int_t r2=1000000){fReqMult = kTRUE; fMultMin=r1;  fMultMax=r2;}
  37     void SetDirectivityRange(Float_t r1=-1e10, Float_t r2=1e10) {fReqDir = kTRUE; fDirMin=r1; fDirMax=r2;}
  38     void SetDirUnitRange(Float_t r1=-1e10, Float_t r2=1e10) {fReqDirUnit = kTRUE; fDirUMin=r1; fDirUMax=r2;}
  39
  40
  41  private:
  42         AliAnalysisCentralCutEvtMC(const AliAnalysisCentralCutEvtMC& ref);
  43         AliAnalysisCentralCutEvtMC& operator=(const AliAnalysisCentralCutEvtMC& ref);
  44
  45     Bool_t fReqMult, fReqDir, fReqDirUnit; //set whether to compute multiplicity, directivity or dir unity
  46     Double_t fMultMin, fMultMax; //stores the multiplicity cut interval
  47     Double_t fDirMin, fDirMax; //stores the directivity cut interval
  48     Double_t fDirUMin, fDirUMax; //stores the directivity unity cut interval
  49
  50
  51     Int_t CalcMult(AliMCEvent* const mcEv);
  52     Double_t CalcDir(AliMCEvent* const mcEv);
  53     Double_t CalcDirUnit(AliMCEvent* const mcEv);
  54
  55
  56     ClassDef(AliAnalysisCentralCutEvtMC, 1);
  57 };
  58
  59 #endif