#include "MagneticField.h"

Public Member Functions
	MagneticField (fhicl::ParameterSet const &pset, art::ActivityRegistry &reg)

	~MagneticField ()

void	reconfigure (fhicl::ParameterSet const &pset)

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

size_t	NumFields () const

MagFieldMode_t const &	UseField (size_t f) const

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

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

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

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

Detailed Description

Definition at line 47 of file MagneticField.h.

Constructor & Destructor Documentation

MagneticField::MagneticField	(	fhicl::ParameterSet const &	pset,
		art::ActivityRegistry &	reg
	)

Definition at line 35 of file MagneticField_service.cc.

References reconfigure().

   {
     this->reconfigure(pset);
   }

mag::MagneticField::~MagneticField ( )

inline

Definition at line 50 of file MagneticField.h.

50 {};

Member Function Documentation

G4ThreeVector const MagneticField::FieldAtPoint ( G4ThreeVector const & p = G4ThreeVector(0) ) const

Definition at line 81 of file MagneticField_service.cc.

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::MagneticField::Fields ( ) const

inline

Definition at line 54 of file MagneticField.h.

Referenced by g4b::DetectorConstruction::Construct().

54 { return fFieldDescriptions; }

mag::MagneticField::fFieldDescriptions

std::vector< MagneticFieldDescription > fFieldDescriptions

Descriptions of the fields.

Definition: MagneticField.h:70

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

inline

Definition at line 57 of file MagneticField.h.

References f.

57 { return fFieldDescriptions[f].fVolume; }

mag::MagneticField::fFieldDescriptions

std::vector< MagneticFieldDescription > fFieldDescriptions

Descriptions of the fields.

Definition: MagneticField.h:70

f

TFile f

Definition: plotHisto.C:6

size_t mag::MagneticField::NumFields ( ) const

inline

Definition at line 55 of file MagneticField.h.

55 { return fFieldDescriptions.size(); }

mag::MagneticField::fFieldDescriptions

std::vector< MagneticFieldDescription > fFieldDescriptions

Descriptions of the fields.

Definition: MagneticField.h:70

void MagneticField::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 41 of file MagneticField_service.cc.

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

Referenced by MagneticField().

   {
     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 MagneticField::UniformFieldInVolume ( std::string const & volName ) const

Definition at line 100 of file MagneticField_service.cc.

References DEFINE_ART_SERVICE, and 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);
   }