gar::geo::seg::SegmentationMuIDAlg Class Reference

#include <SegmentationMuIDAlg.h>

Inheritance diagram for gar::geo::seg::SegmentationMuIDAlg:

Public Member Functions
	SegmentationMuIDAlg (fhicl::ParameterSet const &pset)
	Default constructor used by derived classes passing the encoding string. More...
	SegmentationMuIDAlg (const BitFieldCoder *decoder, fhicl::ParameterSet const &pset)
	Default constructor used by derived classes passing an existing decoder. More...
	~SegmentationMuIDAlg ()
void	reconfigure (fhicl::ParameterSet const &pset) override
void	Initialize (const gar::geo::GeometryCore &geo) override
std::array< double, 3 >	GetPosition (const gar::geo::GeometryCore &geo, const gar::raw::CellID_t &cID) const override
gar::raw::CellID_t	GetCellID (const gar::geo::GeometryCore &geo, const unsigned int &det_id, const unsigned int &stave, const unsigned int &module, const unsigned int &layer, const unsigned int &slice, const std::array< double, 3 > &localPosition) const override
	determine the cell ID based on the position More...
bool	isTile (const gar::raw::CellID_t &cID) const override
bool	isBarrel (const gar::raw::CellID_t &cID) const override
double	getStripLength (const gar::geo::GeometryCore &geo, const std::array< double, 3 > &local, const gar::raw::CellID_t &cID) const override
std::pair< TVector3, TVector3 >	getStripEnds (const gar::geo::GeometryCore &geo, const std::array< double, 3 > &local, const gar::raw::CellID_t &cID) const override
std::pair< float, float >	CalculateLightPropagation (const gar::geo::GeometryCore &geo, const std::array< double, 3 > &local, const gar::raw::CellID_t &cID) const override
std::array< double, 3 >	ReconstructStripHitPosition (const gar::geo::GeometryCore &geo, const std::array< double, 3 > &local, const float &xlocal, const gar::raw::CellID_t &cID) const override
const double &	stripSizeX () const override
const double &	stripSizeY () const
const double &	layerDimX () const
const double &	layerDimY () const
const std::string &	fieldNameX () const
const std::string &	fieldNameY () const
const std::string &	fieldNameLayer () const
const std::string &	fieldNameSlice () const
const unsigned int &	nLayers () const override
void	setStripSizeY (double stripSize)
void	setFieldNameX (const std::string &fieldName)
void	setFieldNameY (const std::string &fieldName)
void	setFieldNameLayer (const std::string &fieldName)
void	setFieldNameSlice (const std::string &fieldName)
void	setLayerDimXY (const double &dimX, const double &dimY) const override
void	setVariables (const double &, const double &) const override
Public Member Functions inherited from gar::geo::seg::SegmentationAlg
virtual	~SegmentationAlg ()
virtual const std::string &	name () const
virtual void	setName (const std::string &value)
virtual const std::string &	type () const
virtual const std::string &	description () const
virtual const BitFieldCoder *	decoder () const
virtual const double &	gridSizeX () const
virtual const unsigned int &	nPlanes () const
virtual const std::string &	cellEncoding () const
virtual void	setDecoder (const BitFieldCoder *decoder)

Protected Member Functions
void	PrintParameters () const override
Protected Member Functions inherited from gar::geo::seg::SegmentationAlg
	SegmentationAlg (fhicl::ParameterSet const &pset)
	SegmentationAlg (const BitFieldCoder *decoder, fhicl::ParameterSet const &pset)

Protected Attributes
std::string	_xId
	the field name used for X More...
std::string	_yId
	the field name used for Y More...
std::string	_layerId
	the field name used for layer More...
std::string	_sliceId
	the field name used for slice More...
std::string	_encoding
	the encoding string More...
double	_stripSizeX
	the strip size in X More...
double	_stripSizeY
	the strip size in Y More...
double	_frac
	fraction of tiles to remove at the edge More...
unsigned int	_nLayers
	number of layers More...
double	_layer_dim_X
	layer dimension in X More...
double	_layer_dim_Y
	layer dimension in Y More...
Protected Attributes inherited from gar::geo::seg::SegmentationAlg
std::string	_name
std::string	_type
std::string	_description
std::string	_encoding
double	_gridSizeX
double	_stripSizeX
unsigned int	_nLayers
unsigned int	_nPlanes
const BitFieldCoder *	_decoder = 0

