#include <Field.h>

Inheritance diagram for sim::Field:

Public Member Functions
	Field (const MagneticField *f, double d)

void	GetFieldValue (const G4double p[4], G4double b[3]) const override

	~Field () override

Private Attributes
double	oldb [3]

double	oldx [3]

const MagneticField *	theCMSMagneticField

double	theDelta

Detailed Description

Definition at line 9 of file Field.h.

Constructor & Destructor Documentation

Field::Field	(	const MagneticField *	f,
		double	d
	)

Definition at line 11 of file Field.cc.

References mps_fire::i, oldb, and oldx.

                                              : G4MagneticField(), theCMSMagneticField(f), theDelta(d) {
   for (int i = 0; i < 3; ++i) {
     oldx[i] = 1.0e12;
     oldb[i] = 0.0;
   }
 }

Field::~Field ( )

override

Definition at line 18 of file Field.cc.

18 {}

Member Function Documentation

void Field::GetFieldValue	(	const G4double	p[4],
		G4double	b[3]
	)		const

override

Definition at line 20 of file Field.cc.

References funct::abs(), objects.autophobj::float, MagneticField::inTesla(), oldb, oldx, theCMSMagneticField, theDelta, findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                          {
   if (std::abs(oldx[0] - xyz[0]) > theDelta || std::abs(oldx[1] - xyz[1]) > theDelta ||
       std::abs(oldx[2] - xyz[2]) > theDelta) {
     static const float lunit = (float)(1.0 / CLHEP::cm);
     GlobalPoint ggg((float)(xyz[0]) * lunit, (float)(xyz[1]) * lunit, (float)(xyz[2]) * lunit);
     GlobalVector v = theCMSMagneticField->inTesla(ggg);
 
     static const float btesla = (float)CLHEP::tesla;
     oldb[0] = (G4double)(v.x() * btesla);
     oldb[1] = (G4double)(v.y() * btesla);
     oldb[2] = (G4double)(v.z() * btesla);
     oldx[0] = xyz[0];
     oldx[1] = xyz[1];
     oldx[2] = xyz[2];
   }
 
   bfield[0] = oldb[0];
   bfield[1] = oldb[1];
   bfield[2] = oldb[2];
 }