QwGeant4
QweakSimTrap.cc
Go to the documentation of this file.
1 #include <iostream>
2 
3 #include "G4Trap.hh"
4 #include "G4SystemOfUnits.hh"
5 
6 #include "QweakSimTrap.hh"
7 
8 int main(int argc,char** argv)
9 {
10  G4double ShieldingWallHousing_FullLength_Z = 80.0*cm;
11 
12  G4double OctantCutOut_Trap_PolarAngle = -21.40*degree;
13  G4double OctantCutOut_Trap_AzimuthalAngle = 90.0*degree;
14 
15  G4double OctantCutOut_Trap_FullHeightFront = 65.8*cm;
16  G4double OctantCutOut_Trap_FullLengthFront_Outer = 79.95*2.0*cm;
17  G4double OctantCutOut_Trap_FullLengthFront_Inner = 67.9495*2.0*cm;
18 
19  G4double OctantCutOut_Trap_FullHeightBack = 59.3*cm;
20  G4double OctantCutOut_Trap_FullLengthBack_Outer = 90.45*2.0*cm;
21  G4double OctantCutOut_Trap_FullLengthBack_Inner = 79.635*2.0*cm;
22 
23 
24  G4double pDz = 0.5 * ShieldingWallHousing_FullLength_Z;
25  G4double pTheta = OctantCutOut_Trap_PolarAngle;
26  G4double pPhi = OctantCutOut_Trap_AzimuthalAngle;
27  G4double pDy1 = 0.5*OctantCutOut_Trap_FullHeightFront;
28  G4double pDx1 = 0.5*OctantCutOut_Trap_FullLengthFront_Outer;
29  G4double pDx2 = 0.5*OctantCutOut_Trap_FullLengthFront_Inner;
30  G4double pAlp1 = 0.0*degree;
31  G4double pDy2 = 0.5*OctantCutOut_Trap_FullHeightBack;
32  G4double pDx3 = 0.5*OctantCutOut_Trap_FullLengthBack_Outer;
33  G4double pDx4 = 0.5*OctantCutOut_Trap_FullLengthBack_Inner;
34  G4double pAlp2 = 0.0*degree;
35 
36  G4Trap* trap1 = new QweakSimTrap("trap1",
37  pDz, pTheta, pPhi,
38  pDy1, pDx1, pDx2, pAlp1,
39  pDy2, pDx3, pDx4, pAlp2);
40  std::cout << *trap1 << std::endl;
41 
42 }
43 
44 
main
int main(int argc, char **argv)
Definition: Hadr00.cc:69