#include <LBNEMagneticField.hh>

Inheritance diagram for LBNEMagneticFieldDecayPipe:

Public Member Functions
	LBNEMagneticFieldDecayPipe (bool isToroidal)

virtual void	GetFieldValue (const double Point[3], double *Bfield) const

void	SetWireCurrent (G4double iC)

G4double	GetWireCurrent () const

void	SetWireRadius (G4double r)

G4double	GetWireRadius () const

Private Attributes
G4bool	isToroid

G4double	fWireCurrent

G4double	fWireRadius

Detailed Description

Definition at line 132 of file LBNEMagneticField.hh.

Constructor & Destructor Documentation

LBNEMagneticFieldDecayPipe::LBNEMagneticFieldDecayPipe ( bool isToroidal )

explicit

Definition at line 479 of file LBNEMagneticField.cc.

                                                                       : 
 isToroid(isToroidal),
 fWireCurrent(0.),
 fWireRadius(0.)
 {
 }

Member Function Documentation

void LBNEMagneticFieldDecayPipe::GetFieldValue	(	const double	Point[3],
		double *	Bfield
	)		const

virtual

Definition at line 485 of file LBNEMagneticField.cc.

 {
 
   
    G4double current = fWireCurrent/CLHEP::ampere/1000.;
    for (size_t k=0; k!=3; ++k) Bfield[k] = 0.;
 //
 // No misalignement of this wire..
 //
    const double r = std::sqrt(Point[0]*Point[0] + Point[1]*Point[1]); 
    if (r < fWireRadius) return;
    
    double magBField = current / (5.*r/CLHEP::cm)/10*CLHEP::tesla; //B(kG)=i(kA)/[5*r(cm)], 1T=10kG
    Bfield[0] = -magBField*Point[1]/r;
    Bfield[1] = magBField*Point[0]/r;
    
 //   std::cerr << " LBNEMagneticFieldWire::GetFieldValue Z " << Point[2] 
 //             << " current " << current << " r " << r << " magfield " << magBField << std::endl;
    
 }