F06DetectorConstruction Class Reference

#include <Doxymodules_field.h>

Inheritance diagram for F06DetectorConstruction:

Public Member Functions
	F06DetectorConstruction ()
	~F06DetectorConstruction () override
G4VPhysicalVolume *	Construct () override
void	ConstructSDandField () override

Private Member Functions
void	DefineMaterials ()

Private Attributes
G4Material *	fVacuum = nullptr

Static Private Attributes
static G4ThreadLocal G4UniformGravityField *	fField = nullptr

Detailed Description

Definition at line 150 of file Doxymodules_field.h.

Constructor & Destructor Documentation

F06DetectorConstruction()

F06DetectorConstruction::F06DetectorConstruction

(

)

Definition at line 70 of file F06DetectorConstruction.cc.

 : fVacuum(nullptr)
{
  // materials
  DefineMaterials();
}

~F06DetectorConstruction()

F06DetectorConstruction::~F06DetectorConstruction ( )

override

Definition at line 79 of file F06DetectorConstruction.cc.

{
  delete fField;
}

Member Function Documentation

Construct()

G4VPhysicalVolume * F06DetectorConstruction::Construct ( )

override

Definition at line 97 of file F06DetectorConstruction.cc.

{
  //
  // World
  //
 
  G4double expHall_x = 1.0*m;
  G4double expHall_y = 1.0*m;
  G4double expHall_z = 1.0*m;
 
  auto solidWorld = new G4Box("World",             //its name
                   expHall_x,expHall_y,expHall_z); //its size
 
  auto logicWorld = new G4LogicalVolume(solidWorld,//its solid
                                    fVacuum,       //its material
                                    "World");      //its name
 
  auto physiWorld = new G4PVPlacement(nullptr,     //no rotation
                                  G4ThreeVector(), //at (0,0,0)
                                  logicWorld,      //its logical volume
                                  "World",         //its name
                                  nullptr,         //its mother  volume
                                  false,           //no boolean operation
                                  0);              //copy number
 
  G4double maxStep = 1.0*mm;
  G4double maxTime = 41.*s;
 
  auto  stepLimit = new G4UserLimits(maxStep,DBL_MAX,maxTime);
 
  logicWorld->SetUserLimits(stepLimit);
 
  //
  // Visualization attributes
  //
  // logicWorld->SetVisAttributes (G4VisAttributes::GetInvisible());
 
  //
  //always return the physical World
  //
  return physiWorld;
}

ConstructSDandField()

void F06DetectorConstruction::ConstructSDandField ( )

override

Definition at line 144 of file F06DetectorConstruction.cc.

{
  using StepperType = G4ClassicalRK4;
 
  if (!fField) {
 
     fField = new G4UniformGravityField();
 
     auto  equation = new G4RepleteEofM(fField);
//     G4EqGravityField* equation = new G4EqGravityField(fField);
 
     G4TransportationManager* transportMgr =
                           G4TransportationManager::GetTransportationManager();
 
     G4FieldManager* fieldManager= transportMgr->GetFieldManager();
     fieldManager->SetDetectorField(fField);
 
     const int nVar= 8;  // 12 for RepleteEofM
     auto  stepper = new StepperType(equation,nVar);
 
     G4double minStep           = 0.01*mm;
     G4ChordFinder* chordFinder = nullptr;
     if( stepper )
     {
        auto intgrDriver = new G4IntegrationDriver<StepperType>(minStep,
                                          stepper,
                                          stepper->GetNumberOfVariables());
        if( intgrDriver ){
           chordFinder = new G4ChordFinder(intgrDriver);
        }
     }
 
     // OLD -- and wrong
     // new G4ChordFinder((G4MagneticField*)fField,minStep,stepper);
 
     // Set accuracy parameters
     G4double deltaChord        = 3.0*mm;
     chordFinder->SetDeltaChord( deltaChord );
 
     G4double deltaIntersection = 0.1*mm;
     fieldManager->SetDeltaIntersection(deltaIntersection);
 
     //  Control accuracy of integration
     //
     G4double deltaOneStep = 0.01*mm;
     fieldManager->SetAccuraciesWithDeltaOneStep(deltaOneStep);
     //
     G4double epsMax       = 1.0e-4;  // Pure number -- maximum relative integration error
     G4double epsMin       = 2.5e-7;  //
     fieldManager->SetMinimumEpsilonStep(epsMin);
     fieldManager->SetMaximumEpsilonStep(epsMax);
     // The acceptable relative accuracy is calculated  as  deltaOneStep / stepsize
     //    but bounded to the interval between these values!
 
     fieldManager->SetChordFinder(chordFinder);
  }
}

DefineMaterials()

void F06DetectorConstruction::DefineMaterials ( )

private

Definition at line 86 of file F06DetectorConstruction.cc.

{
  G4NistManager* nistMan = G4NistManager::Instance();
 
  fVacuum = nistMan->FindOrBuildMaterial("G4_Galactic");
 
  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
}

Member Data Documentation

fVacuum

G4Material* F06DetectorConstruction::fVacuum = nullptr

private

Definition at line 64 of file F06DetectorConstruction.hh.

fField

G4ThreadLocal G4UniformGravityField * F06DetectorConstruction::fField = nullptr

staticprivate

Definition at line 66 of file F06DetectorConstruction.hh.

---

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

• .doxygen/Doxymodules_field.h
• extended/field/field06/include/F06DetectorConstruction.hh
• extended/field/field06/src/F06DetectorConstruction.cc