TFCSBinnedShower Class Reference

#include <TFCSBinnedShower.h>

Inheritance diagram for TFCSBinnedShower:

Collaboration diagram for TFCSBinnedShower:

Classes
struct	event_t
struct	layer_bins_t
struct	layer_t

Public Types
typedef std::vector< event_t >	eventvector_t
typedef std::vector< layer_bins_t >	event_bins_t
enum	FCSEnergyInitializationStatusBits { kOnlyScaleEnergy = }
	Status bit for energy initialization. More...
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
	TFCSBinnedShower (const char *name=nullptr, const char *title=nullptr)
virtual	~TFCSBinnedShower ()
virtual bool	is_match_Ekin_bin (int) const override
void	enable_event_matching ()
void	disable_event_matching ()
void	enable_event_cherry_picking ()
void	disable_event_cherry_picking ()
void	enable_eta_matching ()
void	disable_eta_matching ()
void	enable_upscaling ()
void	disable_upscaling ()
void	set_default_hit_energy (float energy)
float	get_default_hit_energy () const
void	set_max_hits_per_voxel (int max_hits)
int	get_max_hits_per_voxel () const
void	load_event_library (const std::string &filename, std::vector< long unsigned int > &layers, bool only_load_meta_data=false)
void	load_sub_bin_distribution (const std::string &filename)
void	set_hdf5_path (const std::string &filename)
void	delete_hdf5_path ()
const std::string &	get_hdf5_path () const
void	set_layer_energy (long unsigned int event_index, long unsigned int layer_index, const std::vector< unsigned int > &bin_index_vector, const std::vector< float > &E_vector)
void	set_bin_boundaries (long unsigned int layer_index, std::vector< float > &R_lower, std::vector< float > &R_size, std::vector< float > &alpha_lower, std::vector< float > &alpha_size)
void	set_shower_center_information (long unsigned int event_index, long unsigned int reference_layer_index, float eta_center, float phi_center)
const eventvector_t &	get_eventlibrary ()
void	set_event_library (eventvector_t &eventlibrary)
const event_bins_t &	get_coordinates ()
void	set_coordinates (event_bins_t &coordinates)
const std::vector< std::vector< std::vector< std::vector< float > > > > &	get_sub_bin_distribution () const
const std::vector< float > &	get_upscaling_energies () const
void	set_sub_bin_distribution_and_energies (std::vector< std::vector< std::vector< std::vector< float > > > > &sub_bin_distribution, std::vector< float > &upscaling_energies)
bool	OnlyScaleEnergy () const
void	set_OnlyScaleEnergy ()
void	reset_OnlyScaleEnergy ()
virtual FCSReturnCode	simulate (TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const override
	Method in all derived classes to do some simulation.
virtual FCSReturnCode	simulate_hit (Hit &hit, TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) override
	simulated one hit position with some energy.
ICaloGeometry *	get_geometry ()
virtual void	set_geometry (ICaloGeometry *geo) override
	Method to set the geometry access pointer.
virtual int	get_number_of_hits (TFCSSimulationState &simulstate, const TFCSTruthState *, const TFCSExtrapolationState *) const override
	Call get_number_of_hits() only once per shower simulation, as it could be calculated with random numbers and give different results each time.
virtual double	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 by the amount of generated hits.
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	getMinWeight () const
	Get minimum and maximum value of weight for hit energy reweighting.
virtual float	getMaxWeight () const
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	Print (Option_t *option="") const override
void	clear ()
virtual bool	is_match_pdgid (int id) const override
virtual bool	is_match_Ekin (float Ekin) const override
virtual bool	is_match_eta (float eta) const override
virtual bool	is_match_all_Ekin () const override
virtual bool	is_match_all_eta () const override
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 unsigned int	size () const
	Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.
virtual const TFCSParametrizationBase *	operator[] (unsigned int) const
	Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.
virtual TFCSParametrizationBase *	operator[] (unsigned int)
	Some derived classes have daughter instances of TFCSParametrizationBase objects The size() and operator[] methods give general access to these daughters.
virtual void	set_daughter (unsigned int, TFCSParametrizationBase *)
	Some derived classes have daughter instances of TFCSParametrizationBase objects The set_daughter method allows to change these daughters - expert use only!
virtual bool	operator== (const TFCSParametrizationBase &ref) const
	The == operator compares the content of instances.
virtual void	CleanAuxInfo (TFCSSimulationState &) const
	Method in all derived classes to delete objects stored in the simulstate AuxInfo.
void	FindDuplicates (FindDuplicateClasses_t &dup)
void	RemoveDuplicates ()
void	RemoveNameTitle ()
bool	msgLvl (const MSG::Level lvl) const
	Check whether the logging system is active at the provided verbosity level.
MsgStream &	msg () const
	Return a stream for sending messages directly (no decoration)
MsgStream &	msg (const MSG::Level lvl) const
	Return a decorated starting stream for sending messages.
MSG::Level	level () const
	Retrieve output level.
virtual void	setLevel (MSG::Level lvl)
	Update outputlevel.

Static Public Member Functions
static std::string	startMsg (MSG::Level lvl, const std::string &file, int line)
	Make a message to decorate the start of logging.

Protected Member Functions
virtual void	get_event (TFCSSimulationState &simulstate, float eta_center, float phi_center, float e_init, long unsigned int reference_layer_index) const override
	do not persistify
virtual long unsigned int	get_n_hits (TFCSSimulationState &simulstate, long unsigned int layer_index) const override
virtual void	compute_n_hits_and_elayer (TFCSSimulationState &simulstate) const
virtual float	get_layer_energy (TFCSSimulationState &simulstate, long unsigned int layer_index) const override
virtual std::tuple< float, float, float >	get_hit_position_and_energy (TFCSSimulationState &simulstate, long unsigned int layer_index, long unsigned int hit_index) const override
virtual void	delete_event (TFCSSimulationState &simulstate) const override
bool	compare (const TFCSParametrizationBase &ref) const

Protected Attributes
ICaloGeometry *	m_geo {}

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 Member Functions
long unsigned int	find_best_match (float eta_center, float phi_center, float e_init, long unsigned int reference_layer_index, bool phi_mod_matching) const
std::tuple< float, float >	get_coordinates (TFCSSimulationState &simulstate, long unsigned int layer_index, int bin_index) const
void	upscale (TFCSSimulationState &simulstate, float &R_min, float &R_max, float &alpha_min, float &alpha_max, long unsigned int layer_index, int bin_index) const
long unsigned int	get_energy_index (TFCSSimulationState &simulstate, long unsigned int layer_index, long unsigned int hit_index) const
std::tuple< std::vector< float >, std::vector< hsize_t >, bool >	load_hdf5_dataset (const std::string &filename, const std::string &datasetname)
void	load_layer_energy (const std::string &filename, long unsigned int layer_index)
void	load_bin_boundaries (const std::string &filename, long unsigned int layer_index)
void	load_shower_center_information (const std::string &filename)

Private Attributes
const long unsigned int	m_n_layers
float	m_default_hit_energy = 4.
int	m_max_hits_per_voxel = 100
std::string	m_hdf5_file
eventvector_t	m_eventlibrary
event_bins_t	m_coordinates
bool	m_use_event_matching = false
bool	m_use_event_cherry_picking = false
bool	m_use_eta_matching = false
bool	m_use_upscaling = false
std::vector< std::vector< std::vector< std::vector< float > > > >	m_sub_bin_distribution
std::vector< float >	m_upscaling_energies
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 23 of file TFCSBinnedShower.h.

Member Typedef Documentation

event_bins_t

typedef std::vector<layer_bins_t> TFCSBinnedShower::event_bins_t

Definition at line 46 of file TFCSBinnedShower.h.

eventvector_t

typedef std::vector<event_t> TFCSBinnedShower::eventvector_t

Definition at line 37 of file TFCSBinnedShower.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

FCSEnergyInitializationStatusBits

enum TFCSBinnedShowerBase::FCSEnergyInitializationStatusBits

inherited

Status bit for energy initialization.

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

Definition at line 24 of file TFCSBinnedShowerBase.h.

                                         {
    kOnlyScaleEnergy =
        BIT(18)  
  };

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

TFCSBinnedShower()

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

Definition at line 50 of file TFCSBinnedShower.cxx.

    : TFCSBinnedShowerBase(name, title) {}

~TFCSBinnedShower()

TFCSBinnedShower::~TFCSBinnedShower ( )

virtual

Definition at line 53 of file TFCSBinnedShower.cxx.

{}

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

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

compute_n_hits_and_elayer()

void TFCSBinnedShower::compute_n_hits_and_elayer ( TFCSSimulationState & simulstate ) const

protectedvirtual

Definition at line 221 of file TFCSBinnedShower.cxx.

                                           {
 
  event_t *event = static_cast<event_t *>(
      simulstate.getAuxInfo<void *>("BSEventData"_FCShash));
  float e_init = simulstate.getAuxInfo<float>("BSEinit"_FCShash);
 
  // Loop over all layers
  long unsigned int n_layers = event->event_data.size();
  std::vector<std::vector<long unsigned int>> hits_per_layer;
  std::vector<float> elayer;
  hits_per_layer.resize(n_layers);
  elayer.resize(n_layers, 0.0f);
 
  for (long unsigned int layer_index = 0; layer_index < n_layers;
       ++layer_index) {
 
    layer_t &layer = event->event_data.at(layer_index);
 
    // Loop over all voxels in the layer
    long unsigned int n_hits = 0;
    for (float e_voxel : layer.E_vector) {
      long unsigned int hits_per_bin;
      if (e_voxel > std::numeric_limits<float>::epsilon()) {
        hits_per_bin =
            std::min(std::max(int(e_voxel * e_init / m_default_hit_energy), 1),
                     m_max_hits_per_voxel);
        elayer.at(layer_index) += e_voxel * e_init;
 
      } else {
        hits_per_bin = 0;
      }
      n_hits += hits_per_bin;
      hits_per_layer.at(layer_index).push_back(n_hits);
    }
  }
  // Store the hits per layer vector
  std::vector<std::vector<long unsigned int>> *hits_per_layer_ptr =
      new std::vector<std::vector<long unsigned int>>(std::move(hits_per_layer));
  simulstate.setAuxInfo<void *>("BSNHits"_FCShash, hits_per_layer_ptr);
 
  // Store the energy per layer
  std::vector<float> *elayer_ptr = new std::vector<float>(std::move(elayer));
  simulstate.setAuxInfo<void *>("BSELayer"_FCShash, elayer_ptr);
}

delete_event()

void TFCSBinnedShower::delete_event ( TFCSSimulationState & simulstate ) const

overrideprotectedvirtual

Implements TFCSBinnedShowerBase.

Definition at line 509 of file TFCSBinnedShower.cxx.

                                                                         {
  // Delete the event data
  void *event_ptr = simulstate.getAuxInfo<void *>("BSEventData"_FCShash);
  if (event_ptr) {
    delete static_cast<event_t *>(event_ptr);
    simulstate.setAuxInfo<void *>("BSEventData"_FCShash, nullptr);
  } else {
    ATH_MSG_ERROR("No event data found to delete.");
  }
 
  void *n_hits_ptr = simulstate.getAuxInfo<void *>("BSNHits"_FCShash);
  if (n_hits_ptr) {
    delete static_cast<std::vector<std::vector<long unsigned int>> *>(
        n_hits_ptr);
    simulstate.setAuxInfo<void *>("BSNHits"_FCShash, nullptr);
  } else {
    ATH_MSG_ERROR("No event hits data found to delete.");
  }
 
  void *elayer_ptr = simulstate.getAuxInfo<void *>("BSELayer"_FCShash);
  if (elayer_ptr) {
    delete static_cast<std::vector<float> *>(elayer_ptr);
    simulstate.setAuxInfo<void *>("BSELayer"_FCShash, nullptr);
  } else {
    ATH_MSG_ERROR("No event layer energy data found to delete.");
  }
 
  return;
}

delete_hdf5_path()

void TFCSBinnedShower::delete_hdf5_path ( )

inline

Definition at line 102 of file TFCSBinnedShower.h.

{ m_hdf5_file.clear(); }

disable_eta_matching()

void TFCSBinnedShower::disable_eta_matching ( )

inline

Definition at line 78 of file TFCSBinnedShower.h.

{ m_use_eta_matching = false; };

disable_event_cherry_picking()

void TFCSBinnedShower::disable_event_cherry_picking ( )

inline

Definition at line 70 of file TFCSBinnedShower.h.

{ m_use_event_cherry_picking = false; };

disable_event_matching()

void TFCSBinnedShower::disable_event_matching ( )

inline

Definition at line 62 of file TFCSBinnedShower.h.

{ m_use_event_matching = false; }

disable_upscaling()

void TFCSBinnedShower::disable_upscaling ( )

inline

Definition at line 81 of file TFCSBinnedShower.h.

{ m_use_upscaling = false; };

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

enable_eta_matching()

void TFCSBinnedShower::enable_eta_matching ( )

inline

Definition at line 72 of file TFCSBinnedShower.h.

                             {
    m_use_eta_matching = true;
    m_use_event_matching = false;
    m_use_event_cherry_picking = false;
  };

enable_event_cherry_picking()

void TFCSBinnedShower::enable_event_cherry_picking ( )

inline

Definition at line 64 of file TFCSBinnedShower.h.

                                     {
    m_use_event_cherry_picking = true;
    m_use_event_matching = false;
    m_use_eta_matching = false;
  };

enable_event_matching()

void TFCSBinnedShower::enable_event_matching ( )

inline

Definition at line 56 of file TFCSBinnedShower.h.

                               {
    m_use_event_matching = true;
    m_use_event_cherry_picking = false;
    m_use_eta_matching = false;
  };

enable_upscaling()

void TFCSBinnedShower::enable_upscaling ( )

inline

Definition at line 80 of file TFCSBinnedShower.h.

{ m_use_upscaling = true; };

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

find_best_match()

long unsigned int TFCSBinnedShower::find_best_match ( float eta_center, float phi_center, float e_init, long unsigned int reference_layer_index, bool phi_mod_matching ) const

private

Definition at line 115 of file TFCSBinnedShower.cxx.

                                                                          {
 
  float phi_cell_size, phi_within_cell;
 
  if (phi_mod_matching) {
    const CaloDetDescrElement *cellele =
        m_geo->getDDE(reference_layer_index, eta_center, phi_center);
    phi_cell_size = TFCSPhiModulationCorrection::get_phi_cell_size(
        reference_layer_index, eta_center);
 
    float phi_cell = cellele->phi();
    phi_within_cell = phi_center - phi_cell;
 
    phi_within_cell = fmod(phi_within_cell, phi_cell_size);
    if (phi_within_cell < 0)
      phi_within_cell += phi_cell_size;
  }
 
  // Find the event with the closest eta_center and phi_mod (L2 distance)
  float best_distance = std::numeric_limits<float>::max();
  long unsigned int best_match = m_eventlibrary.size() + 1;
 
  for (long unsigned int event_index = 0; event_index < m_eventlibrary.size();
       ++event_index) {
 
    // Absolute eta difference. Normalized by the (usual) eta range of 0.05...
    float eta_diff =
        (m_eventlibrary.at(event_index).center_eta - eta_center) / 0.05;
 
    // Logarithmic energy difference. Use power of two (distance between usual
    // energy points) as reasonable distance measure.
    float e_diff =
        std::log(m_eventlibrary.at(event_index).e_init / e_init) / std::log(2);
 
    // float dist2 = eta_diff * eta_diff;
    float dist2 = eta_diff * eta_diff + e_diff * e_diff;
 
    if (phi_mod_matching) {
      float phi_mod = m_eventlibrary.at(event_index).phi_mod;
      float phi_diff = (phi_mod - phi_within_cell) / phi_cell_size;
      dist2 = dist2 + phi_diff * phi_diff;
    }
    float distance = TMath::Sqrt(dist2);
 
    if (distance < best_distance) {
      best_distance = distance;
      best_match = event_index;
    }
  }
 
  if (best_match == m_eventlibrary.size() + 1) {
    ATH_MSG_ERROR("No best match found");
  } else {
    ATH_MSG_INFO("Best match found for eta "
                 << eta_center << " / "
                 << m_eventlibrary.at(best_match).center_eta << " with energy "
                 << e_init << " / " << m_eventlibrary.at(best_match).e_init);
  }
 
  return best_match;
}

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_coordinates() [1/2]

const event_bins_t & TFCSBinnedShower::get_coordinates ( )

inline

Definition at line 130 of file TFCSBinnedShower.h.

{ return m_coordinates; }

get_coordinates() [2/2]

std::tuple< float, float > TFCSBinnedShower::get_coordinates ( TFCSSimulationState & simulstate, long unsigned int layer_index, int bin_index ) const

private

Definition at line 333 of file TFCSBinnedShower.cxx.

                         {
  float R_min = m_coordinates.at(layer_index).R_lower.at(bin_index);
  float R_max = m_coordinates.at(layer_index).R_size.at(bin_index) +
                m_coordinates.at(layer_index).R_lower.at(bin_index);
 
  float alpha_min = m_coordinates.at(layer_index).alpha_lower.at(bin_index);
  float alpha_max = m_coordinates.at(layer_index).alpha_size.at(bin_index) +
                    m_coordinates.at(layer_index).alpha_lower.at(bin_index);
 
  if (m_use_upscaling) {
    upscale(simulstate, R_min, R_max, alpha_min, alpha_max, layer_index,
            bin_index);
  }
  float R;
  if (TMath::Abs(R_max - R_min) > std::numeric_limits<float>::epsilon()) {
    R = CLHEP::RandFlat::shoot(simulstate.randomEngine(), R_min, R_max);
  } else {
    R = R_min;  // If the range is too small, just use the minimum value
  }
  float alpha =
      CLHEP::RandFlat::shoot(simulstate.randomEngine(), alpha_min, alpha_max);
 
  return std::make_tuple(R, alpha);
}

get_default_hit_energy()

float TFCSBinnedShower::get_default_hit_energy ( ) const

inline

Definition at line 84 of file TFCSBinnedShower.h.

{ return m_default_hit_energy; }

get_E_hit()

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

virtualinherited

Get hit energy from layer energy and number of hits.

Definition at line 31 of file TFCSLateralShapeParametrizationHitBase.cxx.

                                                  {
  const int nhits = get_number_of_hits(simulstate, truth, extrapol);
  const int sample = calosample();
  if (nhits <= 0 || sample < 0)
    return -1.;
  else
    return simulstate.E(sample) / nhits;
}

get_energy_index()

long unsigned int TFCSBinnedShower::get_energy_index ( TFCSSimulationState & simulstate, long unsigned int layer_index, long unsigned int hit_index ) const

private

Definition at line 296 of file TFCSBinnedShower.cxx.

                                       {
  std::vector<std::vector<long unsigned int>> *hits_per_layer_ptr =
      static_cast<std::vector<std::vector<long unsigned int>> *>(
          simulstate.getAuxInfo<void *>("BSNHits"_FCShash));
  if (!hits_per_layer_ptr) {
    ATH_MSG_ERROR("Invalid hits per layer information");
    return 0;
  }
 
  if (layer_index >= hits_per_layer_ptr->size()) {
    ATH_MSG_ERROR("Layer index out of bounds: " << layer_index << " >= "
                                                << hits_per_layer_ptr->size());
    return 0;
  }
 
  // Find the hit index in the hit vector for the given layer
  const std::vector<long unsigned int> &hits =
      hits_per_layer_ptr->at(layer_index);
  auto it = std::upper_bound(hits.begin(), hits.end(), hit_index);
  long unsigned int energy_index = std::distance(hits.begin(), it);
 
  if (energy_index >= hits.size()) {
    ATH_MSG_ERROR("Energy index out of bounds: " << energy_index
                                                 << " >= " << hits.size());
    // Print full hits for debugging
    ATH_MSG_ERROR("Hits per layer: ");
    for (const auto &hit : hits) {
      ATH_MSG_ERROR("Hit: " << hit);
    }
    return 0;
  }
 
  return energy_index;
}

get_event()

void TFCSBinnedShower::get_event ( TFCSSimulationState & simulstate, float eta_center, float phi_center, float e_init, long unsigned int reference_layer_index ) const

overrideprotectedvirtual

do not persistify

Implements TFCSBinnedShowerBase.

Definition at line 179 of file TFCSBinnedShower.cxx.

                                                                 {
 
  if (m_eventlibrary.empty()) {
    ATH_MSG_ERROR(
        "No event library loaded. Please load an event libray for "
        "TFCSBinnedShower::get_event.");
    return;
  }
 
  simulstate.setAuxInfo<float>("BSEinit"_FCShash, e_init);
 
  long unsigned int event_index;
  if (simulstate.hasAuxInfo("EventNr"_FCShash) && m_use_event_matching) {
    event_index = simulstate.getAuxInfo<int>("EventNr"_FCShash);
  } else if (m_use_event_cherry_picking || m_use_eta_matching) {
    // dphi/deta
    event_index =
        find_best_match(eta_center, phi_center, e_init, reference_layer_index,
                        m_use_event_cherry_picking);
 
    if (event_index >= m_eventlibrary.size()) {
      event_index = std::floor(CLHEP::RandFlat::shoot(
          simulstate.randomEngine(), 0, m_eventlibrary.size()));
    }
 
  } else {
    event_index = std::floor(CLHEP::RandFlat::shoot(simulstate.randomEngine(),
                                                    0, m_eventlibrary.size()));
  }
 
  ATH_MSG_DEBUG("Using event index " << event_index << " for eta " << eta_center
                                     << " and phi " << phi_center);
 
  // Store a pointer to the event
  event_t *event_ptr = new event_t(m_eventlibrary.at(event_index));
  simulstate.setAuxInfo<void *>("BSEventData"_FCShash, event_ptr);
 
  compute_n_hits_and_elayer(simulstate);
}

get_eventlibrary()

const eventvector_t & TFCSBinnedShower::get_eventlibrary ( )

inline

Definition at line 124 of file TFCSBinnedShower.h.

{ return m_eventlibrary; }

get_geometry()

ICaloGeometry * TFCSBinnedShowerBase::get_geometry ( )

inlineinherited

Definition at line 46 of file TFCSBinnedShowerBase.h.

{ return m_geo; };

get_hdf5_path()

const std::string & TFCSBinnedShower::get_hdf5_path ( ) const

inline

Definition at line 103 of file TFCSBinnedShower.h.

{ return m_hdf5_file; }

get_hit_position_and_energy()

std::tuple< float, float, float > TFCSBinnedShower::get_hit_position_and_energy ( TFCSSimulationState & simulstate, long unsigned int layer_index, long unsigned int hit_index ) const

overrideprotectedvirtual

Implements TFCSBinnedShowerBase.

Definition at line 467 of file TFCSBinnedShower.cxx.

                                       {
 
  event_t *event = static_cast<event_t *>(
      simulstate.getAuxInfo<void *>("BSEventData"_FCShash));
 
  float e_init = simulstate.getAuxInfo<float>("BSEinit"_FCShash);
 
  if (layer_index >= event->event_data.size()) {
    ATH_MSG_ERROR("Layer index out of bounds: " << layer_index << " >= "
                                                << event->event_data.size());
    return std::make_tuple(0.0f, 0.0f, 0.0f);
  }
 
  long unsigned int energy_index = get_energy_index(
      simulstate, layer_index, hit_index);  // Get the bin index for the hit
 
  std::vector<std::vector<long unsigned int>> *hits_per_layer_ptr =
      static_cast<std::vector<std::vector<long unsigned int>> *>(
          simulstate.getAuxInfo<void *>("BSNHits"_FCShash));
 
  long unsigned int hits_per_bin;
  if (energy_index == 0) {
    hits_per_bin = hits_per_layer_ptr->at(layer_index).at(energy_index);
  } else {
    hits_per_bin = hits_per_layer_ptr->at(layer_index).at(energy_index) -
                   hits_per_layer_ptr->at(layer_index).at(energy_index - 1);
  }
 
  float r, alpha;
 
  layer_t &layer = event->event_data.at(layer_index);
 
  std::tie(r, alpha) = get_coordinates(simulstate, layer_index,
                                       layer.bin_index_vector.at(energy_index));
 
  float E = layer.E_vector.at(energy_index) * e_init / hits_per_bin;
 
  return std::make_tuple(r, alpha, E);
}

get_layer_energy()

float TFCSBinnedShower::get_layer_energy ( TFCSSimulationState & simulstate, long unsigned int layer_index ) const

overrideprotectedvirtual

Implements TFCSBinnedShowerBase.

Definition at line 282 of file TFCSBinnedShower.cxx.

                                                                              {
  std::vector<float> *elayer_ptr = static_cast<std::vector<float> *>(
      simulstate.getAuxInfo<void *>("BSELayer"_FCShash));
  if (!elayer_ptr) {
    ATH_MSG_ERROR("Invalid layer energy information");
    return 0.0f;
  }
  if (layer_index >= elayer_ptr->size()) {
    return 0.0f;
  }
  return elayer_ptr->at(layer_index);
}

get_max_hits_per_voxel()

int TFCSBinnedShower::get_max_hits_per_voxel ( ) const

inline

Definition at line 86 of file TFCSBinnedShower.h.

{ return m_max_hits_per_voxel; }

get_n_hits()

long unsigned int TFCSBinnedShower::get_n_hits ( TFCSSimulationState & simulstate, long unsigned int layer_index ) const

overrideprotectedvirtual

Implements TFCSBinnedShowerBase.

Definition at line 267 of file TFCSBinnedShower.cxx.

                                                                          {
 
  std::vector<std::vector<long unsigned int>> *hits_per_layer_ptr =
      static_cast<std::vector<std::vector<long unsigned int>> *>(
          simulstate.getAuxInfo<void *>("BSNHits"_FCShash));
 
  if (!hits_per_layer_ptr) {
    ATH_MSG_ERROR("Invalid hits per layer information");
    return 0;
  }
 
  return hits_per_layer_ptr->at(layer_index).back();
}

get_number_of_hits()

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

inlineoverridevirtualinherited

Call get_number_of_hits() only once per shower simulation, as it could be calculated with random numbers and give different results each time.

Return a value of -1 if this instance can't determine

Reimplemented from TFCSLateralShapeParametrizationHitBase.

Definition at line 53 of file TFCSBinnedShowerBase.h.

                                                       {
    long unsigned int nhits = get_n_hits(simulstate, calosample());
    if (nhits == 0)
      return -1;
    return static_cast<int>(nhits);
  }

get_sigma2_fluctuation()

double TFCSLateralShapeParametrizationHitBase::get_sigma2_fluctuation ( TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol ) const

virtualinherited

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

Return a value <0 if this instance can't determine

Reimplemented in TFCS2DFunctionLateralShapeParametrization, TFCSHistoLateralShapeParametrization, and TFCSLateralShapeParametrizationHitNumberFromE.

Definition at line 19 of file TFCSLateralShapeParametrizationHitBase.cxx.

                                            {
  return -1;
}

get_sub_bin_distribution()

const std::vector< std::vector< std::vector< std::vector< float > > > > & TFCSBinnedShower::get_sub_bin_distribution ( ) const

inline

Definition at line 137 of file TFCSBinnedShower.h.

                                   {
    return m_sub_bin_distribution;
  }

get_upscaling_energies()

const std::vector< float > & TFCSBinnedShower::get_upscaling_energies ( ) const

inline

Definition at line 141 of file TFCSBinnedShower.h.

                                                          {
    return m_upscaling_energies;
  }

getMaxWeight()

float TFCSLateralShapeParametrizationHitBase::getMaxWeight ( ) const

virtualinherited

Reimplemented in TFCSHistoLateralShapeWeight.

Definition at line 46 of file TFCSLateralShapeParametrizationHitBase.cxx.

                                                                 {
  return -1.;
}

getMinWeight()

float TFCSLateralShapeParametrizationHitBase::getMinWeight ( ) const

virtualinherited

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

Reimplemented in TFCSHistoLateralShapeWeight.

Definition at line 42 of file TFCSLateralShapeParametrizationHitBase.cxx.

                                                                 {
  return -1.;
}

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

virtual bool TFCSBinnedShower::is_match_Ekin_bin ( int ) const

inlineoverridevirtual

Reimplemented from TFCSParametrizationBase.

Definition at line 52 of file TFCSBinnedShower.h.

                                                       {
    return true;
  };

is_match_eta()

virtual bool TFCSParametrization::is_match_eta ( float eta ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 21 of file TFCSParametrization.h.

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

is_match_pdgid()

virtual bool TFCSParametrization::is_match_pdgid ( int id ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 15 of file TFCSParametrization.h.

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

level()

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

inlineinherited

Retrieve output level.

Definition at line 201 of file MLogging.h.

{ return msg().level(); }

load_bin_boundaries()

void TFCSBinnedShower::load_bin_boundaries ( const std::string & filename, long unsigned int layer_index )

private

Definition at line 652 of file TFCSBinnedShower.cxx.

                                                                          {
 
  // Assert that the layer index is valid
  if (layer_index >= m_coordinates.size()) {
    m_coordinates.resize(layer_index + 1);
  }
 
  std::vector<std::string> datasetnames = {
      "binstart_radius_layer_", "binsize_radius_layer_",
      "binstart_alpha_layer_", "binsize_alpha_layer_"};
 
  for (long unsigned int i = 0; i < datasetnames.size(); i++) {
    std::string datasetname = datasetnames.at(i) + std::to_string(layer_index);
    std::vector<float> data;
    std::vector<hsize_t> dims;
    bool success;
    std::tie(data, dims, success) = load_hdf5_dataset(filename, datasetname);
    if (!success) {
      ATH_MSG_ERROR("Error while extracting the bin boundaries for layer "
                    << layer_index << " from " << filename << "."
                    << "Specifically, the key " << datasetname
                    << " could not be loaded.");
      return;
    }
 
    // Fill the corresponding vector in the layer_bins_t structure
    auto &event_bins = m_coordinates.at(layer_index);
    switch (i) {
      case 0:
        event_bins.R_lower = std::move(data);
        break;
      case 1:
        event_bins.R_size = std::move(data);
        break;
      case 2:
        event_bins.alpha_lower = std::move(data);
        break;
      case 3:
        event_bins.alpha_size = std::move(data);
        break;
    }
  }
}

load_event_library()

void TFCSBinnedShower::load_event_library	(	const std::string &	filename,
		std::vector< long unsigned int > &	layers,
		bool	only_load_meta_data = false )

Definition at line 539 of file TFCSBinnedShower.cxx.

                              {
 
  if (!only_load_meta_data) {
    m_eventlibrary.clear();
  }
 
  m_coordinates.clear();
 
  if (m_use_upscaling) {
    m_sub_bin_distribution.clear();
  }
 
  // layer dependent variables
  for (long unsigned int layer_index : layers) {
 
    ATH_MSG_INFO("Loading layer " << layer_index << " from file: " << filename);
 
    // Load the bin boundaries for this layer
    load_bin_boundaries(filename, layer_index);
 
    if (!only_load_meta_data) {
      // Load the layer energies
      load_layer_energy(filename, layer_index);
    }
 
  }
 
  if ((m_use_event_cherry_picking || m_use_eta_matching) &&
      !only_load_meta_data) {
    load_shower_center_information(filename);
  }
 
  return;
}

load_hdf5_dataset()

std::tuple< std::vector< float >, std::vector< hsize_t >, bool > TFCSBinnedShower::load_hdf5_dataset ( const std::string & filename, const std::string & datasetname )

private

Definition at line 577 of file TFCSBinnedShower.cxx.

                                                                  {
 
  // Open the HDF5 file and dataset
  H5::H5File file(filename, H5F_ACC_RDONLY);
 
  // check if the dataset exists
  if (!file.exists(datasetname)) {
    return std::make_tuple(std::vector<float>{}, std::vector<hsize_t>{}, false);
  }
 
  H5::DataSet dataset = file.openDataSet(datasetname);
 
  // Get the dataspace of the dataset
  H5::DataSpace dataspace = dataset.getSpace();
 
  // Get the number of dimensions and the size of each dimension
  int rank = dataspace.getSimpleExtentNdims();
  std::vector<hsize_t> dims_out(rank);
  dataspace.getSimpleExtentDims(dims_out.data(), NULL);
 
  // Calculate the total number of elements
  hsize_t totalSize = 1;
  for (const auto &dim : dims_out) {
    totalSize *= dim;
  }
 
  // Read the dataset into a buffer
  std::vector<float> data(totalSize);
  dataset.read(data.data(), H5::PredType::NATIVE_FLOAT);
  file.close();
  return std::make_tuple(data, dims_out, true);
}

load_layer_energy()

void TFCSBinnedShower::load_layer_energy ( const std::string & filename, long unsigned int layer_index )

private

Definition at line 611 of file TFCSBinnedShower.cxx.

                                                                        {
 
  std::string datasetname = "energy_layer_" + std::to_string(layer_index);
 
  std::vector<float> data;
  std::vector<hsize_t> dims;
  bool success;
  std::tie(data, dims, success) = load_hdf5_dataset(filename, datasetname);
  if (!success) {
    ATH_MSG_ERROR("Error while extracting the layer energy for layer "
                  << layer_index << " from " << filename << ".");
    return;
  }
 
  // Store the data in the event library
  std::vector<unsigned int> bin_index_vector;
  std::vector<float> E_vector;
 
  for (size_t i = 0; i < data.size(); ++i) {
    long unsigned int event_index = i / dims.at(1);
    float bin_index = i % dims.at(1);
 
    if (bin_index == 0) {
      bin_index_vector.clear();
      E_vector.clear();
    }
 
    if (data.at(i) != 0.0) {
      bin_index_vector.push_back(bin_index);
      E_vector.push_back(data.at(i));
    }
 
    if (bin_index ==
        dims.at(1) - 1)  // True for the last voxel of the event in this layer
    {
      set_layer_energy(event_index, layer_index, bin_index_vector, E_vector);
    }
  }
}

load_shower_center_information()

void TFCSBinnedShower::load_shower_center_information ( const std::string & filename )

private

Definition at line 697 of file TFCSBinnedShower.cxx.

                               {
 
  // Open the HDF5 file
  H5::H5File file(filename, H5F_ACC_RDONLY);
 
  // Open the dataset
  std::vector<std::string> datasetnames = {"phi_mod", "center_eta",
                                           "incident_energy"};
 
  for (long unsigned int i = 0; i < datasetnames.size(); i++) {
    std::string datasetname = datasetnames.at(i);
    std::vector<float> data;
    std::vector<hsize_t> dims;
    bool success;
    std::tie(data, dims, success) = load_hdf5_dataset(filename, datasetname);
    if (!success) {
      if (datasetname == "phi_mod" && m_use_eta_matching) {
        // We do not necessarly need the phi_mod for eta matching, so we can
        // just skip this dataset
        continue;
      } else {
        ATH_MSG_ERROR(
            "Error while extracting the shower center information from "
            << filename << "."
            << "Specifically, the key " << datasetname
            << " could not be loaded.");
      }
      return;
    }
 
    for (long unsigned int event_index = 0; event_index < data.size();
         ++event_index) {
      switch (i) {
        case 0:
          m_eventlibrary.at(event_index).phi_mod = data.at(event_index);
          break;
        case 1:
          m_eventlibrary.at(event_index).center_eta = data.at(event_index);
          break;
        case 2:
          m_eventlibrary.at(event_index).e_init = data.at(event_index);
          break;
      }
    }
  }
}

load_sub_bin_distribution()

void TFCSBinnedShower::load_sub_bin_distribution

(

const std::string &

filename

)

Definition at line 776 of file TFCSBinnedShower.cxx.

                                                                          {
  m_use_upscaling = true;
  TFile *file = TFile::Open(filename.c_str(), "READ");
  if (!file || file->IsZombie()) {
    std::cerr << "Failed to open file: " << filename << std::endl;
    return;
  }
 
  std::regex pattern(R"(probabilities_layer_(\d+)_energy_([0-9.]+))");
  std::map<float, std::map<int, std::vector<std::vector<float>>>> temp_storage;
 
  TIter next(file->GetListOfKeys());
  TKey *key;
 
  while ((key = (TKey *)next())) {
    std::string keyname = key->GetName();
    std::smatch match;
    if (std::regex_match(keyname, match, pattern)) {
      int layer = std::stoi(match[1].str());
      float energy = std::stod(match[2].str());
 
      TMatrixD *matrix = dynamic_cast<TMatrixD *>(file->Get(keyname.c_str()));
      if (matrix) {
        std::vector<std::vector<float>> mat_vec(
            matrix->GetNrows(), std::vector<float>(matrix->GetNcols()));
        for (int i = 0; i < matrix->GetNrows(); ++i) {
          for (int j = 0; j < matrix->GetNcols(); ++j) {
            mat_vec[i][j] = static_cast<float>((*matrix)(i, j));
          }
        }
        temp_storage[energy][layer] = std::move(mat_vec);
      }
    }
  }
 
  file->Close();
  delete file;
 
  // Output containers
  std::vector<float> energies;
  std::vector<std::vector<std::vector<std::vector<float>>>>
      data;  // [energy][layer][row][col]
 
  for (const auto &[energy, layer_map] : temp_storage) {
    energies.push_back(energy);
    int max_layer = 0;
    for (const auto &[l, _] : layer_map)
      max_layer = std::max(max_layer, l);
 
    std::vector<std::vector<std::vector<float>>> layer_vec(max_layer + 1);
    for (const auto &[layer_idx, mat] : layer_map) {
      layer_vec[layer_idx] = mat;
    }
    data.push_back(std::move(layer_vec));
  }
 
  // Example output
  for (size_t i = 0; i < energies.size(); ++i) {
    std::cout << "Energy index " << i << ": " << energies[i] << " GeV\n";
    for (size_t j = 0; j < data[i].size(); ++j) {
      if (!data[i][j].empty()) {
        std::cout << "  Layer " << j << " Shape: (" << data[i][j].size() << ", "
                  << data[i][j][0].size() << ")\n";
      }
    }
  }
 
  m_upscaling_energies = std::move(energies);
  m_sub_bin_distribution = std::move(data);
}

msg() [1/2]

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

inlineinherited

Return a stream for sending messages directly (no decoration)

Definition at line 231 of file MLogging.h.

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

msg() [2/2]

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

inlineinherited

Return a decorated starting stream for sending messages.

Definition at line 240 of file MLogging.h.

                                                        {
  return msg() << lvl;
}

msgLvl()

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

inlineinherited

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

Definition at line 222 of file MLogging.h.

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

OnlyScaleEnergy()

bool TFCSBinnedShowerBase::OnlyScaleEnergy ( ) const

inlineinherited

Definition at line 30 of file TFCSBinnedShowerBase.h.

{ return TestBit(kOnlyScaleEnergy); };

operator==()

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

inlinevirtualinherited

The == operator compares the content of instances.

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

Reimplemented in TFCSHitCellMapping, TFCSHitCellMappingWiggle, TFCSLateralShapeParametrizationHitNumberFromE, and TFCSPredictExtrapWeights.

Definition at line 119 of file TFCSParametrizationBase.h.

                                                                    {
    return compare(ref);
  };

operator[]() [1/2]

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

inlinevirtualinherited

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

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 106 of file TFCSParametrizationBase.h.

                                                              {
    return nullptr;
  };

operator[]() [2/2]

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

inlinevirtualinherited

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

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 99 of file TFCSParametrizationBase.h.

                                   {
    return nullptr;
  };

pdgid()

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

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 34 of file TFCSParametrization.h.

{ return m_pdgid; };

Print()

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

overrideinherited

Definition at line 53 of file TFCSLateralShapeParametrization.cxx.

                                                                  {
  TString opt(option);
  bool shortprint = opt.Index("short") >= 0;
  bool longprint = msgLvl(MSG::DEBUG) || (msgLvl(MSG::INFO) && !shortprint);
  TString optprint = opt;
  optprint.ReplaceAll("short", "");
  TFCSParametrization::Print(option);
  if (longprint) {
    if (Ekin_bin() == -1)
      ATH_MSG_INFO(optprint << "  Ekin_bin=all ; calosample=" << calosample());
    else
      ATH_MSG_INFO(optprint << "  Ekin_bin=" << Ekin_bin()
                            << " ; calosample=" << calosample());
  }
}

RemoveDuplicates()

void TFCSParametrizationBase::RemoveDuplicates ( )

inherited

Definition at line 152 of file TFCSParametrizationBase.cxx.

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

RemoveNameTitle()

void TFCSParametrizationBase::RemoveNameTitle ( )

inherited

Definition at line 210 of file TFCSParametrizationBase.cxx.

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

reset_match_all_pdgid()

virtual void TFCSParametrizationBase::reset_match_all_pdgid ( )

inlinevirtualinherited

Definition at line 84 of file TFCSParametrizationBase.h.

{ ResetBit(kMatchAllPDGID); };

reset_OnlyScaleEnergy()

void TFCSBinnedShowerBase::reset_OnlyScaleEnergy ( )

inlineinherited

Definition at line 34 of file TFCSBinnedShowerBase.h.

{ ResetBit(kOnlyScaleEnergy); };

set_bin_boundaries()

void TFCSBinnedShower::set_bin_boundaries	(	long unsigned int	layer_index,
		std::vector< float > &	R_lower,
		std::vector< float > &	R_size,
		std::vector< float > &	alpha_lower,
		std::vector< float > &	alpha_size )

Definition at line 76 of file TFCSBinnedShower.cxx.

                                                                        {
  if (layer_index >= m_coordinates.size()) {
    m_coordinates.resize(layer_index + 1);
  }
  m_coordinates.at(layer_index).R_lower = R_lower;
  m_coordinates.at(layer_index).R_size = R_size;
  m_coordinates.at(layer_index).alpha_lower = alpha_lower;
  m_coordinates.at(layer_index).alpha_size = alpha_size;
}

set_calosample()

void TFCSLateralShapeParametrization::set_calosample ( int cs )

inherited

Definition at line 21 of file TFCSLateralShapeParametrization.cxx.

                                                           {
  m_calosample = cs;
}

set_coordinates()

void TFCSBinnedShower::set_coordinates ( event_bins_t & coordinates )

inline

Definition at line 132 of file TFCSBinnedShower.h.

                                                  {
    m_coordinates = coordinates;
  }

set_daughter()

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

inlinevirtualinherited

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

The original element at this position is not deleted

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 113 of file TFCSParametrizationBase.h.

                                                        {};

set_default_hit_energy()

void TFCSBinnedShower::set_default_hit_energy ( float energy )

inline

Definition at line 83 of file TFCSBinnedShower.h.

{ m_default_hit_energy = energy; }

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

void TFCSBinnedShower::set_event_library ( eventvector_t & eventlibrary )

inline

Definition at line 126 of file TFCSBinnedShower.h.

                                                      {
    m_eventlibrary = eventlibrary;
  }

set_geometry()

virtual void TFCSBinnedShowerBase::set_geometry ( ICaloGeometry * geo )

inlineoverridevirtualinherited

Method to set the geometry access pointer.

Loops over daughter objects if present

Reimplemented from TFCSParametrizationBase.

Definition at line 48 of file TFCSBinnedShowerBase.h.

                                                         {
    m_geo = geo;
    TFCSParametrizationBase::set_geometry(geo);
  };

set_hdf5_path()

void TFCSBinnedShower::set_hdf5_path ( const std::string & filename )

inline

Definition at line 101 of file TFCSBinnedShower.h.

{ m_hdf5_file = filename; }

set_layer_energy()

void TFCSBinnedShower::set_layer_energy	(	long unsigned int	event_index,
		long unsigned int	layer_index,
		const std::vector< unsigned int > &	bin_index_vector,
		const std::vector< float > &	E_vector )

Definition at line 55 of file TFCSBinnedShower.cxx.

                                      {
 
  // Assert that the event index is valid
  if (event_index >= m_eventlibrary.size()) {
    m_eventlibrary.resize(event_index + 1);
  }
 
  // Assert that the layer index is valid
  if (layer_index >= m_eventlibrary.at(event_index).event_data.size()) {
    m_eventlibrary.at(event_index).event_data.resize(layer_index + 1);
  }
 
  // Set the layer energy
  layer_t &layer = m_eventlibrary.at(event_index).event_data.at(layer_index);
  layer.bin_index_vector = bin_index_vector;
  layer.E_vector = E_vector;
}

set_match_all_pdgid()

virtual void TFCSParametrizationBase::set_match_all_pdgid ( )

inlinevirtualinherited

Definition at line 83 of file TFCSParametrizationBase.h.

{ SetBit(kMatchAllPDGID); };

set_max_hits_per_voxel()

void TFCSBinnedShower::set_max_hits_per_voxel ( int max_hits )

inline

Definition at line 85 of file TFCSBinnedShower.h.

{ m_max_hits_per_voxel = max_hits; }

set_OnlyScaleEnergy()

void TFCSBinnedShowerBase::set_OnlyScaleEnergy ( )

inlineinherited

Definition at line 32 of file TFCSBinnedShowerBase.h.

{ SetBit(kOnlyScaleEnergy); };

set_pdgid() [1/2]

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

virtualinherited

Reimplemented in TFCSEnergyBinParametrization.

Definition at line 33 of file TFCSParametrization.cxx.

{ m_pdgid = ids; }

set_pdgid() [2/2]

void TFCSParametrization::set_pdgid ( int id )

virtualinherited

Reimplemented in TFCSEnergyBinParametrization.

Definition at line 28 of file TFCSParametrization.cxx.

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

set_pdgid_Ekin_eta()

void TFCSParametrization::set_pdgid_Ekin_eta ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 72 of file TFCSParametrization.cxx.

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

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

set_shower_center_information()

void TFCSBinnedShower::set_shower_center_information	(	long unsigned int	event_index,
		long unsigned int	reference_layer_index,
		float	eta_center,
		float	phi_center )

Definition at line 90 of file TFCSBinnedShower.cxx.

                                        {
 
  // Compute phi_mod
  float phi_cell_size = TFCSPhiModulationCorrection::get_phi_cell_size(
      reference_layer_index, eta_center);
  float phi_cell =
      m_geo->getDDE(reference_layer_index, eta_center, phi_center)->phi();
  float phi_within_cell = fmod(phi_center - phi_cell, phi_cell_size);
  if (phi_within_cell < 0) {
    phi_within_cell += phi_cell_size;
  }
 
  if (m_eventlibrary.size() <= event_index) {
    m_eventlibrary.resize(event_index + 1);
  }
  if (m_eventlibrary.at(event_index).event_data.size() <=
      reference_layer_index) {
    m_eventlibrary.at(event_index).event_data.resize(reference_layer_index + 1);
  }
  m_eventlibrary.at(event_index).center_eta = eta_center;
  m_eventlibrary.at(event_index).phi_mod = phi_within_cell;
}

set_sub_bin_distribution_and_energies()

void TFCSBinnedShower::set_sub_bin_distribution_and_energies ( std::vector< std::vector< std::vector< std::vector< float > > > > & sub_bin_distribution, std::vector< float > & upscaling_energies )

inline

Definition at line 145 of file TFCSBinnedShower.h.

                                            {
    m_sub_bin_distribution = sub_bin_distribution;
    m_upscaling_energies = upscaling_energies;
    m_use_upscaling = true;
  }

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

overridevirtualinherited

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 29 of file TFCSBinnedShowerBase.cxx.

                                                  {
 
  // select a random event from the library
  float eta_center, phi_center;
  long unsigned int reference_layer_index = CaloCell_ID_FCS::CaloSample_FCS::EMB2;
  eta_center =
      extrapol->eta(reference_layer_index, TFCSExtrapolationState::SUBPOS_MID);
  if (eta_center > 1.4) {  // Endcap becomes more relevant
    reference_layer_index = CaloCell_ID_FCS::CaloSample_FCS::EME2;
  }
  // TODO: What about the endcap?
  eta_center =
      extrapol->eta(reference_layer_index, TFCSExtrapolationState::SUBPOS_MID);
  phi_center =
      extrapol->phi(reference_layer_index, TFCSExtrapolationState::SUBPOS_MID);
 
  // Fill the total energy and layer energies into simulstate
  float Einit;
  const float Ekin = truth->Ekin();
 
  if (OnlyScaleEnergy())
    Einit = simulstate.E();
  else
    Einit = Ekin;
 
  // Reset the total energy
  simulstate.set_E(0);
 
  get_event(simulstate, eta_center, phi_center, Einit, reference_layer_index);
 
  for (long unsigned int layer_index = 0;
       layer_index < CaloCell_ID_FCS::MaxSample; ++layer_index) {
 
    float layer_energy = get_layer_energy(simulstate, layer_index);
 
    // Reset and set the layer energy
    simulstate.set_E(layer_index, 0);
    simulstate.add_E(layer_index, layer_energy);
  }
 
  if (simulstate.E() > std::numeric_limits<double>::epsilon()) {
    for (int ilayer = 0; ilayer < CaloCell_ID_FCS::MaxSample; ++ilayer) {
      simulstate.set_Efrac(ilayer, simulstate.E(ilayer) / simulstate.E());
    }
  }
  return FCSSuccess;
}

simulate_hit()

FCSReturnCode TFCSBinnedShowerBase::simulate_hit ( Hit & hit, TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol )

overridevirtualinherited

simulated one hit position with some energy.

As last step in TFCSLateralShapeParametrizationHitChain::simulate, the hit should be mapped into a cell and this cell recorded in simulstate. All hits/cells should be resacled such that their final sum is simulstate->E(sample)

Reimplemented from TFCSLateralShapeParametrizationHitBase.

Definition at line 79 of file TFCSBinnedShowerBase.cxx.

                                     {
 
  const int pdgId = truth->pdgid();
  const float charge = HepPDT::ParticleID(pdgId).charge();
  long unsigned int layer_index = calosample();
 
  const double center_eta = hit.center_eta();
  const double center_phi = hit.center_phi();
  const double center_r = hit.center_r();
  const double center_z = hit.center_z();
 
  ATH_MSG_VERBOSE(" Layer " << layer_index << " Extrap eta " << center_eta
                            << " phi " << center_phi << " R " << center_r);
                            
  //next MR: change to std::abs, std::sqrt functions
  const float dist000 = TMath::Sqrt(center_r * center_r + center_z * center_z);
  const float eta_jakobi = TMath::Abs(2.0 * TMath::Exp(-center_eta) /
                                      (1.0 + TMath::Exp(-2 * center_eta)));
 
  long unsigned int hit_index = hit.idx();
 
  // Get necessary the hit information
  float r, alpha, E;
  std::tie(r, alpha, E) =
      get_hit_position_and_energy(simulstate, layer_index, hit_index);
 
  hit.reset();
  hit.E() = E;
 
  if (layer_index <= CaloCell_ID_FCS::CaloSample_FCS::FCAL0) {
    float delta_eta_mm = r * std::cos(alpha);
    float delta_phi_mm = r * std::sin(alpha);
 
    // Particles with negative eta are expected to have the same shape
    // as those with positive eta after transformation: delta_eta -->
    // -delta_eta
    if (center_eta < 0.) {
      delta_eta_mm = -delta_eta_mm;
    }
 
    // We derive the shower shapes for electrons and positively charged
    // hadrons. Particle with the opposite charge are expected to have the
    // same shower shape after the transformation: delta_phi -->
    // -delta_phi
    if ((charge < 0. && pdgId != 11) || pdgId == -11)
      delta_phi_mm = -delta_phi_mm;
 
    const float delta_eta = delta_eta_mm / eta_jakobi / dist000;
    const float delta_phi = delta_phi_mm / center_r;
 
    hit.eta() = center_eta + delta_eta;
    hit.phi() = TVector2::Phi_mpi_pi(center_phi + delta_phi);
 
    ATH_MSG_VERBOSE(" Hit eta " << hit.eta() << " phi " << hit.phi()
                                << " layer " << layer_index);
 
  } else {  // FCAL is in (x,y,z)
    const float hit_r = r * std::cos(alpha) + center_r;
    float delta_phi = r * std::sin(alpha) / center_r;
    // We derive the shower shapes for electrons and positively charged
    // hadrons. Particle with the opposite charge are expected to have the
    // same shower shape after the transformation: delta_phi -->
    // -delta_phi
    if ((charge < 0. && pdgId != 11) || pdgId == -11)
      delta_phi = -delta_phi;
    const float hit_phi = TVector2::Phi_mpi_pi(center_phi + delta_phi);
    hit.x() = hit_r * std::cos(hit_phi);
    hit.y() = hit_r * std::sin(hit_phi);
    hit.z() = center_z;
    ATH_MSG_VERBOSE(" Hit x " << hit.x() << " y " << hit.y() << " layer "
                              << layer_index);
  }
 
  return FCSSuccess;
}

size()

virtual unsigned int TFCSParametrizationBase::size ( ) const

inlinevirtualinherited

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

Reimplemented in TFCSLateralShapeParametrizationHitChain, and TFCSParametrizationChain.

Definition at line 93 of file TFCSParametrizationBase.h.

{ return 0; };

startMsg()

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

staticinherited

Make a message to decorate the start of logging.

Print a message for the start of logging.

Definition at line 116 of file MLogging.cxx.

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

upscale()

void TFCSBinnedShower::upscale ( TFCSSimulationState & simulstate, float & R_min, float & R_max, float & alpha_min, float & alpha_max, long unsigned int layer_index, int bin_index ) const

private

Definition at line 360 of file TFCSBinnedShower.cxx.

                                                    {
 
  float p = CLHEP::RandFlat::shoot(simulstate.randomEngine(), 0, 1);
 
  float e_init = simulstate.getAuxInfo<float>("BSEinit"_FCShash);
  std::vector<float> available_energies = m_upscaling_energies;
 
  unsigned int e_index = 0;
 
  std::vector<float> probabilities = {0.25f, 0.5f, 0.75f};
 
  if (available_energies.size() > 1) {
    // find closest energy index using binary search
    auto it = std::upper_bound(available_energies.begin(),
                               available_energies.end(), e_init);
    if (it != available_energies.end()) {
      e_index = std::distance(available_energies.begin(), it);
    } else {
      e_index = available_energies.size() - 1;
    }
 
    float e_high = available_energies.at(e_index);
    if (e_high < e_init || e_index == 0) {
      if (m_sub_bin_distribution.at(e_index).size() > layer_index) {
        probabilities =
            m_sub_bin_distribution.at(e_index).at(layer_index).at(bin_index);
      }
    } else {
      if (m_sub_bin_distribution.at(e_index).size() > layer_index &&
          m_sub_bin_distribution.at(e_index - 1).size() > layer_index) {
        float e_low = available_energies.at(e_index - 1);
        float f_low = std::log(e_high / e_init) / (std::log(e_high / e_low));
        float f_high = 1 - f_low;
        for (unsigned int i = 0; i < 3; ++i) {
          float p_low = m_sub_bin_distribution.at(e_index - 1)
                            .at(layer_index)
                            .at(bin_index)
                            .at(i);
          float p_high = m_sub_bin_distribution.at(e_index)
                             .at(layer_index)
                             .at(bin_index)
                             .at(i);
          probabilities[i] = f_low * p_low + f_high * p_high;
        }
      }
    }
  }
 
  else if (available_energies.size() == 1) {
    probabilities = m_sub_bin_distribution.at(0).at(layer_index).at(bin_index);
  }
 
  float p_alpha_low = probabilities[2] - probabilities[1] + probabilities[0];
  float p_r;
 
  if (p < p_alpha_low) {
    alpha_max = (alpha_min + alpha_max) / 2.;
    p_r = probabilities[0] / (p_alpha_low);
  } else {
    alpha_min = (alpha_min + alpha_max) / 2.;
    p_r = (probabilities[1] - probabilities[0]) / (1 - p_alpha_low);
  }
 
  p = CLHEP::RandFlat::shoot(simulstate.randomEngine(), 0, 1);
  if (layer_index != 2){
  // if ((layer_index != 2) && (layer_index != 1)) {
    if (p > p_r) {
      R_min = (R_min + R_max) / 2.;
      return;
    } else {
      R_max = (R_min + R_max) / 2.;
      return;
    }
  }
 
 
  // Use linear interpolation for layer 2
  // It works better than uniform sampling for the second layer...
  if (p_r < 0.25) {
    p_r = 0.25;  // Values below 0.25 are not allowed for linear pdf
  } else if (p_r > 0.75) {
    p_r = 0.75;  // Values above 0.75 are not allowed for linear pdf
  } else if (TMath::Abs(p_r - 0.5) < std::numeric_limits<float>::epsilon()) {
    return;  // Best upscaling is uniform sampling. Nothing to do here.
  }
 
  // Inverse CDF for linear pdf
  float r = (1. - 4. * p_r) / (2. - 4. * p_r) -
            TMath::Sqrt(((1. - 4. * p_r) / (2. - 4. * p_r)) *
                            ((1. - 4. * p_r) / (2. - 4. * p_r)) +
                        p / ((1. / 2.) - p_r));
  if (r < 0) {
    r = (1. - 4. * p_r) / (2. - 4. * p_r) +
        TMath::Sqrt(((1. - 4. * p_r) / (2. - 4. * p_r)) *
                        ((1. - 4. * p_r) / (2. - 4. * p_r)) +
                    p / ((1. / 2.) - p_r));
  }
 
  R_min = R_min + r / 2 * (R_max - R_min);
  R_max = R_min;
 
  return;
}

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_coordinates

event_bins_t TFCSBinnedShower::m_coordinates

private

Definition at line 198 of file TFCSBinnedShower.h.

m_default_hit_energy

float TFCSBinnedShower::m_default_hit_energy = 4.

private

Definition at line 189 of file TFCSBinnedShower.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_eventlibrary

eventvector_t TFCSBinnedShower::m_eventlibrary

private

Definition at line 197 of file TFCSBinnedShower.h.

m_geo

ICaloGeometry* TFCSBinnedShowerBase::m_geo {}

protectedinherited

Definition at line 63 of file TFCSBinnedShowerBase.h.

{};  

m_hdf5_file

std::string TFCSBinnedShower::m_hdf5_file

private

Definition at line 194 of file TFCSBinnedShower.h.

m_max_hits_per_voxel

int TFCSBinnedShower::m_max_hits_per_voxel = 100

private

Definition at line 191 of file TFCSBinnedShower.h.

m_n_layers

const long unsigned int TFCSBinnedShower::m_n_layers

private

Initial value:

=
      CaloCell_ID_FCS::CaloSample_FCS::MaxSample

Definition at line 185 of file TFCSBinnedShower.h.

m_nm

std::string ISF_FCS::MLogging::m_nm

privateinherited

Message source name.

Definition at line 211 of file MLogging.h.

m_pdgid

std::set<int> TFCSParametrization::m_pdgid

privateinherited

Definition at line 63 of file TFCSParametrization.h.

m_sub_bin_distribution

std::vector<std::vector<std::vector<std::vector<float> > > > TFCSBinnedShower::m_sub_bin_distribution

private

Definition at line 233 of file TFCSBinnedShower.h.

m_upscaling_energies

std::vector<float> TFCSBinnedShower::m_upscaling_energies

private

Definition at line 234 of file TFCSBinnedShower.h.

m_use_eta_matching

bool TFCSBinnedShower::m_use_eta_matching = false

private

Definition at line 216 of file TFCSBinnedShower.h.

m_use_event_cherry_picking

bool TFCSBinnedShower::m_use_event_cherry_picking = false

private

Definition at line 211 of file TFCSBinnedShower.h.

m_use_event_matching

bool TFCSBinnedShower::m_use_event_matching = false

private

Definition at line 205 of file TFCSBinnedShower.h.

m_use_upscaling

bool TFCSBinnedShower::m_use_upscaling = false

private

Definition at line 228 of file TFCSBinnedShower.h.

---

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

• TFCSBinnedShower.h
• TFCSBinnedShower.cxx