BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
BesDetectorConstruction.cc
Go to the documentation of this file.
1//---------------------------------------------------------------------------//
2// BOOST --- BESIII Object_Oriented Simulation Tool //
3//---------------------------------------------------------------------------//
4// Description: Using gdml in G4
5// Author: Liuhm
6// Created: May 25, 2003
7// Modified:
8// Comment:
9//---------------------------------------------------------------------------//
10//
11#include "BesDetectorConstruction.hh"
12#include "BesMagneticField.hh"
13
14#include "EmcSim/BesEmcConstruction.hh"
15#include "MdcSim/BesMdcConstruction.hh"
16#include "MucSim/BesMucConstruction.hh"
17#include "TofSim/BesTofConstruction.hh"
18
19#include "G4Box.hh"
20#include "G4Colour.hh"
21#include "G4Evaluator/GDMLExpressionEvaluator.hh"
22#include "G4Processor/GDMLProcessor.hh"
23#include "G4Tubs.hh"
24#include "G4VisAttributes.hh"
25#include "SimUtil/ReadBoostRoot.hh"
26#include "globals.hh"
27#include <cstdlib>
28#include <fstream>
29// Added here just to help resolve properly dependencies
30#include "G4BooleanSolid.hh"
31#include "G4Geo/BesG4Geo.h"
32
33BesDetectorConstruction::BesDetectorConstruction() {
34 // create commands for interactive definition of the calorimeter
35 // detectorMessenger = new BesDetectorMessenger(this);
36
37 if ( ReadBoostRoot::GetField() ) besField = new BesMagneticField();
38
39 mdc = new BesMdcConstruction();
40 tof = new BesTofConstruction();
41 emc = new BesEmcConstruction();
42 muc = new BesMucConstruction();
43 m_SCM = new BesSCM();
44 m_Pip = new BesPip();
45 m_pipeSCM = 2;
46 lWorld = 0;
47 fWorld = 0;
48}
49
50BesDetectorConstruction::~BesDetectorConstruction() {
51 G4cout << "BesDetectorConstruction::~BesDetectorConstruction()" << G4endl;
52 if ( ReadBoostRoot::GetField() && besField ) delete besField;
53 if ( mdc ) delete mdc;
54 if ( tof ) delete tof;
55 if ( emc ) delete emc;
56 if ( muc ) delete muc;
57 if ( m_SCM ) delete m_SCM;
58 if ( m_Pip ) delete m_Pip;
59}
60
61G4VPhysicalVolume* BesDetectorConstruction::Construct() {
62 if ( m_pipeSCM != 2 )
63 {
64 // define materials
65 G4double a, z, density, abundance, fractionmass;
66 G4double temperature, pressure;
67 G4int iz, n, nel, ncomponents, natoms;
68 G4String name, symbol;
69 G4Isotope* U5 = new G4Isotope( name = "U235", iz = 92, n = 235, a = 235.01 * g / mole );
70 G4Isotope* U8 = new G4Isotope( name = "U238", iz = 92, n = 238, a = 238.03 * g / mole );
71
72 G4Element* U = new G4Element( name = "enriched Uranium", symbol = "U", ncomponents = 2 );
73 U->AddIsotope( U5, abundance = 90. * perCent );
74 U->AddIsotope( U8, abundance = 10. * perCent );
75
76 a = 1.01 * g / mole;
77 G4Element* H = new G4Element( name = "Hydrogen", symbol = "H", z = 1., a );
78
79 a = 12.01 * g / mole;
80 G4Element* C = new G4Element( name = "Carbon", symbol = "C", z = 6., a );
81
82 a = 14.01 * g / mole;
83 G4Element* N = new G4Element( name = "Nitrogen", symbol = "N", z = 7., a );
84
85 a = 16.00 * g / mole;
86 G4Element* O = new G4Element( name = "Oxygen", symbol = "O", z = 8., a );
87
88 a = 28.09 * g / mole;
89 G4Element* Si = new G4Element( name = "Silicon", symbol = "Si", z = 14., a );
90
91 a = 126.90447 * g / mole;
92 G4Element* I = new G4Element( name = "Iodine", symbol = "I", z = 53., a );
93
94 a = 132.90545 * g / mole;
95 G4Element* Cs = new G4Element( name = "Cesium", symbol = "Cs", z = 55., a );
96
97 density = 1.848 * g / cm3;
98 a = 9.012182 * g / mole;
99 G4Material* Beryllium = new G4Material( name = "Beryllium", z = 4.0, a, density );
100
101 density = 2.70 * g / cm3;
102 a = 26.98 * g / mole;
103 G4Material* Aluminium = new G4Material( name = "Aluminium", z = 13., a, density );
104
105 density = 7.84 * g / cm3;
106 a = 55.845 * g / mole;
107 G4Material* Iron = new G4Material( name = "Iron", z = 26.0, a, density );
108
109 density = 8.96 * g / cm3;
110 a = 63.546 * g / mole;
111 G4Material* Copper = new G4Material( name = "Copper", z = 29.0, a, density );
112
113 density = 0.0001664 * g / cm3;
114 a = 4.0026 * g / mole;
115 G4Material* Hegas = new G4Material( name = "Hegas", z = 2.0, a, density );
116
117 density = 0.001879 * g / cm3;
118 G4Material* Propane = new G4Material( name = "Propane", density, nel = 2 );
119 Propane->AddElement( C, natoms = 3 );
120 Propane->AddElement( H, natoms = 8 );
121
122 density = 4.53 * g / cm3;
123 G4Material* Cesiumiodide = new G4Material( name = "Cesiumiodide", density, nel = 2 );
124 Cesiumiodide->AddElement( I, natoms = 1 );
125 Cesiumiodide->AddElement( Cs, natoms = 1 );
126
127 density = 0.00085144 * g / cm3;
128 G4Material* Mdcgas = new G4Material( name = "Mdcgas", density, nel = 2 );
129 Mdcgas->AddMaterial( Hegas, fractionmass = 0.1173 );
130 Mdcgas->AddMaterial( Propane, fractionmass = 0.8827 );
131
132 density = 1.390 * g / cm3;
133 a = 39.95 * g / mole;
134 G4Material* liquidArgon = new G4Material( name = "liquidArgon", z = 18.0, a, density );
135
136 density = 11.35 * g / cm3;
137 a = 207.19 * g / mole;
138 G4Material* Lead = new G4Material( name = "Lead", z = 82., a, density );
139
140 density = 1.0 * g / cm3;
141 G4Material* Water = new G4Material( name = "Water", density, nel = 2 );
142 Water->AddElement( H, natoms = 2 );
143 Water->AddElement( O, natoms = 1 );
144
145 density = 1.032 * g / cm3;
146 G4Material* Scintillator = new G4Material( name = "Scintillator", density, nel = 2 );
147 Scintillator->AddElement( C, natoms = 9 );
148 Scintillator->AddElement( H, natoms = 10 );
149
150 density = 2.200 * g / cm3;
151 G4Material* SiO2 = new G4Material( name = "SiO2", density, nel = 2 );
152 SiO2->AddElement( Si, natoms = 1 );
153 SiO2->AddElement( O, natoms = 2 );
154
155 density = 1.290 * mg / cm3;
156 G4Material* Air = new G4Material( name = "Air", density, nel = 2 );
157 Air->AddElement( N, fractionmass = 0.7 );
158 Air->AddElement( O, fractionmass = 0.3 );
159
160 density = 0.200 * g / cm3;
161 G4Material* Aerogel = new G4Material( name = "Aerogel", density, nel = 3 );
162 Aerogel->AddMaterial( SiO2, fractionmass = 0.625 );
163 Aerogel->AddMaterial( Water, fractionmass = 0.374 );
164 Aerogel->AddElement( C, fractionmass = 0.001 );
165
166 density = 27.0 * mg / cm3;
167 pressure = 50.0 * atmosphere;
168 temperature = 325.0 * kelvin;
169 G4Material* CarbonicGas = new G4Material( name = "CarbonicGas", density, nel = 2,
170 kStateGas, temperature, pressure );
171 CarbonicGas->AddElement( C, natoms = 1 );
172 CarbonicGas->AddElement( O, natoms = 2 );
173
174 density = 0.3 * mg / cm3;
175 pressure = 2.0 * atmosphere;
176 temperature = 500.0 * kelvin;
177 G4Material* WaterSteam = new G4Material( name = "WaterSteam", density, nel = 1, kStateGas,
178 temperature, pressure );
179 WaterSteam->AddMaterial( Water, fractionmass = 1 );
180
181 density = universe_mean_density;
182 pressure = 3.e-18 * pascal;
183 temperature = 2.73 * kelvin;
184 G4Material* Galactic = new G4Material( name = "Galactic", z = 1., a = 1.01 * g / mole,
185 density, kStateGas, temperature, pressure );
186
187 density = 1.0e-5 * g / cm3;
188 pressure = 2.e-2 * bar;
189 G4Material* Beam = new G4Material( name = "Beam", density, nel = 1, kStateGas,
190 STP_Temperature, pressure );
191 Beam->AddMaterial( Air, fractionmass = 1.0 );
192
193 // construct world volume
194 G4Box* sWorld = new G4Box( "solidWorld", 263.5 * cm, 263.5 * cm, 287.5 * cm );
195 lWorld = new G4LogicalVolume( sWorld, Air, "logicalWorld", 0, 0, 0 );
196 fWorld =
197 new G4PVPlacement( 0, G4ThreeVector( 0, 0, 0 ), lWorld, "physicWorld", 0, false, 0 );
198 // set visual attributes
199 lWorld->SetVisAttributes( G4VisAttributes::Invisible );
200
201 if ( m_pipeSCM == 1 )
202 {
203 G4cout << "construct beampipe and SCM with Geant4 classes" << G4endl;
204 // construct beam pipe
205 /*
206 G4Tubs* sPipe = new G4Tubs("solidPipe",0.00*cm,3.58*cm,15.0*cm,0.*deg,360*deg);
207 G4LogicalVolume* lPipe = new G4LogicalVolume(sPipe,Air,"logicalPipe",0,0,0);
208 G4VPhysicalVolume* fPipe = new
209 G4PVPlacement(0,G4ThreeVector(0,0,0),lPipe,"physicPipe",lWorld,false,0);
210
211 G4Tubs* sPip1 = new G4Tubs("solidPip1",3.15*cm,3.23*cm,15.0*cm,0.*deg,360*deg);
212 G4LogicalVolume* lPip1 = new G4LogicalVolume(sPip1,Beryllium,"logicalPip1",0,0,0);
213 G4VPhysicalVolume* fPip1 = new
214 G4PVPlacement(0,G4ThreeVector(0,0,0),lPip1,"physicPip1",lPipe,false,0);
215
216 G4Tubs* sPip2 = new G4Tubs("solidPip2",3.23*cm,3.53*cm,15.0*cm,0.*deg,360*deg);
217 G4LogicalVolume* lPip2 = new G4LogicalVolume(sPip2,Hegas,"logicalPip2",0,0,0);
218 G4VPhysicalVolume* fPip2 = new
219 G4PVPlacement(0,G4ThreeVector(0,0,0),lPip2,"physicPip2",lPipe,false,0);
220
221 G4Tubs* sPip3 = new G4Tubs("solidPip3",3.53*cm,3.58*cm,15.0*cm,0.*deg,360*deg);
222 G4LogicalVolume* lPip3 = new G4LogicalVolume(sPip3,Beryllium,"logicalPip3",0,0,0);
223 G4VPhysicalVolume* fPip3 = new
224 G4PVPlacement(0,G4ThreeVector(0,0,0),lPip3,"physicPip3",lPipe,false,0);
225 //set visual attributes
226 lPipe -> SetVisAttributes(G4VisAttributes::Invisible);
227 lPip1 -> SetVisAttributes(G4VisAttributes::Invisible);
228 lPip2 -> SetVisAttributes(G4VisAttributes::Invisible);
229 lPip3 -> SetVisAttributes(G4VisAttributes::Invisible);
230 */
231 m_Pip->Construct( lWorld );
232 m_SCM->Construct( lWorld );
233 /*
234 //construct SCM
235 G4Tubs* sCoil = new G4Tubs("solidCoil",137.5*cm,170.0*cm,195.5*cm,0.*deg,360*deg);
236 G4LogicalVolume* lCoil = new G4LogicalVolume(sCoil,Air,"logicalCoil",0,0,0);
237 G4VPhysicalVolume* fCoil = new
238 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoil,"physicCoil",lWorld,false,0);
239
240 G4Tubs* sCoilin = new
241 G4Tubs("solidCoilin",137.5*cm,138.1*cm,195.5*cm,0.*deg,360*deg); G4LogicalVolume*
242 lCoilin = new G4LogicalVolume(sCoilin,Iron,"logicalCoilin",0,0,0); G4VPhysicalVolume*
243 fCoilin = new
244 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoilin,"physicCoilin",lCoil,false,0);
245
246 G4Tubs* sCoilout = new
247 G4Tubs("solidCoilout",168.4*cm,170.0*cm,195.5*cm,0.*deg,360*deg); G4LogicalVolume*
248 lCoilout = new G4LogicalVolume(sCoilout,Iron,"logicalCoilout",0,0,0);
249 G4VPhysicalVolume* fCoilout = new
250 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoilout,"physicCoilout",lCoil,false,0);
251
252 G4Tubs* sCoil1 = new
253 G4Tubs("solidCoil1",144.5*cm,145.0*cm,179.5*cm,0.*deg,360*deg); G4LogicalVolume*
254 lCoil1 = new G4LogicalVolume(sCoil1,Aluminium,"logicalCoil1",0,0,0);
255 G4VPhysicalVolume* fCoil1 = new
256 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoil1,"physicCoil1",lCoil,false,0);
257
258 G4Tubs* sCoil2 = new
259 G4Tubs("solidCoil2",157.5*cm,158.0*cm,179.5*cm,0.*deg,360*deg); G4LogicalVolume*
260 lCoil2 = new G4LogicalVolume(sCoil2,Aluminium,"logicalCoil2",0,0,0);
261 G4VPhysicalVolume* fCoil2 = new
262 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoil2,"physicCoil2",lCoil,false,0);
263
264 G4Tubs* sCoil3 = new
265 G4Tubs("solidCoil3",148.0*cm,150.0*cm,179.5*cm,0.*deg,360*deg); G4LogicalVolume*
266 lCoil3 = new G4LogicalVolume(sCoil3,Aluminium,"logicalCoil3",0,0,0);
267 G4VPhysicalVolume* fCoil3 = new
268 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoil3,"physicCoil3",lCoil,false,0);
269
270 G4Tubs* sCoil4 = new
271 G4Tubs("solidCoil4",148.8*cm,149.2*cm,179.5*cm,0.*deg,360*deg); G4LogicalVolume*
272 lCoil4 = new G4LogicalVolume(sCoil4,Copper,"logicalCoil4",0,0,0); G4VPhysicalVolume*
273 fCoil4 = new
274 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoil4,"physicCoil4",lCoil3,false,0);
275
276 G4Tubs* sCoil5 = new
277 G4Tubs("solidCoil5",150.5*cm,152.7*cm,179.5*cm,0.*deg,360*deg); G4LogicalVolume*
278 lCoil5 = new G4LogicalVolume(sCoil5,Aluminium,"logicalCoil5",0,0,0);
279 G4VPhysicalVolume* fCoil5 = new
280 G4PVPlacement(0,G4ThreeVector(0,0,0),lCoil5,"physicCoil5",lCoil,false,0);
281 //set visual attributes
282 lCoil -> SetVisAttributes(G4VisAttributes::Invisible);
283 lCoilin -> SetVisAttributes(G4VisAttributes::Invisible);
284 lCoilout -> SetVisAttributes(G4VisAttributes::Invisible);
285 lCoil1 -> SetVisAttributes(G4VisAttributes::Invisible);
286 lCoil2 -> SetVisAttributes(G4VisAttributes::Invisible);
287 lCoil3 -> SetVisAttributes(G4VisAttributes::Invisible);
288 lCoil4 -> SetVisAttributes(G4VisAttributes::Invisible);
289 lCoil5 -> SetVisAttributes(G4VisAttributes::Invisible);
290 */
291 }
292 }
293
294 else if ( m_pipeSCM == 2 )
295 {
296 /*
297 BesG4Geo* aBesG4Geo = new BesG4Geo();
298
299 //construct with gdml file
300 //G4String GeometryPath = ReadBoostRoot::GetBoostRoot();
301 //G4String GeometryPath = getenv("GDMLMANAGEMENTDATAROOT");
302 if(!GeometryPath){
303 G4Exception("BOOST environment not set!");
304 }
305 GeometryPath += "/dat/Bes.gdml";
306
307 G4cout<<"construc beampipe and SCM from: "<<GeometryPath<<G4endl;
308
309 config.SetURI(GeometryPath);
310 config.SetSetupName( "Default" );
311 config.SetSetupVersion( "1.0" );
312 config.SetType( "G4" );
313 sxp.Configure( &config );
314 sxp.Initialize();
315 sxp.Run();
316
317
318 fWorld=(G4VPhysicalVolume *)GDMLProcessor::GetInstance()->GetWorldVolume();
319 if( fWorld == 0 ) {
320 G4Exception(
321 "World volume not set properly check your setup selection criteria or GDML input!");
322 }
323
324 //make the world volume invisible
325 lWorld = FindLogicalVolume("logicalWorld");
326 lWorld -> SetVisAttributes(G4VisAttributes::Invisible);
327
328 //make the material inside magnet coil invisible
329
330 G4LogicalVolume* lCoil = FindLogicalVolume("logicalCoil" );
331 lCoil -> SetVisAttributes(G4VisAttributes::Invisible);
332
333 G4LogicalVolume* lCoilin = FindLogicalVolume("logicalCoilin" );
334 lCoilin -> SetVisAttributes(G4VisAttributes::Invisible);
335
336 G4LogicalVolume* lCoilout = FindLogicalVolume("logicalCoilout" );
337 lCoilout -> SetVisAttributes(G4VisAttributes::Invisible);
338
339 G4LogicalVolume* lCoil1 = FindLogicalVolume("logicalCoil1");
340 lCoil1 -> SetVisAttributes(G4VisAttributes::Invisible);
341
342 G4LogicalVolume* lCoil2 = FindLogicalVolume("logicalCoil2");
343 lCoil2 -> SetVisAttributes(G4VisAttributes::Invisible);
344
345 G4LogicalVolume* lCoil3 = FindLogicalVolume("logicalCoil3");
346 lCoil3 -> SetVisAttributes(G4VisAttributes::Invisible);
347
348 G4LogicalVolume* lCoil4 = FindLogicalVolume("logicalCoil4");
349 lCoil4 -> SetVisAttributes(G4VisAttributes::Invisible);
350
351 G4LogicalVolume* lCoil5 = FindLogicalVolume("logicalCoil5");
352 lCoil5 -> SetVisAttributes(G4VisAttributes::Invisible);
353
354 delete aBesG4Geo;
355 */
356 SetVisAttributes( "Pip_SCM" );
357 }
358
359 // construct subdetectors
360 if ( ReadBoostRoot::GetMdc() ) mdc->Construct( lWorld );
361 if ( ReadBoostRoot::GetTof() ) tof->Construct( lWorld );
362 if ( ReadBoostRoot::GetEmc() ) emc->Construct( lWorld );
363 if ( ReadBoostRoot::GetMuc() ) muc->Construct( lWorld );
364
365 G4cout << "get world " << fWorld->GetLogicalVolume()->GetNoDaughters() << endl;
366 for ( int nd = 0; nd < fWorld->GetLogicalVolume()->GetNoDaughters(); nd++ )
367 {
368 G4cout << "daughter n: name :" << nd << " "
369 << fWorld->GetLogicalVolume()->GetDaughter( nd )->GetName() << G4endl;
370 }
371
372 return fWorld;
373}
374
375G4LogicalVolume* BesDetectorConstruction::FindLogicalVolume( const G4String& vn ) {
376 return const_cast<G4LogicalVolume*>( GDMLProcessor::GetInstance()->GetLogicalVolume( vn ) );
377}
378
379void BesDetectorConstruction::SetVisAttributes( G4String volumeName ) {
380 if ( volumeName = "Pip_SCM" )
381 {
382 BesG4Geo* aBesG4Geo = new BesG4Geo();
383 fWorld = (G4VPhysicalVolume*)GDMLProcessor::GetInstance()->GetWorldVolume();
384 if ( fWorld == 0 )
385 {
386 G4cout
387 << "World volume not set properly check your setup selection criteria or GDML input!"
388 << G4endl;
389 exit( -1 );
390 }
391
392 // make the world volume invisible
393 lWorld = FindLogicalVolume( "logicalWorld" );
394 lWorld->SetVisAttributes( G4VisAttributes::Invisible );
395
396 // beam pipe
397
398 G4LogicalVolume* logicalPip = FindLogicalVolume( "logicalPip" );
399 G4LogicalVolume* logicalgoldLayer = FindLogicalVolume( "logicalgoldLayer" );
400 G4LogicalVolume* logicalinnerBe = FindLogicalVolume( "logicalinnerBe" );
401 G4LogicalVolume* logicaloilLayer = FindLogicalVolume( "logicaloilLayer" );
402 G4LogicalVolume* logicalouterBe = FindLogicalVolume( "logicalouterBe" );
403 G4LogicalVolume* logicalinnerBeSide = FindLogicalVolume( "logicalinnerBeSide" );
404 G4LogicalVolume* logicalinnerAl = FindLogicalVolume( "logicalinnerAl" );
405 G4LogicalVolume* logicalAlRing = FindLogicalVolume( "logicalAlRing" );
406 G4LogicalVolume* logicalouterAl = FindLogicalVolume( "logicalouterAl" );
407 G4LogicalVolume* logicalAgLayer = FindLogicalVolume( "logicalAgLayer" );
408 G4LogicalVolume* logicalCuLayer = FindLogicalVolume( "logicalCuLayer" );
409 G4LogicalVolume* logicalAlCover = FindLogicalVolume( "logicalAlCover" );
410
411 logicalPip->SetVisAttributes( G4VisAttributes::Invisible );
412 logicalgoldLayer->SetVisAttributes( G4VisAttributes::Invisible );
413 logicalinnerBe->SetVisAttributes( G4VisAttributes::Invisible );
414 logicaloilLayer->SetVisAttributes( G4VisAttributes::Invisible );
415 logicalouterBe->SetVisAttributes( G4VisAttributes::Invisible );
416 logicalinnerBeSide->SetVisAttributes( G4VisAttributes::Invisible );
417 logicalinnerAl->SetVisAttributes( G4VisAttributes::Invisible );
418 logicalAlRing->SetVisAttributes( G4VisAttributes::Invisible );
419 logicalouterAl->SetVisAttributes( G4VisAttributes::Invisible );
420 logicalAgLayer->SetVisAttributes( G4VisAttributes::Invisible );
421 logicalCuLayer->SetVisAttributes( G4VisAttributes::Invisible );
422 logicalAlCover->SetVisAttributes( G4VisAttributes::Invisible );
423
424 // SCM
425 G4LogicalVolume* logicalSCM = FindLogicalVolume( "logicalSCM" );
426 G4LogicalVolume* logicalinnCryo = FindLogicalVolume( "logicalinnCryo" );
427 G4LogicalVolume* logicalouterCryo = FindLogicalVolume( "logicalouterCryo" );
428 G4LogicalVolume* logicalendCryo = FindLogicalVolume( "logicalendCryo" );
429 G4LogicalVolume* logicall1Adia = FindLogicalVolume( "logicall1Adia" );
430 G4LogicalVolume* logicall2Adia = FindLogicalVolume( "logicall2Adia" );
431 G4LogicalVolume* logicall3Adia = FindLogicalVolume( "logicall3Adia" );
432 G4LogicalVolume* logicall4Adia = FindLogicalVolume( "logicall4Adia" );
433 G4LogicalVolume* logicalinnShield = FindLogicalVolume( "logicalinnShield" );
434 G4LogicalVolume* logicaloutShield = FindLogicalVolume( "logicaloutShield" );
435 G4LogicalVolume* logicalendShield = FindLogicalVolume( "logicalendShield" );
436 G4LogicalVolume* logicall1insu = FindLogicalVolume( "logicall1insu" );
437 G4LogicalVolume* logicall2insu = FindLogicalVolume( "logicall2insu" );
438 G4LogicalVolume* logicalcoil = FindLogicalVolume( "logicalcoil" );
439 G4LogicalVolume* logicalendCoil = FindLogicalVolume( "logicalendCoil" );
440 G4LogicalVolume* logicalsupp = FindLogicalVolume( "logicalsupp" );
441 G4LogicalVolume* logicalpipe1 = FindLogicalVolume( "logicalpipe1" );
442 G4LogicalVolume* logicalpipe2 = FindLogicalVolume( "logicalpipe2" );
443 G4LogicalVolume* logicalpipe3 = FindLogicalVolume( "logicalpipe3" );
444 G4LogicalVolume* logicaltub17 = FindLogicalVolume( "logicaltub17" );
445 G4LogicalVolume* logicaltub19 = FindLogicalVolume( "logicaltub19" );
446 G4LogicalVolume* logicaltub20 = FindLogicalVolume( "logicaltub20" );
447 G4LogicalVolume* logicaltub22 = FindLogicalVolume( "logicaltub22" );
448 G4LogicalVolume* logicaltub23 = FindLogicalVolume( "logicaltub23" );
449 G4LogicalVolume* logicaltub25 = FindLogicalVolume( "logicaltub25" );
450 G4LogicalVolume* logicalsub1 = FindLogicalVolume( "logicalsub1" );
451 G4LogicalVolume* logicalsub2 = FindLogicalVolume( "logicalsub2" );
452 G4LogicalVolume* logicalsub3 = FindLogicalVolume( "logicalsub3" );
453 G4LogicalVolume* logicalhole1 = FindLogicalVolume( "logicalhole1" );
454 G4LogicalVolume* logicalhole2 = FindLogicalVolume( "logicalhole2" );
455 G4LogicalVolume* logicalhole3 = FindLogicalVolume( "logicalhole3" );
456 G4LogicalVolume* logicalhole4 = FindLogicalVolume( "logicalhole4" );
457 G4LogicalVolume* logicalhole5 = FindLogicalVolume( "logicalhole5" );
458
459 // G4VisAttributes* visSCM1 = new G4VisAttributes(G4Colour(0.,0.,1.));
460 logicalSCM->SetVisAttributes( G4VisAttributes::Invisible );
461 // logicalinnCryo->SetVisAttributes(visSCM1);
462 logicalinnCryo->SetVisAttributes( G4VisAttributes::Invisible );
463 // logicalouterCryo->SetVisAttributes(visSCM1);
464 logicalouterCryo->SetVisAttributes( G4VisAttributes::Invisible );
465 // logicalendCryo->SetVisAttributes(visSCM1);
466 logicalendCryo->SetVisAttributes( G4VisAttributes::Invisible );
467 // G4VisAttributes* visSCM2 = new G4VisAttributes(G4Colour(0.,1.,0.));
468 // logicall1Adia->SetVisAttributes(visSCM2);
469 logicall1Adia->SetVisAttributes( G4VisAttributes::Invisible );
470 // logicall2Adia->SetVisAttributes(visSCM2);
471 logicall2Adia->SetVisAttributes( G4VisAttributes::Invisible );
472 // logicall3Adia->SetVisAttributes(visSCM2);
473 logicall3Adia->SetVisAttributes( G4VisAttributes::Invisible );
474 // logicall4Adia->SetVisAttributes(visSCM2);
475 logicall4Adia->SetVisAttributes( G4VisAttributes::Invisible );
476 // G4VisAttributes* visSCM3 = new G4VisAttributes(G4Colour(1.,0.,0.));
477 // logicalinnShield->SetVisAttributes(visSCM3);
478 logicalinnShield->SetVisAttributes( G4VisAttributes::Invisible );
479 // logicaloutShield->SetVisAttributes(visSCM3);
480 logicaloutShield->SetVisAttributes( G4VisAttributes::Invisible );
481 // logicalendShield->SetVisAttributes(visSCM3);
482 logicalendShield->SetVisAttributes( G4VisAttributes::Invisible );
483
484 // G4VisAttributes* visSCM4 = new G4VisAttributes(G4Colour(0.,1.,1.));
485 // logicall1insu->SetVisAttributes(visSCM4);
486 logicall1insu->SetVisAttributes( G4VisAttributes::Invisible );
487 // logicall2insu->SetVisAttributes(visSCM4);
488 logicall2insu->SetVisAttributes( G4VisAttributes::Invisible );
489 // G4VisAttributes* visSCM5 = new G4VisAttributes(G4Colour(1.,0.,1.));
490 // logicalcoil->SetVisAttributes(visSCM5);
491 logicalcoil->SetVisAttributes( G4VisAttributes::Invisible );
492 // logicalendCoil->SetVisAttributes(visSCM5);
493 logicalendCoil->SetVisAttributes( G4VisAttributes::Invisible );
494 // G4VisAttributes* visSCM6 = new G4VisAttributes(G4Colour(1.,1.,0.));
495 // logicalsupp->SetVisAttributes(visSCM6);
496 logicalsupp->SetVisAttributes( G4VisAttributes::Invisible );
497 // G4VisAttributes* visSCM7 = new G4VisAttributes(G4Colour(0.,0.,0.));
498 // logicalrein->SetVisAttributes(G4VisAttributes::Invisible);
499 logicalpipe1->SetVisAttributes( G4VisAttributes::Invisible );
500 // logicalpipe2->SetVisAttributes(visSCM2);
501 logicalpipe2->SetVisAttributes( G4VisAttributes::Invisible );
502 logicalpipe3->SetVisAttributes( G4VisAttributes::Invisible );
503 // logicaltub17->SetVisAttributes(visSCM2);
504 logicaltub17->SetVisAttributes( G4VisAttributes::Invisible );
505 // logicaltub19->SetVisAttributes(visSCM5);
506 logicaltub19->SetVisAttributes( G4VisAttributes::Invisible );
507 // logicaltub20->SetVisAttributes(visSCM1);
508 logicaltub20->SetVisAttributes( G4VisAttributes::Invisible );
509 // logicaltub22->SetVisAttributes(visSCM1);
510 logicaltub22->SetVisAttributes( G4VisAttributes::Invisible );
511 // logicaltub23->SetVisAttributes(visSCM2);
512 logicaltub23->SetVisAttributes( G4VisAttributes::Invisible );
513 // logicaltub25->SetVisAttributes(visSCM5);
514 logicaltub25->SetVisAttributes( G4VisAttributes::Invisible );
515 // logicalsub1->SetVisAttributes(visSCM4);
516 logicalsub1->SetVisAttributes( G4VisAttributes::Invisible );
517 // logicalsub2->SetVisAttributes(visSCM4);
518 logicalsub2->SetVisAttributes( G4VisAttributes::Invisible );
519
520 // logicalsub3->SetVisAttributes(visSCM4);
521 logicalsub3->SetVisAttributes( G4VisAttributes::Invisible );
522 // logicalhole1->SetVisAttributes(visSCM3);
523 logicalhole1->SetVisAttributes( G4VisAttributes::Invisible );
524 // logicalhole2->SetVisAttributes(visSCM3);
525 logicalhole2->SetVisAttributes( G4VisAttributes::Invisible );
526 // logicalhole3->SetVisAttributes(visSCM3);
527 logicalhole3->SetVisAttributes( G4VisAttributes::Invisible );
528 // logicalhole4->SetVisAttributes(visSCM3);
529 logicalhole4->SetVisAttributes( G4VisAttributes::Invisible );
530 // logicalhole5->SetVisAttributes(visSCM3);
531 logicalhole5->SetVisAttributes( G4VisAttributes::Invisible );
532
533 delete aBesG4Geo;
534 }
535}
n
const Int_t n
Definition DataBase/tau_mode.c:57
I
const DifComplex I
C
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA all others are from PDG *Only resonances with known widths into electron pairs are sept ************************************************************************C Declarations C
Definition RRes.h:29
BesDetectorConstruction::FindLogicalVolume
G4LogicalVolume * FindLogicalVolume(const G4String &vn)
Definition BesDetectorConstruction.cc:375
BesDetectorConstruction::Construct
G4VPhysicalVolume * Construct()
Definition BesDetectorConstruction.cc:61
BesDetectorConstruction::SetVisAttributes
void SetVisAttributes(G4String volumeName)
Definition BesDetectorConstruction.cc:379
H
IMPLICIT REAL *A H
Definition myXsection.h:1