Updating conv correlation software

[u/mrichter/AliRoot.git] / PWG4 / GammaConv / ConvCorrelations / AliAnaConvCorrBase.h

/* This file is property of and copyright                                 *\r
 * ALICE Experiment at CERN, All rights reserved.                         *\r
 * See cxx source for full Copyright notice                               */\r
\r
/// @file   AliAnaConvCorrPhoton.h\r
/// @author Svein Lindal\r
/// @brief  Base class for analysis of correlations between conversion particles and charged tracks\r
\r
#ifndef ALIANACONVCORRBASE_CXX\r
#define ALIANACONVCORRBASE_CXX\r
\r
#include "Rtypes.h"\r
#include "TMath.h"\r
#include "TList.h"\r
#include "TH1.h"\r
\r
class TH1F;\r
class TH3F;\r
class TH2F;\r
class AliAODConversionParticle;\r
class TClonesArray;\r
class TNtuple;\r
class TString;\r
\r
\r
class AliAnaConvCorrBase : public TObject {\r
\r
public:\r
\r
  \r
\r
  //Constructor / desctructor\r
  AliAnaConvCorrBase(TString name); \r
  virtual ~AliAnaConvCorrBase();\r
  \r
  //Set and get min pt for triggers\r
  // void SetTriggerPt(Float_t pt) { fTriggerPt = pt; }\r
  // inline Float_t GetTriggerPt() const {return fTriggerPt; }\r
\r
\r
  // //Set and get min pt for correlation particles\r
  // void SetCorrelatedPt(Float_t pt) { fCorrelatedPt = pt; }\r
  // inline Float_t GetCorrelatedPt() const {return fCorrelatedPt; }\r
\r
  //CreateHistograms\r
  void CreateBaseHistograms();\r
  //To be overrriden by children. Should always call CreateBaseHistograms()\r
  virtual void CreateHistograms();\r
  \r
  //Get list of histograms\r
  TList * GetHistograms() const { return fHistograms;}\r
\r
  //Add histogram to list\r
  void AddHistogram(TH1 * histogram) { fHistograms->Add(dynamic_cast<TObject*>(histogram));}\r
\r
  //Set and get number of bins in phi direction\r
  void SetNPhiBins(Int_t bins) { fNPhiBins = bins; }\r
  Int_t GetNPhiBins() const { return fNPhiBins;}\r
\r
  ///Get the distance in phi between trigger particle and correlated particle\r
  Float_t GetDPhi(Float_t dPhi) { \r
    if ( dPhi < 3*TMath::PiOver2() && dPhi > - TMath::PiOver2() ) return dPhi;\r
    else return ( (dPhi>0)? dPhi - TMath::TwoPi() : dPhi + TMath::TwoPi() ); \r
  }\r
\r
\r
  TArrayD * GetTriggerBins() const { return fPtBins; }\r
\r
  //Print statistics for histograms\r
  void PrintStatistics();\r
\r
\r
protected:\r
\r
  //Fill histograms\r
  void FillHistograms(Float_t tPt, Float_t cPt, Float_t dPhi, Float_t dEta, Bool_t isolated);\r
  //Fill trigger counter histograms\r
  void FillTriggerCounters(Float_t tPt, Bool_t isolated);\r
\r
 private:\r
  \r
  TString fName; //name of analysis\r
  TList * fHistograms; //List of histograms\r
\r
  Int_t fNPhiBins;  //Nbins in phi direction\r
  TArrayD * fdPhiBins; //!transient phi bin array\r
  TArrayD * fPtBins; //!Array of trigger bins\r
  \r
\r
  TH3F * fHdPhi[2]; //dPhi pt histogram\r
  TH1F * fHNTriggers[2]; //Histograms containing number of triggers in various bins\r
\r
\r
  //Default constructor prohibited\r
  AliAnaConvCorrBase(); //not implemented\r
  AliAnaConvCorrBase(const AliAnaConvCorrBase&); // not implemented\r
  AliAnaConvCorrBase& operator=(const AliAnaConvCorrBase&); // not implemented\r
  ClassDef(AliAnaConvCorrBase, 3); // example of analysis\r
};\r
\r
#endif\r