Removing GetDebug and SetDebug from AliRun and AliModule. Using AliLog for the messages
[u/mrichter/AliRoot.git] / START / AliSTART.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id$ */
17
18 ///////////////////////////////////////////////////////////////////////////////
19 //                                                                           //
20 //  START (T-Zero) Detector                                            //
21 //  This class contains the base procedures for the START     //
22 //  detector                                                                 //
23 //                                                                           //
24 //Begin_Html
25 /*
26 <img src="gif/AliSTARTClass.gif">
27 </pre>
28 <br clear=left>
29 <font size=+2 color=red>
30 <p>The responsible person for this module is
31 <a href="mailto:Alla.Maevskaia@cern.ch">Alla Maevskaia</a>.
32 </font>
33 <pre>
34 */
35 //End_Html
36 //                                                                           //
37 //                                                                           //
38 ///////////////////////////////////////////////////////////////////////////////
39
40 #include <Riostream.h>
41
42 #include <TFile.h>
43 #include <TGeometry.h>
44 #include <TMath.h>
45 #include <TNode.h>
46 #include <TParticle.h>
47 #include <TRandom.h>
48 #include <TTUBE.h>
49 #include <TVirtualMC.h>
50 #include <AliESD.h>
51
52 #include "AliLog.h"
53 #include "AliMC.h"
54 #include "AliLoader.h"
55 #include "AliRun.h"
56 #include "TClonesArray.h"
57 #include "AliSTART.h"
58 #include "AliSTARTLoader.h"
59 #include "AliSTARTdigit.h"
60 #include "AliSTARThit.h"
61 #include "AliSTARTDigitizer.h"
62 #include "AliSTARTRawData.h"
63 #include "AliSTARTRecPoint.h"
64 #include "AliLog.h"
65
66 ClassImp(AliSTART)
67
68   //static  AliSTARTdigit *digits; 
69
70 //_____________________________________________________________________________
71 AliSTART::AliSTART()
72 {
73   //
74   // Default constructor for class AliSTART
75   //
76   fIshunt   = 1;
77   fHits     = 0;
78   fDigits   = 0;
79   fRecPoints = 0;
80 }
81  
82 //_____________________________________________________________________________
83 AliSTART::AliSTART(const char *name, const char *title)
84        : AliDetector(name,title)
85 {
86   //
87   // Standard constructor for START Detector
88   //
89
90   
91   //
92   // Initialise Hit array
93   fHits       = new TClonesArray("AliSTARThit",  405);
94   gAlice->GetMCApp()->AddHitList(fHits);
95   fDigits    = new AliSTARTdigit();
96   fRecPoints = new AliSTARTRecPoint();
97   fIshunt     =  1;
98   fIdSens   =  0;
99   SetMarkerColor(kRed);
100 }
101
102 //_____________________________________________________________________________
103 AliSTART::~AliSTART() {
104   
105   if (fHits) {
106     fHits->Delete();
107     delete fHits;
108     cout<<"  delete fHits; "<<endl;
109   }
110   /*
111   if (fDigits) {
112     fDigits->Delete();
113     delete fDigits;
114     cout<<" delete fDigits; "<<endl;
115   }
116   if (fRecPoints) {
117    fRecPoints ->Delete();
118     delete fRecPoints;
119     cout<<" delete fRecPoints; "<<endl;
120   }
121   */ 
122 }
123  
124 //_____________________________________________________________________________
125 void AliSTART::AddHit(Int_t track, Int_t *vol, Float_t *hits)
126 {
127   //
128   // Add a START hit
129   //
130   TClonesArray &lhits = *fHits;
131   new(lhits[fNhits++]) AliSTARThit(fIshunt,track,vol,hits);
132 }
133
134
135 //_____________________________________________________________________________
136
137 void AliSTART::AddDigit(Int_t besttimeright, Int_t besttimeleft, Int_t meantime, 
138                         Int_t timediff, TArrayI *sumMult,
139                         TArrayI *time, TArrayI *adc, TArrayI *timeAmp, TArrayI *adcAmp)
140 {
141   
142   //  Add a START digit to the list.
143  //
144   
145   if (!fDigits) {
146     fDigits = new AliSTARTdigit();
147   }
148   fDigits-> SetTimeBestRight(besttimeright);
149   fDigits->SetTimeBestLeft(besttimeleft);
150   fDigits-> SetMeanTime(meantime);
151   fDigits-> SetDiffTime(timediff);
152   fDigits-> SetSumMult(*sumMult);
153   fDigits->SetTime(*time);
154   fDigits->SetTimeAmp(*timeAmp);
155   fDigits->SetADC(*adc);
156   fDigits->SetADCAmp(*adcAmp);
157 }
158
159
160 //_____________________________________________________________________________
161 void AliSTART::BuildGeometry()
162 {
163   //
164   // Build simple ROOT TNode geometry for event display
165   //
166   TNode *node, *top;
167   const int kColorSTART  = 19;
168
169   top=gAlice->GetGeometry()->GetNode("alice");
170
171   // START define the different volumes
172   new TRotMatrix("rotx999","rot999",  90,0,90,90,180,0);
173
174   new TTUBE("S_0ST1","START  volume 1","void",5.,10.7,5.3);
175   top->cd();
176   node = new TNode("0ST1","0ST01","S_0ST1",0,0,-69.7,"");
177   node->SetLineColor(kColorSTART);
178   fNodes->Add(node);
179
180   new TTUBE("S_0ST2","START volume 2","void",5.,10.7,5.3);
181   top->cd();
182   node = new TNode("0ST2","0ST2","S_0ST2",0,0,350,"rotx999");
183   node->SetLineColor(kColorSTART);
184   fNodes->Add(node);
185 }
186  
187 //_____________________________________________________________________________
188 Int_t AliSTART::DistanceToPrimitive(Int_t /*px*/, Int_t /*py*/)
189 {
190   //
191   // Calculate the distance from the mouse to the START on the screen
192   // Dummy routine
193   //
194   return 9999;
195 }
196  
197 //-------------------------------------------------------------------------
198 void AliSTART::Init()
199 {
200   //
201   // Initialis the START after it has been built
202   Int_t i;
203   //
204   if(AliLog::GetGlobalDebugLevel()>0) {
205     printf("\n%s: ",ClassName());
206     for(i=0;i<35;i++) printf("*");
207     printf(" START_INIT ");
208     for(i=0;i<35;i++) printf("*");
209     printf("\n%s: ",ClassName());
210     //
211     // Here the START initialisation code (if any!)
212     for(i=0;i<80;i++) printf("*");
213     printf("\n");
214   }
215 }
216
217 //---------------------------------------------------------------------------
218 void AliSTART::MakeBranch(Option_t* option)
219 {
220   //
221 // Create Tree branches for the START.
222
223  // Options:
224   //
225   //    H          Make a branch of TClonesArray of AliSTARTHit's
226   //    D          Make a branch of TClonesArray of AliSTARTDigit's
227   //
228   //    R         Make a branch of  AliSTARTRecPoints
229   //
230   char branchname[20];
231   sprintf(branchname,"%s",GetName());
232
233   const char *cH = strstr(option,"H");
234   const char *cD = strstr(option,"D");
235   const char *cR = strstr(option,"R");
236
237     if (cH && fLoader->TreeH())
238   {
239      if (fHits == 0x0) fHits  = new TClonesArray("AliSTARThit",  405);
240      AliDetector::MakeBranch(option);
241   } 
242     
243     
244   if (cD && fLoader->TreeD())
245     {
246       if (fDigits == 0x0) fDigits  = new AliSTARTdigit();
247       //     MakeBranchInTree(fLoader->TreeD(), branchname,
248       //                       &fDigits, 405, 0);
249       fLoader->TreeD()->Branch(branchname,"AliSTARTdigit",&fDigits,405,1);
250     } 
251   if (cR && fLoader->TreeR())
252     {
253       if (fRecPoints == 0x0) fRecPoints  = new AliSTARTRecPoint();
254       MakeBranchInTree(fLoader->TreeR(), branchname,
255                        &fRecPoints, 405, 0);
256     } 
257   
258 }    
259
260 //_____________________________________________________________________________
261 void AliSTART::ResetHits()
262 {
263   AliDetector::ResetHits();
264   
265 }
266 //____________________________________________________________________
267 void AliSTART::ResetDigits()
268 {
269   //
270   // Reset number of digits and the digits array for this detector
271   //
272   if (fDigits) fDigits->Clear();
273 }
274
275 //_____________________________________________________________________________
276 void AliSTART::SetTreeAddress()
277 {
278
279   TTree    *treeH;
280   treeH = TreeH();
281   
282   if (treeH)
283     {
284       if (fHits == 0x0) fHits  = new TClonesArray("AliSTARThit",  405);
285     }
286     
287   AliDetector::SetTreeAddress();
288   TTree *treeD = fLoader->TreeD();
289   if (treeD) {
290     if (fDigits == 0x0)  fDigits  = new AliSTARTdigit();
291     TBranch* branch = treeD->GetBranch ("START");
292     if (branch) branch->SetAddress(&fDigits);
293   }
294
295   TTree *treeR = fLoader->TreeR();
296   if (treeR) {
297     if (fRecPoints == 0x0) fRecPoints  = new  AliSTARTRecPoint()  ;
298     TBranch* branch = treeR->GetBranch ("START");
299     if (branch) branch->SetAddress(&fRecPoints);
300   }
301  
302 }
303
304
305 //_____________________________________________________________________________
306 void AliSTART::MakeBranchInTreeD(TTree *treeD, const char *file)
307 {
308     //
309     // Create TreeD branches for the FMD
310     //
311     const Int_t kBufferSize = 4000;
312     char branchname[20];
313     sprintf(branchname,"%s",GetName());
314     if(treeD)
315      {
316        MakeBranchInTree(treeD,  branchname,&fDigits, kBufferSize, file);
317      }
318 }
319
320 //_____________________________________________________________________________
321 AliDigitizer* AliSTART::CreateDigitizer(AliRunDigitizer* manager) const
322 {
323   return new AliSTARTDigitizer(manager);
324 }
325 //____________________________________________________________________________
326 void AliSTART::Digits2Raw()
327 {
328 //
329 // Starting from the START digits, writes the Raw Data objects
330 //
331 //  AliSTARTLoader* pStartLoader = (AliSTARTLoader*)fLoader;
332   fLoader ->LoadDigits("read");
333   TTree* treeD = fLoader->TreeD();
334   if (!treeD) {
335     AliError("no digits tree");
336     return;
337   }
338   if (fDigits == 0x0)  fDigits  = new AliSTARTdigit();
339   
340   TBranch *branch = treeD->GetBranch("START");
341   if (branch) {
342     branch->SetAddress(&fDigits);
343   }else{
344     AliError("Branch START DIGIT not found");
345     exit(111);
346   } 
347   AliSTARTRawData rawWriter;
348   rawWriter.SetVerbose(0);
349   
350   AliDebug(2,Form(" Formatting raw data for START "));
351   branch->GetEntry(0);
352   //  rawWriter.RawDataSTART(treeD->GetBranch("START"));
353   rawWriter.RawDataSTART(fDigits);
354   
355   
356   fLoader->UnloadDigits();
357   
358 }