Skip to content
Snippets Groups Projects
Commit 78deb55a authored by clantz's avatar clantz
Browse files

Finished implementing modifiable geometry via UI command.

parent af55661c
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
include/DetectorConstruction.hh
+ 9
− 0
View file @ 78deb55a
...@@ -40,6 +40,10 @@ ...@@ -40,6 +40,10 @@
#include "G4GDMLParser.hh" #include "G4GDMLParser.hh"
#include "G4RunManager.hh" #include "G4RunManager.hh"
#ifdef CADMESH
#include "CADMesh.hh"
#endif
#include "Materials.hh" #include "Materials.hh"
#include "DetectorMessenger.hh" #include "DetectorMessenger.hh"
...@@ -57,6 +61,7 @@ class DetectorConstruction : public G4VUserDetectorConstruction ...@@ -57,6 +61,7 @@ class DetectorConstruction : public G4VUserDetectorConstruction
void OutputToGDML (G4String name); void OutputToGDML (G4String name);
void SetRotation (G4ThreeVector arg); void SetRotation (G4ThreeVector arg);
void SetTranslation (G4ThreeVector arg); void SetTranslation (G4ThreeVector arg);
void SetPMTTranslation (G4ThreeVector arg);
void SetSurfaceSigmaAlpha (G4double v); void SetSurfaceSigmaAlpha (G4double v);
void SetSurfaceFinish (const G4OpticalSurfaceFinish finish); void SetSurfaceFinish (const G4OpticalSurfaceFinish finish);
...@@ -75,6 +80,10 @@ class DetectorConstruction : public G4VUserDetectorConstruction ...@@ -75,6 +80,10 @@ class DetectorConstruction : public G4VUserDetectorConstruction
G4LogicalVolume* m_logicHalfWorld; G4LogicalVolume* m_logicHalfWorld;
G4VPhysicalVolume* m_physHalfWorld; G4VPhysicalVolume* m_physHalfWorld;
#ifdef CADMESH
CADMesh* m_mesh;
#endif
G4LogicalVolume* m_logicLightGuide; G4LogicalVolume* m_logicLightGuide;
G4VPhysicalVolume* m_physLightGuide; G4VPhysicalVolume* m_physLightGuide;
......
include/DetectorMessenger.hh
+ 1
− 0
View file @ 78deb55a
...@@ -71,6 +71,7 @@ class DetectorMessenger: public G4UImessenger{ ...@@ -71,6 +71,7 @@ class DetectorMessenger: public G4UImessenger{
G4UIcmdWithAString* fModelCmd; G4UIcmdWithAString* fModelCmd;
G4UIcmdWith3VectorAndUnit* fModelRotationCmd; G4UIcmdWith3VectorAndUnit* fModelRotationCmd;
G4UIcmdWith3VectorAndUnit* fModelTranslationCmd; G4UIcmdWith3VectorAndUnit* fModelTranslationCmd;
G4UIcmdWith3VectorAndUnit* fPMTTranslationCmd;
}; };
......
run1.mac
+ 68
− 29
View file @ 78deb55a
...@@ -15,60 +15,99 @@ ...@@ -15,60 +15,99 @@
/run/particle/verbose 0 /run/particle/verbose 0
/hits/verbose 0 /hits/verbose 0
/gps/particle opticalphoton
#Set reflective surface properties
/lightGuide/surfaceModel unified #----- Set the light guide to be used and position it -----
/lightGuide/surfaceType dielectric_metal #If no model is selected, the program will default to a
/lightGuide/surfaceFinish ground #simplified version of the 2018 testbeam light guide
/lightGuide/surfaceSigmaAlpha .1
/lightGuide/surfaceProperty SPECULARLOBECONSTANT 0.000002 1.0 0.000008 1.0 #/lightGuide/model/CADmodel models/LightGuide2018TB.stl
/lightGuide/surfaceProperty SPECULARSPIKECONSTANT 0.000002 0.3 0.000008 0.3 #/lightGuide/model/rotate 0 90 0
/lightGuide/surfaceProperty BACKSCATTERCONSTANT 0.000002 .05 0.000008 .05 #/lightGuide/model/translate 0 -250 200 mm
/lightGuide/surfaceProperty REFLECTIVITY 0.000002 .99 0.000008 .99 #/lightGuide/model/translatePMT 0 339 0 mm
# the beam energy is in gaussian profile #/lightGuide/model/CADmodel models/Winston_cone.stl
#/lightGuide/model/rotate 0 0 90
#/lightGuide/model/translate 0 -111 0 mm
#/lightGuide/model/translatePMT 0 339 0 mm
#/lightGuide/model/CADmodel models/LightGuide2007BigPMT.stl
#/lightGuide/model/rotate 0 -90 0
#/lightGuide/model/translate 0 -250 -200 mm
#/lightGuide/model/translatePMT 0 131 0 mm
#/lightGuide/model/CADmodel models/LightGuide2007SmallPMT.stl
#/lightGuide/model/rotate 0 -90 0
#/lightGuide/model/translate 0 -250 -200 mm
#/lightGuide/model/translatePMT 0 131 0 mm
#----- Set surface properties of the light guide -----
/lightGuide/surface/Model unified
/lightGuide/surface/Type dielectric_metal
/lightGuide/surface/Finish ground
/lightGuide/surface/SigmaAlpha .2
/lightGuide/surface/Property SPECULARLOBECONSTANT 0.000002 1.0 0.000008 1.0
/lightGuide/surface/Property SPECULARSPIKECONSTANT 0.000002 0.3 0.000008 0.3
/lightGuide/surface/Property BACKSCATTERCONSTANT 0.000002 .05 0.000008 .05
#/lightGuide/surface/Property REFLECTIVITY 0.000002 .89 0.000008 .89
#----- Set the refractive index of the gas the light guide is inside of -----
#/lightGuide/gasProperty RINDEX 0.000002 1.0 0.000008 1.0
#----- Set the beam energy profile -----
/gps/particle opticalphoton
/gps/ene/type Gauss /gps/ene/type Gauss
/gps/ene/mono 3.4 eV /gps/ene/mono 3.4 eV
/gps/ene/sigma 1 eV /gps/ene/sigma 1 eV
/gps/polarization 0 0 0 /gps/polarization 0 0 0
# General particle source, beam
#----- Set the beam geometry -----
# Circular beam source. Simulates a single fiber optic
#/gps/pos/type Beam #/gps/pos/type Beam
#/gps/pos/shape Circle #/gps/pos/shape Circle
#/gps/pos/centre 0. 0. 0. mm
#/gps/pos/radius 0.75 mm #/gps/pos/radius 0.75 mm
#/gps/pos/sigma_r 0.002 mm #/gps/pos/sigma_r 0.002 mm
# General particle source, rectangular plane # Rectangular plane for infinite resolution
/gps/pos/type Plane /gps/pos/type Plane
/gps/pos/shape Rectangle /gps/pos/shape Rectangle
/gps/pos/centre 0. 0. 0. mm
/gps/pos/halfx 44.875 mm /gps/pos/halfx 44.875 mm
/gps/pos/halfy 82 mm /gps/pos/halfy 82 mm
# # the incident surface is in the x-z plane just below the gas
# the incident surface is in the x-y plane
/gps/pos/rot1 1 0 0 /gps/pos/rot1 1 0 0
/gps/pos/rot2 0 1 0 /gps/pos/rot2 0 0 1
# /gps/pos/centre 0. -0.1 0. mm
#
# the beam is travelling along the z-axis # the photons are emitted around the y-axis
/gps/ang/rot1 0 1 0 /gps/ang/rot1 1 0 0
/gps/ang/rot2 1 0 0 /gps/ang/rot2 0 0 1
/gps/ang/type beam1d /gps/ang/type beam1d
# the beam angular dispersion
#----- Set beam angular dispersion for fiber optic -----
# N.A. is 0.22 so acceptance cone is ~25.4 degrees # N.A. is 0.22 so acceptance cone is ~25.4 degrees
/gps/ang/sigma_r 25.40 deg /gps/ang/sigma_r 25.40 deg
# 1.0milliradians = 0.057 deg
#
#if using these commands in interactive mode to visualize, use this to allow more events to be displayed (do not use in batch mode) #if using these commands in interactive mode to visualize, use this to allow more events to be displayed (do not use in batch mode)
/vis/ogl/set/displayListLimit 1000000 /vis/ogl/set/displayListLimit 10000
############################################################## ##############################################################
############################################################## ##############################################################
# number of events # number of events
/run/beamOn 5 /run/beamOn 100
############################################################## ##############################################################
############################################################## ##############################################################
src/DetectorConstruction.cc
+ 39
− 18
View file @ 78deb55a
...@@ -60,10 +60,6 @@ ...@@ -60,10 +60,6 @@
#include "G4GDMLParser.hh" #include "G4GDMLParser.hh"
#ifdef CADMESH
#include "CADMesh.hh"
#endif
/* /*
* *
*/ */
...@@ -79,6 +75,10 @@ DetectorConstruction::DetectorConstruction() ...@@ -79,6 +75,10 @@ DetectorConstruction::DetectorConstruction()
m_filler = G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR"); m_filler = G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR");
m_GasMPT = new G4MaterialPropertiesTable(); m_GasMPT = new G4MaterialPropertiesTable();
m_filler->SetMaterialPropertiesTable(m_GasMPT); m_filler->SetMaterialPropertiesTable(m_GasMPT);
#ifdef CADMESH
m_mesh = 0;
#endif
} }
...@@ -158,8 +158,7 @@ G4VPhysicalVolume* DetectorConstruction::Construct(){ ...@@ -158,8 +158,7 @@ G4VPhysicalVolume* DetectorConstruction::Construct(){
0, //copy number 0, //copy number
checkOverlaps); //overlaps checking checkOverlaps); //overlaps checking
G4VisAttributes* boxVisAtt_half_world = new G4VisAttributes(G4Colour(0.0,0.0,1.0,0.9)); G4VisAttributes* boxVisAtt_half_world = new G4VisAttributes(G4Colour(0.0,0.0,1.0,0.1)); // or G4Colour(0.0,0.0,1.0,0.1)
m_logicHalfWorld ->SetVisAttributes(boxVisAtt_half_world); m_logicHalfWorld ->SetVisAttributes(boxVisAtt_half_world);
...@@ -257,6 +256,9 @@ G4VPhysicalVolume* DetectorConstruction::Construct(){ ...@@ -257,6 +256,9 @@ G4VPhysicalVolume* DetectorConstruction::Construct(){
false, false,
0); 0);
G4VisAttributes* VisAtt_PMT = new G4VisAttributes(G4Colour(1.0,1.0,0.6,0.7));
m_logicPMT->SetVisAttributes(VisAtt_PMT);
m_logicPMT->SetSensitiveDetector( PMT ); m_logicPMT->SetSensitiveDetector( PMT );
return m_physWorld; return m_physWorld;
...@@ -287,13 +289,14 @@ void DetectorConstruction::UseCADModel(G4String fileName){ ...@@ -287,13 +289,14 @@ void DetectorConstruction::UseCADModel(G4String fileName){
#ifdef CADMESH #ifdef CADMESH
if(fileType == "stl"){ if(fileType == "stl"){
CADMesh* mesh = new CADMesh((char*) fileName.c_str()); if(m_mesh != 0) delete m_mesh;
mesh->SetScale(mm); m_mesh = new CADMesh((char*) fileName.c_str());
mesh->SetOffset( G4ThreeVector(-20*cm, 0, 0) ); m_mesh->SetScale(mm);
mesh->SetReverse(false); m_mesh->SetOffset( G4ThreeVector(-20*cm, 0, 0) );
m_mesh->SetReverse(false);
m_logicLightGuide = m_logicLightGuide =
new G4LogicalVolume(mesh->TessellatedMesh(), //solid new G4LogicalVolume(m_mesh->TessellatedMesh(), //solid
materials->Al, //material materials->Al, //material
"LightGuide"); //name "LightGuide"); //name
} }
...@@ -307,13 +310,21 @@ void DetectorConstruction::UseCADModel(G4String fileName){ ...@@ -307,13 +310,21 @@ void DetectorConstruction::UseCADModel(G4String fileName){
} }
if(m_logicLightGuide !=0 ){ if(m_logicLightGuide !=0 ){
//Print the extent to the console to help the user position the light guide
char message[64];
G4VisExtent extent = m_logicLightGuide->GetSolid()->GetExtent(); G4VisExtent extent = m_logicLightGuide->GetSolid()->GetExtent();
G4cout << "Xmin, Xmax " << extent.GetXmin() << ", " << extent.GetXmax() << G4endl; G4cout << "===== Light guide extent =====" << G4endl;
G4cout << "Ymin, Ymax " << extent.GetYmin() << ", " << extent.GetYmax() << G4endl; sprintf(message,"Xmin = %04.2f, Xmax = %04.2f: Width = %04.2f", extent.GetXmin(), extent.GetXmax(), extent.GetXmax() - extent.GetXmin() );
G4cout << "Zmin, Zmax " << extent.GetZmin() << ", " << extent.GetZmax() << G4endl; G4cout << message << G4endl;
sprintf(message,"Ymin = %04.2f, Ymax = %04.2f: Height = %04.2f", extent.GetYmin(), extent.GetYmax(), extent.GetYmax() - extent.GetYmin() );
G4cout << message << G4endl;
sprintf(message,"Zmin = %04.2f, Zmax = %04.2f: Depth = %04.2f", extent.GetZmin(), extent.GetZmax(), extent.GetZmax() - extent.GetZmin() );
G4cout << message << G4endl;
m_physLightGuide = m_physLightGuide =
new G4PVPlacement(0, new G4PVPlacement(0,
G4ThreeVector(0,-1*m_WorldSizeY/2,0), G4ThreeVector(0,0,0),
m_logicLightGuide, m_logicLightGuide,
"physLightGuide", "physLightGuide",
m_logicHalfWorld, m_logicHalfWorld,
...@@ -410,9 +421,9 @@ void DetectorConstruction::SetRotation(G4ThreeVector arg){ ...@@ -410,9 +421,9 @@ void DetectorConstruction::SetRotation(G4ThreeVector arg){
if(m_rotation) delete m_rotation; if(m_rotation) delete m_rotation;
m_rotation = new G4RotationMatrix(); m_rotation = new G4RotationMatrix();
m_rotation->rotateX(arg.x()*deg); m_rotation->rotateX(arg.x());
m_rotation->rotateY(arg.y()*deg); m_rotation->rotateY(arg.y());
m_rotation->rotateZ(arg.z()*deg); m_rotation->rotateZ(arg.z());
m_physLightGuide->SetRotation(m_rotation); m_physLightGuide->SetRotation(m_rotation);
m_runMan->GeometryHasBeenModified(); m_runMan->GeometryHasBeenModified();
...@@ -428,3 +439,13 @@ void DetectorConstruction::SetTranslation(G4ThreeVector arg){ ...@@ -428,3 +439,13 @@ void DetectorConstruction::SetTranslation(G4ThreeVector arg){
m_runMan->GeometryHasBeenModified(); m_runMan->GeometryHasBeenModified();
G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/rebuild"); G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/rebuild");
} }
/*
*
*/
void DetectorConstruction::SetPMTTranslation(G4ThreeVector arg){
m_physPMT->SetTranslation(arg);
m_runMan->GeometryHasBeenModified();
G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/rebuild");
}
src/DetectorMessenger.cc
+ 19
− 5
View file @ 78deb55a
...@@ -107,17 +107,26 @@ DetectorMessenger::DetectorMessenger(DetectorConstruction * Det) ...@@ -107,17 +107,26 @@ DetectorMessenger::DetectorMessenger(DetectorConstruction * Det)
fModelRotationCmd->AvailableForStates(G4State_PreInit, G4State_Idle); fModelRotationCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
fModelRotationCmd->SetToBeBroadcasted(false); fModelRotationCmd->SetToBeBroadcasted(false);
fModelRotationCmd->SetParameterName("rotationX","rotationY","rotationZ",true); fModelRotationCmd->SetParameterName("rotationX","rotationY","rotationZ",true);
//fModelRotationCmd->SetDefaultValue(G4ThreeVector(0.,0.,0.)); fModelRotationCmd->SetDefaultValue(G4ThreeVector(0.,0.,0.));
//fModelRotationCmd->SetDefaultUnit("deg"); fModelRotationCmd->SetDefaultUnit("deg");
fModelTranslationCmd = fModelTranslationCmd =
new G4UIcmdWith3VectorAndUnit("/lightGuide/model/translate", this); new G4UIcmdWith3VectorAndUnit("/lightGuide/model/translate", this);
fModelTranslationCmd->SetGuidance("Set light guide translation"); fModelTranslationCmd->SetGuidance("Set light guide translation");
fModelTranslationCmd->AvailableForStates(G4State_PreInit, G4State_Idle); fModelTranslationCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
fModelTranslationCmd->SetToBeBroadcasted(false); fModelTranslationCmd->SetToBeBroadcasted(false);
fModelRotationCmd->SetParameterName("X","Y","Z",true); fModelTranslationCmd->SetParameterName("X","Y","Z",true);
fModelRotationCmd->SetDefaultValue(G4ThreeVector(0.,0.,0.)); fModelTranslationCmd->SetDefaultValue(G4ThreeVector(0.,0.,0.));
fModelRotationCmd->SetDefaultUnit("mm"); fModelTranslationCmd->SetDefaultUnit("mm");
fPMTTranslationCmd =
new G4UIcmdWith3VectorAndUnit("/lightGuide/model/translatePMT", this);
fPMTTranslationCmd->SetGuidance("Set PMT translation");
fPMTTranslationCmd->AvailableForStates(G4State_PreInit, G4State_Idle);
fPMTTranslationCmd->SetToBeBroadcasted(false);
fPMTTranslationCmd->SetParameterName("X","Y","Z",true);
fPMTTranslationCmd->SetDefaultValue(G4ThreeVector(0.,0.,0.));
fPMTTranslationCmd->SetDefaultUnit("mm");
} }
...@@ -134,6 +143,7 @@ DetectorMessenger::~DetectorMessenger(){ ...@@ -134,6 +143,7 @@ DetectorMessenger::~DetectorMessenger(){
delete fModelCmd; delete fModelCmd;
delete fModelRotationCmd; delete fModelRotationCmd;
delete fModelTranslationCmd; delete fModelTranslationCmd;
delete fPMTTranslationCmd;
} }
/* /*
...@@ -374,4 +384,8 @@ void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue) ...@@ -374,4 +384,8 @@ void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
else if(command == fModelTranslationCmd){ else if(command == fModelTranslationCmd){
fDetector->SetTranslation(fModelTranslationCmd->GetNew3VectorValue(newValue)); fDetector->SetTranslation(fModelTranslationCmd->GetNew3VectorValue(newValue));
} }
// PMT TRANSLATION
else if(command == fPMTTranslationCmd){
fDetector->SetPMTTranslation(fPMTTranslationCmd->GetNew3VectorValue(newValue));
}
} }
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment