genie::geometry::GeomVolSelectorRockBox Class Reference

GENIE Interface for limiting vertex selection in the rock to a volume that depends (in part) on the neutrino p4. Uses GeomVolSelectorFiducial to possibly exclude an inner region. More...

#include <GeomVolSelectorRockBox.h>

Inheritance diagram for genie::geometry::GeomVolSelectorRockBox:

Public Member Functions
	GeomVolSelectorRockBox ()
virtual	~GeomVolSelectorRockBox ()
void	TrimSegment (PathSegment &segment) const
void	BeginPSList (const PathSegmentList *untrimmed) const
void	EndPSList () const
void	SetRockBoxMinimal (Double_t *xyzmin, Double_t *xyzmax)
void	SetRockBoxInclusion (Double_t *xyzmin, Double_t *xyzmax)
void	SetMinimumWall (Double_t w)
void	SetDeDx (Double_t dedx)
void	SetExpandFromInclusion (bool how=false)
virtual void	ConvertShapeMaster2Top (const ROOTGeomAnalyzer *rgeom)
Public Member Functions inherited from genie::geometry::GeomVolSelectorFiducial
	GeomVolSelectorFiducial ()
virtual	~GeomVolSelectorFiducial ()
void	SetReverseFiducial (Bool_t reverse=true)
void	AdoptFidShape (FidShape *shape)
void	MakeSphere (Double_t x0, Double_t y0, Double_t z0, Double_t radius)
void	MakeXCylinder (Double_t y0, Double_t z0, Double_t radius, Double_t xmin, Double_t xmax)
void	MakeYCylinder (Double_t x0, Double_t z0, Double_t radius, Double_t ymin, Double_t ymax)
void	MakeZCylinder (Double_t x0, Double_t y0, Double_t radius, Double_t zmin, Double_t zmax)
void	MakeCylinder (Double_t *base, Double_t *axis, Double_t radius, Double_t *cap1, Double_t *cap2)
void	MakeBox (Double_t *xyzmin, Double_t *xyzmax)
void	MakeZPolygon (Int_t n, Double_t x0, Double_t y0, Double_t inradius, Double_t phi0deg, Double_t zmin, Double_t zmax)
Public Member Functions inherited from genie::geometry::GeomVolSelectorBasic
	GeomVolSelectorBasic ()
virtual	~GeomVolSelectorBasic ()
void	SetVolumeSelection (string volstr)
void	SetMediumSelection (string medstr)
void	SetMaterialSelection (string matstr)
void	SetPathSelection (string pathstr)
Public Member Functions inherited from genie::geometry::GeomVolSelectorI
virtual	~GeomVolSelectorI ()
virtual PathSegmentList *	GenerateTrimmedList (const PathSegmentList *untrimmed) const
void	SetCurrentRay (const TLorentzVector &x4, const TLorentzVector &p4)
	configure for individual neutrino ray More...
void	SetSI2Local (double scale)
	set scale factor for SI to "raydist" units of PathSegmentList More...
void	SetRemoveEntries (bool rmset)
bool	GetRemoveEntries ()
void	SetNeedPath ()
bool	GetNeedPath () const
	allow toggle on only More...
std::string	GetName () const

Protected Member Functions
void	MakeRockBox () const
Protected Member Functions inherited from genie::geometry::GeomVolSelectorBasic
void	ParseSelection (const string &str, vector< string > &required, vector< string > &forbidden)
bool	RejectString (const string &str, const vector< string > &required, const vector< string > &forbidden) const
Protected Member Functions inherited from genie::geometry::GeomVolSelectorI
	GeomVolSelectorI ()
	GeomVolSelectorI (std::string name)

Protected Attributes
Double_t	fMinimalXYZMin [3]
Double_t	fMinimalXYZMax [3]
	interior box lower corner More...
Double_t	fMinimumWall
	interior box upper corner More...
Double_t	fInclusionXYZMin [3]
	minimum distance around (XYZmin,XYZmax) More...
Double_t	fInclusionXYZMax [3]
	box within which events are always More...
Double_t	fDeDx
	accepted More...
Bool_t	fExpandInclusion
	how to scale from energy to distance More...
FidShape *	fRockBoxShape
	expand from minimal or inclusion box? More...
const ROOTGeomAnalyzer *	fROOTGeom
	shape changes for every nu ray More...
RayIntercept	fInterceptRock
Protected Attributes inherited from genie::geometry::GeomVolSelectorFiducial
Bool_t	fSelectReverse
FidShape *	fShape
	select for "outside" fiducial? More...
const PathSegmentList *	fCurrPathSegmentList
	shape More...
