Utils Class Reference

#include <Utils.h>

Public Member Functions
	Utils ()
	Utils (const Utils &)=default
	~Utils ()
void	SetSeed (int seed)
void	SetOrigin (double *origin)
bool	PointInFiducial (const TVector3 &point)
bool	PointInTPC (const TVector3 &point)
bool	PointInCalo (const TVector3 &point)
bool	PointStopBetween (const TVector3 &point)
bool	isThroughCalo (const TVector3 &point)
bool	hasDecayedInCalo (const TVector3 &point)
bool	isBarrel (const TVector3 &point)
bool	isEndcap (const TVector3 &point)
bool	isBackscatter (const TVector3 &spoint, const TVector3 &epoint)
bool	isBremsstrahlung (const TVector3 &spoint, const int &pdg, const int &motherpdg)
float	GetRamdomNumber ()
float	GaussianSmearing (const float &mean, const float &sigma)
double *	GetOrigin ()

Private Attributes
double	_origin [3]
	coordinates of the origin More...
unsigned long int	_seed
	seed More...
TRandom3 *	_rando
	random generator More...
double	_TPCFidRadius = 222.5
double	_TPCFidLength = 215.
double	_TPCRadius = 273.
double	_TPCLength = 259.
double	_ECALInnerRadius = 278.02
double	_ECALOuterRadius = 322.25
double	_ECALEndX = 374.4
double	_ECALStartX = _ECALEndX - (_ECALOuterRadius - _ECALInnerRadius)

Detailed Description

Definition at line 7 of file Utils.h.

Constructor & Destructor Documentation

Utils::Utils

(

)

Definition at line 3 of file Utils.cpp.

: _seed(0), _rando(new TRandom3(0))
{
    _origin[0] = 0.;
    _origin[1] = 0.;
    _origin[2] = 0.;
}

Utils::Utils ( const Utils &  )

default

Utils::~Utils

(

)

Definition at line 11 of file Utils.cpp.

{
    delete _rando;
}

Member Function Documentation

float Utils::GaussianSmearing ( const float &  mean, const float &  sigma  )

inline

Definition at line 58 of file Utils.h.

{ return _rando->Gaus(mean, sigma); }

double* Utils::GetOrigin ( )

inline

Definition at line 60 of file Utils.h.

{ return &_origin[0]; }

float Utils::GetRamdomNumber ( )

inline

Definition at line 56 of file Utils.h.

{ return _rando->Rndm(); }

bool Utils::hasDecayedInCalo

(

const TVector3 &

point

)

Definition at line 91 of file Utils.cpp.

{
    return PointInCalo(point);
}

bool Utils::isBackscatter	(	const TVector3 &	spoint,
		const TVector3 &	epoint
	)

Definition at line 114 of file Utils.cpp.

{
    bool isBackscatter = false;

    //check if started in the ECAL but made it to the tracker
    float r_spoint = std::sqrt( spoint.Y()*spoint.Y() + spoint.Z()*spoint.Z() );
    float r_epoint = std::sqrt( epoint.Y()*epoint.Y() + epoint.Z()*epoint.Z() );

    //in the Barrel
    if( r_spoint > _ECALInnerRadius && r_epoint < _TPCRadius ) isBackscatter = true;
    //in the Endcap
    if( (r_spoint < _ECALInnerRadius && r_epoint < _TPCRadius) && ( std::abs(spoint.X()) > _ECALStartX && std::abs(epoint.X()) < _TPCLength ) ) isBackscatter = true;

    return isBackscatter;
}

bool Utils::isBarrel

(

const TVector3 &

point

)

Definition at line 96 of file Utils.cpp.

{
    bool isBarrel = false;
    float theta = std::atan(_ECALInnerRadius / std::abs(_ECALStartX) ); //angle for barrel/endcap transition
    float r_point = std::sqrt( point.Y()*point.Y() + point.Z()*point.Z() );
    float theta_point = std::atan(r_point / std::abs(point.X()) ); //angle for barrel/endcap transition for the point

    if( theta_point > theta ) isBarrel = true;
    return isBarrel;
}

bool Utils::isBremsstrahlung	(	const TVector3 &	spoint,
		const int &	pdg,
		const int &	motherpdg
	)

Definition at line 130 of file Utils.cpp.

{
    bool isBremsstrahlung = false;

    //Check if it has origin in the tracker and that the pdg is photon and mother is electron/positron (most probable)
    if(PointInFiducial(spoint) && pdg == 22 && std::abs(motherpdg) == 11) isBremsstrahlung = true;

    return isBremsstrahlung;
}

bool Utils::isEndcap

(

const TVector3 &

point

)

