AliEmcalTrackingQATask.h

Go to the documentation of this file.

#ifndef ALIEMCALTRACKINGQATASK_H
#define ALIEMCALTRACKINGQATASK_H

#include "AliAnalysisTaskEmcal.h"

class AliParticleContainer;
class THnSparse;
class TH3;

class AliEmcalTrackingQATask : public AliAnalysisTaskEmcal {

 public:
  AliEmcalTrackingQATask();
  AliEmcalTrackingQATask(const char *name); 
  virtual ~AliEmcalTrackingQATask();

  void                   UserCreateOutputObjects();
  void                   SetGeneratorLevelName(const char* name);
  void                   SetDetectorLevelName(const char* name);
  void                   SetDoSigma1OverPt(Bool_t s)      {fDoSigma1OverPt     = s; }
  void                   SetDoSigmaPtOverPtGen(Bool_t s)  {fDoSigmaPtOverPtGen = s; }

 protected:
  Bool_t                 FillHistograms()                               ;
  void                   ExecOnce()                                     ;
  void                   GenerateHistoBins()                            ;
  void                   AllocateDetectorLevelTHnSparse()               ;
  void                   AllocateGeneratorLevelTHnSparse()              ;
  void                   AllocateMatchedParticlesTHnSparse()            ;
  void                   FillDetectorLevelTHnSparse(Double_t cent, Double_t trackEta, Double_t trackPhi, Double_t trackPt, 
                                                    Double_t sigma1OverPt, Int_t mcGen, Byte_t trackType);
  void                   FillGeneratorLevelTHnSparse(Double_t cent, Double_t partEta, Double_t partPhi, Double_t partPt, Int_t mcGen, Byte_t findable);
  void                   FillMatchedParticlesTHnSparse(Double_t cent, Double_t partEta, Double_t partPhi, Double_t partPt,
                   Double_t trackEta, Double_t trackPhi, Double_t trackPt, Byte_t trackType);

  // Task configuration
  Bool_t                fDoSigma1OverPt        ; //  add sigma(1/pt), if false add sigma(pt)/pt instead
  Bool_t                fDoSigmaPtOverPtGen    ; //  MC: if true do sigma((ptgen - ptdet) / ptgen), otherwise do sigma((ptgen - ptdet) / ptdet)

  // Service fields (non-streamed)
  AliMCParticleContainer* fGeneratorLevel      ; 
  AliTrackContainer*    fDetectorLevel         ; 
  Int_t                 fNPtHistBins           ; 
  Double_t*             fPtHistBins            ; 
  Int_t                 fNEtaHistBins          ; 
  Double_t*             fEtaHistBins           ; 
  Int_t                 fNPhiHistBins          ; 
  Double_t*             fPhiHistBins           ; 
  Int_t                 fNCentHistBins         ; 
  Double_t*             fCentHistBins          ; 
  Int_t                 fNPtRelDiffHistBins    ; 
  Double_t*             fPtRelDiffHistBins     ; 
  Int_t                 fNPtResHistBins        ; 
  Double_t*             fPtResHistBins         ; 
  Double_t*             f1OverPtResHistBins    ; 
  Int_t                 fN1OverPtResHistBins   ; 
  Int_t                 fNIntegerHistBins      ; 
  Double_t*             fIntegerHistBins       ; 


  // Histograms
  THnSparse*            fTracks                ; 
  THnSparse*            fParticlesPhysPrim     ; 
  THnSparse*            fParticlesMatched      ; 
  
 private:
  AliEmcalTrackingQATask(const AliEmcalTrackingQATask&);            // not implemented
  AliEmcalTrackingQATask &operator=(const AliEmcalTrackingQATask&); // not implemented

  ClassDef(AliEmcalTrackingQATask, 3) // Track QA task (efficiency and pt resolution)
};
#endif