RayIntercept	fIntercept
Protected Attributes inherited from genie::geometry::GeomVolSelectorBasic
vector< string >	fRequiredVol
vector< string >	fRequiredMed
vector< string >	fRequiredMat
vector< string >	fRequiredPath
vector< string >	fForbiddenVol
vector< string >	fForbiddenMed
vector< string >	fForbiddenMat
vector< string >	fForbiddenPath
Protected Attributes inherited from genie::geometry::GeomVolSelectorI
TLorentzVector	fX4
	current neutrino ray's start position (global) More...
TLorentzVector	fP4
	current neutrino ray's momentum (global) More...
double	fScale
	SI->raydist scale factor. More...
bool	fRemoveEntries
	whether selector should remove entries or set hi=lo More...
bool	fNeedPath
	selector needs PathSegment "path" string More...
std::string	fName
	volume selector name More...

Additional Inherited Members
Static Protected Member Functions inherited from genie::geometry::GeomVolSelectorFiducial
static Bool_t	NewStepPairs (Bool_t selectReverse, Double_t raydist, Double_t slo, Double_t shi, const RayIntercept &intercept, Bool_t &split, StepRange &step1, StepRange &step2)

Detailed Description

GENIE Interface for limiting vertex selection in the rock to a volume that depends (in part) on the neutrino p4. Uses GeomVolSelectorFiducial to possibly exclude an inner region.

Author

Robert Hatcher rhatc.nosp@m.her@.nosp@m.fnal..nosp@m.gov FNAL

August 5, 2010

Copyright (c) 2003-2020, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 37 of file GeomVolSelectorRockBox.h.

Constructor & Destructor Documentation

GeomVolSelectorRockBox::GeomVolSelectorRockBox

(

)

Definition at line 22 of file GeomVolSelectorRockBox.cxx.

  : GeomVolSelectorFiducial(), fMinimumWall(0.), fDeDx(1.)
  , fExpandInclusion(false)
  , fRockBoxShape(0), fROOTGeom(0)
{
  fName = "RockBox";
  // base class' fiducial volume always treated as reverse (if even exists)
  this->SetReverseFiducial(true);

  for (int i=0; i<3; ++i) {
    fMinimalXYZMin[i] = 0;
    fMinimalXYZMax[i] = 0;
    fInclusionXYZMin[i] = 0;
    fInclusionXYZMax[i] = 0;
  }
}

GeomVolSelectorRockBox::~GeomVolSelectorRockBox ( )

virtual

Definition at line 40 of file GeomVolSelectorRockBox.cxx.

{
  if ( fRockBoxShape ) delete fRockBoxShape;
  fRockBoxShape = 0;
  fROOTGeom = 0;  // was reference only
}

Member Function Documentation

void GeomVolSelectorRockBox::BeginPSList ( const PathSegmentList *  untrimmed ) const

virtual

Every derived version must also respond to a signal that starts a new path segment list processing and ends it. In general they can simply ignore the signal. If the derived class needs to cache something, make it mutable

Reimplemented from genie::geometry::GeomVolSelectorFiducial.

Definition at line 91 of file GeomVolSelectorRockBox.cxx.

{
  // A new neutrino ray has been set, calculate the entrance/exit distances.

  GeomVolSelectorFiducial::BeginPSList(untrimmed);

  fCurrPathSegmentList = untrimmed;

  MakeRockBox();

  if ( ! fRockBoxShape ) {
    LOG("GeomVolSel", pFATAL) << "no shape defined";
    fInterceptRock = RayIntercept();
  } else {
    fInterceptRock =
      fRockBoxShape->Intercept(fCurrPathSegmentList->GetStartPos(),
                               fCurrPathSegmentList->GetDirection());
  }

  //cout << "BeginPSList: " << endl
  //     << " fid:  " << fIntercept << endl
  //     << " rock: " << fInterceptRock << endl;

}

void GeomVolSelectorRockBox::ConvertShapeMaster2Top ( const ROOTGeomAnalyzer *  rgeom )

virtual

Reimplemented from genie::geometry::GeomVolSelectorFiducial.

Definition at line 123 of file GeomVolSelectorRockBox.cxx.

{
  GeomVolSelectorFiducial::ConvertShapeMaster2Top(rgeom);
  fROOTGeom = rgeom;   // stash away a copy
}

void GeomVolSelectorRockBox::EndPSList ( ) const

virtual

Reimplemented from genie::geometry::GeomVolSelectorFiducial.

Definition at line 117 of file GeomVolSelectorRockBox.cxx.

{
  // Completed current path segment list processsing
}

void GeomVolSelectorRockBox::MakeRockBox ( ) const

protected

Definition at line 162 of file GeomVolSelectorRockBox.cxx.

