Detector Class Reference

Construct detector geometry. More...

Inheritance diagram for Detector:

Collaboration diagram for Detector:

Public Member Functions
	Detector ()
	~Detector ()
G4VPhysicalVolume *	Construct ()
	called at /run/initialize More...
void	SetNewValue (G4UIcommand *cmd, G4String value)
	for G4UI More...

Detailed Description

Construct detector geometry.

This uses two types of instructions to construct a detector:

• Geant4 text geometry
• GDML

It won't work together with HP neutron simulation if Geant4 version is lower than 10 because of this bug: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/hadronprocess/1242.html?inline=-1

Definition at line 413 of file gears.cc.

Constructor & Destructor Documentation

Detector()

Detector::Detector

(

)

Definition at line 430 of file gears.cc.

                  : G4VUserDetectorConstruction(), G4UImessenger(), fWorld(0)
{
#ifdef hasGDML
  fCmdOut = new G4UIcmdWithAString("/geometry/export",this);
  fCmdOut->SetGuidance("Export geometry gdml file name");
  fCmdOut->SetParameterName("gdml geometry output",false);
  fCmdOut->AvailableForStates(G4State_Idle);
#else
  fCmdOut=0;
#endif

  fCmdSrc = new G4UIcmdWithAString("/geometry/source",this);
  fCmdSrc->SetGuidance("Set geometry source file name");
  fCmdSrc->SetParameterName("text geometry input",false);
  fCmdSrc->AvailableForStates(G4State_PreInit);

  fCmdLmt = new G4UIcommand("/tracking/setStepLimit",this);
  fCmdLmt->SetGuidance("set max step length for a logical volume");
  fCmdLmt->AvailableForStates(G4State_Idle);
  G4UIparameter *p0 = new G4UIparameter("logical volume name", 's', false);
  G4UIparameter *p1 = new G4UIparameter("step length in mm", 'd', false);
  fCmdLmt->SetParameter(p0); fCmdLmt->SetParameter(p1);

  fCmdSetB = new G4UIcmdWith3VectorAndUnit("/geometry/SetB",this);
  fCmdSetB->SetGuidance("Set uniform magnetic field value.");
  fCmdSetB->SetParameterName("Bx", "By", "Bz", false);
  fCmdSetB->SetUnitCategory("Magnetic flux density");
}

~Detector()

Detector::~Detector ( )

inline

Definition at line 417 of file gears.cc.

{ delete fCmdSetB; delete fCmdSrc; delete fCmdOut; }

Member Function Documentation

Construct()

G4VPhysicalVolume * Detector::Construct

(

)

called at /run/initialize

Definition at line 510 of file gears.cc.

{
  if (fWorld==NULL) {
    G4cout<<"GEARS: no detector specified, set to a 10x10x10 m^3 box."<<G4endl;
    auto box = new G4Box("hall", 5*CLHEP::m, 5*CLHEP::m, 5*CLHEP::m);
    G4NistManager *nist = G4NistManager::Instance();
    G4Material *vacuum = nist->FindOrBuildMaterial("G4_Galactic");
    auto v = new G4LogicalVolume(box, vacuum, "hall");
    fWorld = new G4PVPlacement(0, G4ThreeVector(), v, "hall", 0, 0, 0);
  }
  return fWorld;
}

SetNewValue()

void Detector::SetNewValue	(	G4UIcommand *	cmd,
		G4String	value
	)

for G4UI

Definition at line 468 of file gears.cc.

{
  if (cmd==fCmdSetB) {
    auto field = new G4UniformMagField(0,0,0);
    field->SetFieldValue(fCmdSetB->GetNew3VectorValue(value));
    G4FieldManager* mgr = 
      G4TransportationManager::GetTransportationManager()->GetFieldManager();
    mgr->SetDetectorField(field);
    mgr->CreateChordFinder(field);
    G4cout<<"GEARS: Magnetic field is set to "<<value<<G4endl;
#ifdef hasGDML
  } else if (cmd==fCmdOut) {
    G4GDMLParser paser;
    paser.Write(value, fWorld);
#endif
  } else if (cmd==fCmdLmt) {
    istringstream iss(value); G4String vlm, limit; iss>>vlm>>limit;
    auto v = G4LogicalVolumeStore::GetInstance()->GetVolume(vlm);
    if (v) {
      v->SetUserLimits(new G4UserLimits(stof(limit)));
      G4cout<<"GEARS: max step length in "<<vlm<<": "<<limit<<" mm"<<G4endl;
    }
  } else { // cmd==fCmdSrc
    if (value.substr(value.length()-4)!="gdml") { // text geometry input
      G4tgbVolumeMgr* mgr = G4tgbVolumeMgr::GetInstance();
      mgr->AddTextFile(value);
      auto tgb = new TextDetectorBuilder;
      mgr->SetDetectorBuilder(tgb); 
      fWorld = mgr->ReadAndConstructDetector();
#ifdef hasGDML
    } else { // GDML input
      G4GDMLParser parser;
      parser.Read(value);
      fWorld=parser.GetWorldVolume();
#endif
    }
  }
}

---

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

• gears.cc