Additional Inherited Members
Static Protected Member Functions inherited from gar::geo::seg::SegmentationAlg
static double	binToPosition (gar::raw::CellID_t bin, double cellSize, double offset=0)
static int	positionToBin (double position, double cellSize, double offset=0)

Detailed Description

Definition at line 18 of file SegmentationMuIDAlg.h.

Constructor & Destructor Documentation

gar::geo::seg::SegmentationMuIDAlg::SegmentationMuIDAlg

(

fhicl::ParameterSet const &

pset

)

Default constructor used by derived classes passing the encoding string.

Definition at line 20 of file SegmentationMuIDAlg.cxx.

            : SegmentationAlg(pset)
            {
                _type = "MuID";
                _description = "Muon system segmentation with cross-strips depending on the slice";

                std::cout << "######### gar::geo::seg::SegmentationMuIDAlg() " << std::endl ;

                this->reconfigure(pset);
            }

gar::geo::seg::SegmentationMuIDAlg::SegmentationMuIDAlg	(	const BitFieldCoder *	decoder,
		fhicl::ParameterSet const &	pset
	)

Default constructor used by derived classes passing an existing decoder.

Definition at line 33 of file SegmentationMuIDAlg.cxx.

            : SegmentationAlg(decode, pset)
            {
                _type = "MuID";
                _description = "Muon system segmentation with cross-strips depending on the slice";

                std::cout << "######### gar::geo::seg::SegmentationMuIDAlg() " << std::endl ;

                this->reconfigure(pset);
            }

gar::geo::seg::SegmentationMuIDAlg::~SegmentationMuIDAlg

(

)

Definition at line 45 of file SegmentationMuIDAlg.cxx.

            {
            }

Member Function Documentation

std::pair< float, float > gar::geo::seg::SegmentationMuIDAlg::CalculateLightPropagation ( const gar::geo::GeometryCore &  geo, const std::array< double, 3 > &  local, const gar::raw::CellID_t &  cID  ) const

overridevirtual

Reimplemented from gar::geo::seg::SegmentationAlg.

Definition at line 220 of file SegmentationMuIDAlg.cxx.

            {
                float c = (CLHEP::c_light * CLHEP::mm / CLHEP::ns) / CLHEP::cm; // in cm/ns
                //time1 is left SiPM
                float time1 = 0.;
                //time2 is right SiPM
                float time2 = 0.;

                int layer = _decoder->get(cID, "layer");

                if(layer == 1 || layer == 3) {
                    time1 = ( _layer_dim_X / 2 + local[0] ) / c;
                    time2 = ( _layer_dim_X / 2 - local[0] ) / c;
                }

                if(layer == 2) {
                    time1 = ( _layer_dim_Y / 2 + local[1] ) / c;
                    time2 = ( _layer_dim_Y / 2 - local[1] ) / c;
                }

                return std::make_pair(time1, time2);
            }

const std::string& gar::geo::seg::SegmentationMuIDAlg::fieldNameLayer ( ) const

inline

Definition at line 59 of file SegmentationMuIDAlg.h.

{ return _layerId; }

const std::string& gar::geo::seg::SegmentationMuIDAlg::fieldNameSlice ( ) const

inline

Definition at line 61 of file SegmentationMuIDAlg.h.

{ return _sliceId; }

const std::string& gar::geo::seg::SegmentationMuIDAlg::fieldNameX ( ) const

inline

Definition at line 55 of file SegmentationMuIDAlg.h.

{ return _xId; }

const std::string& gar::geo::seg::SegmentationMuIDAlg::fieldNameY ( ) const

inline

Definition at line 57 of file SegmentationMuIDAlg.h.

{ return _yId; }

gar::raw::CellID_t gar::geo::seg::SegmentationMuIDAlg::GetCellID ( const gar::geo::GeometryCore &  geo, const unsigned int &  det_id, const unsigned int &  stave, const unsigned int &  module, const unsigned int &  layer, const unsigned int &  slice, const std::array< double, 3 > &  localPosition  ) const

overridevirtual

determine the cell ID based on the position

Implements gar::geo::seg::SegmentationAlg.

