189{
191 TGeoMedium *HPTmat =gGeoManager->GetMedium("HPTgas");
192
194 TGeoMedium *Conc =gGeoManager->GetMedium("Concrete");
195
196
198 TGeoVolume *volMudetBox = gGeoManager->GetVolume("volNuTauMudet");
199 TGeoBBox *HPT =
new TGeoBBox(
"HPT",
DimX/2,
DimY/2,
DimZ/2);
200 TGeoVolume *volHPT = new TGeoVolume("volHPT",HPT,HPTmat);
201 volHPT->SetLineColor(kBlue-5);
202 AddSensitiveVolume(volHPT);
203
205 TGeoVolume *volCbase = new TGeoVolume("volCbase",Cbase,Conc);
206 volCbase->SetLineColor(kOrange-7);
207
208
209 volMudetBox->AddNode(volHPT,1,
new TGeoTranslation(0,0,-
zSizeMudet/2 +
DimZ/2));
211
212
213
214 volMudetBox->AddNode(volHPT,2,
new TGeoTranslation(0,0,-
zSizeMudet/2 + 3*
DimZ/2 +55*
cm));
216
217
218 volMudetBox->AddNode(volHPT,3,
new TGeoTranslation(0,0,-72*
cm/2 -
DimZ/2));
219
220
221 volMudetBox->AddNode(volHPT,4,
new TGeoTranslation(0,0,72*
cm/2 +
DimZ/2));
222
223
224
225 volMudetBox->AddNode(volHPT,5,
new TGeoTranslation(0,0,
zSizeMudet/2 - 3*
DimZ/2 - 55*
cm));
227
228 volMudetBox->AddNode(volHPT,6,
new TGeoTranslation(0,0,
zSizeMudet/2 -
DimZ/2));
230
231 }
233
234 TGeoVolume *volMagRegion=gGeoManager->GetVolume("volMagRegion");
235 TGeoVolume *volTarget =gGeoManager->GetVolume("volTarget");
236 TGeoVolume *tTauNuDet = gGeoManager->GetVolume("tTauNuDet");
237
238 Double_t DZMagnetizedRegion = ((TGeoBBox*) volMagRegion->GetShape())->GetDZ() *2;
239 Double_t DYMagnetizedRegion = ((TGeoBBox*) volMagRegion->GetShape())->GetDY() *2;
240 Double_t DXMagnetizedRegion = ((TGeoBBox*) volMagRegion->GetShape())->GetDX() *2;
241
242 Double_t DZTarget = ((TGeoBBox*) volTarget->GetShape())->GetDZ() *2;
243
244
246 TGeoVolume *volDT = new TGeoVolume("volDT",DT,HPTmat);
247 volDT->SetLineColor(kBlue-5);
248
249
251 TGeoMedium *CarbonComposite = gGeoManager->GetMedium("CarbonComposite");
252
254 TGeoMedium *SciFiMat = gGeoManager->GetMedium("SciFiMat");
255
257 TGeoMedium *Airex = gGeoManager->GetMedium("Airex");
258
259
261 TGeoVolume* HPT_support_volume = new TGeoVolume("HPT_support", HPT_support_box, CarbonComposite);
262 HPT_support_volume->SetLineColor(kGray - 2);
263 HPT_support_volume->SetVisibility(1);
264
265
267 TGeoVolume* HPT_honeycomb_volume = new TGeoVolume("HPT_honeycomb", HPT_honeycomb_box, Airex);
268 HPT_honeycomb_volume->SetLineColor(kYellow);
269 HPT_honeycomb_volume->SetVisibility(1);
270
271
273 TGeoVolume* HPT_scifi_plane_hor_volume = new TGeoVolume("HPT_scifi_plane_hor", HPT_scifi_plane_hor_box, SciFiMat);
274 HPT_scifi_plane_hor_volume->SetVisibility(1);
275
277 TGeoVolume* HPT_scifi_plane_vert_volume = new TGeoVolume("HPT_scifi_plane_vert", HPT_scifi_plane_vert_box, SciFiMat);
278 HPT_scifi_plane_vert_volume->SetVisibility(1);
279
280
282 TGeoVolume* HPT_scifimat_hor_volume = new TGeoVolume("HPT_scifimat_hor", HPT_scifimat_hor_box, SciFiMat);
283 HPT_scifimat_hor_volume->SetLineColor(kCyan);
284
286 TGeoVolume* HPT_scifimat_vert_volume = new TGeoVolume("HPT_scifimat_vert", HPT_scifimat_vert_box, SciFiMat);
287 HPT_scifimat_vert_volume->SetLineColor(kGreen);
288
289 AddSensitiveVolume(HPT_scifimat_hor_volume);
290 AddSensitiveVolume(HPT_scifimat_vert_volume);
291
292
294 HPT_scifi_plane_hor_volume->AddNode(HPT_scifimat_hor_volume, i+1,
new TGeoTranslation(0, (-(
n_hor_planes-1)/2.0 + i)*
scifimat_width, 0));
295 }
297 HPT_scifi_plane_vert_volume->AddNode(HPT_scifimat_vert_volume, 100+i+1 ,
new TGeoTranslation((-(
n_vert_planes-1)/2.0 + i)*
scifimat_width, 0, 0));
298 }
299
300 volDT->AddNode(HPT_support_volume, 0,
new TGeoTranslation(0, 0, -
DimZ / 2 +
support_z / 2));
301 volDT->AddNode(HPT_scifi_plane_hor_volume, 0,
new TGeoTranslation(0, 0, -
DimZ / 2 +
support_z +
scifimat_z / 2));
306
307 Double_t first_DT_position = -DZMagnetizedRegion/2 + DZTarget +
DimZ/2;
311 }
312 }
313 }
314}
Int_t InitMedium(const char *name)