Heed::HeedFieldMap Class Reference

Retrieve electric and magnetic field from Sensor. More...

#include <HeedFieldMap.h>

Public Member Functions
	HeedFieldMap ()=default
void	SetSensor (Garfield::Sensor *sensor)
void	SetCentre (const double x, const double y, const double z)
void	UseEfield (const bool flag)
void	UseBfield (const bool flag)
void	field_map (const point &pt, vec &efield, vec &bfield, vfloat &mrange) const
bool	inside (const point &pt)

Detailed Description

Retrieve electric and magnetic field from Sensor.

Definition at line 15 of file HeedFieldMap.h.

Constructor & Destructor Documentation

HeedFieldMap()

Heed::HeedFieldMap::HeedFieldMap ( )

default

Member Function Documentation

field_map()

void Heed::HeedFieldMap::field_map	(	const point &	pt,
		vec &	efield,
		vec &	bfield,
		vfloat &	mrange
	)		const

Definition at line 11 of file HeedFieldMap.cpp.

                                                   {
 
  const double x = pt.v.x * conv + m_x;
  const double y = pt.v.y * conv + m_y;
  const double z = pt.v.z * conv + m_z;
 
  // Initialise the electric and magnetic field.
  efield.x = bfield.x = 0.;
  efield.y = bfield.y = 0.;
  efield.z = bfield.z = 0.;
  mrange = DBL_MAX;
 
  if (!m_sensor) {
    std::cerr << "HeedFieldMap::field_map: Sensor not defined.\n";
    return;
  }
 
  if (m_useEfield) {
    double ex = 0., ey = 0., ez = 0.;
    int status = 0;
    Garfield::Medium* m = nullptr;
    m_sensor->ElectricField(x, y, z, ex, ey, ez, m, status);
    efield.x = ex * 1.e-7;
    efield.y = ey * 1.e-7;
    efield.z = ez * 1.e-7;
  }
 
  if (m_useBfield) {
    double bx = 0., by = 0., bz = 0.;
    int status = 0;
    m_sensor->MagneticField(x, y, z, bx, by, bz, status);
    bfield.x = bx * 1.e-3;
    bfield.y = by * 1.e-3;
    bfield.z = bz * 1.e-3;
  }
}

Referenced by Heed::eparticle::force().

inside()

bool Heed::HeedFieldMap::inside

(

const point &

pt

)

Definition at line 49 of file HeedFieldMap.cpp.

                                         {
 
  const double x = pt.v.x * conv + m_x;
  const double y = pt.v.y * conv + m_y;
  const double z = pt.v.z * conv + m_z;
  // Check if the point is inside the drift area.
  if (!m_sensor->IsInArea(x, y, z)) return false;
  // Check if the point is inside a medium.
  Garfield::Medium* m = nullptr;
  if (!m_sensor->GetMedium(x, y, z, m) || !m) return false;
  return m->IsIonisable();
}

Referenced by Heed::HeedDeltaElectron::physics_after_new_speed().

SetCentre()

void Heed::HeedFieldMap::SetCentre ( const double  x, const double  y, const double  z  )

inline

Definition at line 20 of file HeedFieldMap.h.

                                                                 {
    m_x = x;
    m_y = y;
    m_z = z;
  }

SetSensor()

void Heed::HeedFieldMap::SetSensor ( Garfield::Sensor *  sensor )

inline

Definition at line 19 of file HeedFieldMap.h.

{ m_sensor = sensor; }

UseBfield()

void Heed::HeedFieldMap::UseBfield ( const bool  flag )

inline

Definition at line 26 of file HeedFieldMap.h.

{ m_useBfield = flag; }

UseEfield()

void Heed::HeedFieldMap::UseEfield ( const bool  flag )

inline

Definition at line 25 of file HeedFieldMap.h.

{ m_useEfield = flag; }

---

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

• HeedFieldMap.h
• HeedFieldMap.cpp