TFCSLateralShapeParametrizationHitChain Class Reference

#include <TFCSLateralShapeParametrizationHitChain.h>

Inheritance diagram for TFCSLateralShapeParametrizationHitChain:

Collaboration diagram for TFCSLateralShapeParametrizationHitChain:

Public Types
typedef std::vector< TFCSLateralShapeParametrizationHitBase * >	Chain_t
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
	TFCSLateralShapeParametrizationHitChain (const char *name=nullptr, const char *title=nullptr)
	TFCSLateralShapeParametrizationHitChain (TFCSLateralShapeParametrizationHitBase *hitsim)
virtual FCSReturnCode	init_hit (TFCSLateralShapeParametrizationHitBase::Hit &hit, TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const
virtual FCSReturnCode	simulate (TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const override
	Method in all derived classes to do some simulation.
virtual unsigned int	size () const override
	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 ind) const override
	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 ind) override
	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 ind, TFCSParametrizationBase *param) override
	Some derived classes have daughter instances of TFCSParametrizationBase objects The set_daughter method allows to change these daughters - expert use only!
const Chain_t &	chain () const
Chain_t &	chain ()
void	push_back (const Chain_t::value_type &value)
void	push_back_init (const Chain_t::value_type &value)
unsigned int	get_nr_of_init () const
void	set_nr_of_init (unsigned int ninit)
virtual void	set_number_of_hits_simul (TFCSLateralShapeParametrizationHitBase *sim)
	set which instance should determine the number of hits
