Construction of detector geometry. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par02DetectorConstruction:

Public Member Functions
	Par02DetectorConstruction ()
	A default constructor.

virtual	~Par02DetectorConstruction ()

virtual G4VPhysicalVolume *	Construct ()
	A method invoked by G4RunManager::Initialize()

virtual void	ConstructSDandField ()
	A method invoked by G4RunManager::Initialize() to construct thread local objects.

Public Attributes
std::vector< G4Region * >	fTrackerList
	A vector of the tracking detector regions.

std::vector< G4Region * >	fECalList
	A vector of the the electromagnetic calorimeter regions.

std::vector< G4Region * >	fHCalList
	A vector of the the hadronic calorimeter regions.

std::vector< G4Region * >	fMuonList
	A vector of the muon detector regions.

G4GlobalMagFieldMessenger *	fMagFieldMessenger
	Messenger of the magnetic field.

Detailed Description

Construction of detector geometry.

A mandatory initialization class of the detector setup. Detector construction allows to use the geometry read from a GDML file. Based on G4 examples/persistency/gdml/G01/include/G01DetectorConstruction.hh .

Author: Anna Zaborowska

Definition at line 37 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

◆ Par02DetectorConstruction()

Par02DetectorConstruction::Par02DetectorConstruction ( )

default

A default constructor.

◆ ~Par02DetectorConstruction()

Par02DetectorConstruction::~Par02DetectorConstruction ( )

virtualdefault

Member Function Documentation

◆ Construct()

G4VPhysicalVolume * Par02DetectorConstruction::Construct ( )

virtual

A method invoked by G4RunManager::Initialize()

Returns: A pointer to the world volume.

Definition at line 49 of file Par02DetectorConstruction.cc.

                                                        {
  G4GDMLParser parser;
  parser.Read( "Par02FullDetector.gdml" );
  G4cout << "Geometry loaded from  file .......Par02FullDetector.gdml " << G4endl;
 
  // This GDML detector description uses the auxiliary information part to store
  // information regarding which Geant4 volumes have a fast simulation model.
 
  const G4GDMLAuxMapType* aAuxMap = parser.GetAuxMap();
  for ( G4GDMLAuxMapType::const_iterator iter = aAuxMap->begin();
        iter != aAuxMap->end(); ++iter ) {
    for ( G4GDMLAuxListType::const_iterator vit = (*iter).second.begin();
          vit != (*iter).second.end(); ++vit ) {
      if ( (*vit).type == "FastSimModel" ) {
        G4LogicalVolume* myvol = (*iter).first;
        if ( ( myvol->GetName() ).find( "Tracker" ) != std::string::npos ) {
          fTrackerList.push_back( new G4Region( myvol->GetName() ) );
          fTrackerList.back()->AddRootLogicalVolume( myvol );
          G4cout << G4endl << "tracker !!!" << G4endl;
        } else if ( ( myvol->GetName() ).find( "HCal" ) != std::string::npos ) {
          fHCalList.push_back( new G4Region( myvol->GetName() ) );
          fHCalList.back()->AddRootLogicalVolume( myvol );
          G4cout << G4endl << "hcal !!!" << G4endl;
        } else if ( ( myvol->GetName() ).find( "ECal" ) != std::string::npos ) {
          fECalList.push_back( new G4Region( myvol->GetName() ) );
          fECalList.back()->AddRootLogicalVolume( myvol );
          G4cout << G4endl << "ecal !!!" << G4endl;
        } else if ( ( myvol->GetName() ).find( "Muon" ) != std::string::npos ) {
          fMuonList.push_back( new G4Region( myvol->GetName() ) );
          fMuonList.back()->AddRootLogicalVolume( myvol );
        } else {
          G4cout << G4endl << "NOT A KNOWN DETECTOR !!!" << G4endl;
        }
      }
    }
  }
  for ( G4int iterTracker = 0; iterTracker < G4int( fTrackerList.size() ); 
        iterTracker++ ) {
    fTrackerList[ iterTracker ]->SetProductionCuts( new G4ProductionCuts() );
    fTrackerList[ iterTracker ]->GetProductionCuts()->SetProductionCut
      ( 1.0* ( ( *fTrackerList[ iterTracker ]->GetRootLogicalVolumeIterator() )->
          GetMaterial()->GetRadlen() ) );
    fTrackerList[ iterTracker ]->GetProductionCuts()->
      SetProductionCut( 1.0*m, idxG4GammaCut );
  }
  for ( G4int iterECal = 0; iterECal < G4int( fECalList.size() ); iterECal++ ) {
    fECalList[ iterECal ]->SetProductionCuts( new G4ProductionCuts() );
    fECalList[ iterECal ]->GetProductionCuts()->SetProductionCut
      ( 0.5* ( ( *fECalList[ iterECal ]->GetRootLogicalVolumeIterator() )->
          GetMaterial()->GetRadlen() ) );
    fECalList[ iterECal ]->GetProductionCuts()->
      SetProductionCut( 0.1*m, idxG4GammaCut );
  }
  for ( G4int iterHCal = 0; iterHCal < G4int( fHCalList.size() ); iterHCal++ ) {
    fHCalList[ iterHCal ]->SetProductionCuts( new G4ProductionCuts() );
    fHCalList[ iterHCal ]->GetProductionCuts()->SetProductionCut(
      0.5* ( ( *fHCalList[iterHCal]->GetRootLogicalVolumeIterator() )->
        GetMaterial()->GetRadlen() ) );
    fHCalList[ iterHCal ]->GetProductionCuts()->
      SetProductionCut( 1.0*m, idxG4GammaCut );
  }
 
  // Returns the pointer to the physical world.
  return parser.GetWorldVolume();
}

