DetectorConstruction Class Reference

#include <DetectorConstruction.hh>

Inheritance diagram for DetectorConstruction:

Public Member Functions
virtual void	CloneF ()
virtual void	CloneSD ()
G4VPhysicalVolume *	Construct ()
G4int	ConstructParallelGeometries ()
void	ConstructParallelSD ()
virtual void	ConstructSDandField ()
	DetectorConstruction ()
G4Material *	GetMaterial ()
G4int	GetNumberOfParallelWorld () const
G4VUserParallelWorld *	GetParallelWorld (G4int i) const
G4double	GetSize ()
const G4VPhysicalVolume *	GetWorld ()
void	PrintParameters ()
void	RegisterParallelWorld (G4VUserParallelWorld *)
void	SetMaterial (G4String)
	~DetectorConstruction ()

Protected Member Functions
void	SetSensitiveDetector (const G4String &logVolName, G4VSensitiveDetector *aSD, G4bool multi=false)
void	SetSensitiveDetector (G4LogicalVolume *logVol, G4VSensitiveDetector *aSD)

Private Member Functions
G4VPhysicalVolume *	ConstructVolumes ()
void	DefineMaterials ()

Private Attributes
G4Material *	aMaterial
G4double	BoxSize
DetectorMessenger *	detectorMessenger
std::vector< G4VUserParallelWorld * >	parallelWorld
G4VPhysicalVolume *	pBox

Detailed Description

Definition at line 45 of file DetectorConstruction.hh.

Constructor & Destructor Documentation

DetectorConstruction()

DetectorConstruction::DetectorConstruction

(

)

Definition at line 53 of file DetectorConstruction.cc.

:pBox(0), aMaterial(0)
{
  BoxSize = 1*mm;
  DefineMaterials();
  SetMaterial("Germanium");  
  detectorMessenger = new DetectorMessenger(this);
}

References BoxSize, DefineMaterials(), detectorMessenger, mm, and SetMaterial().

~DetectorConstruction()

DetectorConstruction::~DetectorConstruction

(

)

Definition at line 64 of file DetectorConstruction.cc.

{ delete detectorMessenger;}

References detectorMessenger.

Member Function Documentation

CloneF()

void G4VUserDetectorConstruction::CloneF ( )

virtualinherited

Definition at line 126 of file G4VUserDetectorConstruction.cc.

{
  using FMtoFMmap = std::map<G4FieldManager*, G4FieldManager*>;
  using FMpair = std::pair<G4FieldManager*, G4FieldManager*>;
 
  FMtoFMmap masterToWorker;
  G4LogicalVolumeStore* const logVolStore = G4LogicalVolumeStore::GetInstance();
  for(auto it = logVolStore->cbegin(); it != logVolStore->cend(); ++it)
  {
    G4LogicalVolume* g4LogicalVolume = *it;
    // Use shadow of master to get instance of FM
    G4FieldManager* masterFM = nullptr;  // g4LogicalVolume->fFieldManager;
    G4FieldManager* clonedFM = nullptr;
    if(masterFM != nullptr)
    {
      auto fmFound = masterToWorker.find(masterFM);
      if(fmFound == masterToWorker.cend())
      {
        // First time we see this SD, let's clone and remember...
        try
        {
          auto insertedEl =
               masterToWorker.insert(FMpair(masterFM, masterFM->Clone()));
          clonedFM = (insertedEl.first)->second;
        } catch(...)
        {
          G4ExceptionDescription msg;
          msg << "Cloning of G4FieldManager failed."
              << " But derived class does not implement cloning. Cannot "
                 "continue.";
          G4Exception("G4VUserDetectorConstruction::CloneSD()", "Run0053",
                      FatalException, msg);
        }
      }
      else
      {
        // We have already seen this SD attached to a different LogicalVolume,
        // let's re-use previous clone
        clonedFM = (*fmFound).second;
      }
    }  // masterFM != 0
    // Note that we do not push FM to daughters (false argument), however, since
    // we area looping on all logical volumes and we implemented the "trick" of
    // the map master<->cloned the final effect is the same as using here the
    // correct Boolean flag: log-volumes that originally were sharing the same
    // FM they will have cloned ones
    g4LogicalVolume->SetFieldManager(clonedFM, false);
  }
}

