Loading...
Searching...
No Matches
BesEmcConstruction Class Reference
#include <BesEmcConstruction.hh>
Inheritance diagram for BesEmcConstruction:
Public Member Functions | |
| BesEmcConstruction () | |
| ~BesEmcConstruction () | |
| void | SetVerboseLevel (G4int val) |
| void | SetCrystalMaterial (G4String) |
| void | SetCasingMaterial (G4String) |
| void | SetCasingThickness (G4ThreeVector) |
| void | SetBSCRmin (G4double) |
| void | SetBSCNbPhi (G4int) |
| void | SetBSCNbTheta (G4int) |
| void | SetStartIDTheta (G4int) |
| void | SetBSCCrystalLength (G4double) |
| void | SetBSCYFront0 (G4double) |
| void | SetBSCYFront (G4double) |
| void | SetBSCPosition0 (G4double) |
| void | SetBSCPosition1 (G4double) |
| void | SetMagField (G4double) |
| void | Construct (G4LogicalVolume *) |
| void | UpdateGeometry () |
| void | ThreeVectorTrans (G4ThreeVector fPnt[8], double x[8], double y[8], double z[8]) |
| void | TransformToArb8 (const G4ThreeVector fPnt[8], G4ThreeVector newfPnt[8], G4ThreeVector ¢er, G4ThreeVector &rotAngle) |
| Hep3Vector | RotAngleFromNewZ (Hep3Vector newZ) |
| void | PrintEMCParameters () |
| G4int | GetVerboseLevel () |
| G4double | GetMagField () |
| G4int | GetStartIDTheta () |
| G4Material * | GetCrystalMaterial () |
| G4Material * | GetCasingMaterial () |
| const G4VPhysicalVolume * | GetEMC () |
| const G4VPhysicalVolume * | GetBSCPhi () |
| const G4VPhysicalVolume * | GetBSCTheta () |
| const G4VPhysicalVolume * | GetBSCCrystal () |
| const G4VPVParameterisation * | GetCrystalParam () |
| void | ConstructSPFrame (G4LogicalVolume *, BesEmcGeometry *) |
| void | ConstructEndGeometry (G4LogicalVolume *) |
| void | GetLogicalVolume () |
| void | SetVisAndSD () |
| G4int | ComputeEndCopyNb (G4int) |
| BesEmcConstruction () | |
| ~BesEmcConstruction () | |
| void | SetVerboseLevel (G4int val) |
| void | SetCrystalMaterial (G4String) |
| void | SetCasingMaterial (G4String) |
| void | SetCasingThickness (G4ThreeVector) |
| void | SetBSCRmin (G4double) |
| void | SetBSCNbPhi (G4int) |
| void | SetBSCNbTheta (G4int) |
| void | SetStartIDTheta (G4int) |
| void | SetBSCCrystalLength (G4double) |
| void | SetBSCYFront0 (G4double) |
| void | SetBSCYFront (G4double) |
| void | SetBSCPosition0 (G4double) |
| void | SetBSCPosition1 (G4double) |
| void | SetMagField (G4double) |
| void | Construct (G4LogicalVolume *) |
| void | UpdateGeometry () |
| void | ThreeVectorTrans (G4ThreeVector fPnt[8], double x[8], double y[8], double z[8]) |
| void | TransformToArb8 (const G4ThreeVector fPnt[8], G4ThreeVector newfPnt[8], G4ThreeVector ¢er, G4ThreeVector &rotAngle) |
| Hep3Vector | RotAngleFromNewZ (Hep3Vector newZ) |
| void | PrintEMCParameters () |
| G4int | GetVerboseLevel () |
| G4double | GetMagField () |
| G4int | GetStartIDTheta () |
| G4Material * | GetCrystalMaterial () |
| G4Material * | GetCasingMaterial () |
| const G4VPhysicalVolume * | GetEMC () |
| const G4VPhysicalVolume * | GetBSCPhi () |
| const G4VPhysicalVolume * | GetBSCTheta () |
| const G4VPhysicalVolume * | GetBSCCrystal () |
| const G4VPVParameterisation * | GetCrystalParam () |
| void | ConstructSPFrame (G4LogicalVolume *, BesEmcGeometry *) |
| void | ConstructEndGeometry (G4LogicalVolume *) |
| void | GetLogicalVolume () |
| void | SetVisAndSD () |
| G4int | ComputeEndCopyNb (G4int) |
| BesEmcConstruction () | |
| ~BesEmcConstruction () | |
| void | SetVerboseLevel (G4int val) |
| void | SetCrystalMaterial (G4String) |
| void | SetCasingMaterial (G4String) |
| void | SetCasingThickness (G4ThreeVector) |
| void | SetBSCRmin (G4double) |
| void | SetBSCNbPhi (G4int) |
| void | SetBSCNbTheta (G4int) |
| void | SetStartIDTheta (G4int) |
| void | SetBSCCrystalLength (G4double) |
| void | SetBSCYFront0 (G4double) |
| void | SetBSCYFront (G4double) |
| void | SetBSCPosition0 (G4double) |
| void | SetBSCPosition1 (G4double) |
| void | SetMagField (G4double) |
| void | Construct (G4LogicalVolume *) |
| void | UpdateGeometry () |
| void | ThreeVectorTrans (G4ThreeVector fPnt[8], double x[8], double y[8], double z[8]) |
| void | TransformToArb8 (const G4ThreeVector fPnt[8], G4ThreeVector newfPnt[8], G4ThreeVector ¢er, G4ThreeVector &rotAngle) |
| Hep3Vector | RotAngleFromNewZ (Hep3Vector newZ) |
| void | PrintEMCParameters () |
| G4int | GetVerboseLevel () |
| G4double | GetMagField () |
| G4int | GetStartIDTheta () |
| G4Material * | GetCrystalMaterial () |
| G4Material * | GetCasingMaterial () |
| const G4VPhysicalVolume * | GetEMC () |
| const G4VPhysicalVolume * | GetBSCPhi () |
| const G4VPhysicalVolume * | GetBSCTheta () |
| const G4VPhysicalVolume * | GetBSCCrystal () |
| const G4VPVParameterisation * | GetCrystalParam () |
| void | ConstructSPFrame (G4LogicalVolume *, BesEmcGeometry *) |
| void | ConstructEndGeometry (G4LogicalVolume *) |
| void | GetLogicalVolume () |
| void | SetVisAndSD () |
| G4int | ComputeEndCopyNb (G4int) |
| Public Member Functions inherited from BesSubdetector | |
| BesSubdetector () | |
| virtual | ~BesSubdetector () |
| G4LogicalVolume * | FindLogicalVolume (const G4String &vn) |
| BesSubdetector () | |
| virtual | ~BesSubdetector () |
| G4LogicalVolume * | FindLogicalVolume (const G4String &vn) |
| BesSubdetector () | |
| virtual | ~BesSubdetector () |
| G4LogicalVolume * | FindLogicalVolume (const G4String &vn) |
Static Public Member Functions | |
| static BesEmcConstruction * | GetBesEmcConstruction () |
| static BesEmcConstruction * | GetBesEmcConstruction () |
| static BesEmcConstruction * | GetBesEmcConstruction () |
Additional Inherited Members | |
| Protected Attributes inherited from BesSubdetector | |
| SAXProcessor | m_sxp |
| ProcessingConfigurator | m_config |
Detailed Description
Definition at line 51 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
Constructor & Destructor Documentation
◆ BesEmcConstruction() [1/3]
| BesEmcConstruction::BesEmcConstruction | ( | ) |
Definition at line 70 of file BesEmcConstruction.cc.
71 : verboseLevel( 0 )
72 , solidEMC( 0 )
73 , logicEMC( 0 )
74 , physiEMC( 0 )
75 , logicBSCWorld( 0 )
76 , solidBSCPhi( 0 )
77 , logicBSCPhi( 0 )
78 , physiBSCPhi( 0 )
79 , solidBSCTheta( 0 )
80 , logicBSCTheta( 0 )
81 , physiBSCTheta( 0 )
82 , solidBSCCrystal( 0 )
83 , logicBSCCrystal( 0 )
84 , physiBSCCrystal( 0 )
85 , magField( 0 )
86 , detectorMessenger( 0 )
87 , besEMCSD( 0 )
88 , crystalParam( 0 )
89 , logicEnd( 0 )
90 , logicEndPhi( 0 )
91 , logicEndCasing( 0 )
92 , logicEndCrystal( 0 )
93 , logicRear( 0 )
94 , logicRearCasing( 0 )
95 , logicOrgGlass( 0 )
96 , logicPD( 0 )
97 , logicAlPlate( 0 )
98 , logicPreAmpBox( 0 )
99 , logicAirInPABox( 0 )
100 , logicHangingPlate( 0 )
101 , logicOCGirder( 0 )
102 , logicCable( 0 )
103 , logicWaterPipe( 0 )
104 , logicSupportBar( 0 )
105 , logicSupportBar1( 0 )
106 , logicEndRing( 0 )
107 , logicGear( 0 )
108 , logicTaperRing1( 0 )
109 , logicTaperRing2( 0 )
110 , logicTaperRing3( 0 ) {
111 if ( fBesEmcConstruction )
112 {
113 G4cout << "BesEmcConstruction constructed twice." << G4endl;
114 exit( -1 );
115 }
116 fBesEmcConstruction = this;
117 // for debug
118 // G4Exception("BesEmcConstruction::BesEmcConstruction() starting........");
119 startID = 1;
120 phiNbCrystals = 0;
121 thetaNbCrystals = 0;
122 besEMCGeometry = new BesEmcGeometry();
123 emcEnd = new BesEmcEndGeometry();
124}
Referenced by GetBesEmcConstruction().
◆ ~BesEmcConstruction() [1/3]
| BesEmcConstruction::~BesEmcConstruction | ( | ) |
Definition at line 126 of file BesEmcConstruction.cc.
126 {
127 if ( detectorMessenger ) delete detectorMessenger;
128 if ( crystalParam ) delete crystalParam;
129 if ( besEMCGeometry ) delete besEMCGeometry;
130 if ( emcEnd ) delete emcEnd;
131
132 BesEmcParameter::Kill();
133}
◆ BesEmcConstruction() [2/3]
| BesEmcConstruction::BesEmcConstruction | ( | ) |
◆ ~BesEmcConstruction() [2/3]
| BesEmcConstruction::~BesEmcConstruction | ( | ) |
◆ BesEmcConstruction() [3/3]
| BesEmcConstruction::BesEmcConstruction | ( | ) |
◆ ~BesEmcConstruction() [3/3]
| BesEmcConstruction::~BesEmcConstruction | ( | ) |
Member Function Documentation
◆ ComputeEndCopyNb() [1/3]
| G4int BesEmcConstruction::ComputeEndCopyNb | ( | G4int | num | ) |
Definition at line 1181 of file BesEmcConstruction.cc.
1181 {
1182 G4int copyNb;
1184 {
1185 case 30: copyNb = 5; break;
1186 case 31: copyNb = 6; break;
1187 case 32: copyNb = 14; break;
1188 case 33: copyNb = 15; break;
1189 case 34: copyNb = 16; break;
1191 }
1192 return copyNb;
1193}
◆ ComputeEndCopyNb() [2/3]
| G4int BesEmcConstruction::ComputeEndCopyNb | ( | G4int | ) |
◆ ComputeEndCopyNb() [3/3]
| G4int BesEmcConstruction::ComputeEndCopyNb | ( | G4int | ) |
◆ Construct() [1/3]
|
virtual |
Implements BesSubdetector.
Definition at line 135 of file BesEmcConstruction.cc.
135 {
136 gErrorIgnoreLevel = kFatal;
137 besEMCGeometry->ComputeEMCParameters();
138 detectorMessenger = new BesEmcDetectorMessenger( this, besEMCGeometry );
139 emcEnd->ComputeParameters();
140
141 G4SDManager* SDman = G4SDManager::GetSDMpointer();
142 if ( !besEMCSD )
143 {
144 besEMCSD = new BesEmcSD( "CalorSD", this, besEMCGeometry );
145 SDman->AddNewDetector( besEMCSD );
146 }
147
148 // Construction
151 {
153
154 if ( logicEMC )
155 {
156 physiEMC = new G4PVPlacement( 0, G4ThreeVector( 0.0, 0.0, 0.0 ), logicEMC, "physicalEMC",
158 G4cout << "logicEmc: === " << logicEMC << " physiEmc " << physiEMC << G4endl;
159
160 GetLogicalVolume();
161 SetVisAndSD();
162 }
163 }
164 else
165 {
166 // for debug
167 // G4Exception("BesEmcConstruction::Construct() starting............");
168 //
169 DefineMaterials();
170 phiNbCrystals = ( *besEMCGeometry ).BSCNbPhi;
171 thetaNbCrystals = ( *besEMCGeometry ).BSCNbTheta * 2;
172
173 G4double da = 0.001 * deg; // delta angle to avoid overlap
174
175 //
176 // BSC
177 //
178 solidBSC = new G4Tubs(
179 "solidBSC", ( *besEMCGeometry ).TaperRingRmin1,
180 ( *besEMCGeometry ).BSCRmax + ( *besEMCGeometry ).SPBarThickness +
181 ( *besEMCGeometry ).SPBarThickness1 + 2.1 * mm, // radius from 942mm to 940 mm
182 ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 +
183 ( *besEMCGeometry ).EndRingDz,
184 0. * deg, 360. * deg );
185
186 solidESC = new G4Cons( "solidESC", ( *emcEnd ).WorldRmin1, ( *emcEnd ).WorldRmax1,
187 ( *emcEnd ).WorldRmin2, ( *emcEnd ).WorldRmax2,
188 ( *emcEnd ).WorldDz / 2, 0. * deg, 360. * deg );
189
190 solidEMC = new G4UnionSolid( "solidEMC0", solidBSC, solidESC, 0,
191 G4ThreeVector( 0, 0, ( *emcEnd ).WorldZPosition ) );
192
193 G4RotationMatrix* rotateESC = new G4RotationMatrix();
194 rotateESC->rotateY( 180. * deg );
195
196 solidEMC = new G4UnionSolid( "solidEMC", solidEMC, solidESC, rotateESC,
197 G4ThreeVector( 0, 0, -( *emcEnd ).WorldZPosition ) );
198
199 logicEMC = new G4LogicalVolume( solidEMC, G4Material::GetMaterial( "Air" ), "logicalEMC" );
200
201 physiEMC = new G4PVPlacement( 0, G4ThreeVector(), logicEMC, "physicalEMC", logicBes, false,
202 0, CHECKLV0 );
203
204 solidBSCWorld =
205 new G4SubtractionSolid( "solidBSCWorld0", solidBSC, solidESC, 0,
206 G4ThreeVector( 0, 0, ( *emcEnd ).WorldZPosition ) );
207
208 solidBSCWorld =
209 new G4SubtractionSolid( "solidBSCWorld", solidBSCWorld, solidESC, rotateESC,
210 G4ThreeVector( 0, 0, -( *emcEnd ).WorldZPosition ) );
211
212 logicBSCWorld = new G4LogicalVolume( solidBSCWorld, G4Material::GetMaterial( "Air" ),
213 "logicalBSCWorld" );
214
215 G4RotationMatrix* rotBSC = new G4RotationMatrix();
216 rotBSC->rotateY( 180. * deg );
217 physiBSCWorld = new G4PVPlacement( rotBSC, G4ThreeVector(), logicBSCWorld,
219
220 G4RotationMatrix* rotateMatrix[200];
221 G4double oOp, ox, oy, oz;
222 G4double delta = 0 * deg;
223 G4ThreeVector axis = G4ThreeVector( 0, 0, 0 );
224 oOp = ( *besEMCGeometry ).BSCRmin /
225 sin( 0.5 * ( *besEMCGeometry ).BSCPhiDphi + 90 * deg ) *
226 sin( ( *besEMCGeometry ).BSCAngleRotat );
227 G4double ll = ( *besEMCGeometry ).BSCCryLength;
228 G4double rr = ( *besEMCGeometry ).BSCRmin;
229 G4double oj = sqrt( ll * ll + rr * rr -
230 2 * ll * rr * cos( 180. * deg - ( *besEMCGeometry ).BSCAngleRotat ) );
231 G4double oij =
232 90. * deg - ( *besEMCGeometry ).BSCPhiDphi / 2. - ( *besEMCGeometry ).BSCAngleRotat;
233 G4double doj = asin( sin( 180. * deg - ( *besEMCGeometry ).BSCAngleRotat ) / oj * ll );
234 G4double ioj = ( *besEMCGeometry ).BSCPhiDphi / 2. + doj;
236 G4double dOp = rr / sin( 90. * deg - ( *besEMCGeometry ).BSCPhiDphi / 2. ) *
237 sin( 90. * deg + ( *besEMCGeometry ).BSCPhiDphi / 2. -
238 ( *besEMCGeometry ).BSCAngleRotat );
239 G4double cOp = rr / sin( 90. * deg + ( *besEMCGeometry ).BSCPhiDphi / 2. ) *
240 sin( 90. * deg - ( *besEMCGeometry ).BSCPhiDphi / 2. -
241 ( *besEMCGeometry ).BSCAngleRotat );
242 G4double ch = ( dOp + ll ) / cos( ( *besEMCGeometry ).BSCPhiDphi ) - cOp;
243 G4double hi = ( dOp + ll ) * tan( ( *besEMCGeometry ).BSCPhiDphi ) - ij;
244 G4double oh = sqrt( ch * ch + rr * rr -
245 2 * ch * rr * cos( 180 * deg - ( *besEMCGeometry ).BSCAngleRotat ) );
246 G4double hoi = asin( sin( 180 * deg - oij ) / oh * hi );
247 G4double dok = asin( sin( 180 * deg - ( *besEMCGeometry ).BSCAngleRotat ) / oh * ch );
248 if ( verboseLevel > 3 )
249 G4cout << "oj=" << oj / cm << G4endl << "oij=" << oij / deg << G4endl
250 << "doj=" << doj / deg << G4endl << "ioj=" << ioj / deg << G4endl
251 << "ij=" << ij / cm << G4endl << "dOp=" << dOp / cm << G4endl
252 << "cOp=" << cOp / cm << G4endl << "ch=" << ch / cm << G4endl << "hi=" << hi / cm
253 << G4endl << "oh=" << oh / cm << G4endl << "hoi=" << hoi / deg << G4endl
254 << "dok=" << dok / deg << G4endl;
255
256 // Phi Cell
257 G4double cmo = asin( sin( 180 * degree - ( *besEMCGeometry ).BSCAngleRotat ) /
258 ( *besEMCGeometry ).BSCRmax2 * dOp );
259 G4double cm = ( *besEMCGeometry ).BSCRmax2 /
260 sin( 180 * degree - ( *besEMCGeometry ).BSCAngleRotat ) *
261 sin( ( *besEMCGeometry ).BSCAngleRotat - cmo );
262 G4ThreeVector Pout( dOp + cm * cos( ( *besEMCGeometry ).BSCAngleRotat ),
263 -cm * sin( ( *besEMCGeometry ).BSCAngleRotat ),
264 ( *besEMCGeometry ).BSCDz );
265
266 G4double rTaperRingOuter1 =
267 ( *besEMCGeometry ).TaperRingRmin1 + ( *besEMCGeometry ).TaperRingThickness1;
268 G4double rTaperRingOuter2 =
269 ( *besEMCGeometry ).TaperRingRmin2 + ( *besEMCGeometry ).TaperRingDr;
270 G4double zTaperRing1 = ( *besEMCGeometry ).BSCDz +
271 ( *besEMCGeometry ).TaperRingThickness3 -
272 ( *besEMCGeometry ).TaperRingDz;
273 G4double zTaperRing2 = ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3;
274
275 G4RotationMatrix* rotIntersection = new G4RotationMatrix();
276 rotIntersection->rotateZ( -( *besEMCGeometry ).BSCAngleRotat -
277 ( *besEMCGeometry ).BSCPhiDphi / 2. - hoi );
279 oOp * sin( -90. * deg + ( *besEMCGeometry ).BSCAngleRotat + hoi ), 0 );
280 G4AffineTransform Td( rotIntersection, P );
281 G4ThreeVector md = Td.Inverse().TransformPoint( G4ThreeVector( 0, 0, 0 ) );
282 G4ThreeVector vd = Td.Inverse().TransformPoint( Pout );
283
284 G4double zPlane0[4] = { -( *besEMCGeometry ).BSCDz, -zTaperRing1, zTaperRing1,
285 ( *besEMCGeometry ).BSCDz };
286 G4double rInner0[4] = { vd.perp(), cOp, cOp, vd.perp() };
287 G4double rOuter0[4] = { ( *besEMCGeometry ).BSCPhiRmax, ( *besEMCGeometry ).BSCPhiRmax,
288 ( *besEMCGeometry ).BSCPhiRmax, ( *besEMCGeometry ).BSCPhiRmax };
289
290 G4double zPlane1[4] = { -zTaperRing2, -zTaperRing1, zTaperRing1, zTaperRing2 };
291 G4double rInner1[4] = { rTaperRingOuter2, rTaperRingOuter1, rTaperRingOuter1,
292 rTaperRingOuter2 };
293 G4double rOuter1[4] = { ( *besEMCGeometry ).BSCRmax, ( *besEMCGeometry ).BSCRmax,
294 ( *besEMCGeometry ).BSCRmax, ( *besEMCGeometry ).BSCRmax };
295
296 G4VSolid* solidBSCPhi1 =
297 new G4Polycone( "solidBSCPhi1", 0, 360 * deg, 4, zPlane1, rInner1, rOuter1 );
298 G4VSolid* solidBSCPhi0 =
299 new G4Polyhedra( "solidBSCPhi0", 360. * deg - ( *besEMCGeometry ).BSCPhiDphi,
300 ( *besEMCGeometry ).BSCPhiDphi, ( *besEMCGeometry ).BSCNbPhi / 2, 4,
301 zPlane0, rInner0, rOuter0 );
302 G4IntersectionSolid* solidBSCPhi =
303 new G4IntersectionSolid( "solidBSCPhi", solidBSCPhi0, solidBSCPhi1, 0, md );
304
305 logicBSCPhi =
306 new G4LogicalVolume( solidBSCPhi, G4Material::GetMaterial( "Air" ), "logicalBSCPhi" );
307
308 G4int i;
309 for ( G4int j = 0; j < ( *besEMCGeometry ).BSCNbPhi; j++ ) //=============
310 {
311 if ( j < ( *besEMCGeometry ).BSCNbPhi / 2 )
312 { // 0~59
313 i = ( *besEMCGeometry ).BSCNbPhi / 2 - j - 1;
314 }
315 else
316 { // 60~119
317 i = ( *besEMCGeometry ).BSCNbPhi * 3 / 2 - j - 1;
318 }
319 rotateMatrix[i] = new G4RotationMatrix();
320 rotateMatrix[i]->rotateZ( -i * ( *besEMCGeometry ).BSCPhiDphi -
321 ( *besEMCGeometry ).BSCAngleRotat -
322 ( *besEMCGeometry ).BSCPhiDphi / 2. - hoi );
323 rotateMatrix[i]->getAngleAxis( delta, axis );
324 // G4cout << "The axis of crystals in the world system is: "
325 // << delta/deg << "(deg)(delta) "
326 //<< axis << "(Z axis)"<< G4endl;
327 ox = oOp * cos( -90. * deg + ( *besEMCGeometry ).BSCAngleRotat + hoi +
328 i * ( *besEMCGeometry ).BSCPhiDphi );
329 oy = oOp * sin( -90. * deg + ( *besEMCGeometry ).BSCAngleRotat + hoi +
330 i * ( *besEMCGeometry ).BSCPhiDphi );
331 oz = 0 * cm;
332
333 ostringstream strPhi;
334 strPhi << "physicalBSCPhi" << j;
335
336 physiBSCPhi =
337 new G4PVPlacement( rotateMatrix[i], G4ThreeVector( ox, oy, oz ), logicBSCPhi,
339 // G4cout << G4ThreeVector(ox/cm,oy/cm,oz/cm) <<"(cm)" << G4endl
340 // << (-(*besEMCGeometry).BSCAngleRotat+(i-1)*(*besEMCGeometry).BSCPhiDphi)/deg
341 // <<"(degree)" << G4endl;
342 }
343
344 //
345 // Crystals
346 //
347 G4double zHalfLength[50];
348 G4double thetaAxis[50];
349 G4double phiAxis[50];
350 G4double yHalfLength1[50];
351 G4double xHalfLength2[50];
352 G4double xHalfLength1[50];
353 G4double tanAlpha1[50];
354 G4double yHalfLength2[50];
355 G4double xHalfLength4[50];
356 G4double xHalfLength3[50];
357 G4double tanAlpha2[50];
358 G4double xPosition[50];
359 G4double yPosition[50];
360 G4double zPosition[50];
361 G4double thetaPosition[50];
362 for ( i = 0; i < ( *besEMCGeometry ).BSCNbTheta; i++ )
363 {
364 zHalfLength[i] = ( *besEMCGeometry ).zHalfLength[i];
365 thetaAxis[i] = ( *besEMCGeometry ).thetaAxis[i];
366 phiAxis[i] = ( *besEMCGeometry ).phiAxis[i];
367 yHalfLength1[i] = ( *besEMCGeometry ).yHalfLength1[i];
368 xHalfLength2[i] = ( *besEMCGeometry ).xHalfLength2[i];
369 xHalfLength1[i] = ( *besEMCGeometry ).xHalfLength1[i];
370 tanAlpha1[i] = ( *besEMCGeometry ).tanAlpha1[i];
371 yHalfLength2[i] = ( *besEMCGeometry ).yHalfLength2[i];
372 xHalfLength4[i] = ( *besEMCGeometry ).xHalfLength4[i];
373 xHalfLength3[i] = ( *besEMCGeometry ).xHalfLength3[i];
374 tanAlpha2[i] = ( *besEMCGeometry ).tanAlpha2[i];
375 xPosition[i] = ( *besEMCGeometry ).xPosition[i];
376 yPosition[i] = ( *besEMCGeometry ).yPosition[i];
377 zPosition[i] = ( *besEMCGeometry ).zPosition[i];
378 thetaPosition[i] = ( *besEMCGeometry ).thetaPosition[i];
379 if ( verboseLevel > 4 )
380 G4cout << "The sizes of the " << i + 1 << " crystal are:" << G4endl
381 << "zHalfLength =" << zHalfLength[i] / cm << "(cm)," << G4endl
382 << "thetaAxis =" << thetaAxis[i] / deg << "(deg)," << G4endl
383 << "phiAxis =" << phiAxis[i] / deg << "(deg)," << G4endl
384 << "yHalfLength1=" << yHalfLength1[i] / cm << "(cm)," << G4endl
385 << "xHalfLength1=" << xHalfLength1[i] / cm << "(cm)," << G4endl
386 << "xHalfLength2=" << xHalfLength2[i] / cm << "(cm)," << G4endl
387 << "tanAlpha1 =" << tanAlpha1[i] << G4endl
388 << "yHalfLength2=" << yHalfLength2[i] / cm << "(cm)," << G4endl
389 << "xHalfLength3=" << xHalfLength3[i] / cm << "(cm)," << G4endl
390 << "xHalfLength4=" << xHalfLength4[i] / cm << "(cm)," << G4endl
391 << "tanAlpha2 =" << tanAlpha2[i] << "." << G4endl;
392 }
393 besEMCGeometry->ModifyForCasing();
394
395 solidBSCCrystal =
396 new G4Trap( "solidCrystal", 100 * cm, 100 * deg, 100 * deg, 100 * cm, 100 * cm,
397 100 * cm, 100 * deg, 100 * cm, 100 * cm, 100 * cm, 100 * deg );
398
399 logicBSCCrystal =
400 new G4LogicalVolume( solidBSCCrystal, fCrystalMaterial, "logicalCrystal" );
401
402 crystalParam = new BesCrystalParameterisation( startID, thetaNbCrystals,
403 ( *besEMCGeometry ).BSCNbTheta * 2,
404 besEMCGeometry, verboseLevel );
405
406 //---------------------------------------------------------------------------------
407 // rear substance
408 solidRear =
409 new G4Box( "solidRearBox", ( *besEMCGeometry ).rearBoxLength / 2,
410 ( *besEMCGeometry ).rearBoxLength / 2, ( *besEMCGeometry ).rearBoxDz / 2 );
411
412 logicRear =
413 new G4LogicalVolume( solidRear, G4Material::GetMaterial( "Air" ), "logicalRearBox" );
414
415 // organic glass
416 solidOrgGlass = new G4Box( "solidOrganicGlass", ( *besEMCGeometry ).orgGlassLengthX / 2,
417 ( *besEMCGeometry ).orgGlassLengthY / 2,
418 ( *besEMCGeometry ).orgGlassLengthZ / 2 );
419
420 logicOrgGlass = new G4LogicalVolume( solidOrgGlass, organicGlass, "logicalOrganicGlass" );
421
422 physiOrgGlass = new G4PVPlacement(
423 0,
424 G4ThreeVector(
425 0, 0,
426 -( ( *besEMCGeometry ).rearBoxDz - ( *besEMCGeometry ).orgGlassLengthZ ) / 2 ),
428
429 // casing
430 solidCasingBox = new G4Box( "solidCasingBox", ( *besEMCGeometry ).rearBoxLength / 2,
431 ( *besEMCGeometry ).rearBoxLength / 2,
432 ( *besEMCGeometry ).rearCasingThickness / 2 );
433
434 solidAirHole =
435 new G4Box( "solidAirHole", ( *besEMCGeometry ).orgGlassLengthX / 2,
436 ( *besEMCGeometry ).orgGlassLengthY / 2,
437 ( *besEMCGeometry ).rearBoxDz / 2 ); // any value more than casing thickness
438
439 solidRearCasing = new G4SubtractionSolid( "solidRearCasing", solidCasingBox, solidAirHole,
440 0, G4ThreeVector() );
441
442 logicRearCasing =
443 new G4LogicalVolume( solidRearCasing, rearCasingMaterial, "logicalRearCasing" );
444
445 physiRearCasing = new G4PVPlacement(
446 0,
447 G4ThreeVector(
448 0, 0,
449 -( ( *besEMCGeometry ).rearBoxDz - ( *besEMCGeometry ).rearCasingThickness ) / 2 ),
451
452 // Al Plate
453 solidAlBox =
454 new G4Box( "solidAlBox", ( *besEMCGeometry ).rearBoxLength / 2,
455 ( *besEMCGeometry ).rearBoxLength / 2, ( *besEMCGeometry ).AlPlateDz / 2 );
456
457 solidAlPlate =
458 new G4SubtractionSolid( "solidAlPlate", solidAlBox, solidAirHole, 0, G4ThreeVector() );
459
460 logicAlPlate = new G4LogicalVolume( solidAlPlate, G4Material::GetMaterial( "Aluminium" ),
461 "logicalAlPlate" );
462
463 physiAlPlate =
464 new G4PVPlacement( 0,
465 G4ThreeVector( 0, 0,
466 -( ( *besEMCGeometry ).rearBoxDz / 2 -
467 ( *besEMCGeometry ).rearCasingThickness -
468 ( *besEMCGeometry ).AlPlateDz / 2 ) ),
470
471 // photodiode
472 solidPD =
473 new G4Box( "solidPD",
474 ( *besEMCGeometry ).PDLengthX, // two PD
475 ( *besEMCGeometry ).PDLengthY / 2, ( *besEMCGeometry ).PDLengthZ / 2 );
476
477 logicPD =
478 new G4LogicalVolume( solidPD, G4Material::GetMaterial( "M_Silicon" ), "logicalPD" );
479
480 physiPD = new G4PVPlacement( 0,
481 G4ThreeVector( 0, 0,
482 -( ( *besEMCGeometry ).rearBoxDz / 2 -
483 ( *besEMCGeometry ).orgGlassLengthZ -
484 ( *besEMCGeometry ).PDLengthZ / 2 ) ),
486
487 // preamplifier box
488 solidPreAmpBox =
489 new G4Box( "solidPreAmpBox", ( *besEMCGeometry ).rearBoxLength / 2,
490 ( *besEMCGeometry ).rearBoxLength / 2, ( *besEMCGeometry ).PABoxDz / 2 );
491
492 logicPreAmpBox = new G4LogicalVolume(
493 solidPreAmpBox, G4Material::GetMaterial( "Aluminium" ), "logicalPreAmpBox" );
494
495 physiPreAmpBox = new G4PVPlacement(
496 0,
497 G4ThreeVector( 0, 0,
498 -( ( *besEMCGeometry ).rearBoxDz / 2 -
499 ( *besEMCGeometry ).rearCasingThickness -
500 ( *besEMCGeometry ).AlPlateDz - ( *besEMCGeometry ).PABoxDz / 2 ) ),
502
503 // air in preamplifier box
504 solidAirInPABox =
505 new G4Box( "solidAirInPABox",
506 ( *besEMCGeometry ).rearBoxLength / 2 - ( *besEMCGeometry ).PABoxThickness,
507 ( *besEMCGeometry ).rearBoxLength / 2 - ( *besEMCGeometry ).PABoxThickness,
508 ( *besEMCGeometry ).PABoxDz / 2 - ( *besEMCGeometry ).PABoxThickness );
509
510 logicAirInPABox = new G4LogicalVolume( solidAirInPABox, G4Material::GetMaterial( "Air" ),
511 "logicalAirInPABox" );
512
513 physiAirInPABox =
514 new G4PVPlacement( 0, G4ThreeVector(), logicAirInPABox, "physicalAirInPABox",
516
517 // stainless steel for hanging the crystal
518 solidHangingPlate = new G4Box( "solidHangingPlate", ( *besEMCGeometry ).rearBoxLength / 2,
519 ( *besEMCGeometry ).rearBoxLength / 2,
520 ( *besEMCGeometry ).HangingPlateDz / 2 );
521
522 logicHangingPlate =
523 new G4LogicalVolume( solidHangingPlate, stainlessSteel, "logicalHangingPlate" );
524
525 physiHangingPlate = new G4PVPlacement(
526 0,
527 G4ThreeVector( 0, 0,
528 -( ( *besEMCGeometry ).rearBoxDz / 2 -
529 ( *besEMCGeometry ).rearCasingThickness -
530 ( *besEMCGeometry ).AlPlateDz - ( *besEMCGeometry ).PABoxDz -
531 ( *besEMCGeometry ).HangingPlateDz / 2 ) ),
533
534 // water pipe
535 solidWaterPipe = new G4Tubs( "solidWaterPipe", 0, ( *besEMCGeometry ).waterPipeDr,
536 ( *besEMCGeometry ).BSCDz, 0. * deg, 360. * deg );
537
538 logicWaterPipe = new G4LogicalVolume( solidWaterPipe, waterPipe, "logicalWaterPipe" );
539
540 physiWaterPipe = new G4PVPlacement(
541 0,
542 G4ThreeVector( ( *besEMCGeometry ).cablePosX[0] - 2 * ( *besEMCGeometry ).cableDr,
543 ( *besEMCGeometry ).cablePosY[0] - ( *besEMCGeometry ).cableDr -
544 ( *besEMCGeometry ).waterPipeDr,
545 0 ),
547 //---------------------------------------------------------------------------------
548
549 //
550 // Theta Cell
551 //
552 G4String nameCrystalAndCasing = "CrystalAndCasing";
553
554 G4int id = 0; // ID of crystals after distinguishing left and right
555 for ( i = startID; i <= thetaNbCrystals; i++ ) //================
556 {
557 ostringstream strSolidCasing;
558 strSolidCasing << "solidBSCCasing" << i - 1;
559 ostringstream strVolumeCasing;
560 strVolumeCasing << "logicalBSCCasing" << i - 1;
561 ostringstream strPhysiCasing;
562 strPhysiCasing << "physicalBSCCasing" << i - 1;
563
564 if ( i > ( *besEMCGeometry ).BSCNbTheta )
565 {
566 id = i - ( *besEMCGeometry ).BSCNbTheta - 1;
567 solidBSCTheta =
568 new G4Trap( strSolidCasing.str(), zHalfLength[id], thetaAxis[id], -phiAxis[id],
569 yHalfLength1[id], xHalfLength2[id], xHalfLength1[id], -tanAlpha1[id],
570 yHalfLength2[id], xHalfLength4[id], xHalfLength3[id], -tanAlpha2[id] );
571
572 // G4cout<<"in EmcConstr1: "<<strSolidCasing.str()<<" x1="<<xHalfLength1[id]<<"
573 // y1="<<yHalfLength1[id]<<" theta="<<thetaAxis[id]
574 //<<" phi="<<-phiAxis[id]<<" a1="<<-tanAlpha1[id]<<G4endl;
575
576 logicBSCTheta =
577 new G4LogicalVolume( solidBSCTheta, fCasingMaterial, strVolumeCasing.str() );
578
579 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1] = new G4RotationMatrix();
580 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1]->rotateZ( -90 * deg );
581 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1]->rotateX( -thetaPosition[id] );
582
583 physiBSCTheta = new G4PVPlacement(
584 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1],
585 G4ThreeVector( xPosition[id], yPosition[id], zPosition[id] ), strPhysiCasing.str(),
587
588 if ( logicBSCTheta )
589 {
590 G4VisAttributes* rightVisAtt = new G4VisAttributes( G4Colour( 1.0, 0., 0. ) );
591 rightVisAtt->SetVisibility( true );
592 logicBSCTheta->SetVisAttributes( rightVisAtt );
593 logicBSCTheta->SetVisAttributes( G4VisAttributes::Invisible );
594 }
595
596 ostringstream strRear;
597 strRear << "physicalRearBox_1_" << i - 1;
598
599 physiRear =
600 new G4PVPlacement( rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1],
601 G4ThreeVector( ( *besEMCGeometry ).rearBoxPosX[id],
602 ( *besEMCGeometry ).rearBoxPosY[id],
603 ( *besEMCGeometry ).rearBoxPosZ[id] ),
605
606 ostringstream strGirder;
607 strGirder << "solidOpenningCutGirder_1_" << i - 1;
608 solidOCGirder =
609 new G4Cons( strGirder.str(), ( *besEMCGeometry ).OCGirderRmin1[id],
610 ( *besEMCGeometry ).BSCPhiRmax, ( *besEMCGeometry ).OCGirderRmin2[id],
611 ( *besEMCGeometry ).BSCPhiRmax, ( *besEMCGeometry ).OCGirderDz[id] / 2,
612 360. * deg - ( *besEMCGeometry ).OCGirderAngle / 2,
613 ( *besEMCGeometry ).OCGirderAngle / 2 - da );
614
615 ostringstream strVGirder;
616 strVGirder << "logicalOpenningCutGirder_1_" << i - 1;
617 logicOCGirder = new G4LogicalVolume( solidOCGirder, stainlessSteel, strVGirder.str() );
618 logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
619
620 ostringstream strPGirder;
621 strPGirder << "physicalOpenningCutGirder_1_" << i - 1;
622 physiOCGirder = new G4PVPlacement(
623 0, G4ThreeVector( 0, 0, ( *besEMCGeometry ).OCGirderPosZ[id] ), logicOCGirder,
625
626 if ( id < ( *besEMCGeometry ).BSCNbTheta - 1 )
627 {
628 G4double zLength = ( *besEMCGeometry ).OCGirderPosZ[id + 1] -
629 ( *besEMCGeometry ).OCGirderPosZ[id] -
630 ( *besEMCGeometry ).OCGirderDz[id + 1] / 2 -
631 ( *besEMCGeometry ).OCGirderDz[id] / 2;
632 G4double zPositionOCGirder = ( *besEMCGeometry ).OCGirderPosZ[id + 1] -
633 ( *besEMCGeometry ).OCGirderDz[id + 1] / 2 -
634 zLength / 2;
635
636 ostringstream strGirder2;
637 strGirder2 << "solidOpenningCutGirder_2_" << i - 1;
638 solidOCGirder = new G4Cons( strGirder2.str(), ( *besEMCGeometry ).OCGirderRmin2[id],
639 ( *besEMCGeometry ).BSCPhiRmax,
640 ( *besEMCGeometry ).OCGirderRmin1[id + 1],
641 ( *besEMCGeometry ).BSCPhiRmax, zLength / 2,
642 360. * deg - ( *besEMCGeometry ).OCGirderAngle / 2,
643 ( *besEMCGeometry ).OCGirderAngle / 2 - da );
644
645 ostringstream strVGirder2;
646 strVGirder2 << "logicalOpenningCutGirder_2_" << i - 1;
647 logicOCGirder =
648 new G4LogicalVolume( solidOCGirder, stainlessSteel, strVGirder2.str() );
649 logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
650
651 ostringstream strPGirder2;
652 strPGirder2 << "physicalOpenningCutGirder_2_" << i - 1;
653 physiOCGirder =
654 new G4PVPlacement( 0, G4ThreeVector( 0, 0, zPositionOCGirder ), logicOCGirder,
656 }
657
658 ostringstream strBSCCable;
659 strBSCCable << "solidBSCCable_1_" << i - 1;
660 solidCable =
661 new G4Tubs( strBSCCable.str(), 0, ( *besEMCGeometry ).cableDr,
662 ( *besEMCGeometry ).cableLength[id] / 2, 0. * deg, 360. * deg );
663
664 ostringstream strVBSCCable;
665 strVBSCCable << "logicalBSCCable_1_" << i - 1;
666 logicCable = new G4LogicalVolume( solidCable, cable, strVBSCCable.str() );
667
668 ostringstream strPBSCCable;
669 strPBSCCable << "physicalBSCCable_1_" << i - 1;
670 physiCable = new G4PVPlacement( 0,
671 G4ThreeVector( ( *besEMCGeometry ).cablePosX[id],
672 ( *besEMCGeometry ).cablePosY[id],
673 ( *besEMCGeometry ).cablePosZ[id] ),
674 logicCable, strPBSCCable.str(), logicBSCPhi, false, 0,
675 CHECKLV3 );
676 logicCable->SetVisAttributes( G4VisAttributes::Invisible );
677 }
678 else
679 {
680 id = ( *besEMCGeometry ).BSCNbTheta - i;
681 solidBSCTheta =
682 new G4Trap( strSolidCasing.str(), zHalfLength[id], thetaAxis[id], phiAxis[id],
683 yHalfLength1[id], xHalfLength1[id], xHalfLength2[id], tanAlpha1[id],
684 yHalfLength2[id], xHalfLength3[id], xHalfLength4[id], tanAlpha2[id] );
685
686 // G4cout<<"in EmcConstr2: "<<strSolidCasing.str()<<"
687 // x1="<<xHalfLength1[id]<<" y1="<<yHalfLength1[id]<<" theta="<<thetaAxis[id]
688 // <<" phi="<<phiAxis[id]<<" a1="<<tanAlpha1[id]<<G4endl;
689
690 logicBSCTheta =
691 new G4LogicalVolume( solidBSCTheta, fCasingMaterial, strVolumeCasing.str() );
692
693 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1] = new G4RotationMatrix();
694 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1]->rotateZ( -90 * deg );
695 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1]->rotateX( -180 * deg +
696 thetaPosition[id] );
697 physiBSCTheta = new G4PVPlacement(
698 rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1],
699 G4ThreeVector( xPosition[id], yPosition[id], -zPosition[id] ),
701 if ( logicBSCTheta )
702 {
703 G4VisAttributes* rightVisAtt = new G4VisAttributes( G4Colour( 1.0, 0., 0. ) );
704 rightVisAtt->SetVisibility( true );
705 logicBSCTheta->SetVisAttributes( rightVisAtt );
706 logicBSCTheta->SetVisAttributes( G4VisAttributes::Invisible );
707 }
708
709 ostringstream strRear;
710 strRear << "physicalRearBox_2_" << i - 1;
711
712 physiRear =
713 new G4PVPlacement( rotateMatrix[( *besEMCGeometry ).BSCNbPhi + i - 1],
714 G4ThreeVector( ( *besEMCGeometry ).rearBoxPosX[id],
715 ( *besEMCGeometry ).rearBoxPosY[id],
716 -( *besEMCGeometry ).rearBoxPosZ[id] ),
718
719 ostringstream strGirder;
720 strGirder << "solidOpenningCutGirder_3_" << i - 1;
721 solidOCGirder =
722 new G4Cons( strGirder.str(), ( *besEMCGeometry ).OCGirderRmin2[id],
723 ( *besEMCGeometry ).BSCPhiRmax, ( *besEMCGeometry ).OCGirderRmin1[id],
724 ( *besEMCGeometry ).BSCPhiRmax, ( *besEMCGeometry ).OCGirderDz[id] / 2,
725 360. * deg - ( *besEMCGeometry ).OCGirderAngle / 2,
726 ( *besEMCGeometry ).OCGirderAngle / 2 - da );
727
728 ostringstream strVGirder;
729 strVGirder << "logicalOpenningCutGirder_3_" << i - 1;
730 logicOCGirder = new G4LogicalVolume( solidOCGirder, stainlessSteel, strVGirder.str() );
731 logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
732
733 ostringstream strPGirder;
734 strPGirder << "physicalOpenningCutGirder_3_" << i - 1;
735 physiOCGirder = new G4PVPlacement(
736 0, G4ThreeVector( 0, 0, -( *besEMCGeometry ).OCGirderPosZ[id] ), logicOCGirder,
738
739 if ( id < ( *besEMCGeometry ).BSCNbTheta - 1 )
740 {
741 G4double zLength = ( *besEMCGeometry ).OCGirderPosZ[id + 1] -
742 ( *besEMCGeometry ).OCGirderPosZ[id] -
743 ( *besEMCGeometry ).OCGirderDz[id + 1] / 2 -
744 ( *besEMCGeometry ).OCGirderDz[id] / 2;
745 G4double zPositionOCGirder = ( *besEMCGeometry ).OCGirderPosZ[id + 1] -
746 ( *besEMCGeometry ).OCGirderDz[id + 1] / 2 -
747 zLength / 2;
748
749 ostringstream strGirder2;
750 strGirder2 << "solidOpenningCutGirder_4_" << i - 1;
751 solidOCGirder = new G4Cons(
752 strGirder2.str(), ( *besEMCGeometry ).OCGirderRmin1[id + 1],
753 ( *besEMCGeometry ).BSCPhiRmax, ( *besEMCGeometry ).OCGirderRmin2[id],
754 ( *besEMCGeometry ).BSCPhiRmax, zLength / 2,
755 360. * deg - ( *besEMCGeometry ).OCGirderAngle / 2,
756 ( *besEMCGeometry ).OCGirderAngle / 2 - da );
757
758 ostringstream strVGirder2;
759 strVGirder2 << "logicalOpenningCutGirder_4_" << i - 1;
760 logicOCGirder =
761 new G4LogicalVolume( solidOCGirder, stainlessSteel, strVGirder2.str() );
762 logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
763
764 ostringstream strPGirder2;
765 strPGirder2 << "physicalOpenningCutGirder_4_" << i - 1;
766 physiOCGirder =
767 new G4PVPlacement( 0, G4ThreeVector( 0, 0, -zPositionOCGirder ), logicOCGirder,
769 }
770
771 ostringstream strBSCCable;
772 strBSCCable << "solidBSCCable_2_" << i - 1;
773 solidCable =
774 new G4Tubs( strBSCCable.str(), 0, ( *besEMCGeometry ).cableDr,
775 ( *besEMCGeometry ).cableLength[id] / 2, 0. * deg, 360. * deg );
776
777 ostringstream strVBSCCable;
778 strVBSCCable << "logicalBSCCable_2_" << i - 1;
779 logicCable = new G4LogicalVolume( solidCable, cable, strVBSCCable.str() );
780
781 ostringstream strPBSCCable;
782 strPBSCCable << "physicalBSCCable_2_" << i - 1;
783 physiCable = new G4PVPlacement( 0,
784 G4ThreeVector( ( *besEMCGeometry ).cablePosX[id],
785 ( *besEMCGeometry ).cablePosY[id],
786 -( *besEMCGeometry ).cablePosZ[id] ),
787 logicCable, strPBSCCable.str(), logicBSCPhi, false, 0,
788 CHECKLV3 );
789 logicCable->SetVisAttributes( G4VisAttributes::Invisible );
790 }
791
792 ostringstream strCrystal;
793 strCrystal << "physicalCrystal" << i - 1;
794 physiBSCCrystal =
795 new G4PVParameterised( strCrystal.str(), logicBSCCrystal, physiBSCTheta, kZAxis,
796 1, // for this method,it must be 1.
797 crystalParam );
798 ( *besEMCGeometry ).physiBSCCrystal[i] = physiBSCCrystal;
799 // G4cout << (*besEMCGeometry).physiBSCCrystal[i] << G4endl;
800 physiBSCCrystal->SetCopyNo( i );
802
803 if ( verboseLevel > 4 )
804 G4cout
805 << "BesEmcConstruction*****************************" << G4endl
806 << "point of crystal =" << physiBSCCrystal
807 << G4endl
808 // << "point of mother =" <<physiBSCCrystal->GetMotherPhysical() << G4endl
809 << "point of excepted=" << physiBSCTheta << G4endl;
810 // G4Exception("BesEMCConstruction::Construct() starting............");
811 }
812 //
813 // always return the physical World
814 //
816 // return physiBSC;
817
818 ConstructSPFrame( logicBSCWorld, besEMCGeometry );
819 ConstructEndGeometry( logicEMC );
820 }
821
822 // Set vis attributes and sensitive detector
823 SetVisAndSD();
824
825 // list geo tree
826 if ( logicEMC && physiEMC && verboseLevel > 4 )
827 {
828 G4cout << "logicEmc " << logicEMC << " physiEmc " << physiEMC << G4endl;
829 G4cout << "list geo tree" << G4endl;
830
831 int NdaughterofEMC = logicEMC->GetNoDaughters();
832
833 for ( int i = 0; i < NdaughterofEMC; i++ )
834 {
835 G4LogicalVolume* daughterofEmc = logicEMC->GetDaughter( i )->GetLogicalVolume();
836 G4cout << i << "/" << NdaughterofEMC << " name: " << daughterofEmc->GetName() << " "
837 << daughterofEmc << " shape: " << daughterofEmc->GetSolid()->GetName() << G4endl;
838 int NdaughterofEmc_2 = daughterofEmc->GetNoDaughters();
839 for ( int j = 0; j < NdaughterofEmc_2; j++ )
840 {
841 G4LogicalVolume* daughterofEmc_2 = daughterofEmc->GetDaughter( j )->GetLogicalVolume();
842 G4cout << " --> " << j << "/" << NdaughterofEmc_2
843 << " name: " << daughterofEmc_2->GetName() << " " << daughterofEmc_2
844 << " shape: " << daughterofEmc_2->GetSolid()->GetName() << G4endl;
845 int NdaughterofEmc_3 = daughterofEmc_2->GetNoDaughters();
846 for ( int k = 0; k < NdaughterofEmc_3; k++ )
847 {
848 G4LogicalVolume* daughterofEmc_3 =
849 daughterofEmc_2->GetDaughter( k )->GetLogicalVolume();
850 G4cout << " --> " << k << "/" << NdaughterofEmc_3
851 << " name: " << daughterofEmc_3->GetName() << " " << daughterofEmc_3
852 << " shape: " << daughterofEmc_3->GetSolid()->GetName() << G4endl;
853 int NdaughterofEmc_4 = daughterofEmc_3->GetNoDaughters();
854 for ( int m = 0; m < NdaughterofEmc_4; m++ )
855 {
856 G4LogicalVolume* daughterofEmc_4 =
857 daughterofEmc_3->GetDaughter( m )->GetLogicalVolume();
858 G4cout << " --> " << m << "/" << NdaughterofEmc_4
859 << " name: " << daughterofEmc_4->GetName() << " " << daughterofEmc_4
860 << " shape: " << daughterofEmc_4->GetSolid()->GetName() << G4endl;
861 if ( daughterofEmc_3->GetSolid()->GetName().contains( "solidBSCCasing" ) )
862 {
863 G4Trap* Crystal = (G4Trap*)daughterofEmc_3->GetSolid();
864 double hz = Crystal->GetZHalfLength();
865 double hx1 = Crystal->GetXHalfLength1();
866 double hx2 = Crystal->GetXHalfLength2();
867 double hx3 = Crystal->GetXHalfLength3();
868 double hx4 = Crystal->GetXHalfLength4();
869 double hy1 = Crystal->GetYHalfLength1();
870 double hy2 = Crystal->GetYHalfLength2();
871 double tanalpha1 = Crystal->GetTanAlpha1();
872 double tanalpha2 = Crystal->GetTanAlpha2();
873 G4cout << " --> " << hx1 << " " << hx2 << " " << hx3 << " "
874 << hx4 << " " << hy1 << " " << hy2 << " " << hz << " " << tanalpha1 << " "
875 << tanalpha2 << G4endl;
876
877 } // if(SolidCrystal)
878 } // 4
879 } // 3
880 } // 2
881 } // 1
882 }
883}
double tan(const BesAngle a)
double sin(const BesAngle a)
double cos(const BesAngle a)
void ConstructSPFrame(G4LogicalVolume *, BesEmcGeometry *)
Definition BesEmcConstruction.cc:1195
void ConstructEndGeometry(G4LogicalVolume *)
Definition BesEmcConstruction.cc:885
void PrintEMCParameters()
Definition BesEmcConstruction.cc:1632
static EmcG4Geo * Instance()
Get a pointer to the single instance of EmcG4Geo.
Definition EmcG4Geo.cxx:43
static G4int GetEmc()
G4LogicalVolume * GetTopVolume()
Get the top(world) volume;.
◆ Construct() [2/3]
|
virtual |
Implements BesSubdetector.
◆ Construct() [3/3]
|
virtual |
Implements BesSubdetector.
◆ ConstructEndGeometry() [1/3]
| void BesEmcConstruction::ConstructEndGeometry | ( | G4LogicalVolume * | logEMC | ) |
Definition at line 885 of file BesEmcConstruction.cc.
885 { // logicEMC is not
886 // necessary, since
887 // it is a member of
888 // this class
889 if ( logEMC != logicEMC )
890 G4cout << "BesEmcConstruction::ConstructEndGeometry parameter transmit error!" << G4endl;
891 // has gotten in DefineMaterials()
892 // G4Material* fCrystalMaterial = G4Material::GetMaterial("Cesiumiodide");
893 G4VisAttributes* crystalVisAtt = new G4VisAttributes( G4Colour( 0.5, 0, 1.0 ) );
894 crystalVisAtt->SetVisibility( false );
895 G4VisAttributes* endPhiVisAtt = new G4VisAttributes( G4Colour( 0, 1.0, 0 ) );
896 endPhiVisAtt->SetVisibility( false );
897 const G4double zoomConst = 0.995;
898 const G4double da = 0.001 * deg;
899
900 // world volume of endcap
901 // east end
902 solidEnd = new G4Cons( "solidEndWorld", ( *emcEnd ).WorldRmin1, ( *emcEnd ).WorldRmax1,
903 ( *emcEnd ).WorldRmin2, ( *emcEnd ).WorldRmax2,
904 ( *emcEnd ).WorldDz / 2, 0. * deg, 360. * deg );
905 logicEnd = new G4LogicalVolume( solidEnd, G4Material::GetMaterial( "Aluminium" ),
906 "logicalEndWorld", 0, 0, 0 );
907 physiEnd = new G4PVPlacement( 0, // no rotation
908 G4ThreeVector( 0, 0, ( *emcEnd ).WorldZPosition ),
909 logicEnd, // its logical volume
910 "physicalEndWorld0", // its name
911 logicEMC, // its mother volume
912 false, // no boolean operations
914 if ( logicEnd ) logicEnd->SetVisAttributes( G4VisAttributes::Invisible );
915
916 // west end
917 G4RotationMatrix* rotateEnd = new G4RotationMatrix();
918 rotateEnd->rotateY( 180. * deg );
919 physiEnd = new G4PVPlacement( rotateEnd, G4ThreeVector( 0, 0, -( *emcEnd ).WorldZPosition ),
921
922 ////////////////////////////////////////////////////////////////////////
923 // emc endcap sectors (east) //
924 //////////////////////////////////////////////////////////////////////////
925 // 20mm gap //
926 // || //
927 // \ 7 || 6 / //
928 // - 8 \ || / 5 - //
929 // - \ || / - //
930 // _ 9 - \ || / - 4 _ //
931 // - _ - \ || / - _ - //
932 // - _ - \||/ - _ - //
933 // 10 - -||- - 3 //
934 // ----------------||---------------- //
935 // 11 - -||- - 2 //
936 // _ - - /||\ - - _ //
937 // _ - - / || \ - - _ //
938 // - 12 - / || \ - 1 - //
939 // - / || \ - //
940 // - 13 / || \ 0 - //
941 // / 14 || 15 \ //
942 // || //
943 ////////////////////////////////////////////////////////////////////////
944
945 // 1/16 of endcap world,which has some symmetry
946 // sector 0-6,8-14
947 solidEndPhi = new G4Cons( "solidEndPhi0", ( *emcEnd ).SectorRmin1, ( *emcEnd ).SectorRmax1,
948 ( *emcEnd ).SectorRmin2, ( *emcEnd ).SectorRmax2,
949 ( *emcEnd ).SectorDz / 2, 0. * deg, 22.5 * deg - da );
950 logicEndPhi = new G4LogicalVolume( solidEndPhi, G4Material::GetMaterial( "Air" ),
951 "logicalEndPhi0", 0, 0, 0 );
952 for ( G4int i = 0; i < 14; i++ )
953 {
954 if ( ( i != 6 ) && ( i != 7 ) )
955 {
956 G4RotationMatrix* rotatePhi = new G4RotationMatrix();
957 rotatePhi->rotateZ( -i * 22.5 * deg + 67.5 * deg );
958 ostringstream strEndPhi;
959 strEndPhi << "physicalEndPhi" << i;
960 physiEndPhi =
961 new G4PVPlacement( rotatePhi, // 0,logicEndPhi,strEndPhi.str(),logicEnd,false,i);
962 G4ThreeVector( 0, 0, ( *emcEnd ).SectorZPosition ), logicEndPhi,
964 }
965 }
966 if ( logicEndPhi ) logicEndPhi->SetVisAttributes( endPhiVisAtt );
967
968 for ( G4int i = 0; i < 35; i++ )
969 {
970 ostringstream strEndCasing;
971 strEndCasing << "solidEndCasing_0_" << i;
972
973 //-************tranform to new coodinate! liangyt 2007.5.7 *******
974 G4ThreeVector newfPnt[8];
975 G4ThreeVector center( 0.0, 0.0, 0.0 );
976 G4ThreeVector rotAngle( 0.0, 0.0, 0.0 );
977
978 TransformToArb8( ( *emcEnd ).fPnt[i], newfPnt, center, rotAngle );
979
980 emcEnd->Zoom( newfPnt, zoomConst ); // change emcEnd.fPnt[i] to newfPnt
981
982 G4RotationMatrix* rotatePhiIrregBox = new G4RotationMatrix();
983 rotatePhiIrregBox->rotateX( rotAngle.x() );
984 rotatePhiIrregBox->rotateY( rotAngle.y() );
985 rotatePhiIrregBox->rotateZ( rotAngle.z() );
986 //-*******************************************************************
987
988 solidEndCasing = new G4IrregBox( strEndCasing.str(), ( *emcEnd ).zoomPoint ); // liangyt
989
990 ostringstream strVEndCasing;
991 strVEndCasing << "logicalEndCasing_0_" << i;
992 logicEndCasing =
993 new G4LogicalVolume( solidEndCasing, fCasingMaterial, strVEndCasing.str() );
994
995 ostringstream strPEndCasing;
996 strPEndCasing << "physicalEndCasing_0_" << i;
997 physiEndCasing =
998 new G4PVPlacement( rotatePhiIrregBox, center, logicEndCasing, strPEndCasing.str(),
1000
1001 ostringstream strEndCrystal;
1002 strEndCrystal << "solidEndCrystal_0_" << i;
1003
1004 emcEnd->ModifyForCasing( ( *emcEnd ).zoomPoint, i );
1005 solidEndCrystal = new G4IrregBox( strEndCrystal.str(), ( *emcEnd ).cryPoint );
1006
1007 ostringstream strVEndCrystal;
1008 strVEndCrystal << "logicalEndCrystal_0_" << i;
1009 logicEndCrystal =
1010 new G4LogicalVolume( solidEndCrystal, fCrystalMaterial, strVEndCrystal.str() );
1011
1012 ostringstream strPEndCrystal;
1013 strPEndCrystal << "physicalEndCrystal_0_" << i;
1014 physiEndCrystal =
1015 new G4PVPlacement( 0, G4ThreeVector(), logicEndCrystal, strPEndCrystal.str(),
1017
1018 logicEndCasing->SetVisAttributes( G4VisAttributes::Invisible );
1019 logicEndCrystal->SetVisAttributes( crystalVisAtt );
1020 logicEndCrystal->SetSensitiveDetector( besEMCSD );
1021 }
1022
1023 // the top area which has 20 mm gap
1024 // sector 6,14
1025 G4double rmin2 =
1026 ( *emcEnd ).WorldRmin1 + ( ( *emcEnd ).WorldRmin2 - ( *emcEnd ).WorldRmin1 ) *
1027 ( *emcEnd ).SectorDz / ( *emcEnd ).WorldDz;
1028 G4double rmax2 =
1029 ( *emcEnd ).WorldRmax1 + ( ( *emcEnd ).WorldRmax2 - ( *emcEnd ).WorldRmax1 ) *
1030 ( *emcEnd ).SectorDz / ( *emcEnd ).WorldDz;
1031 solidEndPhi =
1032 new G4Cons( "solidEndPhi1", ( *emcEnd ).WorldRmin1, ( *emcEnd ).WorldRmax1, rmin2, rmax2,
1033 ( *emcEnd ).SectorDz / 2, 67.5 * deg, 22.5 * deg - da );
1034 logicEndPhi = new G4LogicalVolume( solidEndPhi, G4Material::GetMaterial( "Air" ),
1035 "logicalEndPhi1", 0, 0, 0 );
1036 for ( G4int i = 0; i < 2; i++ )
1037 {
1038 G4RotationMatrix* rotatePhi = new G4RotationMatrix();
1039 rotatePhi->rotateZ( -i * 180. * deg );
1040 ostringstream strEndPhi;
1041 strEndPhi << "physicalEndPhi" << i * 8 + 6;
1042 physiEndPhi = new G4PVPlacement(
1043 rotatePhi, G4ThreeVector( 0, 0, ( *emcEnd ).SectorZPosition ), logicEndPhi,
1045 }
1046 if ( logicEndPhi ) logicEndPhi->SetVisAttributes( endPhiVisAtt );
1047
1048 for ( G4int i = 0; i < 35; i++ )
1049 {
1050 ostringstream strEndCasing;
1051 strEndCasing << "solidEndCasing_1_" << i;
1052
1053 //-************tranform to new coodinate! liangyt 2007.5.7 *******
1054 G4ThreeVector newfPnt[8];
1055 G4ThreeVector center( 0.0, 0.0, 0.0 );
1056 G4ThreeVector rotAngle( 0.0, 0.0, 0.0 );
1057
1058 TransformToArb8( ( *emcEnd ).fPnt1[i], newfPnt, center, rotAngle );
1059
1060 emcEnd->Zoom( newfPnt, zoomConst ); // change emcEnd.fPnt[i] to newfPnt
1061
1062 G4RotationMatrix* rotatePhiIrregBox = new G4RotationMatrix();
1063 rotatePhiIrregBox->rotateX( rotAngle.x() );
1064 rotatePhiIrregBox->rotateY( rotAngle.y() );
1065 rotatePhiIrregBox->rotateZ( rotAngle.z() );
1066 //-*******************************************************************
1067
1068 solidEndCasing = new G4IrregBox( strEndCasing.str(), ( *emcEnd ).zoomPoint );
1069
1070 ostringstream strVEndCasing;
1071 strVEndCasing << "logicalEndCasing_1_" << i;
1072 logicEndCasing =
1073 new G4LogicalVolume( solidEndCasing, fCasingMaterial, strVEndCasing.str() );
1074
1075 ostringstream strPEndCasing;
1076 strPEndCasing << "physicalEndCasing_1_" << i;
1077 physiEndCasing =
1078 new G4PVPlacement( rotatePhiIrregBox, center, logicEndCasing, strPEndCasing.str(),
1080
1081 ostringstream strEndCrystal;
1082 strEndCrystal << "solidEndCrystal_1_" << i;
1083
1084 emcEnd->ModifyForCasing( ( *emcEnd ).zoomPoint, i );
1085 solidEndCrystal = new G4IrregBox( strEndCrystal.str(), ( *emcEnd ).cryPoint );
1086
1087 ostringstream strVEndCrystal;
1088 strVEndCrystal << "logicalEndCrystal_1_" << i;
1089 logicEndCrystal =
1090 new G4LogicalVolume( solidEndCrystal, fCrystalMaterial, strVEndCrystal.str() );
1091
1092 ostringstream strPEndCrystal;
1093 strPEndCrystal << "physicalEndCrystal_1_" << i;
1094 physiEndCrystal =
1095 new G4PVPlacement( 0, G4ThreeVector(), logicEndCrystal, strPEndCrystal.str(),
1097
1098 logicEndCasing->SetVisAttributes( G4VisAttributes::Invisible );
1099 logicEndCrystal->SetVisAttributes( crystalVisAtt );
1100 logicEndCrystal->SetSensitiveDetector( besEMCSD );
1101 }
1102
1103 ( *emcEnd ).ReflectX();
1104
1105 // sector 7,15
1106 for ( G4int i = 0; i < 35; i++ )
1107 for ( G4int j = 0; j < 8; j++ ) ( *emcEnd ).fPnt1[i][j].rotateZ( -90. * deg );
1108
1109 solidEndPhi = new G4Cons( "solidEndPhi2", ( *emcEnd ).WorldRmin1, ( *emcEnd ).WorldRmax1,
1110 rmin2, rmax2, ( *emcEnd ).SectorDz / 2, 0 * deg, 22.5 * deg - da );
1111 logicEndPhi = new G4LogicalVolume( solidEndPhi, G4Material::GetMaterial( "Air" ),
1112 "logicalEndPhi2", 0, 0, 0 );
1113 for ( G4int i = 0; i < 2; i++ )
1114 {
1115 G4RotationMatrix* rotatePhi = new G4RotationMatrix();
1116 rotatePhi->rotateZ( -i * 180. * deg - 90. * deg );
1117 ostringstream strEndPhi;
1118 strEndPhi << "physicalEndPhi" << i * 8 + 7;
1119 physiEndPhi = new G4PVPlacement(
1120 rotatePhi, G4ThreeVector( 0, 0, ( *emcEnd ).SectorZPosition ), logicEndPhi,
1122 }
1123 if ( logicEndPhi ) logicEndPhi->SetVisAttributes( endPhiVisAtt );
1124
1125 for ( G4int i = 0; i < 35; i++ )
1126 {
1127 ostringstream strEndCasing;
1128 strEndCasing << "solidEndCasing_2_" << i;
1129
1130 //-************tranform to new coodinate! liangyt 2007.5.7 *******
1131 G4ThreeVector newfPnt[8];
1132 G4ThreeVector center( 0.0, 0.0, 0.0 );
1133 G4ThreeVector rotAngle( 0.0, 0.0, 0.0 );
1134
1135 TransformToArb8( ( *emcEnd ).fPnt1[i], newfPnt, center, rotAngle );
1136
1137 emcEnd->Zoom( newfPnt, zoomConst ); // change emcEnd.fPnt[i] to newfPnt
1138
1139 G4RotationMatrix* rotatePhiIrregBox = new G4RotationMatrix();
1140 rotatePhiIrregBox->rotateX( rotAngle.x() );
1141 rotatePhiIrregBox->rotateY( rotAngle.y() );
1142 rotatePhiIrregBox->rotateZ( rotAngle.z() );
1143 //-*******************************************************************
1144
1145 solidEndCasing = new G4IrregBox( strEndCasing.str(), ( *emcEnd ).zoomPoint );
1146
1147 ostringstream strVEndCasing;
1148 strVEndCasing << "logicalEndCasing_2_" << i;
1149 logicEndCasing =
1150 new G4LogicalVolume( solidEndCasing, fCasingMaterial, strVEndCasing.str() );
1151
1152 ostringstream strPEndCasing;
1153 strPEndCasing << "physicalEndCasing_2_" << i;
1154 physiEndCasing =
1155 new G4PVPlacement( rotatePhiIrregBox, center, logicEndCasing, strPEndCasing.str(),
1157
1158 ostringstream strEndCrystal;
1159 strEndCrystal << "solidEndCrystal_2_" << i;
1160
1161 emcEnd->ModifyForCasing( ( *emcEnd ).zoomPoint, i );
1162 solidEndCrystal = new G4IrregBox( strEndCrystal.str(), ( *emcEnd ).cryPoint );
1163
1164 ostringstream strVEndCrystal;
1165 strVEndCrystal << "logicalEndCrystal_2_" << i;
1166 logicEndCrystal =
1167 new G4LogicalVolume( solidEndCrystal, fCrystalMaterial, strVEndCrystal.str() );
1168
1169 ostringstream strPEndCrystal;
1170 strPEndCrystal << "physicalEndCrystal_2_" << i;
1171 physiEndCrystal =
1172 new G4PVPlacement( 0, G4ThreeVector(), logicEndCrystal, strPEndCrystal.str(),
1174
1175 logicEndCasing->SetVisAttributes( G4VisAttributes::Invisible );
1176 logicEndCrystal->SetVisAttributes( crystalVisAtt );
1177 logicEndCrystal->SetSensitiveDetector( besEMCSD );
1178 }
1179}
void TransformToArb8(const G4ThreeVector fPnt[8], G4ThreeVector newfPnt[8], G4ThreeVector ¢er, G4ThreeVector &rotAngle)
Definition BesEmcConstruction.cc:1759
Referenced by Construct().
◆ ConstructEndGeometry() [2/3]
| void BesEmcConstruction::ConstructEndGeometry | ( | G4LogicalVolume * | ) |
◆ ConstructEndGeometry() [3/3]
| void BesEmcConstruction::ConstructEndGeometry | ( | G4LogicalVolume * | ) |
◆ ConstructSPFrame() [1/3]
| void BesEmcConstruction::ConstructSPFrame | ( | G4LogicalVolume * | logMother, |
| BesEmcGeometry * | emcGeometry ) |
Definition at line 1195 of file BesEmcConstruction.cc.
1196 {
1197 G4double rmax = ( *besEMCGeometry ).BSCRmax + 2. * mm; // radius from 942mm to 940mm
1198 solidSupportBar = new G4Tubs(
1199 "solidSupportBar0", rmax + ( *besEMCGeometry ).SPBarThickness1,
1200 rmax + ( *besEMCGeometry ).SPBarThickness + ( *besEMCGeometry ).SPBarThickness1,
1201 ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 +
1202 ( *besEMCGeometry ).EndRingDz,
1203 0. * deg, 360. * deg );
1204
1205 logicSupportBar =
1206 new G4LogicalVolume( solidSupportBar, stainlessSteel, "logicalSupportBar0" );
1207
1208 physiSupportBar = new G4PVPlacement( 0, G4ThreeVector(), logicSupportBar,
1210
1211 solidSupportBar1 =
1212 new G4Tubs( "solidSupportBar1", rmax, rmax + ( *besEMCGeometry ).SPBarThickness1,
1213 ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3,
1214 ( *besEMCGeometry ).BSCPhiDphi - ( *besEMCGeometry ).SPBarDphi / 2,
1215 ( *besEMCGeometry ).SPBarDphi );
1216
1217 logicSupportBar1 =
1218 new G4LogicalVolume( solidSupportBar1, stainlessSteel, "logicalSupportBar1" );
1219
1220 for ( G4int i = 0; i < ( *besEMCGeometry ).BSCNbPhi / 2; i++ )
1221 {
1222 G4RotationMatrix* rotateSPBar = new G4RotationMatrix();
1223 rotateSPBar->rotateZ( ( *besEMCGeometry ).BSCPhiDphi -
1224 i * 2 * ( *besEMCGeometry ).BSCPhiDphi );
1225 ostringstream strSupportBar1;
1226 strSupportBar1 << "physicalSupportBar1_" << i;
1227 physiSupportBar1 =
1228 new G4PVPlacement( rotateSPBar, G4ThreeVector(), logicSupportBar1,
1230 }
1231
1232 // end ring
1233 solidEndRing =
1234 new G4Tubs( "solidEndRing", ( *besEMCGeometry ).EndRingRmin,
1235 ( *besEMCGeometry ).EndRingRmin + ( *besEMCGeometry ).EndRingDr / 2,
1236 ( *besEMCGeometry ).EndRingDz / 2, 0. * deg, 360. * deg );
1237
1238 solidGear = new G4Tubs(
1239 "solidGear", ( *besEMCGeometry ).EndRingRmin + ( *besEMCGeometry ).EndRingDr / 2,
1240 ( *besEMCGeometry ).EndRingRmin + ( *besEMCGeometry ).EndRingDr,
1241 ( *besEMCGeometry ).EndRingDz / 2, 0. * deg, ( *besEMCGeometry ).BSCPhiDphi );
1242
1243 // taper ring
1244 solidTaperRing1 =
1245 new G4Tubs( "solidTaperRing1", ( *besEMCGeometry ).TaperRingRmin1,
1246 ( *besEMCGeometry ).TaperRingRmin1 + ( *besEMCGeometry ).TaperRingThickness1,
1247 ( *besEMCGeometry ).TaperRingInnerLength / 2, 0. * deg, 360. * deg );
1248
1249 solidTaperRing2 =
1250 new G4Cons( "solidTaperRing2", ( *besEMCGeometry ).TaperRingRmin1,
1251 ( *besEMCGeometry ).TaperRingRmin1 + ( *besEMCGeometry ).TaperRingThickness1,
1252 ( *besEMCGeometry ).TaperRingRmin2,
1253 ( *besEMCGeometry ).TaperRingRmin2 + ( *besEMCGeometry ).TaperRingDr,
1254 ( *besEMCGeometry ).TaperRingDz / 2, 0. * deg, 360. * deg );
1255
1256 solidTaperRing3 =
1257 new G4Cons( "solidTaperRing3", ( *besEMCGeometry ).BSCRmax2,
1258 ( *besEMCGeometry ).BSCRmax2 + ( *besEMCGeometry ).TaperRingOuterLength1,
1259 ( *besEMCGeometry ).TaperRingRmin2 + ( *besEMCGeometry ).TaperRingDr,
1260 ( *besEMCGeometry ).TaperRingRmin2 + ( *besEMCGeometry ).TaperRingDr +
1261 ( *besEMCGeometry ).TaperRingOuterLength,
1262 ( *besEMCGeometry ).TaperRingThickness3 / 2, 0. * deg, 360. * deg );
1263
1264 logicEndRing = new G4LogicalVolume( solidEndRing, stainlessSteel, "logicalEndRing" );
1265 logicGear = new G4LogicalVolume( solidGear, stainlessSteel, "logicalGear" );
1266 logicTaperRing1 =
1267 new G4LogicalVolume( solidTaperRing1, stainlessSteel, "logicalTaperRing1" );
1268 logicTaperRing2 =
1269 new G4LogicalVolume( solidTaperRing2, stainlessSteel, "logicalTaperRing2" );
1270 logicTaperRing3 =
1271 new G4LogicalVolume( solidTaperRing3, stainlessSteel, "logicalTaperRing3" );
1272
1273 for ( G4int i = 0; i < 2; i++ )
1274 {
1275 G4RotationMatrix* rotateSPRing = new G4RotationMatrix();
1276 G4double zEndRing, z1, z2, z3;
1277 if ( i == 0 )
1278 {
1279 zEndRing = ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 +
1280 ( *besEMCGeometry ).EndRingDz / 2;
1281 z1 = ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 -
1282 ( *besEMCGeometry ).TaperRingDz - ( *besEMCGeometry ).TaperRingInnerLength / 2;
1283 z2 = ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 -
1284 ( *besEMCGeometry ).TaperRingDz / 2;
1285 z3 = ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 / 2;
1286 }
1287 else
1288 {
1289 rotateSPRing->rotateY( 180. * deg );
1290 zEndRing = -( ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 +
1291 ( *besEMCGeometry ).EndRingDz / 2 );
1292 z1 = -( ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 -
1293 ( *besEMCGeometry ).TaperRingDz - ( *besEMCGeometry ).TaperRingInnerLength / 2 );
1294 z2 = -( ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 -
1295 ( *besEMCGeometry ).TaperRingDz / 2 );
1296 z3 = -( ( *besEMCGeometry ).BSCDz + ( *besEMCGeometry ).TaperRingThickness3 / 2 );
1297 }
1298
1299 ostringstream strEndRing;
1300 strEndRing << "physicalEndRing_" << i;
1301 physiEndRing =
1302 new G4PVPlacement( rotateSPRing, G4ThreeVector( 0, 0, zEndRing ), logicEndRing,
1304
1305 for ( G4int j = 0; j < ( *besEMCGeometry ).BSCNbPhi / 2; j++ )
1306 {
1307 G4RotationMatrix* rotateGear = new G4RotationMatrix();
1308 rotateGear->rotateZ( ( *besEMCGeometry ).BSCPhiDphi / 2 -
1309 j * 2 * ( *besEMCGeometry ).BSCPhiDphi );
1310
1311 ostringstream strGear;
1312 strGear << "physicalGear_" << i << "_" << j;
1313 physiGear = new G4PVPlacement( rotateGear, G4ThreeVector( 0, 0, zEndRing ), logicGear,
1315 }
1316
1317 ostringstream strTaperRing1;
1318 strTaperRing1 << "physicalTaperRing1_" << i;
1319 physiTaperRing1 =
1320 new G4PVPlacement( rotateSPRing, G4ThreeVector( 0, 0, z1 ), logicTaperRing1,
1322
1323 ostringstream strTaperRing2;
1324 strTaperRing2 << "physicalTaperRing2_" << i;
1325 physiTaperRing2 =
1326 new G4PVPlacement( rotateSPRing, G4ThreeVector( 0, 0, z2 ), logicTaperRing2,
1328
1329 ostringstream strTaperRing3;
1330 strTaperRing3 << "physicalTaperRing3_" << i;
1331 physiTaperRing3 =
1332 new G4PVPlacement( rotateSPRing, G4ThreeVector( 0, 0, z3 ), logicTaperRing3,
1334 }
1335}
Referenced by Construct().
◆ ConstructSPFrame() [2/3]
| void BesEmcConstruction::ConstructSPFrame | ( | G4LogicalVolume * | , |
| BesEmcGeometry * | ) |
◆ ConstructSPFrame() [3/3]
| void BesEmcConstruction::ConstructSPFrame | ( | G4LogicalVolume * | , |
| BesEmcGeometry * | ) |
◆ GetBesEmcConstruction() [1/3]
|
static |
Definition at line 68 of file BesEmcConstruction.cc.
68{ return fBesEmcConstruction; }
Referenced by BesCrystalParameterisation::ComputeIDAndSide(), and BesCrystalParameterisation::ComputeMaterial().
◆ GetBesEmcConstruction() [2/3]
|
static |
◆ GetBesEmcConstruction() [3/3]
|
static |
◆ GetBSCCrystal() [1/3]
|
inline |
Definition at line 109 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
109{ return physiBSCCrystal; };
◆ GetBSCCrystal() [2/3]
|
inline |
Definition at line 109 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
109{ return physiBSCCrystal; };
◆ GetBSCCrystal() [3/3]
|
inline |
Definition at line 109 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
109{ return physiBSCCrystal; };
◆ GetBSCPhi() [1/3]
|
inline |
Definition at line 107 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
107{ return physiBSCPhi; };
◆ GetBSCPhi() [2/3]
|
inline |
Definition at line 107 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
107{ return physiBSCPhi; };
◆ GetBSCPhi() [3/3]
|
inline |
Definition at line 107 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
107{ return physiBSCPhi; };
◆ GetBSCTheta() [1/3]
|
inline |
Definition at line 108 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
108{ return physiBSCTheta; };
◆ GetBSCTheta() [2/3]
|
inline |
Definition at line 108 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
108{ return physiBSCTheta; };
◆ GetBSCTheta() [3/3]
|
inline |
Definition at line 108 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
108{ return physiBSCTheta; };
◆ GetCasingMaterial() [1/3]
|
inline |
Definition at line 102 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
102{ return fCasingMaterial; };
Referenced by BesCrystalParameterisation::ComputeMaterial().
◆ GetCasingMaterial() [2/3]
|
inline |
Definition at line 102 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
102{ return fCasingMaterial; };
◆ GetCasingMaterial() [3/3]
|
inline |
Definition at line 102 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
102{ return fCasingMaterial; };
◆ GetCrystalMaterial() [1/3]
|
inline |
Definition at line 101 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
101{ return fCrystalMaterial; };
Referenced by BesCrystalParameterisation::ComputeMaterial().
◆ GetCrystalMaterial() [2/3]
|
inline |
Definition at line 101 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
101{ return fCrystalMaterial; };
◆ GetCrystalMaterial() [3/3]
|
inline |
Definition at line 101 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
101{ return fCrystalMaterial; };
◆ GetCrystalParam() [1/3]
|
inline |
Definition at line 110 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
110{ return crystalParam; };
◆ GetCrystalParam() [2/3]
|
inline |
Definition at line 110 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
110{ return crystalParam; };
◆ GetCrystalParam() [3/3]
|
inline |
Definition at line 110 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
110{ return crystalParam; };
◆ GetEMC() [1/3]
|
inline |
Definition at line 106 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
106{ return physiEMC; }
◆ GetEMC() [2/3]
|
inline |
Definition at line 106 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
106{ return physiEMC; }
◆ GetEMC() [3/3]
|
inline |
Definition at line 106 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
106{ return physiEMC; }
◆ GetLogicalVolume() [1/3]
| void BesEmcConstruction::GetLogicalVolume | ( | ) |
Definition at line 1396 of file BesEmcConstruction.cc.
1396 {
1397 //-------------------------------------------------------------
1398 // Barrel
1411
1412 for ( int i = 0; i < 44; i++ )
1413 {
1414 std::ostringstream osnameBSCCasing;
1415 osnameBSCCasing << "logicalBSCCasing" << i;
1416 logicBSCTheta = FindLogicalVolume( osnameBSCCasing.str() );
1417 if ( logicBSCTheta )
1418 {
1419 G4VisAttributes* rightVisAtt = new G4VisAttributes( G4Colour( 1.0, 0.0, 0.0 ) );
1420 rightVisAtt->SetVisibility( false );
1421 logicBSCTheta->SetVisAttributes( rightVisAtt );
1422 }
1423
1424 std::ostringstream osnameBSCCable1;
1425 osnameBSCCable1 << "logicalBSCCable_1_" << i;
1426 logicCable = FindLogicalVolume( osnameBSCCable1.str() );
1427 if ( logicCable ) logicCable->SetVisAttributes( G4VisAttributes::Invisible );
1428
1429 std::ostringstream osnameBSCCable2;
1430 osnameBSCCable2 << "logicalBSCCable_2_" << i;
1431 logicCable = FindLogicalVolume( osnameBSCCable2.str() );
1432 if ( logicCable ) logicCable->SetVisAttributes( G4VisAttributes::Invisible );
1433
1434 std::ostringstream osnameOCGirder1;
1435 osnameOCGirder1 << "logicalOpenningCutGirder_1_" << i;
1436 logicOCGirder = FindLogicalVolume( osnameOCGirder1.str() );
1437 if ( logicOCGirder ) logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
1438
1439 std::ostringstream osnameOCGirder2;
1440 osnameOCGirder2 << "logicalOpenningCutGirder_2_" << i;
1441 logicOCGirder = FindLogicalVolume( osnameOCGirder2.str() );
1442 if ( logicOCGirder ) logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
1443
1444 std::ostringstream osnameOCGirder3;
1445 osnameOCGirder3 << "logicalOpenningCutGirder_3_" << i;
1446 logicOCGirder = FindLogicalVolume( osnameOCGirder3.str() );
1447 if ( logicOCGirder ) logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
1448
1449 std::ostringstream osnameOCGirder4;
1450 osnameOCGirder4 << "logicalOpenningCutGirder_4_" << i;
1451 logicOCGirder = FindLogicalVolume( osnameOCGirder4.str() );
1452 if ( logicOCGirder ) logicOCGirder->SetVisAttributes( G4VisAttributes::Invisible );
1453 }
1454
1455 //-------------------------------------------------------------
1456 // Support system
1464
1465 //-------------------------------------------------------------
1466 // Endcap
1468
1469 for ( G4int sector = 0; sector < 3; sector++ )
1470 {
1471 std::ostringstream osnameEndPhi;
1472 osnameEndPhi << "logicalEndPhi" << sector;
1473 logicEndPhi = FindLogicalVolume( osnameEndPhi.str() );
1474 if ( logicEndPhi ) { logicEndPhi->SetVisAttributes( G4VisAttributes::Invisible ); }
1475 else { G4cout << "Can't find logicEndPhi!" << G4endl; }
1476
1477 for ( G4int cryNb = 0; cryNb < 35; cryNb++ )
1478 {
1479
1480 std::ostringstream osnameEndCrystal;
1481 osnameEndCrystal << "logicalEndCrystal_" << sector << "_" << cryNb;
1482 logicEndCrystal = FindLogicalVolume( osnameEndCrystal.str() );
1483 if ( logicEndCrystal )
1484 {
1485 logicEndCrystal->SetSensitiveDetector( besEMCSD );
1486 G4VisAttributes* crystalVisAtt = new G4VisAttributes( G4Colour( 0.5, 0, 1.0 ) );
1487 crystalVisAtt->SetVisibility( false );
1488 logicEndCrystal->SetVisAttributes( crystalVisAtt );
1489 }
1490 else { G4cout << "Can't find: " << osnameEndCrystal.str() << G4endl; }
1491
1492 std::ostringstream osnameEndCasing;
1493 osnameEndCasing << "logicalEndCasing_" << sector << "_" << cryNb;
1494 logicEndCasing = FindLogicalVolume( osnameEndCasing.str() );
1495 if ( logicEndCasing ) { logicEndCasing->SetVisAttributes( G4VisAttributes::Invisible ); }
1496 else { G4cout << "Can't find: " << osnameEndCasing.str() << G4endl; }
1497 }
1498 }
1499}
G4LogicalVolume * FindLogicalVolume(const G4String &vn)
Referenced by Construct().
◆ GetLogicalVolume() [2/3]
| void BesEmcConstruction::GetLogicalVolume | ( | ) |
◆ GetLogicalVolume() [3/3]
| void BesEmcConstruction::GetLogicalVolume | ( | ) |
◆ GetMagField() [1/3]
|
inline |
Definition at line 97 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
97{ return fmagField; };
◆ GetMagField() [2/3]
|
inline |
Definition at line 97 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
97{ return fmagField; };
◆ GetMagField() [3/3]
|
inline |
Definition at line 97 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
97{ return fmagField; };
◆ GetStartIDTheta() [1/3]
|
inline |
Definition at line 99 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
99{ return startID; };
◆ GetStartIDTheta() [2/3]
|
inline |
Definition at line 99 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
99{ return startID; };
◆ GetStartIDTheta() [3/3]
|
inline |
Definition at line 99 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
99{ return startID; };
◆ GetVerboseLevel() [1/3]
|
inline |
Definition at line 96 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
96{ return verboseLevel; };
Referenced by BesCrystalParameterisation::ComputeIDAndSide().
◆ GetVerboseLevel() [2/3]
|
inline |
Definition at line 96 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
96{ return verboseLevel; };
◆ GetVerboseLevel() [3/3]
|
inline |
Definition at line 96 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
96{ return verboseLevel; };
◆ PrintEMCParameters() [1/3]
| void BesEmcConstruction::PrintEMCParameters | ( | ) |
Definition at line 1632 of file BesEmcConstruction.cc.
1632 {
1633 G4cout << "-------------------------------------------------------" << G4endl
1634 << "---> There are " << phiNbCrystals << "(max=" << ( *besEMCGeometry ).BSCNbPhi
1635 << ") crystals along phi direction and " << thetaNbCrystals
1636 << "(max=" << ( *besEMCGeometry ).BSCNbTheta << ") crystals along theta direction."
1637 << G4endl << "The crystals have sizes of " << ( *besEMCGeometry ).BSCCryLength / cm
1638 << "cm(L) and " << ( *besEMCGeometry ).BSCYFront / cm << "cm(Y) with "
1639 << fCrystalMaterial->GetName() << "." << G4endl << "The casing is layer of "
1640 << ( *besEMCGeometry ).fTyvekThickness / mm << "mm tyvek,"
1641 << ( *besEMCGeometry ).fAlThickness / mm << "mm aluminum and"
1642 << ( *besEMCGeometry ).fMylarThickness / mm << "mm mylar." << G4endl
1643 << "-------------------------------------------------------" << G4endl;
1644 G4cout << G4Material::GetMaterial( "PolyethyleneTerephthlate" ) << G4endl
1645 << G4Material::GetMaterial( "Casing" ) << G4endl
1646 << G4Material::GetMaterial( "Polyethylene" ) << G4endl
1647 << "-------------------------------------------------------" << G4endl;
1648}
Referenced by Construct(), SetCasingMaterial(), and SetCrystalMaterial().
◆ PrintEMCParameters() [2/3]
| void BesEmcConstruction::PrintEMCParameters | ( | ) |
◆ PrintEMCParameters() [3/3]
| void BesEmcConstruction::PrintEMCParameters | ( | ) |
◆ RotAngleFromNewZ() [1/3]
| Hep3Vector BesEmcConstruction::RotAngleFromNewZ | ( | Hep3Vector | newZ | ) |
Definition at line 1797 of file BesEmcConstruction.cc.
1797 {
1798 newZ.setMag( 1.0 );
1799 Hep3Vector x0( 1, 0, 0 ), y0( 0, 1, 0 ), z0( 0, 0, 1 );
1800 double dx, dy, dz = 0.0;
1801
1802 Hep3Vector a( 0.0, newZ.y(), newZ.z() );
1803 // three rotated angles, rotate dx angles(rad) by X axis,
1804 // then rotate dy by new Y axis, then dz by new Z axis;
1805 // all rotate in left hand system;
1806
1807 // a, project of final newZ on original YZ plane, 0 <= theta < pi;
1808 if ( a.mag() != 0.0 ) a.setMag( 1.0 );
1809 else cout << "newZ on X axis, a=(0,0,0)" << endl;
1810 dx = acos( a.dot( z0 ) );
1811 if ( a.dot( z0 ) == -1.0 ) dx = 0.0;
1812
1813 // b, rotate dx angle of z0 around original X axis(x0),
1815 dy = acos( b.dot( newZ ) );
1816 if ( newZ.x() > 0.0 ) dy = -dy;
1817
1818 Hep3Vector v( dx, dy, dz );
1820}
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
Definition KarLud.h:35
Referenced by TransformToArb8().
◆ RotAngleFromNewZ() [2/3]
| Hep3Vector BesEmcConstruction::RotAngleFromNewZ | ( | Hep3Vector | newZ | ) |
◆ RotAngleFromNewZ() [3/3]
| Hep3Vector BesEmcConstruction::RotAngleFromNewZ | ( | Hep3Vector | newZ | ) |
◆ SetBSCCrystalLength() [1/3]
| void BesEmcConstruction::SetBSCCrystalLength | ( | G4double | val | ) |
Definition at line 1701 of file BesEmcConstruction.cc.
1701 {
1702 ( *besEMCGeometry ).BSCCryLength = val;
1703}
◆ SetBSCCrystalLength() [2/3]
| void BesEmcConstruction::SetBSCCrystalLength | ( | G4double | ) |
◆ SetBSCCrystalLength() [3/3]
| void BesEmcConstruction::SetBSCCrystalLength | ( | G4double | ) |
◆ SetBSCNbPhi() [1/3]
| void BesEmcConstruction::SetBSCNbPhi | ( | G4int | val | ) |
Definition at line 1691 of file BesEmcConstruction.cc.
1691{ ( *besEMCGeometry ).BSCNbPhi = val; }
◆ SetBSCNbPhi() [2/3]
| void BesEmcConstruction::SetBSCNbPhi | ( | G4int | ) |
◆ SetBSCNbPhi() [3/3]
| void BesEmcConstruction::SetBSCNbPhi | ( | G4int | ) |
◆ SetBSCNbTheta() [1/3]
| void BesEmcConstruction::SetBSCNbTheta | ( | G4int | val | ) |
Definition at line 1695 of file BesEmcConstruction.cc.
1695{ ( *besEMCGeometry ).BSCNbTheta = val; }
◆ SetBSCNbTheta() [2/3]
| void BesEmcConstruction::SetBSCNbTheta | ( | G4int | ) |
◆ SetBSCNbTheta() [3/3]
| void BesEmcConstruction::SetBSCNbTheta | ( | G4int | ) |
◆ SetBSCPosition0() [1/3]
| void BesEmcConstruction::SetBSCPosition0 | ( | G4double | val | ) |
Definition at line 1717 of file BesEmcConstruction.cc.
1717 {
1718 ( *besEMCGeometry ).BSCPosition0 = val;
1719}
◆ SetBSCPosition0() [2/3]
| void BesEmcConstruction::SetBSCPosition0 | ( | G4double | ) |
◆ SetBSCPosition0() [3/3]
| void BesEmcConstruction::SetBSCPosition0 | ( | G4double | ) |
◆ SetBSCPosition1() [1/3]
| void BesEmcConstruction::SetBSCPosition1 | ( | G4double | val | ) |
Definition at line 1723 of file BesEmcConstruction.cc.
1723 {
1724 ( *besEMCGeometry ).BSCPosition1 = val;
1725}
◆ SetBSCPosition1() [2/3]
| void BesEmcConstruction::SetBSCPosition1 | ( | G4double | ) |
◆ SetBSCPosition1() [3/3]
| void BesEmcConstruction::SetBSCPosition1 | ( | G4double | ) |
◆ SetBSCRmin() [1/3]
| void BesEmcConstruction::SetBSCRmin | ( | G4double | val | ) |
Definition at line 1687 of file BesEmcConstruction.cc.
1687{ ( *besEMCGeometry ).BSCRmin = val; }
◆ SetBSCRmin() [2/3]
| void BesEmcConstruction::SetBSCRmin | ( | G4double | ) |
◆ SetBSCRmin() [3/3]
| void BesEmcConstruction::SetBSCRmin | ( | G4double | ) |
◆ SetBSCYFront() [1/3]
| void BesEmcConstruction::SetBSCYFront | ( | G4double | val | ) |
Definition at line 1713 of file BesEmcConstruction.cc.
1713{ ( *besEMCGeometry ).BSCYFront = val; }
◆ SetBSCYFront() [2/3]
| void BesEmcConstruction::SetBSCYFront | ( | G4double | ) |
◆ SetBSCYFront() [3/3]
| void BesEmcConstruction::SetBSCYFront | ( | G4double | ) |
◆ SetBSCYFront0() [1/3]
| void BesEmcConstruction::SetBSCYFront0 | ( | G4double | val | ) |
Definition at line 1707 of file BesEmcConstruction.cc.
1707 {
1708 ( *besEMCGeometry ).BSCYFront0 = val;
1709}
◆ SetBSCYFront0() [2/3]
| void BesEmcConstruction::SetBSCYFront0 | ( | G4double | ) |
◆ SetBSCYFront0() [3/3]
| void BesEmcConstruction::SetBSCYFront0 | ( | G4double | ) |
◆ SetCasingMaterial() [1/3]
| void BesEmcConstruction::SetCasingMaterial | ( | G4String | materialChoice | ) |
Definition at line 1665 of file BesEmcConstruction.cc.
1665 {
1666 // search the material by its name
1667 G4Material* pttoMaterial = G4Material::GetMaterial( materialChoice );
1668 if ( pttoMaterial )
1669 {
1670 fCasingMaterial = pttoMaterial;
1671 logicBSCTheta->SetMaterial( pttoMaterial );
1672 PrintEMCParameters();
1673 }
1674}
◆ SetCasingMaterial() [2/3]
| void BesEmcConstruction::SetCasingMaterial | ( | G4String | ) |
◆ SetCasingMaterial() [3/3]
| void BesEmcConstruction::SetCasingMaterial | ( | G4String | ) |
◆ SetCasingThickness() [1/3]
| void BesEmcConstruction::SetCasingThickness | ( | G4ThreeVector | val | ) |
Definition at line 1678 of file BesEmcConstruction.cc.
1678 {
1679 // change Gap thickness and recompute the calorimeter parameters
1680 ( *besEMCGeometry ).fTyvekThickness = val( 'X' );
1681 ( *besEMCGeometry ).fAlThickness = val( 'Y' );
1682 ( *besEMCGeometry ).fMylarThickness = val( 'Z' );
1683}
◆ SetCasingThickness() [2/3]
| void BesEmcConstruction::SetCasingThickness | ( | G4ThreeVector | ) |
◆ SetCasingThickness() [3/3]
| void BesEmcConstruction::SetCasingThickness | ( | G4ThreeVector | ) |
◆ SetCrystalMaterial() [1/3]
| void BesEmcConstruction::SetCrystalMaterial | ( | G4String | materialChoice | ) |
Definition at line 1652 of file BesEmcConstruction.cc.
1652 {
1653 // search the material by its name
1654 G4Material* pttoMaterial = G4Material::GetMaterial( materialChoice );
1655 if ( pttoMaterial )
1656 {
1657 fCrystalMaterial = pttoMaterial;
1658 logicBSCCrystal->SetMaterial( pttoMaterial );
1659 PrintEMCParameters();
1660 }
1661}
◆ SetCrystalMaterial() [2/3]
| void BesEmcConstruction::SetCrystalMaterial | ( | G4String | ) |
◆ SetCrystalMaterial() [3/3]
| void BesEmcConstruction::SetCrystalMaterial | ( | G4String | ) |
◆ SetMagField() [1/3]
| void BesEmcConstruction::SetMagField | ( | G4double | fieldValue | ) |
Definition at line 1729 of file BesEmcConstruction.cc.
1729 {
1730 // apply a global uniform magnetic field along Z axis
1731 G4FieldManager* fieldMgr =
1732 G4TransportationManager::GetTransportationManager()->GetFieldManager();
1733
1734 if ( magField ) delete magField; // delete the existing magn field
1735
1736 if ( fieldValue != 0. ) // create a new one if non nul
1737 {
1738 magField = new G4UniformMagField( G4ThreeVector( 0., 0., fieldValue ) );
1739 fieldMgr->SetDetectorField( magField );
1740 fieldMgr->CreateChordFinder( magField );
1741 fmagField = fieldValue;
1742 }
1743 else
1744 {
1745 magField = 0;
1746 fieldMgr->SetDetectorField( magField );
1747 fmagField = 0.;
1748 }
1749}
◆ SetMagField() [2/3]
| void BesEmcConstruction::SetMagField | ( | G4double | ) |
◆ SetMagField() [3/3]
| void BesEmcConstruction::SetMagField | ( | G4double | ) |
◆ SetStartIDTheta() [1/3]
| void BesEmcConstruction::SetStartIDTheta | ( | G4int | val | ) |
Definition at line 1697 of file BesEmcConstruction.cc.
1697{ startID = val; }
◆ SetStartIDTheta() [2/3]
| void BesEmcConstruction::SetStartIDTheta | ( | G4int | ) |
◆ SetStartIDTheta() [3/3]
| void BesEmcConstruction::SetStartIDTheta | ( | G4int | ) |
◆ SetVerboseLevel() [1/3]
|
inline |
Definition at line 65 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcSim/BesEmcConstruction.hh.
65{ verboseLevel = val; }
◆ SetVerboseLevel() [2/3]
|
inline |
Definition at line 65 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcSim/BesEmcConstruction.hh.
65{ verboseLevel = val; }
◆ SetVerboseLevel() [3/3]
|
inline |
Definition at line 65 of file Simulation/BOOST/EmcSim/include/EmcSim/BesEmcConstruction.hh.
65{ verboseLevel = val; }
◆ SetVisAndSD() [1/3]
| void BesEmcConstruction::SetVisAndSD | ( | ) |
Definition at line 1338 of file BesEmcConstruction.cc.
1338 {
1339 //-------------------------------------------------------------
1340 // Barrel
1341 G4VisAttributes* bscVisAtt = new G4VisAttributes( G4Colour( 0.5, 0.5, 0.5 ) );
1342 bscVisAtt->SetVisibility( false );
1343 logicEMC->SetVisAttributes( bscVisAtt );
1344 if ( logicBSCWorld ) logicBSCWorld->SetVisAttributes( G4VisAttributes::Invisible );
1345
1346 if ( logicBSCCrystal )
1347 {
1348 // G4cout<<"find BSCCrystal "<<G4endl;
1349 G4VisAttributes* crystalVisAtt = new G4VisAttributes( G4Colour( 0, 0, 1.0 ) );
1350 crystalVisAtt->SetVisibility( true );
1351 logicBSCCrystal->SetVisAttributes( crystalVisAtt );
1352 logicBSCCrystal->SetSensitiveDetector( besEMCSD );
1353 }
1354
1355 if ( logicBSCPhi )
1356 {
1357 G4VisAttributes* rightVisAtt = new G4VisAttributes( G4Colour( 1.0, 0., 1.0 ) );
1358 rightVisAtt->SetVisibility( false );
1359 logicBSCPhi->SetVisAttributes( rightVisAtt );
1360 }
1361
1362 if ( logicRear ) logicRear->SetVisAttributes( G4VisAttributes::Invisible );
1363 if ( logicOrgGlass ) logicOrgGlass->SetVisAttributes( G4VisAttributes::Invisible );
1364 if ( logicRearCasing ) logicRearCasing->SetVisAttributes( G4VisAttributes::Invisible );
1365 if ( logicAlPlate ) logicAlPlate->SetVisAttributes( G4VisAttributes::Invisible );
1366 if ( logicPD )
1367 {
1368 logicPD->SetVisAttributes( G4VisAttributes::Invisible );
1369 logicPD->SetSensitiveDetector( besEMCSD );
1370 }
1371 if ( logicPreAmpBox ) logicPreAmpBox->SetVisAttributes( G4VisAttributes::Invisible );
1372 if ( logicAirInPABox ) logicAirInPABox->SetVisAttributes( G4VisAttributes::Invisible );
1373 if ( logicHangingPlate ) logicHangingPlate->SetVisAttributes( G4VisAttributes::Invisible );
1374 if ( logicWaterPipe ) logicWaterPipe->SetVisAttributes( G4VisAttributes::Invisible );
1375
1376 //-------------------------------------------------------------
1377 // Support system
1378 G4VisAttributes* ringVisAtt = new G4VisAttributes( G4Colour( 0.5, 0.25, 0. ) );
1379 ringVisAtt->SetVisibility( false );
1380 if ( logicSupportBar ) logicSupportBar->SetVisAttributes( ringVisAtt );
1381 if ( logicSupportBar1 ) logicSupportBar1->SetVisAttributes( ringVisAtt );
1382 if ( logicEndRing ) logicEndRing->SetVisAttributes( ringVisAtt );
1383 if ( logicGear ) logicGear->SetVisAttributes( ringVisAtt );
1384 if ( logicTaperRing1 ) logicTaperRing1->SetVisAttributes( ringVisAtt );
1385 if ( logicTaperRing2 ) logicTaperRing2->SetVisAttributes( ringVisAtt );
1386 if ( logicTaperRing3 ) logicTaperRing3->SetVisAttributes( ringVisAtt );
1387
1388 //-------------------------------------------------------------
1389 // Endcap
1390 G4VisAttributes* endPhiVisAtt = new G4VisAttributes( G4Colour( 0, 1.0, 0 ) );
1391 endPhiVisAtt->SetVisibility( false );
1392 if ( logicEnd ) logicEnd->SetVisAttributes( endPhiVisAtt );
1393}
Referenced by Construct().
◆ SetVisAndSD() [2/3]
| void BesEmcConstruction::SetVisAndSD | ( | ) |
◆ SetVisAndSD() [3/3]
| void BesEmcConstruction::SetVisAndSD | ( | ) |
◆ ThreeVectorTrans() [1/3]
| void BesEmcConstruction::ThreeVectorTrans | ( | G4ThreeVector | fPnt[8], |
| double | x[8], | ||
| double | y[8], | ||
| double | z[8] ) |
Definition at line 1787 of file BesEmcConstruction.cc.
1788 {
1789 for ( int i = 0; i < 8; i++ )
1790 {
1791 x[i] = fPnt[i].x();
1792 y[i] = fPnt[i].y();
1793 z[i] = fPnt[i].z();
1794 }
1795}
◆ ThreeVectorTrans() [2/3]
| void BesEmcConstruction::ThreeVectorTrans | ( | G4ThreeVector | fPnt[8], |
| double | x[8], | ||
| double | y[8], | ||
| double | z[8] ) |
◆ ThreeVectorTrans() [3/3]
| void BesEmcConstruction::ThreeVectorTrans | ( | G4ThreeVector | fPnt[8], |
| double | x[8], | ||
| double | y[8], | ||
| double | z[8] ) |
◆ TransformToArb8() [1/3]
| void BesEmcConstruction::TransformToArb8 | ( | const G4ThreeVector | fPnt[8], |
| G4ThreeVector | newfPnt[8], | ||
| G4ThreeVector & | center, | ||
| G4ThreeVector & | rotAngle ) |
Definition at line 1759 of file BesEmcConstruction.cc.
1761 {
1762 HepPoint3D point[8];
1763 center = G4ThreeVector( 0.0, 0.0, 0.0 );
1764 for ( int i = 0; i < 8; i++ )
1765 {
1767 center += point[i];
1768 }
1769 center /= 8.0;
1770
1771 HepPlane3D bottomPlane( point[4], point[5], point[6] );
1772 HepPoint3D centerProject = bottomPlane.point( center );
1773 Hep3Vector newZ = center - centerProject;
1774
1775 rotAngle = RotAngleFromNewZ( newZ );
1776 G4RotationMatrix* g4Rot = new G4RotationMatrix();
1777 g4Rot->rotateX( rotAngle.x() );
1778 g4Rot->rotateY( rotAngle.y() );
1779
1780 G4AffineTransform* transform = new G4AffineTransform( g4Rot, center );
1781 transform->Invert();
1782 for ( int i = 0; i < 8; i++ ) { newfPnt[i] = transform->TransformPoint( fPnt[i] ); }
1783 delete g4Rot;
1784 delete transform;
1785}
HepGeom::Point3D< double > HepPoint3D
HepGeom::Plane3D< double > HepPlane3D
Hep3Vector RotAngleFromNewZ(Hep3Vector newZ)
Definition BesEmcConstruction.cc:1797
Referenced by ConstructEndGeometry().
◆ TransformToArb8() [2/3]
| void BesEmcConstruction::TransformToArb8 | ( | const G4ThreeVector | fPnt[8], |
| G4ThreeVector | newfPnt[8], | ||
| G4ThreeVector & | center, | ||
| G4ThreeVector & | rotAngle ) |
◆ TransformToArb8() [3/3]
| void BesEmcConstruction::TransformToArb8 | ( | const G4ThreeVector | fPnt[8], |
| G4ThreeVector | newfPnt[8], | ||
| G4ThreeVector & | center, | ||
| G4ThreeVector & | rotAngle ) |
◆ UpdateGeometry() [1/3]
| void BesEmcConstruction::UpdateGeometry | ( | ) |
Definition at line 1753 of file BesEmcConstruction.cc.
1753 {
1754 ; // G4RunManager::GetRunManager()->DefineWorldVolume(BesDetectorConstruction::Construct());
1755}
◆ UpdateGeometry() [2/3]
| void BesEmcConstruction::UpdateGeometry | ( | ) |
◆ UpdateGeometry() [3/3]
| void BesEmcConstruction::UpdateGeometry | ( | ) |
The documentation for this class was generated from the following files:
