Coding conventions + removal of redundant methods
[u/mrichter/AliRoot.git] / ITS / AliITSv11.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 //  Inner Traking System version 11                                         //
21 //  This class contains the base procedures for the Inner Tracking System   //
22 //                                                                          //
23 // Authors: R. Barbera                                                      //
24 // version 6.                                                               //
25 // Created  2000.                                                           //
26 //                                                                          //
27 //  NOTE: THIS IS THE  SYMMETRIC PPR geometry of the ITS.                   //
28 // THIS WILL NOT WORK                                                       //
29 // with the geometry or module classes or any analysis classes. You are     //
30 // strongly encouraged to uses AliITSv5.                                    //
31 //                                                                          //
32 //////////////////////////////////////////////////////////////////////////////
33 // See AliITSv11::StepManager().
34 // General C/C++ includes
35 #include <stdio.h>
36 #include <stdlib.h>
37 // General Root includes
38 #include <Riostream.h>
39 #include <TMath.h>
40 #include <float.h>
41 #include <TObjArray.h>
42 #include <TClonesArray.h>
43 #include <TLorentzVector.h>
44 #include <TObjString.h>
45 // Root Geometry includes
46 #include <TGeoManager.h>
47 #include <TGeoVolume.h>
48 #include <TGeoPcon.h>
49 #include <TGeoCone.h>
50 #include <TGeoTube.h> // contaings TGeoTubeSeg
51 #include <TGeoArb8.h>
52 #include <TGeoCompositeShape.h>
53 #include <TGeoMatrix.h>
54 #include <TGeoNode.h>
55 #include <TGeoMaterial.h>
56 #include <TGeoMedium.h>
57 #include "AliITSv11GeometrySupport.h"
58 // General AliRoot includes
59 #include "AliRun.h"
60 #include "AliMagF.h"
61 #include "AliConst.h"
62 // ITS specific includes
63 #include "AliITSgeomSPD.h"
64 #include "AliITSgeomSDD.h"
65 #include "AliITSDetType.h"
66 #include "AliITSresponseSPD.h"
67 #include "AliITSresponseSDD.h"
68 #include "AliITSsegmentationSPD.h"
69 #include "AliITSsegmentationSDD.h"
70 #include "AliITSsimulationSPD.h"
71 #include "AliITSsimulationSDD.h"
72 #include "AliITSClusterFinderSSD.h"
73 #include "AliITSv11.h"
74 #include "AliITSv11GeometrySupport.h"
75
76 ClassImp(AliITSv11)
77
78 /*
79   Some temparary #define's used untill ROOT has addoppted the proper
80   Getter in it's classes.
81   These Below are for TGeoPcon functions.
82 */
83
84 //______________________________________________________________________
85 AliITSv11::AliITSv11() : AliITS(),
86 fGeomDetOut(kFALSE),
87 fGeomDetIn(kFALSE),
88 fMajorVersion(11),
89 fMinorVersion(0),
90 fDet1(0.0),
91 fDet2(0.0),
92 fChip1(0.0),
93 fChip2(0.0),
94 fRails(0),
95 fFluid(1){
96     // Standard default constructor for the ITS version 11.
97     // Inputs:
98     //   none.
99     // Outputs:
100     //   none.
101     // Return
102     //   A default constructed AliITSv11 class.
103
104     //fITSV = 0;
105     //fcS = 0;
106 //   fcD = 0;
107
108    SetEUCLID(kFALSE);
109 }
110 //______________________________________________________________________
111 AliITSv11::AliITSv11(const char *title) : AliITS("ITS", title),
112 fGeomDetOut(kFALSE),
113 fGeomDetIn(kFALSE),
114 fMajorVersion(11),
115 fMinorVersion(0),
116 fDet1(0.0),
117 fDet2(0.0),
118 fChip1(0.0),
119 fChip2(0.0),
120 fRails(0),
121 fFluid(1){
122     // Standard constructor for the ITS version 11.
123     // Inputs:
124     //   const char *title  The title of for this geometry.
125     // Outputs:
126     //   none.
127     // Return
128     //   A Standard constructed AliITSv11 class.
129    SetEUCLID(kFALSE);
130 }
131 //______________________________________________________________________
132 AliITSv11::~AliITSv11() {
133     // Standard destructor for the ITS version 11.
134     // Inputs:
135     //   none.
136     // Outputs:
137     //   none.
138     // Return
139     //   none.
140
141 //    if(fITSV!=0) delete fITSV;
142 //    if(fcS!=0) delete fcS;
143 //    if(fcD!=0) delete fcD;
144 }
145 //______________________________________________________________________
146 AliITSv11::AliITSv11(const AliITSv11 &source) : AliITS(source){
147     //     Copy Constructor for ITS version 11.
148     // Inputs:
149     //   AliITSv11 &source  class to be copied from.
150     // Outputs:
151     //   none.
152     // Return
153     //   none.
154
155     if(&source == this) return;
156     Error("Copy Constructor","Not allowed to copy AliITSv11");
157     return;
158 }
159 //______________________________________________________________________
160 AliITSv11& AliITSv11::operator=(const AliITSv11 &source){
161     //    Assignment operator for the ITS version 11.
162     // Inputs:
163     //   AliITSv11 &source  class to be copied from.
164     // Outputs:
165     //   none.
166     // Return
167     //   none.
168
169     if(&source == this) return *this;
170     Error("= operator","Not allowed to copy AliITSv11");
171     return *this;
172 }
173 //______________________________________________________________________
174 void AliITSv11::BuildGeometry(){
175     // This routine defines and Creates the geometry for version 11 of 
176     // the ITS for use in the simulation display routines. This is a 
177     // very simplified geometry for speed of viewing.
178     // Inputs:
179     //   none.
180     // Outputs:
181     //   none.
182     // Return
183     //   none.
184     //TVector3 t(0.0,0.0,0.0);
185
186     //if(fITSV==0) fITSV = new AliITSGeometryITSV(this,"ALIC");
187     //if(fcS==0) fcS = new AliITSGeometrySSDCone(this,t,"TSV",1);
188
189     //fcS->BuildDisplayGeometry();
190 }
191 //______________________________________________________________________
192 void AliITSv11::CreateGeometry(){
193     // This routine defines and Creates the geometry for version 11 of 
194     // the ITS. The geometry is used by the particle trasport routines,
195     // and therefore, is very detailed.
196     // Inputs:
197     //   none.
198     // Outputs:
199     //   none.
200     // Return
201     //   none.
202     TVector3 t(0.0,0.0,0.0);
203     const Double_t kcm = 1.0;
204
205     TGeoManager *mgr = gGeoManager;
206     TGeoVolume *vALIC = mgr->GetTopVolume();
207
208     TGeoPcon *sITS = new TGeoPcon("ITS Top Volume, Daughter of ALIC",
209                                   0.0,360.0,2);
210     // DefineSection(section number, Z, Rmin, Rmax).
211     sITS->DefineSection(0,-300.0*kcm,0.01*kcm,50.0*kcm);
212     sITS->DefineSection(1,+300.0*kcm,0.01*kcm,50.0*kcm);
213     TGeoVolume *vITS = new TGeoVolume("ITSV",sITS,0);
214     mgr->AddVolume(vITS);
215     vITS->SetVisibility(kFALSE);
216     vALIC->AddNode(vITS,1,0);
217     //
218     AliITSv11GeometrySupport *sup = new AliITSv11GeometrySupport(GetDebug());
219     sup->SPDCone(vITS);
220     sup->SPDThermalSheald(vITS);
221     sup->SDDCone(vITS);
222     sup->SSDCone(vITS);
223     sup->ServicesCableSupport(vITS);
224 }
225 //______________________________________________________________________
226 void AliITSv11::CreateMaterials(){
227     // Create ITS materials
228     //     This function defines the default materials used in the Geant
229     // Monte Carlo simulations for the geometries AliITSv11.
230     // In general it is automatically replaced by
231     // Inputs:
232     //   none.
233     // Outputs:
234     //   none.
235     // Return
236     //   none.
237
238     //TGeoMaterial *C  = new TGeoMaterial("ITSCarbon",12.0,6.0,2.265);
239     TGeoMaterial *matAl = new TGeoMaterial("ITSAluminum",26.981539,13.0,2.07);
240     TGeoMixture *matCfiber = new TGeoMixture("ITSCarbonFiber",6,1.930);
241     TGeoMixture *matRohacell = new TGeoMixture("ITSRohacell",6,1.930);
242     TGeoMixture *matStaselite = new TGeoMixture("ITSStaselite4411w",6,1.930);
243     TGeoMixture *matAir = new TGeoMixture("ITSAir",6,1.205*1.E-3);
244     TGeoMixture *matStainless = new TGeoMixture("ITSStainless",6,1.930);
245     //
246     Double_t medSPDcone[20];
247     medSPDcone[0] = 1.0; // imat
248     medSPDcone[1] = 0.0; // isvol
249     medSPDcone[2] = gAlice->Field()->Integ(); // ifield
250     medSPDcone[3] = gAlice->Field()->Max(); // fieldm
251     medSPDcone[4] = 1.0; // tmaxfd [degrees]
252     medSPDcone[5] = 1.0; // stemax [cm]
253     medSPDcone[6] = 0.5; // deemax [fraction]
254     medSPDcone[7] = 1.0E-3; // epsil [cm]
255     medSPDcone[8] = 0.0; // stmin [cm]
256     new TGeoMedium("ITSspdCarbonFiber",1,matCfiber,medSPDcone);
257     medSPDcone[0] += 1.0;
258     new TGeoMedium("ITSspdStaselite4411w",2,matStaselite,medSPDcone);
259     medSPDcone[0] += 1.0;
260     new TGeoMedium("ITSspdRohacell50A",3,matRohacell,medSPDcone);
261     medSPDcone[0] += 1.0;
262     new TGeoMedium("ITSspdStainlesSteal",4,matStainless,medSPDcone);
263     medSPDcone[0] += 1.0;
264     new TGeoMedium("ITSspdAir",5,matAir,medSPDcone);
265     medSPDcone[0] += 1.0;
266     new TGeoMedium("ITSspdAl",6,matAl,medSPDcone);
267     //
268     Double_t medSSDcone[20];
269     medSSDcone[0] = 1.0; // imat
270     medSSDcone[1] = 0.0; // isvol
271     medSSDcone[2] = gAlice->Field()->Integ(); // ifield
272     medSSDcone[3] = gAlice->Field()->Max(); // fieldm
273     medSSDcone[4] = 1.0; // tmaxfd [degrees]
274     medSSDcone[5] = 1.0; // stemax [cm]
275     medSSDcone[6] = 0.5; // deemax [fraction]
276     medSSDcone[7] = 1.0E-3; // epsil [cm]
277     medSSDcone[8] = 0.0; // stmin [cm]
278     new TGeoMedium("ITSssdCarbonFiber",1,matCfiber,medSSDcone);
279     medSSDcone[0] += 1.0;
280     new TGeoMedium("ITSssdStaselite4411w",2,matStaselite,medSSDcone);
281     medSSDcone[0] += 1.0;
282     new TGeoMedium("ITSssdRohacell50A",3,matRohacell,medSSDcone);
283     medSSDcone[0] += 1.0;
284     new TGeoMedium("ITSssdStainlesSteal",4,matStainless,medSSDcone);
285     medSSDcone[0] += 1.0;
286     new TGeoMedium("ITSssdAir",5,matAir,medSSDcone);
287     medSSDcone[0] += 1.0;
288     new TGeoMedium("ITSssdAl",6,matAl,medSSDcone);
289 }
290 //______________________________________________________________________
291 void AliITSv11::InitAliITSgeom(){
292     // Based on the geometry tree defined in Geant 3.21, this
293     // routine initilizes the Class AliITSgeom from the Geant 3.21 ITS 
294     // geometry sturture.
295     // Inputs:
296     //   none.
297     // Outputs:
298     //   none.
299     // Return
300     //   none.
301 }
302 //______________________________________________________________________
303 void AliITSv11::Init(){
304     // Initialise the ITS after it has been created.
305     // Inputs:
306     //   none.
307     // Outputs:
308     //   none.
309     // Return
310     //   none.
311 }
312 //______________________________________________________________________
313 void AliITSv11::SetDefaults(){
314     // Sets the default segmentation, response, digit and raw cluster 
315     // classes to be used. These defaults can be overwritten in the
316     // macros that do these later steps. Defaults are give hear for the
317     // general user.
318     // Inputs:
319     //   none.
320     // Outputs:
321     //   none.
322     // Return
323     //   none.
324 }
325 //______________________________________________________________________
326 void AliITSv11::DrawModule(){
327     // Draw a standard set of shaded view of the ITS version 11.
328     // Inputs:
329     //   none.
330     // Outputs:
331     //   none.
332     // Return
333     //   none.
334 }
335 //______________________________________________________________________
336 void AliITSv11::StepManager(){
337     // Called for every step in the ITS, then calles the AliITShit class
338     // creator with the information to be recoreded about that hit.
339     //  The value of the macro ALIITSPRINTGEOM if set to 1 will allow the
340     // printing of information to a file which can be used to create a .det
341     // file read in by the routine CreateGeometry(). If set to 0 or any other
342     // value except 1, the default behavior, then no such file is created nor
343     // is the extra variables and the like used in the printing allocated.
344 }
345