References G4FieldManager::Clone(), FatalException, G4Exception(), G4LogicalVolumeStore::GetInstance(), second, and G4LogicalVolume::SetFieldManager().

CloneSD()

void G4VUserDetectorConstruction::CloneSD ( )

virtualinherited

Definition at line 177 of file G4VUserDetectorConstruction.cc.

{
  // Loop on ALL logial volumes to search for attached SD
  G4LogicalVolumeStore* const logVolStore = G4LogicalVolumeStore::GetInstance();
 
  using SDtoSDmap = std::map<G4VSensitiveDetector*, G4VSensitiveDetector*>;
  using SDpair = std::pair<G4VSensitiveDetector*, G4VSensitiveDetector*>;
  SDtoSDmap masterToWorker;
 
  for(auto it = logVolStore->cbegin(); it != logVolStore->cend(); ++it)
  {
    G4LogicalVolume* g4LogicalVolume = *it;
    // Use shadow of master to get the instance of SD
    G4VSensitiveDetector* masterSD = nullptr;
    G4VSensitiveDetector* clonedSD = nullptr;
    if(masterSD != nullptr)
    {
      auto sdFound = masterToWorker.find(masterSD);
      if(sdFound == masterToWorker.cend())
      {
        // First time we see this SD, let's clone and remember...
        try
        {
          auto insertedEl =
               masterToWorker.insert(SDpair(masterSD, masterSD->Clone()));
          clonedSD = (insertedEl.first)->second;
        } catch(...)
        {
          G4ExceptionDescription msg;
          msg << "Cloning of G4VSensitiveDetector requested for:"
              << masterSD->GetName() << "\n"
#ifndef WIN32
              << " (full path name: " << masterSD->GetFullPathName() << ").\n"
#endif
              << " But derived class does not implement cloning. Cannot "
                 "continue.";
          G4Exception("G4VUserDetectorConstruction::CloneSD()", "Run0053",
                      FatalException, msg);
        }
      }
      else
      {
        // We have already seen this SD attached to a different LogicalVolume,
        // let's re-use previous clone
        clonedSD = (*sdFound).second;
      }
    }  // masterSD!=0
    g4LogicalVolume->SetSensitiveDetector(clonedSD);
  }
}

References G4VSensitiveDetector::Clone(), FatalException, G4Exception(), G4VSensitiveDetector::GetFullPathName(), G4LogicalVolumeStore::GetInstance(), G4VSensitiveDetector::GetName(), second, and G4LogicalVolume::SetSensitiveDetector().

Construct()

G4VPhysicalVolume * DetectorConstruction::Construct ( )

virtual

Implements G4VUserDetectorConstruction.

Definition at line 69 of file DetectorConstruction.cc.

{
  return ConstructVolumes();
}

References ConstructVolumes().

ConstructParallelGeometries()

G4int G4VUserDetectorConstruction::ConstructParallelGeometries ( )

inherited

Definition at line 73 of file G4VUserDetectorConstruction.cc.

{
  G4int nP = 0;
  for(auto pwItr = parallelWorld.cbegin();
           pwItr != parallelWorld.cend(); ++pwItr)
  {
    (*pwItr)->Construct();
    ++nP;
  }
  return nP;
}

References G4VUserDetectorConstruction::parallelWorld.

Referenced by G4RunManager::InitializeGeometry().

ConstructParallelSD()

void G4VUserDetectorConstruction::ConstructParallelSD ( )

inherited

Definition at line 86 of file G4VUserDetectorConstruction.cc.

{
  for(auto pwItr = parallelWorld.cbegin();
           pwItr != parallelWorld.cend(); ++pwItr)
  {
    (*pwItr)->ConstructSD();
  }
}

References G4VUserDetectorConstruction::parallelWorld.

