AliAnalysisTaskSESelectHF4Prong.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKSESELECTHF4PRONG_H
#define ALIANALYSISTASKSESELECTHF4PRONG_H

/* Copyright(c) 1998-2008, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

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

#include <TROOT.h>
#include <TSystem.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TList.h>
#include <TClonesArray.h>
#include <TChain.h>

#include "AliAnalysisTaskSE.h"
#include "AliRDHFCutsD0toKpipipi.h"
#include "AliAODMCParticle.h"
class AliAODEvent;
class AliAnalysisManager;
class AliAnalysisVertexingHF;
class AliRDHFCuts;

class AliAnalysisTaskSESelectHF4Prong : public AliAnalysisTaskSE
{
 public:

  AliAnalysisTaskSESelectHF4Prong();
  AliAnalysisTaskSESelectHF4Prong(const char *name,AliRDHFCutsD0toKpipipi* cuts);
  virtual ~AliAnalysisTaskSESelectHF4Prong();

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

  void SetPtBinH(Double_t* ptlimits);
  void PrintPtBinHandMCFlag();
  void SetMCTruth(int flag) {fMCTruth = flag;};
  Int_t GetMCTruth() const {return fMCTruth;}
  
  void AnalysisReflection(AliAODEvent* aodIn, AliAODRecoDecayHF4Prong* d, Int_t hyp01, Int_t hyp02, Int_t hyp03, Int_t hyp04);
  Int_t StudyMCTruth(TClonesArray* mcArray, AliAODRecoDecayHF4Prong* d);
  void FillReflHistos(AliAODRecoDecayHF4Prong* d, Int_t hyp01, Int_t hyp02, Int_t hyp03, Int_t hyp04, Int_t flagLayer1, Int_t flagLayer2, Int_t flagLayer3, Int_t flagLayer4);
  void FillReflD0Histos(AliAODRecoDecayHF4Prong* d, Int_t hyp01, Int_t hyp02, Int_t hyp03, Int_t hyp04, Int_t flagLayforD0_1, Int_t flagLayforD0_2, Int_t flagLayforD0_3, Int_t flagLayforD0_4);
  void FillReflD0barHistos(AliAODRecoDecayHF4Prong* d, Int_t hyp01, Int_t hyp02, Int_t hyp03, Int_t hyp04, Int_t flagLayforD0bar_1, Int_t flagLayforD0bar_2, Int_t flagLayforD0bar_3, Int_t flagLayforD0bar_4);
 
 private:

  AliAnalysisTaskSESelectHF4Prong(const AliAnalysisTaskSESelectHF4Prong &source);
  AliAnalysisTaskSESelectHF4Prong& operator=(const AliAnalysisTaskSESelectHF4Prong& source); 
  TClonesArray *fVerticesHFTClArr;     
  TClonesArray *fCharm4ProngTClArr;        

  Double_t fPtBinH[6];            
  Double_t fmassD0[2];            
  Double_t fmassD0bar[2];         
  Int_t fSelected;            
  Int_t fMCTruth;           

  TList *fOutput;                             
  TList *fOutput2;                            
  TList *fOutput3;                            
  TList *fOutput4;                            
  TList *fOutput5;                            
  TList *fOutputC;                            

  //output histograms
  TH1F *fhInvMassD0Sum10MevBin1;                        
  TH1F *fhInvMassD0barSum10MevBin1;                     
  TH1F *fhInvMassSumAll10MevBin1;                       
  TH1F *fhInvMassD0Sum5MevBin1;                         
  TH1F *fhInvMassD0barSum5MevBin1;                      
  TH1F *fhInvMassSumAll5MevBin1;                        

  TH1F *fhInvMassD0Sum10MevBin2;                        
  TH1F *fhInvMassD0barSum10MevBin2;                     
  TH1F *fhInvMassSumAll10MevBin2;                       
  TH1F *fhInvMassD0Sum5MevBin2;                         
  TH1F *fhInvMassD0barSum5MevBin2;                      
  TH1F *fhInvMassSumAll5MevBin2;                        

  TH1F *fhInvMassD0Sum10MevBin3;                        
  TH1F *fhInvMassD0barSum10MevBin3;                     
  TH1F *fhInvMassSumAll10MevBin3;                       
  TH1F *fhInvMassD0Sum5MevBin3;                         
  TH1F *fhInvMassD0barSum5MevBin3;                      
  TH1F *fhInvMassSumAll5MevBin3;                        

  TH1F *fhInvMassD0Sum10MevBin4;                        
  TH1F *fhInvMassD0barSum10MevBin4;                     
  TH1F *fhInvMassSumAll10MevBin4;                       
  TH1F *fhInvMassD0Sum5MevBin4;                         
  TH1F *fhInvMassD0barSum5MevBin4;                      
  TH1F *fhInvMassSumAll5MevBin4;                        

  TH1F *fhInvMassD0Sum10MevBin5;                        
  TH1F *fhInvMassD0barSum10MevBin5;                     
  TH1F *fhInvMassSumAll10MevBin5;                       
  TH1F *fhInvMassD0Sum5MevBin5;                         
  TH1F *fhInvMassD0barSum5MevBin5;                      
  TH1F *fhInvMassSumAll5MevBin5;                        

  TH2F *fhReflBin1;         
  TH2F *fhReflBin2;           
  TH2F *fhReflBin3;           
  TH2F *fhReflBin4;           
  TH2F *fhReflBin5;           
  TH2F *fhReflD0Bin1;         
  TH2F *fhReflD0Bin2;           
  TH2F *fhReflD0Bin3;           
  TH2F *fhReflD0Bin4;           
  TH2F *fhReflD0Bin5;           
  TH2F *fhReflD0barBin1;        
  TH2F *fhReflD0barBin2;        
  TH2F *fhReflD0barBin3;        
  TH2F *fhReflD0barBin4;        
  TH2F *fhReflD0barBin5;        

  TH2F *fScatterP4PID;            
  TH2F *fPtVsY;             
  TH2F *fPtVsYAll;            

  TH1F *fEventCounter;               
  TH1F *fCutDCA;               
  TH1F *fCutDCA3;              
  TH1F *fCutDCA2;              
  TH1F *fCutDCA5;              
  TH1F *fCutVertexDist2;             
  TH1F *fCutVertexDist3;             
  TH1F *fCutVertexDist4;             
  TH1F *fCutCosinePoint;             
  TH1F *fCutPt;                      
  TH1F *fCutY;                       
  TH1F *fPIDSel;               
  TH1F *fPIDSelBin1;                     
  TH1F *fPIDSelBin2;                     
  TH1F *fPIDSelBin3;                     
  TH1F *fPIDSelBin4;                     
  TH1F *fPIDSelBin5;                     
  TH2F *fMultipleHyps;               
  TH2F *fMultipleHypsType;             

  TH1F *fPtSel;                      
 
  AliRDHFCutsD0toKpipipi *fCuts;      

  ClassDef(AliAnalysisTaskSESelectHF4Prong,2); 
};

#endif