Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 /// \file exoticphysics/dmparticle/src/Detecto 26 /// \file exoticphysics/dmparticle/src/DetectorConstruction.cc 27 /// \brief Implementation of the DetectorConst 27 /// \brief Implementation of the DetectorConstruction class 28 // 28 // 29 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 32 33 #include "DetectorConstruction.hh" 33 #include "DetectorConstruction.hh" 34 << 35 #include "DetectorMessenger.hh" 34 #include "DetectorMessenger.hh" 36 #include "TargetSD.hh" << 37 35 >> 36 #include "G4Material.hh" 38 #include "G4Box.hh" 37 #include "G4Box.hh" 39 #include "G4FieldManager.hh" << 40 #include "G4GeometryManager.hh" << 41 #include "G4LogicalVolume.hh" 38 #include "G4LogicalVolume.hh" 42 #include "G4Material.hh" << 43 #include "G4NistManager.hh" << 44 #include "G4PVPlacement.hh" 39 #include "G4PVPlacement.hh" 45 #include "G4RunManager.hh" << 40 #include "TargetSD.hh" 46 #include "G4SDManager.hh" 41 #include "G4SDManager.hh" 47 #include "G4SystemOfUnits.hh" << 42 48 #include "G4ThreeVector.hh" << 43 #include "G4GeometryManager.hh" >> 44 49 #include "G4UnitsTable.hh" 45 #include "G4UnitsTable.hh" >> 46 #include "G4NistManager.hh" >> 47 >> 48 #include "G4FieldManager.hh" >> 49 #include "G4ThreeVector.hh" >> 50 #include "G4RunManager.hh" >> 51 #include "G4SystemOfUnits.hh" 50 52 51 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 52 54 53 DetectorConstruction::DetectorConstruction() 55 DetectorConstruction::DetectorConstruction() 54 : G4VUserDetectorConstruction(), fAbsorMater << 56 : G4VUserDetectorConstruction(), >> 57 fAbsorMaterial(nullptr), >> 58 fLogAbsor(nullptr), >> 59 fWorld(nullptr) 55 { 60 { 56 // default parameter values 61 // default parameter values 57 fAbsorSizeZ = fAbsorSizeXY = 10 * cm; 62 fAbsorSizeZ = fAbsorSizeXY = 10 * cm; 58 fWorldSizeZ = fWorldSizeXY = 1.2 * fAbsorSiz 63 fWorldSizeZ = fWorldSizeXY = 1.2 * fAbsorSizeZ; 59 64 60 SetMaterial("G4_Al"); 65 SetMaterial("G4_Al"); 61 fWorldMaterial = G4NistManager::Instance()-> << 66 fWorldMaterial = >> 67 G4NistManager::Instance()->FindOrBuildMaterial("G4_Galactic"); 62 68 63 // create commands for interactive definitio 69 // create commands for interactive definition of the detector 64 fDetectorMessenger = new DetectorMessenger(t 70 fDetectorMessenger = new DetectorMessenger(this); 65 } 71 } 66 72 67 //....oooOO0OOooo........oooOO0OOooo........oo 73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 68 74 69 DetectorConstruction::~DetectorConstruction() 75 DetectorConstruction::~DetectorConstruction() 70 { 76 { 71 delete fDetectorMessenger; 77 delete fDetectorMessenger; 72 } 78 } 73 79 74 //....oooOO0OOooo........oooOO0OOooo........oo 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 75 81 76 G4VPhysicalVolume* DetectorConstruction::Const 82 G4VPhysicalVolume* DetectorConstruction::Construct() 77 { << 83 { 78 if (fWorld) { << 84 if(fWorld) { return fWorld; } 79 return fWorld; << 80 } << 81 85 82 /**************************** World *** 86 /**************************** World *****************************/ 83 G4Box* sWorld = new G4Box("world", fWorldSiz << 87 G4Box * sWorld = new G4Box("world", >> 88 fWorldSizeXY / 2, fWorldSizeXY / 2, fWorldSizeZ / 2); 84 89 85 G4LogicalVolume* lWorld = new G4LogicalVolum << 90 G4LogicalVolume * lWorld = new G4LogicalVolume(sWorld, >> 91 fWorldMaterial, >> 92 "world"); >> 93 >> 94 fWorld = new G4PVPlacement(0, //no rotation >> 95 G4ThreeVector(), //at (0,0,0) >> 96 lWorld, //logical volume >> 97 "world", //name >> 98 0, //mother volume >> 99 false, //no boolean operation >> 100 0); //copy number 86 101 87 fWorld = new G4PVPlacement(0, // no rotatio << 88 G4ThreeVector(), << 89 lWorld, // logic << 90 "world", // name << 91 0, // mother vo << 92 false, // no boo << 93 0); // copy numb << 94 102 95 /************************** Absorber * 103 /************************** Absorber ***************************/ 96 104 97 G4Box* sAbsor = new G4Box("Absorber", fAbsor << 105 G4Box * sAbsor = new G4Box("Absorber", 98 << 106 fAbsorSizeXY / 2, fAbsorSizeXY / 2, fAbsorSizeZ / 2); 99 fLogAbsor = new G4LogicalVolume(sAbsor, fAbs << 100 107 101 new G4PVPlacement(0, // no rotation << 108 fLogAbsor = new G4LogicalVolume(sAbsor, 102 G4ThreeVector(), // at (0 << 109 fAbsorMaterial, 103 fLogAbsor, // logical vol << 110 "Absorber"); 104 "Absorber", // name << 111 105 lWorld, // mother volume << 112 new G4PVPlacement(0, //no rotation 106 false, // no boolean oper << 113 G4ThreeVector(), //at (0,0,0) 107 0); // copy number << 114 fLogAbsor, //logical volume >> 115 "Absorber", //name >> 116 lWorld, //mother volume >> 117 false, //no boolean operation >> 118 0); //copy number 108 119 109 PrintParameters(); 120 PrintParameters(); 110 121 111 /************ always return the World vo 122 /************ always return the World volume *****************/ 112 return fWorld; 123 return fWorld; 113 } 124 } 114 125 115 //....oooOO0OOooo........oooOO0OOooo........oo 126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 116 127 117 void DetectorConstruction::PrintParameters() 128 void DetectorConstruction::PrintParameters() 118 { 129 { 119 G4cout << "\n------------------------------- 130 G4cout << "\n---------------------------------------------------------\n"; 120 G4cout << "---> The Absorber is " << G4BestU << 131 G4cout << "---> The Absorber is " << G4BestUnit(fAbsorSizeZ, "Length") 121 << fAbsorMaterial->GetName() << G4end << 132 << " of " << fAbsorMaterial->GetName() << G4endl; 122 G4cout << "\n------------------------------- 133 G4cout << "\n---------------------------------------------------------\n"; >> 134 123 } 135 } 124 136 125 //....oooOO0OOooo........oooOO0OOooo........oo 137 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 126 138 127 void DetectorConstruction::SetSizeZ(G4double v 139 void DetectorConstruction::SetSizeZ(G4double value) 128 { 140 { 129 if (value > 0.0) { << 141 if(value > 0.0) { 130 fAbsorSizeZ = value; << 142 fAbsorSizeZ = value; 131 fWorldSizeZ = 1.2 * fAbsorSizeZ; 143 fWorldSizeZ = 1.2 * fAbsorSizeZ; 132 } 144 } 133 } 145 } 134 << 146 135 //....oooOO0OOooo........oooOO0OOooo........oo 147 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 136 148 137 void DetectorConstruction::SetSizeXY(G4double 149 void DetectorConstruction::SetSizeXY(G4double value) 138 { 150 { 139 if (value > 0.0) { << 151 if(value > 0.0) 140 fAbsorSizeXY = value; << 152 { >> 153 fAbsorSizeXY = value; 141 fWorldSizeXY = 1.2 * fAbsorSizeXY; 154 fWorldSizeXY = 1.2 * fAbsorSizeXY; 142 } 155 } 143 } << 156 } 144 157 145 //....oooOO0OOooo........oooOO0OOooo........oo 158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 146 159 147 void DetectorConstruction::SetMaterial(const G 160 void DetectorConstruction::SetMaterial(const G4String& namemat) 148 { 161 { 149 // search the material by its name << 162 // search the material by its name 150 163 151 G4Material* mat = G4NistManager::Instance()- 164 G4Material* mat = G4NistManager::Instance()->FindOrBuildMaterial(namemat); 152 if (!mat) { << 165 if(!mat) { 153 G4cout << "!!! DetectorConstruction::SetMa << 166 G4cout << "!!! DetectorConstruction::SetMaterial: WARNING Material <" 154 << "> does not exist in DB" << G4en << 167 << namemat << "> does not exist in DB" << G4endl; 155 return; 168 return; 156 } 169 } 157 // new material is found out 170 // new material is found out 158 if (mat != fAbsorMaterial) { 171 if (mat != fAbsorMaterial) { 159 fAbsorMaterial = mat; 172 fAbsorMaterial = mat; 160 if (fLogAbsor) { << 173 if(fLogAbsor) { fLogAbsor->SetMaterial(mat); } 161 fLogAbsor->SetMaterial(mat); << 162 } << 163 G4RunManager::GetRunManager()->PhysicsHasB 174 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 164 } 175 } 165 } 176 } 166 177 167 //....oooOO0OOooo........oooOO0OOooo........oo 178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 168 179 169 void DetectorConstruction::ConstructSDandField 180 void DetectorConstruction::ConstructSDandField() 170 { 181 { 171 auto sd = new TargetSD("Target"); 182 auto sd = new TargetSD("Target"); 172 G4SDManager::GetSDMpointer()->AddNewDetector 183 G4SDManager::GetSDMpointer()->AddNewDetector(sd); 173 SetSensitiveDetector(fLogAbsor, sd); << 184 SetSensitiveDetector(fLogAbsor, sd); 174 } 185 } 175 186 176 //....oooOO0OOooo........oooOO0OOooo........oo 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 177 188