AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKSEXICPLUS2XIPIPIFROMAODTRACKS_H
#define ALIANALYSISTASKSEXICPLUS2XIPIPIFROMAODTRACKS_H

/**************************************************************************
 * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */ 

#include "TROOT.h"
#include "TSystem.h"

#include "AliAnalysisTaskSE.h"
#include "AliAODEvent.h"
#include "AliPID.h"
#include "AliRDHFCutsXicPlustoXiPiPifromAODtracks.h"


class THnSparse;
class TH1F;
class TH2F;
class TClonesArray;
class AliAODRecoCascadeHF3Prong;
class AliAODPidHF;
class AliESDtrackCuts;
class AliESDVertex;
class AliAODMCParticle;

class AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks : public AliAnalysisTaskSE 
{
 public:
  AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks();
  AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks(const Char_t* name, AliRDHFCutsXicPlustoXiPiPifromAODtracks* cuts, Bool_t writeVariableTree=kTRUE);
  virtual ~AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks();

  virtual void UserCreateOutputObjects();
  virtual void Init();
  virtual void LocalInit() {Init();}
  virtual void UserExec(Option_t *option);
  virtual void Terminate(Option_t *option);

  void FillROOTObjects(AliAODRecoCascadeHF3Prong *xicobj, AliAODMCParticle *mcpart, AliAODMCParticle *mcdau1, AliAODMCParticle *mcdau2, AliAODMCParticle *mcdauxi, Int_t mcnused);
  void MakeAnalysis(AliAODEvent *aod, TClonesArray *mcArray);

  
  void SetMC(Bool_t theMCon) {fUseMCInfo = theMCon;}
  Bool_t GetMC() const {return fUseMCInfo;}
  
  void   SetReconstructPrimVert(Bool_t a) { fReconstructPrimVert=a; }
  void SelectCascade( const AliVEvent *event,Int_t nCascades,Int_t &nSeleCasc, Bool_t *seleCascFlags);
  void SelectTrack( const AliVEvent *event, Int_t trkEntries, Int_t &nSeleTrks,Bool_t *seleFlags);
  Bool_t SelectLikeSign(AliAODTrack *trk1, AliAODTrack *trk2);
  AliAODRecoCascadeHF3Prong* MakeCascadeHF3Prong(AliAODcascade *casc, AliAODTrack *trk1, AliAODTrack *trk2, AliAODEvent *aod, AliAODVertex *secvert, Double_t dispersion);
  
 private:
  
  AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks(const AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks &source);
  AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks& operator=(const AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks& source); 
  
  void DefineTreeVariables();
  void DefineGeneralHistograms();
  void DefineAnalysisHistograms();

  AliAODVertex *CallPrimaryVertex(AliAODcascade *casc, AliAODTrack *trk1, AliAODTrack *trk2, AliAODEvent *evt);
  AliAODVertex* PrimaryVertex(const TObjArray *trkArray,AliVEvent *event);
  AliAODVertex* CallReconstructSecondaryVertex(AliAODTrack *trk1, AliAODTrack *trk2,Double_t &disp);
  AliAODVertex* ReconstructSecondaryVertex(TObjArray *trkArray, Double_t &dispersion,Bool_t useTRefArray=kTRUE);
  
  Bool_t fUseMCInfo;          
  TList *fOutput;             
  TList *fOutputAll;          
  TList *fListCuts;           
  TH1F *fCEvents;             
  TH1F *fHTrigger;            
  TH1F *fHCentrality;         
  AliRDHFCutsXicPlustoXiPiPifromAODtracks *fAnalCuts;      
  Bool_t fIsEventSelected;           
  Bool_t    fWriteVariableTree;       
  TTree    *fVariablesTree;           
  Bool_t fReconstructPrimVert;            
  Bool_t fIsMB;            
  Bool_t fIsSemi;          
  Bool_t fIsCent;          
  Bool_t fIsINT7;          
  Bool_t fIsEMC7;          
  Float_t *fCandidateVariables;     
  AliAODVertex *fVtx1;              
  AliESDVertex *fV1;                
  Double_t fBzkG;                   
  Float_t  fCentrality;             
  Float_t  fTriggerCheck;           
  
  //--------------------- My histograms ------------------
  THnSparse*  fHistoXicMass;        
  
  TH1F*  fHistoDcaPi1Pi2;                    
  TH1F*  fHistoDcaPiCasc;                    
  TH1F*  fHistoLikeDecayLength;              
  TH1F*  fHistoLikeDecayLengthXY;            
  TH1F*  fHistoXicCosPAXY;                   
  
  TH1F*  fHistoXiMass;                       
  TH1F*  fHistoCascDcaXiDaughters;           
  TH1F*  fHistoCascDcaV0Daughters;           
  TH1F*  fHistoCascDcaV0ToPrimVertex;        
  TH1F*  fHistoCascDcaPosToPrimVertex;       
  TH1F*  fHistoCascDcaNegToPrimVertex;       
  TH1F*  fHistoCascDcaBachToPrimVertex;      
  TH1F*  fHistoCascCosPAXiPrim;              
  TH1F*  fHistoXiPt;                         
  
  TH1F*  fHistoPiPt;                         
  TH1F*  fHistoPid0;                         
  TH1F*  fHistonSigmaTPCpi;                  
  TH1F*  fHistonSigmaTOFpi;                  
  TH1F*  fHistoProbPion;                     

  TH2F*  fHistoXiMassvsPtRef;                    
  TH2F*  fHistoXiMassvsPtRef2;                    
  TH2F*  fHistoXiMassvsPtRef3;                    
  TH2F*  fHistoXiMassvsPtRef4;                    
  TH2F*  fHistoXiMassvsPtRef5;                    
  TH2F*  fHistoXiMassvsPtRef6;                    
  TH1F*  fHistoPiPtRef;                      

  ClassDef(AliAnalysisTaskSEXicPlus2XiPiPifromAODtracks,5); 
};
#endif