#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