#include <TFCSEnergyBinParametrization.h>

Inheritance diagram for TFCSEnergyBinParametrization:

Collaboration diagram for TFCSEnergyBinParametrization:

Public Types
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
	TFCSEnergyBinParametrization (const char name=nullptr, const char title=nullptr)

virtual void	set_pdgid (int id) override

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

virtual void	add_pdgid (int id) override

virtual void	clear_pdgid () override

virtual int	n_bins () const override

void	set_number_of_Ekin_bins (int n_Ekin_bin)
	current convention is to start Ekin_bin counting at 1, to be updated to start counting with 0 More...

virtual void	set_pdgid_Ekin_bin_probability (int id, std::vector< float > prob)
	set_pdgid_Ekin_bin_probability(): format of probability should be a vector of float with probability values, e.g. More...

virtual bool	load_pdgid_Ekin_bin_probability_from_file (int id, TFile *file, std::string prob_object_name)

virtual FCSReturnCode	simulate (TFCSSimulationState &simulstate, const TFCSTruthState truth, const TFCSExtrapolationState extrapol) const override
	Method in all derived classes to do some simulation. More...

virtual bool	is_match_Ekin_bin (int Ekin_bin) const override

void	Print (Option_t *option="") const override

virtual bool	is_match_calosample (int) 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_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 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 Member Functions
void	resize ()

Private Attributes
int	m_number_of_Ekin_bins

std::map< int, std::vector< float > >	m_pdgid_Ebin_probability

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 13 of file TFCSEnergyBinParametrization.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

◆ 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

◆ TFCSEnergyBinParametrization()

TFCSEnergyBinParametrization::TFCSEnergyBinParametrization	(	const char *	name = `nullptr`,
		const char *	title = `nullptr`
	)

Definition at line 18 of file TFCSEnergyBinParametrization.cxx.

     : TFCSEnergyParametrization(name, title) {
   m_number_of_Ekin_bins = 1;
 }

Member Function Documentation

◆ add_pdgid()

void TFCSEnergyBinParametrization::add_pdgid ( int id )

overridevirtual

Reimplemented from TFCSParametrization.

Definition at line 54 of file TFCSEnergyBinParametrization.cxx.

                                                    {
   TFCSEnergyParametrization::add_pdgid(id);
   resize();
 }

◆ 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 TFCSEnergyBinParametrization::clear_pdgid ( )

overridevirtual

Reimplemented from TFCSParametrization.

Definition at line 59 of file TFCSEnergyBinParametrization.cxx.

                                                {
   TFCSEnergyParametrization::clear_pdgid();
   resize();
 }

◆ 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);
     }
 }

◆ 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 TFCSEnergyParametrization::is_match_calosample ( int ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Reimplemented in TFCSPCAEnergyParametrization.

Definition at line 18 of file TFCSEnergyParametrization.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()

bool TFCSEnergyBinParametrization::is_match_Ekin_bin ( int Ekin_bin ) const

overridevirtual

Reimplemented from TFCSEnergyParametrization.

Definition at line 24 of file TFCSEnergyBinParametrization.cxx.

                                                                        {
   if (Ekin_bin >= 1 && Ekin_bin <= n_bins())
     return true;
   return false;
 }

◆ 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(); }

◆ load_pdgid_Ekin_bin_probability_from_file()

bool TFCSEnergyBinParametrization::load_pdgid_Ekin_bin_probability_from_file	(	int	id,
		TFile *	file,
		std::string	prob_object_name
	)

virtual

Definition at line 84 of file TFCSEnergyBinParametrization.cxx.

                                                      {
   add_pdgid(id);
  
   file->cd();
   auto *probFromFile =
       dynamic_cast<TVectorF *>(gDirectory->Get(prob_object_name.c_str()));
  
   std::vector<float> prob;
   prob.reserve(m_pdgid_Ebin_probability[id].size());
   if (probFromFile == nullptr) {
     ATH_MSG_INFO("TFCSEnergyBinParametrization::load_pdgid_Ekin_bin_"
                  "probability_from_file(): "
                  << prob_object_name
                  << " is null. Using equal PCA probabilities.");
     prob.push_back(0.0);
     for (size_t i = 1; i < m_pdgid_Ebin_probability[id].size(); i++) {
       prob.push_back(1.0);
     }
   } else {
     auto size = static_cast<size_t>(probFromFile->GetNoElements());
     if (size != m_pdgid_Ebin_probability[id].size()) {
       ATH_MSG_ERROR(
           "TFCSEnergyBinParametrization::load_pdgid_Ekin_bin_probability_from_"
           "file(): size of prob array does not match! in.size()="
           << size << " instance=" << m_pdgid_Ebin_probability[id].size());
       return false;
     }
  
     const float *probArray = probFromFile->GetMatrixArray();
     for (size_t i = 0; i < size; i++) {
       prob.push_back(probArray[i]);
     }
   }
  
   float ptot{};
   for (int iEbin = 0; iEbin <= n_bins(); ++iEbin)
     ptot += prob[iEbin];
   float p{};
   for (int iEbin = 0; iEbin <= n_bins(); ++iEbin) {
     p += prob[iEbin] / ptot;
     m_pdgid_Ebin_probability[id][iEbin] = p;
   }
  
   return true;
 }

