BOSS 7.0.1
BESIII Offline Software System
Loading...
Searching...
No Matches
BesTofConstruction.cc
Go to the documentation of this file.
1//---------------------------------------------------------------------------//
2// BOOST --- BESIII Object_Oriented Simulation Tool //
3//---------------------------------------------------------------------------//
4//Description: Geometry Construction for barrel and endcap Tof
5//Author: Dengzy
6//Created: Mar, 2004
7//Modified: Way of reading txt file in dat/ is changed.
8//The material for BC404 and BC408 has been exchanged. The original
9//one is opposite. (Liu Yong @ Dec 2008)
10//$Id: BesTofConstruction.cc
11//
12
13#include "GaudiKernel/MsgStream.h"
14#include "GaudiKernel/Bootstrap.h"
15#include "GaudiKernel/PropertyMgr.h"
16#include "GaudiKernel/IJobOptionsSvc.h"
17#include "GaudiKernel/ISvcLocator.h"
18#include "GaudiKernel/IDataProviderSvc.h"
19
20#include "BesTofConstruction.hh"
21#include "BesTofGeoParameter.hh"
22#include "PartProduce.hh"
23#include "BesTofSD.hh"
24#include "G4Geo/TofG4Geo.h"
25#include "G4Geo/MRPCG4Geo.h"
26#include "ReadBoostRoot.hh"
27
28#include "G4Material.hh"
29#include "G4Trap.hh"
30#include "G4Tubs.hh"
31#include "G4Box.hh"
32#include "G4Trd.hh"
33#include "G4Cons.hh"
34#include "G4UnionSolid.hh"
35#include "G4SubtractionSolid.hh"
36#include "G4IntersectionSolid.hh"
37
38#include "G4LogicalVolume.hh"
39#include "G4Transform3D.hh"
40#include "G4PVPlacement.hh"
41#include "G4AssemblyVolume.hh"
42#include "G4SDManager.hh"
43#include "G4UnitsTable.hh"
44#include "G4ios.hh"
45#include "G4Colour.hh"
46#include "G4VisAttributes.hh"
47#include "G4UnionSolid.hh"
48#include "G4UserLimits.hh"
49
50
51
52BesTofConstruction::BesTofConstruction() //Constructor, set variables (G4logical or G4physical volumes) declared in BesTofConstruction.hh
53{
54 PropertyMgr m_propMgr1;
55 m_propMgr1.declareProperty("UserLimits", m_userLimits = 0.03);
56 m_propMgr1.declareProperty("IonE", m_ionE = 100);
57 IJobOptionsSvc* jobSvc1;
58 Gaudi::svcLocator()->service("JobOptionsSvc", jobSvc1);
59 jobSvc1->setMyProperties("BesTofConstruction", &m_propMgr1);
60
61 logicalTof = 0;
62 physicalTof = 0;
63
64 logicalBrTof = 0;
65 logicalEcTofWest=0;
66 logicalEcTofEast=0;
67
68 logicalScinBr1 = 0;
69 logicalScinBr2 = 0;
70 physicalScinBr1 = 0;
71 physicalScinBr2 = 0;
72
73 logicalAlBr1 = 0;
74 logicalAlBr2 = 0;
75 physicalAlBr1 = 0;
76 physicalAlBr2 = 0;
77
78 logicalPVFBr1 = 0;
79 logicalPVFBr2 = 0;
80 physicalPVFBr1 = 0;
81 physicalPVFBr2 = 0;
82
83 logicalBucketBr1= 0;
84 logicalBucketBr2= 0;
85
86 physicalBucket1Br1 = 0;
87 physicalBucket2Br1 = 0;
88 physicalBucket1Br2 = 0;
89 physicalBucket2Br2 = 0;
90
91 logicalScinEcWest = 0;
92 logicalScinEcEast = 0;
93 physicalScinEcWest = 0;
94 physicalScinEcEast = 0;
95
96 logicalAlEcWest = 0;
97 logicalAlEcEast = 0;
98 physicalAlEcWest = 0;
99 physicalAlEcEast = 0;
100
101 logicalPVFEcWest = 0;
102 logicalPVFEcEast = 0;
103 physicalPVFEcWest = 0;
104 physicalPVFEcEast = 0;
105
106 logicalBucketEc = 0;
107 physicalBucketEc = 0;
108
109 BC404 = 0;
110 BC408 = 0;
111 PVF = 0;
112 PMTmaterial = 0;
113
114 //MRPC
115 logical_gasLayer=0;
116}
117
118BesTofConstruction::~BesTofConstruction()
119{}
120
121void BesTofConstruction::Construct(G4LogicalVolume* logicalBes)
122{
123 DefineMaterial(); //This function just defines my material, used for the TOF
124
125 if ( ReadBoostRoot::GetTof() == 2 )
126 {
127 TofG4Geo* aTofG4Geo = new TofG4Geo();
128 logicalTof = aTofG4Geo->GetTopVolume();
129 if (!logicalTof)
130 G4cout<<"BesTofConstruction::Construct(), logicalTof not found"<<G4endl;
131 else
132 {
133 //produce the logical Tof
134 physicalTof = new G4PVPlacement(0,G4ThreeVector(0,0,0),
135 logicalTof,"physicalTof",logicalBes,false,0); // Class representing a single volume positioned within and relative to a mother volume.
136
137 //for visual attributes
138 logicalScinBr1 = FindLogicalVolume("logicalScinBr1");
139 logicalScinBr2 = FindLogicalVolume("logicalScinBr2");
140
141 logicalAlBr1 = FindLogicalVolume("logicalAlBr1");
142 logicalAlBr2 = FindLogicalVolume("logicalAlBr2");
143 logicalPVFBr1 = FindLogicalVolume("logicalPVFBr1");
144 logicalPVFBr2 = FindLogicalVolume("logicalPVFBr2");
145 logicalBucketBr1 = FindLogicalVolume("logicalBucketBr1");
146 logicalBucketBr2 = FindLogicalVolume("logicalBucketBr2");
147
148 logicalScinEcWest = FindLogicalVolume("logicalScinEcWest");
149 logicalScinEcEast = FindLogicalVolume("logicalScinEcEast");
150 logicalBucketEc = FindLogicalVolume("logicalBucketEc");
151 }
152 //m_sxp.Finalize();
153 delete aTofG4Geo;
154 } //close (ReadBoostRoot::GetTof() == 2)
155
156 else if(ReadBoostRoot::GetTof() == 3 ) //start TOF with MRPC Endcaps
157 {
158 //This part intialize a function where the Parameter of the TOF are saved!
159 BesTofGeoParameter* tofPara = BesTofGeoParameter::GetInstance();
160
161 //Get the coordinates of the TOF
162 G4double r1=380; //Newly changed from 399 for MRPC to avoid overlap
163 G4double r2=tofPara->GetBucketPosR()+0.5*tofPara->GetBucketDEc()+1; //=474
164 G4double r3=810;
165 G4double r4=925;
166 G4double a1=1382+tofPara->GetBucketLEc()+1; //1463
167 G4double a2=1381; // Changed from 1382 to avoid overlapping with EMC
168 //to make the center of tub6 is tofPara.GetzPosEastEc()=1356
169 G4double a3=1330;
170 G4cout<<"Tof Volume: "<<r1<<" "<<r2<<" "<<r3<<" "<<r4<<" "<<a1<<" "<<a2<<" "<<a3<<G4endl;
171
172 //Tof: produce a cylinder with inner radius r1 and outer radius r2, z height: 2a1, starting angle 0 untill 360
173 G4Tubs* tub1 = new G4Tubs("tub1",r1,r2,a1,0,360);
174 G4Tubs* tub2 = new G4Tubs("tub2",r2-4,r4,a2,0,360);
175 G4Tubs* tub3 = new G4Tubs("tub3",0,r3,a3,0,360);
176 G4UnionSolid* tub4 = new G4UnionSolid("tub4",tub1,tub2,0,G4ThreeVector(0,0,0) ); //Vereinigt tub 1 und 2 in tub4
177 G4SubtractionSolid* solidTof = new G4SubtractionSolid("solidTof",tub4,tub3,0, G4ThreeVector(0,0,0) ); //Subtracts tub4 - tub3
178 logicalTof = new G4LogicalVolume(solidTof, G4Material::GetMaterial("Air"),"logicalTof"); //book logical Tof VOlume
179 physicalTof = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalTof,"physicalTof",logicalBes,false,0);
180
181 //Logical + Physical volumes for the barrel
182 G4Tubs* tub5 = new G4Tubs("tub5",r3,r4,a3,0,360);
183 logicalBrTof = new G4LogicalVolume(tub5, G4Material::GetMaterial("Air"),"logicalBrTof"); //book logical Barrel TOF
184 G4VPhysicalVolume* physicalBrTof = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalBrTof,"physicalBrTof",logicalTof,false,0);
185
186
187 //Logical + Physical volumes for the endcaps
188 G4Tubs* tub6 = new G4Tubs("tub6",r1,r4,(a2-a3)/2,0,360); //ring: dz/2=2.55cm, xy=52.6 cm
189 G4Tubs* tub7 = new G4Tubs("tub7",r1,r2,(a1-a2)/2,0,360); //ring: dz/2=0.1cm, xy=7.5 cm
190 G4UnionSolid* tub8 = new G4UnionSolid("tub8",tub6,tub7,0,G4ThreeVector(0,0,(a3-a1)/2)); //Unify both tubes,translate z-axes -6.15 cm to the originally one
191 logicalEcTofWest = new G4LogicalVolume(tub8, G4Material::GetMaterial("Air"),"logicalEcTofWest");
192 G4VPhysicalVolume* physicalEcTofWest = new G4PVPlacement(0,G4ThreeVector(0,0,tofPara->GetzPosWestEc()),logicalEcTofWest,"physicalEcTofWest",logicalTof,false,0);
193
194 G4UnionSolid* tub9 = new G4UnionSolid("tub9",tub6,tub7,0,G4ThreeVector(0,0,(a1-a3)/2) ); //Unify both tubes,translate z-axes +6.15 cm to the originally one
195 logicalEcTofEast = new G4LogicalVolume(tub9, G4Material::GetMaterial("Air"),"logicalEcTofEast");
196 G4VPhysicalVolume* physicalEcTofEast = new G4PVPlacement(0,G4ThreeVector(0,0,tofPara->GetzPosEastEc()),logicalEcTofEast,"physicalEcTofEast",logicalTof,false,0);
197
198
199 //Functions which construct the different parts of the TOF -->See below
200 //Construct Tape
201 //put these lines before ConstructBr1Tof();
202 //so in BesTofSD.cc, no need to change the tofid caculation method
203 G4Tubs* tubTape = new G4Tubs("tubTape",866, 866.3, 1150, 0, 360);
204 G4LogicalVolume* logicalTape = new G4LogicalVolume(tubTape,G4Material::GetMaterial("tape"),"logicalTape");
205 G4VPhysicalVolume* physicalTape = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalTape,"physicalTape",logicalBrTof,false,0);
206 logicalTape->SetVisAttributes(G4VisAttributes::Invisible);
207
208 ConstructBr1Tof();
209 ConstructBr2Tof();
210 ConstructEcTof_mrpc();
211 }
212
213 else if (ReadBoostRoot::GetTof() == 4 ) //start TOF with MRPC Endcaps GDML Construction
214 {
215 MRPCG4Geo* aTofG4Geo = new MRPCG4Geo();
216 logicalTof = aTofG4Geo->GetTopVolume();
217 if (!logicalTof)
218 G4cout<<"BesTofConstruction::Construct(), logicalTof not found"<<G4endl;
219 else
220 {
221 physicalTof = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalTof,"physicalTof",logicalBes,false,0);
222
223 logicalScinBr1 = FindLogicalVolume("logicalScinBr1");
224 logicalScinBr2 = FindLogicalVolume("logicalScinBr2");
225
226 logicalAlBr1 = FindLogicalVolume("logicalAlBr1");
227 logicalAlBr2 = FindLogicalVolume("logicalAlBr2");
228 logicalPVFBr1 = FindLogicalVolume("logicalPVFBr1");
229 logicalPVFBr2 = FindLogicalVolume("logicalPVFBr2");
230 logicalBucketBr1 = FindLogicalVolume("logicalBucketBr1");
231 logicalBucketBr2 = FindLogicalVolume("logicalBucketBr2");
232
233 logical_gasLayer=FindLogicalVolume("logical_gasLayer");
234 }
235 delete aTofG4Geo;
236 } //close else if(==4)
237
238 else
239 {
240 //-----------------logicalTof added here
241 BesTofGeoParameter* tofPara = BesTofGeoParameter::GetInstance();
242
243 //Get the coordinates of the TOF
244 G4double r1=tofPara->GetEcR1()-1; //=399
245 G4double r2=tofPara->GetBucketPosR()+0.5*tofPara->GetBucketDEc()+1; //=474
246 G4double r3=810;;
247 G4double r4=925;
248 G4double a1=1382+tofPara->GetBucketLEc()+1; //1463
249 G4double a2=1381; // Changed from 1382 to avoid overlapping with EMC
250 //to make the center of tub6 is tofPara.GetzPosEastEc()=1356
251 G4double a3=1330;
252 G4cout<<"Tof Volume: "<<r1<<" "<<r2<<" "<<r3<<" "<<r4<<" "<<a1<<" "<<a2<<" "<<a3<<G4endl;
253
254 G4Tubs* tub1 = new G4Tubs("tub1",r1,r2,a1,0,360); //Just produce a cylinder with inner radius r1 and outer radius r2, z height: 2a1, starting angle 0 untill 360
255 G4Tubs* tub2 = new G4Tubs("tub2",r2-4,r4,a2,0,360);
256 G4Tubs* tub3 = new G4Tubs("tub3",0,r3,a3,0,360);
257 G4UnionSolid* tub4 = new G4UnionSolid("tub4",tub1,tub2,0,G4ThreeVector(0,0,0) ); //Vereinigt tub 1 und 2 in tub4
258 G4SubtractionSolid* solidTof = new G4SubtractionSolid("solidTof",tub4,tub3,0, G4ThreeVector(0,0,0) ); //Subtracts tub4 - tub3
259 logicalTof = new G4LogicalVolume(solidTof, G4Material::GetMaterial("Air"),"logicalTof"); //book logical Tof VOlume
260 physicalTof = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalTof,"physicalTof",logicalBes,false,0);
261
262 G4Tubs* tub5 = new G4Tubs("tub5",r3,r4,a3,0,360);
263 logicalBrTof = new G4LogicalVolume(tub5, G4Material::GetMaterial("Air"),"logicalBrTof"); //book logical Barrel TOF
264 G4VPhysicalVolume* physicalBrTof = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalBrTof,"physicalBrTof",logicalTof,false,0);
265
266 //Logical + Physical volumes for the endcaps
267 G4Tubs* tub6 = new G4Tubs("tub6",r1,r4,(a2-a3)/2,0,360);//r1=39.9cm, r4=92.5cm, zlength=2.55 cm
268 G4Tubs* tub7 = new G4Tubs("tub7",r1,r2,(a1-a2)/2,0,360);//r1=39.9cm, r2=47.4cm, zlength=4.05 cm
269 G4UnionSolid* tub8 = new G4UnionSolid("tub8",tub6,tub7,0,G4ThreeVector(0,0,(a3-a1)/2) );// -6.65cm
270 logicalEcTofWest = new G4LogicalVolume(tub8, G4Material::GetMaterial("Air"),"logicalEcTofWest");
271 G4VPhysicalVolume* physicalEcTofWest = new G4PVPlacement(0,G4ThreeVector(0,0,tofPara->GetzPosWestEc()),logicalEcTofWest,"physicalEcTofWest",logicalTof,false,0);
272
273 G4UnionSolid* tub9 = new G4UnionSolid("tub9",tub6,tub7,0,G4ThreeVector(0,0,(a1-a3)/2) );
274 logicalEcTofEast = new G4LogicalVolume(tub9, G4Material::GetMaterial("Air"),"logicalEcTofEast");
275 G4VPhysicalVolume* physicalEcTofEast = new G4PVPlacement(0,G4ThreeVector(0,0,tofPara->GetzPosEastEc()),logicalEcTofEast,"physicalEcTofEast",logicalTof,false,0);
276
277
278 //Construct Tape
279 //put these lines before ConstructBr1Tof();
280 //so in BesTofSD.cc, no need to change the tofid caculation method
281 G4Tubs* tubTape = new G4Tubs("tubTape",866, 866.3, 1150, 0, 360);
282 G4LogicalVolume* logicalTape = new G4LogicalVolume(tubTape,G4Material::GetMaterial("tape"),"logicalTape");
283 G4VPhysicalVolume* physicalTape = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalTape,"physicalTape",logicalBrTof,false,0);
284 logicalTape->SetVisAttributes(G4VisAttributes::Invisible);
285
286 //Functions which construct the different parts of the TOF -->See below
287 ConstructBr1Tof();
288 ConstructBr2Tof();
289 ConstructEcTof();
290 }
291
292
293 //sensitive detector
294 G4SDManager* SDman = G4SDManager::GetSDMpointer();
295 BesTofSD* tofSD = new BesTofSD("BesTofSD");
296 SDman->AddNewDetector(tofSD);
297
298 if (logicalScinBr1) logicalScinBr1->SetSensitiveDetector(tofSD);
299 if (logicalScinBr2) logicalScinBr2->SetSensitiveDetector(tofSD);
300 if (logicalScinEcWest) logicalScinEcWest->SetSensitiveDetector(tofSD);
301 if (logicalScinEcEast) logicalScinEcEast->SetSensitiveDetector(tofSD);
302 if(logical_gasLayer)
303 {
304 logical_gasLayer->SetSensitiveDetector(tofSD);
305 G4cout<<"!!! Find logical_gasLayer !!!"<<G4endl;
306
307 //For limits smaller than 0.02 one has to use the PAI-Model for ionisation!
308 logical_gasLayer->SetUserLimits (new G4UserLimits(m_userLimits*mm));
309 //logical_gasLayer->SetUserLimits (new G4UserLimits(0.03*mm));
310 }
311
312
313 G4VisAttributes* visBrTof = new G4VisAttributes(G4Colour(1.,0.,0.));
314 G4VisAttributes* visEcTof = new G4VisAttributes(G4Colour(0.,1.,0.));
315
316 if (logicalTof)
317 logicalTof->SetVisAttributes(G4VisAttributes::Invisible);
318 if (logicalBrTof)
319 logicalBrTof->SetVisAttributes(G4VisAttributes::Invisible);
320 //logicalBrTof->SetVisAttributes(visBrTof);
321 if (logicalEcTofWest)
322 logicalEcTofWest->SetVisAttributes(G4VisAttributes::Invisible);
323 //logicalEcTofWest->SetVisAttributes(visEcTof);
324 if (logicalEcTofEast)
325 logicalEcTofEast->SetVisAttributes(G4VisAttributes::Invisible);
326 //logicalEcTofEast->SetVisAttributes(visEcTof);
327
328
329 if (logicalScinBr1)
330 logicalScinBr1->SetVisAttributes(G4VisAttributes::Invisible);
331 if (logicalScinBr2)
332 logicalScinBr2->SetVisAttributes(G4VisAttributes::Invisible);
333
334 if (logicalAlBr1)
335 logicalAlBr1->SetVisAttributes(G4VisAttributes::Invisible);
336 if (logicalAlBr2)
337 logicalAlBr2->SetVisAttributes(G4VisAttributes::Invisible);
338
339 if (logicalPVFBr1)
340 logicalPVFBr1->SetVisAttributes(G4VisAttributes::Invisible);
341 if (logicalPVFBr2)
342 logicalPVFBr2->SetVisAttributes(G4VisAttributes::Invisible);
343
344 G4VisAttributes* visAttBrBuck = new G4VisAttributes(G4Colour(1.,1.,0.));
345 if (logicalBucketBr1)
346 logicalBucketBr1->SetVisAttributes(visAttBrBuck);
347 //logicalBucketBr1->SetVisAttributes(G4VisAttributes::Invisible);
348 if (logicalBucketBr2)
349 logicalBucketBr2->SetVisAttributes(visAttBrBuck);
350 //logicalBucketBr2->SetVisAttributes(G4VisAttributes::Invisible);
351
352 G4VisAttributes* visAttEcTof = new G4VisAttributes(G4Colour(0.,1.,1.));
353 if (logicalScinEcWest)
354 logicalScinEcWest->SetVisAttributes(visAttEcTof);
355 //logicalScinEc->SetVisAttributes(G4VisAttributes::Invisible);
356 if (logicalScinEcEast)
357 logicalScinEcEast->SetVisAttributes(visAttEcTof);
358
359 G4VisAttributes* visAttEcBuck = new G4VisAttributes(G4Colour(1.,1.,0.));
360 if (logicalBucketEc)
361 logicalBucketEc->SetVisAttributes(visAttEcBuck);
362 //logicalBucketEc->SetVisAttributes(G4VisAttributes::Invisible);
363
364 if (logicalAlEcWest)
365 logicalAlEcWest->SetVisAttributes(G4VisAttributes::Invisible);
366 if (logicalAlEcEast)
367 logicalAlEcEast->SetVisAttributes(G4VisAttributes::Invisible);
368 if (logicalPVFEcWest)
369 logicalPVFEcWest->SetVisAttributes(G4VisAttributes::Invisible);
370 if (logicalPVFEcEast)
371 logicalPVFEcEast->SetVisAttributes(G4VisAttributes::Invisible);
372}
373
374void BesTofConstruction::ConstructBr1Tof()
375{
376 BesTofGeoParameter* tofPara = BesTofGeoParameter::GetInstance();
377
378 //Get barrel tof layer1 geometry data
379 G4int nScinBr=tofPara->GetnScinBr(); //number of barrel scintillators
380 G4double br1L=tofPara->GetBr1L();
381 G4double br1TrapW1=tofPara->GetBr1TrapW1();
382 G4double br1TrapW2=tofPara->GetBr1TrapW2();
383 G4double br1TrapH=tofPara->GetBr1TrapH();
384 G4double br1R1=tofPara->GetBr1R1();
385 G4double AlThickness=tofPara->GetAlThickness();
386 G4double PVFThickness=tofPara->GetPVFThickness();
387 //barrel PMT bucket geometry data
388 G4double bucketDBr=tofPara->GetBucketDBr(); //diameter of barrel PMT bucket
389 G4double bucketLBr=tofPara->GetBucketLBr(); //length of barrel PMT bucket
390
391 //computer from the original data
392 G4double angleBr=360./nScinBr*deg;
393 G4double scinTrapdx1,scinTrapdx2,scinTrapdx3,scinTrapdx4;
394 scinTrapdx1=scinTrapdx3=br1TrapW1/2.;
395 scinTrapdx2=scinTrapdx4=br1TrapW2/2.;
396 G4double scinTrapdy1,scinTrapdy2;
397 scinTrapdy1=scinTrapdy2=0.5*br1TrapH;
398 G4double scinTrapdz=br1L/2.;
399 G4double scinPosR=br1R1+scinTrapdy1;
400 G4double theta = atan((br1TrapW2/2.-br1TrapW1/2.)/br1TrapH);
401 G4double delta13 = AlThickness*(1/cos(theta)-tan(theta));
402 G4double delta24 = AlThickness*(1/cos(theta)+tan(theta));
403 G4double delta13P = (PVFThickness+AlThickness)*(1/cos(theta)-tan(theta));
404 G4double delta24P = (PVFThickness+AlThickness)*(1/cos(theta)+tan(theta));
405
406 //construct barrel tof scintillator
407 G4Trap* solidScinBr1=new G4Trap("solidScinBr1",scinTrapdz,0*deg,0*deg,
408 scinTrapdy1,scinTrapdx1,scinTrapdx2,0*deg,
409 scinTrapdy2,scinTrapdx3,scinTrapdx4,0*deg);
410 logicalScinBr1 = new G4LogicalVolume(solidScinBr1,BC408,"logicalScinBr1");
411
412 //construct barrel Al film and PVF film
413 G4Trap* solidAlBr1 = new G4Trap("solidAlBr1", scinTrapdz+0.001,0*deg,0*deg,
414 scinTrapdy1+AlThickness, scinTrapdx1+delta13, scinTrapdx2+delta24 ,0*deg,
415 scinTrapdy2+AlThickness, scinTrapdx3+delta13, scinTrapdx4+delta24 ,0*deg);
416
417 logicalAlBr1 = new G4LogicalVolume(solidAlBr1,G4Material::GetMaterial("Aluminium"),"logicalAlBr1");
418 G4Trap* solidPVFBr1 = new G4Trap("solidPVFBr1",scinTrapdz+0.002,0*deg,0*deg,
419 scinTrapdy1+AlThickness+PVFThickness,scinTrapdx1+delta13P, scinTrapdx2+delta24P,0*deg,
420 scinTrapdy2+AlThickness+PVFThickness,scinTrapdx3+delta13P, scinTrapdx4+delta24P,0*deg);
421
422 logicalPVFBr1 = new G4LogicalVolume(solidPVFBr1,PVF,"logicalPVFBr1");
423 //put daughter in mother logical volume
424 physicalAlBr1 = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalAlBr1, "physicalAlBr1",logicalPVFBr1, false,0);
425 physicalScinBr1 = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalScinBr1, "physicalScinBr1",logicalAlBr1, false, 0);
426
427 //construct barrel PMT bucket (bucket=Behaelter)
428 G4Tubs* solidBucketBr1= new G4Tubs("solidBucketBr1",0,bucketDBr/2.,
429 bucketLBr/2., 0, 360.*deg);
430 logicalBucketBr1=new G4LogicalVolume(solidBucketBr1,PMTmaterial,"logicalBucketBr1");
431
432
433 //Rotate the crystalls:
434 G4RotationMatrix trapRot;
435 trapRot.rotateZ(0.5*angleBr-90.*deg);
436 // /| y
437 //------- / | |
438 //\ / | | rotate->(0.5*angleBr-90.*deg): |
439 // \ / ->>>> | | -----x
440 // --- \ |
441 // \|
442 G4ThreeVector pos(0,0,0);
443 G4double x,y;
444 for (int i=0;i<nScinBr;i++)
445 {
446 ostringstream osnamePVFBr1;
447 osnamePVFBr1 << "physicalPVFBr1_"<<i;
448 ostringstream osnameBucket1Br1;
449 osnameBucket1Br1 << "physicalBucket1Br1_"<<i;
450 ostringstream osnameBucket2Br1;
451 osnameBucket2Br1 << "physicalBucket2Br1_"<<i;
452
453 //scintillator and casing
454 x=scinPosR*cos((i+0.5)*angleBr); //+0.5: start from phi=0
455 y=scinPosR*sin((i+0.5)*angleBr);
456 pos.setX(x);
457 pos.setY(y);
458 pos.setZ(0);
459
460 //to make code consistent with gdml
461 double a, b, c;
462 getXYZ(trapRot, a, b ,c); //Function is defined below
463 G4RotationMatrix trapRotTemp;
464 trapRotTemp.rotateZ(c);
465 G4Transform3D transform(trapRotTemp,pos);
466
467 //G4PVPlacement: Class representing a single volume positioned within and relative to a mother volume.
468 //G4Transform3D transform(trapRot,pos);
469 physicalPVFBr1 = new G4PVPlacement(transform,logicalPVFBr1,osnamePVFBr1.str(), logicalBrTof,false,i);
470 //bucket --> Two PMT, one on each side
471 pos.setZ(-(bucketLBr+br1L)/2.-0.01);
472 physicalBucket1Br1 = new G4PVPlacement(0,pos,logicalBucketBr1,osnameBucket1Br1.str(),logicalBrTof,false,0);
473 pos.setZ((bucketLBr+br1L)/2.+0.01);
474 physicalBucket2Br1 = new G4PVPlacement(0,pos,logicalBucketBr1,osnameBucket2Br1.str(),logicalBrTof,false,0);
475
476 //rotate for next scintillator
477 trapRot.rotateZ(angleBr);
478 }
479}
480
481void BesTofConstruction::ConstructBr2Tof()
482{
483 BesTofGeoParameter* tofPara = BesTofGeoParameter::GetInstance();
484
485 //barrel tof layer2 geometry data
486 G4int nScinBr=tofPara->GetnScinBr(); //number of barrel scintillators
487 G4double br2L=tofPara->GetBr2L();
488 G4double br2TrapW1=tofPara->GetBr2TrapW1();
489 G4double br2TrapW2=tofPara->GetBr2TrapW2();
490 G4double br2TrapH=tofPara->GetBr2TrapH();
491 G4double br2R1=tofPara->GetBr2R1();
492 G4double AlThickness=tofPara->GetAlThickness();
493 G4double PVFThickness=tofPara->GetPVFThickness();
494 //barrel PMT bucket geometry data
495 G4double bucketDBr=tofPara->GetBucketDBr(); //diameter of barrel PMT bucket
496 G4double bucketLBr=tofPara->GetBucketLBr(); //length of barrel PMT bucket
497
498 //computer from the original data
499 G4double angleBr=360./nScinBr*deg;
500 G4double scinTrapdx1,scinTrapdx2,scinTrapdx3,scinTrapdx4;
501 scinTrapdx1=scinTrapdx3=br2TrapW1/2.;
502 scinTrapdx2=scinTrapdx4=br2TrapW2/2.;
503 G4double scinTrapdy1,scinTrapdy2;
504 scinTrapdy1=scinTrapdy2=0.5*br2TrapH;
505 G4double scinTrapdz=br2L/2.;
506 G4double scinPosR=br2R1+scinTrapdy1;
507 G4double theta = atan((br2TrapW2/2.-br2TrapW1/2.)/br2TrapH);
508 G4double delta13 = AlThickness*(1/cos(theta)-tan(theta));
509 G4double delta24 = AlThickness*(1/cos(theta)+tan(theta));
510 G4double delta13P = (PVFThickness+AlThickness)*(1/cos(theta)-tan(theta));
511 G4double delta24P = (PVFThickness+AlThickness)*(1/cos(theta)+tan(theta));
512
513 //construct barrel tof scintillator
514 G4Trap* solidScinBr2=new G4Trap("solidScinBr2",scinTrapdz,0*deg,0*deg,
515 scinTrapdy1,scinTrapdx1,scinTrapdx2,0*deg,
516 scinTrapdy2,scinTrapdx3,scinTrapdx4,0*deg);
517 logicalScinBr2 = new G4LogicalVolume(solidScinBr2,BC408,"logicalScinBr2");
518
519 //construct barrel Al film and PVF film
520 G4Trap* solidAlBr2 = new G4Trap("solidAlBr2", scinTrapdz+0.001,0*deg,0*deg,
521 scinTrapdy1+AlThickness, scinTrapdx1+delta13, scinTrapdx2+delta24 ,0*deg,
522 scinTrapdy2+AlThickness, scinTrapdx3+delta13, scinTrapdx4+delta24 ,0*deg);
523
524 logicalAlBr2 = new G4LogicalVolume(solidAlBr2,G4Material::GetMaterial("Aluminium"),"logicalAlBr2");
525 G4Trap* solidPVFBr2 = new G4Trap("solidPVFBr2",scinTrapdz+0.002,0*deg,0*deg,
526 scinTrapdy1+AlThickness+PVFThickness,scinTrapdx1+delta13P, scinTrapdx2+delta24P,0*deg,
527 scinTrapdy2+AlThickness+PVFThickness,scinTrapdx3+delta13P, scinTrapdx4+delta24P,0*deg);
528
529 logicalPVFBr2 = new G4LogicalVolume(solidPVFBr2,PVF,"logicalPVFBr2");
530 //put daughter in mother logical volume
531 physicalAlBr2 = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalAlBr2, "physicalAlBr2",logicalPVFBr2, false,0);
532 physicalScinBr2 = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalScinBr2, "physicalScinBr2",logicalAlBr2, false, 0);
533
534 //construct barrel PMT bucket
535 G4Tubs* solidBucketBr2= new G4Tubs("solidBucketBr2",0,bucketDBr/2.,
536 bucketLBr/2., 0, 360.*deg);
537 logicalBucketBr2=new G4LogicalVolume(solidBucketBr2,PMTmaterial,"logicalBucketBr2");
538
539 G4RotationMatrix trapRot;
540 trapRot.rotateZ(-90.*deg);
541 // /| y
542 //------- / | |
543 //\ / | | rotate->(-90.*deg): |
544 // \ / ->>>> | | -----x
545 // --- \ |
546 // \|
547 G4ThreeVector pos(0,0,0);
548 G4double x,y;
549 for (int i=0;i<nScinBr;i++)
550 {
551 ostringstream osnamePVFBr2;
552 osnamePVFBr2 << "physicalPVFBr2_"<<i;
553 ostringstream osnameBucket1Br2;
554 osnameBucket1Br2 << "physicalBucket1Br2_"<<i;
555 ostringstream osnameBucket2Br2;
556 osnameBucket2Br2 << "physicalBucket2Br2_"<<i;
557
558 //scintillator and casing
559 x=scinPosR*cos(i*angleBr); //+0.5: start from phi=0
560 y=scinPosR*sin(i*angleBr);
561 pos.setX(x);
562 pos.setY(y);
563 pos.setZ(0);
564 //to make code consistent with gdml
565 double a, b, c;
566 getXYZ(trapRot, a, b ,c);
567 G4RotationMatrix trapRotTemp;
568 trapRotTemp.rotateZ(c);
569 G4Transform3D transform(trapRotTemp,pos);
570 //G4Transform3D transform(trapRot,pos);
571 physicalPVFBr2 = new G4PVPlacement(transform,logicalPVFBr2,osnamePVFBr2.str(), logicalBrTof,false,i);
572 //bucket
573 pos.setZ(-(bucketLBr+br2L)/2.-0.01);
574 physicalBucket1Br2 = new G4PVPlacement(0,pos,logicalBucketBr2,osnameBucket1Br2.str(),logicalBrTof,false,0);
575 pos.setZ((bucketLBr+br2L)/2.+0.01);
576 physicalBucket2Br2 = new G4PVPlacement(0,pos,logicalBucketBr2,osnameBucket2Br2.str(),logicalBrTof,false,0);
577
578 //rotate for next scintillator
579 trapRot.rotateZ(angleBr);
580 }
581
582}
583
584void BesTofConstruction::ConstructEcTof()
585{
586 BesTofGeoParameter* tofPara = BesTofGeoParameter::GetInstance();
587
588 //Get all interesting parameters for the TOF construction
589 //endcap geometry data
590 G4int nScinEc = tofPara->GetnScinEc(); //number of endcap scintillators
591 G4double ecL = tofPara->GetEcL();
592 G4double ecTrapW1 = tofPara->GetEcTrapW1();
593 G4double ecTrapW2 = tofPara->GetEcTrapW2();
594 G4double ecTrapH = tofPara->GetEcTrapH();
595 G4double ecTrapH1 = tofPara->GetEcTrapH1();
596 G4double zPosEastEc = tofPara->GetzPosEastEc(); //z position of east endcap
597 G4double zPosWestEc = tofPara->GetzPosWestEc(); //z position of west endcap
598 G4double ecR1 = tofPara->GetEcR1();
599 //G4double ecR2 = tofPara.GetEcR2();
600
601 //endcap PMT bucket geometry data
602 G4double bucketDEc = tofPara->GetBucketDEc(); //diameter of endcap PMT bucket
603 G4double bucketLEc = tofPara->GetBucketLEc(); //length of endcap PMT bucket
604 G4double bucketPosR = tofPara->GetBucketPosR(); //R of bucket center
605
606 G4double AlThickness=tofPara->GetAlThickness();
607 G4double PVFThickness=tofPara->GetPVFThickness();
608
609 G4double angleEc=360./nScinEc*deg;
610 G4double ecTrapW3 = ecTrapW1 + (ecTrapW2-ecTrapW1)*ecTrapH1/ecTrapH;
611 G4double ecTrapH2 = ecTrapH-ecTrapH1;
612 G4double pdz = ecL/2;
613 G4double ptheta = atan(ecTrapH1/(2*ecL));
614 G4double pdy1 = ecTrapH2/2;
615 G4double pdx1 = ecTrapW3/2;
616 G4double pdx2 = ecTrapW2/2;
617 G4double pdy2 = ecTrapH/2;
618 G4double pdx3 = ecTrapW1/2;
619 G4double pdx4 = ecTrapW2/2;
620
621 //because of removing a heighth of H1,
622 //plus initial center position ecR=ecR1+ecTrapH/2 with ecTrapH1/4
623 G4double ecR = ecR1+ecTrapH/2+ecTrapH1/4;
624
625 //construct endcap scintillator
626 G4Trap* solidScinEc=new G4Trap("solidScinEc",pdz,ptheta, 270*deg,
627 pdy1,pdx1,pdx2,0*deg,pdy2,pdx3,pdx4,0*deg);
628
629 logicalScinEcWest=new G4LogicalVolume(solidScinEc,BC404,"logicalScinEcWest");
630 logicalScinEcEast=new G4LogicalVolume(solidScinEc,BC404,"logicalScinEcEast");
631
632 //construct endcap PMT bucket
633 G4Tubs* solidBucketEc = new G4Tubs("solidBucketEc", 0, bucketDEc/2.,
634 bucketLEc/2., 0, 360.*deg);
635 logicalBucketEc = new G4LogicalVolume(solidBucketEc, PMTmaterial, "logicalBucketEc");
636
637 //construct Al and PVF film
638 G4double pthetaAl = atan(ecTrapH1/(2*(ecL+AlThickness*2)));
639 G4double theta1 = atan((ecTrapW2/2.-ecTrapW3/2.)/ecTrapH2);
640 G4double theta2 = atan((ecTrapW2/2.-ecTrapW1/2.)/ecTrapH);
641
642 G4double delta1 = AlThickness*(1/cos(theta1)-tan(theta1));
643 G4double delta2 = AlThickness*(1/cos(theta1)+tan(theta1));
644 G4double delta3 = AlThickness*(1/cos(theta2)-tan(theta2));
645 G4double delta4 = AlThickness*(1/cos(theta2)+tan(theta2));
646
647 G4double thick = AlThickness+PVFThickness;
648 G4double pthetaPVF = atan(ecTrapH1/(2*(ecL+thick*2)));
649 G4double delta1P = thick*(1/cos(theta1)-tan(theta1));
650 G4double delta2P = thick*(1/cos(theta1)+tan(theta1));
651 G4double delta3P = thick*(1/cos(theta2)-tan(theta2));
652 G4double delta4P = thick*(1/cos(theta2)+tan(theta2));
653
654 G4Trap* solidAlEc= new G4Trap("solidAlEc",pdz+AlThickness,pthetaAl,270*deg,pdy1+AlThickness,pdx1+delta1,pdx2+delta2,0*deg,pdy2+AlThickness,pdx3+delta3,pdx4+delta4,0*deg);
655
656 logicalAlEcWest = new G4LogicalVolume(solidAlEc,G4Material::GetMaterial("Aluminium"),"logicalAlEcWest");
657 logicalAlEcEast = new G4LogicalVolume(solidAlEc,G4Material::GetMaterial("Aluminium"),"logicalAlEcEast");
658
659 G4Trap* solidPVFEc = new G4Trap("solidPVFEc",pdz+thick,pthetaPVF,270*deg,
660 pdy1+thick,pdx1+delta1P,pdx2+delta2P,0*deg,
661 pdy2+thick,pdx3+delta3P,pdx4+delta4P,0*deg);
662
663 logicalPVFEcWest = new G4LogicalVolume(solidPVFEc,PVF,"logicalPVFEcWest");
664 logicalPVFEcEast = new G4LogicalVolume(solidPVFEc,PVF,"logicalPVFEcEast");
665
666 //put scintilator in Al, then put Al in PVF
667 physicalAlEcWest = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalAlEcWest, "physicalAlEcWest",logicalPVFEcWest, false,0);
668 physicalScinEcWest = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalScinEcWest, "physicalScinEcWest",logicalAlEcWest, false, 0);
669
670
671 //put scintilator in Al, then put Al in PVF
672 physicalAlEcEast = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalAlEcEast, "physicalAlEcEast",logicalPVFEcEast, false,0);
673 physicalScinEcEast = new G4PVPlacement(0,G4ThreeVector(0,0,0),logicalScinEcEast, "physicalScinEcEast",logicalAlEcEast, false, 0);
674
675 //zPosWestEc<0, zPosEastEc>0, direction of axis z: west--->east
676 //construct west endcap first
677 //when no rotation, the first tub is start from phi=90
678 //after rotation, the first tub is start from 0, rotate angle = 0.5*angleEc
679 G4ThreeVector posW(0,0,0);
680 G4ThreeVector posE(0,0,0);
681 G4RotationMatrix trapRot;
682 trapRot.rotateY(180.*deg);
683 trapRot.rotateZ(0.5*angleEc-90*deg);
684 for (int i=0;i<nScinEc;i++)
685 {
686 ostringstream osnamePVFEc;
687 osnamePVFEc << "physicalPVFEcWest_"<<i;
688 ostringstream osnameBucketEc;
689 osnameBucketEc << "physicalBucketEcWest_"<<i;
690
691 posW.setX( ecR*cos((0.5+i)*angleEc) );
692 posW.setY( ecR*sin((0.5+i)*angleEc) );
693 //posW.setZ( zPosWestEc );
694 posW.setZ(0);
695
696 //to make code consistent with gdml
697 double a, b, c;
698 getXYZ(trapRot, a, b ,c);
699 G4RotationMatrix trapRotTemp;
700 trapRotTemp.rotateX(a);
701 trapRotTemp.rotateY(b);
702 trapRotTemp.rotateZ(c);
703 G4Transform3D transform(trapRotTemp,posW);
704 //G4Transform3D transform(trapRot,posW);
705 physicalPVFEcWest = new G4PVPlacement(transform,logicalPVFEcWest,osnamePVFEc.str(), logicalEcTofWest,false,i);
706 posW.setX( bucketPosR*cos((0.5+i)*angleEc) );
707 posW.setY( bucketPosR*sin((0.5+i)*angleEc) );
708 //posW.setZ( zPosWestEc-ecL/2-thick-bucketLEc/2-0.01 );
709 posW.setZ(-ecL/2-thick-bucketLEc/2-0.01 );
710 physicalBucketEc = new G4PVPlacement(0,posW,logicalBucketEc,osnameBucketEc.str(),logicalEcTofWest,false,0);
711 trapRot.rotateZ(angleEc); //pay attention: not i*angleEc
712 }
713 //to make east scintillator start from phi=0
714 // _ _
715 // _ - _ -
716 // - - - - ====> - - - - ====> - - - -
717 // ~ - _
718 //
719 trapRot.rotateZ(-angleEc);
720 trapRot.rotateX(180.*deg); //make east bucket point to the center
721 for (int i=0;i<nScinEc;i++)
722 {
723 ostringstream osnamePVFEc;
724 osnamePVFEc << "physicalPVFEcEast_"<<i;
725 ostringstream osnameBucketEc;
726 osnameBucketEc << "physicalBucketEcEast_"<<i;
727
728 posE.setX( ecR*cos((0.5+i)*angleEc) );
729 posE.setY( ecR*sin((0.5+i)*angleEc) );
730 //posE.setZ( zPosEastEc );
731 posE.setZ(0);
732
733 //to make code consistent with gdml
734 double a, b, c;
735 getXYZ(trapRot, a, b ,c);
736 G4RotationMatrix trapRotTemp;
737 trapRotTemp.rotateX(a);
738 trapRotTemp.rotateY(b);
739 trapRotTemp.rotateZ(c);
740 G4Transform3D transform(trapRotTemp,posE);
741 //G4Transform3D transform(trapRot,posE);
742 physicalPVFEcEast = new G4PVPlacement(transform,logicalPVFEcEast,osnamePVFEc.str(), logicalEcTofEast,false,i);
743 posE.setX( bucketPosR*cos((0.5+i)*angleEc) );
744 posE.setY( bucketPosR*sin((0.5+i)*angleEc) );
745 //posE.setZ( zPosEastEc+ecL/2+thick+bucketLEc/2+0.01 );
746 posE.setZ(ecL/2+thick+bucketLEc/2+0.01 );
747 physicalBucketEc = new G4PVPlacement(0,posE,logicalBucketEc,osnameBucketEc.str(),logicalEcTofEast,false,0);
748 trapRot.rotateZ(angleEc); //pay attention: not i*angleEc
749 }
750}
751
752
753//This function construct the TOF_Endcap_MRPC
754void BesTofConstruction::ConstructEcTof_mrpc()
755{
756 //I construct my detector in the Geant4 coordinate system. It will be rotated into the BES coordinatesystem later!
757 // |z ___________
758 // | /y \ /
759 // | / \ /
760 // |/______x \_____/
761
762
763 G4double smallL = 0.01*mm;
764 G4RotationMatrix rot_dummy(0*deg,0*deg,0*deg);
765 IniParam_mrpc();
766 PartProduce* partProduce = new PartProduce();
767 partProduce->IniParam();
768 bool checkOverlap = false;
769
770 //Get the parts and place them in the container
771 //container
772 std::ostringstream name;
773 //Create 4 containers and gasContainers
774 //0: top module, closer to IP, longer bottom cover;
775 //1: same to [0], with gas sensitive detector overturned;
776 //2: bottom module, next to EMC, shorter bottom cover;
777 //3: same to [2], with gas sensitive detector overturned
778 G4LogicalVolume* logical_container[4];
779 G4LogicalVolume* logical_gasContainer[4];
780 for(int kk=0; kk<4; kk++)
781 {
782 name.str("");
783 name<<"logical_container_m"<<kk;
784 logical_container[kk] = partProduce->lg_container((int)kk/2, name.str());
785
786 name.str("");
787 name<<"logical_gasContainer_m"<<kk;
788 logical_gasContainer[kk] = partProduce->lg_gasContainer(name.str());
789 }
790
791
792 //The parts that compose the container
793 G4LogicalVolume* logical_containerFrame = partProduce->lg_containerFrame("logical_containerFrame");
794 G4LogicalVolume* logical_upCover = partProduce->lg_upCover("logical_upCover");
795 G4LogicalVolume* logical_lowCover1 = partProduce->lg_lowCover1("logical_lowCover1");
796 G4LogicalVolume* logical_lowCover2 = partProduce->lg_lowCover2("logical_lowCover2");
797 G4LogicalVolume* logical_upFEE = partProduce->lg_upFEE("logical_upFEE");
798 G4LogicalVolume* logical_sideFEE = partProduce->lg_sideFEE("logical_sideFEE");
799 G4LogicalVolume* logical_castingDie = partProduce->lg_castingDie("logical_castingDie");
800 G4LogicalVolume* logical_bareChamber = partProduce->lg_bareChamber("logical_bareChamber");
801 G4LogicalVolume* logical_bracket = partProduce->lg_bracket("logical_bracket");
802 G4LogicalVolume* logical_sideStopBlock = partProduce->lg_sideStopBlock("logical_sideStopBlock");
803 G4LogicalVolume* logical_upStopBlock = partProduce->lg_upStopBlock("logical_upStopBlock");
804
805
806 //Place the parts in the container
807 for(int kk=0; kk<4; kk++)
808 {
809 //container frame
810 name.str("");
811 name<<"physical_containerFrame_m"<<kk;
812 new G4PVPlacement(0, G4ThreeVector(0,0,0), logical_containerFrame, name.str(), logical_container[kk], false, 0, checkOverlap);
813
814 //top cover
815 G4ThreeVector translation_upCover(0, (containerOut_y+upCover_y)/2+smallL, (upCover_z-containerOut_z)/2-upCover_posz);
816 name.str("");
817 name<<"physical_upCover_m"<<kk;
818 new G4PVPlacement(0, translation_upCover, logical_upCover, name.str(), logical_container[kk], false, 0, checkOverlap);
819
820 //bottom cover
821 name.str("");
822 name<<"physical_lowCover_m"<<kk;
823 if(kk<2)
824 {
825 G4ThreeVector translation_lowCover(0, -1*(containerOut_y+upCover_y)/2, (lowCover1_z-containerOut_z)/2-lowCover1_posz);
826 new G4PVPlacement(0, translation_lowCover, logical_lowCover1, name.str(), logical_container[kk], false, 0, checkOverlap);
827 }
828 else
829 {
830 G4ThreeVector translation_lowCover(0, -1*(containerOut_y+upCover_y)/2, (lowCover2_z-containerOut_z)/2-lowCover2_posz);
831 new G4PVPlacement(0, translation_lowCover, logical_lowCover2, name.str(), logical_container[kk], false, 0, checkOverlap);
832 }
833
834 //FEE on the top cover
835 G4ThreeVector translation_upFEE(0, containerOut_y/2+upCover_y+upFEE_y/2+smallL, upCover_z-upFEE_posz-upFEE_z/2-upCover_posz-containerOut_z/2);
836 name.str("");
837 name<<"physical_upFEE_m"<<kk;
838 new G4PVPlacement(0, translation_upFEE, logical_upFEE, name.str(), logical_container[kk], false, 0, checkOverlap);
839
840 //FEE at the side
841 double angle_containerFrameSide = atan((containerOut_xup-containerOut_xdown)/2/containerOut_z) *rad;
842 double l_containerFrameSide = containerOut_z/cos(angle_containerFrameSide);
843 double x_sideFEE = (l_containerFrameSide/2-sideFEE_posl)*sin(angle_containerFrameSide)+(containerOut_xup+containerOut_xdown)/4+sideFEE_x/2*cos(angle_containerFrameSide);
844 double z_sideFEE = (l_containerFrameSide/2-sideFEE_posl)*cos(angle_containerFrameSide)-sideFEE_x/2*sin(angle_containerFrameSide);
845 G4ThreeVector translation_sideFEE_left(-x_sideFEE, 0, z_sideFEE);
846 G4RotationMatrix rot_sideFEE_left;
847 rot_sideFEE_left.rotateY(-angle_containerFrameSide);
848 G4Transform3D transform_sideFEE_left(rot_sideFEE_left, translation_sideFEE_left);
849 name.str("");
850 name<<"physical_sideFEE_left_m"<<kk;
851 new G4PVPlacement(transform_sideFEE_left, logical_sideFEE, name.str(), logical_container[kk], false, 0, checkOverlap);
852
853 G4ThreeVector translation_sideFEE_right(x_sideFEE, 0, z_sideFEE);
854 G4RotationMatrix rot_sideFEE_right;
855 rot_sideFEE_right.rotateZ(180*deg);
856 rot_sideFEE_right.rotateY(angle_containerFrameSide);
857 G4Transform3D transform_sideFEE_right(rot_sideFEE_right, translation_sideFEE_right);
858 name.str("");
859 name<<"physical_sideFEE_right_m"<<kk;
860 new G4PVPlacement(transform_sideFEE_right, logical_sideFEE, name.str(), logical_container[kk], false, 0, checkOverlap);
861
862 //casting die
863 G4RotationMatrix rot_castingDie;
864 rot_castingDie.rotateX(-90*deg);
865 G4ThreeVector translation_castingDie_1(0, containerOut_y/2+upCover_y+castingDie_z/2+smallL, -containerOut_z/2+upCover_posz+castingDie_posz1);
866 G4Transform3D transform_castingDie_1(rot_castingDie, translation_castingDie_1);
867 name.str("");
868 name<<"physical_castingDie_1_m"<<kk;
869 new G4PVPlacement(transform_castingDie_1, logical_castingDie, name.str(), logical_container[kk], false, 0, checkOverlap);
870
871 G4ThreeVector translation_castingDie_2(0, containerOut_y/2+upCover_y+castingDie_z/2+smallL, -containerOut_z/2+upCover_posz+castingDie_posz2);
872 G4Transform3D transform_castingDie_2(rot_castingDie, translation_castingDie_2);
873 name.str("");
874 name<<"physical_castingDie_2_m"<<kk;
875 new G4PVPlacement(transform_castingDie_2, logical_castingDie, name.str(), logical_container[kk], false, 0, checkOverlap);
876
877 //gas container
878 G4ThreeVector translation_gasContainer(0, 0, (container_lowBorder-container_thickness)/2);
879 name.str("");
880 name<<"physical_gasContainer_m"<<kk;
881 new G4PVPlacement(0, translation_gasContainer, logical_gasContainer[kk], name.str(), logical_container[kk], false, 0, checkOverlap);
882 }
883
884
885 //Fill in the inside of the container
886 //1mm from the bottom cover in design. Considering the modules's position are not real, change this number to make the bare chamber in the right place
887 double l_chamber = honeyComb_y*2+tape_mrpc_y*2+pcbBoard_y0*2+pcbBoard_y1+mylar_y*4+glass0_y*4+glass1_y*10+gasLayer_y*12;
888 double y_chamber;
889 double z_chamber = (containerIn_z-pcbBoard_z)/2-smallL;
890 for(int kk=0; kk<4; kk++)
891 {
892 if(kk<2) y_chamber = -(containerIn_y-l_chamber)/2+chamber_in_posy;
893 else y_chamber = -(containerIn_y-l_chamber)/2+chamber_out_posy;
894 G4Transform3D transform_bareChamber;
895 if(kk==0 || kk==2)
896 {
897 transform_bareChamber = G4Transform3D(rot_dummy, G4ThreeVector(0,y_chamber,z_chamber));
898 }
899 else
900 {
901 G4RotationMatrix rot_bareChamber;
902 rot_bareChamber.rotateZ(180*deg);
903 transform_bareChamber = G4Transform3D(rot_bareChamber, G4ThreeVector(0,y_chamber,z_chamber));
904 }
905 new G4PVPlacement(transform_bareChamber, logical_bareChamber, "physical_bareChamber", logical_gasContainer[kk], false, 0, checkOverlap);
906
907 //lower bracket
908 G4double z_bracket = -(containerIn_z-bracket_z)/2+smallL;
909 new G4PVPlacement(0, G4ThreeVector(-bracket_posx, 0, z_bracket), logical_bracket, "physical_bracket_0", logical_gasContainer[kk], false, 0, checkOverlap);
910 new G4PVPlacement(0, G4ThreeVector(bracket_posx, 0, z_bracket), logical_bracket, "physical_bracket_1", logical_gasContainer[kk], false, 0, checkOverlap);
911
912 //side stop block
913 G4Transform3D transform_sideStopBlock;
914 G4double angle_gasContainerSide = (atan((containerIn_xup-containerIn_xdown)/2/(containerIn_z))) *rad;
915 G4double x_sideStopBlock = (containerIn_xup+containerIn_xdown)/4+sideStopBlock_posz0*tan(angle_gasContainerSide)-sideStopBlock_x/2*cos(angle_gasContainerSide);
916 G4double z_sideStopBlock = sideStopBlock_posz0+sideStopBlock_x/2*sin(angle_gasContainerSide);
917
918 G4RotationMatrix rot_sideStopBlock;
919 rot_sideStopBlock.rotateY(angle_gasContainerSide);
920 rot_sideStopBlock.rotateZ(180*deg);
921 transform_sideStopBlock = G4Transform3D(rot_sideStopBlock, G4ThreeVector(-x_sideStopBlock, 0, z_sideStopBlock));
922 new G4PVPlacement(transform_sideStopBlock, logical_sideStopBlock, "physical_sideStopBlock_0", logical_gasContainer[kk], false, 0, checkOverlap);
923
924 rot_sideStopBlock = G4RotationMatrix(0,0,0);
925 rot_sideStopBlock.rotateY(angle_gasContainerSide);
926 transform_sideStopBlock = G4Transform3D(rot_sideStopBlock, G4ThreeVector(x_sideStopBlock, 0, z_sideStopBlock));
927 new G4PVPlacement(transform_sideStopBlock, logical_sideStopBlock, "physical_sideStopBlock_1", logical_gasContainer[kk], false, 0, checkOverlap);
928
929 x_sideStopBlock = (containerIn_xup+containerIn_xdown)/4+sideStopBlock_posz1*tan(angle_gasContainerSide)-sideStopBlock_x/2*cos(angle_gasContainerSide);
930 z_sideStopBlock = sideStopBlock_posz1+sideStopBlock_x/2*sin(angle_gasContainerSide);
931 rot_sideStopBlock = G4RotationMatrix(0,0,0);
932 rot_sideStopBlock.rotateY(angle_gasContainerSide);
933 rot_sideStopBlock.rotateZ(180*deg);
934 transform_sideStopBlock = G4Transform3D(rot_sideStopBlock, G4ThreeVector(-x_sideStopBlock, 0, z_sideStopBlock));
935 new G4PVPlacement(transform_sideStopBlock, logical_sideStopBlock, "physical_sideStopBlock_2", logical_gasContainer[kk], false, 0, checkOverlap);
936
937 rot_sideStopBlock = G4RotationMatrix(0,0,0);
938 rot_sideStopBlock.rotateY(angle_gasContainerSide);
939 transform_sideStopBlock = G4Transform3D(rot_sideStopBlock, G4ThreeVector(x_sideStopBlock, 0, z_sideStopBlock));
940 new G4PVPlacement(transform_sideStopBlock, logical_sideStopBlock, "physical_sideStopBlock_3", logical_gasContainer[kk], false, 0, checkOverlap);
941 }
942
943 //stop block
944 //upper stop block
945 G4double x_upStopBlock = containerIn_xup/4;
946 G4double y_upStopBlock = pcbBoard_y1/2.+mylar_y+glass0_y+glass1_y*2.5+gasLayer_y*3;
947 G4double z_upStopBlock = (pcbBoard_z-upStopBlock_z)/2-smallL;
948 new G4PVPlacement(0, G4ThreeVector(-x_upStopBlock, -y_upStopBlock, z_upStopBlock ), logical_upStopBlock, "physical_upStopBlock_0", logical_bareChamber, false,0, checkOverlap);
949 new G4PVPlacement(0, G4ThreeVector(-x_upStopBlock, y_upStopBlock, z_upStopBlock ), logical_upStopBlock, "physical_upStopBlock_1", logical_bareChamber, false,0, checkOverlap);
950 new G4PVPlacement(0, G4ThreeVector(x_upStopBlock, -y_upStopBlock, z_upStopBlock ), logical_upStopBlock, "physical_upStopBlock_2", logical_bareChamber, false,0, checkOverlap);
951 new G4PVPlacement(0, G4ThreeVector(x_upStopBlock, y_upStopBlock, z_upStopBlock ), logical_upStopBlock, "physical_upStopBlock_3", logical_bareChamber, false,0, checkOverlap);
952
953
954 //honeyComb
955 G4LogicalVolume* logical_honeyComb = partProduce->lg_honeyComb("logical_honeyComb");
956 G4double y_honeyComb = (honeyComb_y+pcbBoard_y1)/2+tape_mrpc_y+pcbBoard_y0+mylar_y*2+glass0_y*2+glass1_y*5+gasLayer_y*6;
957 G4double z_honeyComb = (pcbBoard_z-honeyComb_z)/2-honeyComb_posz;
958 new G4PVPlacement(0,G4ThreeVector(0,y_honeyComb,z_honeyComb), logical_honeyComb,"physical_honeyComb_0",logical_bareChamber,false,0, checkOverlap);
959 new G4PVPlacement(0,G4ThreeVector(0,-1*y_honeyComb,z_honeyComb), logical_honeyComb,"physical_honeyComb_1",logical_bareChamber,false,1, checkOverlap);
960
961 //tape_mrpc
962 G4LogicalVolume* logical_tape_mrpc = partProduce->lg_tape_mrpc("logical_tape_mrpc");
963 G4double y_tape_mrpc = (tape_mrpc_y+pcbBoard_y1)/2+pcbBoard_y0+mylar_y*2+glass0_y*2+glass1_y*5+gasLayer_y*6;
964 G4double z_tape_mrpc = (pcbBoard_z-tape_mrpc_z)/2-tape_mrpc_posz;
965 new G4PVPlacement(0,G4ThreeVector(0,y_tape_mrpc,z_tape_mrpc), logical_tape_mrpc,"physical_tape_mrpc_0",logical_bareChamber,false,0, checkOverlap);
966 new G4PVPlacement(0,G4ThreeVector(0,-1*y_tape_mrpc,z_tape_mrpc), logical_tape_mrpc,"physical_tape_mrpc_1",logical_bareChamber,false,1, checkOverlap);
967
968 //pcbBoard
969 G4LogicalVolume* logical_pcbBoard0 = partProduce->lg_pcbBoard(0, "logical_pcbBoard0");
970 G4double y_pcbBoard = (pcbBoard_y0+pcbBoard_y1)/2+mylar_y*2+glass0_y*2+glass1_y*5+gasLayer_y*6;
971 new G4PVPlacement(0,G4ThreeVector(0,y_pcbBoard,0),logical_pcbBoard0,"physical_pcbBoard0_0",logical_bareChamber,false,0, checkOverlap);
972 //Notice!!! I rotate 180 deg in order to make the strips near the PCB surface,
973 //but keep in mind that this will make the readout strips not mirror symmetry, but Z axial symmetry
974 G4RotationMatrix* rot_pcbBoard = new G4RotationMatrix();
975 rot_pcbBoard->rotateZ(180*degree);
976 new G4PVPlacement(rot_pcbBoard,G4ThreeVector(0,-1*y_pcbBoard,0),logical_pcbBoard0,"physical_pcbBoard0_1",logical_bareChamber,false,1, checkOverlap);
977
978 G4LogicalVolume* logical_pcbBoard1 = partProduce->lg_pcbBoard(1, "logical_pcbBoard1");
979 new G4PVPlacement(0,G4ThreeVector(0,0,0),logical_pcbBoard1,"physical_pcbBoard1",logical_bareChamber,false,0, checkOverlap);
980
981 //readout strip
982 G4AssemblyVolume* logical_strip = partProduce->lg_strip("logical_strip");
983 G4double z_strip = (pcbBoard_z-12*strip_z-11*strip_gap)/2-strip_posz;
984 G4Transform3D transform_strip(rot_dummy, G4ThreeVector(0, pcbBoard_y0/2.-strip_y/2-smallL, z_strip));
985 logical_strip->MakeImprint(logical_pcbBoard0, transform_strip);
986 transform_strip = G4Transform3D(rot_dummy, G4ThreeVector(0,0,z_strip));
987 logical_strip->MakeImprint(logical_pcbBoard1, transform_strip);
988
989 //mylar
990 G4LogicalVolume* logical_mylar = partProduce->lg_mylar("logical_mylar");
991 G4double y_mylar = (mylar_y+pcbBoard_y1)/2+mylar_y+glass0_y*2+glass1_y*5+gasLayer_y*6;
992 G4double z_mylar = (pcbBoard_z-mylar_z)/2-mylar_posz;
993 new G4PVPlacement(0,G4ThreeVector(0,y_mylar,z_mylar),logical_mylar,"physical_mylar_0",logical_bareChamber,false,0, checkOverlap);
994 new G4PVPlacement(0,G4ThreeVector(0,-y_mylar,z_mylar),logical_mylar,"physical_mylar_3",logical_bareChamber,false,3, checkOverlap);
995
996 y_mylar = (mylar_y+pcbBoard_y1)/2;
997 new G4PVPlacement(0,G4ThreeVector(0,y_mylar,z_mylar),logical_mylar,"physical_mylar_1",logical_bareChamber,false,1, checkOverlap);
998 new G4PVPlacement(0,G4ThreeVector(0,-y_mylar,z_mylar),logical_mylar,"physical_mylar_2",logical_bareChamber,false,2, checkOverlap);
999
1000 //glass
1001 G4LogicalVolume* logical_glass0 = partProduce->lg_glass(0, "logical_glass0");
1002 G4double y_glass = (glass0_y+pcbBoard_y1)/2.+mylar_y+glass0_y+glass1_y*5+gasLayer_y*6;
1003 G4double z_glass = (pcbBoard_z-glass0_z)/2.-glass0_posz;
1004 new G4PVPlacement(0, G4ThreeVector(0, y_glass, z_glass), logical_glass0, "physical_glass0_0", logical_bareChamber, false, 0, checkOverlap);
1005 new G4PVPlacement(0, G4ThreeVector(0, -y_glass, z_glass), logical_glass0, "physical_glass0_3", logical_bareChamber, false, 3, checkOverlap);
1006 y_glass = pcbBoard_y1/2.+mylar_y+glass0_y/2.;
1007 new G4PVPlacement(0, G4ThreeVector(0, y_glass, z_glass), logical_glass0, "physical_glass0_1", logical_bareChamber, false, 1, checkOverlap);
1008 new G4PVPlacement(0, G4ThreeVector(0, -y_glass, z_glass), logical_glass0, "physical_glass0_2", logical_bareChamber, false, 2, checkOverlap);
1009
1010 G4LogicalVolume* logical_glass1 = partProduce->lg_glass(1, "logical_glass1");
1011 z_glass = (pcbBoard_z-glass1_z)/2.-glass1_posz;
1012 for(G4int i=0; i<5; i++)
1013 {
1014 y_glass = pcbBoard_y1/2.+mylar_y+glass0_y+glass1_y*(4.5-i)+gasLayer_y*(5-i);
1015 name.str("");
1016 name<<"physical_glass1_"<<i;
1017 new G4PVPlacement(0, G4ThreeVector(0, y_glass, z_glass), logical_glass1, name.str(), logical_bareChamber, false, i, checkOverlap);
1018 name.str("");
1019 name<<"physical_glass1_"<<9-i;
1020 new G4PVPlacement(0, G4ThreeVector(0, -y_glass, z_glass), logical_glass1, name.str(), logical_bareChamber, false, 9-i, checkOverlap);
1021 }
1022
1023 //gas
1024 logical_gasLayer = partProduce->lg_gasLayer("logical_gasLayer");
1025 G4double y_gasLayer;
1026 G4double z_gasLayer = (pcbBoard_z-gasLayer_z)/2.-gasLayer_posz;
1027 G4VPhysicalVolume* physical_gasLayer[12];
1028 for(G4int i=0; i<6; i++) //y->larger, gasNp->larger
1029 {
1030 y_gasLayer = pcbBoard_y1/2.+mylar_y+glass0_y+glass1_y*(5-i)+gasLayer_y*(5.5-i);
1031 name.str("");
1032 name<<"physical_gasLayer_"<<11-i;
1033 physical_gasLayer[i] = new G4PVPlacement(0, G4ThreeVector(0, -y_gasLayer, z_gasLayer), logical_gasLayer, name.str(), logical_bareChamber, false, 11-i, checkOverlap);
1034 }
1035 for(G4int i=6; i<12; i++)
1036 {
1037 y_gasLayer = pcbBoard_y1/2.+mylar_y+glass0_y+glass1_y*(i-6)+gasLayer_y*(-5.5+i);
1038 name.str("");
1039 name<<"physical_gasLayer_"<<11-i;
1040 physical_gasLayer[i] = new G4PVPlacement(0, G4ThreeVector(0, y_gasLayer, z_gasLayer), logical_gasLayer, name.str(), logical_bareChamber, false, 11-i, checkOverlap);
1041 }
1042
1043
1044 //arrange the 72 modules in the endcap, it's assumed to form two semicircles with a gap of about 2mm
1045
1046 // Z: East; y: Up; x: North
1047 // /y |z y y|/|
1048 // |z/ ___________ rotateX | / rotateZ /| |
1049 // |/_____x \ / ___|_/_ ||_|____x
1050 // \ / ->>>> / |/__\___x ->>>> | |
1051 // small end \_____/ bigger end /_________\ \ |
1052 // \|
1053
1054 const G4int n_module = 36;
1055 G4double angle = 360.0*deg/n_module;
1056 G4double z_layerIn = endcap_length/2-containerOut_y/2-lowCover1_y-layer_posz;
1057 G4double z_layerOut = endcap_length/2-containerOut_y/2-lowCover2_y-layer_posz;
1058 //Adjusting parameters
1059 //rOffset_east: along R direction
1060 //angle_east: along phi direction
1061 //angleOffset_east: along self asymmetric axis
1062 rOffset_east = tofPara->GetVec("rOffset_east");
1063 angle_east = tofPara->GetVec("angle_east");
1064 angleOffset_east = tofPara->GetVec("angleOffset_east");
1065 rOffset_west = tofPara->GetVec("rOffset_west");
1066 angle_west = tofPara->GetVec("angle_west");
1067 angleOffset_west = tofPara->GetVec("angleOffset_west");
1068
1069
1070 //arrange east endcap
1071 for(int i=0; i<n_module; i++)
1072 {
1073 G4double angle_module = startAngle_east+angle_east[i]+i*angle;
1074 G4double r_module = endcap_r+rOffset_east[i];
1075
1076 G4RotationMatrix rot_layerIn_east;
1077 rot_layerIn_east.rotateX(90.*deg);
1078 rot_layerIn_east.rotateZ(90.*deg);
1079 rot_layerIn_east.rotateZ(angle_module+angleOffset_east[i]);
1080
1081 G4ThreeVector translation_layerIn_east = G4ThreeVector(r_module*cos(angle_module), r_module*sin(angle_module), -z_layerIn);
1082 G4Transform3D transform_layerIn_east = G4Transform3D(rot_layerIn_east, translation_layerIn_east);
1083
1084
1085 G4RotationMatrix rot_layerOut_east;
1086 rot_layerOut_east.rotateZ(180.*deg);
1087 rot_layerOut_east.rotateX(90.*deg);
1088 rot_layerOut_east.rotateZ(90.*deg);
1089 rot_layerOut_east.rotateZ(angle_module+angleOffset_east[i]);
1090
1091 G4ThreeVector translation_layerOut_east = G4ThreeVector(r_module*cos(angle_module), r_module*sin(angle_module), z_layerOut);
1092 G4Transform3D transform_layerOut_east = G4Transform3D(rot_layerOut_east, translation_layerOut_east);
1093
1094 name.str("");
1095 name<<"physical_mrpc_east_"<<i;
1096 if(i%2==0)
1097 {
1098 new G4PVPlacement(transform_layerOut_east, logical_container[3], name.str(), logicalEcTofEast, false, i, checkOverlap);
1099 }
1100 else
1101 {
1102 new G4PVPlacement(transform_layerIn_east, logical_container[0], name.str(), logicalEcTofEast, false, i, checkOverlap);
1103 }
1104 }
1105
1106 //arrange west endcap
1107 for(int i=0; i<n_module; i++)
1108 {
1109 G4double angle_module = startAngle_west+angle_west[i]+i*angle;
1110 G4double r_module = endcap_r+rOffset_west[i];
1111
1112 G4RotationMatrix rot_layerIn_west;
1113 rot_layerIn_west.rotateZ(180.*deg);
1114 rot_layerIn_west.rotateX(90.*deg);
1115 rot_layerIn_west.rotateZ(90.*deg);
1116 rot_layerIn_west.rotateZ(angle_module+angleOffset_west[i]);
1117
1118 G4ThreeVector translation_layerIn_west = G4ThreeVector(r_module*cos(angle_module), r_module*sin(angle_module), z_layerIn);
1119 G4Transform3D transform_layerIn_west = G4Transform3D(rot_layerIn_west, translation_layerIn_west);
1120
1121
1122 G4RotationMatrix rot_layerOut_west;
1123 rot_layerOut_west.rotateX(90.*deg);
1124 rot_layerOut_west.rotateZ(90.*deg);
1125 rot_layerOut_west.rotateZ(angle_module+angleOffset_west[i]);
1126
1127 G4ThreeVector translation_layerOut_west = G4ThreeVector(r_module*cos(angle_module), r_module*sin(angle_module), -z_layerOut);
1128 G4Transform3D transform_layerOut_west = G4Transform3D(rot_layerOut_west, translation_layerOut_west);
1129
1130 name.str("");
1131 name<<"physical_mrpc_west_"<<i;
1132 if(i%2==0)
1133 {
1134 new G4PVPlacement(transform_layerOut_west, logical_container[2], name.str(), logicalEcTofWest, false, i, checkOverlap);
1135 }
1136 else
1137 {
1138 new G4PVPlacement(transform_layerIn_west, logical_container[1], name.str(), logicalEcTofWest, false, i, checkOverlap);
1139 }
1140 }
1141}
1142
1143
1144
1145void BesTofConstruction::DefineMaterial() //This fucntion just "build" my materials
1146{
1147 G4double a,z,density,fraction;
1148 G4int nel, natoms, ncomponents;
1149 G4String name, symbol;
1150
1151
1152 G4Element* H=G4Element::GetElement("Hydrogen");
1153 if (!H)
1154 H = new G4Element(name="Hydrogen",symbol="H" , z= 1., a=1.01*g/mole);
1155
1156 G4Element* C=G4Element::GetElement("Carbon");
1157 if (!C)
1158 C = new G4Element(name="Carbon" ,symbol="C" , z= 6., a=12.01*g/mole);
1159
1160 G4Element* F=G4Element::GetElement("Fluorin");
1161 if (!F)
1162 F = new G4Element(name="Fluorin", symbol="F", z=9., a=18.01*g/mole);
1163
1164 G4Element* O=G4Element::GetElement("Oxygen");
1165 if (!O)
1166 O = new G4Element(name="Oxygen",symbol="O",z= 8., a= 16.00*g/mole);
1167
1168 G4Element* N=G4Element::GetElement("Nitrogen");
1169 if (!N)
1170 N = new G4Element(name="Nitrogen",symbol="N",z= 7., a= 14.01*g/mole);
1171
1172 G4Element* S=G4Element::GetElement("Sulfur");
1173 if(!S)
1174 S = new G4Element(name="Sulfur",symbol="S",z= 16., a= 32.06*g/mole);
1175
1176
1177 BC404 = new G4Material("BC404",density=1.032*g/cm3,nel=2);
1178 BC404->AddElement(C,10);
1179 BC404->AddElement(H,11);
1180
1181 BC408 = new G4Material("BC408",density=1.032*g/cm3,nel=2);
1182 BC408->AddElement(C,1000);
1183 BC408->AddElement(H,1104);
1184
1185 PVF = new G4Material("PVF",density=1.45*g/cm3,nel=3);
1186 PVF->AddElement(C,2);
1187 PVF->AddElement(H,3);
1188 PVF->AddElement(F,1);
1189
1190 //PMT mixed material
1191 G4Material* Cu = G4Material::GetMaterial("Copper");
1192 G4Material* Al = G4Material::GetMaterial("Aluminium");
1193
1194 //vacuum
1195 //G4Material* Vacuum = new G4Material("Galactic", z=1., a=1.01*g/mole,
1196 //density= universe_mean_density, kStateGas, 3.e-18*pascal, 2.73*kelvin);
1197
1198 density = 1.4618815*g/cm3;
1199 PMTmaterial = new G4Material(name="PMTmaterial", density, 4);
1200 PMTmaterial->AddMaterial(Al, 0.4495);
1201 PMTmaterial->AddMaterial(Cu,0.35);
1202 PMTmaterial->AddMaterial(G4Material::GetMaterial("SiO2"),0.2);
1203 PMTmaterial->AddMaterial(G4Material::GetMaterial("Air"),0.0005);
1204 //PMTmaterial->AddMaterial(Vacuum,0.1);
1205 G4cout<<PMTmaterial;
1206
1207 density = 1.002*g/cm3;
1208 G4Material* tape = new G4Material(name="tape",density,nel=2);
1209 tape->AddElement(C,2);
1210 tape->AddElement(H,5);
1211
1212
1213 //MRPC material
1214 //The stop block in the gas container
1215 G4Material* Nylon = new G4Material("Nylon", density=1.15*g/cm3, ncomponents=4);
1216 Nylon->AddElement(C, natoms=2);
1217 Nylon->AddElement(H, natoms=3);
1218 Nylon->AddElement(O, natoms=1);
1219 Nylon->AddElement(N, natoms=1);
1220 G4cout<<Nylon<<G4endl;
1221
1222 //The core of honeyComb is thought to be cellulose
1223 G4Material* honeycombCore = new G4Material("honeycombCore", density=0.024*g/cm3, ncomponents=3);
1224 honeycombCore->AddElement(C,natoms=6);
1225 honeycombCore->AddElement(H,natoms=10);
1226 honeycombCore->AddElement(O,natoms=5);
1227 G4cout<<honeycombCore<<G4endl;
1228
1229 //The honeyComb surface and the PCB board are thought to be composed of SiO2(60%) and epoxy resins(40%)
1230 G4Material* epoxy = new G4Material("epoxy", density=1.2*g/cm3, ncomponents=3);
1231 epoxy->AddElement(C,natoms=11);
1232 epoxy->AddElement(H,natoms=12);
1233 epoxy->AddElement(O,natoms=3);
1234 G4cout<<epoxy<<G4endl;
1235
1236 G4Material* insulationBoard = new G4Material("insulationBoard", density=1.85*g/cm3, 2);
1237 insulationBoard->AddMaterial(G4Material::GetMaterial("SiO2"), 0.6);
1238 insulationBoard->AddMaterial(epoxy, 0.4);
1239 G4cout<<insulationBoard<<G4endl;
1240
1241 //Mylar and tape are thought to be PET
1242 G4Material* PET = new G4Material("PET", density= 1.39*g/cm3, ncomponents=3);
1243 PET->AddElement(C,natoms=10);
1244 PET->AddElement(H,natoms=8);
1245 PET->AddElement(O,natoms=4);
1246 G4cout<<PET<<G4endl;
1247
1248 //MRPC gas: 90% r134a, 5% isobutan, 5% SF6
1249 G4Material* FreonR134A = new G4Material("FreonR134A", density =4.241*mg/cm3,ncomponents=3);
1250 FreonR134A->AddElement(C,natoms=2);
1251 FreonR134A->AddElement(H,natoms=2);
1252 FreonR134A->AddElement(F,natoms=4);
1253 G4cout<<FreonR134A<<G4endl;
1254
1255 G4Material* SF6 = new G4Material("SF6", density =6.14*mg/cm3,ncomponents=2);
1256 SF6->AddElement(S,natoms=1);
1257 SF6->AddElement(F,natoms=6);
1258 G4cout<<SF6<<G4endl;
1259
1260 G4Material* Isobutan = new G4Material("Isobutan", density =2.487*mg/cm3,ncomponents=2);
1261 Isobutan->AddElement(C, natoms=4);
1262 Isobutan->AddElement(H, natoms=10);
1263 G4cout<<Isobutan<<G4endl;
1264
1265 G4Material* MRPCGas = new G4Material(name="MRPCGas", density=4.17*mg/cm3,ncomponents=3);
1266 MRPCGas->AddMaterial(FreonR134A, fraction= 89.69*perCent);
1267 MRPCGas->AddMaterial(SF6, fraction=7.34*perCent);
1268 MRPCGas->AddMaterial(Isobutan, fraction=2.97*perCent);
1269
1270 //This setting is to produce ion pairs, default is 0
1271 if(0.0 == MRPCGas->GetIonisation()->GetMeanEnergyPerIonPair())
1272 {
1273 //MRPCGas->GetIonisation()->SetMeanEnergyPerIonPair(20*eV);
1274 MRPCGas->GetIonisation()->SetMeanEnergyPerIonPair(m_ionE*eV);
1275 }
1276 G4cout<<MRPCGas<<G4endl;
1277}
1278
1279void BesTofConstruction::getXYZ( const G4RotationMatrix r, double& a, double& b, double& c) const
1280{
1281 double cosb = sqrt( r.xx()*r.xx() + r.yx()*r.yx() );
1282
1283 if (cosb > 16*FLT_EPSILON)
1284 {
1285 a = atan2( r.zy(), r.zz() );
1286 b = atan2( -r.zx(), cosb );
1287 c = atan2( r.yx(), r.xx() );
1288 }
1289 else
1290 {
1291 a = atan2( -r.yz(), r.yy() );
1292 b = atan2( -r.zx(), cosb );
1293 c = 0.;
1294 }
1295 //std::cout.precision(20);
1296 //std::cout<<"in getXYZ :( "<<r.xx()<<" "<<r.xy()<<" "<<r.xz()<<std::endl;
1297 //std::cout<<" "<<r.yx()<<" "<<r.yy()<<" "<<r.yz()<<std::endl;
1298 //std::cout<<" "<<r.zx()<<" "<<r.zy()<<" "<<r.zz()<<std::endl;
1299 //std::cout<<"details: cosb="<<cosb<<" a="<<a<<" b="<<b<<" c="<<c<<" min="<<16*FLT_EPSILON<<std::endl;
1300
1301
1302 G4RotationMatrix temp;
1303 temp.rotateZ((c/deg)*deg);
1304 //std::cout.precision(20);
1305 //std::cout<<"in detail temp2:( "<<temp.xx()<<" "<<temp.xy()<<" "<<temp.xz()<<std::endl;
1306 //std::cout<<" "<<temp.yx()<<" "<<temp.yy()<<" "<<temp.yz()<<std::endl;
1307 //std::cout<<" "<<temp.zx()<<" "<<temp.zy()<<" "<<temp.zz()<<std::endl;
1308}
1309
1310
1311void BesTofConstruction::IniParam_mrpc()
1312{
1313 tofPara = BesTofGeoParameter::GetInstance();
1314
1315 //All the values are in units of mm or deg
1316 containerOut_xdown = tofPara->Get("containerOut_xdown");
1317 containerOut_xup = tofPara->Get("containerOut_xup");
1318 containerOut_y = tofPara->Get("containerOut_y");
1319 containerOut_z = tofPara->Get("containerOut_z");
1320 container_thickness = tofPara->Get("container_thickness");
1321 container_lowBorder = tofPara->Get("container_lowBorder");
1322
1323 containerIn_xdown = tofPara->Get("containerIn_xdown");
1324 containerIn_xup = tofPara->Get("containerIn_xup");
1325 containerIn_y = tofPara->Get("containerIn_y");
1326 containerIn_z = tofPara->Get("containerIn_z");
1327
1328 upCover_xdown = tofPara->Get("upCover_xdown");
1329 upCover_xup = tofPara->Get("upCover_xup");
1330 upCover_y = tofPara->Get("upCover_y");
1331 upCover_z = tofPara->Get("upCover_z");
1332 upCover_posz = tofPara->Get("upCover_posz");
1333
1334 lowCover1_xdown = tofPara->Get("lowCover1_xdown");
1335 lowCover1_xup = tofPara->Get("lowCover1_xup");
1336 lowCover1_y = tofPara->Get("lowCover1_y");
1337 lowCover1_z = tofPara->Get("lowCover1_z");
1338 lowCover1_addx = tofPara->Get("lowCover1_addx");
1339 lowCover1_addz = tofPara->Get("lowCover1_addz");
1340 lowCover1_posz = tofPara->Get("lowCover1_posz");
1341
1342 lowCover2_xdown = tofPara->Get("lowCover2_xdown");
1343 lowCover2_xup = tofPara->Get("lowCover2_xup");
1344 lowCover2_y = tofPara->Get("lowCover2_y");
1345 lowCover2_z = tofPara->Get("lowCover2_z");
1346 lowCover2_posz = tofPara->Get("lowCover2_posz");
1347
1348 upFEE_xdown = tofPara->Get("upFEE_xdown");
1349 upFEE_xup = tofPara->Get("upFEE_xup");
1350 upFEE_y = tofPara->Get("upFEE_y");
1351 upFEE_z = tofPara->Get("upFEE_z");
1352 FEE_thickness = tofPara->Get("FEE_thickness");
1353 upFEE_posz = tofPara->Get("upFEE_posz");
1354
1355 upFEE_pcb_xdown = tofPara->Get("upFEE_pcb_xdown");
1356 upFEE_pcb_xup = tofPara->Get("upFEE_pcb_xup");
1357 upFEE_pcb_y = tofPara->Get("upFEE_pcb_y");
1358 upFEE_pcb_z = tofPara->Get("upFEE_pcb_z");
1359
1360 sideFEE_x = tofPara->Get("sideFEE_x");
1361 sideFEE_y = tofPara->Get("sideFEE_y");
1362 sideFEE_z = tofPara->Get("sideFEE_z");
1363 sideFEE_posl = tofPara->Get("sideFEE_posl");
1364
1365 sideFEE_pcb_x = tofPara->Get("sideFEE_pcb_x");
1366 sideFEE_pcb_y = tofPara->Get("sideFEE_pcb_y");
1367 sideFEE_pcb_z = tofPara->Get("sideFEE_pcb_z");
1368
1369 castingDie_x = tofPara->Get("castingDie_x");
1370 castingDie_y = tofPara->Get("castingDie_y");
1371 castingDie_z = tofPara->Get("castingDie_z");
1372 castingDie_addr = tofPara->Get("castingDie_addr");
1373 castingDie_subtractx = tofPara->Get("castingDie_subtractx");
1374 castingDie_subtract_posz = tofPara->Get("castingDie_subtract_posz");
1375 castingDie_posz1 = tofPara->Get("castingDie_posz1");
1376 castingDie_posz2 = tofPara->Get("castingDie_posz2");
1377
1378 chamber_in_posy = tofPara->Get("chamber_in_posy");
1379 chamber_out_posy = tofPara->Get("chamber_out_posy");
1380
1381 bracket_x = tofPara->Get("bracket_x");
1382 bracket_y = tofPara->Get("bracket_y");
1383 bracket_z = tofPara->Get("bracket_z");
1384 bracket_r = tofPara->Get("bracket_r");
1385 bracket_posx = tofPara->Get("bracket_posx");
1386
1387 upStopBlock_x = tofPara->Get("upStopBlock_x");
1388 upStopBlock_y = tofPara->Get("upStopBlock_y");
1389 upStopBlock_z = tofPara->Get("upStopBlock_z");
1390 upStopBlock_subtractr = tofPara->Get("upStopBlock_subtractr");
1391 upStopBlock_subtractposz = tofPara->Get("upStopBlock_subtractposz");
1392
1393 sideStopBlock_x = tofPara->Get("sideStopBlock_x");
1394 sideStopBlock_y = tofPara->Get("sideStopBlock_y");
1395 sideStopBlock_z = tofPara->Get("sideStopBlock_z");
1396 sideStopBlock_thickness = tofPara->Get("sideStopBlock_thickness");
1397 sideStopBlock_posz0 = tofPara->Get("sideStopBlock_posz0");
1398 sideStopBlock_posz1 = tofPara->Get("sideStopBlock_posz1");
1399
1400
1401 honeyComb_xdown = tofPara->Get("honeyComb_xdown");
1402 honeyComb_xup = tofPara->Get("honeyComb_xup");
1403 honeyComb_y = tofPara->Get("honeyComb_y");
1404 honeyComb_z = tofPara->Get("honeyComb_z");
1405 honeyComb_surface_y = tofPara->Get("honeyComb_surface_y");
1406 honeyComb_posz = tofPara->Get("honeyComb_posz");
1407
1408 tape_mrpc_xdown = tofPara->Get("tape_mrpc_xdown");
1409 tape_mrpc_xup = tofPara->Get("tape_mrpc_xup");
1410 tape_mrpc_y = tofPara->Get("tape_mrpc_y");
1411 tape_mrpc_z = tofPara->Get("tape_mrpc_z");
1412 tape_mrpc_posz = tofPara->Get("tape_mrpc_posz");
1413
1414 pcbBoard_xdown = tofPara->Get("pcbBoard_xdown");
1415 pcbBoard_xup = tofPara->Get("pcbBoard_xup");
1416 pcbBoard_y0 = tofPara->Get("pcbBoard_y0");
1417 pcbBoard_y1 = tofPara->Get("pcbBoard_y1");
1418 pcbBoard_z = tofPara->Get("pcbBoard_z");
1419
1420 strip_x[0] = tofPara->Get("strip_x[0]");
1421 strip_x[1] = tofPara->Get("strip_x[1]");
1422 strip_x[2] = tofPara->Get("strip_x[2]");
1423 strip_x[3] = tofPara->Get("strip_x[3]");
1424 strip_x[4] = tofPara->Get("strip_x[4]");
1425 strip_x[5] = tofPara->Get("strip_x[5]");
1426 strip_x[6] = tofPara->Get("strip_x[6]");
1427 strip_x[7] = tofPara->Get("strip_x[7]");
1428 strip_x[8] = tofPara->Get("strip_x[8]");
1429 strip_x[9] = tofPara->Get("strip_x[9]");
1430 strip_x[10] = tofPara->Get("strip_x[10]");
1431 strip_x[11] = tofPara->Get("strip_x[11]");
1432 strip_y = tofPara->Get("strip_y");
1433 strip_z = tofPara->Get("strip_z");
1434 strip_gap = tofPara->Get("strip_gap");
1435 strip_posz = tofPara->Get("strip_posz");
1436
1437 mylar_xdown = tofPara->Get("mylar_xdown");
1438 mylar_xup = tofPara->Get("mylar_xup");
1439 mylar_y = tofPara->Get("mylar_y");
1440 mylar_z = tofPara->Get("mylar_z");
1441 mylar_posz = tofPara->Get("mylar_posz");
1442
1443 glass0_xdown = tofPara->Get("glass0_xdown");
1444 glass0_xup = tofPara->Get("glass0_xup");
1445 glass0_y = tofPara->Get("glass0_y");
1446 glass0_z = tofPara->Get("glass0_z");
1447 glass0_posz = tofPara->Get("glass0_posz");
1448
1449 glass1_xdown = tofPara->Get("glass1_xdown");
1450 glass1_xup = tofPara->Get("glass1_xup");
1451 glass1_y = tofPara->Get("glass1_y");
1452 glass1_z = tofPara->Get("glass1_z");
1453 glass1_posz = tofPara->Get("glass1_posz");
1454
1455 gasLayer_xdown = tofPara->Get("gasLayer_xdown");
1456 gasLayer_xup = tofPara->Get("gasLayer_xup");
1457 gasLayer_y = tofPara->Get("gasLayer_y");
1458 gasLayer_z = tofPara->Get("gasLayer_z");
1459 gasLayer_posz = tofPara->Get("gasLayer_posz");
1460
1461 endcap_r = tofPara->Get("endcap_r");
1462 endcap_length = tofPara->Get("endcap_length");
1463 layer_posz = tofPara->Get("layer_posz");
1464 startAngle_east = tofPara->Get("startAngle_east");
1465 startAngle_west = tofPara->Get("startAngle_west");
1466}
x
Double_t x[10]
Definition: DataBase/tau_mode.c:57
tan
double tan(const BesAngle a)
Definition: InstallArea/include/MdcGeom/MdcGeom/BesAngle.h:216
sin
double sin(const BesAngle a)
Definition: InstallArea/include/MdcGeom/MdcGeom/BesAngle.h:210
cos
double cos(const BesAngle a)
Definition: InstallArea/include/MdcGeom/MdcGeom/BesAngle.h:213
smallL
G4double smallL
Definition: PartProduce.cc:19
rot_dummy
G4RotationMatrix rot_dummy(0 *deg, 0 *deg, 0 *deg)
checkOverlap
bool checkOverlap
Definition: PartProduce.cc:21
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
BesSubdetector::FindLogicalVolume
G4LogicalVolume * FindLogicalVolume(const G4String &vn)
Definition: InstallArea/include/SimUtil/SimUtil/BesSubdetector.hh:38
BesTofConstruction::getXYZ
void getXYZ(const G4RotationMatrix r, double &a, double &b, double &c) const
Definition: BesTofConstruction.cc:1279
BesTofConstruction::Construct
void Construct(G4LogicalVolume *logicalBes)
Definition: BesTofConstruction.cc:121
BesTofGeoParameter::Get
double Get(std::string key)
Definition: BesTofGeoParameter.cc:284
BesTofGeoParameter::GetInstance
static BesTofGeoParameter * GetInstance()
Definition: BesTofGeoParameter.cc:21
BesTofGeoParameter::GetVec
VEC GetVec(std::string key)
Definition: BesTofGeoParameter.cc:299
PartProduce::lg_sideFEE
G4LogicalVolume * lg_sideFEE(string name="")
Definition: PartProduce.cc:390
PartProduce::IniParam
void IniParam()
Definition: PartProduce.cc:34
PartProduce::lg_mylar
G4LogicalVolume * lg_mylar(string name="")
Definition: PartProduce.cc:732
PartProduce::lg_lowCover1
G4LogicalVolume * lg_lowCover1(string name="")
Definition: PartProduce.cc:270
PartProduce::lg_glass
G4LogicalVolume * lg_glass(G4int flag, string name="")
Definition: PartProduce.cc:747
PartProduce::lg_container
G4LogicalVolume * lg_container(int flag, string name="")
Definition: PartProduce.cc:485
PartProduce::lg_gasContainer
G4LogicalVolume * lg_gasContainer(string name="")
Definition: PartProduce.cc:589
PartProduce::lg_sideStopBlock
G4LogicalVolume * lg_sideStopBlock(string name="")
Definition: PartProduce.cc:848
PartProduce::lg_strip
G4AssemblyVolume * lg_strip(string name="")
Definition: PartProduce.cc:700
PartProduce::lg_honeyComb
G4LogicalVolume * lg_honeyComb(string name="")
Definition: PartProduce.cc:627
PartProduce::lg_containerFrame
G4LogicalVolume * lg_containerFrame(string name="")
Definition: PartProduce.cc:204
PartProduce::lg_lowCover2
G4LogicalVolume * lg_lowCover2(string name="")
Definition: PartProduce.cc:285
PartProduce::lg_pcbBoard
G4LogicalVolume * lg_pcbBoard(G4int flag, string name="")
Definition: PartProduce.cc:680
PartProduce::lg_bracket
G4LogicalVolume * lg_bracket(string name="")
Definition: PartProduce.cc:796
PartProduce::lg_bareChamber
G4LogicalVolume * lg_bareChamber(string name="")
Definition: PartProduce.cc:605
PartProduce::lg_upStopBlock
G4LogicalVolume * lg_upStopBlock(string name="")
Definition: PartProduce.cc:820
PartProduce::lg_tape_mrpc
G4LogicalVolume * lg_tape_mrpc(string name="")
Definition: PartProduce.cc:664
PartProduce::lg_gasLayer
G4LogicalVolume * lg_gasLayer(string name="")
Definition: PartProduce.cc:781
PartProduce::lg_castingDie
G4LogicalVolume * lg_castingDie(string name="")
Definition: PartProduce.cc:460
PartProduce::lg_upFEE
G4LogicalVolume * lg_upFEE(string name="")
Definition: PartProduce.cc:321
PartProduce::lg_upCover
G4LogicalVolume * lg_upCover(string name="")
Definition: PartProduce.cc:255
SubDetectorG4Geo::GetTopVolume
G4LogicalVolume * GetTopVolume()
Get the top(world) volume;.
Definition: DetectorDescription/G4Geo/G4Geo-00-00-13/G4Geo/SubDetectorG4Geo.h:50
H
IMPLICIT REAL *A H
Definition: myXsection.h:1