#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

Field *	theField

G4Mag_UsualEqRhs *	theFieldEquation

edm::ParameterSet	thePSet

G4LogicalVolume *	theTopVolume

Detailed Description

Definition at line 15 of file FieldBuilder.h.

Constructor & Destructor Documentation

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

Definition at line 23 of file FieldBuilder.cc.

References edm::ParameterSet::getParameter(), theDelta, theField, and theFieldEquation.

                                                                            : theTopVolume(nullptr), thePSet(p) {
   theDelta = p.getParameter<double>("delta") * CLHEP::mm;
   theField = new Field(f, theDelta);
   theFieldEquation = nullptr;
 }

FieldBuilder::~FieldBuilder ( )

Definition at line 29 of file FieldBuilder.cc.

29 {}

Member Function Documentation

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

Definition at line 31 of file FieldBuilder.cc.

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

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

                                                                      {
   edm::ParameterSet thePSetForGMFM = thePSet.getParameter<edm::ParameterSet>("ConfGlobalMFM");
   std::string volName = thePSetForGMFM.getParameter<std::string>("Volume");
   edm::ParameterSet volPSet = thePSetForGMFM.getParameter<edm::ParameterSet>(volName);
 
   configureForVolume(volName, volPSet, fM, fP);
 
   edm::LogVerbatim("SimG4CoreMagneticField") << " FieldBuilder::build: Global magnetic field is used";
 }

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

Definition at line 41 of file FieldBuilder.cc.

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

Referenced by build().

                                                                {
   G4LogicalVolumeStore *theStore = G4LogicalVolumeStore::GetInstance();
   for (auto vol : *theStore) {
     if ((std::string)vol->GetName() == volName) {
       theTopVolume = vol;
       break;
     }
   }
 
   std::string fieldType = volPSet.getParameter<std::string>("Type");
   std::string stepper = volPSet.getParameter<std::string>("Stepper");
 
   edm::ParameterSet stpPSet = volPSet.getParameter<edm::ParameterSet>("StepperParam");
   double minStep = stpPSet.getParameter<double>("MinStep") * CLHEP::mm;
 
   if (stepper == "G4TDormandPrince45") {
     theFieldEquation = new G4TMagFieldEquation<Field>(theField);
   } else {
     theFieldEquation = new G4Mag_UsualEqRhs(theField);
   }
 
   FieldStepper *dStepper = new FieldStepper(theFieldEquation, theDelta, stepper);
   G4ChordFinder *cf = new G4ChordFinder(theField, minStep, dStepper);
 
   MonopoleEquation *monopoleEquation = new MonopoleEquation(theField);
   G4MagIntegratorStepper *mStepper = new G4ClassicalRK4(monopoleEquation, 8);
   G4ChordFinder *cfmon = new G4ChordFinder(theField, minStep, mStepper);
 
   fM->InitialiseForVolume(stpPSet, theField, cf, cfmon, volName, fieldType, stepper, theDelta, fP);
 }