#include <box.h>

Inheritance diagram for Heed::box:

Public Member Functions
	box (void)
	Constructors.

	box (vfloat fdx, vfloat fdy, vfloat fdz, const String &fname)

	box (vfloat fdx, vfloat fdy, vfloat fdz, vfloat fprec, const String &fname)

	box (box &fb)

	box (const box &fb)

	macro_copy_header (box)

virtual	~box ()
	Destructor.

void	init_prec (void)

void	init_planes (void)

virtual int	check_point_inside (const point &fpt, const vec &dir) const

virtual int	range_ext (trajestep &fts, int s_ext) const

virtual void	income (gparticle *gp)

virtual void	chname (char *nm) const

virtual void	print (std::ostream &file, int l) const

Public Member Functions inherited from absvol
virtual	~absvol ()

virtual int	check_point_inside (const point &fpt, const vec &dir) const =0

virtual int	find_embed_vol (const point &fpt, const vec &dir, manip_absvol_treeid *atid) const

virtual int	range (trajestep &fts, int s_ext, int &sb, manip_absvol_eid *faeid) const

virtual int	range_ext (trajestep &fts, int s_ext) const =0

	macro_copy_header (absvol)

virtual void	income (gparticle *)

virtual void	chname (char *nm) const

virtual void	print (std::ostream &file, int l) const

virtual DynLinArr< manip_absvol * >	Gamanip_embed (void) const

Public Member Functions inherited from absref
virtual	~absref ()

virtual void	down (const abssyscoor *fasc)

virtual void	up (const abssyscoor *fasc)

virtual void	turn (const vec &dir, vfloat angle)

virtual void	shift (const vec &dir)

Public Attributes
vfloat	dx

vfloat	dy

vfloat	dz

vfloat	dxh

vfloat	dyh

vfloat	dzh

ulsvolume	ulsv

String	name

Public Attributes inherited from absvol
vfloat	prec

Protected Member Functions
virtual void	get_components (ActivePtr< absref_transmit > &aref_tran)

Detailed Description

Definition at line 24 of file box.h.

Constructor & Destructor Documentation

◆ box() [1/5]

Heed::box::box ( void )

Constructors.

Definition at line 26 of file box.cpp.

             : dx(0), dy(0), dz(0), dxh(0), dyh(0), dzh(0), name("none") {
  mfunname("box::box(void)");
  init_prec();
  init_planes();
}

◆ box() [2/5]

Heed::box::box	(	vfloat	fdx,
		vfloat	fdy,
		vfloat	fdz,
		const String &	fname
	)

Definition at line 32 of file box.cpp.

                                                                {
  pvecerror("box(vfloat fdx, vfloat fdy, vfloat fdz, const String &fname)");
  dx = abslt(fdx);
  dy = abslt(fdy);
  dz = abslt(fdz);
  dxh = 0.5 * dx;
  dyh = 0.5 * dy;
  dzh = 0.5 * dz;
  name = fname;
  init_prec();
  init_planes();
}

◆ box() [3/5]

Heed::box::box	(	vfloat	fdx,
		vfloat	fdy,
		vfloat	fdz,
		vfloat	fprec,
		const String &	fname
	)

Definition at line 45 of file box.cpp.

                              {
  pvecerror("box(vfloat fdx, vfloat fdy, vfloat fdz, const String &fname)");
  dx = abslt(fdx);
  dy = abslt(fdy);
  dz = abslt(fdz);
  dxh = 0.5 * dx;
  dyh = 0.5 * dy;
  dzh = 0.5 * dz;
  name = fname;
  prec = fprec;
  init_planes();
}

◆ box() [4/5]

Heed::box::box ( box & fb )

Definition at line 59 of file box.cpp.

                : absref(fb), absvol(fb) {
  pvecerror("box(const box& fb)");
  dx = fb.dx;
  dy = fb.dy;
  dz = fb.dz;
  dxh = 0.5 * dx;
  dyh = 0.5 * dy;
  dzh = 0.5 * dz;
  prec = fb.prec;
  name = fb.name;
  init_planes();
}

