gar::CAFHelper Class Reference

Public Member Functions
	CAFHelper (const geo::GeometryCore *fGeo, CLHEP::HepRandomEngine &engine)
	~CAFHelper ()
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)
float	GetRamdomNumber ()
float	GaussianSmearing (const float mean, const float sigma)

Private Member Functions
void	PrintParameters ()

Private Attributes
CLHEP::HepRandomEngine &	fEngine
const double	fTPCFidRadius = 222.5
const double	fTPCFidLength = 215.
double	fTPCRadius
double	fTPCLength
double	fECALBarrelInnerRadius
double	fECALBarrelOuterRadius
double	fECALEndcapInnerRadius
double	fECALEndcapOuterRadius
double	fECALStartX
double	fECALEndX

Detailed Description

Definition at line 51 of file ParamSim_module.cc.

Constructor & Destructor Documentation

gar::CAFHelper::CAFHelper	(	const geo::GeometryCore *	fGeo,
		CLHEP::HepRandomEngine &	engine
	)

Definition at line 120 of file ParamSim_module.cc.

    : fEngine(engine)
  {
    fTPCRadius = fGeo->TPCRadius();
    fTPCLength = fGeo->TPCLength()/2.;
    fECALBarrelInnerRadius = fGeo->GetECALInnerBarrelRadius();
    fECALBarrelOuterRadius = fGeo->GetECALOuterBarrelRadius();
    fECALEndcapInnerRadius = fGeo->GetECALInnerEndcapRadius();
    fECALEndcapOuterRadius = fGeo->GetECALOuterEndcapRadius();
    fECALStartX = fGeo->GetECALEndcapStartX();
    fECALEndX = fGeo->GetECALEndcapOuterX();

    PrintParameters();
  }

gar::CAFHelper::~CAFHelper

(

)

Definition at line 136 of file ParamSim_module.cc.

  {

  }

Member Function Documentation

float gar::CAFHelper::GaussianSmearing ( const float  mean, const float  sigma  )

inline

Definition at line 114 of file ParamSim_module.cc.

                                                                              {
    CLHEP::RandGauss GausRand(fEngine);
    return GausRand.fire(mean, sigma);;
  }

float gar::CAFHelper::GetRamdomNumber ( )

inline

Definition at line 108 of file ParamSim_module.cc.

                                          {
    CLHEP::RandFlat FlatRand(fEngine);
    return FlatRand.fire();;
  }

bool gar::CAFHelper::hasDecayedInCalo ( const TVector3 &  point )

inline

Definition at line 224 of file ParamSim_module.cc.

  {
    return PointInCalo(point);
  }

bool gar::CAFHelper::isBarrel ( const TVector3 &  point )

inline

Definition at line 230 of file ParamSim_module.cc.

  {
    bool isBarrel = false;
    float theta = std::atan(fECALBarrelInnerRadius / std::abs(fECALStartX) ); //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 gar::CAFHelper::isEndcap ( const TVector3 &  point )

inline

Definition at line 242 of file ParamSim_module.cc.

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

bool gar::CAFHelper::isThroughCalo ( const TVector3 &  point )

inline

Definition at line 218 of file ParamSim_module.cc.

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

bool gar::CAFHelper::PointInCalo ( const TVector3 &  point )

inline

Definition at line 188 of file ParamSim_module.cc.

  {
    //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 > fECALBarrelInnerRadius && r_point < fECALBarrelOuterRadius && std::abs(point.X()) < fECALStartX ) isInCalo = true;
    //in the Endcap
    if( r_point < fECALEndcapOuterRadius && std::abs(point.X()) > fECALStartX && std::abs(point.X()) < fECALEndX ) isInCalo = true;

    return isInCalo;
  }

bool gar::CAFHelper::PointInFiducial ( const TVector3 &  point )

inline

Definition at line 157 of file ParamSim_module.cc.

  {
    //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 > fTPCFidRadius ) isInFiducial = false;
    if( r_point < fTPCFidRadius && std::abs(point.X()) > fTPCFidLength ) isInFiducial = false;

    return isInFiducial;
  }

bool gar::CAFHelper::PointInTPC ( const TVector3 &  point )

inline

Definition at line 172 of file ParamSim_module.cc.

  {
    //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 > fTPCRadius ) isInTPC = false;
    if( r_point < fTPCRadius && std::abs(point.X()) > fTPCLength ) isInTPC = false;

    return isInTPC;
  }

bool gar::CAFHelper::PointStopBetween ( const TVector3 &  point )

inline

Definition at line 203 of file ParamSim_module.cc.

  {
    //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 < fECALBarrelInnerRadius && r_point > fTPCRadius && std::abs(point.X()) < fTPCLength ) isStopBetween = true;
    //in the Endcap
    if( r_point < fECALEndcapOuterRadius && std::abs(point.X()) > fTPCLength && std::abs(point.X()) < fECALStartX ) isStopBetween = true;

    return isStopBetween;
  }

void gar::CAFHelper::PrintParameters ( )

private

Definition at line 142 of file ParamSim_module.cc.

  {
    std::cout << " ==== CAFHelper Parameters ==== " << std::endl;
    std::cout << "TPC Radius " << fTPCRadius << " cm" << std::endl;
    std::cout << "TPC Length " << fTPCLength << " cm" << std::endl;
    std::cout << "ECAL Barrel Inner Radius " << fECALBarrelInnerRadius << " cm" << std::endl;
    std::cout << "ECAL Barrel Outer Radius " << fECALBarrelOuterRadius << " cm" << std::endl;
    std::cout << "ECAL Endcap Inner Radius  " << fECALEndcapInnerRadius << " cm" << std::endl;
    std::cout << "ECAL Endcap Outer Radius  " << fECALEndcapOuterRadius << " cm" << std::endl;
    std::cout << "ECAL Endcap Start X " << fECALStartX << " cm" << std::endl;
    std::cout << "ECAL Endcap End X " << fECALEndX << " cm" << std::endl;
    std::cout << " ==== CAFHelper Parameters ==== " << std::endl;
  }

Member Data Documentation

double gar::CAFHelper::fECALBarrelInnerRadius

private

Definition at line 99 of file ParamSim_module.cc.

double gar::CAFHelper::fECALBarrelOuterRadius

private

Definition at line 100 of file ParamSim_module.cc.

double gar::CAFHelper::fECALEndcapInnerRadius

private

Definition at line 101 of file ParamSim_module.cc.

double gar::CAFHelper::fECALEndcapOuterRadius

private

Definition at line 102 of file ParamSim_module.cc.

double gar::CAFHelper::fECALEndX

private

Definition at line 104 of file ParamSim_module.cc.

double gar::CAFHelper::fECALStartX

private

Definition at line 103 of file ParamSim_module.cc.

CLHEP::HepRandomEngine& gar::CAFHelper::fEngine

private

Definition at line 91 of file ParamSim_module.cc.

const double gar::CAFHelper::fTPCFidLength = 215.

private

Definition at line 95 of file ParamSim_module.cc.

const double gar::CAFHelper::fTPCFidRadius = 222.5

private

Definition at line 94 of file ParamSim_module.cc.

double gar::CAFHelper::fTPCLength

private

Definition at line 98 of file ParamSim_module.cc.

double gar::CAFHelper::fTPCRadius

private

Definition at line 97 of file ParamSim_module.cc.

---

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

• garsoft/Ana/ParamSim_module.cc