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 506 of file gears.cc.

Constructor & Destructor Documentation

Detector()

Detector::Detector

(

)

Definition at line 526 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 509 of file gears.cc.

              {
    delete fCmdSetB;
    delete fCmdSrc;
    delete fCmdOut;
  }

Member Function Documentation

Construct()

G4VPhysicalVolume * Detector::Construct

(

)

called at /run/initialize

Definition at line 609 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 565 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