◆ box() [5/5]

Heed::box::box ( const box & fb )

Definition at line 72 of file box.cpp.

                      : absref(fb), absvol(fb) {
  pvecerror("box(const box& fb)");
  dx = fb.dx;
  dy = fb.dy;
  dz = fb.dz;
  dxh = 0.5 * dx;
  dyh = 0.5 * dy;
  dzh = 0.5 * dz;
  name = fb.name;
  prec = fb.prec;
  init_planes();
}

◆ ~box()

virtual Heed::box::~box ( )

inlinevirtual

Destructor.

Definition at line 43 of file box.h.

43{}

Member Function Documentation

◆ check_point_inside()

int Heed::box::check_point_inside	(	const point &	fpt,
		const vec &	dir
	)		const

virtual

Implements absvol.

Definition at line 105 of file box.cpp.

                                                                  {
  mfunname("virtual int check_point_inside(const point& fpt, const vec& dir)");
#ifdef TRACE_find_embed_vol
  mcout << "box::check_point_inside: \n";
  print(mcout, 1);
  mcout << "fpt=" << fpt << "dir=" << dir;
#endif
  if (dir == dv0) {
    if (abslt(fpt.v.x) <= dxh && abslt(fpt.v.y) <= dyh && 
        abslt(fpt.v.z) <= dzh) {
      return 1;
    }
    return 0;
  } else {
    if (abslt(fpt.v.x) <= dxh - prec && abslt(fpt.v.y) <= dyh - prec &&
        abslt(fpt.v.z) <= dzh - prec) {
#ifdef TRACE_find_embed_vol
      mcout << "cond 1, returning 1\n";
#endif
      return 1;
    }
    if (abslt(fpt.v.x) > dxh + prec || abslt(fpt.v.y) > dyh + prec ||
        abslt(fpt.v.z) > dzh + prec) {
#ifdef TRACE_find_embed_vol
      if (abslt(fpt.v.x) > dxh + prec) mcout << "cond 2.1 satisfied\n";
      if (abslt(fpt.v.y) > dyh + prec) mcout << "cond 2.2 satisfied\n";
      if (abslt(fpt.v.z) > dzh + prec) mcout << "cond 2.3 satisfied\n";
      mcout << "cond 2, returning 0\n";
#endif
      return 0;
    }
    // What remains is point belonging to border
#ifdef IMPROVED_BOUNDARY
    // Below we detect cases when particle is exiting, leaving the
    // case when it is entering
    if (abslt(fpt.v.x) > dxh - prec) {
      if (dir.x == 0.0) return 0;
      if ((fpt.v.x > 0 && dir.x > 0) || (fpt.v.x < 0 && dir.x < 0)) {
#ifdef TRACE_find_embed_vol
        mcout << "cond 3, returning 0\n";
#endif
        return 0;
      }
    }
    if (abslt(fpt.v.y) > dyh - prec) {
      if (dir.y == 0.0) return 0;
      if ((fpt.v.y > 0 && dir.y > 0) || (fpt.v.y < 0 && dir.y < 0)) {
#ifdef TRACE_find_embed_vol
        mcout << "cond 4, returning 0\n";
#endif
        return 0;
      }
    }
    if (abslt(fpt.v.z) > dzh - prec) {
      if (dir.z == 0.0) return 0;
      if ((fpt.v.z > 0 && dir.z > 0) || (fpt.v.z < 0 && dir.z < 0)) {
#ifdef TRACE_find_embed_vol
        mcout << "cond 5, returning 0\n";
#endif
        return 0;
      }
    }
#ifdef TRACE_find_embed_vol
    mcout << "finish, returning 1\n";
#endif
    return 1;
 
#else
    // for IMPROVED_BOUNDARY
    // In the old version, which is below,
    // if the track is parallel to a boundary, it was interpreted as
    // signature of being inside volume.
    // The general principal, saying that
    // if particle is located at the boundary and headed
    // outside the volume,
    // it is considered in the volume to which it is headed,
    // cannot help to choose between two versions.
    // Of course, the box is convex and the particle flying along the border
    // enters the external space sooner or later and it can be regarded as
    // already outside. But the same can be said about the particle,
    // which is completely indise the volume.
    // So this principle does not work here.
    // The other principle should be applied.
    // If we allow the absolutely thing volumes (control surfaces)
    // we should use the algorithm, which, in particular stops the particle
    // crossing exactly the corner of volume without entering its inside.
    // But this does not allow to choose. So now the old (this) variant
    // is used, untill other arguments appear.
 
    if (abslt(fpt.v.x) > dxh - prec &&
        ((fpt.v.x > 0 && dir.x > 0) || (fpt.v.x < 0 && dir.x < 0))) {
#ifdef TRACE_find_embed_vol
      mcout << "cond 3, returning 0\n";
#endif
      return 0;  // exiting
    }
    if (abslt(fpt.v.y) > dyh - prec &&
        ((fpt.v.y > 0 && dir.y > 0) || (fpt.v.y < 0 && dir.y < 0))) {
#ifdef TRACE_find_embed_vol
      mcout << "cond 4, returning 0\n";
#endif
      return 0;
    }
    if (abslt(fpt.v.z) > dzh - prec &&
        ((fpt.v.z > 0 && dir.z > 0) || (fpt.v.z < 0 && dir.z < 0))) {
#ifdef TRACE_find_embed_vol
      mcout << "cond 5, returning 0\n";
#endif
      return 0;
    }
#ifdef TRACE_find_embed_vol
    mcout << "finish, returning 1\n";
#endif
    return 1;
#endif
  }
}