Definition at line 111 of file SegmentationMuIDAlg.cxx.

            {
                gar::raw::CellID_t cID = 0;

                _decoder->set(cID, "system", det_id);
                _decoder->set(cID, "module", module);
                _decoder->set(cID, "stave", stave);
                _decoder->set(cID, "layer", layer);
                _decoder->set(cID, "slice", slice);

                double localX = localPosition[0];
                double localY = localPosition[1];

                //Segmentation depending on the slice in the layer (1 layer with Sc/Fe/Sc/Fe/Sc)
                //Slices 1/3/5 active
                int nCellsX = 1;
                int nCellsY = 1;

                if(layer % 2 == 0) {
                    //Segmentation in Y
                    nCellsX = 1;
                    nCellsY = int(_layer_dim_Y / _stripSizeY);

                }

                if(layer % 2 != 0) {
                    //Segmentation in X
                    nCellsY = 1;
                    nCellsX = int(_layer_dim_X / _stripSizeX);
                }

                int _cellIndexX = int ( localX / ( _layer_dim_X / nCellsX ) );
                int _cellIndexY = int ( localY / ( _layer_dim_Y / nCellsY ) );

                _decoder->set(cID, _xId, _cellIndexX);
                _decoder->set(cID, _yId, _cellIndexY);

                return cID;
            }

std::array< double, 3 > gar::geo::seg::SegmentationMuIDAlg::GetPosition ( const gar::geo::GeometryCore &  geo, const gar::raw::CellID_t &  cID  ) const

overridevirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 77 of file SegmentationMuIDAlg.cxx.

            {
                std::array<double, 3> cellPosition;

                int cellIndexX = _decoder->get(cID, _xId);
                int cellIndexY = _decoder->get(cID, _yId);
                int layer = _decoder->get(cID, "layer");

                if (layer % 2 == 0)
                {
                    //Segmentation in Y
                    int nCellsX = 1;
                    // int nCellsY = int(_layer_dim_Y / _stripSizeY);

                    cellPosition[0] = ( cellIndexX + 0.5 ) * (_layer_dim_X / nCellsX ) - (_layer_dim_X / 2.);
                    cellPosition[1] = ( cellIndexY + 0.5 ) * _stripSizeY;
                    cellPosition[2] = 0.;
                }

                if(layer % 2 != 0) {
                    //Segmentation in X
                    int nCellsY = 1;
                    // int nCellsX = int(_layer_dim_Y / _stripSizeY);

                    cellPosition[0] = ( cellIndexX + 0.5 ) * _stripSizeX;
                    cellPosition[1] = ( cellIndexY + 0.5 ) * (_layer_dim_Y / nCellsY ) - (_layer_dim_Y / 2.);
                    cellPosition[2] = 0.;
                }

                return cellPosition;
            }

std::pair< TVector3, TVector3 > gar::geo::seg::SegmentationMuIDAlg::getStripEnds ( const gar::geo::GeometryCore &  geo, const std::array< double, 3 > &  local, const gar::raw::CellID_t &  cID  ) const

overridevirtual

Reimplemented from gar::geo::seg::SegmentationAlg.

Definition at line 190 of file SegmentationMuIDAlg.cxx.

            {
                int layer = _decoder->get(cID, "layer");

                TVector3 stripEnd1(0., 0., 0.);
                TVector3 stripEnd2(0., 0., 0.);

                if(layer == 1 || layer == 3) {
                    stripEnd1.SetX(-_layer_dim_X/2.);
                    stripEnd2.SetX(_layer_dim_X/2.);
                    stripEnd1.SetY(local[1]);
                    stripEnd2.SetY(local[1]);
                    stripEnd1.SetZ(local[2]);
                    stripEnd2.SetZ(local[2]);
                }

                if(layer == 2) {
                    stripEnd1.SetY(-_layer_dim_Y/2.);
                    stripEnd2.SetY(_layer_dim_Y/2.);
                    stripEnd1.SetX(local[0]);
                    stripEnd2.SetX(local[0]);
                    stripEnd1.SetZ(local[2]);
                    stripEnd2.SetZ(local[2]);
                }


                return std::make_pair(stripEnd1, stripEnd2);
            }

double gar::geo::seg::SegmentationMuIDAlg::getStripLength ( const gar::geo::GeometryCore &  geo, const std::array< double, 3 > &  local, const gar::raw::CellID_t &  cID  ) const

overridevirtual

