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.3  2007/10/29 13:48:42  gustavo
  21  * Corrected coding violations
  22  *
  23  * Revision 1.2  2007/08/17 12:40:04  schutz
  24  * New analysis classes by Gustavo Conesa
  25  *
  26  * Revision 1.1.2.1  2007/07/26 10:32:09  schutz
  27  * new analysis classes in the the new analysis framework
  28  *
  29  *
  30  */
  31
  32 //_________________________________________________________________________
  33 // Class for the analysis of gamma-jet (standard jet finder) correlations
  34 //*-- Author: Gustavo Conesa (LNF-INFN)
  35 //////////////////////////////////////////////////////////////////////////////
  36
  37
  38 // --- ROOT system ---
  39 #include "Riostream.h"
  40
  41 //---- AliRoot system ----
  42 #include "AliAnaGammaJetFinder.h"
  43 #include "AliLog.h"
  44
  45   ClassImp(AliAnaGammaJetFinder)
  46
  47
  48 //____________________________________________________________________________
  49   AliAnaGammaJetFinder::AliAnaGammaJetFinder() :
  50     AliAnaGammaCorrelation(),
  51     fhDeltaEtaJet(0), fhDeltaPhiJet(0),
  52     fhDeltaPtJet(0), fhPtRatJet(0)
  53 {
  54   //Default Ctor
  55
  56   SetCorrelationType(kJetFinder);
  57   //Initialize parameters
  58   InitParameters();
  59 }
  60
  61 //____________________________________________________________________________
  62 AliAnaGammaJetFinder::AliAnaGammaJetFinder(const AliAnaGammaJetFinder & g) :
  63   AliAnaGammaCorrelation(g),
  64   fhDeltaEtaJet(g.fhDeltaEtaJet), fhDeltaPhiJet(g.fhDeltaPhiJet),
  65   fhDeltaPtJet(g.fhDeltaPtJet), fhPtRatJet(g.fhPtRatJet)
  66 {
  67   // cpy ctor
  68
  69 }
  70
  71 //_________________________________________________________________________
  72 AliAnaGammaJetFinder & AliAnaGammaJetFinder::operator = (const AliAnaGammaJetFinder & source)
  73 {
  74   // assignment operator
  75
  76   if(this == &source)return *this;
  77   ((AliAnaGammaCorrelation *)this)->operator=(source);
  78   fhDeltaEtaJet = source.fhDeltaEtaJet;
  79   fhDeltaPhiJet = source.fhDeltaPhiJet;
  80   fhDeltaPtJet = source.fhDeltaPtJet;
  81   fhPtRatJet = source.fhPtRatJet;
  82
  83   return *this;
  84
  85 }
  86
  87 //____________________________________________________________________________
  88 AliAnaGammaJetFinder::~AliAnaGammaJetFinder()
  89 {
  90    // Remove all pointers except analysis output pointers.
  91
  92 }
  93
  94
  95 //________________________________________________________________________
  96 TList *  AliAnaGammaJetFinder::GetCreateOutputObjects()
  97 {
  98   // Create histograms to be saved in output file and
  99   // store them in fOutputContainer
 100
 101   AliDebug(1,"Init jet histograms");
 102
 103   TList * outputContainer = new TList() ;
 104   outputContainer->SetName("GammaJetHistos") ;
 105
 106   //---kJet---
 107   fhDeltaPhiJet  = new TH2F
 108     ("DeltaPhiJet","#phi_{#gamma} - #phi_{jet} vs p_{T #gamma}",
 109      200,0,120,200,0,6.4);
 110   fhDeltaPhiJet->SetYTitle("#Delta #phi");
 111   fhDeltaPhiJet->SetXTitle("p_{T #gamma} (GeV/c)");
 112   outputContainer->Add(fhDeltaPhiJet);
 113
 114   fhDeltaEtaJet  = new TH2F
 115     ("DeltaEtaJet","#eta_{#gamma} - #eta_{jet} vs p_{T #gamma}",
 116      200,0,120,200,-2,2);
 117   fhDeltaEtaJet->SetYTitle("#Delta #eta");
 118   fhDeltaEtaJet->SetXTitle("p_{T #gamma} (GeV/c)");
 119   outputContainer->Add(fhDeltaEtaJet);
 120
 121   fhDeltaPtJet  = new TH2F
 122     ("DeltaPtJet","#p_{T #gamma} - #p_{T jet} vs p_{T #gamma}",
 123      200,0,120,100,-10,10);
 124   fhDeltaPtJet->SetYTitle("#Delta #p_{T}");
 125   fhDeltaPtJet->SetXTitle("p_{T #gamma} (GeV/c)");
 126   outputContainer->Add(fhDeltaPtJet);
 127
 128   fhPtRatJet  = new TH2F
 129     ("PtRatJet","#p_{T jet} / #p_{T #gamma} vs p_{T #gamma}",
 130      200,0,120,200,0,5);
 131   fhPtRatJet->SetYTitle("ratio");
 132   fhPtRatJet->SetXTitle("p_{T #gamma} (GeV/c)");
 133   outputContainer->Add(fhPtRatJet);
 134
 135   SetOutputContainer(outputContainer);
 136
 137   return outputContainer;
 138 }
 139
 140  //____________________________________________________________________________
 141 void AliAnaGammaJetFinder::InitParameters()
 142 {
 143
 144   //Initialize the parameters of the analysis.
 145
 146   ;
 147
 148 }
 149
 150 //__________________________________________________________________
 151 void AliAnaGammaJetFinder::Print(const Option_t * opt) const
 152 {
 153
 154   //Print some relevant parameters set for the analysis
 155   if(! opt)
 156     return;
 157
 158 }
 159
 160 //__________________________________________________________________
 161 void  AliAnaGammaJetFinder::MakeGammaCorrelation(TParticle * pGamma, TClonesArray *pl, TClonesArray *)
 162 {
 163   //Gamma -Jet  Correlation Analysis
 164   AliDebug(2, "Begin jet analysis");
 165   cout<<pGamma<<" "<<pl<<endl;
 166   AliInfo("Not implemented");
 167 }