◆ 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;
   }
 }

◆ n_bins()

virtual int TFCSEnergyBinParametrization::n_bins ( ) const

inlineoverridevirtual

Reimplemented from TFCSEnergyParametrization.

Definition at line 23 of file TFCSEnergyBinParametrization.h.

23 { return m_number_of_Ekin_bins; };

◆ 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 TFCSEnergyBinParametrization::Print ( Option_t * option = "" ) const

override

Definition at line 131 of file TFCSEnergyBinParametrization.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", "");
   TFCSEnergyParametrization::Print(option);
   if (longprint) {
     for (std::set<int>::iterator it = pdgid().begin(); it != pdgid().end();
          ++it) {
       ATH_MSG(INFO) << optprint << "  PDGID=" << *it << " : ";
       float p = 0;
       for (int iEbin = 0; iEbin <= n_bins(); ++iEbin) {
         if (iEbin > 0)
           ATH_MSG(INFO) << ", ";
         auto mapit = m_pdgid_Ebin_probability.find(*it);
         ATH_MSG(INFO) << "b" << iEbin << "="
                       << (mapit->second[iEbin] - p) / mapit->second.back() * 100
                       << "%";
         p = mapit->second[iEbin];
       }
       ATH_MSG(INFO) << END_MSG(INFO);
     }
   }
 }

◆ 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); };

◆ resize()

void TFCSEnergyBinParametrization::resize ( )

private

Definition at line 30 of file TFCSEnergyBinParametrization.cxx.

                                           {
   for (int id : pdgid()) {
     ATH_MSG_VERBOSE("PDGid=" << id << " resize to " << n_bins() + 1);
     m_pdgid_Ebin_probability[id].resize(n_bins() + 1);
   }
   for (auto it = m_pdgid_Ebin_probability.begin();
        it != m_pdgid_Ebin_probability.end(); ++it) {
     if (!is_match_pdgid(it->first)) {
       m_pdgid_Ebin_probability.erase(it);
       it = m_pdgid_Ebin_probability.begin();
     }
   }
 }

◆ 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_number_of_Ekin_bins()

void TFCSEnergyBinParametrization::set_number_of_Ekin_bins ( int n_Ekin_bin )

inline

current convention is to start Ekin_bin counting at 1, to be updated to start counting with 0

Definition at line 27 of file TFCSEnergyBinParametrization.h.

                                                {
     m_number_of_Ekin_bins = n_Ekin_bin;
     resize();
   };

◆ set_pdgid() [1/2]

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

overridevirtual

Reimplemented from TFCSParametrization.

Definition at line 49 of file TFCSEnergyBinParametrization.cxx.

                                                                    {
   TFCSEnergyParametrization::set_pdgid(ids);
   resize();
 }

◆ set_pdgid() [2/2]

void TFCSEnergyBinParametrization::set_pdgid ( int id )

overridevirtual

Reimplemented from TFCSParametrization.

Definition at line 44 of file TFCSEnergyBinParametrization.cxx.

                                                    {
   TFCSEnergyParametrization::set_pdgid(id);
   resize();
 }

◆ set_pdgid_Ekin_bin_probability()

void TFCSEnergyBinParametrization::set_pdgid_Ekin_bin_probability	(	int	id,
		std::vector< float >	prob
	)

virtual

set_pdgid_Ekin_bin_probability(): format of probability should be a vector of float with probability values, e.g.

for equal probability of 5 Ekin bins specify a vector with {0,1,1,1,1,1} or {0,0.2,0.2,0.2,0.2,0.2} This would give a probability of 0 for Ekin_bin=0 and a probability of 20% for 1<=Ekin_bin<=5 the function will normalize probabilities automatically, if the sum of values is not 1 current convention is to start Ekin_bin counting at 1, to be updated to start counting with 0

Definition at line 64 of file TFCSEnergyBinParametrization.cxx.

                                    {
   add_pdgid(id);
   if (prob.size() != m_pdgid_Ebin_probability[id].size()) {
     ATH_MSG_ERROR("TFCSEnergyBinParametrization::set_pdgid_Ekin_bin_"
                   "probability(): size of vectors does not match! in.size()="
                   << prob.size()
                   << " instance=" << m_pdgid_Ebin_probability[id].size());
     return;
   }
   float ptot = 0;
   for (int iEbin = 0; iEbin <= n_bins(); ++iEbin)
     ptot += prob[iEbin];
   float p = 0;
   for (int iEbin = 0; iEbin <= n_bins(); ++iEbin) {
     p += prob[iEbin] / ptot;
     m_pdgid_Ebin_probability[id][iEbin] = p;
   }
 }

