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