Definition at line 107 of file Utils.cpp.

{
    bool isEndcap = false;
    if( !isBarrel(point) ) isEndcap = true;
    return isEndcap;
}

bool Utils::isThroughCalo

(

const TVector3 &

point

)

Definition at line 86 of file Utils.cpp.

{
    return !PointInTPC(point) && !PointStopBetween(point) && !PointInCalo(point);
}

bool Utils::PointInCalo

(

const TVector3 &

point

)

Definition at line 58 of file Utils.cpp.

{
    //Barrel Radius 278 cm
    //Endcap starts at 364 cm
    bool isInCalo = false;
    float r_point = std::sqrt( point.Y()*point.Y() + point.Z()*point.Z() );
    //in the Barrel
    if( r_point > _ECALInnerRadius && r_point < _ECALOuterRadius && std::abs(point.X()) < _ECALStartX ) isInCalo = true;
    //in the Endcap
    if( r_point < _ECALInnerRadius && std::abs(point.X()) > _ECALStartX && std::abs(point.X()) < _ECALEndX ) isInCalo = true;

    return isInCalo;
}

bool Utils::PointInFiducial

(

const TVector3 &

point

)

Definition at line 29 of file Utils.cpp.

{
    //TPC Fiducial volume defined as
    //R < 260 cm
    //|X| < 250 cm
    bool isInFiducial = true;

    float r_point = std::sqrt( point.Y()*point.Y() + point.Z()*point.Z() );
    if( r_point > _TPCFidRadius ) isInFiducial = false;
    if( r_point < _TPCFidRadius && std::abs(point.X()) > _TPCFidLength ) isInFiducial = false;

    return isInFiducial;
}

bool Utils::PointInTPC

(

const TVector3 &

point

)

Definition at line 43 of file Utils.cpp.

{
    //TPC volume defined as
    //R < 260 cm
    //|X| < 250 cm
    if(PointInFiducial(point)) return true;
    bool isInTPC = true;

    float r_point = std::sqrt( point.Y()*point.Y() + point.Z()*point.Z() );
    if( r_point > _TPCRadius ) isInTPC = false;
    if( r_point < _TPCRadius && std::abs(point.X()) > _TPCLength ) isInTPC = false;

    return isInTPC;
}

bool Utils::PointStopBetween

(

const TVector3 &

point

)

Definition at line 72 of file Utils.cpp.

{
    //Barrel Radius 278 cm
    //Endcap starts at 364 cm
    bool isStopBetween = false;
    float r_point = std::sqrt( point.Y()*point.Y() + point.Z()*point.Z() );
    //in the Barrel
    if( r_point < _ECALInnerRadius && r_point > _TPCRadius && std::abs(point.X()) < _TPCLength ) isStopBetween = true;
    //in the Endcap
    if( r_point < _ECALInnerRadius && std::abs(point.X()) > _TPCLength && std::abs(point.X()) < _ECALStartX ) isStopBetween = true;

    return isStopBetween;
}

void Utils::SetOrigin

(

double *

origin

)

Definition at line 22 of file Utils.cpp.

{
    _origin[0] = origin[0];
    _origin[1] = origin[1];
    _origin[2] = origin[2];
}

void Utils::SetSeed

(

int

seed

)

Definition at line 16 of file Utils.cpp.

{
    _seed = seed;
    _rando->SetSeed(_seed);
}

Member Data Documentation

double Utils::_ECALEndX = 374.4

private

Definition at line 84 of file Utils.h.

double Utils::_ECALInnerRadius = 278.02

private

Definition at line 82 of file Utils.h.

double Utils::_ECALOuterRadius = 322.25

private

Definition at line 83 of file Utils.h.

double Utils::_ECALStartX = _ECALEndX - (_ECALOuterRadius - _ECALInnerRadius)

private

Definition at line 85 of file Utils.h.

double Utils::_origin[3]

private

coordinates of the origin

Definition at line 63 of file Utils.h.

TRandom3* Utils::_rando

private

random generator

Definition at line 65 of file Utils.h.

unsigned long int Utils::_seed

private

seed

Definition at line 64 of file Utils.h.

double Utils::_TPCFidLength = 215.

private

Definition at line 79 of file Utils.h.

double Utils::_TPCFidRadius = 222.5

private

Definition at line 78 of file Utils.h.

double Utils::_TPCLength = 259.

private

Definition at line 81 of file Utils.h.

double Utils::_TPCRadius = 273.

private

Definition at line 80 of file Utils.h.

---

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

• garsoft/Ana/ParamSim/include/Utils.h
• garsoft/Ana/ParamSim/src/Utils.cpp