LUTs mapping symnames and original global matrices removed from AliGeomManager, which...
[u/mrichter/AliRoot.git] / PWG4 / AliAnaParticleJetFinderCorrelation.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: AliAnaGammaJetFinder.cxx 22232 2007-11-17 16:39:49Z gustavo $ */
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