Updated for coding convenstions compilation errors and warnings and the like.
[u/mrichter/AliRoot.git] / ITS / AliITSv5.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 /*
17 $Log$
18 Revision 1.14.4.9  2000/06/12 19:14:51  barbera
19 Remove partical transision to new Config.C calling convension. Bug Found.
20
21 Revision 1.14.4.8  2000/06/12 18:15:52  barbera
22 fixed posible compilation errors on HP unix. Modifided default constructor
23 for use with new calling requirements.
24
25 Revision 1.14.4.7  2000/06/11 20:37:59  barbera
26 coding convenstion update.
27
28 Revision 1.14.4.4  2000/05/19 10:10:21  nilsen
29 fix for bug with HP and Sun unix + fix for event display in ITS-working branch
30
31 Revision 1.14.4.3  2000/03/04 23:46:38  nilsen
32 Fixed up the comments/documentation.
33
34 Revision 1.14.4.2  2000/03/02 21:53:02  nilsen
35 To make it compatable with the changes in AliRun/AliModule.
36
37 Revision 1.14.4.1  2000/01/12 19:03:33  nilsen
38 This is the version of the files after the merging done in December 1999.
39 See the ReadMe110100.txt file for details
40
41 Revision 1.14  1999/10/22 08:16:49  fca
42 Correct destructors, thanks to I.Hrivnacova
43
44 Revision 1.13  1999/10/06 10:15:19  fca
45 Correct bug in allocation of layer name and add destructor
46
47 Revision 1.12  1999/10/05 08:05:09  fca
48 Minor corrections for uninitialised variables.
49
50 Revision 1.11  1999/09/29 09:24:20  fca
51 Introduction of the Copyright and cvs Log
52
53 */
54
55 ///////////////////////////////////////////////////////////////////////////////
56 //
57 //  Inner Traking System version 5
58 //  This class contains the base procedures for the Inner Tracking System
59 //
60 // Authors: R. Barbera, B. S. Nilsen.
61 // version 5.
62 // Created September 17 1999.
63 //
64 ///////////////////////////////////////////////////////////////////////////////
65
66 // See AliITSv5::StepManager().
67 #define ALIITSPRINTGEOM 0 // default. don't print out gemetry information
68 //#define ALIITSPRINTGEOM 1 // print out geometry information
69
70 #include <stdio.h>
71 #include <stdlib.h>
72 #include <TMath.h>
73 #include <TGeometry.h>
74 #include <TNode.h>
75 #include <TTUBE.h>
76 #include <TFile.h>    // only required for Tracking function?
77 #include <TCanvas.h>
78 #include <TObjArray.h>
79 #include <TClonesArray.h>
80
81 #include "AliRun.h"
82 #include "TSystem.h"
83 #if ALIITSPRINTGEOM==1
84 #include "../TGeant3/TGeant3.h"
85 #endif
86 #include "AliITShit.h"
87 #include "AliITS.h"
88 #include "AliITSv5.h"
89 #include "AliITSgeom.h"
90
91 ClassImp(AliITSv5)
92  
93 //_____________________________________________________________________________
94   AliITSv5::AliITSv5()/*:AliITS("ITS","TDR version")*/{
95 ////////////////////////////////////////////////////////////////////////
96 //    Standard default constructor for the ITS version 5.
97 ////////////////////////////////////////////////////////////////////////
98     fId5N = 6;
99     fId5Name = new char*[fId5N];
100     fId5Name[0] = "ITS1";
101     fId5Name[1] = "ITS2";
102     fId5Name[2] = "ITS3";
103     fId5Name[3] = "ITS4";
104     fId5Name[4] = "ITS5";
105     fId5Name[5] = "ITS6";
106
107     fEuclidMaterial = "$ALICE_ROOT/Euclid/ITSgeometry_5.tme";
108     fEuclidGeometry = "$ALICE_ROOT/Euclid/ITSgeometry_5.euc";
109     printf("Created ITS TDR Detailed version reading geometry from file\n");
110 }
111 //____________________________________________________________________________
112 AliITSv5::AliITSv5(const AliITSv5 &source){
113 ////////////////////////////////////////////////////////////////////////
114 //     Copy Constructor for ITS version 5.
115 ////////////////////////////////////////////////////////////////////////
116     if(&source == this) return;
117     this->fId5N = source.fId5N;
118     this->fId5Name = new char*[fId5N];
119     Int_t i;
120     for(i=0;i<6;i++) strcpy(this->fId5Name[i],source.fId5Name[i]);
121     return;
122 }
123 //_____________________________________________________________________________
124 AliITSv5& AliITSv5::operator=(const AliITSv5 &source){
125 ////////////////////////////////////////////////////////////////////////
126 //    Assignment operator for the ITS version 1.
127 ////////////////////////////////////////////////////////////////////////
128   if(&source == this) return *this;
129   this->fId5N = source.fId5N;
130   this->fId5Name = new char*[fId5N];
131   Int_t i;
132   for(i=0;i<6;i++) strcpy(this->fId5Name[i],source.fId5Name[i]);
133   return *this;
134 }
135 //_____________________________________________________________________________
136 AliITSv5::~AliITSv5() {
137 ////////////////////////////////////////////////////////////////////////
138 //    Standard destructor for the ITS version 5.
139 ////////////////////////////////////////////////////////////////////////
140   delete [] fId5Name;
141 }
142 //_____________________________________________________________________________
143 AliITSv5::AliITSv5(const char *name, const char *title) : AliITS(name, title){
144 ////////////////////////////////////////////////////////////////////////
145 //    Standard constructor for the ITS version 5.
146 ////////////////////////////////////////////////////////////////////////
147     fId5N = 6;
148     fId5Name = new char*[fId5N];
149     fId5Name[0] = "ITS1";
150     fId5Name[1] = "ITS2";
151     fId5Name[2] = "ITS3";
152     fId5Name[3] = "ITS4";
153     fId5Name[4] = "ITS5";
154     fId5Name[5] = "ITS6";
155
156     fEuclidMaterial = "$ALICE_ROOT/Euclid/ITSgeometry_5.tme";
157     fEuclidGeometry = "$ALICE_ROOT/Euclid/ITSgeometry_5.euc";
158     printf("Created ITS TDR Detailed version reading geometry from file\n");
159 }
160 //_____________________________________________________________________________
161 void AliITSv5::BuildGeometry(){
162 ////////////////////////////////////////////////////////////////////////
163 //    Geometry builder for the ITS version 5.
164 ////////////////////////////////////////////////////////////////////////
165   //
166   // Build ITS TNODE geometry for event display using detailed geometry.
167   // This function builds a simple ITS geometry used by the ROOT macro
168   // ITSdisplay.C.
169
170   TNode *top;
171   TNode *nd;
172   //const Int_t kColorITSSPD=kRed;
173   //const Int_t kColorITSSDD=kGreen;
174   const Int_t kColorITSSSD=kBlue;
175   //
176   top=gAlice->GetGeometry()->GetNode("alice");
177   AliITSgeom  *gm = this->GetITSgeom();
178
179   Int_t       lay,lad,det,i;
180   Text_t      name[10];
181   Float_t     xg[3];
182   Float_t     rt[9];
183   Double_t    rtd[9];
184   TBRIK       *box;
185   TRotMatrix  *rm;
186   //TCanvas     *c1 = new TCanvas("c1","ITS");
187
188   for(lay=1;lay<=2;lay++)
189    for(lad=1;lad<=gm->GetNladders(lay);lad++)
190     for(det=1;det<=gm->GetNdetectors(lay);det++){
191           try {
192               box  = new TBRIK ("ActiveSPD","Active volume of SPD","SPD SI DET",
193                                     0.64,0.0075,4.19); 
194           } catch (...) {
195               cout << "EXCEPTION in box = new TBRIK" << endl;
196               return;
197           }
198           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
199           gm->GetRotMatrix(lay,lad,det,rt);
200           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
201           for(i=0;i<9;i++) rtd[i] = rt[i];
202           try {
203                 rm  = new TRotMatrix(name,name,rtd);
204           } catch (...) {
205                 cout << "EXCEPTION in   new TRotMatrix" << endl;
206                 return;
207           }
208          top->cd();
209           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); 
210          try {
211               nd  = new TNode("SPD"," ",box,xg[0],xg[1],xg[2],rm);
212          } catch (...) {
213               cout << "EXCEPTION in new TNode" << endl;
214               return;
215          }
216          nd->SetLineColor(kColorITSSSD);
217          fNodes->Add(nd);
218     }
219
220   for(lay=3;lay<=3;lay++)
221    for(lad=1;lad<=gm->GetNladders(lay);lad++)
222     for(det=1;det<=gm->GetNdetectors(lay);det++){
223           try {
224               box  = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
225                                     3.5,0.014,3.763); 
226           } catch (...) {
227               cout << "EXCEPTION in box = new TBRIK" << endl;
228               return;
229           }
230           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
231           gm->GetRotMatrix(lay,lad,det,rt);
232           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
233           for(i=0;i<9;i++) rtd[i] = rt[i];
234           try {
235                 rm  = new TRotMatrix(name,name,rtd);
236           } catch (...) {
237                 cout << "EXCEPTION in   new TRotMatrix" << endl;
238                 return;
239           }
240          top->cd();
241           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); 
242          try {
243               nd  = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
244          } catch (...) {
245               cout << "EXCEPTION in new TNode" << endl;
246               return;
247          }
248          nd->SetLineColor(kColorITSSSD);
249          fNodes->Add(nd);
250     }
251
252   for(lay=4;lay<=4;lay++)
253    for(lad=1;lad<=gm->GetNladders(lay);lad++)
254     for(det=1;det<=gm->GetNdetectors(lay);det++){
255           try {
256               box  = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
257                                     3.5,0.014,3.763); 
258           } catch (...) {
259               cout << "EXCEPTION in box = new TBRIK" << endl;
260               return;
261           }
262           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
263           gm->GetRotMatrix(lay,lad,det,rt);
264           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
265           for(i=0;i<9;i++) rtd[i] = rt[i];
266           try {
267                 rm  = new TRotMatrix(name,name,rtd);
268           } catch (...) {
269                 cout << "EXCEPTION in   new TRotMatrix" << endl;
270                 return;
271           }
272          top->cd();
273           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); 
274          try {
275               nd  = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
276          } catch (...) {
277               cout << "EXCEPTION in new TNode" << endl;
278               return;
279          }
280          nd->SetLineColor(kColorITSSSD);
281          fNodes->Add(nd);
282     }
283  for(lay=5;lay<=5;lay++)
284    for(lad=1;lad<=gm->GetNladders(lay);lad++)
285     for(det=1;det<=gm->GetNdetectors(lay);det++){
286           try {
287               box  = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
288                                     3.65,0.015,2.0); 
289           } catch (...) {
290               cout << "EXCEPTION in box = new TBRIK" << endl;
291               return;
292           }
293           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
294           gm->GetRotMatrix(lay,lad,det,rt);
295           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
296           for(i=0;i<9;i++) rtd[i] = rt[i];
297           try {
298                 rm  = new TRotMatrix(name,name,rtd);
299           } catch (...) {
300                 cout << "EXCEPTION in   new TRotMatrix" << endl;
301                 return;
302           }
303          top->cd();
304           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); 
305          try {
306               nd  = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
307          } catch (...) {
308               cout << "EXCEPTION in new TNode" << endl;
309               return;
310          }
311          nd->SetLineColor(kColorITSSSD);
312          fNodes->Add(nd);
313     }
314
315  for(lay=6;lay<=6;lay++)
316    for(lad=1;lad<=gm->GetNladders(lay);lad++)
317     for(det=1;det<=gm->GetNdetectors(lay);det++){
318           try {
319               box  = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
320                                     3.65,0.015,2.0); 
321           } catch (...) {
322               cout << "EXCEPTION in box = new TBRIK" << endl;
323               return;
324           }
325
326           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]); 
327           gm->GetRotMatrix(lay,lad,det,rt);
328           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
329           for(i=0;i<9;i++) rtd[i] = rt[i];
330           try {
331                 rm  = new TRotMatrix(name,name,rtd);
332           } catch (...) {
333                 cout << "EXCEPTION in   new TRotMatrix" << endl;
334                 return;
335           }
336          top->cd();
337           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); 
338          try {
339               nd  = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
340          } catch (...) {
341               cout << "EXCEPTION in new TNode" << endl;
342               return;
343          }
344          nd->SetLineColor(kColorITSSSD);
345          fNodes->Add(nd);
346     }
347
348
349
350 //_____________________________________________________________________________
351 void AliITSv5::CreateMaterials(){
352 ////////////////////////////////////////////////////////////////////////
353 //     Read a file containing the materials for the ITS version 5.
354 ////////////////////////////////////////////////////////////////////////
355     char *filtmp;
356
357   filtmp = gSystem->ExpandPathName(fEuclidMaterial.Data());
358
359   FILE *file = fopen(filtmp,"r");
360   if(file) {
361     fclose(file);
362     ReadEuclidMedia(filtmp);
363   } else {
364     Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",filtmp);
365     exit(1);
366   } // end if(file)
367 }
368 //_____________________________________________________________________________
369 void AliITSv5::CreateGeometry(){
370 //////////////////////////////////////////////////////////////////////
371 //    This is the geometry used for the ITS Pre-TDR and comes from an 
372 // Euclid to Geant conversion. The only difference
373 // is in the details of the ITS supports. The detectors elements, 
374 // detector numbering, and local and global reference frames are shown in
375 //  the following figures.
376 //Begin_Html
377 /*
378 <img src="picts/ITS/its1+2_convention_front_5.gif">
379 </pre>
380 <br clear=left>
381 <font size=+2 color=red>
382 <p>This shows the front view of the SPDs.
383 </font>
384 <pre>
385 <img src="picts/ITS/its1+2_convention_side_5.gif">
386 </pre>
387 <br clear=left>
388 <font size=+2 color=red>
389 <p>This shows the perspective view of the SPDs.
390 </font>
391 <img src="picts/ITS/its1+2_tree.gif">
392 </pre>
393 <br clear=left>
394 <font size=+2 color=red>
395 <p>This shows the geometry Tree for the SPDs.
396 </font>
397 <pre>
398
399 <pre>
400 <img src="picts/ITS/its3+4_convention_front_5.gif">
401 </pre>
402 <br clear=left>
403 <font size=+2 color=red>
404 <p>This shows the front view of the SDDs.
405 </font>
406 <pre>
407 <img src="picts/ITS/its3+4_convention_side_5.gif">
408 </pre>
409 <br clear=left>
410 <font size=+2 color=red>
411 <p>This shows the perspective view of the SDDs.
412 </font>
413 <img src="picts/ITS/its3+4_tree.gif">
414 </pre>
415 <br clear=left>
416 <font size=+2 color=red>
417 <p>This shows the geometry Tree for the SDDs.
418 </font>
419 <pre>
420
421 <pre>
422 <img src="picts/ITS/its5+6_convention_front_5.gif">
423 </pre>
424 <br clear=left>
425 <font size=+2 color=red>
426 <p>This shows the front view of the SSDs.
427 </font>
428 <pre>
429 <img src="picts/ITS/its5+6_convention_side_5.gif">
430 </pre>
431 <br clear=left>
432 <font size=+2 color=red>
433 <p>This shows the perspective view of the SSDs.
434 </font>
435 <pre>
436 <img src="picts/ITS/its5+6_tree.gif">
437 </pre>
438 <br clear=left>
439 <font size=+2 color=red>
440 <p>This shows the geometry Tree for the SSDs.
441 </font>
442 <pre>
443
444
445 <img src="picts/ITS/its_layer1-6_2.gif">
446 </pre>
447 <br clear=left>
448 <font size=+2 color=red>
449 <p>This shows the front view of the whole ITS..
450 </font>
451 <pre>
452
453 <img src="picts/ITS/its_layer1-6_1.gif">
454 </pre>
455 <br clear=left>
456 <font size=+2 color=red>
457 <p>This shows the perspective view of the whole ITS..
458 </font>
459 <pre>
460
461 <img src="picts/ITS/its1-6_tree.gif">
462 </pre>
463 <br clear=left>
464 <font size=+2 color=red>
465 <p>This shows the geometry Tree for the whole ITS.
466 </font>
467 <pre>
468 */
469 //End_Html
470 //
471 //
472 //      Here are shown the details of the ITS support cones and services.
473 // First is a GEANT tree showing the organization of all of the volumes
474 // that make up the ITS supports and services.
475 //Begin_Html
476 /*
477 <img src="picts/ITS/supports_tree.gif">
478  */
479 //End_Html
480 //     What follows are a number of figures showing what these support
481 // structures look like.
482 //Begin_Html
483 /*
484
485 <img src="picts/ITS/supports_3.gif">
486 </pre>
487 <br clear=left>
488 <font size=+2 color=red>
489 <p>This shows the geometry of the supports for the Drift and Strip layers only.
490 </font>
491 <pre>
492
493 <img src="picts/ITS/supports_2.gif">
494 </pre>
495 <br clear=left>
496 <font size=+2 color=red>
497 <p>This shows the geometry of the supports for the Drift and Strip layers in front cut out.
498 </font>
499 <pre>
500
501 <img src="picts/ITS/supports_1.gif">
502 </pre>
503 <br clear=left>
504 <font size=+2 color=red>
505 <p>This shows the geometry of the supports for the Drift and Strip layers in a back cut out..
506 </font>
507 <pre>
508
509 <img src="picts/ITS/suppssd.gif">
510 </pre>
511 <br clear=left>
512 <font size=+2 color=red>
513 <p>This shows the geometry for the Strip layers supports.
514 </font>
515 <pre>
516
517 <img src="picts/ITS/suppsdd.gif">
518 </pre>
519 <br clear=left>
520 <font size=+2 color=red>
521 <p>This shows the geometry for the Drift layers supports.
522 </font>
523 <pre>
524  */
525 //End_Html
526 //
527 //    Read a file containing the geometry for the ITS version 5.
528 ////////////////////////////////////////////////////////////////////////
529
530     Int_t size;
531     char topvol[5];
532     char *filtmp;
533
534   filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
535   FILE *file = fopen(filtmp,"r");
536   delete [] filtmp;
537   if(file) {
538     fclose(file);
539     printf("Ready to read Euclid geometry file\n");
540     ReadEuclid(fEuclidGeometry.Data(),topvol);
541     printf("Read in euclid geometries\n");
542   } else {
543     Error("CreateGeometry"," THE GEOM FILE %s DOES NOT EXIST !",
544           fEuclidGeometry.Data());
545     exit(1);
546   } // end if(file)
547   //
548   // Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
549   // invisible
550   //
551   gMC->Gspos("ITSV",1,"ALIC",0,0,0,0,"ONLY");
552   //
553   // Outputs the geometry tree in the EUCLID/CAD format if requested to do so
554   
555     if (fEuclidOut) {
556       gMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
557     } // end if (fEuclidOut)
558
559     // read in the file containing the transformations for the active
560     // volumes for the ITS version 5. This is expected to be in a file
561     // ending in .det. This geometry is kept in the AliITSgeom class.
562     filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
563     size = strlen(filtmp);
564     if(size>4){
565         filtmp[size-3] = 'd'; // change from .euc to .det
566         filtmp[size-2] = 'e';
567         filtmp[size-1] = 't';
568         file = fopen(filtmp,"r");
569         if(file){ // if file exists use it to fill AliITSgeom structure.
570             fclose(file);
571             printf("ready to read .det file %s\n",filtmp);
572             fITSgeom = new AliITSgeom(filtmp);
573         }else{
574             fITSgeom = 0;
575             // fill AliITSgeom structure from geant structure just filled above
576         }// end if(file)
577         delete [] filtmp;
578     }// end if(size>4)
579     printf("finished with euclid geometrys\n");
580 }
581 //_____________________________________________________________________________
582 void AliITSv5::Init(){
583 ////////////////////////////////////////////////////////////////////////
584 //     Initialise the ITS after it has been created.
585 ////////////////////////////////////////////////////////////////////////
586     Int_t i,j,l;
587
588     fIdN    = fId5N;
589     fIdName = new char*[fId5N];
590     fIdSens = new Int_t[fId5N];
591     for(i=0;i<fId5N;i++) {
592         l = strlen(fId5Name[i]);
593         fIdName[i] = new char[l+1];
594         for(j=0;j<l;j++) fIdName[i][j] = fId5Name[i][j];
595         fIdName[i][l] = '\0'; // Null terminate this string.
596     } // end for i
597
598     AliITS::Init();
599     fMajorVersion = 5;
600     fMinorVersion = 0;
601
602 //_____________________________________________________________________________
603 void AliITSv5::StepManager(){
604 ////////////////////////////////////////////////////////////////////////
605 //    Called for every step in the ITS, then calles the AliITShit class
606 // creator with the information to be recoreded about that hit.
607 //     The value of the macro ALIITSPRINTGEOM if set to 1 will allow the
608 // printing of information to a file which can be used to create a .det
609 // file read in by the routine CreateGeometry(). If set to 0 or any other
610 // value except 1, the default behavior, then no such file is created nor
611 // it the extra variables and the like used in the printing allocated.
612 ////////////////////////////////////////////////////////////////////////
613   Int_t          copy, id;
614   Int_t          copy1,copy2;
615   Float_t        hits[8];
616   Int_t          vol[4];
617   TLorentzVector position, momentum;
618   TClonesArray   &lhits = *fHits;
619 #if ALIITSPRINTGEOM==1
620   FILE          *fp;
621   Int_t         i;
622   Float_t       xl[3],xt[3],angl[6];
623 //  Float_t       par[20],att[20];
624   Float_t      mat[9];
625   static Bool_t first=kTRUE,printit[6][50][50];
626   if(first){ for(copy1=0;copy1<6;copy1++)for(copy2=0;copy2<50;copy2++)
627       for(id=0;id<50;id++) printit[copy1][copy2][id] = kTRUE;
628   first = kFALSE;
629   }
630   // end if first
631 #endif
632   //
633   // Track status
634   vol[3] = 0;
635   if(gMC->IsTrackInside())      vol[3] +=  1;
636   if(gMC->IsTrackEntering())    vol[3] +=  2;
637   if(gMC->IsTrackExiting())     vol[3] +=  4;
638   if(gMC->IsTrackOut())         vol[3] +=  8;
639   if(gMC->IsTrackDisappeared()) vol[3] += 16;
640   if(gMC->IsTrackStop())        vol[3] += 32;
641   if(gMC->IsTrackAlive())       vol[3] += 64;
642   //
643   // Fill hit structure.
644   if(!(gMC->TrackCharge())) return;
645   //
646   // Only entering charged tracks
647   if((id = gMC->CurrentVolID(copy)) == fIdSens[0]) {
648       vol[0] = 1;
649       id = gMC->CurrentVolOffID(0,copy);
650       //detector copy in the ladder = 1<->4  (ITS1)
651       vol[1] = copy;
652       gMC->CurrentVolOffID(1,copy1);
653       //ladder copy in the module   = 1<->2  (I186)
654       gMC->CurrentVolOffID(2,copy2);
655       //module copy in the layer    = 1<->10 (I132)
656       vol[2] = copy1+(copy2-1)*2;//# of ladders in one module  = 2
657   } else if(id == fIdSens[1]){
658       vol[0] = 2;
659       id = gMC->CurrentVolOffID(0,copy);
660       //detector copy in the ladder = 1<->4  (ITS2)
661       vol[1] = copy;
662       gMC->CurrentVolOffID(1,copy1);
663       //ladder copy in the module   = 1<->4  (I131)
664       gMC->CurrentVolOffID(2,copy2);
665       //module copy in the layer    = 1<->10 (I132)
666       vol[2] = copy1+(copy2-1)*4;//# of ladders in one module  = 4
667   } else if(id == fIdSens[2]){
668       vol[0] = 3;
669       id = gMC->CurrentVolOffID(1,copy);
670       //detector copy in the ladder = 1<->5  (ITS3 is inside I314)
671       vol[1] = copy;
672       id = gMC->CurrentVolOffID(2,copy);
673       //ladder copy in the layer    = 1<->12 (I316)
674       vol[2] = copy;
675   } else if(id == fIdSens[3]){
676       vol[0] = 4;
677       id = gMC->CurrentVolOffID(1,copy);
678       //detector copy in the ladder = 1<->8  (ITS4 is inside I414)
679       vol[1] = copy;
680       id = gMC->CurrentVolOffID(2,copy);
681       //ladder copy in the layer    = 1<->22 (I417)
682       vol[2] = copy;
683   }else if(id == fIdSens[4]){
684       vol[0] = 5;
685       id = gMC->CurrentVolOffID(1,copy);
686       //detector copy in the ladder = 1<->23  (ITS5 is inside I562)
687       vol[1] = copy;
688       id = gMC->CurrentVolOffID(2,copy);
689      //ladder copy in the layer    = 1<->34 (I565)
690       vol[2] = copy;
691   }else if(id == fIdSens[5]){
692       vol[0] = 6;
693       id = gMC->CurrentVolOffID(1,copy);
694       //detector copy in the ladder = 1<->26  (ITS6 is inside I566)
695       vol[1] = copy;
696       id = gMC->CurrentVolOffID(2,copy);
697       //ladder copy in the layer = 1<->38 (I569)
698       vol[2] = copy;
699   } else {
700       return; // not an ITS volume?
701   } // end if/else if (gMC->CurentVolID(copy) == fIdSens[i])
702 //
703   gMC->TrackPosition(position);
704   gMC->TrackMomentum(momentum);
705   hits[0]=position[0];
706   hits[1]=position[1];
707   hits[2]=position[2];
708   hits[3]=momentum[0];
709   hits[4]=momentum[1];
710   hits[5]=momentum[2];
711   hits[6]=gMC->Edep();
712   hits[7]=gMC->TrackTime();
713   // Fill hit structure with this new hit.
714   new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
715 #if ALIITSPRINTGEOM==1
716   if(printit[vol[0]][vol[2]][vol[1]]){
717       printit[vol[0]][vol[2]][vol[1]] = kFALSE;
718       xl[0] = xl[1] = xl[2] = 0.0;
719       gMC->Gdtom(xl,xt,1);
720       for(i=0;i<9;i++) mat[i] = 0.0;
721       mat[0] = mat[4] = mat[8] = 1.0;  // default with identity matrix
722       xl[0] = 1.0;
723       xl[1] = xl[2] =0.0;
724       gMC->Gdtom(xl,&(mat[0]),2);
725       xl[1] = 1.0;
726       xl[0] = xl[2] =0.0;
727       gMC->Gdtom(xl,&(mat[3]),2);
728       xl[2] = 1.0;
729       xl[1] = xl[0] =0.0;
730       gMC->Gdtom(xl,&(mat[6]),2);
731
732       angl[0] = TMath::ACos(mat[2]);
733       if(mat[2]==1.0) angl[0] = 0.0;
734       angl[1] = TMath::ATan2(mat[1],mat[0]);
735       if(angl[1]<0.0) angl[1] += 2.0*TMath::Pi();
736
737       angl[2] = TMath::ACos(mat[5]);
738       if(mat[5]==1.0) angl[2] = 0.0;
739       angl[3] = TMath::ATan2(mat[4],mat[3]);
740       if(angl[3]<0.0) angl[3] += 2.0*TMath::Pi();
741
742       angl[4] = TMath::ACos(mat[8]);
743       if(mat[8]==1.0) angl[4] = 0.0;
744       angl[5] = TMath::ATan2(mat[7],mat[6]);
745       if(angl[5]<0.0) angl[5] += 2.0*TMath::Pi();
746
747       for(i=0;i<6;i++) angl[i] *= 180.0/TMath::Pi(); // degrees
748       fp = fopen("ITSgeometry_v5.det","a");
749       fprintf(fp,"%2d %2d %2d %9e %9e %9e %9e %9e %9e %9e %9e %9e ",
750               vol[0],vol[2],vol[1], // layer ladder detector
751               xt[0],xt[1],xt[2],    // Translation vector
752               angl[0],angl[1],angl[2],angl[3],angl[4],angl[5] // Geant rotaion
753                                                            // angles (degrees)
754               );
755       fprintf(fp,"%9e %9e %9e %9e %9e %9e %9e %9e %9e",
756              mat[0],mat[1],mat[2],mat[3],mat[4],mat[5],mat[6],mat[7],mat[8]
757           );  // Adding the rotation matrix.
758       fprintf(fp,"\n");
759       fclose(fp);
760   } // end if printit[layer][ladder][detector]
761 #endif
762   return;
763 }
764 //____________________________________________________________________________
765 void AliITSv5::Streamer(TBuffer &R__b){
766 ////////////////////////////////////////////////////////////////////////
767 //    A dummy Streamer function for this class AliITSv5. By default it
768 // only streams the AliITS class as it is required. Since this class
769 // dosen't contain any "real" data to be saved, it doesn't.
770 ////////////////////////////////////////////////////////////////////////
771
772   //    printf("AliITSv5Streamer Starting\n");
773    if (R__b.IsReading()) {
774       Version_t R__v = R__b.ReadVersion(); 
775       if (R__v==1) {
776           AliITS::Streamer(R__b);
777           // This information does not need to be read. It is "hard wired"
778           // into this class via its creators.
779           //R__b >> fId5N;
780           //R__b.ReadArray(fId5Name);
781       }else{
782       } // end if R__v==1
783    } else {
784       R__b.WriteVersion(AliITSv5::IsA());
785       AliITS::Streamer(R__b);
786       // This information does not need to be saved. It is "hard wired"
787       // into this class via its creators.
788       //R__b << fId5N;
789       //R__b.WriteArray(fId5Name, __COUNTER__);
790    } // end if R__b.IsReading()
791    //    printf("AliITSv5Streamer Finishing\n");
792 }
793