Par02DetectorConstruction Class Reference

#include <Par02DetectorConstruction.hh>

Inheritance diagram for Par02DetectorConstruction:

Collaboration diagram for Par02DetectorConstruction:

Public Member Functions
	Par02DetectorConstruction ()
	A default constructor. More...
virtual	~Par02DetectorConstruction ()
virtual G4VPhysicalVolume *	Construct ()
virtual void	ConstructSDandField ()
	A method invoked by G4RunManager::Initialize() to construct thread local objects. More...
Public Member Functions inherited from G4VUserDetectorConstruction
	G4VUserDetectorConstruction ()
virtual	~G4VUserDetectorConstruction ()
virtual void	CloneSD ()
virtual void	CloneF ()
void	RegisterParallelWorld (G4VUserParallelWorld *)
G4int	ConstructParallelGeometries ()
void	ConstructParallelSD ()
G4int	GetNumberOfParallelWorld () const
G4VUserParallelWorld *	GetParallelWorld (G4int i) const

Public Attributes
std::vector< G4Region *>	fTrackerList
	A vector of the tracking detector regions. More...
std::vector< G4Region *>	fECalList
	A vector of the the electromagnetic calorimeter regions. More...
std::vector< G4Region *>	fHCalList
	A vector of the the hadronic calorimeter regions. More...
std::vector< G4Region *>	fMuonList
	A vector of the muon detector regions. More...
G4GlobalMagFieldMessenger *	fMagFieldMessenger
	Messenger of the magnetic field. More...

Additional Inherited Members
Protected Member Functions inherited from G4VUserDetectorConstruction
void	SetSensitiveDetector (const G4String &logVolName, G4VSensitiveDetector *aSD, G4bool multi=false)
void	SetSensitiveDetector (G4LogicalVolume *logVol, G4VSensitiveDetector *aSD)

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 51 of file Par02DetectorConstruction.hh.

Constructor & Destructor Documentation

Par02DetectorConstruction()

Par02DetectorConstruction::Par02DetectorConstruction

(

)

A default constructor.

Definition at line 42 of file Par02DetectorConstruction.cc.

{}

~Par02DetectorConstruction()

Par02DetectorConstruction::~Par02DetectorConstruction ( )

virtual

Definition at line 46 of file Par02DetectorConstruction.cc.

{}

Member Function Documentation

Construct()

G4VPhysicalVolume * Par02DetectorConstruction::Construct ( void  )

virtual

A method invoked by G4RunManager::Initialize()

Returns

A pointer to the world volume.

Implements G4VUserDetectorConstruction.

Definition at line 50 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();
}

Here is the call graph for this function:

ConstructSDandField()

void Par02DetectorConstruction::ConstructSDandField ( )

virtual

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

Reimplemented from G4VUserDetectorConstruction.

Definition at line 118 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);
}

Here is the call graph for this function:

Member Data Documentation

fECalList

std::vector< G4Region* > Par02DetectorConstruction::fECalList

A vector of the the electromagnetic calorimeter regions.

Definition at line 70 of file Par02DetectorConstruction.hh.

fHCalList

std::vector< G4Region* > Par02DetectorConstruction::fHCalList

A vector of the the hadronic calorimeter regions.

Definition at line 73 of file Par02DetectorConstruction.hh.

fMagFieldMessenger

G4GlobalMagFieldMessenger* Par02DetectorConstruction::fMagFieldMessenger

Messenger of the magnetic field.

Definition at line 79 of file Par02DetectorConstruction.hh.

fMuonList

std::vector< G4Region* > Par02DetectorConstruction::fMuonList

A vector of the muon detector regions.

Definition at line 76 of file Par02DetectorConstruction.hh.

fTrackerList

std::vector< G4Region* > Par02DetectorConstruction::fTrackerList

A vector of the tracking detector regions.

Definition at line 67 of file Par02DetectorConstruction.hh.

---

The documentation for this class was generated from the following files:

• Par02DetectorConstruction.hh
• Par02DetectorConstruction.cc