method Ranmar is removed
[u/mrichter/AliRoot.git] / PMD / AliPMDesdanal.C
1 //********************************************************************
2 //     Example (very naive for the moment) of the data analysis 
3 //                    using the ESD classes
4 //     Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
5 //     Modified by Basanta K. Nandi for PMD analysis
6 //********************************************************************
7
8 #if !defined( __CINT__) || defined(__MAKECINT__)
9   #include <Riostream.h>
10   #include "TTree.h"
11   #include "TFile.h"
12   #include "TH1F.h"
13   #include "TCanvas.h"
14   #include "TStyle.h"
15   #include "TStopwatch.h"
16
17   #include "AliESD.h"
18 #endif
19
20 extern TStyle *gStyle;
21
22 Int_t AliPMDesdanal() { 
23    TStopwatch timer;
24
25    gStyle->SetOptStat(111110);
26    gStyle->SetOptFit(1);
27
28 //****** File with the ESD
29    TFile *ef=TFile::Open("AliESDcheck.root");
30    if (!ef || !ef->IsOpen()) {cerr<<"Can't AliESDs.root !\n"; return 1;}
31    AliESD* event = new AliESD;
32    TTree* tree = (TTree*) ef->Get("esdTree");
33    if (!tree) {cerr<<"no ESD tree found\n"; return 1;};
34    tree->SetBranchAddress("ESD", &event);
35
36    Int_t n=0;
37
38 //******* The loop over events
39    while (tree->GetEvent(n)) {
40      cout<<endl<<"Processing event number : "<<n++<<endl;
41
42
43      Int_t npmdcl=event->GetNumberOfPmdTracks();
44      cout<<"Number of PMD tracks : "<<npmdcl<<endl; 
45
46      //****** The loop over PMD clusters
47      while (npmdcl--) {
48        AliESDPmdTrack *pmdtr = event->GetPmdTrack(npmdcl);
49        
50        Int_t   det   = pmdtr->GetDetector();
51        Float_t theta = pmdtr->GetTheta();
52        Float_t phi   = pmdtr->GetPhi();
53        Float_t adc   = pmdtr->GetClusterADC();
54        Float_t pid   = pmdtr->GetClusterPID();
55        
56      }
57    }
58    timer.Stop();
59    timer.Print();
60
61    return 0;
62 }