AliAnalysisTaskSEHFQA.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKSEHFQUALITYASSURANCE_H
#define ALIANALYSISTASKSEHFQUALITYASSURANCE_H

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

/* $Id$ */ 

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

#include <TROOT.h>
#include <TSystem.h>
#include "TH1F.h"
#include "TH2F.h"
#include "TH3F.h"
#include "TH1I.h"
#include "TProfile2D.h"

#include "AliAnalysisTaskSE.h"

class AliRDHFCuts;
class TH1F;
class AliAODEvent;
class AliFlowEvent;
class AliFlowTrackCuts;

class AliAnalysisTaskSEHFQA : public AliAnalysisTaskSE
{

 public:

  enum DecChannel {kDplustoKpipi,kD0toKpi,kDstartoKpipi,kDstoKKpi,kD0toKpipipi,kLambdactopKpi,kLambdactoV0};

  AliAnalysisTaskSEHFQA();
  AliAnalysisTaskSEHFQA(const char *name, DecChannel ch, AliRDHFCuts* cuts);
  virtual ~AliAnalysisTaskSEHFQA();

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

  void SetReadMC(Bool_t mcflag){fReadMC=mcflag;}
  void SetSimpleMode(Bool_t flag){fSimpleMode=flag;}
  void SetTrackOn(Bool_t trackon=kTRUE){fOnOff[0]=trackon;}
  void SetPIDOn(Bool_t pidon=kTRUE){fOnOff[1]=pidon;}
  void SetCentralityOn(Bool_t centron=kTRUE){fOnOff[2]=centron;}
  void SetEvSelectionOn(Bool_t evselon=kTRUE){fOnOff[3]=evselon;}
  void SetFlowObsOn(Bool_t flowobson=kTRUE){fOnOff[4]=flowobson;}
  void SetUseSelectionBit(Bool_t selectionbiton=kTRUE){fUseSelectionBit=selectionbiton;}
  void SetSecondCentralityEstimator(AliRDHFCuts::ECentrality est){fEstimator = est;}
  void SetFillDistributionsForTrackEffChecks(Bool_t filldistrtrackeffcheckson=kFALSE){fFillDistrTrackEffChecks = filldistrtrackeffcheckson;}
  void SetAODMismatchProtection(Int_t opt=1) {fAODProtection=opt;}

  AliRDHFCuts* GetCutObject() const {return fCuts;}
  DecChannel GetDecayChannel()const {return fDecayChannel;}
  Bool_t GetTrackStatus() const {return fOnOff[0];}
  Bool_t GetPIDStatus() const {return fOnOff[1];}
  Bool_t GetCentralityStatus() const {return fOnOff[2];}
  Bool_t GetEvSelStatus() const {return fOnOff[3];}
  Bool_t GetFlowObsStatus() const {return fOnOff[4];}
  Bool_t GetUseSelectionBit() const {return fUseSelectionBit;}
  AliRDHFCuts::ECentrality GetSecondCentralityEstimator()const {return fEstimator;}
  Bool_t GetFillDistributionsForTrackEffChecks()const {return fFillDistrTrackEffChecks;}

 private:
  AliAnalysisTaskSEHFQA(const AliAnalysisTaskSEHFQA &source);
  AliAnalysisTaskSEHFQA operator=(const AliAnalysisTaskSEHFQA &source);
  void FillFlowObs(AliAODEvent *aod);

 TList* fOutputEntries;    
 TList* fOutputPID;        
 TList* fOutputTrack;      
 TList* fOutputCounters;   
 TList* fOutputCheckCentrality;   
 TList* fOutputEvSelection; 
 TList* fOutputFlowObs;    
 DecChannel fDecayChannel; //identify the decay channel
 AliRDHFCuts* fCuts;       // object containing cuts 
 AliFlowEvent *fFlowEvent; 
 AliFlowTrackCuts *fRFPcuts; 
 AliRDHFCuts::ECentrality fEstimator; //2nd estimator for centrality
 Bool_t fReadMC;           
 Bool_t fSimpleMode;       
 Bool_t fUseSelectionBit;  
 Bool_t fOnOff[5];         
 Bool_t fFillDistrTrackEffChecks;
 Int_t fAODProtection;     

