1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 $ */
17 /* History of cvs commits:
20 * Revision 1.3 2007/10/29 13:48:42 gustavo
21 * Corrected coding violations
23 * Revision 1.2 2007/08/17 12:40:04 schutz
24 * New analysis classes by Gustavo Conesa
26 * Revision 1.1.2.1 2007/07/26 10:32:09 schutz
27 * new analysis classes in the the new analysis framework
32 //_________________________________________________________________________
33 // Class for the analysis of gamma-jet (standard jet finder) correlations
34 //*-- Author: Gustavo Conesa (LNF-INFN)
35 //////////////////////////////////////////////////////////////////////////////
38 // --- ROOT system ---
39 #include "Riostream.h"
41 //---- AliRoot system ----
42 #include "AliAnaGammaJetFinder.h"
45 ClassImp(AliAnaGammaJetFinder)
48 //____________________________________________________________________________
49 AliAnaGammaJetFinder::AliAnaGammaJetFinder() :
50 AliAnaGammaCorrelation(),
51 fhDeltaEtaJet(0), fhDeltaPhiJet(0),
52 fhDeltaPtJet(0), fhPtRatJet(0)
56 SetCorrelationType(kJetFinder);
57 //Initialize parameters
61 //____________________________________________________________________________
62 AliAnaGammaJetFinder::AliAnaGammaJetFinder(const AliAnaGammaJetFinder & g) :
63 AliAnaGammaCorrelation(g),
64 fhDeltaEtaJet(g.fhDeltaEtaJet), fhDeltaPhiJet(g.fhDeltaPhiJet),
65 fhDeltaPtJet(g.fhDeltaPtJet), fhPtRatJet(g.fhPtRatJet)
71 //_________________________________________________________________________
72 AliAnaGammaJetFinder & AliAnaGammaJetFinder::operator = (const AliAnaGammaJetFinder & source)
74 // assignment operator
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;
87 //____________________________________________________________________________
88 AliAnaGammaJetFinder::~AliAnaGammaJetFinder()
90 // Remove all pointers except analysis output pointers.
95 //________________________________________________________________________
96 TList * AliAnaGammaJetFinder::GetCreateOutputObjects()
98 // Create histograms to be saved in output file and
99 // store them in fOutputContainer
101 AliDebug(1,"Init jet histograms");
103 TList * outputContainer = new TList() ;
104 outputContainer->SetName("GammaJetHistos") ;
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);
114 fhDeltaEtaJet = new TH2F
115 ("DeltaEtaJet","#eta_{#gamma} - #eta_{jet} vs p_{T #gamma}",
117 fhDeltaEtaJet->SetYTitle("#Delta #eta");
118 fhDeltaEtaJet->SetXTitle("p_{T #gamma} (GeV/c)");
119 outputContainer->Add(fhDeltaEtaJet);
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);
128 fhPtRatJet = new TH2F
129 ("PtRatJet","#p_{T jet} / #p_{T #gamma} vs p_{T #gamma}",
131 fhPtRatJet->SetYTitle("ratio");
132 fhPtRatJet->SetXTitle("p_{T #gamma} (GeV/c)");
133 outputContainer->Add(fhPtRatJet);
135 SetOutputContainer(outputContainer);
137 return outputContainer;
140 //____________________________________________________________________________
141 void AliAnaGammaJetFinder::InitParameters()
144 //Initialize the parameters of the analysis.
150 //__________________________________________________________________
151 void AliAnaGammaJetFinder::Print(const Option_t * opt) const
154 //Print some relevant parameters set for the analysis
160 //__________________________________________________________________
161 void AliAnaGammaJetFinder::MakeGammaCorrelation(TParticle * pGamma, TClonesArray *pl, TClonesArray *)
163 //Gamma -Jet Correlation Analysis
164 AliDebug(2, "Begin jet analysis");
165 cout<<pGamma<<" "<<pl<<endl;
166 AliInfo("Not implemented");