#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)> More...

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. More...

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. More...

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. More...

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. More...

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. More...

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! The original element at this position is not deleted. More...

virtual bool	operator== (const TFCSParametrizationBase &ref) const
	The == operator compares the content of instances. More...

virtual void	CleanAuxInfo (TFCSSimulationState &) const
	Method in all derived classes to delete objects stored in the simulstate AuxInfo. More...

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. More...

MsgStream &	msg () const
	Return a stream for sending messages directly (no decoration) More...

MsgStream &	msg (const MSG::Level lvl) const
	Return a decorated starting stream for sending messages. More...

MSG::Level	level () const
	Retrieve output level. More...

virtual void	setLevel (MSG::Level lvl)
	Update outputlevel. More...

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. More...

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! More...

static constexpr double	init_Ekin_max = 14000000
	Do not persistify! More...

static constexpr double	init_eta_nominal = 0
	Do not persistify! More...

static constexpr double	init_eta_min = -100
	Do not persistify! More...

static constexpr double	init_eta_max = 100
	Do not persistify! More...

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. More...

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

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.

34 : 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.

35 { 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.

130 {};

◆ 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.

37 { 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.

37 { return m_Ekin_max; };

◆ Ekin_min()

double TFCSParametrization::Ekin_min ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 36 of file TFCSParametrization.h.

36 { return m_Ekin_min; };

◆ Ekin_nominal()

double TFCSParametrization::Ekin_nominal ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 35 of file TFCSParametrization.h.

35 { return m_Ekin_nominal; };

◆ eta_max()

double TFCSParametrization::eta_max ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 40 of file TFCSParametrization.h.

40 { return m_eta_max; };

◆ eta_min()

double TFCSParametrization::eta_min ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 39 of file TFCSParametrization.h.

39 { return m_eta_min; };

◆ eta_nominal()

double TFCSParametrization::eta_nominal ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 38 of file TFCSParametrization.h.

38 { 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_DEBUG("Found duplicate pointer for: " << param << "="
                                                       << param->GetName());
         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.

37 { 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.

35 { m_hist = hist; };

◆ is_match_all_calosample()

virtual bool TFCSParametrization::is_match_all_calosample ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

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

Definition at line 32 of file TFCSParametrization.h.

32 { 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 TFCSPCAEnergyParametrization, TFCSEnergyAndHitGAN, TFCSEnergyAndHitGANV2, TFCSLateralShapeParametrization, and TFCSInitWithEkin.

Definition at line 31 of file TFCSParametrization.h.

31 { 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.

66 { 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.

201 { 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.

22 { 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 TFCSHitCellMappingWiggle, TFCSLateralShapeParametrizationHitNumberFromE, TFCSHitCellMapping, 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.

34 { 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 142 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 200 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.

84 { ResetBit(kMatchAllPDGID); };

◆ reset_OnlyScaleEnergy()

void TFCSEnergyInterpolationHistogram::reset_OnlyScaleEnergy ( )

inline

Definition at line 24 of file TFCSEnergyInterpolationHistogram.h.

24 { 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 TFCSParametrizationChain, and TFCSLateralShapeParametrizationHitChain.

Definition at line 113 of file TFCSParametrizationBase.h.

114 {};

◆ 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.

45 { m_Ekin_max = max; }

◆ set_Ekin_min()

void TFCSParametrization::set_Ekin_min ( double min )

virtualinherited

Definition at line 43 of file TFCSParametrization.cxx.

43 { 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.

53 { m_eta_max = max; }

◆ set_eta_min()

void TFCSParametrization::set_eta_min ( double min )

virtualinherited

Definition at line 51 of file TFCSParametrization.cxx.

51 { 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 TFCSHistoLateralShapeParametrization, TFCSVoxelHistoLateralCovarianceFluctuations, and TFCSHitCellMapping.

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.

83 { SetBit(kMatchAllPDGID); };

◆ set_OnlyScaleEnergy()

void TFCSEnergyInterpolationHistogram::set_OnlyScaleEnergy ( )

inline

Definition at line 23 of file TFCSEnergyInterpolationHistogram.h.

23 { 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.

33 { 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.

93 { 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

constexpr double TFCSParametrizationBase::init_Ekin_max = 14000000

staticconstexprprotectedinherited

Do not persistify!

Definition at line 155 of file TFCSParametrizationBase.h.

◆ init_Ekin_min

constexpr double TFCSParametrizationBase::init_Ekin_min = 0

staticconstexprprotectedinherited

Do not persistify!

Definition at line 154 of file TFCSParametrizationBase.h.

◆ init_Ekin_nominal

constexpr double TFCSParametrizationBase::init_Ekin_nominal = 0

staticconstexprprotectedinherited

Definition at line 153 of file TFCSParametrizationBase.h.

◆ init_eta_max

constexpr double TFCSParametrizationBase::init_eta_max = 100

staticconstexprprotectedinherited

Do not persistify!

Definition at line 158 of file TFCSParametrizationBase.h.

◆ init_eta_min

constexpr double TFCSParametrizationBase::init_eta_min = -100

staticconstexprprotectedinherited

Do not persistify!

Definition at line 157 of file TFCSParametrizationBase.h.

◆ init_eta_nominal

constexpr 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:

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Static Protected Attributes

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ FindDuplicateClasses_t

◆ FindDuplicates_t

Member Enumeration Documentation

◆ FCSEnergyInitializationStatusBits

◆ FCSStatusBits

Constructor & Destructor Documentation

◆ TFCSEnergyInterpolationHistogram()

Member Function Documentation

◆ add_pdgid()

◆ CleanAuxInfo()

◆ clear()

◆ clear_pdgid()

◆ compare()

◆ Ekin_max()

◆ Ekin_min()

◆ Ekin_nominal()

◆ eta_max()

◆ eta_min()

◆ eta_nominal()

◆ FindDuplicates()

◆ hist()

◆ InitFromHist()

◆ is_match_all_calosample()

◆ is_match_all_Ekin()

◆ is_match_all_Ekin_bin()

◆ is_match_all_eta()

◆ is_match_all_pdgid()

◆ is_match_calosample()

◆ is_match_Ekin()

◆ is_match_Ekin_bin()

◆ is_match_eta()

◆ is_match_pdgid()

◆ level()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ OnlyScaleEnergy()

◆ operator==()

◆ operator[]() [1/2]

◆ operator[]() [2/2]

◆ pdgid()

◆ Print()

◆ RemoveDuplicates()

◆ RemoveNameTitle()

◆ reset_match_all_pdgid()

◆ reset_OnlyScaleEnergy()

◆ set_daughter()

◆ set_Ekin()

◆ set_Ekin_eta()

◆ set_Ekin_max()

◆ set_Ekin_min()

◆ set_Ekin_nominal()

◆ set_eta()

◆ set_eta_max()

◆ set_eta_min()

◆ set_eta_nominal()

◆ set_geometry()

◆ set_match_all_pdgid()

◆ set_OnlyScaleEnergy()

◆ set_pdgid() [1/2]

◆ set_pdgid() [2/2]

◆ set_pdgid_Ekin_eta()

◆ setLevel()

◆ simulate()

◆ size()

◆ startMsg()

◆ unit_test()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ init_Ekin_max

◆ init_Ekin_min