]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/AliITSv5asymm.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Access function to local momenta renamed.
[u/mrichter/AliRoot.git] / ITS / AliITSv5asymm.cxx
CommitLineData
d683f3d5 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$
202679a1 18Revision 1.9 2001/02/14 14:12:43 nilsen
19Forgot about the nonexisting v in ITSgeometry_5.det.
20
262d891d 21Revision 1.8 2001/02/14 14:06:23 nilsen
22Wrong default .det file name for old .det file format.
23
88cc8e8a 24Revision 1.7 2001/02/13 16:53:35 nilsen
25Fixed a but when trying to use GEANT4. Needed to replace
26if(!((TGeant3*)gMC)) with if(!(dynamic_casst<TGeant3*>(gMC)))
27because just casting gMC to be TGeant3* even when it realy is a TGeant3 pointer
28did not result in a zero value. For AliITSv5asymm and AliITSv5symm, needed
29to fix a bug in the initilizers and a bug in BuildGeometry. This is now done
30in the same way as in AliITSv5.cxx.
31
9a4c6ea3 32Revision 1.6 2001/02/09 20:06:26 nilsen
33Fixed bug in distructor. Can't distroy fixxed length arrays. Thanks Peter.
34
58e13d99 35Revision 1.5 2001/02/09 00:05:31 nilsen
36Added fMajor/MinorVersion variables and made other changes to better make
37use of the new code changes in AliITSgeom related classes.
38
bae7e562 39Revision 1.4 2001/02/02 23:57:29 nilsen
40Added include file that are no londer included in AliITSgeom.h
41
d953664a 42Revision 1.3 2001/01/30 09:23:13 hristov
43Streamers removed (R.Brun)
44
a8a6107b 45Revision 1.2 2000/11/30 11:13:11 barbera
46 Added changes suggested by Federico Carminati on nov, 30, 2000
47
a07ced09 48Revision 1.1 2000/10/07 15:46:38 barbera
49Version 5 of the geometry with asymmetric services
50
d683f3d5 51Revision 1.25 2000/10/05 20:50:00 nilsen
52Now using root generated streamers.
53
54Revision 1.14.4.12 2000/10/02 16:04:03 barbera
55Forward declarations added
56
57Revision 1.22 2000/07/10 16:07:19 fca
58Release version of ITS code
59
60Revision 1.14.4.4 2000/05/19 10:10:21 nilsen
61fix for bug with HP and Sun unix + fix for event display in ITS-working branch
62
63Revision 1.14.4.3 2000/03/04 23:46:38 nilsen
64Fixed up the comments/documentation.
65
66Revision 1.14.4.2 2000/03/02 21:53:02 nilsen
67To make it compatable with the changes in AliRun/AliModule.
68
69Revision 1.14.4.1 2000/01/12 19:03:33 nilsen
70This is the version of the files after the merging done in December 1999.
71See the ReadMe110100.txt file for details
72
73Revision 1.14 1999/10/22 08:16:49 fca
74Correct destructors, thanks to I.Hrivnacova
75
76Revision 1.13 1999/10/06 10:15:19 fca
77Correct bug in allocation of layer name and add destructor
78
79Revision 1.12 1999/10/05 08:05:09 fca
80Minor corrections for uninitialised variables.
81
82Revision 1.11 1999/09/29 09:24:20 fca
83Introduction of the Copyright and cvs Log
84
85*/
86
87///////////////////////////////////////////////////////////////////////////////
88//
89// Inner Traking System version 5 with asymmetric services
90// This class contains the base procedures for the Inner Tracking System
91//
92// Authors: R. Barbera, B. S. Nilsen.
93// version 5.
94// Created October 7 2000.
95//
96///////////////////////////////////////////////////////////////////////////////
d953664a 97#include <iostream.h>
98#include <iomanip.h>
d683f3d5 99#include <stdio.h>
100#include <stdlib.h>
101#include <TMath.h>
102#include <TGeometry.h>
103#include <TNode.h>
104#include <TTUBE.h>
105#include <TFile.h> // only required for Tracking function?
106#include <TCanvas.h>
107#include <TObjArray.h>
108#include <TObjString.h>
109#include <TClonesArray.h>
a07ced09 110#include <TLorentzVector.h>
d683f3d5 111#include <TBRIK.h>
112#include <TSystem.h>
113
114#include "AliMC.h"
115#include "AliRun.h"
d683f3d5 116#include "../TGeant3/TGeant3.h"
d683f3d5 117#include "AliITShit.h"
bae7e562 118#include "AliITSGeant3Geometry.h"
d683f3d5 119#include "AliITS.h"
120#include "AliITSv5asymm.h"
121#include "AliITSgeom.h"
bae7e562 122#include "AliITSgeomSPD.h"
123#include "AliITSgeomSDD.h"
124#include "AliITSgeomSSD.h"
d683f3d5 125
126ClassImp(AliITSv5asymm)
127
128//_____________________________________________________________________________
129AliITSv5asymm::AliITSv5asymm() {
130////////////////////////////////////////////////////////////////////////
131// Standard default constructor for the ITS version 5.
132////////////////////////////////////////////////////////////////////////
bae7e562 133 Int_t i;
d683f3d5 134
135 fIdN = 0;
136 fIdName = 0;
137 fIdSens = 0;
bae7e562 138 fEuclidOut = kFALSE; // Don't write Euclide file
139 fGeomDetOut = kFALSE; // Don't write .det file
9a4c6ea3 140 fGeomDetIn = kTRUE; // Read .det file
bae7e562 141 fGeomOldDetIn = kTRUE; // Read old formatted .det file
142 fMajorVersion = IsVersion();
143 fMinorVersion = 3;
144 for(i=0;i<60;i++) fRead[i] = '\0';
145 for(i=0;i<60;i++) fWrite[i] = '\0';
146 for(i=0;i<60;i++) fEuclidGeomDet[i] = '\0';
147}
148//_____________________________________________________________________________
149AliITSv5asymm::AliITSv5asymm(const char *name, const char *title) : AliITS(name, title){
150/////////////////////////////////////////////////////////////////////////////
151// Standard constructor for the ITS version 5 with symmetrical services.
152/////////////////////////////////////////////////////////////////////////////
153 Int_t i;
154
155 fIdN = 6;
156 fIdName = new TString[fIdN];
157 fIdName[0] = "ITS1";
158 fIdName[1] = "ITS2";
159 fIdName[2] = "ITS3";
160 fIdName[3] = "ITS4";
161 fIdName[4] = "ITS5";
162 fIdName[5] = "ITS6";
163 fIdSens = new Int_t[fIdN];
164 for (i=0;i<fIdN;i++) fIdSens[i] = 0;
165 fEuclidOut = kFALSE; // Don't write Euclide file
166 fGeomDetOut = kFALSE; // Don't write .det file
9a4c6ea3 167 fGeomDetIn = kTRUE; // Read .det file
bae7e562 168 fGeomOldDetIn = kTRUE; // Read old formatted .det file
169 fMajorVersion = IsVersion();
170 fMinorVersion = 3;
171 for(i=0;i<60;i++) fRead[i] = '\0';
172 for(i=0;i<60;i++) fWrite[i] = '\0';
173
174 fEuclidMaterial = "$ALICE_ROOT/Euclid/ITSgeometry_5asymm.tme";
175 fEuclidGeometry = "$ALICE_ROOT/Euclid/ITSgeometry_5asymm.euc";
202679a1 176 strncpy(fEuclidGeomDet,"$ALICE_ROOT/Euclid/ITSgeometry_5.det",60);
bae7e562 177 strncpy(fRead,fEuclidGeomDet,60);
178 strncpy(fWrite,fEuclidGeomDet,60);
d683f3d5 179}
180//____________________________________________________________________________
181AliITSv5asymm::AliITSv5asymm(const AliITSv5asymm &source){
182////////////////////////////////////////////////////////////////////////
183// Copy Constructor for ITS version 5.
184////////////////////////////////////////////////////////////////////////
185 if(&source == this) return;
bae7e562 186 Warning("Copy Constructor","Not allowed to copy AliITSv5asymm");
d683f3d5 187 return;
188}
189//_____________________________________________________________________________
190AliITSv5asymm& AliITSv5asymm::operator=(const AliITSv5asymm &source){
191////////////////////////////////////////////////////////////////////////
bae7e562 192// Assignment operator for the ITS version 5.
d683f3d5 193////////////////////////////////////////////////////////////////////////
194
195 if(&source == this) return *this;
bae7e562 196 Warning("= operator","Not allowed to copy AliITSv5asymm");
d683f3d5 197 return *this;
198
199}
200//_____________________________________________________________________________
201AliITSv5asymm::~AliITSv5asymm() {
202////////////////////////////////////////////////////////////////////////
203// Standard destructor for the ITS version 5.
204////////////////////////////////////////////////////////////////////////
d683f3d5 205}
206//_____________________________________________________________________________
207void AliITSv5asymm::BuildGeometry(){
208////////////////////////////////////////////////////////////////////////
209// Geometry builder for the ITS version 5.
210////////////////////////////////////////////////////////////////////////
211 //
212 // Build ITS TNODE geometry for event display using detailed geometry.
213 // This function builds a simple ITS geometry used by the ROOT macro
214 // ITSdisplay.C.
215
216 TNode *top;
217 TNode *nd;
218 //const int kColorITSSPD=kRed;
219 //const int kColorITSSDD=kGreen;
220 const int kColorITSSSD=kBlue;
221 //
d683f3d5 222 AliITSgeom *gm = this->GetITSgeom();
9a4c6ea3 223 if(gm==0) return;
224 top=gAlice->GetGeometry()->GetNode("alice");
d683f3d5 225
226 Int_t lay,lad,det,i;
227 Text_t name[10];
228 Float_t xg[3];
229 Float_t rt[9];
230 Double_t rtd[9];
231 TBRIK *box;
232 TRotMatrix *rm;
233 //TCanvas *c1 = new TCanvas("c1","ITS");
234
235 for(lay=1;lay<=2;lay++)
236 for(lad=1;lad<=gm->GetNladders(lay);lad++)
237 for(det=1;det<=gm->GetNdetectors(lay);det++){
238 try {
239 box = new TBRIK ("ActiveSPD","Active volume of SPD","SPD SI DET",
240 0.64,0.0075,4.19);
241 } catch (...) {
242 cout << "EXCEPTION in box = new TBRIK" << endl;
243 return;
244 }
245 gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
246 gm->GetRotMatrix(lay,lad,det,rt);
247 //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
248 for(i=0;i<9;i++) rtd[i] = rt[i];
249 try {
250 rm = new TRotMatrix(name,name,rtd);
251 } catch (...) {
252 cout << "EXCEPTION in new TRotMatrix" << endl;
253 return;
254 }
255 top->cd();
256 //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
257 try {
258 nd = new TNode("SPD"," ",box,xg[0],xg[1],xg[2],rm);
259 } catch (...) {
260 cout << "EXCEPTION in new TNode" << endl;
261 return;
262 }
263 nd->SetLineColor(kColorITSSSD);
264 fNodes->Add(nd);
265 }
266
267 for(lay=3;lay<=3;lay++)
268 for(lad=1;lad<=gm->GetNladders(lay);lad++)
269 for(det=1;det<=gm->GetNdetectors(lay);det++){
270 try {
271 box = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
272 3.5,0.014,3.763);
273 } catch (...) {
274 cout << "EXCEPTION in box = new TBRIK" << endl;
275 return;
276 }
277 gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
278 gm->GetRotMatrix(lay,lad,det,rt);
279 //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
280 for(i=0;i<9;i++) rtd[i] = rt[i];
281 try {
282 rm = new TRotMatrix(name,name,rtd);
283 } catch (...) {
284 cout << "EXCEPTION in new TRotMatrix" << endl;
285 return;
286 }
287 top->cd();
288 //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
289 try {
290 nd = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
291 } catch (...) {
292 cout << "EXCEPTION in new TNode" << endl;
293 return;
294 }
295 nd->SetLineColor(kColorITSSSD);
296 fNodes->Add(nd);
297 }
298
299 for(lay=4;lay<=4;lay++)
300 for(lad=1;lad<=gm->GetNladders(lay);lad++)
301 for(det=1;det<=gm->GetNdetectors(lay);det++){
302 try {
303 box = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
304 3.5,0.014,3.763);
305 } catch (...) {
306 cout << "EXCEPTION in box = new TBRIK" << endl;
307 return;
308 }
309 gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
310 gm->GetRotMatrix(lay,lad,det,rt);
311 //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
312 for(i=0;i<9;i++) rtd[i] = rt[i];
313 try {
314 rm = new TRotMatrix(name,name,rtd);
315 } catch (...) {
316 cout << "EXCEPTION in new TRotMatrix" << endl;
317 return;
318 }
319 top->cd();
320 //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
321 try {
322 nd = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
323 } catch (...) {
324 cout << "EXCEPTION in new TNode" << endl;
325 return;
326 }
327 nd->SetLineColor(kColorITSSSD);
328 fNodes->Add(nd);
329 }
330 for(lay=5;lay<=5;lay++)
331 for(lad=1;lad<=gm->GetNladders(lay);lad++)
332 for(det=1;det<=gm->GetNdetectors(lay);det++){
333 try {
334 box = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
335 3.65,0.015,2.0);
336 } catch (...) {
337 cout << "EXCEPTION in box = new TBRIK" << endl;
338 return;
339 }
340 gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
341 gm->GetRotMatrix(lay,lad,det,rt);
342 //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
343 for(i=0;i<9;i++) rtd[i] = rt[i];
344 try {
345 rm = new TRotMatrix(name,name,rtd);
346 } catch (...) {
347 cout << "EXCEPTION in new TRotMatrix" << endl;
348 return;
349 }
350 top->cd();
351 //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
352 try {
353 nd = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
354 } catch (...) {
355 cout << "EXCEPTION in new TNode" << endl;
356 return;
357 }
358 nd->SetLineColor(kColorITSSSD);
359 fNodes->Add(nd);
360 }
361
362 for(lay=6;lay<=6;lay++)
363 for(lad=1;lad<=gm->GetNladders(lay);lad++)
364 for(det=1;det<=gm->GetNdetectors(lay);det++){
365 try {
366 box = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
367 3.65,0.015,2.0);
368 } catch (...) {
369 cout << "EXCEPTION in box = new TBRIK" << endl;
370 return;
371 }
372
373 gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
374 gm->GetRotMatrix(lay,lad,det,rt);
375 //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
376 for(i=0;i<9;i++) rtd[i] = rt[i];
377 try {
378 rm = new TRotMatrix(name,name,rtd);
379 } catch (...) {
380 cout << "EXCEPTION in new TRotMatrix" << endl;
381 return;
382 }
383 top->cd();
384 //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
385 try {
386 nd = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
387 } catch (...) {
388 cout << "EXCEPTION in new TNode" << endl;
389 return;
390 }
391 nd->SetLineColor(kColorITSSSD);
392 fNodes->Add(nd);
393 }
394
395
396}
397//_____________________________________________________________________________
398void AliITSv5asymm::CreateMaterials(){
399////////////////////////////////////////////////////////////////////////
400// Read a file containing the materials for the ITS version 5.
401////////////////////////////////////////////////////////////////////////
402 char *filtmp;
403
404 filtmp = gSystem->ExpandPathName(fEuclidMaterial.Data());
405
406 FILE *file = fopen(filtmp,"r");
407 if(file) {
408 fclose(file);
409 ReadEuclidMedia(filtmp);
410 } else {
411 Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",filtmp);
412 exit(1);
413 } // end if(file)
414}
415//_____________________________________________________________________________
416void AliITSv5asymm::CreateGeometry(){
417//////////////////////////////////////////////////////////////////////
418// This is the geometry used for the ITS Pre-TDR and comes from an
419// Euclid to Geant conversion. The only difference
420// is in the details of the ITS supports. The detectors elements,
421// detector numbering, and local and global reference frames are shown in
422// the following figures.
423//Begin_Html
424/*
425<img src="picts/ITS/its1+2_convention_front_5.gif">
426</pre>
427<br clear=left>
428<font size=+2 color=red>
429<p>This shows the front view of the SPDs.
430</font>
431<pre>
432<img src="picts/ITS/its1+2_convention_side_5.gif">
433</pre>
434<br clear=left>
435<font size=+2 color=red>
436<p>This shows the perspective view of the SPDs.
437</font>
438<img src="picts/ITS/its1+2_tree.gif">
439</pre>
440<br clear=left>
441<font size=+2 color=red>
442<p>This shows the geometry Tree for the SPDs.
443</font>
444<pre>
445
446<pre>
447<img src="picts/ITS/its3+4_convention_front_5.gif">
448</pre>
449<br clear=left>
450<font size=+2 color=red>
451<p>This shows the front view of the SDDs.
452</font>
453<pre>
454<img src="picts/ITS/its3+4_convention_side_5.gif">
455</pre>
456<br clear=left>
457<font size=+2 color=red>
458<p>This shows the perspective view of the SDDs.
459</font>
460<img src="picts/ITS/its3+4_tree.gif">
461</pre>
462<br clear=left>
463<font size=+2 color=red>
464<p>This shows the geometry Tree for the SDDs.
465</font>
466<pre>
467
468<pre>
469<img src="picts/ITS/its5+6_convention_front_5.gif">
470</pre>
471<br clear=left>
472<font size=+2 color=red>
473<p>This shows the front view of the SSDs.
474</font>
475<pre>
476<img src="picts/ITS/its5+6_convention_side_5.gif">
477</pre>
478<br clear=left>
479<font size=+2 color=red>
480<p>This shows the perspective view of the SSDs.
481</font>
482<pre>
483<img src="picts/ITS/its5+6_tree.gif">
484</pre>
485<br clear=left>
486<font size=+2 color=red>
487<p>This shows the geometry Tree for the SSDs.
488</font>
489<pre>
490
491
492<img src="picts/ITS/its_layer1-6_2.gif">
493</pre>
494<br clear=left>
495<font size=+2 color=red>
496<p>This shows the front view of the whole ITS..
497</font>
498<pre>
499
500<img src="picts/ITS/its_layer1-6_1.gif">
501</pre>
502<br clear=left>
503<font size=+2 color=red>
504<p>This shows the perspective view of the whole ITS..
505</font>
506<pre>
507
508<img src="picts/ITS/its1-6_tree.gif">
509</pre>
510<br clear=left>
511<font size=+2 color=red>
512<p>This shows the geometry Tree for the whole ITS.
513</font>
514<pre>
515*/
516//End_Html
517//
518//
519// Here are shown the details of the ITS support cones and services.
520// First is a GEANT tree showing the organization of all of the volumes
521// that make up the ITS supports and services.
522//Begin_Html
523/*
524<img src="picts/ITS/supports_tree.gif">
525 */
526//End_Html
527// What follows are a number of figures showing what these support
528// structures look like.
529//Begin_Html
530/*
531
532<img src="picts/ITS/supports_3.gif">
533</pre>
534<br clear=left>
535<font size=+2 color=red>
536<p>This shows the geometry of the supports for the Drift and Strip layers only.
537</font>
538<pre>
539
540<img src="picts/ITS/supports_2.gif">
541</pre>
542<br clear=left>
543<font size=+2 color=red>
544<p>This shows the geometry of the supports for the Drift and Strip layers in front cut out.
545</font>
546<pre>
547
548<img src="picts/ITS/supports_1.gif">
549</pre>
550<br clear=left>
551<font size=+2 color=red>
552<p>This shows the geometry of the supports for the Drift and Strip layers in a back cut out..
553</font>
554<pre>
555
556<img src="picts/ITS/suppssd.gif">
557</pre>
558<br clear=left>
559<font size=+2 color=red>
560<p>This shows the geometry for the Strip layers supports.
561</font>
562<pre>
563
564<img src="picts/ITS/suppsdd.gif">
565</pre>
566<br clear=left>
567<font size=+2 color=red>
568<p>This shows the geometry for the Drift layers supports.
569</font>
570<pre>
571 */
572//End_Html
573//
574// Read a file containing the geometry for the ITS version 5.
575////////////////////////////////////////////////////////////////////////
d683f3d5 576 char topvol[5];
577 char *filtmp;
578
579 filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
580 FILE *file = fopen(filtmp,"r");
581 delete [] filtmp;
582 if(file) {
583 fclose(file);
9a4c6ea3 584 cout << "Ready to read Euclid geometry file" << endl;
d683f3d5 585 ReadEuclid(fEuclidGeometry.Data(),topvol);
9a4c6ea3 586 cout << "Read in euclid geometries" << endl;
d683f3d5 587 } else {
588 Error("CreateGeometry"," THE GEOM FILE %s DOES NOT EXIST !",
589 fEuclidGeometry.Data());
590 exit(1);
591 } // end if(file)
592 //
593 // Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
594 // invisible
595 //
596 gMC->Gspos("ITSV",1,"ALIC",0,0,0,0,"ONLY");
597 //
598 // Outputs the geometry tree in the EUCLID/CAD format if requested to do so
599
600 if (fEuclidOut) {
601 gMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
602 } // end if (fEuclidOut)
603
9a4c6ea3 604 cout << "finished with euclid geometrys" << endl;
d683f3d5 605}
bae7e562 606
607//______________________________________________________________________
608void AliITSv5asymm::ReadOldGeometry(const char *filename){
609 // read in the file containing the transformations for the active
610 // volumes for the ITS version 5. This is expected to be in a file
611 // ending in .det. This geometry is kept in the AliITSgeom class.
612 Int_t size;
613 char *filtmp;
614 FILE *file;
615
616 filtmp = gSystem->ExpandPathName(filename);
617 size = strlen(filtmp);
618 if(size>4 && fGeomDetIn){
619 filtmp[size-3] = 'd'; // change from .euc to .det
620 filtmp[size-2] = 'e';
621 filtmp[size-1] = 't';
622 file = fopen(filtmp,"r");
623 if(file){ // if file exists use it to fill AliITSgeom structure.
624 fclose(file);
625 fITSgeom = new AliITSgeom(filtmp);
626 fITSgeom->DefineShapes(3); // if fShape isn't defined define it.
627 // Now define the detector types/shapes.
628 fITSgeom->ReSetShape(kSPD,(TObject *) new AliITSgeomSPD300());
629 fITSgeom->ReSetShape(kSDD,(TObject *) new AliITSgeomSDD300());
630 fITSgeom->ReSetShape(kSSD,(TObject *) new AliITSgeomSSD());
631 }else{
632 fITSgeom = 0;
633 // fill AliITSgeom structure from geant structure just filled above
634 }// end if(file)
635 delete [] filtmp;
636 }// end if(size>4)
637}
638//______________________________________________________________________
639void AliITSv5asymm::InitAliITSgeom(){
640// Based on the geometry tree defined in Geant 3.21, this
641// routine initilizes the Class AliITSgeom from the Geant 3.21 ITS geometry
642// sturture.
9a4c6ea3 643 if(!(dynamic_cast<TGeant3*>(gMC))) {
bae7e562 644 Error("InitAliITSgeom",
645 "Wrong Monte Carlo. InitAliITSgeom uses TGeant3 calls");
646 return;
647 } // end if
648 cout << "Reading Geometry transformation directly from Geant 3." << endl;
649 const Int_t nlayers = 6;
650 const Int_t ndeep = 7;
651 Int_t itsGeomTreeNames[nlayers][ndeep],lnam[20],lnum[20];
652 Int_t nlad[nlayers],ndet[nlayers];
653 Double_t t[3],r[10];
654 Float_t par[20],att[20];
655 Int_t npar,natt,idshape,imat,imed;
656 AliITSGeant3Geometry *ig = new AliITSGeant3Geometry();
657 Int_t mod,lay,lad,det,i,j,k;
658 char *names[nlayers][ndeep] = {
659 {"ALIC","ITSV","ITSD","IT12","I132","I186","ITS1"}, // lay=1
660 {"ALIC","ITSV","ITSD","IT12","I132","I131","ITS2"}, // lay=2
661 {"ALIC","ITSV","ITSD","IT34","I004","I302","ITS3"}, // lay=3
662 {"ALIC","ITSV","ITSD","IT34","I005","I402","ITS4"}, // lay=4
663 {"ALIC","ITSV","ITSD","IT56","I565","I562","ITS5"}, // lay=5
664 {"ALIC","ITSV","ITSD","IT56","I569","I566","ITS6"}};// lay=6
665 Int_t itsGeomTreeCopys[nlayers][ndeep] = {{1,1,1,1,10, 2,4}, // lay=1
666 {1,1,1,1,10, 4,4}, // lay=2
667 {1,1,1,1,14, 6,1}, // lay=3
668 {1,1,1,1,22, 8,1}, // lay=4
669 {1,1,1,1,34,23,1}, // lay=5
670 {1,1,1,1,38,26,1}};// lay=6
671
672 // Sorry, but this is not very pritty code. It should be replaced
673 // at some point with a version that can search through the geometry
674 // tree its self.
675 cout << "Reading Geometry informaton from Geant3 common blocks" << endl;
676 for(i=0;i<20;i++) lnam[i] = lnum[i] = 0;
677 for(i=0;i<nlayers;i++)for(j=0;j<ndeep;j++)
678 itsGeomTreeNames[i][j] = ig->StringToInt(names[i][j]);
679 mod = 0;
680 for(i=0;i<nlayers;i++){
681 k = 1;
682 for(j=0;j<ndeep;j++) if(itsGeomTreeCopys[i][j]!=0)
683 k *= TMath::Abs(itsGeomTreeCopys[i][j]);
684 mod += k;
685 } // end for i
686
687 if(fITSgeom!=0) delete fITSgeom;
688 nlad[0]=20;nlad[1]=40;nlad[2]=14;nlad[3]=22;nlad[4]=34;nlad[5]=38;
689 ndet[0]=4;ndet[1]=4;ndet[2]=6;ndet[3]=8;ndet[4]=22;ndet[5]=25;
690 fITSgeom = new AliITSgeom(0,6,nlad,ndet,mod);
691 mod = -1;
692 for(lay=1;lay<=nlayers;lay++){
693 for(j=0;j<ndeep;j++) lnam[j] = itsGeomTreeNames[lay-1][j];
694 for(j=0;j<ndeep;j++) lnum[j] = itsGeomTreeCopys[lay-1][j];
695 switch (lay){
696 case 1: case 2: // layers 1 and 2 are a bit special
697 lad = 0;
698 for(j=1;j<=itsGeomTreeCopys[lay-1][4];j++){
699 lnum[4] = j;
700 for(k=1;k<=itsGeomTreeCopys[lay-1][5];k++){
701 lad++;
702 lnum[5] = k;
703 for(det=1;det<=itsGeomTreeCopys[lay-1][6];det++){
704 lnum[6] = det;
705 mod++;
706 ig->GetGeometry(ndeep,lnam,lnum,t,r,idshape,npar,natt,
707 par,att,imat,imed);
708 fITSgeom->CreatMatrix(mod,lay,lad,det,kSPD,t,r);
709 if(!(fITSgeom->IsShapeDefined((Int_t)kSPD)))
710 if(fMinorVersion==1){
711 fITSgeom->ReSetShape(kSPD,
712 new AliITSgeomSPD300());
713 } else if(fMinorVersion==2)
714 fITSgeom->ReSetShape(kSPD,
715 new AliITSgeomSPD300());
716 } // end for det
717 } // end for k
718 } // end for j
719 break;
720 case 3: case 4: case 5: case 6: // layers 3-6
721 lnum[6] = 1;
722 for(lad=1;lad<=itsGeomTreeCopys[lay-1][4];lad++){
723 lnum[4] = lad;
724 for(det=1;det<=itsGeomTreeCopys[lay-1][5];det++){
725 lnum[5] = det;
726 mod++;
727 ig->GetGeometry(7,lnam,lnum,t,r,idshape,npar,natt,
728 par,att,imat,imed);
729 switch (lay){
730 case 3: case 4:
731 fITSgeom->CreatMatrix(mod,lay,lad,det,kSDD,t,r);
732 if(!(fITSgeom->IsShapeDefined(kSDD)))
733 fITSgeom->ReSetShape(kSDD,new AliITSgeomSDD300());
734 break;
735 case 5: case 6:
736 fITSgeom->CreatMatrix(mod,lay,lad,det,kSSD,t,r);
737 if(!(fITSgeom->IsShapeDefined(kSSD)))
738 fITSgeom->ReSetShape(kSSD,new AliITSgeomSSD175());
739 break;
740 } // end switch
741 } // end for det
742 } // end for lad
743 break;
744 } // end switch
745 } // end for lay
746 return;
747}
d683f3d5 748//_____________________________________________________________________________
749void AliITSv5asymm::Init(){
750////////////////////////////////////////////////////////////////////////
751// Initialise the ITS after it has been created.
752////////////////////////////////////////////////////////////////////////
bae7e562 753 Int_t i;
9a4c6ea3 754 Bool_t bg = kFALSE;
d683f3d5 755
bae7e562 756 cout << endl;
757 for(i=0;i<30;i++) cout << "*";cout << " ITSv5_Init ";
758 for(i=0;i<30;i++) cout << "*";cout << endl;
759//
9a4c6ea3 760 if(fITSgeom==0) bg = kTRUE;
bae7e562 761 if(fRead[0]=='\0') strncpy(fRead,fEuclidGeomDet,60);
762 if(fWrite[0]=='\0') strncpy(fWrite,fEuclidGeomDet,60);
763 if(fITSgeom!=0) delete fITSgeom;
764 fITSgeom = new AliITSgeom();
765 if(fGeomDetIn && !fGeomOldDetIn) fITSgeom->ReadNewFile(fRead);
766 if(fGeomDetIn && fGeomOldDetIn) this->ReadOldGeometry(fRead);
767
768 if(!fGeomDetIn) this->InitAliITSgeom();
769 if(fGeomDetOut) fITSgeom->WriteNewFile(fWrite);
9a4c6ea3 770 if(bg) BuildGeometry(); // call BuildGeometry here if fITSgeom was not
771 // defined ealier.
d683f3d5 772 AliITS::Init();
bae7e562 773//
774 for(i=0;i<72;i++) cout << "*";
775 cout << endl;
d683f3d5 776}
777//_____________________________________________________________________________
778void AliITSv5asymm::StepManager(){
779////////////////////////////////////////////////////////////////////////
780// Called for every step in the ITS, then calles the AliITShit class
781// creator with the information to be recoreded about that hit.
782// The value of the macro ALIITSPRINTGEOM if set to 1 will allow the
783// printing of information to a file which can be used to create a .det
784// file read in by the routine CreateGeometry(). If set to 0 or any other
785// value except 1, the default behavior, then no such file is created nor
786// it the extra variables and the like used in the printing allocated.
787////////////////////////////////////////////////////////////////////////
788 Int_t copy, id;
789 Int_t copy1,copy2;
790 Float_t hits[8];
791 Int_t vol[4];
792 TLorentzVector position, momentum;
793 TClonesArray &lhits = *fHits;
d683f3d5 794 //
795 // Track status
796 vol[3] = 0;
797 if(gMC->IsTrackInside()) vol[3] += 1;
798 if(gMC->IsTrackEntering()) vol[3] += 2;
799 if(gMC->IsTrackExiting()) vol[3] += 4;
800 if(gMC->IsTrackOut()) vol[3] += 8;
801 if(gMC->IsTrackDisappeared()) vol[3] += 16;
802 if(gMC->IsTrackStop()) vol[3] += 32;
803 if(gMC->IsTrackAlive()) vol[3] += 64;
804 //
805 // Fill hit structure.
806 if(!(gMC->TrackCharge())) return;
807 //
808 // Only entering charged tracks
809 if((id = gMC->CurrentVolID(copy)) == fIdSens[0]) {
810 vol[0] = 1;
811 id = gMC->CurrentVolOffID(0,copy);
812 //detector copy in the ladder = 1<->4 (ITS1)
813 vol[1] = copy;
814 gMC->CurrentVolOffID(1,copy1);
815 //ladder copy in the module = 1<->2 (I186)
816 gMC->CurrentVolOffID(2,copy2);
817 //module copy in the layer = 1<->10 (I132)
818 vol[2] = copy1+(copy2-1)*2;//# of ladders in one module = 2
819 } else if(id == fIdSens[1]){
820 vol[0] = 2;
821 id = gMC->CurrentVolOffID(0,copy);
822 //detector copy in the ladder = 1<->4 (ITS2)
823 vol[1] = copy;
824 gMC->CurrentVolOffID(1,copy1);
825 //ladder copy in the module = 1<->4 (I131)
826 gMC->CurrentVolOffID(2,copy2);
827 //module copy in the layer = 1<->10 (I132)
828 vol[2] = copy1+(copy2-1)*4;//# of ladders in one module = 4
829 } else if(id == fIdSens[2]){
830 vol[0] = 3;
831 id = gMC->CurrentVolOffID(1,copy);
832 //detector copy in the ladder = 1<->5 (ITS3 is inside I314)
833 vol[1] = copy;
834 id = gMC->CurrentVolOffID(2,copy);
835 //ladder copy in the layer = 1<->12 (I316)
836 vol[2] = copy;
837 } else if(id == fIdSens[3]){
838 vol[0] = 4;
839 id = gMC->CurrentVolOffID(1,copy);
840 //detector copy in the ladder = 1<->8 (ITS4 is inside I414)
841 vol[1] = copy;
842 id = gMC->CurrentVolOffID(2,copy);
843 //ladder copy in the layer = 1<->22 (I417)
844 vol[2] = copy;
845 }else if(id == fIdSens[4]){
846 vol[0] = 5;
847 id = gMC->CurrentVolOffID(1,copy);
848 //detector copy in the ladder = 1<->23 (ITS5 is inside I562)
849 vol[1] = copy;
850 id = gMC->CurrentVolOffID(2,copy);
851 //ladder copy in the layer = 1<->34 (I565)
852 vol[2] = copy;
853 }else if(id == fIdSens[5]){
854 vol[0] = 6;
855 id = gMC->CurrentVolOffID(1,copy);
856 //detector copy in the ladder = 1<->26 (ITS6 is inside I566)
857 vol[1] = copy;
858 id = gMC->CurrentVolOffID(2,copy);
859 //ladder copy in the layer = 1<->38 (I569)
860 vol[2] = copy;
861 } else {
862 return; // not an ITS volume?
863 } // end if/else if (gMC->CurentVolID(copy) == fIdSens[i])
864//
865 gMC->TrackPosition(position);
866 gMC->TrackMomentum(momentum);
867 hits[0]=position[0];
868 hits[1]=position[1];
869 hits[2]=position[2];
870 hits[3]=momentum[0];
871 hits[4]=momentum[1];
872 hits[5]=momentum[2];
873 hits[6]=gMC->Edep();
874 hits[7]=gMC->TrackTime();
875 // Fill hit structure with this new hit.
876 new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
d683f3d5 877 return;
878}