Referenced by G4RunManager::InitializeGeometry(), and G4WorkerRunManager::InitializeGeometry().

ConstructSDandField()

void G4VUserDetectorConstruction::ConstructSDandField ( )

virtualinherited

Definition at line 111 of file G4VUserDetectorConstruction.cc.

{
  // G4RunManager::RMType rmtype =
  // G4RunManager::GetRunManager()->GetRunManagerType(); if(rmtype !=
  // G4RunManager::sequentialRM)
  // {
  //  G4cout
  //  << "User-derived detector construction class does not implement \n"
  //  << "ConstructSDandFiled method: workers will not have SD and fields!\n"
  //  << "The user can safely ignore this message if (s)he has no sensitive\n"
  //  << "detector or field in her/his application." << G4endl;
  // }
}

Referenced by G4RunManager::InitializeGeometry(), and G4WorkerRunManager::InitializeGeometry().

ConstructVolumes()

G4VPhysicalVolume * DetectorConstruction::ConstructVolumes ( )

private

Definition at line 140 of file DetectorConstruction.cc.

{
  // Cleanup old geometry
  G4GeometryManager::GetInstance()->OpenGeometry();
  G4PhysicalVolumeStore::GetInstance()->Clean();
  G4LogicalVolumeStore::GetInstance()->Clean();
  G4SolidStore::GetInstance()->Clean();
 
  G4Box* 
    sBox = new G4Box("Container",           //its name
             BoxSize/2,BoxSize/2,BoxSize/2);    //its dimensions
           
  G4LogicalVolume*
    lBox = new G4LogicalVolume(sBox,            //its shape
                   aMaterial,       //its material
                   aMaterial->GetName());   //its name
 
  pBox = new G4PVPlacement(0,               //no rotation
               G4ThreeVector(),     //at (0,0,0)
                           lBox,            //its logical volume
                           aMaterial->GetName(),    //its name
                           0,               //its mother  volume
                           false,           //no boolean operation
                           0);              //copy number
                 
  //always return the root volume
  //
  return pBox;
}

References aMaterial, BoxSize, G4LogicalVolumeStore::Clean(), G4PhysicalVolumeStore::Clean(), G4SolidStore::Clean(), G4GeometryManager::GetInstance(), G4LogicalVolumeStore::GetInstance(), G4PhysicalVolumeStore::GetInstance(), G4SolidStore::GetInstance(), G4Material::GetName(), G4GeometryManager::OpenGeometry(), and pBox.

Referenced by Construct(), and SetMaterial().

DefineMaterials()

void DetectorConstruction::DefineMaterials ( )

private

H2O->SetChemicalFormula("H_2O");

Definition at line 76 of file DetectorConstruction.cc.