◆ ConstructSDandField()

void Par02DetectorConstruction::ConstructSDandField ( )

virtual

A method invoked by G4RunManager::Initialize() to construct thread local objects.

Definition at line 117 of file Par02DetectorConstruction.cc.

                                                    {
  for ( G4int iterTracker = 0; iterTracker < G4int( fTrackerList.size() ); 
        iterTracker++ ) {
    // Bound the fast simulation model for the tracker subdetector
    // to all the corresponding Geant4 regions
    Par02FastSimModelTracker* fastSimModelTracker
      = new Par02FastSimModelTracker( "fastSimModelTracker", fTrackerList[ iterTracker ],
                                      Par02DetectorParametrisation::eCMS );
                                     
    // Register the fast simulation model for deleting
    G4AutoDelete::Register(fastSimModelTracker);
  }
  for ( G4int iterECal = 0; iterECal < G4int( fECalList.size() ); iterECal++ ) {
    // Bound the fast simulation model for the electromagnetic calorimeter
    // to all the corresponding Geant4 regions
    Par02FastSimModelEMCal* fastSimModelEMCal
      = new Par02FastSimModelEMCal( "fastSimModelEMCal", fECalList[ iterECal ],
                                    Par02DetectorParametrisation::eCMS );
                                     
    // Register the fast simulation model for deleting
    G4AutoDelete::Register(fastSimModelEMCal);
  }
  for ( G4int iterHCal = 0; iterHCal < G4int( fHCalList.size() ); iterHCal++ ) {
    // Bound the fast simulation model for the hadronic calorimeter
    // to all the corresponding Geant4 regions
    Par02FastSimModelHCal* fastSimModelHCal
      = new Par02FastSimModelHCal( "fastSimModelHCal", fHCalList[ iterHCal ],
                                   Par02DetectorParametrisation::eCMS );
                                     
    // Register the fast simulation model for deleting
    G4AutoDelete::Register( fastSimModelHCal );
  }
  // Currently we don't have a fast muon simulation model to be bound
  // to all the corresponding Geant4 regions.
  // But it could be added in future, in a similar way as done above for
  // the tracker subdetector and the electromagnetic and hadronic calorimeters.
 
  // Add global magnetic field
  G4ThreeVector fieldValue = G4ThreeVector();
  fMagFieldMessenger = new G4GlobalMagFieldMessenger( fieldValue );
  fMagFieldMessenger->SetVerboseLevel(1);
}