CCalMagneticField Class Reference

#include <CCalMagneticField.hh>

Inheritance 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 ()

Detailed Description

Definition at line 38 of file CCalMagneticField.hh.

Constructor & Destructor Documentation

CCalMagneticField::CCalMagneticField

(

const G4String &

name

)

Definition at line 42 of file CCalMagneticField.cc.

References findDO(), G4cout, G4endl, python.hepunit::mm, openGeomFile(), and tab().

                                                             :
  fval(0), pos(0), slope(0), intercept(0) {

  //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;
#ifdef debug
    G4cout << "Field value " << fval << " # points " << npts << " offset in x "
     << xoff*mm << G4endl;
#endif

    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];
#ifdef debug
    G4cout << tab << "Position " << i << " " << pos[i] << " Slope "
         << slope[i] << " Intercept " << intercept[i] << G4endl;
#endif
      }
    }

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

CCalMagneticField::~CCalMagneticField

(

)

Definition at line 82 of file CCalMagneticField.cc.

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

Member Function Documentation

G4double CCalMagneticField::GetConstantFieldvalue ( ) const

inline

Definition at line 47 of file CCalMagneticField.hh.

Referenced by CCalDetectorConstruction::Construct().

{return fval;}

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

virtual

Definition at line 141 of file CCalMagneticField.cc.

References MagneticField().

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

G4FieldManager* CCalMagneticField::GetGlobalFieldManager ( )

protected

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

Definition at line 94 of file CCalMagneticField.cc.

References plottest35::c1, G4cout, G4endl, python.hepunit::kilogauss, python.hepunit::m, python.hepunit::mm, and python.hepunit::tesla.

Referenced by GetFieldValue(), and MagneticField().

                                                                          {

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

  G4double m1=0, 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;
#ifdef ddebug
  G4cout << "Field at x: " << x[0]/mm << "mm (" << xnew << ") = " << B[2]/tesla
     << "T (m = " << m1 << ", c = " << c1 << ", scale = " << scor << ")"
     << G4endl;
#endif
}

CLHEP::Hep3Vector CCalMagneticField::MagneticField

(

const CLHEP::Hep3Vector

Point

)

const

Definition at line 125 of file CCalMagneticField.cc.

References MagneticField(), CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), CLHEP::Hep3Vector::setZ(), test::v, test::x, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

                                                 {

  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;
}

---

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

• CCalMagneticField.hh
• CCalMagneticField.cc