o First Version of TRDnSigma implementation (Xianguo) o still requires some catching...
[u/mrichter/AliRoot.git] / JETAN / AliJetHistos.cxx
CommitLineData
22003a60 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
16/* $Id$ */
17
139cbd96 18//---------------------------------------------------------------------
19// Jet Histos class
20// Creates and fills a few cummon histograms for jet analysis
21//
22//---------------------------------------------------------------------
23
22003a60 24#include <TList.h>
25#include <TClonesArray.h>
26#include <TH1I.h>
27#include <TH1F.h>
28#include <TMath.h>
29
30#include "AliAODJet.h"
31#include "AliJetHistos.h"
32
33ClassImp(AliJetHistos)
34
139cbd96 35///////////////////////////////////////////////////////////////////////
36
22003a60 37AliJetHistos::AliJetHistos():
38 fNJetsH(0x0),
39 fPtH(0x0),
40 fEtaH(0x0),
41 fEneH(0x0),
42 fPhiH(0x0)
43{
44 // Default constructor
139cbd96 45}
46
47//-----------------------------------------------------------------------
48void AliJetHistos::CreateHistos()
49{
50 // create histos
22003a60 51
52 fNJetsH = new TH1I("NJetsH","Number of Jets",12,0,11);
53 SetProperties(fNJetsH,"Number of jets","Entries");
54
55 fPtH = new TH1F("PtH","Pt of Jets",50,0.,200.);
56 SetProperties(fPtH,"P_{#perp} [GeV]","Entries");
57
58 fEtaH = new TH1F("EtaH","Pseudorapidity of Jets",30,-1.5,1.5);
59 SetProperties(fEtaH,"#eta","Entries");
60
61 fEneH = new TH1F("EneH","Energy of Jets",50,0.,200.);
62 SetProperties(fEneH,"Energy [GeV]","Entries");
63
64 fPhiH = new TH1F("PhiH","Azimuthal angle of Jets",
65 60,0.,2.0*TMath::Pi());
66 SetProperties(fPhiH,"#phi","Entries");
67}
68
139cbd96 69//-----------------------------------------------------------------------
22003a60 70AliJetHistos::~AliJetHistos()
71{
139cbd96 72 // Destructor
22003a60 73 delete fNJetsH;
74 delete fPtH;
75 delete fEtaH;
76 delete fEneH;
77 delete fPhiH;
78}
79
139cbd96 80//-----------------------------------------------------------------------
22003a60 81void AliJetHistos::SetProperties(TH1* h,const char* x, const char* y) const
82{
139cbd96 83 // Sets the histogram style properties
22003a60 84 h->SetMarkerStyle(20);
85 h->SetMarkerSize(.5);
86 h->SetMarkerColor(2);
87 h->SetXTitle(x);
88 h->SetYTitle(y);
89 h->Sumw2();
90}
91
139cbd96 92//-----------------------------------------------------------------------
93void AliJetHistos::AddHistosToList(TList *list) const
22003a60 94{
139cbd96 95 // Add histos to the list
22003a60 96 list->Add(fNJetsH);
97 list->Add(fPtH);
98 list->Add(fEtaH);
99 list->Add(fEneH);
100 list->Add(fPhiH);
101}
102
139cbd96 103//-----------------------------------------------------------------------
22003a60 104void AliJetHistos::FillHistos(TClonesArray *jets)
105{
139cbd96 106 // Fill histograms
8bca421f 107 if(!jets)return;
22003a60 108 Int_t nj = jets->GetEntries();
109 fNJetsH->Fill(nj,1);
110
111 if (nj == 0 ) return;
112
113 AliAODJet *j;
114 for (Int_t i=0;i<nj;i++) {
115 j = (AliAODJet *) jets->At(i);
116 fPtH->Fill(j->Pt(),1);
117 fEtaH->Fill(j->Eta(),1);
118 fEneH->Fill(j->E(),1);
119 fPhiH->Fill(j->Phi(),1);
120 }
121
122}