]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PHOS/PetiteMacrodeNuit.C
Bayesian weight for PID calculated for time of flight info
[u/mrichter/AliRoot.git] / PHOS / PetiteMacrodeNuit.C
1 void PetiteMacroDeNuit()
2 {
3 Int_t evt = 0 ; 
4 RecAna * t = new RecAna("junk.root"); 
5 t->GetEvent(evt);   
6 TObjArray * lp = t->PHOSPpsdRP ; 
7  cout << "Tree macro = " << lp << endl ; 
8 for (int i = 0 ; i < lp->GetEntries() ; i++ ) {
9   AliPHOSPpsdRecPoint * rpp = (AliPHOSPpsdRecPoint *)lp->At(i)   ; 
10   rpp.Print(); 
11 }
12 TObjArray * le = t->PHOSEmcRP ; 
13 for (int i = 0 ; i < le->GetEntries() ; i++ ) {
14   AliPHOSEmcRecPoint * rp = (AliPHOSEmcRecPoint *)le->At(i)   ; 
15   rp->Print(); 
16 }
17
18 AliPHOSIndexToObject * please = AliPHOSIndexToObject::GetInstance() ;
19
20 for (int i = 0 ; i < (t->PHOSTS_-1) ; i++) {
21   cout << "TrackSegment # " << i << endl 
22        << "====================" << endl ; 
23   int index = t->PHOSTS_fEmcRecPoint[i] ;
24   AliPHOSEmcRecPoint * emrp = (AliPHOSEmcRecPoint *) ( please->GimeRecPoint(index, TString("emc") ) ) ; 
25   emrp->Print() ; 
26   index = t->PHOSTS_fPpsdLowRecPoint[i] ;
27   AliPHOSPpsdRecPoint * ppsdl = (AliPHOSPpsdRecPoint *) ( please->GimeRecPoint(index, TString("ppsd") ) ) ; 
28   if (ppsdl)
29     ppsdl->Print() ; 
30   index = t->PHOSTS_fPpsdUpRecPoint[i] ;
31   AliPHOSPpsdRecPoint * ppsdu = (AliPHOSPpsdRecPoint *) ( please->GimeRecPoint(index, TString("ppsd") ) ) ; 
32   if (ppsdu)
33     ppsdu->Print() ; 
34 }
35 for (int i = 0 ; i < (t->PHOSRP_-1) ; i++) {
36   cout << "RecParticles # " << i << endl 
37        << "====================" << endl ; 
38   cout << "type = " << t->PHOSRP_fType[i] << " energy = " << t->PHOSRP_fE[i] << endl ;
39 }  
40 delete t ; 
41 }