Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
B2bDetectorConstruction.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id$
27 //
30 
31 #include "B2bDetectorConstruction.hh"
32 #include "B2bDetectorMessenger.hh"
33 #include "B2bChamberParameterisation.hh"
34 #include "B2MagneticField.hh"
35 #include "B2TrackerSD.hh"
36 
37 #include "G4Material.hh"
38 #include "G4NistManager.hh"
39 
40 #include "G4Box.hh"
41 #include "G4Tubs.hh"
42 #include "G4LogicalVolume.hh"
43 #include "G4PVPlacement.hh"
44 #include "G4PVParameterised.hh"
45 
46 #include "G4SDManager.hh"
47 
48 #include "G4GeometryTolerance.hh"
49 #include "G4GeometryManager.hh"
50 
51 #include "G4UserLimits.hh"
52 
53 #include "G4VisAttributes.hh"
54 #include "G4Colour.hh"
55 
56 #include "G4SystemOfUnits.hh"
57 
58 //#include "G4ios.hh"
59 
60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
61 
62 B2bDetectorConstruction::B2bDetectorConstruction()
63 :
64  fLogicTarget(NULL), fLogicChamber(NULL),
65  fTargetMaterial(NULL), fChamberMaterial(NULL),
66  fStepLimit(NULL),
67  fCheckOverlaps(true)
68 {
69  fMessenger = new B2bDetectorMessenger(this);
70  fMagField = new B2MagneticField();
71 }
72 
73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
74 
75 B2bDetectorConstruction::~B2bDetectorConstruction()
76 {
77  delete fMagField;
78  delete fStepLimit;
79  delete fMessenger;
80 }
81 
82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
83 
84 G4VPhysicalVolume* B2bDetectorConstruction::Construct()
85 {
86  // Define materials
87  DefineMaterials();
88 
89  // Define volumes
90  return DefineVolumes();
91 }
92 
93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
94 
95 void B2bDetectorConstruction::DefineMaterials()
96 {
97  // Material definition
98 
99  G4NistManager* nistManager = G4NistManager::Instance();
100 
101  G4bool fromIsotopes = false;
102 
103  // Air defined using NIST Manager
104  nistManager->FindOrBuildMaterial("G4_AIR", fromIsotopes);
105 
106  // Lead defined using NIST Manager
107  fTargetMaterial = nistManager->FindOrBuildMaterial("G4_Pb", fromIsotopes);
108 
109  // Xenon gas defined using NIST Manager
110  fChamberMaterial = nistManager->FindOrBuildMaterial("G4_Xe", fromIsotopes);
111 
112  // Print materials
113  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
114 }
115 
116 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
117 
118 G4VPhysicalVolume* B2bDetectorConstruction::DefineVolumes()
119 {
120  G4Material* air = G4Material::GetMaterial("G4_AIR");
121 
122  // Sizes of the principal geometrical components (solids)
123 
124  G4int NbOfChambers = 5;
125  G4double chamberSpacing = 80*cm; // from chamber center to center!
126 
127  G4double chamberWidth = 20.0*cm; // width of the chambers
128  G4double targetLength = 5.0*cm; // full length of Target
129 
130  G4double trackerLength = (NbOfChambers+1)*chamberSpacing;
131 
132  G4double worldLength = 1.2 * (2*targetLength + trackerLength);
133 
134  G4double targetRadius = 0.5*targetLength; // Radius of Target
135  targetLength = 0.5*targetLength; // Half length of the Target
136  G4double trackerSize = 0.5*trackerLength; // Half length of the Tracker
137 
138  // Definitions of Solids, Logical Volumes, Physical Volumes
139 
140  // World
141 
142  G4GeometryManager::GetInstance()->SetWorldMaximumExtent(worldLength);
143 
144  G4cout << "Computed tolerance = "
145  << G4GeometryTolerance::GetInstance()->GetSurfaceTolerance()/mm
146  << " mm" << G4endl;
147 
148  G4Box* worldS
149  = new G4Box("world", //its name
150  worldLength/2,worldLength/2,worldLength/2); //its size
151  G4LogicalVolume* worldLV
152  = new G4LogicalVolume(
153  worldS, //its solid
154  air, //its material
155  "World"); //its name
156 
157  // Must place the World Physical volume unrotated at (0,0,0).
158  //
159  G4VPhysicalVolume* worldPV
160  = new G4PVPlacement(
161  0, // no rotation
162  G4ThreeVector(), // at (0,0,0)
163  worldLV, // its logical volume
164  "World", // its name
165  0, // its mother volume
166  false, // no boolean operations
167  0, // copy number
168  fCheckOverlaps); // checking overlaps
169 
170  // Target
171 
172  G4ThreeVector positionTarget = G4ThreeVector(0,0,-(targetLength+trackerSize));
173 
174  G4Tubs* targetS
175  = new G4Tubs("target",0.,targetRadius,targetLength,0.*deg,360.*deg);
176  fLogicTarget
177  = new G4LogicalVolume(targetS, fTargetMaterial,"Target",0,0,0);
178  new G4PVPlacement(0, // no rotation
179  positionTarget, // at (x,y,z)
180  fLogicTarget, // its logical volume
181  "Target", // its name
182  worldLV, // its mother volume
183  false, // no boolean operations
184  0, // copy number
185  fCheckOverlaps); // checking overlaps
186 
187  G4cout << "Target is " << 2*targetLength/cm << " cm of "
188  << fTargetMaterial->GetName() << G4endl;
189 
190  // Tracker
191 
192  G4ThreeVector positionTracker = G4ThreeVector(0,0,0);
193 
194  G4Tubs* trackerS
195  = new G4Tubs("tracker",0,trackerSize,trackerSize, 0.*deg, 360.*deg);
196  G4LogicalVolume* trackerLV
197  = new G4LogicalVolume(trackerS, air, "Tracker",0,0,0);
198  new G4PVPlacement(0, // no rotation
199  positionTracker, // at (x,y,z)
200  trackerLV, // its logical volume
201  "Tracker", // its name
202  worldLV, // its mother volume
203  false, // no boolean operations
204  0, // copy number
205  fCheckOverlaps); // checking overlaps
206 
207  // Tracker segments
208 
209  // An example of Parameterised volumes
210  // Dummy values for G4Tubs -- modified by parameterised volume
211 
212  G4Tubs* chamberS
213  = new G4Tubs("tracker",0, 100*cm, 100*cm, 0.*deg, 360.*deg);
214  fLogicChamber
215  = new G4LogicalVolume(chamberS,fChamberMaterial,"Chamber",0,0,0);
216 
217  G4double firstPosition = -trackerSize + chamberSpacing;
218  G4double firstLength = trackerLength/10;
219  G4double lastLength = trackerLength;
220 
221  G4VPVParameterisation* chamberParam =
222  new B2bChamberParameterisation(
223  NbOfChambers, // NoChambers
224  firstPosition, // Z of center of first
225  chamberSpacing, // Z spacing of centers
226  chamberWidth, // chamber width
227  firstLength, // initial length
228  lastLength); // final length
229 
230  // dummy value : kZAxis -- modified by parameterised volume
231 
232  new G4PVParameterised("Chamber", // their name
233  fLogicChamber, // their logical volume
234  trackerLV, // Mother logical volume
235  kZAxis, // Are placed along this axis
236  NbOfChambers, // Number of chambers
237  chamberParam, // The parametrisation
238  fCheckOverlaps); // checking overlaps
239 
240  G4cout << "There are " << NbOfChambers << " chambers in the tracker region. "
241  << "\nThe chambers are " << chamberWidth/cm << " cm of "
242  << fChamberMaterial->GetName() << "\nThe distance between chamber is "
243  << chamberSpacing/cm << " cm" << G4endl;
244 
245  // Sensitive detectors
246 
247  G4String trackerChamberSDname = "B2/TrackerChamberSD";
248  B2TrackerSD* aTrackerSD = new B2TrackerSD(trackerChamberSDname,
249  "TrackerHitsCollection");
250  G4SDManager::GetSDMpointer()->AddNewDetector( aTrackerSD );
251  fLogicChamber->SetSensitiveDetector( aTrackerSD );
252 
253  // Visualization attributes
254 
255  G4VisAttributes* boxVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
256  worldLV ->SetVisAttributes(boxVisAtt);
257  fLogicTarget ->SetVisAttributes(boxVisAtt);
258  trackerLV ->SetVisAttributes(boxVisAtt);
259 
260  G4VisAttributes* chamberVisAtt = new G4VisAttributes(G4Colour(1.0,1.0,0.0));
261  fLogicChamber->SetVisAttributes(chamberVisAtt);
262 
263  // Example of User Limits
264  //
265  // Below is an example of how to set tracking constraints in a given
266  // logical volume
267  //
268  // Sets a max step length in the tracker region, with G4StepLimiter
269 
270  G4double maxStep = 0.5*chamberWidth;
271  fStepLimit = new G4UserLimits(maxStep);
272  trackerLV->SetUserLimits(fStepLimit);
273 
281 
282  // Always return the physical world
283 
284  return worldPV;
285 }
286 
287 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
288 
289 void B2bDetectorConstruction::SetTargetMaterial(G4String materialName)
290 {
291  G4NistManager* nistManager = G4NistManager::Instance();
292  G4bool fromIsotopes = false;
293 
294  G4Material* pttoMaterial =
295  nistManager->FindOrBuildMaterial(materialName, fromIsotopes);
296 
297  if (fTargetMaterial != pttoMaterial) {
298  if ( pttoMaterial ) {
299  fTargetMaterial = pttoMaterial;
300  if (fLogicTarget) fLogicTarget->SetMaterial(fTargetMaterial);
301  G4cout << "\n----> The target is made of " << materialName << G4endl;
302  } else {
303  G4cout << "\n--> WARNING from SetTargetMaterial : "
304  << materialName << " not found" << G4endl;
305  }
306  }
307 }
308 
309 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
310 
311 void B2bDetectorConstruction::SetChamberMaterial(G4String materialName)
312 {
313  G4NistManager* nistManager = G4NistManager::Instance();
314  G4bool fromIsotopes = false;
315 
316  G4Material* pttoMaterial =
317  nistManager->FindOrBuildMaterial(materialName, fromIsotopes);
318 
319  if (fChamberMaterial != pttoMaterial) {
320  if ( pttoMaterial ) {
321  fChamberMaterial = pttoMaterial;
322  if (fLogicChamber) fLogicChamber->SetMaterial(fChamberMaterial);
323  G4cout << "\n----> The chambers are made of " << materialName << G4endl;
324  } else {
325  G4cout << "\n--> WARNING from SetChamberMaterial : "
326  << materialName << " not found" << G4endl;
327  }
328  }
329 }
330 
331 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
332 
333 void B2bDetectorConstruction::SetMagField(G4double fieldValue)
334 {
335  fMagField->SetMagFieldValue(fieldValue);
336 }
337 
338 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
339 
340 void B2bDetectorConstruction::SetMaxStep(G4double maxStep)
341 {
342  if ((fStepLimit)&&(maxStep>0.)) fStepLimit->SetMaxAllowedStep(maxStep);
343 }
344 
345 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Generated on Sat May 25 2013 14:32:16 for Geant4 by   doxygen 1.8.4