Extended pt range; corrected pi0/gamma weights

[u/mrichter/AliRoot.git] / PWGGA / PHOSTasks / PHOS_embedding / AliAnalysisTaskPi0DiffEfficiency.h

#ifndef ALIANALYSISTASKPI0DIFFEFFICIENCY_H\r
#define ALIANALYSISTASKPI0DIFFEFFICIENCY_H\r
\r
// Task calculating effciency of pi0 registration \r
// as a difference between spectrum after and before \r
// embedding\r
\r
class TClonesArray;\r
class AliAODCaloCluster ;\r
class AliPHOSAodCluster ;\r
class AliPHOSCalibData ;\r
class AliPHOSGeometry ;\r
\r
#include "AliAnalysisTaskSE.h"\r
\r
class AliAnalysisTaskPi0DiffEfficiency : public AliAnalysisTaskSE {\r
public:\r
  AliAnalysisTaskPi0DiffEfficiency(const char *name = "AliAnalysisTaskPi0DiffEfficiency");\r
  virtual ~AliAnalysisTaskPi0DiffEfficiency() {}\r
  \r
  virtual void   UserCreateOutputObjects();\r
  virtual void   UserExec(Option_t *option);\r
\r
  void SetPHOSBadMap(Int_t mod,TH2I * h)\r
  {\r
    if(fPHOSBadMap[mod]) delete fPHOSBadMap[mod] ;\r
    fPHOSBadMap[mod]=new TH2I(*h) ;\r
    printf("Set %s \n",fPHOSBadMap[mod]->GetName());\r
  }\r
\r
  protected:\r
  AliAnalysisTaskPi0DiffEfficiency(const AliAnalysisTaskPi0DiffEfficiency& a):AliAnalysisTaskSE(a){} // not implemented\r
  AliAnalysisTaskPi0DiffEfficiency& operator=(const AliAnalysisTaskPi0DiffEfficiency& ){return *this;} // not implemented\r
  Bool_t IsGoodChannel(const char * det, Int_t mod,Int_t ix, Int_t iz); //Use addisional bad map for PHOS\r
  void FillHistogram(const char * key,Double_t x) const ; //Fill 1D histogram witn name key\r
  void FillHistogram(const char * key,Double_t x, Double_t y) const ; //Fill 2D histogram witn name key\r
  void FillHistogram(const char * key,Double_t x, Double_t y, Double_t z) const ; //Fill 3D histogram witn name key\r
  void FillHistogram(const char * key,Double_t x, Double_t y, Double_t z, Double_t weight) const ; //Fill 3D histogram witn name key\r
  Bool_t TestLambda(Double_t e,Double_t l1,Double_t l2) ;  //Evaluate Dispersion cuts for photons\r
  Bool_t TestLambda2(Double_t e,Double_t l1,Double_t l2) ;  //Evaluate Dispersion cuts for photons\r
  Double_t TestCPV(Double_t dx, Double_t dz, Double_t pt, Int_t charge);\r
  Bool_t TestTOF(Double_t t, Double_t e) ;\r
  \r
  void ProcessMC() ;\r
  Double_t CoreEnergy(AliPHOSAodCluster * clu); \r
   \r
private:\r
  Bool_t IsSameCluster(AliAODCaloCluster * c1,AliAODCaloCluster * c2)const ;\r
  void EvalLambdas(AliAODCaloCluster * clu, Int_t iR,Double_t &m02, Double_t &m20) ;\r
  Double_t PrimaryWeight(Double_t x) ;\r
 \r
private:\r
  AliStack * fStack ;               //! Stack of primary\r
  THashList * fOutputContainer;     //! final histogram container\r
  TList * fPHOSEvents[1][10][11] ;  //! Containers for events with PHOS photons\r
  TClonesArray * fPHOSEvent ;       //!PHOS photons in current event\r
  TClonesArray * fPHOSEvent1 ;      //!PHOS photons in current event\r
  TClonesArray * fPHOSEvent2 ;      //!PHOS photons in current event\r
  AliPHOSCalibData *fPHOSCalibData; //! PHOS calibration object\r
  TF1 *fNonLinCorr;                 // Non-linearity correction\r
 \r
  //Reaction plane for v2\r
  Float_t fRPfull ; //!Reaction plain calculated with full TPC \r
  Float_t fRPA ;    //!Reaction plain calculated with A-side TPC: eta>0.15 \r
  Float_t fRPC ;    //!Reaction plain calculated with C-side TPC: eta<-0.15\r
  Float_t fRPFar ;  //!Reaction plain calculated with TPC: eta>0.6 \r
  Float_t fRPAFar ; //!Reaction plain calculated with A-side TPC: eta>0.6 \r
  Float_t fRPCFar ; //!Reaction plain calculated with C-side TPC: eta<-0.6\r
\r
  Float_t fCentrality ; //!Centrality of the currecnt event\r
\r
  Int_t fCenBin ;       //! Current centrality bin\r
\r
  TH2I *fPHOSBadMap[6] ;    //!Container for PHOS bad channels map\r
\r
  AliPHOSGeometry  *fPHOSGeo;  //! PHOS geometry\r
  Int_t fEventCounter;         //! number of analyzed events\r
 \r
  ClassDef(AliAnalysisTaskPi0DiffEfficiency, 2); // PHOS analysis task\r
};\r
\r
#endif\r