AliRDHFCutsOmegactoeleOmegafromAODtracks.h

Go to the documentation of this file.

#ifndef ALIRDHFCUTSOMEGACTOELEOMEGAFROMAODTRACKS_H
#define ALIRDHFCUTSOMEGACTOELEOMEGAFROMAODTRACKS_H
/* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

//***********************************************************
//***********************************************************

#include "AliRDHFCuts.h"

class AliRDHFCutsOmegactoeleOmegafromAODtracks : public AliRDHFCuts
{
 public:

  enum EPIDStrategy{
    kNSigmaCuts,
    kNSigmaCustomizedCuts,
    kCombinedCuts
  };

  AliRDHFCutsOmegactoeleOmegafromAODtracks(const char* name="CutsOmegactoeleOmega");
  virtual ~AliRDHFCutsOmegactoeleOmegafromAODtracks();
  AliRDHFCutsOmegactoeleOmegafromAODtracks(const AliRDHFCutsOmegactoeleOmegafromAODtracks& source);
  AliRDHFCutsOmegactoeleOmegafromAODtracks& operator=(const AliRDHFCutsOmegactoeleOmegafromAODtracks& source);

  using AliRDHFCuts::GetCutVarsForOpt;
  virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters);

  using AliRDHFCuts::IsSelected;
  virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel);
  using AliRDHFCuts::IsSelectedPID;
  virtual Int_t IsSelectedPID(AliAODRecoDecayHF* obj);
  Int_t IsSelectedCombinedPID(AliAODRecoDecayHF* obj);
  Bool_t IsSelectedeID(AliAODTrack* trk);
  Bool_t IsSelectedCustomizedeID(AliAODTrack* trk);
  Bool_t IsSelectedCombinedeID(AliAODTrack* trk);

  void SetPIDStrategy(EPIDStrategy pidStrategy){fPIDStrategy=pidStrategy;}
  EPIDStrategy GetPIDStrategy() const {return fPIDStrategy;}
  void SetCombinedPIDThreshold(Double_t a){fCombinedPIDThreshold=a;}
  Double_t GetCombinedPIDThreshold(){return fCombinedPIDThreshold;}


  Bool_t SingleTrkCuts(AliAODTrack *trk, AliAODVertex *primvert);
  Bool_t SingleCascadeCuts(AliAODcascade *casc, Double_t *vert);
  Bool_t SelectWithRoughCuts(AliAODcascade *casc, AliAODTrack *trk1);

  void SetProdTrackTPCNclsPIDMin(Int_t a){fProdTrackTPCNclsPIDMin=a;}
  void SetProdTrackTPCNclsRatioMin(Double_t a){fProdTrackTPCNclsRatioMin=a;}
  void SetProdUseAODFilterBit(Bool_t a){fProdUseAODFilterBit=a;}
  void SetProdMassTolLambda(Double_t a){fProdMassTolLambda=a;}
  void SetProdMassTolOmega(Double_t a){fProdMassTolOmega=a;}
  void SetProdMassRejXi(Double_t a){fProdMassRejXi=a;}
  void SetProdRfidMinV0(Double_t a){fProdRfidMinV0=a;}
  void SetProdRfidMaxV0(Double_t a){fProdRfidMaxV0=a;}
  void SetProdRfidMinOmega(Double_t a){fProdRfidMinOmega=a;}
  void SetProdRfidMaxOmega(Double_t a){fProdRfidMaxOmega=a;}
  void SetProdCascProperDecayLengthMax(Double_t a){fProdCascProperDecayLengthMax=a;}
  void SetProdDcaOmegaDaughtersMax(Double_t a){fProdDcaOmegaDaughtersMax=a;}
  void SetProdDcaV0DaughtersMax(Double_t a){fProdDcaV0DaughtersMax=a;}
  void SetProdDcaBachToPrimVertexMin(Double_t a){fProdDcaBachToPrimVertexMin=a;}
  void SetProdDcaV0ToPrimVertexMin(Double_t a){fProdDcaV0ToPrimVertexMin=a;}
  void SetProdDcaV0PrToPrimVertexMin(Double_t a){fProdDcaV0PrToPrimVertexMin=a;}
  void SetProdDcaV0PiToPrimVertexMin(Double_t a){fProdDcaV0PiToPrimVertexMin=a;}
  void SetProdXiCosineOfPoiningAngleMin(Double_t a){fProdXiCosineOfPoiningAngleMin=a;}
  void SetProdV0CosineOfPoiningAngleXiMin(Double_t a){fProdV0CosineOfPoiningAngleXiMin=a;}
  void SetProdCascNTPCClustersMin(Double_t a){fProdCascNTPCClustersMin=a;}
  void SetProdCascEtaRange(Double_t a, Double_t b){fProdCascEtaMin=a;fProdCascEtaMax=b;}
  void SetProdCascRapRange(Double_t a, Double_t b){fProdCascRapMin=a;fProdCascRapMax=b;}
  void SetProdRoughMassTol(Double_t a){fProdRoughMassTol=a;}
  void SetProdRoughPtMin(Double_t a){fProdRoughPtMin=a;}
  
  Int_t GetProdTrackTPCNclsPIDMin(){return fProdTrackTPCNclsPIDMin;}
  Double_t GetProdTrackTPCNclsRatioMin(){return fProdTrackTPCNclsRatioMin;}
  Bool_t   GetProdUseAODFilterBit(){return fProdUseAODFilterBit;}
  Double_t GetProdMassTolLambda(){return fProdMassTolLambda;}
  Double_t GetProdMassTolOmega(){return fProdMassTolOmega;}
  Double_t GetProdMassRejXi(){return fProdMassRejXi;}
  Double_t GetProdRfidMinV0(){return fProdRfidMinV0;}
  Double_t GetProdRfidMaxV0(){return fProdRfidMaxV0;}
  Double_t GetProdRfidMinOmega(){return fProdRfidMinOmega;}
  Double_t GetProdRfidMaxOmega(){return fProdRfidMaxOmega;}
  Double_t GetProdCascProperDecayLengthMax(){return fProdCascProperDecayLengthMax;}
  Double_t GetProdDcaOmegaDaughtersMax(){return fProdDcaOmegaDaughtersMax;}
  Double_t GetProdDcaV0DaughtersMax(){return fProdDcaV0DaughtersMax;}
  Double_t GetProdDcaBachToPrimVertexMin(){return fProdDcaBachToPrimVertexMin;}
  Double_t GetProdDcaV0ToPrimVertexMin(){return fProdDcaV0ToPrimVertexMin;}
  Double_t GetProdDcaV0PrToPrimVertexMin(){return fProdDcaV0PrToPrimVertexMin;}
  Double_t GetProdDcaV0PiToPrimVertexMin(){return fProdDcaV0PiToPrimVertexMin;}
  Double_t GetProdXiCosineOfPoiningAngleMin(){return fProdXiCosineOfPoiningAngleMin;}
  Double_t GetProdV0CosineOfPoiningAngleXiMin(){return fProdV0CosineOfPoiningAngleXiMin;}
  Double_t GetProdRoughMassTol(){return fProdRoughMassTol;}
  Double_t GetProdRoughPtMin(){return fProdRoughPtMin;}
  Double_t GetProdCascNTPCClustersMin(){return fProdCascNTPCClustersMin;}
  void GetProdCascEtaRange(Double_t &a, Double_t &b){a=fProdCascEtaMin;b=fProdCascEtaMax;}
  void GetProdCascRapRange(Double_t &a, Double_t &b){a=fProdCascRapMin;b=fProdCascRapMax;}

  Bool_t GetUseCascadePID(){return fUseCascadePID;}
  void SetUseCascadePID(Bool_t a){fUseCascadePID=a;}
  void SetPidCascPi(AliAODPidHF* pidPion) { 
      if(fPidObjCascPi) delete fPidObjCascPi;
      fPidObjCascPi=new AliAODPidHF(*pidPion);
      }
  AliAODPidHF* GetPidCascPi() const {return fPidObjCascPi;}
  void SetPidCascPr(AliAODPidHF* pidProton) { 
      if(fPidObjCascPr) delete fPidObjCascPr;
      fPidObjCascPr=new AliAODPidHF(*pidProton);
      }
  AliAODPidHF* GetPidCascPr() const {return fPidObjCascPr;}
  void SetPidCascKa(AliAODPidHF* pidKaon) { 
      if(fPidObjCascKa) delete fPidObjCascKa;
      fPidObjCascKa=new AliAODPidHF(*pidKaon);
      }
  AliAODPidHF* GetPidCascKa() const {return fPidObjCascKa;}

  void SetExcludePionTPC(Bool_t a){fExcludePionTPC=a;}
  void SetExcludeProtonTPC(Bool_t a){fExcludeProtonTPC=a;}
  void SetExcludeKaonTPC(Bool_t a){fExcludeKaonTPC=a;}
  void SetExcludenSigmaPionTPC(Double_t a){fExcludenSigmaPionTPC=a;}
  void SetExcludenSigmaProtonTPC(Double_t a){fExcludenSigmaProtonTPC=a;}
  void SetExcludenSigmaKaonTPC(Double_t a){fExcludenSigmaKaonTPC=a;}
  void SetSigmaElectronTPCRange(Double_t a,Double_t b){fSigmaElectronTPCMin=a;fSigmaElectronTPCMax=b;}
  void SetSigmaElectronTOFRange(Double_t a,Double_t b){fSigmaElectronTOFMin=a;fSigmaElectronTOFMax=b;}
  

 protected:
  
 private:

  EPIDStrategy fPIDStrategy;        
  Double_t fCombinedPIDThreshold;   
  Bool_t fUseCascadePID;            
  AliAODPidHF *fPidObjCascPi;         
  AliAODPidHF *fPidObjCascPr;         
  AliAODPidHF *fPidObjCascKa;         
//  Bool_t   fUseOnTheFlyV0;          /// Flag to check if we use on-the-fly v0
  
  Int_t fProdTrackTPCNclsPIDMin;      
  Double_t fProdTrackTPCNclsRatioMin;      
  Bool_t   fProdUseAODFilterBit;    
  Double_t fProdMassTolLambda;      
  Double_t fProdMassTolOmega;          
  Double_t fProdMassRejXi;          
  Double_t fProdRfidMinV0;          
  Double_t fProdRfidMaxV0;          
  Double_t fProdRfidMinOmega;          
  Double_t fProdRfidMaxOmega;          
  Double_t fProdCascProperDecayLengthMax;        
  Double_t fProdDcaOmegaDaughtersMax;  
  Double_t fProdDcaV0DaughtersMax;  
  Double_t fProdDcaBachToPrimVertexMin;  
  Double_t fProdDcaV0ToPrimVertexMin;  
  Double_t fProdDcaV0PrToPrimVertexMin;  
  Double_t fProdDcaV0PiToPrimVertexMin;  
  Double_t fProdXiCosineOfPoiningAngleMin;  
  Double_t fProdV0CosineOfPoiningAngleXiMin;  
  Double_t fProdCascNTPCClustersMin;         
  Double_t fProdCascEtaMin; 
  Double_t fProdCascEtaMax; 
  Double_t fProdCascRapMin; 
  Double_t fProdCascRapMax; 
  Double_t fProdRoughMassTol;       
  Double_t fProdRoughPtMin;         

  Bool_t fExcludePionTPC; 
  Bool_t fExcludeProtonTPC; 
  Bool_t fExcludeKaonTPC; 
  Double_t fExcludenSigmaPionTPC; 
  Double_t fExcludenSigmaProtonTPC; 
  Double_t fExcludenSigmaKaonTPC; 
  Double_t fSigmaElectronTPCMin; 
  Double_t fSigmaElectronTPCMax; 
  Double_t fSigmaElectronTOFMin; 
  Double_t fSigmaElectronTOFMax; 
  
  ClassDef(AliRDHFCutsOmegactoeleOmegafromAODtracks,2);
};

#endif