Reimplemented from gar::geo::seg::SegmentationAlg.

Definition at line 178 of file SegmentationMuIDAlg.cxx.

            {
                double stripLength = 0.;
                int layer = _decoder->get(cID, "layer");

                if(layer == 1 || layer == 3) stripLength = _layer_dim_X;
                if(layer == 2) stripLength = _layer_dim_Y;

                return stripLength;
            }

void gar::geo::seg::SegmentationMuIDAlg::Initialize ( const gar::geo::GeometryCore &  geo )

overridevirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 71 of file SegmentationMuIDAlg.cxx.

            {

            }

bool gar::geo::seg::SegmentationMuIDAlg::isBarrel ( const gar::raw::CellID_t &  cID ) const

overridevirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 172 of file SegmentationMuIDAlg.cxx.

            {
                return true;
            }

bool gar::geo::seg::SegmentationMuIDAlg::isTile ( const gar::raw::CellID_t &  cID ) const

overridevirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 166 of file SegmentationMuIDAlg.cxx.

            {
                return false;
            }

const double& gar::geo::seg::SegmentationMuIDAlg::layerDimX ( ) const

inline

Definition at line 51 of file SegmentationMuIDAlg.h.

{ return _layer_dim_X; }

const double& gar::geo::seg::SegmentationMuIDAlg::layerDimY ( ) const

inline

Definition at line 53 of file SegmentationMuIDAlg.h.

{ return _layer_dim_Y; }

const unsigned int& gar::geo::seg::SegmentationMuIDAlg::nLayers ( ) const

inlineoverridevirtual

Reimplemented from gar::geo::seg::SegmentationAlg.

Definition at line 63 of file SegmentationMuIDAlg.h.

{ return _nLayers; }

void gar::geo::seg::SegmentationMuIDAlg::PrintParameters ( ) const

overrideprotectedvirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 152 of file SegmentationMuIDAlg.cxx.

            {
                std::cout << "cell encoding: " << _encoding << std::endl;
                std::cout << "identifier_x: " << _xId << std::endl;
                std::cout << "identifier_y: " << _yId << std::endl;
                std::cout << "strip_size_x: " << _stripSizeX << " cm" << std::endl;
                std::cout << "strip_size_y: " << _stripSizeY << " cm" << std::endl;

                std::cout << "identifier_layer: " << _layerId << std::endl;
                std::cout << "identifier_slice: " << _sliceId << std::endl;
                std::cout << "nLayers: " << _nLayers << std::endl;
            }

void gar::geo::seg::SegmentationMuIDAlg::reconfigure ( fhicl::ParameterSet const &  pset )

overridevirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 50 of file SegmentationMuIDAlg.cxx.

            {
                _stripSizeY = pset.get<double>("strip_size_y");
                _stripSizeX = _stripSizeY;
                _encoding = pset.get<std::string>("cellEncoding");

                _xId = pset.get<std::string>("identifier_x");
                _yId = pset.get<std::string>("identifier_y");
                _layerId = pset.get<std::string>("identifier_layer");
                _sliceId = pset.get<std::string>("identifier_slice");

                _nLayers = pset.get<unsigned int>("nlayers");

                _frac = 1./3.;

                this->PrintParameters();

                return;
            }

std::array< double, 3 > gar::geo::seg::SegmentationMuIDAlg::ReconstructStripHitPosition ( const gar::geo::GeometryCore &  geo, const std::array< double, 3 > &  local, const float &  xlocal, const gar::raw::CellID_t &  cID  ) const

overridevirtual

Reimplemented from gar::geo::seg::SegmentationAlg.

Definition at line 244 of file SegmentationMuIDAlg.cxx.

            {
                int layer = _decoder->get(cID, "layer");
                float pos = xlocal;
                std::array<double, 3> newlocal;

                //Need to check if the local position is bigger than the strip length!
                // bool isBarrel = this->isBarrel(cID);
                float stripLength = this->getStripLength(geo, local, cID);

                if( std::abs(pos) > stripLength / 2. ) pos = (pos > 0) ? stripLength / 2. : -stripLength / 2.;

                if(layer == 1 || layer == 3) {
                    newlocal = {pos, local[1], local[2]};
                }

                if(layer == 2) {
                    newlocal = {local[0], pos, local[2]};
                }

                return newlocal;
            }

