#include <CMSFieldManager.h>

Inheritance diagram for CMSFieldManager:

Public Member Functions
	CMSFieldManager ()

void	ConfigureForTrack (const G4Track *) override

void	InitialiseForVolume (const edm::ParameterSet &, sim::Field , G4ChordFinder , G4ChordFinder *, const std::string &vol, const std::string &type, const std::string &stepper, double delta, double minstep)

void	SetMonopoleTracking (G4bool)

	~CMSFieldManager () override

Private Member Functions
	CMSFieldManager (const CMSFieldManager &)=delete

CMSFieldManager &	operator= (const CMSFieldManager &)=delete

Private Attributes
G4ChordFinder *	chordFinder

G4ChordFinder *	chordFinderMonopole

G4ChordFinder *	currChordFinder

double	dChord

double	dChordSimple

double	dIntersection

double	dIntersectionSimple

double	dOneStep

double	dOneStepSimple

double	energyThreshold

std::unique_ptr< sim::Field >	theField

Detailed Description

Definition at line 17 of file CMSFieldManager.h.

Constructor & Destructor Documentation

CMSFieldManager::CMSFieldManager ( )

explicit

Definition at line 9 of file CMSFieldManager.cc.

   : G4FieldManager(), currChordFinder(nullptr), chordFinder(nullptr),
     chordFinderMonopole(nullptr), dChord(0.001), dOneStep(0.001),
     dIntersection(0.0001), energyThreshold(0.0), dChordSimple(0.1),
     dOneStepSimple(0.1), dIntersectionSimple(0.01)
 {}

CMSFieldManager::~CMSFieldManager ( )

override

Definition at line 16 of file CMSFieldManager.cc.

17 {}

CMSFieldManager::CMSFieldManager ( const CMSFieldManager & )

privatedelete

Member Function Documentation

void CMSFieldManager::ConfigureForTrack ( const G4Track * track )

override

Definition at line 71 of file CMSFieldManager.cc.

References chordFinder, dChord, dChordSimple, dIntersection, dIntersectionSimple, dOneStep, dOneStepSimple, and energyThreshold.

Referenced by MonopoleTransportation::AlongStepGetPhysicalInteractionLength().

 {
   // run time parameters per track
   if(track->GetKineticEnergy() <= energyThreshold && track->GetParentID() > 0) {
     chordFinder->SetDeltaChord(dChordSimple);
     SetDeltaOneStep(dOneStepSimple);
     SetDeltaIntersection(dIntersectionSimple);
   } else {
     chordFinder->SetDeltaChord(dChord);
     SetDeltaOneStep(dOneStep);
     SetDeltaIntersection(dIntersection);
   }
 } 

void CMSFieldManager::InitialiseForVolume	(	const edm::ParameterSet &	p,
		sim::Field *	field,
		G4ChordFinder *	cf,
		G4ChordFinder *	cfmon,
		const std::string &	vol,
		const std::string &	type,
		const std::string &	stepper,
		double	delta,
		double	minstep
	)

Definition at line 19 of file CMSFieldManager.cc.

References chordFinder, chordFinderMonopole, dChord, dChordSimple, dIntersection, dIntersectionSimple, dOneStep, dOneStepSimple, energyThreshold, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), GeV, createfilelist::int, SetMonopoleTracking(), and theField.

Referenced by sim::FieldBuilder::configureForVolume().

 {
   dChord = p.getParameter<double>("DeltaChord")*CLHEP::mm;
   dOneStep = p.getParameter<double>("DeltaOneStep")*CLHEP::mm;
   dIntersection = p.getParameter<double>("DeltaIntersection")*CLHEP::mm;
   energyThreshold = p.getParameter<double>("EnergyThSimple")*CLHEP::GeV;
   dChordSimple = p.getParameter<double>("DeltaChordSimple")*CLHEP::mm;
   dOneStepSimple = p.getParameter<double>("DeltaOneStepSimple")*CLHEP::mm;
   dIntersectionSimple = p.getParameter<double>("DeltaIntersectionSimple")*CLHEP::mm;
   int maxLC = (int)p.getUntrackedParameter<double>("MaximumLoopCounts",1000);
   double minEpsStep = 
     p.getUntrackedParameter<double>("MinimumEpsilonStep",0.00001)*CLHEP::mm;
   double maxEpsStep = 
     p.getUntrackedParameter<double>("MaximumEpsilonStep",0.01)*CLHEP::mm;
   edm::LogInfo("SimG4CoreApplication") 
     << " New CMSFieldManager: LogicalVolume:      <" << vol << ">\n"
     << "               Stepper:                   <" << stepper << ">\n" 
     << "               Field type                 <" << type<< ">\n"
     << "               Field const delta(mm)       " << delta << "\n"
     << "               MinStep(mm)                 " << minstep<< "\n"
     << "               DeltaChord(mm)              " << dChord<< "\n"
     << "               DeltaOneStep(mm)            " << dOneStep<< "\n"
     << "               DeltaIntersection(mm)       " << dIntersection<< "\n"
     << "               EnergyThreshold(GeV)        " << energyThreshold<< "\n"
     << "               DeltaChordSimple(mm)        " << dChord<< "\n"
     << "               DeltaOneStepSimple(mm)      " << dOneStep<< "\n"
     << "               DeltaIntersectionSimple(mm) " << dIntersection<< "\n"
     << "               MaximumLoopCounts           " << maxLC<< "\n"
     << "               MinimumEpsilonStep          " << minEpsStep<< "\n"
     << "               MaximumEpsilonStep          " << maxEpsStep;
 
   // initialisation of chord finders
   chordFinder = cf;
   chordFinderMonopole = cfmon;
   chordFinder->SetDeltaChord(dChord);
   chordFinderMonopole->SetDeltaChord(dChord);
 
   // initialisation of field manager
   theField.reset(field);
   SetDetectorField(field);
   SetDeltaOneStep(dOneStep);
   SetDeltaIntersection(dIntersection);
   SetMinimumEpsilonStep(minEpsStep);
   SetMaximumEpsilonStep(maxEpsStep);
 
   SetMonopoleTracking(false);
 }

CMSFieldManager& CMSFieldManager::operator= ( const CMSFieldManager & )

privatedelete

void CMSFieldManager::SetMonopoleTracking ( G4bool flag )

Definition at line 85 of file CMSFieldManager.cc.

References chordFinder, chordFinderMonopole, and currChordFinder.

Referenced by MonopoleTransportation::AlongStepGetPhysicalInteractionLength(), InitialiseForVolume(), and MonopoleTransportation::PostStepDoIt().

 {
   if(flag) {
     currChordFinder = chordFinderMonopole;
     SetFieldChangesEnergy(true);
   } else { 
     currChordFinder = chordFinder;
     SetFieldChangesEnergy(false);
   }
   SetChordFinder(currChordFinder);
 }