LinearityADC Class Reference

parameterization of correction for ADC non linearity More...

#include <LinearityADC.h>

Inheritance diagram for LinearityADC:

Collaboration diagram for LinearityADC:

Public Member Functions
	LinearityADC (const std::string &filename, const std::string &name="ADCNonLinearityCorrection")
	constructor (initialization done there reading a root file for the HV maps per period
	~LinearityADC ()
float	getCorr (float etaCalo, float et, PATCore::ParticleType::Type ptype=PATCore::ParticleType::Electron) const
	get relative correction factor to apply to energy: corrected energy should be original energy * correction factor
void	setLevel (MSG::Level lvl)
	Change the current logging level.
Functions providing the same interface as AthMessaging
bool	msgLvl (const MSG::Level lvl) const
	Test the output level of the object.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
std::map< PATCore::ParticleType::Type, std::unique_ptr< TProfile2D > >	m_hcorr
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels)
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging)

Detailed Description

parameterization of correction for ADC non linearity

deltaADC vs ADC parameterized from calibration runs in different eta bins for HG/MG/HG by Luca applied to single particle samples to parameterize vs eta and Et and particle type

Definition at line 25 of file LinearityADC.h.

Constructor & Destructor Documentation

LinearityADC()

LinearityADC::LinearityADC	(	const std::string &	filename,
		const std::string &	name = "ADCNonLinearityCorrection" )

constructor (initialization done there reading a root file for the HV maps per period

Definition at line 9 of file LinearityADC.cxx.

    : asg::AsgMessaging(name) {
  ATH_MSG_INFO("Initialize LinearityADC with file " << filename);
  std::unique_ptr<TFile> adcFile(TFile::Open(filename.c_str(), "READ"));
  static const std::map<TString, PATCore::ParticleType::Type> pN = {
      {"elec", PATCore::ParticleType::Electron},
      {"unco", PATCore::ParticleType::UnconvertedPhoton},
      {"conv", PATCore::ParticleType::ConvertedPhoton}};
  for (const auto & [ptclName, ptclType] : pN) {
    m_hcorr[ptclType].reset(dynamic_cast<TProfile2D*>(
        adcFile->Get(Form("linearity_%s", ptclName.Data()))));
    m_hcorr[ptclType]->SetDirectory(nullptr);
  }
  adcFile->Close();
}

~LinearityADC()

LinearityADC::~LinearityADC ( )

inline

Definition at line 32 of file LinearityADC.h.

{};

Member Function Documentation

getCorr()

float LinearityADC::getCorr	(	float	etaCalo,
		float	et,
		PATCore::ParticleType::Type	ptype = PATCore::ParticleType::Electron ) const

get relative correction factor to apply to energy: corrected energy should be original energy * correction factor

inputs: etaCalo and Et in MeV, pType (enum)

Definition at line 26 of file LinearityADC.cxx.

                                                                 {
  float corr = 0.;
  float etGeV = et / 1000.;
  if (etGeV < 6.)
    etGeV = 6.;
  if (etGeV > 2999.)
    etGeV = 2999.;
 
  auto it = m_hcorr.find(pType);
  if (it == m_hcorr.end()) {
    ATH_MSG_ERROR("unknown particle type " << pType);
    return 1.;
  }
  TProfile2D* hh = it->second.get();
 
  float aeta = std::abs(etaCalo);
  if (aeta > 2.49)
    aeta = 2.49;
 
  int ix = hh->GetXaxis()->FindBin(aeta);
  int iy = hh->GetYaxis()->FindBin(etGeV);
  float y = hh->GetYaxis()->GetBinCenter(iy);
  int iy0, iy1;
  if (y > etGeV) {
    iy0 = iy - 1;
    iy1 = iy;
  } else {
    iy0 = iy;
    iy1 = iy + 1;
  }
 
  int iyc = hh->GetYaxis()->FindBin(40.);
 
  if (iy0 < 1)
    corr = hh->GetBinContent(ix, iy1);
  if (iy1 > hh->GetYaxis()->GetNbins())
    corr = hh->GetBinContent(ix, iy0);
 
  float y0 = hh->GetYaxis()->GetBinCenter(iy0);
  float y1 = hh->GetYaxis()->GetBinCenter(iy1);
  // subtract effect for 40 GeV Et electrons
  it = m_hcorr.find(PATCore::ParticleType::Electron);
  if (it == m_hcorr.end()){
    ATH_MSG_ERROR("Unknown particle type " << pType<<" in LinearityADC::getCorr");
    return 1.;
  }
  float c0 = hh->GetBinContent(ix, iy0) - it->second->GetBinContent(ix, iyc);
  float c1 = hh->GetBinContent(ix, iy1) - it->second->GetBinContent(ix, iyc);
 
  corr = (c1 * (etGeV - y0) + c0 * (y1 - etGeV)) / (y1 - y0);
  ATH_MSG_VERBOSE("correction from ADC non linearity = " << corr);
 
  return (1. + corr);
}

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

msg() [1/2]

MsgStream & asg::AsgMessaging::msg ( ) const

inherited

The standard message stream.

Returns

A reference to the default message stream of this object.

Definition at line 49 of file AsgMessaging.cxx.

                                      {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg();
#else // not XAOD_STANDALONE
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msg() [2/2]

MsgStream & asg::AsgMessaging::msg ( const MSG::Level lvl ) const

inherited

The standard message stream.

Parameters

lvl	The message level to set the stream to

Returns

A reference to the default message stream, set to level "lvl"

Definition at line 57 of file AsgMessaging.cxx.

                                                          {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg( lvl );
#else // not XAOD_STANDALONE
      m_msg << lvl;
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msgLvl()

bool asg::AsgMessaging::msgLvl ( const MSG::Level lvl ) const

inherited

Test the output level of the object.

Parameters

lvl	The message level to test against

Returns

boolean Indicting if messages at given level will be printed

true If messages at level "lvl" will be printed

Definition at line 41 of file AsgMessaging.cxx.

                                                       {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msgLvl( lvl );
#else // not XAOD_STANDALONE
      return m_msg.msgLevel( lvl );
#endif // not XAOD_STANDALONE
   }

setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_hcorr

std::map<PATCore::ParticleType::Type, std::unique_ptr<TProfile2D> > LinearityADC::m_hcorr

private

Definition at line 42 of file LinearityADC.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

---

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

• LinearityADC.h
• LinearityADC.cxx