Log replaced by Id
[u/mrichter/AliRoot.git] / ITS / AliITSGeometryITSV.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 #include <Riostream.h>
19 #include <stdio.h>
20 #include <stdlib.h>
21 #include <TMath.h>
22 #include <TGeometry.h>
23 #include <TNode.h>
24 #include <TTUBE.h>
25 #include <TTUBS.h>
26 #include <TPCON.h>
27 #include <TFile.h>    // only required for Tracking function?
28 #include <TCanvas.h>
29 #include <TObjArray.h>
30 #include <TLorentzVector.h>
31 #include <TObjString.h>
32 #include <TClonesArray.h>
33 #include <TBRIK.h>
34 #include <TSystem.h>
35 #include <TVector3.h>
36 #include <AliRun.h>
37 #include <AliITS.h>
38
39 #include "AliITSGeometryITSV.h"
40
41 ClassImp(AliITSGeometryITSV)
42
43 //______________________________________________________________________
44 AliITSGeometryITSV::AliITSGeometryITSV() : AliITSBaseGeometry(){
45     //Default Constructor for SSD Cone geometry
46     // Inputs:
47     //   none.
48     // Outputs:
49     //   none.
50     // Return:
51     //   none.
52
53     SetScalemm();
54 }
55 //______________________________________________________________________
56 AliITSGeometryITSV::AliITSGeometryITSV(AliITS *its,const char *moth):
57     AliITSBaseGeometry(its,0){
58     //Standard Constructor for SSD Cone geometry
59     // Inputs:
60     //   const char *moth, The volume name into which the ITS mother volume
61     //                     will reside.
62     // Outputs:
63     //   none.
64     // Return:
65     //   none.
66
67     fAir = 5; // Air material number.
68
69     fA.Size(16,"ITS Mother Volume");
70     fA.SetVid(ITSG3VnameToIndex("ITSV")); // ITSV is a special name.
71     Double_t rlim  = 50.;
72     Double_t zmax  = 74.;
73     Double_t ztpc = 284.;
74     //
75     fA.P0() = 0.;
76     fA.dP() = 360.;
77     //
78     fA.Z(0)  = -ztpc-5.0-0.1;
79     fA.Rn(0) = 46.0;   
80     fA.Rx(0) = 85.0;
81     //
82     fA.Z(1) = -ztpc;
83     fA.Rn(1) = 46;   
84     fA.Rx(1) = 85.;
85     //
86     fA.Z(2) = -ztpc;
87     fA.Rn(2) = 46;  
88     fA.Rx(2) = rlim+6;
89     //
90     fA.Z(3) = -97.5;
91     fA.Rn(3) = 46;  
92     fA.Rx(3) = rlim+6;
93     //
94     fA.Z(4) = -zmax;
95     fA.Rn(4) = 46;  
96     fA.Rx(4) = rlim+6;
97     //
98     fA.Z(5) = -48;   
99     fA.Rn(5) = 6;
100     fA.Rx(5) = rlim+6;
101     //
102     fA.Z(6) = -28.6;   
103     fA.Rn(6) = 6;
104     fA.Rx(6) = rlim+6;
105     //
106     fA.Z(7) = -27.6;  
107     fA.Rn(7) = 3.295;
108     fA.Rx(7) = rlim+6; 
109     //
110     fA.Z(8) = 27.6;
111     fA.Rn(8) = 3.295;
112     fA.Rx(8) = rlim+6;
113     //
114     fA.Z(9) = 28.6;
115     fA.Rn(9) = 6;
116     fA.Rx(9) = rlim+6;
117     //
118     fA.Z(10) = 48;   
119     fA.Rn(10) = 6;
120     fA.Rx(10) = rlim+6;
121     //
122     fA.Z(11) = zmax;
123     fA.Rn(11) = 46;
124     fA.Rx(11) = rlim+6;
125     //
126     fA.Z(12) = 97.5;
127     fA.Rn(12) = 46;  
128     fA.Rx(12) = rlim+6;
129     //
130     fA.Z(13) = ztpc;
131     fA.Rn(13) = 62;     
132     fA.Rx(13) = 62+4.;
133     //
134     fA.Z(14) = ztpc;
135     fA.Rn(14) = 62;     
136     fA.Rx(14) = 85.;
137     //
138     fA.Z(15) = ztpc+4.+0.1;
139     fA.Rn(15) = 62.4;
140     fA.Rx(15) = 85.;
141 }
142 //______________________________________________________________________
143 void AliITSGeometryITSV::CreateG3Geometry(){
144     // Calls Geant 3 geometry inilization routines with the information
145     // stored in this class.
146     // Inputs:
147     //    none.
148     // Outputs:
149     //    none.
150     // Return:
151     //    none.
152
153     PolyCone(fA,fAir);
154     return;
155 }
156 //______________________________________________________________________
157 void AliITSGeometryITSV::PositionGeometry(const char *moth,Int_t cn,
158                                           TVector3 &t,Int_t irot){
159     // Positions the ITSV geometry in the way needed by Geant3.
160     // Inputs:
161     //   const char *moth, The volume name into which the ITS mother volume
162     //                     will reside.
163     //   Int_t copy        copy number.
164     //   TVector &t        Translection vector
165     //   Int_t   irot      Rotation matrix index number
166     // Outputs:
167     //    none.
168     // Return:
169     //    none.
170
171     gMC->Gspos("ITSV",1,moth,t.X(),t.Y(),t.Z(),irot,"ONLY");
172     return;
173 }
174 //______________________________________________________________________
175 void AliITSGeometryITSV::CreateG3Materials(){
176     // Fills the Geant 3 banks with Material and Medium definisions.
177     // Inputs:
178     //   none.
179     // Outputs:
180     //   none.
181     // Returns:
182     //   none.
183
184     // Dry Air = 78.084% N2, 20.946% O2, 0.934% Ar, 0.033% CO2 by Volume
185     // "Handbook of Chemistry and Physics" 60th eddition, CRC-Press (1980)
186     // page F-211. Others components <0.003% and measured in parts per 
187     // million.
188     Int_t i;
189     Int_t    Z[4]={6,7,8,18}; //C, N, O, Ar
190     Double_t W[4],f[4]={0.00033,0.78084,0.20946,0.00934};//C02,N2,O2,Ar by vol
191     Double_t T = 293.15; // 20 degrees C.
192     Double_t P = 760.0; // mm of Hg.
193     // op. cite, page  F-9.
194     Double_t dens = 1.2929*(273.13/T)*(P/760.0); // Density in grams/leter
195
196     dens *= 1.0E-3; // Convert 1 leter = 1E3 cm^3.
197     // Covert fraction of compounds by volume to graction of atoms by weight.
198     W[0] = GetA(Z[0])*f[0]; // C
199     W[1] = GetA(Z[1])*2.0*f[1]; // N
200     W[2] = GetA(Z[2])*2.0*(f[0]+f[2]); // O
201     W[3] = GetA(Z[3])*f[3];
202     // Renormilize the weights.
203     f[0] = 0.0;
204     for(i=0;i<4;i++) f[0] += W[i];
205     for(i=0;i<4;i++) W[i] /= f[0];
206     MixtureByWeight(fAir,"Dry Standard ITS Air",Z,W,dens,4,0);
207 }
208 //______________________________________________________________________
209 void AliITSGeometryITSV::BuildDisplayGeometry(){
210     // Fill Root geometry banks for fast simple ITS simulation event
211     // display. See Display.C, and related code, for more details.
212     // Inputs:
213     //    none.
214     // Outputs:
215     //   none.
216     // Return:
217     //  none.
218
219     // No need to display ITS cones.
220 }
221
222 //______________________________________________________________________
223 void AliITSGeometryITSV::PolyCone(AliITSPConeData &d,Int_t med){
224     // Interface to TMC->Gsvolu() for ITS PCON geometry. Poly-cone It has 9 
225     // parameters or more. See SetScale() for units. Default units are geant
226     // 3 [cm].
227     // Inputs:
228     //    AliITSPConeData &d  Object with poly cone data stored in it.
229     //    Int_t    med        media index number.
230     // Output:
231     //    none.
232     // Return.
233     //    none.
234     Float_t *param;
235     Int_t n,i;
236
237     n = 3+3*d.Nz();
238     param = new Float_t[n];
239     param[0] = d.Phi0();
240     param[1] = d.DPhi();
241     param[2] = (Float_t) d.Nz();
242     for(i=0;i<d.Nz();i++){
243         param[3+3*i] = fScale*d.ZAt(i);
244         param[4+3*i] = fScale*d.Rmin(i);
245         param[5+3*i] = fScale*d.Rmax(i);
246     } // end for if
247     if(fVolName==0){ // must create array.
248         fVolNameSize = 38624;
249         fVolName = new TString[fVolNameSize];
250         fVolNameLast = 0;
251     } // end if
252     d.SetVid(ITSG3VnameToIndex("ITSV"));
253     fVolName[d.GetVid()] = (d.GetName())->Data();
254     gMC->Gsvolu("ITSV","PCON",GetMed(med),param,n);
255
256     delete[] param;
257 }
258