 //___ Declaration of histograms
 TH1F* fHisNentries;                         
 TH2F* fHisNentriesSelBit;                   
 TH1F* fHisTOFflags;                         
 TH1F* fHisTOFsig;                           
 TH1F* fHisTOFstartTimeMask;                 
 TH1F* fHisTOFstartTimeRes;                  
 TH1F* fHisTOFstartTimeDistrib;              
 TH1F* fHisTOFtime;                          
 TH2F* fHisTOFtimeKaonHyptime;               
 TH2F* fHisTOFtimeKaonHyptimeAC;             
 TH2F* fHisTOFsigmaKSigPid;                  
 TH2F* fHisTOFsigmaPionSigPid;               
 TH2F* fHisTOFsigmaProtonSigPid;             
 TH1F* fHisTOFsigPid3sigPion;                
 TH1F* fHisTOFsigPid3sigKaon;                
 TH1F* fHisTOFsigPid3sigProton;              
 TH1F* fHisTPCsig;                           
 TH2F* fHisTPCsigvsp;                        
 TH2F* fHisTPCsigvspAC;                      
 TH2F* fHisTPCsigmaK;                        
 TH2F* fHisTPCsigmaPion;                     
 TH2F* fHisTPCsigmaProton;                   
 TH2F* fHisTPCsigNvsPtAllTracks;             
 TH2F* fHisTPCsigNvsPhiAllTracks;            
 TH2F* fHisTPCsigNvsEtaAllTracks;            
 TH2F* fHisTPCsigNvsPtDaughters;             
 TH2F* fHisTPCsigNvsPhiDaughters;            
 TH2F* fHisTPCsigNvsEtaDaughters;            
 TH2F* fHisTOFsigmaMCKSigPid;                
 TH2F* fHisTOFsigmaMCPionSigPid;             
 TH2F* fHisTOFsigmaMCProtonSigPid;           
 TH2F* fHisTPCsigmaMCK;                      
 TH2F* fHisTPCsigmaMCPion;                   
 TH2F* fHisTPCsigmaMCProton;                 
 TH1F* fHisnClsITS;                          
 TH1F* fHisnClsITSselTr;                     
 TH1F* fHisnClsITSSA;                        
 TH1F* fHisnClsITSSAspdAny;                  
 TH1F* fHisnClsITSSAspdIn;                   
 TH1F* fHisnClsITSSAspdOut;                  
 TH1F* fHisnLayerITS;                        
 TH1F* fHisnLayerITSselTr;                   
 TH1F* fHisnLayerITSsa;                      
 TH1F* fHisnClsSPD;                          
 TH1F* fHisptGoodTr;                         
 TH1F* fHisptGoodTrFromDaugh;                
 TH1F* fHisptGoodTrFromDaugh_filt;           
 TH1F* fHisdistrGoodTr;                      
 TH1F* fHisdistrSelTr;                       
 TH1F* fHisd0dau;                            
 TH1F* fHisd0dau_filt;                       
 TH2F* fHisd0dauphi;                         
 TH2F* fHisd0dauphi_filt;                    
 TH1F* fHisd0zdau;                           
 TH1F* fHisd0zdau_filt;                      
 TH2F* fHisd0zdauphi;                        
 TH2F* fHisd0zdauphi_filt;                   
 TH1F* fHisd0TracksSPDin;                    
 TH1F* fHisd0TracksSPDany;                   
 TH1F* fHisd0TracksFilterBit4;               
 TH1F* fHisd0TracksTPCITSSPDany;             
 TH1F* fHisPtDaughters;                      
 TH1F* fHisPhiDaughters;                     
 TH1F* fHisEtaDaughters;                     
 TH2F* fHisEtavsPhiDaughters;                
 TH2F* fHisNTPCclsvsPtDaughters;             
 TH2F* fHisNTPCclsvsPhiDaughters;            
 TH2F* fHisNTPCclsvsEtaDaughters;            
 TH2F* fHisNTPCCrossedRowsvsPtDaughters;     
 TH2F* fHisNTPCCrossedRowsvsPhiDaughters;    
 TH2F* fHisNTPCCrossedRowsvsEtaDaughters;    
 TH2F* fHisRatioCRowsOverFclsvsPtDaughters;  
 TH2F* fHisRatioCRowsOverFclsvsPhiDaughters; 
 TH2F* fHisRatioCRowsOverFclsvsEtaDaughters; 
 TH2F* fHisNITSclsvsPtDaughters;             
 TH2F* fHisNITSclsvsPhiDaughters;            
 TH2F* fHisNITSclsvsEtaDaughters;            
 TH1F* fHisSPDclsDaughters;                  
 TH1F* fHisPtAllTracks;                      
 TH1F* fHisPhiAllTracks;                     
 TH1F* fHisEtaAllTracks;                     
 TH2F* fHisEtavsPhiAllTracks;                
 TH2F* fHisNTPCclsvsPtAllTracks;             
 TH2F* fHisNTPCclsvsPhiAllTracks;            
 TH2F* fHisNTPCclsvsEtaAllTracks;            
 TH2F* fHisNTPCCrossedRowsvsPtAllTracks;     
 TH2F* fHisNTPCCrossedRowsvsPhiAllTracks;    
 TH2F* fHisNTPCCrossedRowsvsEtaAllTracks;    
 TH2F* fHisRatioCRowsOverFclsvsPtAllTracks;  
 TH2F* fHisRatioCRowsOverFclsvsPhiAllTracks; 
 TH2F* fHisRatioCRowsOverFclsvsEtaAllTracks; 
 TH2F* fHisNITSclsvsPtAllTracks;             
 TH2F* fHisNITSclsvsPhiAllTracks;            
 TH2F* fHisNITSclsvsEtaAllTracks;            
 TH1F* fHisSPDclsAllTracks;                  
 TH1F* fHisdistrFakeTr;                      
 TH1F* fHisd0f;                              
 TH1F* fHisd0f_filt;                         
 TH1F* fHisptFakeTr;                         
 TH1F* fHisptFakeTrFromDaugh;                
 TH1F* fHisptFakeTrFromDaughFilt;            
 TH1F* fHisNtracklets;                       
 TH1F* fHisNtracklets01;                     
 TH1F* fHisNtracklets01AllEv;                
 TH1F* fHisMult;                             
 TH1F* fHisMultFBit4;                        
 TH1F* fHisMultComb05;                       
 TH1F* fHisMultComb08;                       
 TH1F* fHisNtrackletsIn;                     
 TH1F* fHisMultIn;                           
 TH1F* fHisNtrackletsOut;                    
 TH1F* fHisMultOut;                          
 TH2F* fHisMultvsPercentile;                 
 TH2F* fHisntrklvsPercentile;                
 TH2F* fHisntrklvsPercentile01;              
 TH2F* fHisntrklvsPercentile01AllEv;         
 TH2F* fHisnTPCTracksvsPercentile;           
 TH2F* fHisnTPCITSTracksvsPercentile;        
 TH2F* fHisnTPCITS1SPDTracksvsPercentile;    
 TH2F* fHisStdEstimSignalPercentile;         
 TH2F* fHisStdEstimSignalNtrackletsIn;       
 TH1F* fHisStdEstimSignal;                   
 TH2F* fHisStdPercentileSecondPercentile;    
 TH2F* fHisStdSignalSecondSignal;            
 TH2F* fHisStdPercentileOldFrwPercentile;    
 TH1F* fHisStdPercentileOldFrwPercentileDev; 
 TH1F* fHisxvtx;                             
 TH1F* fHisyvtx;                             
 TH1F* fHiszvtx;                             
 TH1F* fHisxvtxSelEv;                        
 TH1F* fHisyvtxSelEv;                        
 TH1F* fHiszvtxSelEv;                        
 TH1F* fHisWhichVert;                        
 TH1F* fHisWhichVertSelEv;                   
 TH2F* fHisTrigCent;                         
 TH2F* fHisTrigMul;                          
 TH2F* fHisTrigCentSel;                      
 TH2F* fHisTrigMulSel;                       
 TH1F* fHisWhyEvRejected;                    
 TH2F* fHisnClsITSvsNtrackletsSel;           
 TH2F* fHiszvtxvsSPDzvtx;                    
 TH2F* fHiszvtxvsSPDzvtxSel;                 
 TH2F* fHisFEvents;                          
 TH3F* fHisTPCVZE_AngleQ;                    
 TH2F* fHisCentVsMultRPS;                    
 TH2F* fHisAngleQ[3];                        
 TH3F* fHisPhiEta[3];                        
 TProfile2D *fHisQ[3];                       

 ClassDef(AliAnalysisTaskSEHFQA,17); 
};

#endif