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