119 G4Material* air = G4Material::GetMaterial(
"G4_AIR");
123 G4int NbOfChambers = 5;
124 G4double chamberSpacing = 80*cm;
126 G4double chamberWidth = 20.0*cm;
127 G4double targetLength = 5.0*cm;
129 G4double trackerLength = (NbOfChambers+1)*chamberSpacing;
131 G4double worldLength = 1.2 * (2*targetLength + trackerLength);
133 G4double targetRadius = 0.5*targetLength;
134 targetLength = 0.5*targetLength;
135 G4double trackerSize = 0.5*trackerLength;
141 G4GeometryManager::GetInstance()->SetWorldMaximumExtent(worldLength);
143 G4cout <<
"Computed tolerance = "
144 << G4GeometryTolerance::GetInstance()->GetSurfaceTolerance()/mm
149 worldLength/2,worldLength/2,worldLength/2);
171 G4ThreeVector positionTarget = G4ThreeVector(0,0,-(targetLength+trackerSize));
174 =
new G4Tubs(
"target",0.,targetRadius,targetLength,0.*deg,360.*deg);
186 G4cout <<
"Target is " << 2*targetLength/cm <<
" cm of "
191 G4ThreeVector positionTracker = G4ThreeVector(0,0,0);
194 =
new G4Tubs(
"tracker",0,trackerSize,trackerSize, 0.*deg, 360.*deg);
212 =
new G4Tubs(
"tracker",0, 100*cm, 100*cm, 0.*deg, 360.*deg);
216 G4double firstPosition = -trackerSize + chamberSpacing;
217 G4double firstLength = trackerLength/10;
218 G4double lastLength = trackerLength;
239 G4cout <<
"There are " << NbOfChambers <<
" chambers in the tracker region. "
240 <<
"\nThe chambers are " << chamberWidth/cm <<
" cm of "
242 << chamberSpacing/cm <<
" cm" << G4endl;
247 worldLV ->SetVisAttributes(boxVisAtt);
249 trackerLV ->SetVisAttributes(boxVisAtt);
261 G4double maxStep = 0.5*chamberWidth;