{
  // This sets parameters for a box

  // expanded box
  double energy = fP4.Energy();
  double boxXYZMin[3], boxXYZMax[3];
  for ( int j = 0; j < 3; ++j ) {
    double dmin = 0, dmax = 0;
    double dircos = fCurrPathSegmentList->GetDirection()[j];
    if ( dircos > 0 ) dmin =  dircos*energy/fDeDx;  // pad upstream
    else              dmax = -dircos*energy/fDeDx;

//#define RWH_DEBUG
#ifdef RWH_DEBUG
    cout << "MakeRockBox [" << j << "] wall " << fMinimumWall
         << " dm " << dmin << " " << dmax   << " dircos " << dircos
         << " e " << energy << " dedx " << fDeDx << endl;
#endif

    if ( fExpandInclusion ) {
      boxXYZMin[j] = fInclusionXYZMin[j] - dmin;
      boxXYZMax[j] = fInclusionXYZMax[j] + dmax;
    } else {
      boxXYZMin[j] = TMath::Min(fMinimalXYZMin[j]-dmin,fInclusionXYZMin[j]);
      boxXYZMax[j] = TMath::Max(fMinimalXYZMax[j]+dmin,fInclusionXYZMax[j]);
    }
  }

  FidPolyhedron* poly = new FidPolyhedron();
  // careful about sign of "d" vs. direction normal
  PlaneParam pln0(-1,0,0, boxXYZMin[0]);  poly->push_back(pln0);
  PlaneParam pln1(0,-1,0, boxXYZMin[1]);  poly->push_back(pln1);
  PlaneParam pln2(0,0,-1, boxXYZMin[2]);  poly->push_back(pln2);
  PlaneParam pln3(+1,0,0,-boxXYZMax[0]);  poly->push_back(pln3);
  PlaneParam pln4(0,+1,0,-boxXYZMax[1]);  poly->push_back(pln4);
  PlaneParam pln5(0,0,+1,-boxXYZMax[2]);  poly->push_back(pln5);

  if ( fRockBoxShape ) delete fRockBoxShape;
  fRockBoxShape = poly;

#ifdef RWH_DEBUG
  static bool first = true;
  if ( first ) {
    cout << "MakeRockBox first Minimal min ["
         << fMinimalXYZMin[0] << ","
         << fMinimalXYZMin[1] << ","
         << fMinimalXYZMin[2] << "]"
         << "  max ["
         << fMinimalXYZMax[0] << ","
         << fMinimalXYZMax[1] << ","
         << fMinimalXYZMax[2] << "]" << endl;
    cout << "MakeRockBox first Inclusion min ["
         << fInclusionXYZMin[0] << ","
         << fInclusionXYZMin[1] << ","
         << fInclusionXYZMin[2] << "]"
         << "  max ["
         << fInclusionXYZMax[0] << ","
         << fInclusionXYZMax[1] << ","
         << fInclusionXYZMax[2] << "]" << endl;
    first = false;
  }
  cout << "MakeRockBox this ray using min ["
       << boxXYZMin[0] << ","
       << boxXYZMin[1] << ","
       << boxXYZMin[2] << "]"
       << "  max ["
       << boxXYZMax[0] << ","
       << boxXYZMax[1] << ","
       << boxXYZMax[2] << "]" << endl;
  cout << "rock before:" << *fRockBoxShape << endl;
#endif

  if ( fROOTGeom ) fRockBoxShape->ConvertMaster2Top(fROOTGeom);

#ifdef RWH_DEBUG
  cout << "rock after: " << *fRockBoxShape << endl;
  cout << "fid  after: " << *fShape << endl;
#endif

}

void genie::geometry::GeomVolSelectorRockBox::SetDeDx ( Double_t  dedx )

inline

Definition at line 59 of file GeomVolSelectorRockBox.h.

{ fDeDx = dedx; }

void genie::geometry::GeomVolSelectorRockBox::SetExpandFromInclusion ( bool  how = false )

inline

Definition at line 60 of file GeomVolSelectorRockBox.h.

{ fExpandInclusion = how; }

void GeomVolSelectorRockBox::SetMinimumWall

(

Double_t

w

)

Definition at line 153 of file GeomVolSelectorRockBox.cxx.

{
  fMinimumWall = w;
  for ( int j = 0; j < 3; ++j ) {
    fInclusionXYZMin[j] = fMinimalXYZMin[j] - fMinimumWall;
    fInclusionXYZMax[j] = fMinimalXYZMax[j] + fMinimumWall;
  }
}

void GeomVolSelectorRockBox::SetRockBoxInclusion	(	Double_t *	xyzmin,
		Double_t *	xyzmax
	)

Definition at line 142 of file GeomVolSelectorRockBox.cxx.

{
  // This sets parameters for the inclusion (outer) box

  for ( int j = 0; j < 3; ++j ) {
    fInclusionXYZMin[j] = TMath::Min(xyzmin[j],xyzmax[j]);
    fInclusionXYZMax[j] = TMath::Max(xyzmax[j],xyzmax[j]);
  }
}

