CCalMagneticField Class Reference

#include <CCalMagneticField.hh>

Inheritance diagram for CCalMagneticField:

Collaboration diagram for CCalMagneticField:

Public Member Functions
	CCalMagneticField (const G4String &name)
	~CCalMagneticField ()
void	MagneticField (const double Point[3], double Bfield[3]) const
CLHEP::Hep3Vector	MagneticField (const CLHEP::Hep3Vector Point) const
virtual void	GetFieldValue (const double Point[3], double *Bfield) const
G4double	GetConstantFieldvalue () const
Public Member Functions inherited from G4MagneticField
	G4MagneticField ()
virtual	~G4MagneticField ()
	G4MagneticField (const G4MagneticField &r)
G4MagneticField &	operator= (const G4MagneticField &p)
G4bool	DoesFieldChangeEnergy () const
virtual void	GetFieldValue (const G4double Point[4], G4double *Bfield) const =0
Public Member Functions inherited from G4ElectroMagneticField
	G4ElectroMagneticField ()
virtual	~G4ElectroMagneticField ()
	G4ElectroMagneticField (const G4ElectroMagneticField &r)
G4ElectroMagneticField &	operator= (const G4ElectroMagneticField &p)
Public Member Functions inherited from G4Field
	G4Field (G4bool gravityOn=false)
	G4Field (const G4Field &)
virtual	~G4Field ()
G4Field &	operator= (const G4Field &p)
G4bool	IsGravityActive () const
void	SetGravityActive (G4bool OnOffFlag)
virtual G4Field *	Clone () const

Protected Member Functions
G4FieldManager *	GetGlobalFieldManager ()

Private Attributes
G4double	fval
G4int	npts
G4double	xoff
G4double *	pos
G4double *	slope
G4double *	intercept
G4int	fVerbosity

Detailed Description

Definition at line 38 of file CCalMagneticField.hh.

Constructor & Destructor Documentation

CCalMagneticField()

CCalMagneticField::CCalMagneticField

(

const G4String &

name

)

Definition at line 42 of file CCalMagneticField.cc.

                                                             :
  fval(0), pos(0), slope(0), intercept(0) 
{
#ifdef debug
  fVerbosity = 1;
#else
  fVerbosity = 0;
#endif

  //Let's open the file
  G4cout << " ==> Opening file " << filename << " to read magnetic field..."
     << G4endl;
  G4String pathName = getenv("CCAL_GLOBALPATH");
  std::ifstream is;
  bool ok = openGeomFile(is, pathName, filename);
  
  if (ok) {
    findDO(is, G4String("FLDM"));
    is >> fval >> npts >> xoff;

    if (fVerbosity)
      G4cout << "Field value " << fval << " # points " << npts << 
    " offset in x "
         << xoff*mm << G4endl;

    if (npts > 0) {
      pos       = new G4double[npts];
      slope     = new G4double[npts];
      intercept = new G4double[npts];

      for (G4int i = 0; i < npts; i++) {
    is >> pos[i] >> slope[i] >> intercept[i];
    if (fVerbosity)
      G4cout << tab << "Position " << i << " " << pos[i] << " Slope "
         << slope[i] << " Intercept " << intercept[i] << G4endl;
      }
    }

    // Close the file
    G4cout << " ==> Closing file " << filename << G4endl;
    is.close();
  }
}

Here is the call graph for this function:

~CCalMagneticField()

CCalMagneticField::~CCalMagneticField

(

)

Definition at line 88 of file CCalMagneticField.cc.

                                      {
  if (pos)
    delete[] pos;
  if (slope)
    delete[] slope;
  if (intercept)
    delete[] intercept;
}

Member Function Documentation

GetConstantFieldvalue()

G4double CCalMagneticField::GetConstantFieldvalue ( ) const

inline

Definition at line 47 of file CCalMagneticField.hh.

{return fval;}

Here is the call graph for this function:

Here is the caller graph for this function:

GetFieldValue()

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

virtual

Definition at line 152 of file CCalMagneticField.cc.

                                                                        {
  CCalMagneticField::MagneticField(x, B);
}

Here is the call graph for this function:

GetGlobalFieldManager()

G4FieldManager* CCalMagneticField::GetGlobalFieldManager ( )

protected

Here is the caller graph for this function:

MagneticField() [1/2]

void CCalMagneticField::MagneticField	(	const double	Point[3],
		double	Bfield[3]
	)		const

Definition at line 100 of file CCalMagneticField.cc.

{
  G4int i=0;
  for (i=0; i<2; i++) {
    B[i]   = 0*kilogauss;
  }

  G4double m1=0;
  G4double c1=0;
  G4double xnew = x[0]/mm + xoff;
  if (npts > 0) {
    for (i=0; i<npts; i++) {
      if (xnew > pos[i]*mm) {
    m1 = slope[i];
        c1 = intercept[i];
      }
    }
  }
  G4double scor = c1 + m*xnew;
  if (scor < 0.) scor = 0.;
  if (scor > 1.) scor = 1.0;

  B[2] = scor*fval*kilogauss;
  if (fVerbosity)
    {

      G4cout << "Field at x: " << x[0]/mm << "mm (" << xnew << ") = " << 
    B[2]/tesla
         << "T (m = " << m1 << ", c = " << 
    c1 << ", scale = " << scor << ")"
         << G4endl;
    }
}

Here is the caller graph for this function:

MagneticField() [2/2]

CLHEP::Hep3Vector CCalMagneticField::MagneticField

(

const CLHEP::Hep3Vector

Point

)

const

Definition at line 136 of file CCalMagneticField.cc.

                                                 {

  G4double x[3],B[3];
  CLHEP::Hep3Vector v;
  
  x[0] = point.x();
  x[1] = point.y();
  x[2] = point.z();
  CCalMagneticField::MagneticField(x, B);
  v.setX(B[0]);   
  v.setY(B[1]);   
  v.setZ(B[2]);   
  return v;
}

Here is the call graph for this function:

Member Data Documentation

fval

G4double CCalMagneticField::fval

private

Definition at line 54 of file CCalMagneticField.hh.

fVerbosity

G4int CCalMagneticField::fVerbosity

private

Definition at line 61 of file CCalMagneticField.hh.

intercept

G4double* CCalMagneticField::intercept

private

Definition at line 59 of file CCalMagneticField.hh.

npts

G4int CCalMagneticField::npts

private

Definition at line 55 of file CCalMagneticField.hh.

pos

G4double* CCalMagneticField::pos

private

Definition at line 57 of file CCalMagneticField.hh.

slope

G4double* CCalMagneticField::slope

private

Definition at line 58 of file CCalMagneticField.hh.

xoff

G4double CCalMagneticField::xoff

private

Definition at line 56 of file CCalMagneticField.hh.

---

The documentation for this class was generated from the following files:

• CCalMagneticField.hh
• CCalMagneticField.cc