#include "MagneticFieldStandard.h"

Inheritance diagram for mag::MagneticFieldStandard:

Public Member Functions
	MagneticFieldStandard (fhicl::ParameterSet const &pset)

	MagneticFieldStandard (MagneticFieldStandard const &)=delete

virtual	~MagneticFieldStandard ()=default

void	reconfigure (fhicl::ParameterSet const &pset)

std::vector< MagneticFieldDescription > const &	Fields () const override

size_t	NumFields () const override

MagFieldMode_t const &	UseField (size_t f) const override

std::string const &	MagnetizedVolume (size_t f) const override

G4ThreeVector const	FieldAtPoint (G4ThreeVector const &p=G4ThreeVector(0)) const override

G4ThreeVector const	UniformFieldInVolume (std::string const &volName) const override

Private Attributes
std::vector< MagneticFieldDescription >	fFieldDescriptions
	Descriptions of the fields. More...

Detailed Description

Definition at line 26 of file MagneticFieldStandard.h.

Constructor & Destructor Documentation

mag::MagneticFieldStandard::MagneticFieldStandard ( fhicl::ParameterSet const & pset )

explicit

Definition at line 26 of file MagneticFieldStandard.cxx.

References reconfigure().

   {
     this->reconfigure(pset);
   }

mag::MagneticFieldStandard::MagneticFieldStandard ( MagneticFieldStandard const & )

delete

virtual mag::MagneticFieldStandard::~MagneticFieldStandard ( )

virtualdefault

Member Function Documentation

G4ThreeVector const mag::MagneticFieldStandard::FieldAtPoint ( G4ThreeVector const & p = G4ThreeVector(0) ) const

overridevirtual

Implements mag::MagneticField.

Definition at line 68 of file MagneticFieldStandard.cxx.

References fFieldDescriptions.

   {
       // check that the input point is in the magnetized volume
       // Use the gGeoManager to determine what node the point
       // is in
       double point[3] = { p.x(), p.y(), p.z() };
       // loop over the field descriptions to see if the point is in any of them
       for(auto fd : fFieldDescriptions){
           // we found a node, see if its name is the same as
           // the volume with the field
           if(fd.fGeoVol->Contains(point)) return fd.fField;
       }
       // if we get here, we can't find a field
       return G4ThreeVector(0);
   }

std::vector<MagneticFieldDescription> const& mag::MagneticFieldStandard::Fields ( ) const

inlineoverridevirtual

Implements mag::MagneticField.

Definition at line 37 of file MagneticFieldStandard.h.

37 { return fFieldDescriptions; }

mag::MagneticFieldStandard::fFieldDescriptions

std::vector< MagneticFieldDescription > fFieldDescriptions

Descriptions of the fields.

Definition: MagneticFieldStandard.h:60

std::string const& mag::MagneticFieldStandard::MagnetizedVolume ( size_t f ) const

inlineoverridevirtual

Implements mag::MagneticField.

Definition at line 46 of file MagneticFieldStandard.h.

References f.

46 { return fFieldDescriptions[f].fVolume; }

f

TFile f

Definition: plotHisto.C:6

mag::MagneticFieldStandard::fFieldDescriptions

std::vector< MagneticFieldDescription > fFieldDescriptions

Descriptions of the fields.

Definition: MagneticFieldStandard.h:60

size_t mag::MagneticFieldStandard::NumFields ( ) const

inlineoverridevirtual

Implements mag::MagneticField.

Definition at line 40 of file MagneticFieldStandard.h.

40 { return fFieldDescriptions.size(); }

mag::MagneticFieldStandard::fFieldDescriptions

std::vector< MagneticFieldDescription > fFieldDescriptions

Descriptions of the fields.

Definition: MagneticFieldStandard.h:60

void mag::MagneticFieldStandard::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 32 of file MagneticFieldStandard.cxx.

References fFieldDescriptions, fhicl::ParameterSet::get(), and mag::kNoBFieldMode.

Referenced by MagneticFieldStandard(), and mag::MagneticFieldServiceStandard::reconfigure().

   {
       auto fieldDescriptions = pset.get<std::vector<fhicl::ParameterSet> >("FieldDescriptions");
       
       MagneticFieldDescription fieldDescription;
       for(auto itr : fieldDescriptions){
         fieldDescription.fMode   = (mag::MagFieldMode_t)(itr.get<int>("UseField"));
         // if the mode is turned to off, no point looking for a volume or 
         // trying to put a description into the fFieldDescriptions data member.
         // if all input field descriptions are set to fMode = kNoBFieldMode, then the 
         // methods to return the fields will not go into the loop over fFieldDescriptions
         // and will just return a 0 field.
         if(fieldDescription.fMode == mag::kNoBFieldMode) continue;
         fieldDescription.fVolume = itr.get<std::string>("MagnetizedVolume");
         fieldDescription.fGeoVol = gGeoManager->FindVolumeFast(fieldDescription.fVolume.c_str());
         // check that we have a good volume
         if( fieldDescription.fGeoVol == nullptr )
           throw cet::exception("MagneticField")
           << "cannot locat volume "
           << fieldDescription.fVolume
           << " in gGeoManager, bail";
         // These need to be read as types that FHICL know about, but they
         // are used by Geant, so I store them in Geant4 types.
         std::vector<double> field = itr.get<std::vector<double> >("ConstantField");
         
         // Force the dimension of the field definition
         field.resize(3);
         for(size_t i = 0; i < 3; ++i) fieldDescription.fField[i] = field[i];
         
         fFieldDescriptions.push_back(fieldDescription);
       }
       
       return;
   }

G4ThreeVector const mag::MagneticFieldStandard::UniformFieldInVolume ( std::string const & volName ) const

overridevirtual

Implements mag::MagneticField.

Definition at line 85 of file MagneticFieldStandard.cxx.

References fFieldDescriptions.

   {
       // if the input volume name is the same as the magnetized volume
       // return the uniform field
 
       for(auto fd : fFieldDescriptions){
           if (fd.fVolume.compare(volName) == 0) return fd.fField;
       }
 
       // if we get here, we can't find a field
       return G4ThreeVector(0);
   }