]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PWG2/UNICOR/AliDEventAliceESD.h
Added seperate methods to write histograms to a file when not using the task framework
[u/mrichter/AliRoot.git] / PWG2 / UNICOR / AliDEventAliceESD.h
1 // Author: Dariusz Miskowiec <mailto:d.miskowiec@gsi.de> 2008
2
3 #ifndef ALIDEVENTALICEESD_H
4 #define ALIDEVENTALICEESD_H
5
6 #include <cmath>
7 #include "TVector2.h"
8 #include "AliDEvent.h"
9 #include "AliESDEvent.h"
10 #include "AliESDVertex.h"
11
12 //=============================================================================
13 class AliDEventAliceESD : public AliDEvent, public AliESDEvent {
14
15  public:
16               AliDEventAliceESD();
17   virtual    ~AliDEventAliceESD();
18   Double_t    Etamin() const {return -0.75;}
19   Double_t    Etamax() const {return  0.75;}
20   void        AttachTree(TTree *tr) {ReadFromTree(tr);}
21   Bool_t      Good() const;
22   Double_t    Centrality() {return 0.9999*exp(-NParticles()/20.0);} // OK for pp
23   void        RP(Double_t &qx, Double_t &qy) const {AliDEvent::RP(qx,qy,2);}
24   Double_t    RPphi() const {Double_t qx,qy; RP(qx,qy); return atan2(qy,qx);}
25   Double_t    Zver() const {return AliESDEvent::GetPrimaryVertex()->GetZv()/10.0;}
26   Int_t       NParticles() const {return AliESDEvent::GetNumberOfTracks();}
27
28   Bool_t      ParticleGood(Int_t i, Int_t pidi=0) const;
29   Double_t    ParticleP(Int_t i)     const {return AliESDEvent::GetTrack(i)->GetTPCInnerParam()->P();}
30   Double_t    ParticleTheta(Int_t i) const {return AliESDEvent::GetTrack(i)->GetTPCInnerParam()->Theta();}
31   Double_t    ParticlePhi(Int_t i)   const {return TVector2::Phi_mpi_pi(AliESDEvent::GetTrack(i)->GetTPCInnerParam()->Phi());}
32   Double_t    ParticleDedx(Int_t i)  const {return AliESDEvent::GetTrack(i)->GetTPCsignal()/47.0;}
33   Bool_t      PairGood(Double_t p0, Double_t the0, Double_t phi0, 
34                        Double_t p1, Double_t the1, Double_t phi1) const;
35   // alternative: GetTPCInnerParam, GetConstrainedParam
36   ClassDef(AliDEventAliceESD,0)
37 };
38 #endif 
39 //=============================================================================