TFCSEnergyInterpolationHistogram Class Reference

#include <TFCSEnergyInterpolationHistogram.h>

Inheritance diagram for TFCSEnergyInterpolationHistogram:

Collaboration diagram for TFCSEnergyInterpolationHistogram:

Public Types
enum	FCSEnergyInitializationStatusBits { kOnlyScaleEnergy = }
enum	FCSStatusBits { kMatchAllPDGID }
	Status bit for FCS needs. More...
typedef std::map< TFCSParametrizationBase *, Duplicate_t >	FindDuplicates_t
typedef std::map< std::string, FindDuplicates_t >	FindDuplicateClasses_t

Public Member Functions
	TFCSEnergyInterpolationHistogram (const char *name=nullptr, const char *title=nullptr)
bool	OnlyScaleEnergy () const
void	set_OnlyScaleEnergy ()
void	reset_OnlyScaleEnergy ()
virtual bool	is_match_Ekin_bin (int) const override
virtual bool	is_match_calosample (int) const override
void	InitFromHist (const TH1F &hist)
	Initialize interpolation from histogram x values should be Ekin, y values should <E(reco)/Ekin(true)>
const TH1F &	hist () const
virtual FCSReturnCode	simulate (TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const override
	Initialize simulstate with the mean reconstructed energy in the calorimater expeted from the true kinetic energy.
void	Print (Option_t *option="") const override
void	clear ()
virtual bool	is_match_pdgid (int id) const override
virtual bool	is_match_Ekin (float Ekin) const override
virtual bool	is_match_eta (float eta) const override
virtual bool	is_match_all_Ekin () const override
virtual bool	is_match_all_eta () const override
virtual bool	is_match_all_Ekin_bin () const override
virtual bool	is_match_all_calosample () const override
const std::set< int > &	pdgid () const override
double	Ekin_nominal () const override
double	Ekin_min () const override
double	Ekin_max () const override
double	eta_nominal () const override
double	eta_min () const override
double	eta_max () const override
virtual void	set_pdgid (int id)
virtual void	set_pdgid (const std::set< int > &ids)
virtual void	add_pdgid (int id)
virtual void	clear_pdgid ()
virtual void	set_Ekin_nominal (double min)
virtual void	set_Ekin_min (double min)
virtual void	set_Ekin_max (double max)
virtual void	set_eta_nominal (double min)
virtual void	set_eta_min (double min)
virtual void	set_eta_max (double max)
virtual void	set_Ekin (const TFCSParametrizationBase &ref)
virtual void	set_eta (const TFCSParametrizationBase &ref)
virtual void	set_Ekin_eta (const TFCSParametrizationBase &ref)
virtual void	set_pdgid_Ekin_eta (const TFCSParametrizationBase &ref)
virtual bool	is_match_all_pdgid () const
virtual void	set_match_all_pdgid ()
virtual void	reset_match_all_pdgid ()
virtual void	set_geometry (ICaloGeometry *geo)
	Method to set the geometry access pointer.
virtual unsigned int	size () const
	Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.
virtual const TFCSParametrizationBase *	operator[] (unsigned int) const
	Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.
virtual TFCSParametrizationBase *	operator[] (unsigned int)
	Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.
virtual void	set_daughter (unsigned int, TFCSParametrizationBase *)
	Some derived classes have daughter instances of TFCSParametrizationBase objects The set_daughter method allows to change these daughters - expert use only!
virtual bool	operator== (const TFCSParametrizationBase &ref) const
	The == operator compares the content of instances.
virtual void	CleanAuxInfo (TFCSSimulationState &) const
	Method in all derived classes to delete objects stored in the simulstate AuxInfo.
void	FindDuplicates (FindDuplicateClasses_t &dup)
void	RemoveDuplicates ()
void	RemoveNameTitle ()
bool	msgLvl (const MSG::Level lvl) const
	Check whether the logging system is active at the provided verbosity level.
MsgStream &	msg () const
	Return a stream for sending messages directly (no decoration)
MsgStream &	msg (const MSG::Level lvl) const
	Return a decorated starting stream for sending messages.
MSG::Level	level () const
	Retrieve output level.
virtual void	setLevel (MSG::Level lvl)
	Update outputlevel.

Static Public Member Functions
static void	unit_test (TFCSSimulationState *simulstate=nullptr, TFCSTruthState *truth=nullptr, const TFCSExtrapolationState *extrapol=nullptr, TH1F *hist=nullptr)
static std::string	startMsg (MSG::Level lvl, const std::string &file, int line)
	Make a message to decorate the start of logging.

Protected Member Functions
bool	compare (const TFCSParametrizationBase &ref) const

Static Protected Attributes
static constexpr double	init_Ekin_nominal = 0
static constexpr double	init_Ekin_min = 0
	Do not persistify!
static constexpr double	init_Ekin_max = 14000000
	Do not persistify!
static constexpr double	init_eta_nominal = 0
	Do not persistify!
static constexpr double	init_eta_min = -100
	Do not persistify!
static constexpr double	init_eta_max = 100
	Do not persistify!

Private Attributes
TH1F	m_hist
std::set< int >	m_pdgid
double	m_Ekin_nominal
double	m_Ekin_min
double	m_Ekin_max
double	m_eta_nominal
double	m_eta_min
double	m_eta_max
std::string	m_nm
	Message source name.

Static Private Attributes
static boost::thread_specific_ptr< MsgStream > m_msg_tls	ATLAS_THREAD_SAFE
	Do not persistify!

Detailed Description

Definition at line 11 of file TFCSEnergyInterpolationHistogram.h.

Member Typedef Documentation

FindDuplicateClasses_t

typedef std::map<std::string, FindDuplicates_t> TFCSParametrizationBase::FindDuplicateClasses_t

inherited

Definition at line 142 of file TFCSParametrizationBase.h.

FindDuplicates_t

typedef std::map<TFCSParametrizationBase *, Duplicate_t> TFCSParametrizationBase::FindDuplicates_t

inherited

Definition at line 141 of file TFCSParametrizationBase.h.

Member Enumeration Documentation

FCSEnergyInitializationStatusBits

enum TFCSEnergyInterpolationHistogram::FCSEnergyInitializationStatusBits

Enumerator
kOnlyScaleEnergy	Set this bit in the TObject bit field the simulated energy. should only be scaled by the spline

Definition at line 16 of file TFCSEnergyInterpolationHistogram.h.

                                         {
    kOnlyScaleEnergy =
        BIT(15) 
  };

FCSStatusBits

enum TFCSParametrizationBase::FCSStatusBits

inherited

Status bit for FCS needs.

Enumerator
kMatchAllPDGID	Set this bit in the TObject bit field if valid for all PDGID.

Definition at line 52 of file TFCSParametrizationBase.h.

                     {
    kMatchAllPDGID = BIT(
        14) 
  };

Constructor & Destructor Documentation

TFCSEnergyInterpolationHistogram()

TFCSEnergyInterpolationHistogram::TFCSEnergyInterpolationHistogram	(	const char *	name = nullptr,
		const char *	title = nullptr )

Definition at line 32 of file TFCSEnergyInterpolationHistogram.cxx.

    : TFCSParametrization(name, title) {}

Member Function Documentation

add_pdgid()

void TFCSParametrization::add_pdgid ( int id )

virtualinherited

Reimplemented in TFCSEnergyBinParametrization.

Definition at line 35 of file TFCSParametrization.cxx.

{ m_pdgid.insert(id); }

CleanAuxInfo()

virtual void TFCSParametrizationBase::CleanAuxInfo ( TFCSSimulationState & ) const

inlinevirtualinherited

Method in all derived classes to delete objects stored in the simulstate AuxInfo.

Definition at line 130 of file TFCSParametrizationBase.h.

{};

clear()

void TFCSParametrization::clear ( )

inherited

Definition at line 18 of file TFCSParametrization.cxx.

                                {
  m_pdgid.clear();
  m_Ekin_nominal = init_Ekin_nominal;
  m_Ekin_min = init_Ekin_min;
  m_Ekin_max = init_Ekin_max;
  m_eta_nominal = init_eta_nominal;
  m_eta_min = init_eta_min;
  m_eta_max = init_eta_max;
}

clear_pdgid()

void TFCSParametrization::clear_pdgid ( )

virtualinherited

Reimplemented in TFCSEnergyBinParametrization.

Definition at line 37 of file TFCSParametrization.cxx.

{ m_pdgid.clear(); }

compare()

bool TFCSParametrization::compare ( const TFCSParametrizationBase & ref ) const

protectedinherited

Definition at line 78 of file TFCSParametrization.cxx.

                                                                          {
  if (IsA() != ref.IsA()) {
    ATH_MSG_DEBUG("compare(): different class types "
                  << IsA()->GetName() << " != " << ref.IsA()->GetName());
    return false;
  }
  std::string name(GetName());
  if (name != ref.GetName()) {
    ATH_MSG_DEBUG("compare(): different names " << GetName()
                                                << " != " << ref.GetName());
    return false;
  }
  std::string title(GetTitle());
  if (title != ref.GetTitle()) {
    ATH_MSG_DEBUG("compare(): different titles " << GetTitle()
                                                 << " != " << ref.GetTitle());
    return false;
  }
  if (is_match_all_pdgid() != ref.is_match_all_pdgid() ||
      pdgid() != ref.pdgid()) {
    ATH_MSG_DEBUG("compare(): different pdgids");
    return false;
  }
  if (Ekin_nominal() != ref.Ekin_nominal() || Ekin_min() != ref.Ekin_min() ||
      Ekin_max() != ref.Ekin_max()) {
    ATH_MSG_DEBUG("compare(): different Ekin range");
    return false;
  }
  if (eta_nominal() != ref.eta_nominal() || eta_min() != ref.eta_min() ||
      eta_max() != ref.eta_max()) {
    ATH_MSG_DEBUG("compare(): different eta range");
    return false;
  }
 
  return true;
}

Ekin_max()

double TFCSParametrization::Ekin_max ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 37 of file TFCSParametrization.h.

{ return m_Ekin_max; };

Ekin_min()

double TFCSParametrization::Ekin_min ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 36 of file TFCSParametrization.h.

{ return m_Ekin_min; };

Ekin_nominal()

double TFCSParametrization::Ekin_nominal ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 35 of file TFCSParametrization.h.

{ return m_Ekin_nominal; };

eta_max()

double TFCSParametrization::eta_max ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 40 of file TFCSParametrization.h.

{ return m_eta_max; };

eta_min()

double TFCSParametrization::eta_min ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 39 of file TFCSParametrization.h.

{ return m_eta_min; };

eta_nominal()

double TFCSParametrization::eta_nominal ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 38 of file TFCSParametrization.h.

{ return m_eta_nominal; };

FindDuplicates()

void TFCSParametrizationBase::FindDuplicates ( FindDuplicateClasses_t & dup )

inherited

Definition at line 90 of file TFCSParametrizationBase.cxx.

                                        {
 
  for (unsigned int i = 0; i < size(); ++i)
    if ((*this)[i]) {
      TFCSParametrizationBase *param = (*this)[i];
      FindDuplicates_t &dup = dupclasses[param->GetName()];
      // If param is already in the duplication list, skip over
      auto checkexist = dup.find(param);
      if (checkexist != dup.end()) {
        ATH_MSG_WARNING(" [TFCSParametrizationBase::FindDuplicates] "
                        "DUPLICATE POINTER DETECTED");
 
        ATH_MSG_WARNING(" - Pointer           : " << param);
        ATH_MSG_WARNING(" - Name              : " << param->GetName());
        ATH_MSG_WARNING(" - Class             : " << param->ClassName());
        ATH_MSG_WARNING(" - Occurs in parent  : " << this
                          << " (" << this->ClassName() << ")"
                          << " index=" << i);
 
        if (checkexist->second.replace) {
          TFCSParametrizationBase *refparam = checkexist->second.replace;
          ATH_MSG_DEBUG("Found duplicate pointer: "
                        << refparam << "=" << refparam->GetName()
                        << ", duplicate is " << param << "=" << param->GetName()
                        << " index " << i << " of " << this);
          dup[refparam].mother.push_back(this);
          dup[refparam].index.push_back(i);
        }
        continue;
      }
      // Add param to duplication list
      dup[param] = Duplicate_t();
      for (auto &ref : dup) {
        TFCSParametrizationBase *refparam = ref.first;
        // skip itself, as this just added above
        if (param == refparam)
          continue;
        // skip nullptr reference
        if (refparam == nullptr)
          continue;
        // skip reference that is itself going to get replaced
        if (ref.second.replace)
          continue;
        // Check for objects with identical content
        if (*param == *refparam) {
          ATH_MSG_DEBUG("Found duplicate: "
                       << refparam << "=" << refparam->GetName()
                        << ", duplicate is " << param << "=" << param->GetName()
                        << " index " << i << " of " << this);
 
          dup[param].replace = refparam;
          dup[refparam].mother.push_back(this);
          dup[refparam].index.push_back(i);
          break;
        }
      }
      // Continue for child objects in param
      param->FindDuplicates(dupclasses);
    }
}

hist()

const TH1F & TFCSEnergyInterpolationHistogram::hist ( ) const

inline

Definition at line 37 of file TFCSEnergyInterpolationHistogram.h.

{ return m_hist; };

InitFromHist()

void TFCSEnergyInterpolationHistogram::InitFromHist ( const TH1F & hist )

inline

Initialize interpolation from histogram x values should be Ekin, y values should <E(reco)/Ekin(true)>

Definition at line 35 of file TFCSEnergyInterpolationHistogram.h.

{ m_hist = hist; };

is_match_all_calosample()

virtual bool TFCSParametrization::is_match_all_calosample ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Reimplemented in TFCSEnergyAndHitGAN, TFCSEnergyAndHitGANV2, TFCSInitWithEkin, TFCSLateralShapeParametrization, and TFCSPCAEnergyParametrization.

Definition at line 32 of file TFCSParametrization.h.

{ return true; };

is_match_all_Ekin()

virtual bool TFCSParametrization::is_match_all_Ekin ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 25 of file TFCSParametrization.h.

                                                  {
    return Ekin_min() == init_Ekin_min && Ekin_max() == init_Ekin_max;
  };

is_match_all_Ekin_bin()

virtual bool TFCSParametrization::is_match_all_Ekin_bin ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Reimplemented in TFCSEnergyAndHitGAN, TFCSEnergyAndHitGANV2, TFCSInitWithEkin, TFCSLateralShapeParametrization, and TFCSPCAEnergyParametrization.

Definition at line 31 of file TFCSParametrization.h.

{ return true; };

is_match_all_eta()

virtual bool TFCSParametrization::is_match_all_eta ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 28 of file TFCSParametrization.h.

                                                 {
    return eta_min() == init_eta_min && eta_max() == init_eta_max;
  };

is_match_all_pdgid()

virtual bool TFCSParametrizationBase::is_match_all_pdgid ( ) const

inlinevirtualinherited

Definition at line 66 of file TFCSParametrizationBase.h.

{ return TestBit(kMatchAllPDGID); };

is_match_calosample()

virtual bool TFCSEnergyInterpolationHistogram::is_match_calosample ( int ) const

inlineoverridevirtual

Reimplemented from TFCSParametrizationBase.

Definition at line 29 of file TFCSEnergyInterpolationHistogram.h.

                                                         {
    return true;
  };

is_match_Ekin()

virtual bool TFCSParametrization::is_match_Ekin ( float Ekin ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 18 of file TFCSParametrization.h.

                                                        {
    return (Ekin >= m_Ekin_min) && (Ekin < m_Ekin_max);
  };

is_match_Ekin_bin()

virtual bool TFCSEnergyInterpolationHistogram::is_match_Ekin_bin ( int ) const

inlineoverridevirtual

Reimplemented from TFCSParametrizationBase.

Definition at line 26 of file TFCSEnergyInterpolationHistogram.h.

                                                       {
    return true;
  };

is_match_eta()

virtual bool TFCSParametrization::is_match_eta ( float eta ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 21 of file TFCSParametrization.h.

                                                      {
    return (eta >= m_eta_min) && (eta < m_eta_max);
  };

is_match_pdgid()

virtual bool TFCSParametrization::is_match_pdgid ( int id ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 15 of file TFCSParametrization.h.

                                                     {
    return TestBit(kMatchAllPDGID) || m_pdgid.find(id) != m_pdgid.end();
  };

level()

MSG::Level ISF_FCS::MLogging::level ( ) const

inlineinherited

Retrieve output level.

Definition at line 201 of file MLogging.h.

{ return msg().level(); }

msg() [1/2]

MsgStream & ISF_FCS::MLogging::msg ( ) const

inlineinherited

Return a stream for sending messages directly (no decoration)

Definition at line 231 of file MLogging.h.

                                      {
  MsgStream *ms = m_msg_tls.get();
  if (!ms) {
    ms = new MsgStream(Athena::getMessageSvc(), m_nm);
    m_msg_tls.reset(ms);
  }
  return *ms;
}

msg() [2/2]

MsgStream & ISF_FCS::MLogging::msg ( const MSG::Level lvl ) const

inlineinherited

Return a decorated starting stream for sending messages.

Definition at line 240 of file MLogging.h.

                                                        {
  return msg() << lvl;
}

msgLvl()

bool ISF_FCS::MLogging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Check whether the logging system is active at the provided verbosity level.

Definition at line 222 of file MLogging.h.

                                                     {
  if (msg().level() <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

OnlyScaleEnergy()

bool TFCSEnergyInterpolationHistogram::OnlyScaleEnergy ( ) const

inline

Definition at line 22 of file TFCSEnergyInterpolationHistogram.h.

{ return TestBit(kOnlyScaleEnergy); };

operator==()

virtual bool TFCSParametrizationBase::operator== ( const TFCSParametrizationBase & ref ) const

inlinevirtualinherited

The == operator compares the content of instances.

The implementation in the base class only returns true for a comparison with itself

Reimplemented in TFCSHitCellMapping, TFCSHitCellMappingWiggle, TFCSLateralShapeParametrizationHitNumberFromE, and TFCSPredictExtrapWeights.

Definition at line 119 of file TFCSParametrizationBase.h.

                                                                    {
    return compare(ref);
  };

operator[]() [1/2]

virtual TFCSParametrizationBase * TFCSParametrizationBase::operator[] ( unsigned int )

inlinevirtualinherited

Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 106 of file TFCSParametrizationBase.h.

                                                              {
    return nullptr;
  };

operator[]() [2/2]

virtual const TFCSParametrizationBase * TFCSParametrizationBase::operator[] ( unsigned int ) const

inlinevirtualinherited

Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 99 of file TFCSParametrizationBase.h.

                                   {
    return nullptr;
  };

pdgid()

const std::set< int > & TFCSParametrization::pdgid ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 34 of file TFCSParametrization.h.

{ return m_pdgid; };

Print()

void TFCSEnergyInterpolationHistogram::Print ( Option_t * option = "" ) const

override

Definition at line 67 of file TFCSEnergyInterpolationHistogram.cxx.

                                                                   {
  TString opt(option);
  bool shortprint = opt.Index("short") >= 0;
  bool longprint = msgLvl(MSG::DEBUG) || (msgLvl(MSG::INFO) && !shortprint);
  TString optprint = opt;
  optprint.ReplaceAll("short", "");
  TFCSParametrization::Print(option);
 
  if (longprint)
    ATH_MSG_INFO(
        optprint << (OnlyScaleEnergy() ? "  E()*" : "  Ekin()*")
                 << "histNbins=" << m_hist.GetNbinsX() << " "
                 << m_hist.GetXaxis()->GetBinLowEdge(1) << "<=Ekin<="
                 << m_hist.GetXaxis()->GetBinUpEdge(m_hist.GetNbinsX()));
}

RemoveDuplicates()

void TFCSParametrizationBase::RemoveDuplicates ( )

inherited

Definition at line 152 of file TFCSParametrizationBase.cxx.

                                               {
  FindDuplicateClasses_t dupclasses;
  FindDuplicates(dupclasses);
 
  std::set<TFCSParametrizationBase *> dellist;
  for (auto &dupiter : dupclasses) {
    FindDuplicates_t &dup = dupiter.second;
    for (auto onedup : dup) {
      if (onedup.second.mother.empty())
        continue;
      TFCSParametrizationBase *ref = onedup.first;
      ATH_MSG_DEBUG("Main object " << ref << "=" << ref->GetName());
      for (unsigned int i = 0; i < onedup.second.mother.size(); ++i) {
        int index = onedup.second.index[i];
        TFCSParametrizationBase *mother = onedup.second.mother[i];
        TFCSParametrizationBase *delparam = mother->operator[](index);
        unsigned int delcount = dup[delparam].mother.size();
        if (delcount == 0) {
          ATH_MSG_DEBUG("  - Delete object "
                        << delparam << "=" << delparam->GetName() << " index "
                        << index << " of " << mother << ", has " << delcount
                        << " other replacements attached. Deleting");
          mother->set_daughter(index, ref);
          dellist.insert(delparam);
        } else {
          ATH_MSG_WARNING("  - Delete object "
                          << delparam << "=" << delparam->GetName() << " index "
                          << index << " of " << mother << ", has " << delcount
                          << " other replacements attached. Skipping");
        }
      }
    }
  }
 
  ATH_MSG_INFO("RERUNNING DUPLICATE FINDING");
  FindDuplicateClasses_t dupclasses2;
  FindDuplicates(dupclasses2);
 
  std::map<std::string, int> ndel;
  for (auto *delparam : dellist) {
    FindDuplicates_t &dup2 = dupclasses2[delparam->GetName()];
    bool present = dup2.find(delparam) != dup2.end();
    if (present) {
      ATH_MSG_WARNING("- Delete object " << delparam << "="
                                         << delparam->GetName()
                                         << " still referenced somewhere!");
    } else {
      ATH_MSG_DEBUG("- Delete object " << delparam << "="
                                       << delparam->GetName());
      ++ndel[delparam->ClassName()];
      delete delparam;
    }
  }
  for (auto &del : ndel)
    ATH_MSG_INFO("Deleted " << del.second << " duplicate objects of class "
                            << del.first);
}

RemoveNameTitle()

void TFCSParametrizationBase::RemoveNameTitle ( )

inherited

Definition at line 210 of file TFCSParametrizationBase.cxx.

                                              {
  for (unsigned int i = 0; i < size(); ++i)
    if ((*this)[i]) {
      TFCSParametrizationBase *param = (*this)[i];
      param->SetName("");
      param->SetTitle("");
 
      // Continue for child objects in param
      param->RemoveNameTitle();
    }
}

reset_match_all_pdgid()

virtual void TFCSParametrizationBase::reset_match_all_pdgid ( )

inlinevirtualinherited

Definition at line 84 of file TFCSParametrizationBase.h.

{ ResetBit(kMatchAllPDGID); };

reset_OnlyScaleEnergy()

void TFCSEnergyInterpolationHistogram::reset_OnlyScaleEnergy ( )

inline

Definition at line 24 of file TFCSEnergyInterpolationHistogram.h.

{ ResetBit(kOnlyScaleEnergy); };

set_daughter()

virtual void TFCSParametrizationBase::set_daughter ( unsigned int , TFCSParametrizationBase *  )

inlinevirtualinherited

Some derived classes have daughter instances of TFCSParametrizationBase objects The set_daughter method allows to change these daughters - expert use only!

The original element at this position is not deleted

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 113 of file TFCSParametrizationBase.h.

                                                        {};

set_Ekin()

void TFCSParametrization::set_Ekin ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 55 of file TFCSParametrization.cxx.

                                                                     {
  set_Ekin_nominal(ref.Ekin_nominal());
  set_Ekin_min(ref.Ekin_min());
  set_Ekin_max(ref.Ekin_max());
}

set_Ekin_eta()

void TFCSParametrization::set_Ekin_eta ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 67 of file TFCSParametrization.cxx.

                                                                         {
  set_Ekin(ref);
  set_eta(ref);
}

set_Ekin_max()

void TFCSParametrization::set_Ekin_max ( double max )

virtualinherited

Definition at line 45 of file TFCSParametrization.cxx.

{ m_Ekin_max = max; }

set_Ekin_min()

void TFCSParametrization::set_Ekin_min ( double min )

virtualinherited

Definition at line 43 of file TFCSParametrization.cxx.

{ m_Ekin_min = min; }

set_Ekin_nominal()

void TFCSParametrization::set_Ekin_nominal ( double min )

virtualinherited

Definition at line 39 of file TFCSParametrization.cxx.

                                                         {
  m_Ekin_nominal = nominal;
}

set_eta()

void TFCSParametrization::set_eta ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 61 of file TFCSParametrization.cxx.

                                                                    {
  set_eta_nominal(ref.eta_nominal());
  set_eta_min(ref.eta_min());
  set_eta_max(ref.eta_max());
}

set_eta_max()

void TFCSParametrization::set_eta_max ( double max )

virtualinherited

Definition at line 53 of file TFCSParametrization.cxx.

{ m_eta_max = max; }

set_eta_min()

void TFCSParametrization::set_eta_min ( double min )

virtualinherited

Definition at line 51 of file TFCSParametrization.cxx.

{ m_eta_min = min; }

set_eta_nominal()

void TFCSParametrization::set_eta_nominal ( double min )

virtualinherited

Definition at line 47 of file TFCSParametrization.cxx.

                                                        {
  m_eta_nominal = nominal;
}

set_geometry()

void TFCSParametrizationBase::set_geometry ( ICaloGeometry * geo )

virtualinherited

Method to set the geometry access pointer.

Loops over daughter objects if present

Reimplemented in TFCSBinnedShowerBase, TFCSHistoLateralShapeParametrization, TFCSHitCellMapping, TFCSPhiModulationCorrection, and TFCSVoxelHistoLateralCovarianceFluctuations.

Definition at line 24 of file TFCSParametrizationBase.cxx.

                                                             {
  for (unsigned int i = 0; i < size(); ++i)
    (*this)[i]->set_geometry(geo);
}

set_match_all_pdgid()

virtual void TFCSParametrizationBase::set_match_all_pdgid ( )

inlinevirtualinherited

Definition at line 83 of file TFCSParametrizationBase.h.

{ SetBit(kMatchAllPDGID); };

set_OnlyScaleEnergy()

void TFCSEnergyInterpolationHistogram::set_OnlyScaleEnergy ( )

inline

Definition at line 23 of file TFCSEnergyInterpolationHistogram.h.

{ SetBit(kOnlyScaleEnergy); };

set_pdgid() [1/2]

void TFCSParametrization::set_pdgid ( const std::set< int > & ids )

virtualinherited

Reimplemented in TFCSEnergyBinParametrization.

Definition at line 33 of file TFCSParametrization.cxx.

{ m_pdgid = ids; }

set_pdgid() [2/2]

void TFCSParametrization::set_pdgid ( int id )

virtualinherited

Reimplemented in TFCSEnergyBinParametrization.

Definition at line 28 of file TFCSParametrization.cxx.

                                          {
  m_pdgid.clear();
  m_pdgid.insert(id);
}

set_pdgid_Ekin_eta()

void TFCSParametrization::set_pdgid_Ekin_eta ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 72 of file TFCSParametrization.cxx.

                                        {
  set_Ekin_eta(ref);
  set_pdgid(ref.pdgid());
}

setLevel()

void ISF_FCS::MLogging::setLevel ( MSG::Level lvl )

virtualinherited

Update outputlevel.

Definition at line 105 of file MLogging.cxx.

                                    {
  lvl = (lvl >= MSG::NUM_LEVELS) ? MSG::ALWAYS
        : (lvl < MSG::NIL)       ? MSG::NIL
                                 : lvl;
  msg().setLevel(lvl);
}

simulate()

FCSReturnCode TFCSEnergyInterpolationHistogram::simulate ( TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol ) const

overridevirtual

Initialize simulstate with the mean reconstructed energy in the calorimater expeted from the true kinetic energy.

Reimplemented from TFCSParametrizationBase.

Definition at line 36 of file TFCSEnergyInterpolationHistogram.cxx.

                                          {
  float Emean;
  float Einit;
  const float Ekin = truth->Ekin();
  if (OnlyScaleEnergy())
    Einit = simulstate.E();
  else
    Einit = Ekin;
  if (Einit < m_hist.GetXaxis()->GetBinLowEdge(1)) {
    Emean = m_hist.GetBinContent(1) * Einit;
  } else {
    if (Einit > m_hist.GetXaxis()->GetBinUpEdge(m_hist.GetNbinsX())) {
      Emean = m_hist.GetBinContent(m_hist.GetNbinsX()) * Einit;
    } else {
      Emean = m_hist.GetBinContent(m_hist.GetXaxis()->FindBin(Einit)) * Einit;
    }
  }
 
  if (OnlyScaleEnergy()) {
    ATH_MSG_DEBUG("set E=" << Emean << " for true Ekin=" << truth->Ekin()
                           << " and E=" << Einit);
  } else {
    ATH_MSG_DEBUG("set E=" << Emean << " for true Ekin=" << truth->Ekin());
  }
  simulstate.set_E(Emean);
 
  return FCSSuccess;
}

size()

virtual unsigned int TFCSParametrizationBase::size ( ) const

inlinevirtualinherited

Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 93 of file TFCSParametrizationBase.h.

{ return 0; };

startMsg()

std::string ISF_FCS::MLogging::startMsg ( MSG::Level lvl, const std::string & file, int line )

staticinherited

Make a message to decorate the start of logging.

Print a message for the start of logging.

Definition at line 116 of file MLogging.cxx.

                                                                          {
  int col1_len = 20;
  int col2_len = 5;
  int col3_len = 10;
  auto last_slash = file.find_last_of('/');
  int path_len = last_slash == std::string::npos ? 0 : last_slash;
  int trim_point = path_len;
  int total_len = file.length();
  if (total_len - path_len > col1_len)
    trim_point = total_len - col1_len;
  std::string trimmed_name = file.substr(trim_point);
  const char *LevelNames[MSG::NUM_LEVELS] = {
      "NIL", "VERBOSE", "DEBUG", "INFO", "WARNING", "ERROR", "FATAL", "ALWAYS"};
  std::string level = LevelNames[lvl];
  std::string level_string = std::string("(") + level + ") ";
  std::stringstream output;
  output << std::setw(col1_len) << std::right << trimmed_name << ":"
         << std::setw(col2_len) << std::left << line << std::setw(col3_len)
         << std::right << level_string;
  return output.str();
}

unit_test()

void TFCSEnergyInterpolationHistogram::unit_test ( TFCSSimulationState * simulstate = nullptr, TFCSTruthState * truth = nullptr, const TFCSExtrapolationState * extrapol = nullptr, TH1F * hist = nullptr )

static

Definition at line 83 of file TFCSEnergyInterpolationHistogram.cxx.

                                                        {
  if (!simulstate)
    simulstate = new TFCSSimulationState();
  if (!truth)
    truth = new TFCSTruthState();
  if (!extrapol)
    extrapol = new TFCSExtrapolationState();
 
  if (!hist) {
    hist = new TH1F("h1", "h1 title", 16, 0., 512.);
    hist->SetBinContent(1, 0.687165);
    hist->SetBinContent(2, 0.756837);
    hist->SetBinContent(3, 0.836673);
    hist->SetBinContent(4, 0.896336);
    hist->SetBinContent(5, 0.889785);
    hist->SetBinContent(6, 0.901266);
    hist->SetBinContent(7, 0.888937);
    hist->SetBinContent(8, 0.919943);
    hist->SetBinContent(9, 0.941806);
    hist->SetBinContent(10, 0.934668);
    hist->SetBinContent(11, 0.939502);
    hist->SetBinContent(12, 0.940796);
    hist->SetBinContent(13, 0.945894);
    hist->SetBinContent(14, 0.955445);
    hist->SetBinContent(15, 0.955593);
    hist->SetBinContent(16, 0.943673);
  }
 
  TH1F *hdraw = (TH1F *)hist->Clone();
  hdraw->SetMarkerColor(46);
  hdraw->SetMarkerStyle(8);
 
  TFCSEnergyInterpolationHistogram test(
      "testTFCSEnergyInterpolationHistogram",
      "test TFCSEnergyInterpolationHistogram");
  test.set_pdgid(22);
  test.set_Ekin_nominal(0.5 *
                        (hdraw->GetXaxis()->GetBinLowEdge(1) +
                         hdraw->GetXaxis()->GetBinUpEdge(hdraw->GetNbinsX())));
  test.set_Ekin_min(hdraw->GetXaxis()->GetBinLowEdge(1));
  test.set_Ekin_max(hdraw->GetXaxis()->GetBinUpEdge(hdraw->GetNbinsX()));
  test.set_eta_nominal(0.225);
  test.set_eta_min(0.2);
  test.set_eta_max(0.25);
  test.InitFromHist(*hist);
  // test.set_OnlyScaleEnergy();
  test.Print();
  test.hist().Dump();
 
  truth->set_pdgid(22);
 
  TGraph *gr = new TGraph();
  gr->SetNameTitle("testTFCSEnergyInterpolationHistogramLogX",
                   "test TFCSEnergyInterpolationHistogram");
  gr->GetXaxis()->SetTitle("Ekin [MeV]");
  gr->GetYaxis()->SetTitle("<E(reco)>/Ekin(true)");
 
  int ip = 0;
  for (float Ekin = std::max(test.Ekin_min() * 0.25, 0.1);
       Ekin <= test.Ekin_max() * 4; Ekin *= 1.05) {
    // Init LorentzVector for truth. For photon Ekin=E
    truth->SetPxPyPzE(Ekin, 0, 0, Ekin);
    simulstate->set_E(Ekin);
    if (test.simulate(*simulstate, truth, extrapol) != FCSSuccess) {
      return;
    }
    gr->SetPoint(ip, Ekin, simulstate->E() / Ekin);
    ++ip;
  }
 
// Drawing doesn't make sense inside athena and necessary libraries not linked
// by default
#if defined(__FastCaloSimStandAlone__)
  TCanvas *c = new TCanvas(hdraw->GetName(), hdraw->GetTitle());
  hdraw->Draw("HIST");
  gr->Draw("same,APL");
  c->SetLogx();
#endif
}

Member Data Documentation

ATLAS_THREAD_SAFE

boost::thread_specific_ptr<MsgStream> m_msg_tls ISF_FCS::MLogging::ATLAS_THREAD_SAFE

inlinestaticprivateinherited

Do not persistify!

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

Definition at line 215 of file MLogging.h.

init_Ekin_max

double TFCSParametrizationBase::init_Ekin_max = 14000000

staticconstexprprotectedinherited

Do not persistify!

Definition at line 155 of file TFCSParametrizationBase.h.

init_Ekin_min

double TFCSParametrizationBase::init_Ekin_min = 0

staticconstexprprotectedinherited

Do not persistify!

Definition at line 154 of file TFCSParametrizationBase.h.

init_Ekin_nominal

double TFCSParametrizationBase::init_Ekin_nominal = 0

staticconstexprprotectedinherited

Definition at line 153 of file TFCSParametrizationBase.h.

init_eta_max

double TFCSParametrizationBase::init_eta_max = 100

staticconstexprprotectedinherited

Do not persistify!

Definition at line 158 of file TFCSParametrizationBase.h.

init_eta_min

double TFCSParametrizationBase::init_eta_min = -100

staticconstexprprotectedinherited

Do not persistify!

Definition at line 157 of file TFCSParametrizationBase.h.

init_eta_nominal

double TFCSParametrizationBase::init_eta_nominal = 0

staticconstexprprotectedinherited

Do not persistify!

Definition at line 156 of file TFCSParametrizationBase.h.

m_Ekin_max

double TFCSParametrization::m_Ekin_max

privateinherited

Definition at line 64 of file TFCSParametrization.h.

m_Ekin_min

double TFCSParametrization::m_Ekin_min

privateinherited

Definition at line 64 of file TFCSParametrization.h.

m_Ekin_nominal

double TFCSParametrization::m_Ekin_nominal

privateinherited

Definition at line 64 of file TFCSParametrization.h.

m_eta_max

double TFCSParametrization::m_eta_max

privateinherited

Definition at line 65 of file TFCSParametrization.h.

m_eta_min

double TFCSParametrization::m_eta_min

privateinherited

Definition at line 65 of file TFCSParametrization.h.

m_eta_nominal

double TFCSParametrization::m_eta_nominal

privateinherited

Definition at line 65 of file TFCSParametrization.h.

m_hist

TH1F TFCSEnergyInterpolationHistogram::m_hist

private

Definition at line 52 of file TFCSEnergyInterpolationHistogram.h.

m_nm

std::string ISF_FCS::MLogging::m_nm

privateinherited

Message source name.

Definition at line 211 of file MLogging.h.

m_pdgid

std::set<int> TFCSParametrization::m_pdgid

privateinherited

Definition at line 63 of file TFCSParametrization.h.

---

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

• TFCSEnergyInterpolationHistogram.h
• TFCSEnergyInterpolationHistogram.cxx