PWG2/UNICOR/AliDAnalGlobal.cxx

   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 //=============================================================================