Fix for FMD DA
[u/mrichter/AliRoot.git] / ITS / AliITSHits2FastRecPoints.C
1 #if !defined(__CINT__) || defined(__MAKECINT__)
2
3 #include <Riostream.h>
4 #include<TClassTable.h>
5 #include <TError.h>
6 #include <TInterpreter.h>
7 #include <TStopwatch.h>
8 #include <TString.h>
9 #include "AliRun.h"
10 #include "AliRunLoader.h"
11 #include "AliITS.h"
12 #include "AliITSInitGeometry.h"
13 #include "AliITSgeom.h"
14 #include "AliITSLoader.h"
15 #include "AliITSsimulationFastPoints.h"
16 #include "AliGeomManager.h"
17 #endif
18
19 /*
20 $Id$ 
21 */
22
23
24 void AliITSHits2FastRecPoints (Int_t evNumber1=0,Int_t evNumber2=0, TString inFile = "galice.root", Int_t nsignal=25, Int_t size=-1) 
25 {
26   /////////////////////////////////////////////////////////////////////////
27   //   
28   //   This macro creates fast recpoints
29   //   
30   /////////////////////////////////////////////////////////////////////////
31
32
33   // Dynamically link some shared libs
34
35   if (gClassTable->GetID("AliRun") < 0) {
36     gInterpreter->ExecuteMacro("loadlibs.C");
37   }
38   else { 
39     if(gAlice){
40       delete AliRunLoader::Instance();
41       delete gAlice;
42       gAlice=0;
43     }
44   }
45
46
47   // Get geometry
48   AliGeomManager::LoadGeometry("geometry.root");
49
50   // Connect the Root Galice file containing Geometry, Kine and Hits
51   AliRunLoader* rl = AliRunLoader::Open(inFile.Data());
52   if (rl == 0x0)
53     {
54       Error("AliITSHits2FastRecPoints.C","Can not open session RL=NULL");
55       return;
56     }
57      
58   Int_t retval = rl->LoadgAlice();
59   if (retval)
60     {
61       Error("AliITSHits2FastRecPoints.C","LoadgAlice returned error");
62       delete rl;
63       return;
64     }
65   gAlice=rl->GetAliRun();
66   rl->LoadHeader();
67   retval = rl->LoadKinematics();
68   if (retval)
69     {
70       Error("AliITSHits2FastRecPoints.C","LoadKinematics returned error");
71       delete rl;
72       return;
73     }
74
75   AliITSInitGeometry initgeom;
76   AliITSgeom *geom = initgeom.CreateAliITSgeom();
77   printf("Geometry name: %s \n",(initgeom.GetGeometryName()).Data());
78     
79   AliITSLoader* gime = (AliITSLoader*)rl->GetLoader("ITSLoader");
80   if (gime == 0x0)
81     {
82       ::Error("AliITSHits2FastRecPoints.C","can not get ITS loader");
83       delete rl;
84       return;
85     }
86   gime->SetITSgeom(geom);
87   retval = gime->LoadHits("read");
88   if (retval)
89     {
90       ::Error("AliITSHits2FastRecPoints.C","LoadHits returned error");
91       delete rl;
92       return;
93     }
94   gime->SetRecPointsFileName("ITS.FastRecPoints.root"); 
95   retval = gime->LoadRecPoints("update");
96   if (retval)
97     {
98       ::Error("AliITSHits2FastRecPoints.C","LoadRecPoints returned error");
99       delete rl;
100       return;
101     }
102     
103    
104
105  
106   AliITS *ITS  = (AliITS*) gAlice->GetModule("ITS");
107   if (!ITS) return;
108
109   // Set the simulation model
110
111    
112
113   //
114   // Event Loop
115   //
116
117   TStopwatch timer;
118
119   cout << "Creating fast reconstructed points from hits for the ITS..." << endl;
120   AliITSDetTypeSim* dettyp = new AliITSDetTypeSim();
121   dettyp->SetLoader(gime);
122   dettyp->SetITSgeom(geom);
123   ITS->SetDetTypeSim(dettyp);
124   for (Int_t i=0;i<3;i++) {
125     ITS->SetSimulationModel(i,new AliITSsimulationFastPoints());
126   }
127
128   
129   for (int ev=evNumber1; ev<= evNumber2; ev++) {
130     cout << "...working on event "<< ev << " ..." << endl;
131     Int_t nparticles = gAlice->GetEvent(ev);
132     cout << "event         " <<ev<<endl;
133     cout << "nparticles  " <<nparticles<<endl;
134     rl->GetEvent(ev);
135     //if(gime->TreeR() == 0x0) gime->MakeTree("R");
136     
137     if (ev < evNumber1) continue;
138     if (nparticles <= 0) return;
139
140     Int_t bgr_ev=Int_t(ev/nsignal);
141     timer.Start();
142     ITS->HitsToFastRecPoints(ev,bgr_ev,size," ","All"," ");
143     timer.Stop(); timer.Print();
144   } // event loop 
145
146   delete rl;
147 }
148