68{
69
70
71
73 G4double z;
74 G4double density;
75
77 new G4Material(
"ArgonGas", z= 18.,
a= 39.95*g/mole, density= 1.782*mg/cm3);
78
80 new G4Material(
"Aluminum", z= 13.,
a= 26.98*g/mole, density= 2.7*g/cm3);
81
83 new G4Material(
"Lead", z= 82.,
a= 207.19*g/mole, density= 11.35*g/cm3);
84
85
86
87
88
89
90 G4double expHall_x = 3.0*m;
91 G4double expHall_y = 1.0*m;
92 G4double expHall_z = 1.0*m;
93 G4Box* experimentalHall_box
94 =
new G4Box(
"expHall_box",expHall_x,expHall_y,expHall_z);
96 Ar,"expHall_log",0,0,0);
99
100
101
102 G4double innerRadiusOfTheTube = 0.*cm;
103 G4double outerRadiusOfTheTube = 60.*cm;
104 G4double hightOfTheTube = 50.*cm;
105 G4double startAngleOfTheTube = 0.*deg;
106 G4double spanningAngleOfTheTube = 360.*deg;
107 G4Tubs* tracker_tube =
new G4Tubs(
"tracker_tube",innerRadiusOfTheTube,
108 outerRadiusOfTheTube,hightOfTheTube,
109 startAngleOfTheTube,spanningAngleOfTheTube);
111 G4double trackerPos_x = -1.0*m;
112 G4double trackerPos_y = 0.*m;
113 G4double trackerPos_z = 0.*m;
115 G4ThreeVector(trackerPos_x,trackerPos_y,trackerPos_z),
117
118
119
120 G4double block_x = 1.0*m;
121 G4double block_y = 50.0*cm;
122 G4double block_z = 50.0*cm;
123 G4Box* calorimeterBlock_box =
new G4Box(
"calBlock_box",block_x,
124 block_y,block_z);
126 Pb,"caloBlock_log",0,0,0);
127 G4double blockPos_x = 1.0*m;
128 G4double blockPos_y = 0.0*m;
129 G4double blockPos_z = 0.0*m;
131 G4ThreeVector(blockPos_x,blockPos_y,blockPos_z),
133
134
135
136 G4double calo_x = 1.*cm;
137 G4double calo_y = 40.*cm;
138 G4double calo_z = 40.*cm;
139 G4Box* calorimeterLayer_box =
new G4Box(
"caloLayer_box",
140 calo_x,calo_y,calo_z);
142 Al,"caloLayer_log",0,0,0);
143 for(G4int i=0;i<19;i++)
144 {
145 G4double caloPos_x = (i-9)*10.*cm;
146 G4double caloPos_y = 0.0*m;
147 G4double caloPos_z = 0.0*m;
149 G4ThreeVector(caloPos_x,caloPos_y,caloPos_z),
151 }
152
153
154
155
156
157
158 TSystem ts;
159
160 gSystem->Load("libExP02ClassesDict");
161
162
163
164 const G4ElementTable* eltab = G4Element::GetElementTable();
165
167
168 TFile fo("geo.root","RECREATE");
169
170 fo.WriteObject(geotree, "my_geo");
171
173}
std::vector< ExP01TrackerHit * > a
Helper class needed for the ROOT I/O.