CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
sim::FieldBuilder Class Reference

#include <FieldBuilder.h>

Public Member Functions

void build (CMSFieldManager *fM, G4PropagatorInField *fP)
 
void configureForVolume (const std::string &volName, edm::ParameterSet &volPSet, CMSFieldManager *fM, G4PropagatorInField *fP)
 
 FieldBuilder (const MagneticField *, const edm::ParameterSet &)
 
 ~FieldBuilder ()
 

Private Attributes

double theDelta
 
FieldtheField
 
G4Mag_UsualEqRhs * theFieldEquation
 
edm::ParameterSet thePSet
 
G4LogicalVolume * theTopVolume
 

Detailed Description

Definition at line 15 of file FieldBuilder.h.

Constructor & Destructor Documentation

◆ FieldBuilder()

FieldBuilder::FieldBuilder ( const MagneticField f,
const edm::ParameterSet p 
)

Definition at line 24 of file FieldBuilder.cc.

References f, AlCaHLTBitMon_ParallelJobs::p, theDelta, theField, and theFieldEquation.

24  : theTopVolume(nullptr), thePSet(p) {
25  theDelta = p.getParameter<double>("delta") * CLHEP::mm;
26  theField = new Field(f, theDelta);
27  theFieldEquation = nullptr;
28 }
G4LogicalVolume * theTopVolume
Definition: FieldBuilder.h:31
double f[11][100]
G4Mag_UsualEqRhs * theFieldEquation
Definition: FieldBuilder.h:30
edm::ParameterSet thePSet
Definition: FieldBuilder.h:32

◆ ~FieldBuilder()

FieldBuilder::~FieldBuilder ( )

Definition at line 30 of file FieldBuilder.cc.

30 {}

Member Function Documentation

◆ build()

void FieldBuilder::build ( CMSFieldManager fM,
G4PropagatorInField *  fP 
)

Definition at line 32 of file FieldBuilder.cc.

References configureForVolume(), edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and thePSet.

Referenced by RunManagerMTWorker::initializeG4(), and GeometryProducer::updateMagneticField().

32  {
33  edm::ParameterSet thePSetForGMFM = thePSet.getParameter<edm::ParameterSet>("ConfGlobalMFM");
34  std::string volName = thePSetForGMFM.getParameter<std::string>("Volume");
35  edm::ParameterSet volPSet = thePSetForGMFM.getParameter<edm::ParameterSet>(volName);
36 
37  configureForVolume(volName, volPSet, fM, fP);
38 
39  edm::LogVerbatim("SimG4CoreMagneticField") << " FieldBuilder::build: Global magnetic field is used";
40 }
Log< level::Info, true > LogVerbatim
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
edm::ParameterSet thePSet
Definition: FieldBuilder.h:32
void configureForVolume(const std::string &volName, edm::ParameterSet &volPSet, CMSFieldManager *fM, G4PropagatorInField *fP)
Definition: FieldBuilder.cc:42

◆ configureForVolume()

void FieldBuilder::configureForVolume ( const std::string &  volName,
edm::ParameterSet volPSet,
CMSFieldManager fM,
G4PropagatorInField *  fP 
)

Definition at line 42 of file FieldBuilder.cc.

References edm::ParameterSet::getParameter(), CMSFieldManager::InitialiseForVolume(), AlCaHLTBitMon_QueryRunRegistry::string, theDelta, theField, theFieldEquation, and theTopVolume.

Referenced by build().

45  {
46  G4LogicalVolumeStore *theStore = G4LogicalVolumeStore::GetInstance();
47  for (auto vol : *theStore) {
48  if ((std::string)vol->GetName() == volName) {
49  theTopVolume = vol;
50  break;
51  }
52  }
53 
54  std::string fieldType = volPSet.getParameter<std::string>("Type");
55  std::string stepper = volPSet.getParameter<std::string>("Stepper");
56 
57  edm::ParameterSet stpPSet = volPSet.getParameter<edm::ParameterSet>("StepperParam");
58  double minStep = stpPSet.getParameter<double>("MinStep") * CLHEP::mm;
59 
60  if (stepper == "CMSTDormandPrince45") {
62  } else if (stepper == "G4TDormandPrince45") {
63  theFieldEquation = new G4TMagFieldEquation<sim::Field>(theField);
64  } else {
65  theFieldEquation = new G4Mag_UsualEqRhs(theField);
66  }
67 
68  FieldStepper *dStepper = new FieldStepper(theFieldEquation, theDelta, stepper);
69  G4ChordFinder *cf = new G4ChordFinder(theField, minStep, dStepper);
70 
71  MonopoleEquation *monopoleEquation = new MonopoleEquation(theField);
72  G4MagIntegratorStepper *mStepper = new G4ClassicalRK4(monopoleEquation, 8);
73  G4ChordFinder *cfmon = new G4ChordFinder(theField, minStep, mStepper);
74 
75  fM->InitialiseForVolume(stpPSet, theField, cf, cfmon, volName, fieldType, stepper, theDelta, fP);
76 }
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
G4LogicalVolume * theTopVolume
Definition: FieldBuilder.h:31
G4Mag_UsualEqRhs * theFieldEquation
Definition: FieldBuilder.h:30
void InitialiseForVolume(const edm::ParameterSet &, sim::Field *, G4ChordFinder *cfDefault, G4ChordFinder *cfMonopole, const std::string &vol, const std::string &fieldType, const std::string &stepperName, double delta, G4PropagatorInField *)

Member Data Documentation

◆ theDelta

double sim::FieldBuilder::theDelta
private

Definition at line 33 of file FieldBuilder.h.

Referenced by configureForVolume(), and FieldBuilder().

◆ theField

Field* sim::FieldBuilder::theField
private

Definition at line 29 of file FieldBuilder.h.

Referenced by configureForVolume(), and FieldBuilder().

◆ theFieldEquation

G4Mag_UsualEqRhs* sim::FieldBuilder::theFieldEquation
private

Definition at line 30 of file FieldBuilder.h.

Referenced by configureForVolume(), and FieldBuilder().

◆ thePSet

edm::ParameterSet sim::FieldBuilder::thePSet
private

Definition at line 32 of file FieldBuilder.h.

Referenced by build().

◆ theTopVolume

G4LogicalVolume* sim::FieldBuilder::theTopVolume
private

Definition at line 31 of file FieldBuilder.h.

Referenced by configureForVolume().