PWG4/AliAnaGammaJetFinder.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15 /* $Id$ */
  16
  17 /* History of cvs commits:
  18  *
  19  * $Log$
  20  * Revision 1.1.2.1  2007/07/26 10:32:09  schutz
  21  * new analysis classes in the the new analysis framework
  22  *
  23  *
  24  */
  25
  26 //_________________________________________________________________________
  27 // Class for the analysis of gamma-jet (standard jet finder) correlations
  28 //*-- Author: Gustavo Conesa (LNF-INFN)
  29 //////////////////////////////////////////////////////////////////////////////
  30
  31
  32 // --- ROOT system ---
  33 #include "Riostream.h"
  34
  35 //---- AliRoot system ----
  36 #include "AliAnaGammaJetFinder.h"
  37 #include "AliLog.h"
  38
  39   ClassImp(AliAnaGammaJetFinder)
  40
  41
  42 //____________________________________________________________________________
  43   AliAnaGammaJetFinder::AliAnaGammaJetFinder() :
  44     AliAnaGammaCorrelation(),
  45     fhDeltaEtaJet(0), fhDeltaPhiJet(0),
  46     fhDeltaPtJet(0), fhPtRatJet(0)
  47 {
  48   //Default Ctor
  49
  50   SetCorrelationType(kJetFinder);
  51   //Initialize parameters
  52   InitParameters();
  53 }
  54
  55 //____________________________________________________________________________
  56 AliAnaGammaJetFinder::AliAnaGammaJetFinder(const AliAnaGammaJetFinder & g) :
  57   AliAnaGammaCorrelation(g),
  58   fhDeltaEtaJet(g.fhDeltaEtaJet), fhDeltaPhiJet(g.fhDeltaPhiJet),
  59   fhDeltaPtJet(g.fhDeltaPtJet), fhPtRatJet(g.fhPtRatJet)
  60 {
  61   // cpy ctor
  62
  63 }
  64
  65 //_________________________________________________________________________
  66 AliAnaGammaJetFinder & AliAnaGammaJetFinder::operator = (const AliAnaGammaJetFinder & source)
  67 {
  68   // assignment operator
  69
  70   if(this == &source)return *this;
  71   ((AliAnaGammaCorrelation *)this)->operator=(source);
  72   fhDeltaEtaJet = source.fhDeltaEtaJet;
  73   fhDeltaPhiJet = source.fhDeltaPhiJet;
  74   fhDeltaPtJet = source.fhDeltaPtJet;
  75   fhPtRatJet = source.fhPtRatJet;
  76
  77   return *this;
  78
  79 }
  80
  81 //____________________________________________________________________________
  82 AliAnaGammaJetFinder::~AliAnaGammaJetFinder()
  83 {
  84
  85   delete fhDeltaEtaJet ;
  86   delete fhDeltaPhiJet ;
  87   delete fhDeltaPtJet ;
  88   delete fhPtRatJet ;
  89
  90 }
  91
  92
  93 //________________________________________________________________________
  94 TList *  AliAnaGammaJetFinder::GetCreateOutputObjects()
  95 {
  96   // Create histograms to be saved in output file and
  97   // store them in fOutputContainer
  98
  99   AliDebug(1,"Init jet histograms");
 100
 101   TList * outputContainer = new TList() ;
 102   outputContainer->SetName("GammaJetHistos") ;
 103
 104   //---kJet---
 105   fhDeltaPhiJet  = new TH2F
 106     ("DeltaPhiJet","#phi_{#gamma} - #phi_{jet} vs p_{T #gamma}",
 107      200,0,120,200,0,6.4);
 108   fhDeltaPhiJet->SetYTitle("#Delta #phi");
 109   fhDeltaPhiJet->SetXTitle("p_{T #gamma} (GeV/c)");
 110   outputContainer->Add(fhDeltaPhiJet);
 111
 112   fhDeltaEtaJet  = new TH2F
 113     ("DeltaEtaJet","#eta_{#gamma} - #eta_{jet} vs p_{T #gamma}",
 114      200,0,120,200,-2,2);
 115   fhDeltaEtaJet->SetYTitle("#Delta #eta");
 116   fhDeltaEtaJet->SetXTitle("p_{T #gamma} (GeV/c)");
 117   outputContainer->Add(fhDeltaEtaJet);
 118
 119   fhDeltaPtJet  = new TH2F
 120     ("DeltaPtJet","#p_{T #gamma} - #p_{T jet} vs p_{T #gamma}",
 121      200,0,120,100,-10,10);
 122   fhDeltaPtJet->SetYTitle("#Delta #p_{T}");
 123   fhDeltaPtJet->SetXTitle("p_{T #gamma} (GeV/c)");
 124   outputContainer->Add(fhDeltaPtJet);
 125
 126   fhPtRatJet  = new TH2F
 127     ("PtRatJet","#p_{T jet} / #p_{T #gamma} vs p_{T #gamma}",
 128      200,0,120,200,0,5);
 129   fhPtRatJet->SetYTitle("ratio");
 130   fhPtRatJet->SetXTitle("p_{T #gamma} (GeV/c)");
 131   outputContainer->Add(fhPtRatJet);
 132
 133   SetOutputContainer(outputContainer);
 134
 135   return outputContainer;
 136 }
 137
 138  //____________________________________________________________________________
 139 void AliAnaGammaJetFinder::InitParameters()
 140 {
 141
 142   //Initialize the parameters of the analysis.
 143
 144   ;
 145
 146 }
 147
 148 //__________________________________________________________________
 149 void AliAnaGammaJetFinder::Print(const Option_t * opt) const
 150 {
 151
 152   //Print some relevant parameters set for the analysis
 153   if(! opt)
 154     return;
 155
 156 }
 157
 158 //__________________________________________________________________
 159 void  AliAnaGammaJetFinder::MakeGammaCorrelation(TParticle * pGamma, TClonesArray *pl, TClonesArray *)
 160 {
 161   //Gamma -Jet  Correlation Analysis
 162   AliDebug(2, "Begin jet analysis");
 163   AliInfo("Not implemented");
 164 }