#include "DetectorConstruction.h"

Inheritance diagram for g4b::DetectorConstruction:

Public Member Functions
	DetectorConstruction (std::string const &gdmlFile, bool const &overlapCheck=false, bool const &validateSchema=true)
	Standard constructor and destructor. More...

virtual	~DetectorConstruction ()

G4VPhysicalVolume *	Construct ()

void	SetCustomMagneticService ()
	Set a custom Magnetic service. More...

Static Public Member Functions
static G4VPhysicalVolume *	GetWorld ()

static G4FieldManager *	GetFieldMgr ()
	Magnetic field. More...

Static Private Attributes
static G4VPhysicalVolume *	fWorld = nullptr
	pointer to the world volume More...

static G4FieldManager *	fFieldMgr = nullptr
	pointer to the field manager More...

Detailed Description

Definition at line 19 of file DetectorConstruction.h.

Constructor & Destructor Documentation

DetectorConstruction::DetectorConstruction	(	std::string const &	gdmlFile,
		bool const &	overlapCheck = `false`,
		bool const &	validateSchema = `true`
	)

explicit

Standard constructor and destructor.

Definition at line 32 of file DetectorConstruction.cxx.

References fWorld, and parser().

   {
     if ( gdmlFile.empty() ) {
       throw cet::exception("DetectorConstruction") << "Supplied GDML filename is empty\n"
                                                    << __FILE__ << ":" << __LINE__ << "\n";
     }
     // Get the path to the GDML file from the Geometry interface.
     const G4String GDMLfile = static_cast<const G4String>( gdmlFile );
 
     // Use Geant4's GDML parser to convert the geometry to Geant4 format.
     G4GDMLParser parser;
     parser.SetOverlapCheck(overlapCheck);
     parser.Read(GDMLfile,validateSchema);
 
     // Fetch the world physical volume from the parser.  This contains
     // the entire detector, not just the outline of the experimental
     // hall.
     fWorld = parser.GetWorldVolume();
     
   }

DetectorConstruction::~DetectorConstruction ( )

virtual

Definition at line 57 of file DetectorConstruction.cxx.

58 {

59 }

Member Function Documentation

G4VPhysicalVolume * DetectorConstruction::Construct ( )

The key method in this class; returns the Geant4 version of the detector geometry. Required of any class that inherits from G4VUserDetectorConstruction.

Definition at line 62 of file DetectorConstruction.cxx.

References fFieldMgr, fWorld, mag::kConstantBFieldMode, mag::kFieldRZMapMode, mag::kFieldXYZMapMode, mag::kNoBFieldMode, MF_LOG_ERROR, MF_LOG_INFO, and mag::MagneticFieldService::provider().

   {
     // Setup the magnetic field situation 
     art::ServiceHandle<mag::MagneticFieldService> bField;
     auto const * pProvider = bField->provider(); //get the provider
 
     // loop over the possible fields
     for(auto fd : pProvider->Fields()){
       switch (fd.fMode) {
         case mag::kNoBFieldMode:
           /* NOP */
           break;
         case mag::kConstantBFieldMode: {
           // Attach this to the magnetized volume only, so get that volume
           G4LogicalVolume *bvol = G4LogicalVolumeStore::GetInstance()->GetVolume(fd.fVolume);
           
           // Define the basic field, using p we should get the uniform field
           G4UniformMagField* magField = new G4UniformMagField( fd.fField * CLHEP::tesla );
           fFieldMgr = new G4FieldManager();
           fFieldMgr->SetDetectorField(magField);
           fFieldMgr->CreateChordFinder(magField);
           
           MF_LOG_INFO("DetectorConstruction")
           << "Setting uniform magnetic field to be "
           << magField->GetConstantFieldValue().x() << " "
           << magField->GetConstantFieldValue().y() << " "
           << magField->GetConstantFieldValue().z() << " "
           << " in " << bvol->GetName();
           
           // Reset the chord finding accuracy
           // fFieldMgr->GetChordFinder()->SetDeltaChord(1.0 * cm);
           
           // the boolean tells the field manager to use local volume
           bvol->SetFieldManager(fFieldMgr, true);
           
           break;
         } // case mag::kConstantBFieldMode
         case mag::kFieldRZMapMode: {
  
           // Attach this to the magnetized volume only, so get that volume
           G4LogicalVolume *bvol = G4LogicalVolumeStore::GetInstance()->GetVolume(fd.fVolume);
 
           mag::GlobalMagneticField *magField = new mag::GlobalMagneticField();
           fFieldMgr = new G4FieldManager();
           fFieldMgr->SetDetectorField(magField);
           fFieldMgr->CreateChordFinder(magField);
 
           MF_LOG_INFO("DetectorConstruction")
           << "Setting magnetic field in kFieldRZMapMode to be"
           << " in " << bvol->GetName();
 
           // the boolean tells the field manager to use local volume
           bvol->SetFieldManager(fFieldMgr, true);
 
           break;
         } // case mag::kFieldRZMapMode
         case mag::kFieldXYZMapMode: {
 
           // Attach this to the magnetized volume only, so get that volume
           G4LogicalVolume *bvol = G4LogicalVolumeStore::GetInstance()->GetVolume(fd.fVolume);
 
           mag::GlobalMagneticField *magField = new mag::GlobalMagneticField();
           fFieldMgr = new G4FieldManager();
           fFieldMgr->SetDetectorField(magField);
           fFieldMgr->CreateChordFinder(magField);
   
           MF_LOG_INFO("DetectorConstruction")
           << "Setting magnetic field in kFieldXYZMapMode to be"
           << " in " << bvol->GetName();
 
            // the boolean tells the field manager to use local volume
           bvol->SetFieldManager(fFieldMgr, true);
           
           break;
         } // case mag::kFieldXYZMapMode
         default: // Complain if the user asks for something not handled
           MF_LOG_ERROR("DetectorConstruction")
           << "Unknown or illegal Magneticfield "
           << "mode specified: " 
           << fd.fMode
           << ". Note that AutomaticBFieldMode is reserved.";
           break;
       } // end switch cases
       
     } // end loop over fields
     
     return fWorld;
   }

static G4FieldManager* g4b::DetectorConstruction::GetFieldMgr ( )

inlinestatic

Magnetic field.

Definition at line 43 of file DetectorConstruction.h.

References fFieldMgr.

43 { return fFieldMgr; }

g4b::DetectorConstruction::fFieldMgr

static G4FieldManager * fFieldMgr

pointer to the field manager

Definition: DetectorConstruction.h:47

static G4VPhysicalVolume* g4b::DetectorConstruction::GetWorld ( )

inlinestatic

Return a pointer to the world volume. This returns the same pointer as Construct(), but Construct() also does other actions such as setting up sensitive detectors. World() just returns the pointer to the top physical volume.

Definition at line 40 of file DetectorConstruction.h.

References fWorld.

Referenced by larg4::OpDetReadoutGeometry::Construct(), larg4::AuxDetReadoutGeometry::Construct(), and larg4::LArVoxelReadoutGeometry::Construct().

40 { return fWorld; }

g4b::DetectorConstruction::fWorld

static G4VPhysicalVolume * fWorld

pointer to the world volume

Definition: DetectorConstruction.h:46

void g4b::DetectorConstruction::SetCustomMagneticService ( )

Set a custom Magnetic service.