void gar::geo::seg::SegmentationMuIDAlg::setFieldNameLayer ( const std::string &  fieldName )

inline

Definition at line 71 of file SegmentationMuIDAlg.h.

{ _layerId = fieldName; }

void gar::geo::seg::SegmentationMuIDAlg::setFieldNameSlice ( const std::string &  fieldName )

inline

Definition at line 73 of file SegmentationMuIDAlg.h.

{ _sliceId = fieldName; }

void gar::geo::seg::SegmentationMuIDAlg::setFieldNameX ( const std::string &  fieldName )

inline

Definition at line 67 of file SegmentationMuIDAlg.h.

{ _xId = fieldName; }

void gar::geo::seg::SegmentationMuIDAlg::setFieldNameY ( const std::string &  fieldName )

inline

Definition at line 69 of file SegmentationMuIDAlg.h.

{ _yId = fieldName; }

void gar::geo::seg::SegmentationMuIDAlg::setLayerDimXY ( const double &  dimX, const double &  dimY  ) const

inlineoverridevirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 75 of file SegmentationMuIDAlg.h.

{ _layer_dim_X = dimX; _layer_dim_Y = dimY; }

void gar::geo::seg::SegmentationMuIDAlg::setStripSizeY ( double  stripSize )

inline

Definition at line 65 of file SegmentationMuIDAlg.h.

{ _stripSizeY = stripSize; _stripSizeX = stripSize; }

void gar::geo::seg::SegmentationMuIDAlg::setVariables ( const double &  , const double &    ) const

inlineoverridevirtual

Implements gar::geo::seg::SegmentationAlg.

Definition at line 78 of file SegmentationMuIDAlg.h.

{ /* no op */ }

const double& gar::geo::seg::SegmentationMuIDAlg::stripSizeX ( ) const

inlineoverridevirtual

Reimplemented from gar::geo::seg::SegmentationAlg.

Definition at line 47 of file SegmentationMuIDAlg.h.

{ return _stripSizeX; }

const double& gar::geo::seg::SegmentationMuIDAlg::stripSizeY ( ) const

inline

Definition at line 49 of file SegmentationMuIDAlg.h.

{ return _stripSizeY; }

Member Data Documentation

std::string gar::geo::seg::SegmentationMuIDAlg::_encoding

protected

the encoding string

Definition at line 93 of file SegmentationMuIDAlg.h.

double gar::geo::seg::SegmentationMuIDAlg::_frac

protected

fraction of tiles to remove at the edge

Definition at line 99 of file SegmentationMuIDAlg.h.

double gar::geo::seg::SegmentationMuIDAlg::_layer_dim_X

mutableprotected

layer dimension in X

Definition at line 103 of file SegmentationMuIDAlg.h.

double gar::geo::seg::SegmentationMuIDAlg::_layer_dim_Y

mutableprotected

layer dimension in Y

Definition at line 105 of file SegmentationMuIDAlg.h.

std::string gar::geo::seg::SegmentationMuIDAlg::_layerId

protected

the field name used for layer

Definition at line 89 of file SegmentationMuIDAlg.h.

unsigned int gar::geo::seg::SegmentationMuIDAlg::_nLayers

protected

number of layers

Definition at line 101 of file SegmentationMuIDAlg.h.

std::string gar::geo::seg::SegmentationMuIDAlg::_sliceId

protected

the field name used for slice

Definition at line 91 of file SegmentationMuIDAlg.h.

double gar::geo::seg::SegmentationMuIDAlg::_stripSizeX

protected

the strip size in X

Definition at line 95 of file SegmentationMuIDAlg.h.

double gar::geo::seg::SegmentationMuIDAlg::_stripSizeY

protected

the strip size in Y

Definition at line 97 of file SegmentationMuIDAlg.h.

std::string gar::geo::seg::SegmentationMuIDAlg::_xId

protected

the field name used for X

Definition at line 85 of file SegmentationMuIDAlg.h.

std::string gar::geo::seg::SegmentationMuIDAlg::_yId

protected

the field name used for Y

Definition at line 87 of file SegmentationMuIDAlg.h.

---

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

• garsoft/Geometry/ChannelMapAlgs/SegmentationMuIDAlg.h
• garsoft/Geometry/ChannelMapAlgs/SegmentationMuIDAlg.cxx