[u/mrichter/AliRoot.git] / JETAN / AliJetHistos.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

#include <TList.h>
#include <TClonesArray.h>
#include <TH1I.h>
#include <TH1F.h>
#include <TMath.h>

#include "AliAODJet.h"
#include "AliJetHistos.h"

ClassImp(AliJetHistos)

AliJetHistos::AliJetHistos():
  fNJetsH(0x0),
  fPtH(0x0),
  fEtaH(0x0),
  fEneH(0x0),
  fPhiH(0x0)
{
  // Default constructor

  fNJetsH = new TH1I("NJetsH","Number of Jets",12,0,11);
  SetProperties(fNJetsH,"Number of jets","Entries");

  fPtH = new TH1F("PtH","Pt of Jets",50,0.,200.);
  SetProperties(fPtH,"P_{#perp} [GeV]","Entries");

  fEtaH = new TH1F("EtaH","Pseudorapidity of Jets",30,-1.5,1.5);
  SetProperties(fEtaH,"#eta","Entries");

  fEneH = new TH1F("EneH","Energy of Jets",50,0.,200.);
  SetProperties(fEneH,"Energy [GeV]","Entries");

  fPhiH = new TH1F("PhiH","Azimuthal angle of Jets",
                   60,0.,2.0*TMath::Pi());
  SetProperties(fPhiH,"#phi","Entries");
}

AliJetHistos::~AliJetHistos()
{
  delete fNJetsH;
  delete fPtH;
  delete fEtaH;
  delete fEneH;
  delete fPhiH;
}

void AliJetHistos::SetProperties(TH1* h,const char* x, const char* y) const
{
// Sets the histogram style properties
  h->SetMarkerStyle(20);
  h->SetMarkerSize(.5);
  h->SetMarkerColor(2);
  h->SetXTitle(x);
  h->SetYTitle(y);
  h->Sumw2();
}

void AliJetHistos::AddHistosToList(TList *list)
{
  list->Add(fNJetsH);
  list->Add(fPtH);
  list->Add(fEtaH);
  list->Add(fEneH);
  list->Add(fPhiH);
}


void AliJetHistos::FillHistos(TClonesArray *jets)
{
  Int_t nj = jets->GetEntries();
  fNJetsH->Fill(nj,1);

  if (nj == 0 ) return;

  AliAODJet *j;
  for (Int_t i=0;i<nj;i++) {
    j = (AliAODJet *) jets->At(i);
    fPtH->Fill(j->Pt(),1);
    fEtaH->Fill(j->Eta(),1);
    fEneH->Fill(j->E(),1);
    fPhiH->Fill(j->Phi(),1);
  }
  
}
Commit	Line	Data
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
	18	#include <TList.h>
	19	#include <TClonesArray.h>
	20	#include <TH1I.h>
	21	#include <TH1F.h>
	22	#include <TMath.h>
	23
	24	#include "AliAODJet.h"
	25	#include "AliJetHistos.h"
	26
	27	ClassImp(AliJetHistos)
	28
	29	AliJetHistos::AliJetHistos():
	30	fNJetsH(0x0),
	31	fPtH(0x0),
	32	fEtaH(0x0),
	33	fEneH(0x0),
	34	fPhiH(0x0)
	35	{
	36	// Default constructor
	37
	38	fNJetsH = new TH1I("NJetsH","Number of Jets",12,0,11);
	39	SetProperties(fNJetsH,"Number of jets","Entries");
	40
	41	fPtH = new TH1F("PtH","Pt of Jets",50,0.,200.);
	42	SetProperties(fPtH,"P_{#perp} [GeV]","Entries");
	43
	44	fEtaH = new TH1F("EtaH","Pseudorapidity of Jets",30,-1.5,1.5);
	45	SetProperties(fEtaH,"#eta","Entries");
	46
	47	fEneH = new TH1F("EneH","Energy of Jets",50,0.,200.);
	48	SetProperties(fEneH,"Energy [GeV]","Entries");
	49
	50	fPhiH = new TH1F("PhiH","Azimuthal angle of Jets",
	51	60,0.,2.0*TMath::Pi());
	52	SetProperties(fPhiH,"#phi","Entries");
	53	}
	54
	55	AliJetHistos::~AliJetHistos()
	56	{
	57	delete fNJetsH;
	58	delete fPtH;
	59	delete fEtaH;
	60	delete fEneH;
	61	delete fPhiH;
	62	}
	63
	64	void AliJetHistos::SetProperties(TH1* h,const char* x, const char* y) const
65	{
66	// Sets the histogram style properties
67	h->SetMarkerStyle(20);
68	h->SetMarkerSize(.5);
69	h->SetMarkerColor(2);
70	h->SetXTitle(x);
71	h->SetYTitle(y);
72	h->Sumw2();
73	}
74
75	void AliJetHistos::AddHistosToList(TList *list)
76	{
77	list->Add(fNJetsH);
78	list->Add(fPtH);
79	list->Add(fEtaH);
80	list->Add(fEneH);
81	list->Add(fPhiH);
82	}
83
84
85	void AliJetHistos::FillHistos(TClonesArray *jets)
86	{
87	Int_t nj = jets->GetEntries();
88	fNJetsH->Fill(nj,1);
89
90	if (nj == 0 ) return;
91
92	AliAODJet *j;
93	for (Int_t i=0;i<nj;i++) {
94	j = (AliAODJet *) jets->At(i);
95	fPtH->Fill(j->Pt(),1);
96	fEtaH->Fill(j->Eta(),1);
97	fEneH->Fill(j->E(),1);
98	fPhiH->Fill(j->Phi(),1);
99	}
100
101	}