/[PAMELA software]/PamVMC/ac/src/PamVMCDetCat.cxx
ViewVC logotype

Contents of /PamVMC/ac/src/PamVMCDetCat.cxx

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.5 - (show annotations) (download)
Fri Jun 12 18:39:02 2009 UTC (15 years, 5 months ago) by pam-rm2
Branch: MAIN
CVS Tags: v1r0, HEAD
Changes since 1.1: +0 -0 lines 
- Introduced user-defined names of output files and random seeds number.
Users can do it use options of PamVMCApplication constructor:
PamVMCApplication(const char* name,  const char *title, const char*
filename="pamtest", Int_t seed=0).
The Random object that I use is TRandom3 object which has astronomical
large period (in case of default initialization 0). All random generators
in the code use this object by calling of gRandom singleton which keeps
it.

- Corrected TOF digitization routine. No problems with TDC hits due to
hadronic interactions anymore.

- Some small changes was done to compile code under Root 5.23. +
geant4_vmc v. 2.6 without any warnings

- Some classes of PamG4RunConfiguartion was changed for geant4_vmc v.
2.6.Some obsolete classes was deleted as soon as developers implemented
regions.

- Navigation was changed from "geomRootToGeant4" to "geomRoot", because on
VMC web page written that as soon as Geant4 has no option ONLY/MANY
translation of overlapped geometry to Geant4 through VGM could be wrong.
I'd like to stay with Root navigation:
http://root.cern.ch/root/vmc/Geant4VMC.html. This should be default
option.

- New Tracker digitization routine written by Sergio was implemented

- PamVMC again became compatible with geant4_vmc v.2.5 and ROOT 5.20.
 The problem was that ROOT developers introduced in TVirtualMC class a new