◆ chname()

void Heed::box::chname ( char * nm ) const

virtual

Reimplemented from absvol.

Reimplemented in Heed::manip_box, and Heed::sh_manip_box.

Definition at line 260 of file box.cpp.

                               {
#ifdef USE_STLSTRING
  strcpy(nm, "box: ");
  strcat(nm, name.c_str());
#else
  strcpy(nm, "box: ");
  strcat(nm, name);
#endif
}

Referenced by print().

◆ get_components()

void Heed::box::get_components ( ActivePtr< absref_transmit > & aref_tran )

protectedvirtual

Reimplemented from absref.

Reimplemented in Garfield::HeedChamber, and Heed::sh_manip_box.

Definition at line 18 of file box.cpp.

                                                      {
  mfunnamep("box::get_components(...)");
  funnw.ehdr(mcerr);
  mcerr << "one should not call this function, since this object cannot be "
           "modified\n";
  spexit(mcerr);
}

◆ income()

virtual void Heed::box::income ( gparticle * gp )

virtual

Reimplemented from absvol.

◆ init_planes()

void Heed::box::init_planes ( void )

Definition at line 90 of file box.cpp.

                          {
  mfunname("void box::init_planes(void)");
  splane spl[6];
  spl[0] = splane(plane(point(dxh, 0, 0), vec(-1, 0, 0)), vec(-1, 0, 0));
  spl[1] = splane(plane(point(-dxh, 0, 0), vec(1, 0, 0)), vec(1, 0, 0));
  spl[2] = splane(plane(point(0, dyh, 0), vec(0, -1, 0)), vec(0, -1, 0));
  spl[3] = splane(plane(point(0, -dyh, 0), vec(0, 1, 0)), vec(0, 1, 0));
  spl[4] = splane(plane(point(0, 0, dzh), vec(0, 0, -1)), vec(0, 0, -1));
  spl[5] = splane(plane(point(0, 0, -dzh), vec(0, 0, 1)), vec(0, 0, 1));
  surface* fsurf[6];
  for (int n = 0; n < 6; ++n)
    fsurf[n] = &spl[n];
  ulsv.ulsvolume_init(fsurf, 6, "ulsv of box", prec);
}

