73 G4Box* areaSolid=
new G4Box(
"AREA",
76 G4Material* vacuum= G4Material::GetMaterial(
"Vacuum");
77 G4LogicalVolume* areaLV=
new G4LogicalVolume(areaSolid, vacuum,
"AREA_LV");
78 G4PVPlacement* area=
new G4PVPlacement(0, G4ThreeVector(),
"AREA_PV",
81 va=
new G4VisAttributes(G4Color(1.,1.,1.));
82 va-> SetForceWireframe(
true);
83 areaLV-> SetVisAttributes(va);
89 const G4double dvoxel= 10.*mm;
90 const G4double dl= 10.*cm;
92 G4Box* svoxel=
new G4Box(
"voxel", dvoxel, dl, dvoxel);
93 mate= G4Material::GetMaterial(
"Vacuum");
94 G4LogicalVolume* lvoxel=
new G4LogicalVolume(svoxel, mate,
"voxel");
95 va=
new G4VisAttributes(G4Color(0.,0.8,0.8));
96 va-> SetVisibility(
false);
97 lvoxel-> SetVisAttributes(va);
101 for (iz=0; iz<5; iz++) {
102 for (ix=-7; ix<=7; ix++) {
103 G4double x0= (2.*ix)*cm;
104 G4double z0= (-13.+2.*iz)*cm;
105 G4PVPlacement* pvoxel=
new 106 G4PVPlacement(0, G4ThreeVector(x0, 0., z0),
107 lvoxel,
"voxel", areaLV,
false, index);
114 G4Tubs* stube=
new G4Tubs(
"tube", 0.*mm, 19./2.*mm, dl,
116 mate= G4Material::GetMaterial(
"Al");
117 G4LogicalVolume* ltube=
new G4LogicalVolume(stube, mate,
"tube");
118 va=
new G4VisAttributes(G4Color(0.,0.8,0.8));
119 ltube-> SetVisAttributes(va);
121 G4RotationMatrix* rmtube=
new G4RotationMatrix;
122 rmtube-> rotateX(-90.*deg);
123 G4PVPlacement* ptube=
new 124 G4PVPlacement(rmtube, G4ThreeVector(),
125 ltube,
"tube", lvoxel,
false, 0);
128 const G4double dxycal= 25.*mm;
129 const G4double dzcal= 3.*cm;
131 G4Box* scal=
new G4Box(
"cal", dxycal, dxycal, dzcal);
132 mate= G4Material::GetMaterial(
"CsI");
133 G4LogicalVolume* lcal=
new G4LogicalVolume(scal, mate,
"cal");
134 va=
new G4VisAttributes(G4Color(0.5,0.5,0.));
135 lcal-> SetVisAttributes(va);
138 for (ix=-2; ix<=2; ix++) {
139 G4double x0= (5.*ix)*cm;
140 G4PVPlacement* pcal=
new 141 G4PVPlacement(0, G4ThreeVector(x0, 0., 2.*cm),
142 lcal,
"cal", areaLV,
false, index);
static const G4double DXYZ_AREA