method SetMagField and new base class:TVirtualMagField from which
user-defined classes shoukd be derived
1 #include "PamVMCDetCat.h"
2
3 ClassImp(PamVMCDetCat)
4 void PamVMCDetCat::DefineGeometry(){
5
6 TGeoMaterial *m=0;
7 TGeoMixture *x=0;
8 TString s = "ALUMINIUM";
9 m = GetMat(s);
10 if(!m){
11 m = new TGeoMaterial(s,26.98,13,2.7);
12 m->SetUniqueID( 9);
13 }
14
15
16 s = "SCINT";
17 x = (TGeoMixture*)GetMat(s);
18 if(!x){
19 x = new TGeoMixture(s,2, 1.03200 );
20 x->SetUniqueID( 17);
21 x->DefineElement(0,12.01,6,0.922427);
22 x->DefineElement(1,1.01,1,0.7757296E-01);
23 }
24
25 s = "N2GAS";
26 m = GetMat(s);
27 if (!m){
28 m = new TGeoMaterial(s,14.01,7,0.1250000E-02);
29 m->SetUniqueID( 19);
30 }
31
32 s = "PLAST";
33 x = (TGeoMixture*)GetMat(s);
34 if(!x){
35 x = new TGeoMixture(s,2, 1.03200 );
36 x->SetUniqueID( 28);
37 x->DefineElement(0,12.01,6,0.9929774);
38 x->DefineElement(1,1.01,1,0.7022570E-02);
39 }
40
41 s = "TIALVA";
42 x = (TGeoMixture*)GetMat(s);
43 if(!x){
44 x = new TGeoMixture(s,3, 4.49234 );
45 x->SetUniqueID( 32);
46 x->DefineElement(0,47.867,22,0.9);
47 x->DefineElement(1,26.9815,13,0.6000000E-01);
48 x->DefineElement(2,50.9415,23,0.4000000E-01);
49 }
50
51 TGeoMedium *n =0;
52
53 s="ALUMINIUM";
54 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,1,9,-1,1,40,1,100,0.5000000E-01,0.136283,0.5000000E-03);
55
56 s = "N2 GAS";
57 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,3,19,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.1000000E-02,0.5000000E-03);
58
59 s = "SCINTILLATOR";
60 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,10,17,1,1,40,1,0.5000000E-01,0.5000000E-01,0.1000000E-01,0.5000000E-03);
61
62 s = "PLASTIC";
63 n = (GetMed(s))?GetMed(s): new TGeoMedium(s,18,28,-1,1,40,1,100,0.5000000E-01,0.1000000E-01,0.5000000E-03);
64
65 s = "TITANIUM Mixture";
66 n = (GetMed(s))?GetMed(s):new TGeoMedium(s,24,32,-1,1,40,1,0.5000000E-01,0.5000000E-01,0.3200000E-01,0.5000000E-03);
67
68
69 TGeoVolume *v = 0;
70
71 TPLAfirst = kFALSE;
72 s = "TPLA";
73 v = GetVol(s);
74 if(!v) {
75 v = gGeoManager->MakeBox("TPLA",GetMed("N2 GAS"),28.,27.165,1.6);
76 TPLAfirst = kTRUE;
77 }
78
79 s = "TPTU";
80 v = (GetVol(s))?GetVol(s):gGeoManager->MakeBox(s,GetMed("TITANIUM Mixture"),2.55,1.95,0.6);
81 s = "TPTM";
82 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("TITANIUM Mixture"),3.45,3.3,0.75);
83 s = "TPTL";
84 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("TITANIUM Mixture"),21.9,21.1,0.25);
85 s = "TPCV";
86 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),8.5,7.,0.25);
87 s = "CATA";
88 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),28.,0.6,27.165);
89 s = "ATBB";
90 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),28.,0.6,27.165);
91 s = "PLA2";
92 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("PLASTIC"),27.65,0.4,26.815);
93 s = "TOP1";
94 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("SCINTILLATOR"),26.85,0.4,26.015);
95 s = "PTBM";
96 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("PLASTIC"),10.11,0.4,8.44);
97 s = "ATBM";
98 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),9.41,0.4,7.74);
99 s = "LTBM";
100 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),9.11,0.6,7.44);
101 s = "PTBC";
102 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("PLASTIC"),4.975,0.4,5.5575);
103 s = "ATBC";
104 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("ALUMINIUM"),4.975,0.4,5.5575);
105 s = "LTBC";
106 v = (GetVol(s))?GetVol(s): gGeoManager->MakeBox(s,GetMed("N2 GAS"),4.975,0.6,5.5575);
107 s = "PTRT";
108 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("PLASTIC"),0.,16.9,0.4,2.405);
109 s = "ATRT";
110 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("ALUMINIUM"),0.,17.7,0.4,2.405);
111 s = "LTRT";
112 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("N2 GAS"),0.,18.05,0.6,2.405);
113 s = "PTRB";
114 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("PLASTIC"),16.9,0.,0.4,2.405);
115 s = "ATRB";
116 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("ALUMINIUM"),17.7,0.,0.4,2.405);
117 s = "LTRB";
118 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("N2 GAS"),18.05,0.,0.6,2.405);
119 s = "PTRL";
120 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("PLASTIC"),0.,14.9,0.4,2.375);
121 s = "ATRL";
122 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("ALUMINIUM"),0.,15.7,0.4,2.375);
123 s = "LTRL";
124 v = (GetVol(s))?GetVol(s): gGeoManager->MakeTrd1(s,GetMed("N2 GAS"),0.,16.05,0.6,2.375);
125
126 GetVol("TPTL")->AddNode(GetVol("TPCV"),1,gGeoIdentity);
127 GetVol("CATA")->AddNode(GetVol("ATBB"),1,gGeoIdentity); //only
128 GetVol("ATBB")->AddNodeOverlap(GetVol("PLA2"),1,gGeoIdentity); //only
129 GetVol("PLA2")->AddNodeOverlap(GetVol("TOP1"),1,gGeoIdentity); //only
130 GetVol("TOP1")->AddNode(GetVol("PTBM") ,1,gGeoIdentity);//only
131 GetVol("TOP1")->AddNode(GetVol("PTBC"),1,new TGeoTranslation(-21.875,0.,20.4575));
132 GetVol("TOP1")->AddNode(GetVol("PTBC"),2,new TGeoTranslation(21.875,0.,20.4575));
133 GetVol("TOP1")->AddNode(GetVol("PTBC"),3,new TGeoTranslation(-21.875,0.,-20.4575));
134 GetVol("TOP1")->AddNode(GetVol("PTBC"),4,new TGeoTranslation(21.875,0.,-20.4575));
135 GetVol("TOP1")->AddNode(GetVol("PTRT"),1,new TGeoTranslation(0.,0.,23.61));
136 GetVol("TOP1")->AddNode(GetVol("PTRB"),1,new TGeoTranslation(0.,0.,-23.61));
137 GetVol("TOP1")->AddNode(GetVol("PTRL"),1,new TGeoCombiTrans(-24.475,0.,0.,GetRot("rot7")));
138 GetVol("TOP1")->AddNode(GetVol("PTRL"),2,new TGeoCombiTrans(24.475,0.,0.,GetRot("rot8")));
139 GetVol("PLA2")->AddNode(GetVol("ATBM"),1,gGeoIdentity);
140 GetVol("PLA2")->AddNode(GetVol("ATBC"),1,new TGeoTranslation(-22.675,0.,21.2575));
141 GetVol("PLA2")->AddNode(GetVol("ATBC"),2,new TGeoTranslation(22.675,0.,21.2575));
142 GetVol("PLA2")->AddNode(GetVol("ATBC"),3,new TGeoTranslation(-22.675,0.,-21.2575));
143 GetVol("PLA2")->AddNode(GetVol("ATBC"),4,new TGeoTranslation(22.675,0.,-21.2575));
144 GetVol("PLA2")->AddNode(GetVol("ATRT"),1,new TGeoTranslation(0.,0.,24.41));
145 GetVol("PLA2")->AddNode(GetVol("ATRB"),1,new TGeoTranslation(0.,0.,-24.41));
146 GetVol("PLA2")->AddNode(GetVol("ATRL"),1,new TGeoCombiTrans(-25.275,0.,0.,GetRot("rot7")));
147 GetVol("PLA2")->AddNode(GetVol("ATRL"),2,new TGeoCombiTrans(25.275,0.,0.,GetRot("rot8")));
148 GetVol("ATBB")->AddNode(GetVol("LTBM"),1,gGeoIdentity);
149 GetVol("ATBB")->AddNode(GetVol("LTBC"),1,new TGeoTranslation(-23.025,0.,21.6075));
150 GetVol("ATBB")->AddNode(GetVol("LTBC"),2,new TGeoTranslation(23.025,0.,21.6075));
151 GetVol("ATBB")->AddNode(GetVol("LTBC"),3,new TGeoTranslation(-23.025,0.,-21.6075));
152 GetVol("ATBB")->AddNode(GetVol("LTBC"),4,new TGeoTranslation(23.025,0.,-21.6075));
153 GetVol("ATBB")->AddNode(GetVol("LTRT"),1,new TGeoTranslation(0.,0.,24.76));
154 GetVol("ATBB")->AddNode(GetVol("LTRB"),1,new TGeoTranslation(0.,0.,-24.76));
155 GetVol("ATBB")->AddNode(GetVol("LTRL"),1,new TGeoCombiTrans(-25.625,0.,0.,GetRot("rot7")));
156 GetVol("ATBB")->AddNode(GetVol("LTRL"),2,new TGeoCombiTrans(25.625,0.,0.,GetRot("rot8")));
157 GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),1,new TGeoTranslation(-20.65,-20.15,1.));
158 GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),2,new TGeoTranslation(-20.65,20.15,1.));
159 GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),3,new TGeoTranslation(20.65,-20.15,1.));
160 GetVol("TPLA")->AddNodeOverlap(GetVol("TPTU"),4,new TGeoTranslation(20.65,20.15,1.));
161
162
163 if(TPLAfirst) SetMotherProp(GetVol("TPLA"),1,new TGeoTranslation(0.,0.,73.489));
164
165 GetVol("TPLA")->AddNode(GetVol("TPTL"),1,new TGeoTranslation(0.,0.,-1.35));
166 GetVol("TPLA")->AddNode(GetVol("TPTM"),1,new TGeoTranslation(-19.75,-18.8,-0.35));
167 GetVol("TPLA")->AddNode(GetVol("TPTM"),2,new TGeoTranslation(-19.75,18.8,-0.35));
168 GetVol("TPLA")->AddNode(GetVol("TPTM"),3,new TGeoTranslation(19.75,-18.8,-0.35));
169 GetVol("TPLA")->AddNode(GetVol("TPTM"),4,new TGeoTranslation(19.75,18.8,-0.35));
170 GetVol("TPLA")->AddNode(GetVol("CATA"),1,new TGeoCombiTrans(0.,0.,1.,GetRot("rot2")));
171 };
172
173 void PamVMCDetCat::DefineCuts(){
174
175 TString s ="ALUMINIUM"; // default GPAMELA parameters
176 if (GetMed(s) && !GetCC(s))
177 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
178 0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
179 s ="N2 GAS";
180 if (GetMed(s) && !GetCC(s))
181 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
182 0.001, 0.01, 1., 1., 1., 0., 0., 1., 1., 4., 1., 1., 2., 1.));
183 s ="SCINTILLATOR";// default GPAMELA parameters
184 if (GetMed(s) && !GetCC(s))
185 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
186 0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
187 s = "PLASTIC";// default GPAMELA parameters
188 if (GetMed(s) && !GetCC(s))
189 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
190 0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
191 s = "TITANIUM Mixture"; // default GPAMELA parameters
192 if (GetMed(s) && !GetCC(s))
193 SetCC(s, new pCutControl(GetMedID(s), 0.0001, 0.001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001,
194 0.001, 0.01, 1., 1., 1., 0., 1., 1., 1., 4., 1., 1., 1., 1.));
195
196 }
  ViewVC Help
Powered by ViewVC 1.1.23