Hadr10/src/DetectorConstruction.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 /// \file DetectorConstruction.cc
0027 /// \brief Implementation of the DetectorConstruction class
0028 //
0029 //
0030
0031 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0032 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0033
0034 #include "DetectorConstruction.hh"
0035
0036 #include "DetectorMessenger.hh"
0037 #include "PrimaryGeneratorAction.hh"
0038
0039 #include "G4Box.hh"
0040 #include "G4FieldManager.hh"
0041 #include "G4GeometryManager.hh"
0042 #include "G4LogicalVolume.hh"
0043 #include "G4LogicalVolumeStore.hh"
0044 #include "G4Material.hh"
0045 #include "G4NistManager.hh"
0046 #include "G4PVPlacement.hh"
0047 #include "G4PhysicalConstants.hh"
0048 #include "G4PhysicalVolumeStore.hh"
0049 #include "G4RunManager.hh"
0050 #include "G4SolidStore.hh"
0051 #include "G4SystemOfUnits.hh"
0052 #include "G4ThreeVector.hh"
0053 #include "G4TransportationManager.hh"
0054 #include "G4Tubs.hh"
0055 #include "G4UniformMagField.hh"
0056 #include "globals.hh"
0057
0058 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0059
0060 DetectorConstruction::DetectorConstruction()
0061   : G4VUserDetectorConstruction(),
0062     fTargetMaterial(nullptr),
0063     fLogicExperimentalHall(nullptr),
0064     fPhysExperimentalHall(nullptr),
0065     fLogicTargetLayer(nullptr),
0066     fPhysTargetLayer(nullptr),
0067     fFieldMgr(nullptr),
0068     fUniformMagField(nullptr),
0069     fDetectorMessenger(nullptr),
0070     fTargetInnerRadius(9.0 * mm),
0071     fTargetOuterRadius(11.0 * mm)  //***LOOKHERE*** Default values
0072 {
0073   fFieldMgr = G4TransportationManager::GetTransportationManager()->GetFieldManager();
0074   //***LOOKHERE*** Default material
0075   fTargetMaterial = G4NistManager::Instance()->FindOrBuildMaterial("G4_Be");
0076   fDetectorMessenger = new DetectorMessenger(this);
0077 }
0078
0079 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0080
0081 DetectorConstruction::~DetectorConstruction()
0082 {
0083   delete fUniformMagField;
0084   delete fDetectorMessenger;
0085 }
0086
0087 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0088
0089 G4VPhysicalVolume* DetectorConstruction::Construct()
0090 {
0091   return ConstructLayer();
0092 }
0093
0094 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0095
0096 G4VPhysicalVolume* DetectorConstruction::ConstructLayer()
0097 {
0098   // Clean old geometry, if any
0099   G4GeometryManager::GetInstance()->OpenGeometry();
0100   G4PhysicalVolumeStore::GetInstance()->Clean();
0101   G4LogicalVolumeStore::GetInstance()->Clean();
0102   G4SolidStore::GetInstance()->Clean();
0103   // The geometry consists of a cylinder with axis along the z-direction, and
0104   // positioned at the center, (0.0, 0.0, 0.0). Its inner and outer radius, and
0105   // its material can be set via UI commands.
0106   // The world volume (experimental hall) is a box slightly bigger than the cylinder
0107   // and it is filled of "G4_Galactic" material.
0108   const G4double halfLength = 1.0 * m;  //***LOOKHERE*** Half-length of the cylinder
0109   const G4double expHall_x = 1.01 * halfLength;  // half dimension along x
0110   const G4double expHall_y = 1.01 * halfLength;  // half dimension along y
0111   const G4double expHall_z = 1.01 * halfLength;  // half dimension along z
0112   G4Material* vacuum = G4NistManager::Instance()->FindOrBuildMaterial("G4_Galactic");
0113   G4Box* experimentalHallBox = new G4Box("expHallBox", expHall_x, expHall_y, expHall_z);
0114   fLogicExperimentalHall = new G4LogicalVolume(experimentalHallBox,  // solid
0115                                                vacuum,  // material
0116                                                "logicExpHall",  // name
0117                                                0,  // field manager
0118                                                0,  // sensitive detector
0119                                                0);  // user limits
0120   fPhysExperimentalHall = new G4PVPlacement(0,  // rotation
0121                                             G4ThreeVector(),  // translation
0122                                             "expHall",  // name
0123                                             fLogicExperimentalHall,  // logical volume
0124                                             0,  // mother physical volume
0125                                             false,  // boolean operation
0126                                             0);  // copy number
0127   // Cylinder along the z-axis, with inner and outer diameter
0128   G4Tubs* solidTargetLayer = new G4Tubs("solidTargetLayer",
0129                                         fTargetInnerRadius,  // inner radius
0130                                         fTargetOuterRadius,  // outer radius
0131                                         halfLength,  // half cylinder length in z
0132                                         0.0,  // starting phi angle in rad
0133                                         2.0 * pi);  // final phi angle in rad
0134   fLogicTargetLayer = new G4LogicalVolume(solidTargetLayer,  // solid
0135                                           fTargetMaterial,  // material
0136                                           "logicTargetLayer",  // name
0137                                           0,  // field manager
0138                                           0,  // sensitive detector
0139                                           0);  // user limits
0140   fPhysTargetLayer = new G4PVPlacement(0,  // rotation
0141                                        G4ThreeVector(),  // translation
0142                                        "physTargetLayer",  // name
0143                                        fLogicTargetLayer,  // logical volume
0144                                        fPhysExperimentalHall,  // mother physical volume
0145                                        false,  // boolean operation
0146                                        0);  // copy number
0147   PrintParameters();
0148   G4cout << G4endl << "DetectorConstruction::ConstructLayer() : " << G4endl
0149          << "\t World (box) size: " << G4endl << "\t \t x : -/+ " << expHall_x << " mm ;"
0150          << "\t y : -/+ " << expHall_y << " mm ;"
0151          << "\t z : -/+ " << expHall_z << " mm ;" << G4endl
0152          << "\t Layer (cylinder) size : " << G4endl << "\t \t radii : " << fTargetInnerRadius
0153          << " , " << fTargetOuterRadius << " mm ;"
0154          << "\t \t length (along z) : " << 2.0 * halfLength << " mm ;" << G4endl << G4endl
0155          << G4endl;
0156   return fPhysExperimentalHall;
0157 }
0158
0159 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0160
0161 void DetectorConstruction::SetTargetMaterial(const G4String name)
0162 {
0163   fTargetMaterial = G4NistManager::Instance()->FindOrBuildMaterial(name);
0164   if (!fTargetMaterial) {
0165     G4cout << G4endl << G4endl << "WARNING: the name of the material has not been recognized!"
0166            << G4endl << "     ===> the default  * G4_Be *  will be used." << G4endl << G4endl;
0167     //***LOOKHERE*** Default material
0168     fTargetMaterial = G4NistManager::Instance()->FindOrBuildMaterial("G4_Be");
0169   }
0170   if (fLogicTargetLayer) fLogicTargetLayer->SetMaterial(fTargetMaterial);
0171 }
0172
0173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0174
0175 void DetectorConstruction::UpdateGeometry()
0176 {
0177   G4RunManager::GetRunManager()->ReinitializeGeometry();
0178   PrintParameters();
0179   // Update also the position of the gun
0180   const PrimaryGeneratorAction* pPrimaryAction = dynamic_cast<const PrimaryGeneratorAction*>(
0181     G4RunManager::GetRunManager()->GetUserPrimaryGeneratorAction());
0182   if (pPrimaryAction) pPrimaryAction->SetGunPosition();
0183 }
0184
0185 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0186
0187 void DetectorConstruction::PrintParameters()
0188 {
0189   G4cout << G4endl << G4endl << " ------  DetectorConstruction::PrintParameters() ------ " << G4endl
0190          << " Material            = " << fTargetMaterial->GetName() << G4endl
0191          << " Target Inner Radius = " << fTargetInnerRadius << " mm" << G4endl
0192          << " Target Outer Radius = " << fTargetOuterRadius << " mm" << G4endl
0193          << " B [T]               = "
0194          << (fUniformMagField ? fUniformMagField->GetConstantFieldValue() / CLHEP::tesla
0195                               : G4ThreeVector(0.0, 0.0, 0.0))
0196          << G4endl << " ------------------------------------------------------ " << G4endl
0197          << G4endl;
0198 }
0199
0200 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0201
0202 void DetectorConstruction::SetMagField(const G4double fieldValue)
0203 {
0204   if (fUniformMagField) delete fUniformMagField;
0205   if (std::abs(fieldValue) > 0.0) {
0206     // Apply a global uniform magnetic field along the Z axis
0207     fUniformMagField = new G4UniformMagField(G4ThreeVector(0.0, 0.0, fieldValue));
0208     fFieldMgr->SetDetectorField(fUniformMagField);
0209     fFieldMgr->CreateChordFinder(fUniformMagField);
0210   }
0211 }
0212
0213 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......