sample_burried_wire.C

/*
MyTRIM - a three dimensional binary collision Monte Carlo library.
Copyright (C) 2008-2018  Daniel Schwen <daniel@schwen.de>

This library is free software; you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as
published by the Free Software Foundation; either version 2.1 of the
License, or (at your option) any later version.

This library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public
License along with this library; if not, write to the Free Software
Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
02110-1301 USA
*/

#include <cmath>
#if !defined(__APPLE__)
#include "malloc.h"
#endif

#include "sample_burried_wire.h"

using namespace MyTRIM_NS;

SampleBurriedWire::SampleBurriedWire(Real x, Real y, Real z) : SampleWire(x, y, z)
{
  bc[0] = INF;
  bc[1] = INF;
  bc[2] = INF;
}

// look if we are within dr of the wire axis
MaterialBase *
SampleBurriedWire::lookupMaterial(Point & pos)
{
  // cover layer
  if (pos(2) < 0.0 && pos(2) >= -250.0)
    return material[1];

  // above sample or inside substrate
  if (pos(2) > w[2] || pos(2) < -250.0)
    return 0;

  // in wire layer
  MaterialBase * ret = SampleWire::lookupMaterial(pos);
  if (ret == 0)
    return material[1];
  else
    return ret;
}