#include <TruthTrackSmearer.h>

Collaboration diagram for L0Muon::TruthTrackSmearer:

Public Member Functions
	TruthTrackSmearer (ATHRNG::RNGWrapper *rndWrapper)
	~TruthTrackSmearer ()=default
	TruthTrackSmearer (const TruthTrackSmearer &other)=delete
TruthTrackSmearer &	operator= (const TruthTrackSmearer &other)=delete
double	effFunc (double pt) const
bool	emulateL0MuonTrack (double curv, float eta, float phi, L0MuonTrack &otrack) const

Private Attributes
ATHRNG::RNGWrapper *	m_rngWrapper
std::array< float, 2 >	m_efficiencyMap

Detailed Description

Definition at line 17 of file TruthTrackSmearer.h.

Constructor & Destructor Documentation

◆ TruthTrackSmearer() [1/2]

L0Muon::TruthTrackSmearer::TruthTrackSmearer ( ATHRNG::RNGWrapper * rndWrapper )

Definition at line 11 of file TruthTrackSmearer.cxx.

: m_rngWrapper(rngWrapper) {
  m_efficiencyMap[0] = 0.0;   // pT < 3 GeV
  m_efficiencyMap[1] = 0.95;  // pT > 3 GeV
}

◆ ~TruthTrackSmearer()

L0Muon::TruthTrackSmearer::~TruthTrackSmearer ( )

default

◆ TruthTrackSmearer() [2/2]

L0Muon::TruthTrackSmearer::TruthTrackSmearer ( const TruthTrackSmearer & other )

delete

Member Function Documentation

◆ effFunc()

double L0Muon::TruthTrackSmearer::effFunc ( double pt ) const

Definition at line 17 of file TruthTrackSmearer.cxx.

                                                 {
  // TODO: Efficiency is changed at only 3 GeV...
  return (pt < 3000.) ? m_efficiencyMap[0] : m_efficiencyMap[1];
}

◆ emulateL0MuonTrack()

bool L0Muon::TruthTrackSmearer::emulateL0MuonTrack	(	double	curv,
		float	eta,
		float	phi,
		L0MuonTrack &	otrack ) const

Definition at line 22 of file TruthTrackSmearer.cxx.

                                                                                                       {
  // input curv(q/pT) is in MeV
  CLHEP::HepRandomEngine* engine = m_rngWrapper->getEngine(Gaudi::Hive::currentContext());
 
  // efficiency emulation
  double abspt = std::abs(1.0 / curv);
  if (CLHEP::RandFlat::shoot(engine) > effFunc(abspt)) return false;
 
  // TODO: extrapolate to the pivot plane
 
  // position efficiency map
  if (std::abs(eta) > 2.41) return false;  // for the time being...
 
  // TODO: The position is NOT smeared here - since the binning of RoIs could be larger. To be checked.
 
  // pt smearing
  double sigma = curv * 0.05;  // TODO: 5% uniform for the time being...
 
  double gencurv = CLHEP::RandGauss::shoot(engine, curv, sigma);
 
  otrack.setTrack(gencurv, eta, phi);
 
  return true;
}

◆ operator=()

TruthTrackSmearer & L0Muon::TruthTrackSmearer::operator= ( const TruthTrackSmearer & other )

delete

Member Data Documentation

◆ m_efficiencyMap

std::array<float, 2> L0Muon::TruthTrackSmearer::m_efficiencyMap

private

Definition at line 31 of file TruthTrackSmearer.h.

◆ m_rngWrapper

ATHRNG::RNGWrapper* L0Muon::TruthTrackSmearer::m_rngWrapper