Referenced by box().

◆ init_prec()

void Heed::box::init_prec ( void )

Definition at line 85 of file box.cpp.

                        {
  prec = (dxh + dyh + dzh) / 3.0;
  prec *= vprecision;
}

Referenced by box().

◆ macro_copy_header()

Heed::box::macro_copy_header ( box )

◆ print()

void Heed::box::print	(	std::ostream &	file,
		int	l
	)		const

virtual

Reimplemented from absvol.

Reimplemented in Heed::manip_box, and Heed::sh_manip_box.

Definition at line 223 of file box.cpp.

                                             {
  if (l <= 0) return;
  char s[1000];
  chname(s);
  Ifile << "box::print(l=" << l << "): " << s << '\n';
  indn.n += 2;
  Ifile << " dx=" << dx << " dy=" << dy << " dz=" << dz << " prec=" << prec
        << '\n';
  Ifile << " dxh=" << dxh << " dyh=" << dyh << " dzh=" << dzh << '\n';
  if (l >= 10) {
    l--;
    indn.n += 2;
    ulsv.print(file, l);
    indn.n -= 2;
  }
  absvol::print(file, l);
  indn.n -= 2;
}

Referenced by check_point_inside().

◆ range_ext()

int Heed::box::range_ext	(	trajestep &	fts,
		int	s_ext
	)		const

virtual

Implements absvol.

Definition at line 242 of file box.cpp.

                                                  {
  mfunname("virtual int box::range_ext(trajestep& fts, int s_ext) const");
  if (s_ext == 0) {
    if (abslt(fts.currpos.v.x) > dxh + fts.mrange) return 0;
    if (abslt(fts.currpos.v.y) > dyh + fts.mrange) return 0;
    if (abslt(fts.currpos.v.z) > dzh + fts.mrange) return 0;
  } else {
    if (abslt(fts.currpos.v.x) < dxh - fts.mrange &&
        abslt(fts.currpos.v.y) < dyh - fts.mrange &&
        abslt(fts.currpos.v.z) < dzh - fts.mrange) {
      return 0;
    }
  }
  return ulsv.range_ext(fts, s_ext);
}

Member Data Documentation

◆ dx

vfloat Heed::box::dx

Definition at line 26 of file box.h.

Referenced by box(), and print().

◆ dxh

vfloat Heed::box::dxh

Definition at line 27 of file box.h.

Referenced by box(), check_point_inside(), init_planes(), init_prec(), print(), and range_ext().

◆ dy

vfloat Heed::box::dy

Definition at line 26 of file box.h.

Referenced by box(), and print().

◆ dyh

vfloat Heed::box::dyh

Definition at line 27 of file box.h.

Referenced by box(), check_point_inside(), init_planes(), init_prec(), print(), and range_ext().

◆ dz

vfloat Heed::box::dz

Definition at line 26 of file box.h.

Referenced by box(), and print().

◆ dzh

vfloat Heed::box::dzh

Definition at line 27 of file box.h.

Referenced by box(), check_point_inside(), init_planes(), init_prec(), print(), and range_ext().

◆ name

String Heed::box::name

Definition at line 29 of file box.h.

Referenced by box(), chname(), and Heed::macro_copy_body().

◆ ulsv

ulsvolume Heed::box::ulsv

Definition at line 28 of file box.h.

Referenced by init_planes(), print(), and range_ext().

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

Public Member Functions

Public Attributes

Protected Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ box() [1/5]

◆ box() [2/5]

◆ box() [3/5]

◆ box() [4/5]

◆ box() [5/5]

◆ ~box()

Member Function Documentation

◆ check_point_inside()

◆ chname()

◆ get_components()

◆ income()

◆ init_planes()

◆ init_prec()

◆ macro_copy_header()

◆ print()

◆ range_ext()

Member Data Documentation

◆ dx

◆ dxh

◆ dy

◆ dyh

◆ dz

◆ dzh

◆ name

◆ ulsv