virtual int	get_number_of_hits (TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const
	Call get_number_of_hits() only once, as it could contain a random number.
virtual float	get_E_hit (TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const
	Get hit energy from layer energy and number of hits.
virtual float	get_sigma2_fluctuation (TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const
	Give the effective size sigma^2 of the fluctuations that should be generated.
virtual float	getMinWeight () const
	Get minimum and maximum value of weight for hit energy reweighting.
virtual float	getMaxWeight () const
void	Print (Option_t *option="") const override
	Do not persistify!
bool	is_match_Ekin_bin (int bin) const override
bool	is_match_calosample (int calosample) const override
virtual bool	is_match_all_Ekin_bin () const override
virtual bool	is_match_all_calosample () const override
int	Ekin_bin () const
void	set_Ekin_bin (int bin)
int	calosample () const
void	set_calosample (int cs)
virtual void	set_pdgid_Ekin_eta_Ekin_bin_calosample (const TFCSLateralShapeParametrization &ref)
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
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 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 std::string	startMsg (MSG::Level lvl, const std::string &file, int line)
	Make a message to decorate the start of logging.

Static Public Attributes
static constexpr float	s_max_sigma2_fluctuation = 1000

Protected Member Functions
void	PropagateMSGLevel (MSG::Level level) const
virtual bool	check_all_hits_simulated (TFCSLateralShapeParametrizationHitBase::Hit &hit, TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol, bool success) const
bool	compare (const TFCSParametrizationBase &ref) const

Protected Attributes
Chain_t	m_chain

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
TFCSLateralShapeParametrizationHitBase *	m_number_of_hits_simul
unsigned int	m_ninit = 0
int	m_Ekin_bin
int	m_calosample
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 14 of file TFCSLateralShapeParametrizationHitChain.h.

Member Typedef Documentation

Chain_t

typedef std::vector<TFCSLateralShapeParametrizationHitBase *> TFCSLateralShapeParametrizationHitChain::Chain_t

Definition at line 43 of file TFCSLateralShapeParametrizationHitChain.h.

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

TFCSLateralShapeParametrizationHitChain() [1/2]

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

Definition at line 47 of file TFCSLateralShapeParametrizationHitChain.cxx.

    : TFCSLateralShapeParametrization(name, title),
      m_number_of_hits_simul(nullptr) {}

TFCSLateralShapeParametrizationHitChain() [2/2]

TFCSLateralShapeParametrizationHitChain::TFCSLateralShapeParametrizationHitChain

(

TFCSLateralShapeParametrizationHitBase *

hitsim

)

Definition at line 53 of file TFCSLateralShapeParametrizationHitChain.cxx.

    : TFCSLateralShapeParametrization(TString("hit_chain_") + hitsim->GetName(),
                                      TString("hit chain for ") +
                                          hitsim->GetTitle()),
      m_number_of_hits_simul(nullptr) {
  set_pdgid_Ekin_eta_Ekin_bin_calosample(*hitsim);
 
  m_chain.push_back(hitsim);
}

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

calosample()

int TFCSLateralShapeParametrization::calosample ( ) const

inlineinherited

Definition at line 34 of file TFCSLateralShapeParametrization.h.

{ return m_calosample; };

chain() [1/2]

Chain_t & TFCSLateralShapeParametrizationHitChain::chain ( )

inline

Definition at line 51 of file TFCSLateralShapeParametrizationHitChain.h.

{ return m_chain; };

chain() [2/2]

const Chain_t & TFCSLateralShapeParametrizationHitChain::chain ( ) const

inline

Definition at line 50 of file TFCSLateralShapeParametrizationHitChain.h.

{ return m_chain; };

check_all_hits_simulated()

bool TFCSLateralShapeParametrizationHitChain::check_all_hits_simulated ( TFCSLateralShapeParametrizationHitBase::Hit & hit, TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol, bool success ) const

protectedvirtual

Reimplemented in TFCSLateralShapeParametrizationFixedHitChain.

Definition at line 457 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                                {
 
  (void)truth;     // unused parameter
  (void)extrapol;  // unused parameter
 
  if (!success) {
    // The hit simulation failed => Nothing changed => Still not done
    return false;
  }
 
  float sumEhit = 0;
  if (simulstate.hasAuxInfo("FCSHitChainEnergySum"_FCShash) && hit.idx() > 1) {
    sumEhit = simulstate.getAuxInfo<float>("FCSHitChainEnergySum"_FCShash);
  }
 
  sumEhit += hit.E();
 
  simulstate.setAuxInfo<float>("FCSHitChainEnergySum"_FCShash, sumEhit);
 
  float Elayer = simulstate.E(calosample());
  return (std::abs(sumEhit) >= std::abs(Elayer));
}

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 TFCSLateralShapeParametrization::compare ( const TFCSParametrizationBase & ref ) const

protectedinherited

Definition at line 32 of file TFCSLateralShapeParametrization.cxx.

                                              {
  if (IsA() != ref.IsA()) {
    ATH_MSG_DEBUG("compare(): different class types "
                  << IsA()->GetName() << " != " << ref.IsA()->GetName());
    return false;
  }
  const TFCSLateralShapeParametrization &ref_typed =
      static_cast<const TFCSLateralShapeParametrization &>(ref);
  if (Ekin_bin() != ref_typed.Ekin_bin()) {
    ATH_MSG_DEBUG("compare(): different Ekin bin");
    return false;
  }
  if (calosample() != ref_typed.calosample()) {
    ATH_MSG_DEBUG("compare(): different calosample");
    return false;
  }
 
  return true;
}

Ekin_bin()

int TFCSLateralShapeParametrization::Ekin_bin ( ) const

inlineinherited

Definition at line 31 of file TFCSLateralShapeParametrization.h.

{ return m_Ekin_bin; };

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

get_E_hit()

float TFCSLateralShapeParametrizationHitChain::get_E_hit ( TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol ) const

virtual

Get hit energy from layer energy and number of hits.

Reimplemented in TFCSLateralShapeParametrizationFluctChain.

Definition at line 162 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                  {
  const int nhits = get_number_of_hits(simulstate, truth, extrapol);
 
  // Store nhits in the simulation state as the get function may only be
  // called once per event
  simulstate.setAuxInfo<int>("FCSHitChainNHits"_FCShash, nhits);
 
  const int sample = calosample();
  if (sample < 0)
    return 0;
  if (nhits <= 0)
    return simulstate.E(sample);
  const float maxWeight =
      getMaxWeight(); // maxWeight = -1 if  shapeWeight class is not in m_chain
 
  if (maxWeight > 0)
    return simulstate.E(sample) /
           (maxWeight *
            nhits); // maxWeight is used only if shapeWeight class is in m_chain
  else
    return simulstate.E(sample) /
           nhits; // Otherwise, old definition of E_hit is used
}

get_nr_of_init()

unsigned int TFCSLateralShapeParametrizationHitChain::get_nr_of_init ( ) const

inline

Definition at line 57 of file TFCSLateralShapeParametrizationHitChain.h.

{ return m_ninit; };

get_number_of_hits()

int TFCSLateralShapeParametrizationHitChain::get_number_of_hits ( TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol ) const

virtual

Call get_number_of_hits() only once, as it could contain a random number.

Definition at line 125 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                  {
  // TODO: should we still do it?
  if (m_number_of_hits_simul) {
    int n =
        m_number_of_hits_simul->get_number_of_hits(simulstate, truth, extrapol);
    if (n < 1)
      n = 1;
    return n;
  }
  for (TFCSLateralShapeParametrizationHitBase *hitsim : m_chain) {
    int n = hitsim->get_number_of_hits(simulstate, truth, extrapol);
    if (n > 0)
      return n;
  }
  return 1;
}

get_sigma2_fluctuation()

float TFCSLateralShapeParametrizationHitChain::get_sigma2_fluctuation ( TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol ) const

virtual

Give the effective size sigma^2 of the fluctuations that should be generated.

Definition at line 188 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                  {
  if (m_number_of_hits_simul) {
    double sigma2 = m_number_of_hits_simul->get_sigma2_fluctuation(
        simulstate, truth, extrapol);
    if (sigma2 > 0)
      return sigma2;
  }
  for (TFCSLateralShapeParametrizationHitBase *hitsim : m_chain) {
    double sigma2 = hitsim->get_sigma2_fluctuation(simulstate, truth, extrapol);
    if (sigma2 > 0)
      return sigma2;
  }
  // Limit to factor s_max_sigma2_fluctuation fluctuations
  return s_max_sigma2_fluctuation;
}

getMaxWeight()

float TFCSLateralShapeParametrizationHitChain::getMaxWeight ( ) const

virtual

Definition at line 153 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                                  {
  for (TFCSLateralShapeParametrizationHitBase *hitsim : m_chain) {
    float weight = hitsim->getMaxWeight();
    if (weight > 0.)
      return weight;
  }
  return -1.;
}

getMinWeight()

float TFCSLateralShapeParametrizationHitChain::getMinWeight ( ) const

virtual

Get minimum and maximum value of weight for hit energy reweighting.

Definition at line 144 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                                  {
  for (TFCSLateralShapeParametrizationHitBase *hitsim : m_chain) {
    float weight = hitsim->getMinWeight();
    if (weight > 0.)
      return weight;
  }
  return -1.;
}

init_hit()

FCSReturnCode TFCSLateralShapeParametrizationHitChain::init_hit ( TFCSLateralShapeParametrizationHitBase::Hit & hit, TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol ) const

virtual

Definition at line 212 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                  {
  hit.reset_center();
  if (get_nr_of_init() > 0) {
    ATH_MSG_DEBUG("E(" << calosample() << ")=" << simulstate.E(calosample())
                       << " before init");
 
    auto initloopend = m_chain.begin() + get_nr_of_init();
    for (auto hititr = m_chain.begin(); hititr != initloopend; ++hititr) {
      TFCSLateralShapeParametrizationHitBase *hitsim = *hititr;
 
      FCSReturnCode status =
          hitsim->simulate_hit(hit, simulstate, truth, extrapol);
 
      if (status != FCSSuccess) {
        ATH_MSG_ERROR("TFCSLateralShapeParametrizationHitChain::simulate(): "
                      "simulate_hit init call failed");
        return FCSFatal;
      }
    }
  }
  return FCSSuccess;
}

is_match_all_calosample()

virtual bool TFCSLateralShapeParametrization::is_match_all_calosample ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrization.

Definition at line 29 of file TFCSLateralShapeParametrization.h.

{ return false; };

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 TFCSLateralShapeParametrization::is_match_all_Ekin_bin ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrization.

Definition at line 24 of file TFCSLateralShapeParametrization.h.

                                                      {
    if (Ekin_bin() == -1)
      return true;
    return false;
  };

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

bool TFCSLateralShapeParametrization::is_match_calosample ( int calosample ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 20 of file TFCSLateralShapeParametrization.h.

                                                          {
    return calosample == m_calosample;
  };

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 TFCSLateralShapeParametrization::is_match_Ekin_bin ( int bin ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 15 of file TFCSLateralShapeParametrization.h.

                                                 {
    if (Ekin_bin() == -1)
      return true;
    return bin == Ekin_bin();
  };

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

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]

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

overridevirtual

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

Reimplemented from TFCSParametrizationBase.

Definition at line 84 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                                          {
  if (m_number_of_hits_simul) {
    if (ind == 0)
      return m_number_of_hits_simul;
    return m_chain[ind - 1];
  } else {
    return m_chain[ind];
  }
}

operator[]() [2/2]

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

overridevirtual

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

Reimplemented from TFCSParametrizationBase.

Definition at line 95 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                                    {
  if (m_number_of_hits_simul) {
    if (ind == 0)
      return m_number_of_hits_simul;
    return m_chain[ind - 1];
  } else {
    return m_chain[ind];
  }
}

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

override

Do not persistify!

Definition at line 483 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                                          {
  TFCSLateralShapeParametrization::Print(option);
  TString opt(option);
  bool shortprint = opt.Index("short") >= 0;
  bool longprint = msgLvl(MSG::DEBUG) || (msgLvl(MSG::INFO) && !shortprint);
  TString optprint = opt;
  optprint.ReplaceAll("short", "");
 
  if (m_number_of_hits_simul) {
    if (longprint)
      ATH_MSG_INFO(optprint << "#:Number of hits simulation:");
    m_number_of_hits_simul->Print(opt + "#:");
  }
  if (longprint && get_nr_of_init() > 0)
    ATH_MSG_INFO(optprint << "> Simulation init chain:");
  auto hitloopstart = m_chain.begin() + get_nr_of_init();
  for (auto hititr = m_chain.begin(); hititr != hitloopstart; ++hititr) {
    TFCSLateralShapeParametrizationHitBase *hitsim = *hititr;
    hitsim->Print(opt + "> ");
  }
  if (longprint)
    ATH_MSG_INFO(optprint << "- Simulation chain:");
  char count = 'A';
  for (auto hititr = hitloopstart; hititr != m_chain.end(); ++hititr) {
    TFCSLateralShapeParametrizationHitBase *hitsim = *hititr;
    hitsim->Print(opt + count + " ");
    count++;
  }
}

PropagateMSGLevel()

void TFCSLateralShapeParametrizationHitChain::PropagateMSGLevel ( MSG::Level level ) const

protected

Definition at line 206 of file TFCSLateralShapeParametrizationHitChain.cxx.

                          {
  for (TFCSLateralShapeParametrizationHitBase *reset : m_chain)
    reset->setLevel(level);
}

push_back()

void TFCSLateralShapeParametrizationHitChain::push_back ( const Chain_t::value_type & value )

inline

Definition at line 52 of file TFCSLateralShapeParametrizationHitChain.h.

                                                 {
    m_chain.push_back(value);
  };

push_back_init()

void TFCSLateralShapeParametrizationHitChain::push_back_init

(

const Chain_t::value_type &

value

)

Definition at line 72 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                    {
  if (m_ninit == size()) {
    chain().push_back(value);
  } else {
    const auto it = chain().begin() + m_ninit;
    chain().insert(it, value);
  }
  ++m_ninit;
}

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

set_calosample()

void TFCSLateralShapeParametrization::set_calosample ( int cs )

inherited

Definition at line 21 of file TFCSLateralShapeParametrization.cxx.

                                                           {
  m_calosample = cs;
}

set_daughter()

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

overridevirtual

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 from TFCSParametrizationBase.

Definition at line 105 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                      {
  TFCSLateralShapeParametrizationHitBase *param_typed = nullptr;
  if (param != nullptr) {
    if (!param->InheritsFrom(TFCSLateralShapeParametrizationHitBase::Class())) {
      ATH_MSG_ERROR("Wrong class type " << param->IsA()->GetName());
      return;
    }
    param_typed = static_cast<TFCSLateralShapeParametrizationHitBase *>(param);
  }
  if (m_number_of_hits_simul) {
    if (ind == 0)
      m_number_of_hits_simul = param_typed;
    else
      m_chain.at(ind - 1) = param_typed;
  } else {
    m_chain.at(ind) = param_typed;
  }
}

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

void TFCSLateralShapeParametrization::set_Ekin_bin ( int bin )

inherited

Definition at line 17 of file TFCSLateralShapeParametrization.cxx.

                                                          {
  m_Ekin_bin = bin;
}

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

void TFCSLateralShapeParametrizationHitChain::set_nr_of_init ( unsigned int ninit )

inline

Definition at line 58 of file TFCSLateralShapeParametrizationHitChain.h.

{ m_ninit = ninit; };

set_number_of_hits_simul()

virtual void TFCSLateralShapeParametrizationHitChain::set_number_of_hits_simul ( TFCSLateralShapeParametrizationHitBase * sim )

inlinevirtual

set which instance should determine the number of hits

Definition at line 62 of file TFCSLateralShapeParametrizationHitChain.h.

                                                                        {
    m_number_of_hits_simul = sim;
  };

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

set_pdgid_Ekin_eta_Ekin_bin_calosample()

void TFCSLateralShapeParametrization::set_pdgid_Ekin_eta_Ekin_bin_calosample ( const TFCSLateralShapeParametrization & ref )

virtualinherited

Definition at line 25 of file TFCSLateralShapeParametrization.cxx.

                                                {
  set_calosample(ref.calosample());
  set_Ekin_bin(ref.Ekin_bin());
  set_pdgid_Ekin_eta(ref);
}

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

Reimplemented in TFCSLateralShapeParametrizationFluctChain.

Definition at line 238 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                  {
  MSG::Level old_level = level();
  const bool debug = msgLvl(MSG::DEBUG);
  const bool verbose = msgLvl(MSG::VERBOSE);
 
  // Execute the first get_nr_of_init() simulate calls only once. Used for
  // example to initialize the center position
  TFCSLateralShapeParametrizationHitBase::Hit hit;
  if (init_hit(hit, simulstate, truth, extrapol) != FCSSuccess) {
    ATH_MSG_ERROR("init_hit() failed");
    return FCSFatal;
  }
 
  int cs = calosample();
  // Initialize hit energy only now, as init loop above might change the layer
  // energy
  const float Elayer = simulstate.E(cs);
  if (Elayer == 0) {
    ATH_MSG_VERBOSE("Elayer=0, nothing to do");
    return FCSSuccess;
  }
  const float Ehit = get_E_hit(simulstate, truth, extrapol);
  if (Ehit * Elayer <= 0) {
    ATH_MSG_ERROR("TFCSLateralShapeParametrizationHitChain::simulate(): Ehit "
                  "and Elayer have different sign! Ehit="
                  << Ehit << " Elayer=" << Elayer);
    return FCSFatal;
  }
 
  // Call get_number_of_hits() only once inside get_E_hit(...),
  // as it could contain a random number
  //TODO: should read simulate value!
  int nhit_signed = TMath::Nint(Elayer / Ehit);
  if (nhit_signed < 1)
    nhit_signed = 1;
  unsigned int nhit = (unsigned int)nhit_signed;
 
  float sumEhit = 0;
 
  if (debug) {
    PropagateMSGLevel(old_level);
    ATH_MSG_DEBUG("E(" << cs << ")=" << simulstate.E(cs) << " #hits~" << nhit);
  }
 
  //
  // simulate the hits in GPU
  //
#ifdef USE_GPU
  int ichn = 0;
  bool our_chainA = false;
  if (cs > 0 && cs < 8 && cs != 4)
    our_chainA = true;
  if (nhit > MIN_GPU_HITS && our_chainA && simulstate.get_geold() != nullptr) {
    ATH_MSG_DEBUG(
        "running FastCaloSim in GPU, nhits = " << nhit << " ; Ehit = " << Ehit);
    GeoLoadGpu *gld = (GeoLoadGpu *)simulstate.get_geold();
 
    Chain0_Args args;
    args.debug = debug;
    args.cs = cs;
    args.extrapol_eta_ent = extrapol->eta(cs, SUBPOS_ENT);
    args.extrapol_eta_ext = extrapol->eta(cs, SUBPOS_EXT);
    args.extrapol_phi_ent = extrapol->phi(cs, SUBPOS_ENT);
    args.extrapol_phi_ext = extrapol->phi(cs, SUBPOS_EXT);
    args.extrapol_r_ent = extrapol->r(cs, SUBPOS_ENT);
    args.extrapol_r_ext = extrapol->r(cs, SUBPOS_EXT);
    args.extrapol_z_ent = extrapol->z(cs, SUBPOS_ENT);
    args.extrapol_z_ext = extrapol->z(cs, SUBPOS_EXT);
    args.pdgId = truth->pdgid();
    args.charge = HepPDT::ParticleID(args.pdgId).charge();
    args.nhits = nhit;
    args.rand = 0;
    args.geo = gld->get_geoPtr();
    args.rd4h = simulstate.get_gpu_rand();
    args.ncells = gld->get_ncells();
    ichn = 0;
    // reweight the cell energies for pion
    bool reweight = false;
    for (auto hitsim : m_chain) {
      if (ichn == 0) {
        args.extrapWeight =
            ((TFCSCenterPositionCalculation *)hitsim)->getExtrapWeight();
      }
 
      if (ichn == 1) {
 
        args.is_phi_symmetric = ((TFCSHistoLateralShapeParametrization *)hitsim)
                                    ->is_phi_symmetric();
        args.fh2d =
            ((TFCSHistoLateralShapeParametrization *)hitsim)->LdFH()->hf2d_d();
        args.fh2d_h = *(
            ((TFCSHistoLateralShapeParametrization *)hitsim)->LdFH()->hf2d_h());
      }
      if (ichn == 2 && m_chain.size() == 4) {
        args.fh1d =
            ((TFCSHistoLateralShapeGausLogWeight *)hitsim)->LdFH()->hf1d_d();
        args.fh1d_h =
            *(((TFCSHistoLateralShapeGausLogWeight *)hitsim)->LdFH()->hf1d_h());
        reweight = true;
      }
      if ((ichn == 2 && m_chain.size() == 3) ||
          (ichn == 3 && m_chain.size() == 4)) {
        //      ATH_MSG_DEBUG("---NumberOfBins :" << ( (TFCSHitCellMappingWi
        //      ggle*)hitsim )->get_number_of_bins() );
        args.fhs = ((TFCSHitCellMappingWiggle *)hitsim)->LdFH()->hf_d();
        args.fhs_h = *(((TFCSHitCellMappingWiggle *)hitsim)->LdFH()->hf_h());
      }
 
      ichn++;
    }
    if (reweight)
      CaloGpuGeneral::simulate_hits(Ehit, 1.2 * nhit, args, reweight);
    else
      CaloGpuGeneral::simulate_hits(Ehit, nhit, args, reweight);
 
    for (unsigned int ii = 0; ii < args.ct; ++ii) {
 
      const CaloDetDescrElement_Gpu *cellele =
          gld->index2cell(args.hitcells_E_h[ii].cellid);
      const CaloDetDescrElement *cell_h =
          m_geo->getDDE(cs, cellele->eta(), cellele->phi());
      sumEhit += args.hitcells_E_h[ii].energy;
      if (reweight) {
        if (std::abs(sumEhit) > std::abs(Elayer))
          simulstate.deposit(cell_h, args.hitcells_E_h[ii].energy);
      }
    }
  } else
#endif
  {
 
    auto hitloopstart = m_chain.begin() + get_nr_of_init();
    hit.idx() = 0;
    int ifail = 0;
    int itotalfail = 0;
    int retry_warning = 1;
    int retry = 0;
    bool failed = false;
    do {
      hit.reset();
      hit.E() = Ehit;
      failed = false;
      if (debug)
        if (hit.idx() == 2)
          if (!verbose) {
            // Switch debug output back to INFO to avoid huge logs, but keep
            // full log in verbose
            PropagateMSGLevel(MSG::INFO);
          }
      for (auto hititr = hitloopstart; hititr != m_chain.end(); ++hititr) {
        TFCSLateralShapeParametrizationHitBase *hitsim = *hititr;
 
        FCSReturnCode status =
            hitsim->simulate_hit(hit, simulstate, truth, extrapol);
 
        if (status == FCSSuccess)
          continue;
        if (status == FCSFatal) {
          if (debug)
            PropagateMSGLevel(old_level);
          return FCSFatal;
        }
        failed = true;
        ++ifail;
        ++itotalfail;
        retry = status - FCSRetry;
        retry_warning = retry >> 1;
        if (retry_warning < 1)
          retry_warning = 1;
        break;
      }
      if (!failed) {
        ifail = 0;
        sumEhit += hit.E();
        ++hit.idx();
 
        // TODO: Different signedness:
        if (((hit.idx() == 20 * nhit) || (hit.idx() == 100 * nhit)) && hit.idx() >= 100) {
          ATH_MSG_DEBUG(
              "TFCSLateralShapeParametrizationHitChain::simulate(): Iterated "
              << hit.idx() << " times, expected " << nhit << " times. Deposited E("
              << cs << ")=" << sumEhit << " expected E=" << Elayer);
        }
        if (hit.idx() >= 1000 * nhit && hit.idx() >= 1000) {
          ATH_MSG_DEBUG("TFCSLateralShapeParametrizationHitChain::simulate():"
                          " Aborting hit chain, iterated "
                          << hit.idx() << " times, expected " << nhit
                          << " times. Deposited E(" << cs << ")=" << sumEhit
                          << " expected E=" << Elayer << ", caused by:");
          if (debug) Print();
          break;
        }
      } else {
        if (ifail >= retry) {
          ATH_MSG_ERROR("TFCSLateralShapeParametrizationHitChain::simulate(): "
                        "simulate_hit call failed after "
                        << ifail << "/" << retry
                        << "retries, total fails=" << itotalfail);
          if (debug)
            PropagateMSGLevel(old_level);
          return FCSFatal;
        }
        if (ifail >= retry_warning) {
          ATH_MSG_WARNING("TFCSLateralShapeParametrizationHitChain::simulate():"
                          " retry simulate_hit calls "
                          << ifail << "/" << retry
                          << ", total fails=" << itotalfail);
        }
      }
    } while (!check_all_hits_simulated(hit, simulstate, truth, extrapol, !failed));
  }
 
  if (debug)
    PropagateMSGLevel(old_level);
  return FCSSuccess;
}

size()

unsigned int TFCSLateralShapeParametrizationHitChain::size ( ) const

overridevirtual

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

Reimplemented from TFCSParametrizationBase.

Definition at line 65 of file TFCSLateralShapeParametrizationHitChain.cxx.

                                                                 {
  if (m_number_of_hits_simul)
    return m_chain.size() + 1;
  else
    return m_chain.size();
}

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

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_calosample

int TFCSLateralShapeParametrization::m_calosample

privateinherited

Definition at line 47 of file TFCSLateralShapeParametrization.h.

m_chain

Chain_t TFCSLateralShapeParametrizationHitChain::m_chain

protected

Definition at line 99 of file TFCSLateralShapeParametrizationHitChain.h.

m_Ekin_bin

int TFCSLateralShapeParametrization::m_Ekin_bin

privateinherited

Definition at line 46 of file TFCSLateralShapeParametrization.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_ninit

unsigned int TFCSLateralShapeParametrizationHitChain::m_ninit = 0

private

Definition at line 109 of file TFCSLateralShapeParametrizationHitChain.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_hits_simul

TFCSLateralShapeParametrizationHitBase* TFCSLateralShapeParametrizationHitChain::m_number_of_hits_simul

private

Definition at line 108 of file TFCSLateralShapeParametrizationHitChain.h.

m_pdgid

std::set<int> TFCSParametrization::m_pdgid

privateinherited

Definition at line 63 of file TFCSParametrization.h.

s_max_sigma2_fluctuation

float TFCSLateralShapeParametrizationHitChain::s_max_sigma2_fluctuation = 1000

staticconstexpr

Definition at line 87 of file TFCSLateralShapeParametrizationHitChain.h.

---

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

• TFCSLateralShapeParametrizationHitChain.h
• TFCSLateralShapeParametrizationHitChain.cxx