1 // Author: Dariusz Miskowiec <mailto:d.miskowiec@gsi.de> 2007
3 //=============================================================================
4 // single particle analyzer
5 //=============================================================================
11 #include <TParticlePDG.h>
13 #include "AliDEvent.h"
14 #include "AliDAnalSingle.h"
16 ClassImp(AliDAnalSingle)
18 //=============================================================================
19 AliDAnalSingle::AliDAnalSingle(Char_t *nam, Double_t emi, Double_t ema, Int_t pid) :
20 AliDAnal(nam), fPid(pid), fMass(0.0)
23 // emi and ema define the rapidity range for histograms
26 TParticlePDG *part = AliDAnal::fgPDG.GetParticle(fPid);
27 fMass = part? part->Mass() : 0;
29 double pi = TMath::Pi();
31 ax[0] = new TAxis(30,-1,1); ax[0]->SetTitle("vertex z");
32 ax[1] = new TAxis(80,emi,ema); ax[1]->SetTitle("eta");
33 ax[2] = new TAxis(90,-pi,pi); ax[2]->SetTitle("phi");
34 AliDHN *zep = new AliDHN("zep",3,ax);
35 for (int i=0; i<3; i++) delete ax[i];
37 ax[0] = new TAxis(20,0,1); ax[0]->SetTitle("centrality");
38 ax[1] = new TAxis(80,emi,ema); ax[1]->SetTitle("y");
39 ax[2] = new TAxis(80,0,2); ax[2]->SetTitle("pt (GeV)");
40 AliDHN *cyp = new AliDHN("cyp",3,ax);
41 for (int i=0; i<3; i++) delete ax[i];
43 ax[0] = new TAxis(10,emi,ema); ax[0]->SetTitle("eta");
44 ax[1] = new TAxis(150,0,3); ax[1]->SetTitle("p (GeV)");
45 ax[2] = new TAxis(150,0.5,3.5);ax[2]->SetTitle("sqrt(dedx (mips))");
46 AliDHN *epd = new AliDHN("epd",3,ax);
47 for (int i=0; i<3; i++) delete ax[i];
53 printf("%s object named %s created\n",ClassName(),GetName());
55 //=============================================================================
56 void AliDAnalSingle::Process(AliDEvent *ev)
58 // fill single particle histograms
60 AliDHN *zep = (AliDHN*) fHistos.At(0);
61 AliDHN *cyp = (AliDHN*) fHistos.At(1);
62 AliDHN *epd = (AliDHN*) fHistos.At(2);
63 for (int i=0; i<ev->NParticles(); i++) {
64 if (!ev->ParticleGood(i,fPid)) continue;
65 zep->Fill(ev->Zver(),ev->ParticleEta(i),ev->ParticlePhi(i),1.0);
66 double y = fMass>0? ev->ParticleY(i,fMass) : ev->ParticleEta(i);
67 cyp->Fill(ev->Centrality(),y,ev->ParticlePt(i),1.0);
68 epd->Fill(ev->ParticleEta(i),ev->ParticleP(i),sqrt(ev->ParticleDedx(i)),1.0);
71 //=============================================================================