coverity defects fixed
[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("AliESDs.root");
30    if (!ef || !ef->IsOpen()) {cerr<<"Can't AliESDs.root !\n"; return 1;}
31    AliESDEvent * event = new AliESDEvent;
32    TTree* tree = (TTree*) ef->Get("esdTree");
33    if (!tree) {cerr<<"no ESD tree found\n"; return 1;};
34    event->ReadFromTree(tree);
35    Int_t n=0;
36
37 //******* The loop over events
38    while (tree->GetEvent(n)) {
39      cout<<endl<<"Processing event number : "<<n++<<endl;
40
41
42      Int_t npmdcl=event->GetNumberOfPmdTracks();
43      cout<<"Number of PMD tracks : "<<npmdcl<<endl; 
44
45      //****** The loop over PMD clusters
46      while (npmdcl--) {
47        AliESDPmdTrack *pmdtr = event->GetPmdTrack(npmdcl);
48        
49        Int_t   det   = pmdtr->GetDetector(); 
50        Float_t clsX  = pmdtr->GetClusterX();
51        Float_t clsY  = pmdtr->GetClusterY();
52        Float_t clsZ  = pmdtr->GetClusterZ();
53        Float_t ncell = pmdtr->GetClusterCells();
54        Float_t adc   = pmdtr->GetClusterADC();
55        Float_t pid   = pmdtr->GetClusterPID();
56        
57      }
58    }
59
60    delete event;
61
62    timer.Stop();
63    timer.Print();
64
65    return 0;
66 }