{
  //
  // define Elements
  //
  G4double z,a;
  
  G4Element* H  = new G4Element("Hydrogen" ,"H" , z= 1., a=   1.01*g/mole);
  G4Element* N  = new G4Element("Nitrogen" ,"N" , z= 7., a=  14.01*g/mole);
  G4Element* O  = new G4Element("Oxygen"   ,"O" , z= 8., a=  16.00*g/mole);
  G4Element* Ge = new G4Element("Germanium","Ge", z=32., a=  72.59*g/mole);
  G4Element* Bi = new G4Element("Bismuth"  ,"Bi", z=83., a= 208.98*g/mole);
  
  //
  // define materials
  //
  G4double density;
  G4int ncomponents, natoms;
  G4double fractionmass;  
  
  G4Material* Air = 
  new G4Material("Air", density= 1.290*mg/cm3, ncomponents=2);
  Air->AddElement(N, fractionmass=70.*perCent);
  Air->AddElement(O, fractionmass=30.*perCent);
 
  G4Material* H2l = 
  new G4Material("H2liquid", density= 70.8*mg/cm3, ncomponents=1);
  H2l->AddElement(H, fractionmass=1.);
  
  G4Material* H2O = 
  new G4Material("Water", density= 1.000*g/cm3, ncomponents=2);
  H2O->AddElement(H, natoms=2);
  H2O->AddElement(O, natoms=1);
  H2O->GetIonisation()->SetMeanExcitationEnergy(75.0*eV);
  
  G4Material* steam = 
  new G4Material("WaterSteam", density= 1.0*mg/cm3, ncomponents=1);
  steam->AddMaterial(H2O, fractionmass=1.);
  steam->GetIonisation()->SetMeanExcitationEnergy(71.6*eV);  
    
  G4Material* BGO = 
  new G4Material("BGO", density= 7.10*g/cm3, ncomponents=3);
  BGO->AddElement(O , natoms=12);
  BGO->AddElement(Ge, natoms= 3);
  BGO->AddElement(Bi, natoms= 4);  
    
  new G4Material("Aluminium"  , z=13., a= 26.98*g/mole, density= 2.700*g/cm3);
  new G4Material("Silicon"    , z=14., a= 28.09*g/mole, density= 2.330*g/cm3);
  new G4Material("liquidArgon", z=18., a= 39.95*g/mole, density= 1.390*g/cm3);  
  new G4Material("Iron"       , z=26., a= 55.85*g/mole, density= 7.870*g/cm3);
  new G4Material("Copper"     , z=29., a= 63.55*g/mole, density= 8.960*g/cm3);
  new G4Material("Germanium"  , z=32., a= 72.61*g/mole, density= 5.323*g/cm3);    
  new G4Material("Silver"     , z=47., a=107.87*g/mole, density= 10.50*g/cm3);
  new G4Material("Tungsten"   , z=74., a=183.85*g/mole, density= 19.30*g/cm3);
  new G4Material("Lead"       , z=82., a=207.19*g/mole, density= 11.35*g/cm3);
  new G4Material("Uranium"    , z=92., a=238.03*g/mole, density= 18.95*g/cm3);
  
 
  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
}

References G4Material::AddElement(), G4Material::AddMaterial(), cm3, eV, g, G4cout, G4endl, G4Material::GetIonisation(), G4Material::GetMaterialTable(), mg, mole, perCent, and G4IonisParamMat::SetMeanExcitationEnergy().

Referenced by DetectorConstruction().

GetMaterial()

G4Material * DetectorConstruction::GetMaterial ( )

inline

Definition at line 64 of file DetectorConstruction.hh.

{return aMaterial;};

References aMaterial.

Referenced by RunAction::BeginOfRunAction(), RunAction::CriticalEnergy(), and RunAction::GetCuts().

GetNumberOfParallelWorld()

G4int G4VUserDetectorConstruction::GetNumberOfParallelWorld ( ) const

inherited

Definition at line 96 of file G4VUserDetectorConstruction.cc.

{
  return parallelWorld.size();
}

References G4VUserDetectorConstruction::parallelWorld.

Referenced by G4VUserDetectorConstruction::GetParallelWorld().

GetParallelWorld()

G4VUserParallelWorld * G4VUserDetectorConstruction::GetParallelWorld ( G4int  i ) const

inherited

Definition at line 103 of file G4VUserDetectorConstruction.cc.

{
  if(i < 0 || i >= GetNumberOfParallelWorld())
    return nullptr;
  return parallelWorld[i];
}

References G4VUserDetectorConstruction::GetNumberOfParallelWorld(), and G4VUserDetectorConstruction::parallelWorld.

GetSize()

G4double DetectorConstruction::GetSize ( )

inline

Definition at line 63 of file DetectorConstruction.hh.

{return BoxSize;};      

References BoxSize.

Referenced by PrimaryGeneratorAction::SetDefaultKinematic().

GetWorld()

const G4VPhysicalVolume * DetectorConstruction::GetWorld ( )

inline

Definition at line 61 of file DetectorConstruction.hh.

{return pBox;};           

References pBox.

PrintParameters()

void DetectorConstruction::PrintParameters

(

)

Definition at line 172 of file DetectorConstruction.cc.

