Double check if SM is running added. Some redundant output removed from SM
[u/mrichter/AliRoot.git] / ITS / AliITSHits2Digits.C
1 #if !defined(__CINT__) || defined(__MAKECINT__)
2
3 #include "TROOT.h"
4 #include "Riostream.h"
5 #include "TDatime.h"
6 #include "TClassTable.h"
7 #include "AliRun.h"
8 #include "AliRunDigitizer.h"
9 #include "AliITS.h"
10 #include "AliITSDetTypeSim.h"
11 #include "AliITSLoader.h"
12 #include "AliITSresponseSDD.h"
13 #include "TStopwatch.h"
14
15 #endif
16
17 //#define DEBUG
18 Int_t AliITSHits2Digits(TString inFile = "galice.root"){
19     // Standard ITS Hits to Digits, excluding creation of SDigits.
20
21     // Dynamically link some shared libs
22    if (gClassTable->GetID("AliRun") < 0) {
23       gROOT->ProcessLine(".x $(ALICE_ROOT)/macros/loadlibs.C");
24     }else if (gAlice){
25       delete AliRunLoader::Instance();
26       delete gAlice;
27       gAlice=0;
28      } 
29
30    AliRunLoader* rl = AliRunLoader::Open(inFile.Data());
31     if (rl == 0x0)
32      {
33       cerr<<"AliITSHits2Digits.C : Can not open session RL=NULL"
34            << endl;
35        return 3;
36      }
37      
38     Int_t retval = rl->LoadgAlice();
39     if (retval)
40      {
41       cerr<<"AliITSHits2Digits.C : LoadgAlice returned error"
42            << endl;
43        return 3;
44      }
45     gAlice=rl->GetAliRun();
46     AliITSLoader* gime = (AliITSLoader*)rl->GetLoader("ITSLoader");
47     if (gime == 0x0)
48      {
49       cerr<<"AliITSHits2Digits.C : can not get ITS loader"
50            << endl;
51      }
52     AliITS *ITS = (AliITS*)gAlice->GetDetector("ITS");      
53     if (!ITS) {
54         cerr<<"AliITSHits2Digit.C : AliITS object not found on file"
55             << endl;
56         return 3;
57     }  // end if !ITS
58
59     if(!(ITS->GetITSgeom())){
60         cerr << " AliITSgeom not found. Can't digitize without it." << endl;
61         return 4;
62     } // end if
63
64     
65
66     TStopwatch timer;
67     Int_t evNumber1 = 0;
68     Int_t evNumber2 = rl->GetNumberOfEvents();
69     timer.Start();
70     retval = gime->LoadHits();
71     if (retval)
72      {
73       cerr<<"AliITSHits2Digits.C : ITSLoader::LoadHits returned error"
74            << endl;
75        return 3;
76      }
77
78     retval = gime->LoadDigits("recreate");
79     if (retval)
80      {
81       cerr<<"AliITSHits2Digits.C : ITSLoader::LoadDigits returned error"
82            << endl;
83        return 3;
84      }
85     for(Int_t nevent = evNumber1; nevent < evNumber2; nevent++){
86         // cout<<"Producing Digits for event n."<<nevent<<endl;
87
88         rl->GetEvent(nevent);
89         if(!gime->TreeD()){ 
90             cout << "Having to create the Digits Tree." << endl;
91             gime->MakeTree("D");
92         } // end if creating digits tree
93         ITS->MakeBranch("D");
94         ITS->SetTreeAddress();   
95         ITS->Hits2Digits();
96     } // end for nevent
97     timer.Stop();
98     timer.Print();
99
100     delete rl;
101     return 0;
102 }