1 /*************************************************************************
2 * Copyright(c) 1998-2048, 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 **************************************************************************/
16 // Author: Dariusz Miskowiec <mailto:d.miskowiec@gsi.de> 2007
18 //=============================================================================
19 // event global variable analyzer
20 // Histogram the event variables like multiplicity, flow vector, vertex z,
22 //=============================================================================
27 #include "AliUnicorEvent.h"
28 #include "AliUnicorAnalGlobal.h"
31 ClassImp(AliUnicorAnalGlobal)
33 //=============================================================================
34 AliUnicorAnalGlobal::AliUnicorAnalGlobal(const char *nam) : AliUnicorAnal(nam)
38 TH1D *mult = new TH1D("mult","mult",5000,-0.5,4999.5);
39 mult->SetXTitle("multiplicity");
40 TH1D *cent = new TH1D("cent","cent",100,0,1);
41 cent->SetXTitle("centrality");
42 TH2D *cemu = new TH2D("cemu","cemu",100,0,1,5000,-0.5,4999.5);
43 cemu->SetXTitle("centrality");
44 cemu->SetYTitle("Nch");
45 TH2D *dire = new TH2D("dire","dire",160,-40,40,160,-40,40);
46 dire->SetXTitle("Qx (GeV)");
47 dire->SetYTitle("Qy (GeV)");
48 TH1D *zver = new TH1D("zver","zver",120,-1.2,1.2);
49 zver->SetXTitle("normalized z-vertex");
57 //=============================================================================
58 void AliUnicorAnalGlobal::Process(AliUnicorEvent *ev) const
60 // fill event variable histograms
62 TH1D *mult = (TH1D*) fHistos.At(0);
63 TH1D *cent = (TH1D*) fHistos.At(1);
64 TH2D *cemu = (TH2D*) fHistos.At(2);
65 TH2D *dire = (TH2D*) fHistos.At(3);
66 TH1D *zver = (TH1D*) fHistos.At(4);
68 double n = ev->NGoodParticles();
70 cent->Fill(ev->Centrality(),1.0);
71 cemu->Fill(ev->Centrality(),n);
74 dire->Fill(qx,qy,1.0);
75 zver->Fill(ev->Zver(),1.0);
77 //=============================================================================