{
  G4cout << "\n The Box is " << G4BestUnit(BoxSize,"Length")
         << " of " << aMaterial->GetName() << G4endl;
}

References aMaterial, BoxSize, G4BestUnit, G4cout, G4endl, and G4Material::GetName().

RegisterParallelWorld()

void G4VUserDetectorConstruction::RegisterParallelWorld ( G4VUserParallelWorld *  aPW )

inherited

Definition at line 57 of file G4VUserDetectorConstruction.cc.

{
  auto pwItr = std::find(parallelWorld.cbegin(), parallelWorld.cend(), aPW);
  if (pwItr != parallelWorld.cend())
  {
    G4String eM = "A parallel world <";
    eM += aPW->GetName();
    eM += "> is already registered to the user detector construction.";
    G4Exception("G4VUserDetectorConstruction::RegisterParallelWorld",
                "Run0051", FatalErrorInArgument, eM);
  }
  parallelWorld.push_back(aPW);
}

References FatalErrorInArgument, G4Exception(), G4VUserParallelWorld::GetName(), and G4VUserDetectorConstruction::parallelWorld.

SetMaterial()

void DetectorConstruction::SetMaterial

(

G4String

materialChoice

)

Definition at line 182 of file DetectorConstruction.cc.

{
  // search the material by its name, or build it from nist data base
  G4Material* pttoMaterial = 
     G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
 
  if (pttoMaterial) {
    aMaterial = pttoMaterial;
    if (pBox) G4RunManager::GetRunManager()
                             ->DefineWorldVolume(ConstructVolumes());
  } else {
    G4cout << "\n--> warning from DetectorConstruction::SetMaterial : "
           << materialChoice << " not found" << G4endl;  
  }  
}

References aMaterial, ConstructVolumes(), G4RunManager::DefineWorldVolume(), G4NistManager::FindOrBuildMaterial(), G4cout, G4endl, G4RunManager::GetRunManager(), G4NistManager::Instance(), and pBox.

Referenced by BOOST_PYTHON_MODULE(), DetectorConstruction(), and DetectorMessenger::SetNewValue().

SetSensitiveDetector() [1/2]

void G4VUserDetectorConstruction::SetSensitiveDetector ( const G4String &  logVolName, G4VSensitiveDetector *  aSD, G4bool  multi = false  )

protectedinherited

Definition at line 229 of file G4VUserDetectorConstruction.cc.

{
  G4bool found                = false;
  G4LogicalVolumeStore* store = G4LogicalVolumeStore::GetInstance();
  auto volmap = store->GetMap();
  auto pos = volmap.find(logVolName);
  if(pos != volmap.cend())
  {
    if ((pos->second.size()>1) && !multi)
    {
      G4String eM = "More than one logical volumes of name <";
      eM += pos->first;
      eM += "> are found and thus the sensitive detector <";
      eM += aSD->GetName();
      eM += "> cannot be uniquely assigned.";
      G4Exception("G4VUserDetectorConstruction::SetSensitiveDetector()",
                  "Run0052", FatalErrorInArgument, eM);
    }
    found = true;
    for (std::size_t i = 0; i < pos->second.size(); ++i)
    {
      SetSensitiveDetector(pos->second[i], aSD);
    }
  }
  if(!found)
  {
    G4String eM2 = "No logical volume of name <";
    eM2 += logVolName;
    eM2 += "> is found. The specified sensitive detector <";
    eM2 += aSD->GetName();
    eM2 += "> couldn't be assigned to any volume.";
    G4Exception("G4VUserDetectorConstruction::SetSensitiveDetector()",
                "Run0053", FatalErrorInArgument, eM2);
  }
}

References FatalErrorInArgument, G4Exception(), G4LogicalVolumeStore::GetInstance(), G4LogicalVolumeStore::GetMap(), G4VSensitiveDetector::GetName(), pos, and G4VUserDetectorConstruction::SetSensitiveDetector().

Referenced by MyDetectorConstruction::SetSDtoScoreVoxel(), and G4VUserDetectorConstruction::SetSensitiveDetector().

