#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 20 of file Field.cc.

20 {}

Member Function Documentation

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

override

Definition at line 22 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];
 }