dNdPt analysis

[u/mrichter/AliRoot.git] / PWG0 / dNdPt / AlidNdPtEventCuts.h

#ifndef ALIDNDPTEVENTCUTS_H\r
#define ALIDNDPTEVENTCUTS_H\r
\r
//------------------------------------------------------------------------------\r
// Class to keep event selection cuts for dNdPt analysis. \r
// \r
// Author: J.Otwinowski 01/11/2008 \r
//------------------------------------------------------------------------------\r
\r
#include "AliAnalysisCuts.h"\r
#include "dNdPt/AlidNdPtHelper.h"\r
\r
class AliESDEvent;\r
class AliESDVertex;\r
class AliMCEvent;\r
class AliHeader;\r
class AliGenEventHeader;\r
\r
class AlidNdPtEventCuts : public AliAnalysisCuts\r
{\r
public:\r
  AlidNdPtEventCuts(const Char_t* name ="AlidNdPtEventCuts", const Char_t *title ="");\r
  virtual ~AlidNdPtEventCuts(); \r
 \r
  // setters \r
  void SetRecVertexRequired(const Bool_t bFlag=kTRUE)  {fRecVertexRequired=bFlag;}\r
  void SetEventProcessType(AlidNdPtHelper::MCProcessType type=AlidNdPtHelper::kInvalidProcess)  {fEventProcessType=type;}\r
  void SetNContributorsRange(const Float_t min=0.,const Float_t max=1e99) {fMinNContributors=min; fMaxNContributors=max;}\r
  void SetMaxR(const Float_t max=1e99) {fMaxR=max;}\r
  void SetZvRange(const Float_t min=-1e99, const Float_t max=1e99) {fMinZv=min; fMaxZv=max;}\r
\r
  void SetMeanXYZv(const Float_t xv=0.0, const Float_t yv=0.0, const Float_t zv=0.0) {\r
    fMeanXv = xv; fMeanYv = yv; fMeanZv = zv;\r
  }\r
\r
  void SetSigmaMeanXYZv(const Float_t sxv=1.0, const Float_t syv=1.0, const Float_t szv=10.0) {\r
    fSigmaMeanXv = sxv; fSigmaMeanYv = syv; fSigmaMeanZv = szv;\r
  }\r
\r
  // getters \r
  Bool_t  IsRecVertexRequired() const {return fRecVertexRequired;}\r
  Int_t   GetEventProcessType() const {return fEventProcessType;}  \r
  Float_t GetMinNContributors() const {return fMinNContributors;}\r
  Float_t GetMaxNContributors() const {return fMaxNContributors;}\r
  Float_t GetMaxR() const {return fMaxR;}\r
  Float_t GetMinZv() const {return fMinZv;}\r
  Float_t GetMaxZv() const {return fMaxZv;}\r
\r
  Float_t GetMeanXv() const {return fMeanXv;}\r
  Float_t GetMeanYv() const {return fMeanYv;}\r
  Float_t GetMeanZv() const {return fMeanZv;}\r
\r
  Float_t GetSigmaMeanXv() const {return fSigmaMeanXv;}\r
  Float_t GetSigmaMeanYv() const {return fSigmaMeanYv;}\r
  Float_t GetSigmaMeanZv() const {return fSigmaMeanZv;}\r
 \r
  // cuts init function\r
  void Init();\r
\r
  // check MC tracks\r
  Bool_t IsSelected(TObject *) {return kTRUE;}\r
  Bool_t IsSelected(TList *) {return kTRUE;}\r
\r
  // accept event\r
  Bool_t AcceptEvent(AliESDEvent *event=0, AliMCEvent *mcEvent=0, const AliESDVertex *vtx=0);\r
  Bool_t AcceptMCEvent(AliMCEvent *mcEvent=0);\r
\r
  // Merge output objects (needed by PROOF) \r
  virtual Long64_t Merge(TCollection* list);\r
\r
private:\r
  Bool_t fRecVertexRequired; // reconstructed event vertex required  \r
  Int_t fEventProcessType;   // select MC event process type (ND, SD, DD)\r
  Float_t fMinNContributors; // min. number of contributing vertex tracks\r
  Float_t fMaxNContributors; // max. number of contributing vertex tracks\r
  Float_t fMaxR;             // max. vertex radii (R = sqrt(Xv^2+Yv^2) \r
  Float_t fMinZv;            // min. Zv vertex\r
  Float_t fMaxZv;            // max. Zv vertex\r
\r
  // interaction spot constraint\r
  Float_t fMeanXv; // mean Xv position\r
  Float_t fMeanYv; // mean Yv position\r
  Float_t fMeanZv; // mean Zv position\r
\r
  Float_t fSigmaMeanXv; // sigma mean Xv position \r
  Float_t fSigmaMeanYv; // sigma mean Yv position\r
  Float_t fSigmaMeanZv; // sigma mean Zv position\r
 \r
  AlidNdPtEventCuts(const AlidNdPtEventCuts&); // not implemented\r
  AlidNdPtEventCuts& operator=(const AlidNdPtEventCuts&); // not implemented\r
\r
  ClassDef(AlidNdPtEventCuts, 1)\r
};\r
\r
#endif // ALIDNDPTEVENTCUTS_H\r