SetSensitiveDetector() [2/2]

void G4VUserDetectorConstruction::SetSensitiveDetector ( G4LogicalVolume *  logVol, G4VSensitiveDetector *  aSD  )

protectedinherited

Definition at line 268 of file G4VUserDetectorConstruction.cc.

{
  assert(logVol != nullptr && aSD != nullptr);
 
  // The aSD has already been added by user to the manager if needed
  // G4SDManager::GetSDMpointer()->AddNewDetector(aSD);
 
  // New Logic: allow for "multiple" SDs being attached to a single LV.
  // To do that we use a special proxy SD called G4MultiSensitiveDetector
 
  // Get existing SD if already set and check if it is of the special type
  G4VSensitiveDetector* originalSD = logVol->GetSensitiveDetector();
  if(originalSD == aSD)
  {
    G4ExceptionDescription msg;
    msg << "Attempting to add multiple times the same sensitive detector (\"";
    msg << originalSD->GetName() << "\") is not allowed, skipping.";
    G4Exception("G4VUserDetectorConstruction::SetSensitiveDetector", "Run0054",
                JustWarning, msg);
    return;
  }
  if(originalSD == nullptr)
  {
    logVol->SetSensitiveDetector(aSD);
  }
  else
  {
    G4MultiSensitiveDetector* msd =
      dynamic_cast<G4MultiSensitiveDetector*>(originalSD);
    if(msd != nullptr)
    {
      msd->AddSD(aSD);
    }
    else
    {
      std::ostringstream mn;
      mn << "/MultiSD_" << logVol->GetName() << "_" << logVol;
      const G4String msdname = mn.str();
      msd                    = new G4MultiSensitiveDetector(msdname);
      // We need to register the proxy to have correct handling of IDs
      G4SDManager::GetSDMpointer()->AddNewDetector(msd);
      msd->AddSD(originalSD);
      msd->AddSD(aSD);
      logVol->SetSensitiveDetector(msd);
    }
  }
}

References G4SDManager::AddNewDetector(), G4MultiSensitiveDetector::AddSD(), G4Exception(), G4VSensitiveDetector::GetName(), G4LogicalVolume::GetName(), G4SDManager::GetSDMpointer(), G4LogicalVolume::GetSensitiveDetector(), JustWarning, and G4LogicalVolume::SetSensitiveDetector().

Field Documentation

aMaterial

G4Material* DetectorConstruction::aMaterial

private

Definition at line 73 of file DetectorConstruction.hh.

Referenced by ConstructVolumes(), GetMaterial(), PrintParameters(), and SetMaterial().

BoxSize

G4double DetectorConstruction::BoxSize

private

Definition at line 72 of file DetectorConstruction.hh.

Referenced by ConstructVolumes(), DetectorConstruction(), GetSize(), and PrintParameters().

detectorMessenger

DetectorMessenger* DetectorConstruction::detectorMessenger

private

Definition at line 75 of file DetectorConstruction.hh.

Referenced by DetectorConstruction(), and ~DetectorConstruction().

parallelWorld

std::vector<G4VUserParallelWorld*> G4VUserDetectorConstruction::parallelWorld

privateinherited

Definition at line 82 of file G4VUserDetectorConstruction.hh.

Referenced by G4VUserDetectorConstruction::ConstructParallelGeometries(), G4VUserDetectorConstruction::ConstructParallelSD(), G4VUserDetectorConstruction::GetNumberOfParallelWorld(), G4VUserDetectorConstruction::GetParallelWorld(), and G4VUserDetectorConstruction::RegisterParallelWorld().

pBox

G4VPhysicalVolume* DetectorConstruction::pBox

private

Definition at line 70 of file DetectorConstruction.hh.

Referenced by ConstructVolumes(), GetWorld(), and SetMaterial().

---

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

• environments/g4py/examples/demos/TestEm0/module/DetectorConstruction.hh
• environments/g4py/examples/demos/TestEm0/module/DetectorConstruction.cc