◆ 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 TFCSEnergyBinParametrization::simulate	(	TFCSSimulationState &	simulstate,
		const TFCSTruthState *	truth,
		const TFCSExtrapolationState *	extrapol
	)		const

overridevirtual

Method in all derived classes to do some simulation.

Result should be returned in simulstate.

Simulate all energies in calo layers for energy parametrizations. Simulate cells for shape simulation.

Reimplemented from TFCSParametrizationBase.

Definition at line 157 of file TFCSEnergyBinParametrization.cxx.

                                             {
   if (!simulstate.randomEngine()) {
     return FCSFatal;
   }
  
   const int truth_pdgid = truth->pdgid();
   int pdgid = -99;
  
   if (is_match_pdgid(truth_pdgid))
     pdgid = truth_pdgid;
   else if (is_match_pdgid(0))
     pdgid = 0;
   else {
     ATH_MSG_ERROR(
         "TFCSEnergyBinParametrization::simulate(): cannot simulate pdgid="
         << truth_pdgid);
     return FCSFatal;
   }
  
   const float searchRand = CLHEP::RandFlat::shoot(simulstate.randomEngine());
   const auto &Ebin_probability = m_pdgid_Ebin_probability.at(pdgid);
   int chosenBin =
       TMath::BinarySearch(n_bins() + 1, Ebin_probability.data(), searchRand) +
       1;
   if (chosenBin < 0) {
     ATH_MSG_WARNING("TFCSEnergyBinParametrization::simulate(): chosenBin<0 "
                     "(will use chosenBin=0)");
     std::string array = "";
     for (const auto &prob : Ebin_probability) {
       array += prob;
       array += " ";
     }
     ATH_MSG_WARNING(" E=" << simulstate.E() << " Ebin=" << chosenBin
                           << " rnd=" << searchRand << " array=" << array);
     chosenBin = 0;
   } else if (chosenBin > n_bins()) {
     ATH_MSG_WARNING("TFCSEnergyBinParametrization::simulate(): "
                     "chosenBin>n_bins() (will use chosenBin=n_bins())");
     std::string array = "";
     for (const auto &prob : Ebin_probability) {
       array += prob;
       array += " ";
     }
     ATH_MSG_WARNING(" E=" << simulstate.E() << " Ebin=" << chosenBin
                           << " rnd=" << searchRand << " array=" << array);
     chosenBin = n_bins();
   }
   simulstate.set_Ebin(chosenBin);
   ATH_MSG_DEBUG("Ebin=" << chosenBin);
  
   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();
 }

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_nm

std::string ISF_FCS::MLogging::m_nm

privateinherited

Message source name.

Definition at line 211 of file MLogging.h.

◆ m_number_of_Ekin_bins

int TFCSEnergyBinParametrization::m_number_of_Ekin_bins

private

Definition at line 53 of file TFCSEnergyBinParametrization.h.

◆ m_pdgid

std::set<int> TFCSParametrization::m_pdgid

privateinherited

Definition at line 63 of file TFCSParametrization.h.

◆ m_pdgid_Ebin_probability

std::map<int, std::vector<float> > TFCSEnergyBinParametrization::m_pdgid_Ebin_probability

private

Definition at line 54 of file TFCSEnergyBinParametrization.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 Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ FindDuplicateClasses_t

◆ FindDuplicates_t

Member Enumeration Documentation

◆ FCSStatusBits

Constructor & Destructor Documentation

◆ TFCSEnergyBinParametrization()

Member Function Documentation

◆ add_pdgid()

◆ CleanAuxInfo()

◆ clear()

◆ clear_pdgid()

◆ compare()

◆ Ekin_max()

◆ Ekin_min()

◆ Ekin_nominal()

◆ eta_max()

◆ eta_min()

◆ eta_nominal()

◆ FindDuplicates()

◆ 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()

◆ load_pdgid_Ekin_bin_probability_from_file()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ n_bins()

◆ operator==()

◆ operator[]() [1/2]

◆ operator[]() [2/2]

◆ pdgid()

◆ Print()

◆ RemoveDuplicates()

◆ RemoveNameTitle()

◆ reset_match_all_pdgid()

◆ resize()

◆ 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_number_of_Ekin_bins()

◆ set_pdgid() [1/2]

◆ set_pdgid() [2/2]

◆ set_pdgid_Ekin_bin_probability()

◆ set_pdgid_Ekin_eta()

◆ setLevel()

◆ simulate()

◆ size()

◆ startMsg()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ init_Ekin_max

◆ init_Ekin_min

◆ init_Ekin_nominal