#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	delta

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 25 of file FieldBuilder.cc.

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

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

FieldBuilder::~FieldBuilder ( )

Definition at line 33 of file FieldBuilder.cc.

34 {}

Member Function Documentation

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

Definition at line 36 of file FieldBuilder.cc.

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

Referenced by RunManagerMT::initG4().

 {    
   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::LogInfo("SimG4CoreApplication") 
     << " FieldBuilder::build: Global magnetic field is used";
 }

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

Definition at line 52 of file FieldBuilder.cc.

References delta, edm::ParameterSet::getParameter(), i, CMSFieldManager::InitialiseForVolume(), sim::FieldStepper::select(), AlCaHLTBitMon_QueryRunRegistry::string, theField, theFieldEquation, and theTopVolume.

Referenced by build().

 {
   G4LogicalVolumeStore* theStore = G4LogicalVolumeStore::GetInstance();
   for (unsigned int i=0; i<(*theStore).size(); ++i ) {
     std::string curVolName = ((*theStore)[i])->GetName();
     if ( curVolName == volName ) {
       theTopVolume = (*theStore)[i] ;
     }
   }
 
   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") ;
   int maxLoopCount = 
     (int)stpPSet.getUntrackedParameter<double>("MaximumLoopCounts",1000);
   double minEpsilonStep = 
     stpPSet.getUntrackedParameter<double>("MinimumEpsilonStep",0.00001);
   double maxEpsilonStep = 
     stpPSet.getUntrackedParameter<double>("MaximumEpsilonStep",0.01);
 
   FieldStepper * theStepper = new FieldStepper(theFieldEquation, delta);
   theStepper->select(stepper);
   G4ChordFinder * cf = new G4ChordFinder(theField,minStep,theStepper);
 
   G4MonopoleEquation* monopoleEquation = new G4MonopoleEquation(theField);
   G4MagIntegratorStepper* theStepperMon = new G4ClassicalRK4(monopoleEquation,8);
   G4ChordFinder * cfmon = new G4ChordFinder(theField,minStep,theStepperMon);
 
   fM->InitialiseForVolume(stpPSet, theField, cf, cfmon, volName, 
               fieldType, stepper, delta, minStep); 
 
   if(fP) {
     fP->SetMaxLoopCount(maxLoopCount);
     fP->SetMinimumEpsilonStep(minEpsilonStep);
     fP->SetMaximumEpsilonStep(maxEpsilonStep);
     //fP->SetVerboseLevel(0);
   }
 }