[u/mrichter/AliRoot.git] / PWG2 / UNICOR / AliDAnalGlobal.cxx

// Author: Dariusz Miskowiec <mailto:d.miskowiec@gsi.de> 2007

//=============================================================================
// event global variable analyzer
//=============================================================================

#include <TROOT.h>
#include <TH1.h>
#include <TH2.h>
#include "AliDEvent.h"
#include "AliDAnalGlobal.h"

ClassImp(AliDAnalGlobal)
  
//=============================================================================
AliDAnalGlobal::AliDAnalGlobal(Char_t *nam) : AliDAnal(nam) 
{
  // constructor

  TH1D *mult = new TH1D("mult","mult",5000,-0.5,4999.5);
  mult->SetXTitle("multiplicity");
  TH1D *cent = new TH1D("cent","cent",100,0,1);
  cent->SetXTitle("centrality");
  TH2D *dire = new TH2D("dire","dire",100,-40,40,100,-40,40);
  dire->SetXTitle("Qx (GeV)");
  dire->SetYTitle("Qy (GeV)");
  TH1D *zver = new TH1D("zver","zver",120,-1.2,1.2);
  zver->SetXTitle("normalized z-vertex");
  fHistos.Add(mult);
  fHistos.Add(cent);
  fHistos.Add(dire);
  fHistos.Add(zver);
  gROOT->cd();
  printf("%s object named %s created\n",ClassName(),GetName());
}
//=============================================================================
void AliDAnalGlobal::Process(AliDEvent *ev) 
{
  // fill event variable histograms

  TH1D *mult = (TH1D*) fHistos.At(0);
  TH1D *cent = (TH1D*) fHistos.At(1);
  TH2D *dire = (TH2D*) fHistos.At(2);
  TH1D *zver = (TH1D*) fHistos.At(3);

  mult->Fill(ev->NParticles(),1.0);
  cent->Fill(ev->Centrality(),1.0);
  Double_t qx,qy;
  ev->RP(qx,qy);
  dire->Fill(qx,qy,1.0);
  zver->Fill(ev->Zver(),1.0);
}
//=============================================================================
Commit	Line	Data
7148817a	1	// Author: Dariusz Miskowiec <mailto:d.miskowiec@gsi.de> 2007
	2
	3	//=============================================================================
	4	// event global variable analyzer
	5	//=============================================================================
	6
	7	#include <TROOT.h>
	8	#include <TH1.h>
	9	#include <TH2.h>
	10	#include "AliDEvent.h"
	11	#include "AliDAnalGlobal.h"
	12
	13	ClassImp(AliDAnalGlobal)
	14
	15	//=============================================================================
	16	AliDAnalGlobal::AliDAnalGlobal(Char_t *nam) : AliDAnal(nam)
	17	{
	18	// constructor
	19
	20	TH1D *mult = new TH1D("mult","mult",5000,-0.5,4999.5);
	21	mult->SetXTitle("multiplicity");
	22	TH1D *cent = new TH1D("cent","cent",100,0,1);
	23	cent->SetXTitle("centrality");
	24	TH2D *dire = new TH2D("dire","dire",100,-40,40,100,-40,40);
	25	dire->SetXTitle("Qx (GeV)");
	26	dire->SetYTitle("Qy (GeV)");
	27	TH1D *zver = new TH1D("zver","zver",120,-1.2,1.2);
	28	zver->SetXTitle("normalized z-vertex");
	29	fHistos.Add(mult);
	30	fHistos.Add(cent);
	31	fHistos.Add(dire);
	32	fHistos.Add(zver);
	33	gROOT->cd();
	34	printf("%s object named %s created\n",ClassName(),GetName());
	35	}
	36	//=============================================================================
	37	void AliDAnalGlobal::Process(AliDEvent *ev)
	38	{
	39	// fill event variable histograms
	40
	41	TH1D mult = (TH1D) fHistos.At(0);
	42	TH1D cent = (TH1D) fHistos.At(1);
	43	TH2D dire = (TH2D) fHistos.At(2);
	44	TH1D zver = (TH1D) fHistos.At(3);
	45
	46	mult->Fill(ev->NParticles(),1.0);
	47	cent->Fill(ev->Centrality(),1.0);
	48	Double_t qx,qy;
	49	ev->RP(qx,qy);
	50	dire->Fill(qx,qy,1.0);
	51	zver->Fill(ev->Zver(),1.0);
	52	}
	53	//=============================================================================