void GeomVolSelectorRockBox::SetRockBoxMinimal	(	Double_t *	xyzmin,
		Double_t *	xyzmax
	)

Definition at line 129 of file GeomVolSelectorRockBox.cxx.

{
  // This sets parameters for a minimal box

  for ( int j = 0; j < 3; ++j ) {
    fMinimalXYZMin[j] = TMath::Min(xyzmin[j],xyzmax[j]);
    fMinimalXYZMax[j] = TMath::Max(xyzmax[j],xyzmax[j]);
  }
  // create the default inner (exclusion) box
  this->MakeBox(fMinimalXYZMin,fMinimalXYZMax);
}

void GeomVolSelectorRockBox::TrimSegment ( PathSegment &  segment ) const

virtual

This is the method every derived version must implement To reject a segment outright: segment.fStepRangeSet.clear()

so a miss means blindly accept all segments

Reimplemented from genie::geometry::GeomVolSelectorFiducial.

Definition at line 48 of file GeomVolSelectorRockBox.cxx.

{
  // First trim the segment based on the ray vs. cylinder or box
  // Then trim futher according to the Basic parameters

  if ( ! fInterceptRock.fIsHit ) {
      // want in rock box, ray misses => reject all segments
      ps.fStepRangeSet.clear();  //
  } else {
    // ray hit rock box volume, some segments steps need rejection, some need splitting...
    // check the steps in this segment
    Double_t dist = ps.fRayDist;
    StepRangeSet::iterator srs_itr = ps.fStepRangeSet.begin();
    StepRangeSet::iterator srs_end = ps.fStepRangeSet.end();
    StepRangeSet modifiedStepRangeSet;
    Bool_t ismod = false;

    // loop over steps within this segement
    for ( ; srs_itr != srs_end; ++srs_itr ) {
      Double_t slo = srs_itr->first;
      Double_t shi = srs_itr->second;
      Bool_t split = false;
      StepRange step1, step2;
      // determine new trimmed or split steps
      ismod |= NewStepPairs(false,dist,slo,shi,
                            fInterceptRock,split,step1,step2);
      // build up new step list
      bool nonzerostep = ( step1.first != step1.second );
      if ( nonzerostep || ! fRemoveEntries ) {
        modifiedStepRangeSet.push_back(step1);
        if (split) {
          modifiedStepRangeSet.push_back(step2);
        }
      }
    } // loop over step range set elements
    if ( ismod ) ps.fStepRangeSet = modifiedStepRangeSet;

  } // fIsHit

  GeomVolSelectorFiducial::TrimSegment(ps);
}

Member Data Documentation

Double_t genie::geometry::GeomVolSelectorRockBox::fDeDx

protected

accepted

Definition at line 77 of file GeomVolSelectorRockBox.h.

Bool_t genie::geometry::GeomVolSelectorRockBox::fExpandInclusion

protected

how to scale from energy to distance

Definition at line 78 of file GeomVolSelectorRockBox.h.

Double_t genie::geometry::GeomVolSelectorRockBox::fInclusionXYZMax[3]

protected

box within which events are always

Definition at line 76 of file GeomVolSelectorRockBox.h.

Double_t genie::geometry::GeomVolSelectorRockBox::fInclusionXYZMin[3]

protected

minimum distance around (XYZmin,XYZmax)

Definition at line 75 of file GeomVolSelectorRockBox.h.

RayIntercept genie::geometry::GeomVolSelectorRockBox::fInterceptRock

mutableprotected

Definition at line 85 of file GeomVolSelectorRockBox.h.

Double_t genie::geometry::GeomVolSelectorRockBox::fMinimalXYZMax[3]

protected

interior box lower corner

Definition at line 73 of file GeomVolSelectorRockBox.h.

Double_t genie::geometry::GeomVolSelectorRockBox::fMinimalXYZMin[3]

protected

Definition at line 72 of file GeomVolSelectorRockBox.h.

Double_t genie::geometry::GeomVolSelectorRockBox::fMinimumWall

protected

interior box upper corner

Definition at line 74 of file GeomVolSelectorRockBox.h.

FidShape* genie::geometry::GeomVolSelectorRockBox::fRockBoxShape

mutableprotected

expand from minimal or inclusion box?

Definition at line 80 of file GeomVolSelectorRockBox.h.

const ROOTGeomAnalyzer* genie::geometry::GeomVolSelectorRockBox::fROOTGeom

protected

shape changes for every nu ray

Definition at line 82 of file GeomVolSelectorRockBox.h.

---

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

• Generator/src/Tools/Geometry/GeomVolSelectorRockBox.h
• Generator/src/Tools/Geometry/GeomVolSelectorRockBox.cxx