77 new G4Material(
"ArgonGas", z= 18.,
a= 39.95*g/mole, density= 1.782*mg/cm3);
80 new G4Material(
"Aluminum", z= 13.,
a= 26.98*g/mole, density= 2.7*g/cm3);
83 new G4Material(
"Lead", z= 82.,
a= 207.19*g/mole, density= 11.35*g/cm3);
90 G4double expHall_x = 3.0*m;
91 G4double expHall_y = 1.0*m;
92 G4double expHall_z = 1.0*m;
93 G4Box* experimentalHall_box
94 =
new G4Box(
"expHall_box",expHall_x,expHall_y,expHall_z);
96 Ar,
"expHall_log",0,0,0);
102 G4double innerRadiusOfTheTube = 0.*cm;
103 G4double outerRadiusOfTheTube = 60.*cm;
104 G4double hightOfTheTube = 50.*cm;
105 G4double startAngleOfTheTube = 0.*deg;
106 G4double spanningAngleOfTheTube = 360.*deg;
107 G4Tubs* tracker_tube =
new G4Tubs(
"tracker_tube",innerRadiusOfTheTube,
108 outerRadiusOfTheTube,hightOfTheTube,
109 startAngleOfTheTube,spanningAngleOfTheTube);
111 G4double trackerPos_x = -1.0*m;
112 G4double trackerPos_y = 0.*m;
113 G4double trackerPos_z = 0.*m;
115 G4ThreeVector(trackerPos_x,trackerPos_y,trackerPos_z),
120 G4double block_x = 1.0*m;
121 G4double block_y = 50.0*cm;
122 G4double block_z = 50.0*cm;
123 G4Box* calorimeterBlock_box =
new G4Box(
"calBlock_box",block_x,
126 Pb,
"caloBlock_log",0,0,0);
127 G4double blockPos_x = 1.0*m;
128 G4double blockPos_y = 0.0*m;
129 G4double blockPos_z = 0.0*m;
131 G4ThreeVector(blockPos_x,blockPos_y,blockPos_z),
136 G4double calo_x = 1.*cm;
137 G4double calo_y = 40.*cm;
138 G4double calo_z = 40.*cm;
139 G4Box* calorimeterLayer_box =
new G4Box(
"caloLayer_box",
140 calo_x,calo_y,calo_z);
142 Al,
"caloLayer_log",0,0,0);
143 for(G4int i=0;i<19;i++)
145 G4double caloPos_x = (i-9)*10.*cm;
146 G4double caloPos_y = 0.0*m;
147 G4double caloPos_z = 0.0*m;
149 G4ThreeVector(caloPos_x,caloPos_y,caloPos_z),
160 gSystem->Load(
"libExP02ClassesDict");
164 const G4ElementTable* eltab = G4Element::GetElementTable();
168 TFile fo(
"geo.root",
"RECREATE");
170 fo.WriteObject(geotree,
"my_geo");