evgen::ldm::RayTraceBox Class Reference

RayTraceBox class definiton. More...

Inheritance diagram for evgen::ldm::RayTraceBox:

Public Member Functions
	RayTraceBox (fhicl::ParameterSet const &pset)
	Constructor. More...
	~RayTraceBox ()
	Destructor. More...
void	configure (const fhicl::ParameterSet &) override
	Interface for configuring the particular algorithm tool. More...
bool	IntersectDetector (MeVPrtlFlux &flux, std::array< TVector3, 2 > &intersection, double &weight) override
double	MaxWeight () override
Public Member Functions inherited from evgen::ldm::IRayTrace
virtual	~IRayTrace () noexcept=default
	Virtual Destructor. More...
Public Member Functions inherited from evgen::ldm::IMeVPrtlStage
virtual	~IMeVPrtlStage () noexcept
	Virtual Destructor. More...
	IMeVPrtlStage (const char *name)
TVector3	RandomUnitVector ()
double	GetRandom ()
const char *	Name ()

Private Attributes
geo::BoxBoundedGeo	fBox

Additional Inherited Members
Protected Attributes inherited from evgen::ldm::IMeVPrtlStage
CLHEP::HepRandomEngine *	fEngine
const char *	fName

Detailed Description

RayTraceBox class definiton.

Definition at line 38 of file RayTraceBox_tool.cc.

Constructor & Destructor Documentation

evgen::ldm::RayTraceBox::RayTraceBox

(

fhicl::ParameterSet const &

pset

)

Constructor.

Definition at line 62 of file RayTraceBox_tool.cc.

                                                     :
  IMeVPrtlStage("RayTraceBox") 
{
    this->configure(pset);
}

evgen::ldm::RayTraceBox::~RayTraceBox

(

)

Destructor.

Definition at line 70 of file RayTraceBox_tool.cc.

{
}

Member Function Documentation

void evgen::ldm::RayTraceBox::configure ( const fhicl::ParameterSet &  )

overridevirtual

Interface for configuring the particular algorithm tool.

Parameters

ParameterSet

The input set of parameters for configuration

Implements evgen::ldm::IMeVPrtlStage.

Definition at line 75 of file RayTraceBox_tool.cc.

{
  if (pset.has_key("Box")) {
    std::array<double, 6> box_config = pset.get<std::array<double, 6>>("Box");
    // xmin, xmax, ymin, ymax, zmin, zmax
    fBox = geo::BoxBoundedGeo(box_config[0], box_config[1], box_config[2], box_config[3], box_config[4], box_config[5]);
  }
  else {
    const geo::GeometryCore *geometry = lar::providerFrom<geo::Geometry>();
    fBox = geometry->DetectorEnclosureBox(pset.get<std::string>("Volume"));
  }

  std::cout << "Detector Box." << std::endl;
  std::cout << "X " << fBox.MinX() << " " << fBox.MaxX() << std::endl;
  std::cout << "Y " << fBox.MinY() << " " << fBox.MaxY() << std::endl;
  std::cout << "Z " << fBox.MinZ() << " " << fBox.MaxZ() << std::endl;

}

bool evgen::ldm::RayTraceBox::IntersectDetector ( MeVPrtlFlux &  flux, std::array< TVector3, 2 > &  intersection, double &  weight  )

overridevirtual

Implements evgen::ldm::IRayTrace.

Definition at line 95 of file RayTraceBox_tool.cc.

                                                                                                          {
  std::vector<TVector3> box_intersections = fBox.GetIntersections(flux.pos.Vect(), flux.mom.Vect().Unit());

  if (box_intersections.size() != 2) return false;

  TVector3 A = box_intersections[0];
  TVector3 B = box_intersections[1];

  // make sure that the flux start lies outside the detector
  if ((flux.pos.Vect() - A).Mag() < (A-B).Mag() && (flux.pos.Vect() - B).Mag() < (A-B).Mag()) {
    throw cet::exception("RayTraceBox Exception", "Input portal flux starts inside detector volume: "
        "MeVPrtl start At: (" + std::to_string(flux.pos.X()) + ", " + std::to_string(flux.pos.Y()) + ", " + std::to_string(flux.pos.Z()) + "). "
        "Intersection A At: " + std::to_string(A.X()) + ", " + std::to_string(A.Y()) + ", " + std::to_string(A.Z()) + "). "
        "Intersection B At: " + std::to_string(B.X()) + ", " + std::to_string(B.Y()) + ", " + std::to_string(B.Z()) + ").\n"
    );
  } 

  // if the ray points the wrong way, it doesn't intersect
  if (flux.mom.Vect().Unit().Dot((A - flux.pos.Vect()).Unit()) < 0.) {
    std::cout << "RAYTRACE: MeVPrtl points wrong way" << std::endl;
    return false;
  }

  if ((flux.pos.Vect() - A).Mag() < (flux.pos.Vect() - B).Mag()) {
    intersection = {A, B}; // A is entry, B is exit
  }
  else {
    intersection = {B, A}; // reversed
  }

  weight = 1.;
  return true;
}

double evgen::ldm::RayTraceBox::MaxWeight ( )

inlineoverridevirtual

Implements evgen::ldm::IMeVPrtlStage.

Definition at line 56 of file RayTraceBox_tool.cc.

{ return -1.; }

Member Data Documentation

geo::BoxBoundedGeo evgen::ldm::RayTraceBox::fBox

private

Definition at line 59 of file RayTraceBox_tool.cc.

---

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

• RayTraceBox_tool.cc