TFCSEnergyAndHitGANV2 Class Reference

#include <TFCSEnergyAndHitGANV2.h>

Inheritance diagram for TFCSEnergyAndHitGANV2:

Collaboration diagram for TFCSEnergyAndHitGANV2:

Public Types
enum	FCSGANfreemem { kGANfreemem = BIT(17) }
	Status bit for chain persistency. More...
enum	FCSEnergyInitializationStatusBits { kOnlyScaleEnergy = }
	Status bit for energy initialization. More...
enum	FCSSplitChainObjects { kSplitChainObjects , kRetryChainFromStart = BIT(17) }
	Status bit for chain persistency. More...
typedef std::vector< TFCSParametrizationBase * >	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
	TFCSEnergyAndHitGANV2 (const char *name=nullptr, const char *title=nullptr)
virtual	~TFCSEnergyAndHitGANV2 ()
virtual bool	is_match_Ekin_bin (int) const override
virtual 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
bool	GANfreemem () const
void	set_GANfreemem ()
void	reset_GANfreemem ()
bool	OnlyScaleEnergy () const
void	set_OnlyScaleEnergy ()
void	reset_OnlyScaleEnergy ()
virtual int	get_bin (TFCSSimulationState &simulstate, const TFCSTruthState *, const TFCSExtrapolationState *) const override
	use the layer to be done as binning of the GAN chain
virtual const std::string	get_variable_text (TFCSSimulationState &simulstate, const TFCSTruthState *, const TFCSExtrapolationState *) const override
unsigned int	get_nr_of_init (unsigned int bin) const
void	set_nr_of_init (unsigned int bin, unsigned int ninit)
const TFCSGANXMLParameters::Binning &	get_Binning () const
const TFCSGANEtaSlice::ExtrapolatorWeights	get_ExtrapolationWeights () const
bool	initializeNetwork (int const &pid, int const &etaMin, const std::string &FastCaloGANInputFolderName)
bool	fillEnergy (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 void	Print (Option_t *option="") const override
virtual void	push_before_first_bin (TFCSParametrizationBase *param)
virtual void	push_back_in_bin (TFCSParametrizationBase *param, unsigned int bin)
virtual unsigned int	get_number_of_bins () const
virtual const std::string	get_bin_text (int bin) const
	print the range of a bin; for bin -1, print the allowed range
bool	SplitChainObjects () const
void	set_SplitChainObjects ()
void	reset_SplitChainObjects ()
bool	RetryChainFromStart () const
void	set_RetryChainFromStart ()
void	reset_RetryChainFromStart ()
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 &param)
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 void	test_path (const std::string &path, TFCSSimulationState *simulstate=nullptr, const TFCSTruthState *truth=nullptr, const TFCSExtrapolationState *extrapol=nullptr, const std::string &outputname="unnamed", int pid=211)
static void	unit_test (TFCSSimulationState *simulstate=nullptr, const TFCSTruthState *truth=nullptr, const TFCSExtrapolationState *extrapol=nullptr)
static std::string	startMsg (MSG::Level lvl, const std::string &file, int line)
	Make a message to decorate the start of logging.

Protected Member Functions
void	SetRegionAndSliceFromXML (int pid, int etaMax, std::string FastCaloGANInputFolderName)
void	recalc_pdgid_intersect ()
void	recalc_pdgid_union ()
void	recalc_Ekin_intersect ()
void	recalc_eta_intersect ()
void	recalc_Ekin_eta_intersect ()
void	recalc_Ekin_union ()
void	recalc_eta_union ()
void	recalc_Ekin_eta_union ()
virtual void	recalc ()
	Default is to call recalc_pdgid_intersect() and recalc_Ekin_eta_intersect().
FCSReturnCode	simulate_and_retry (TFCSParametrizationBase *parametrization, TFCSSimulationState &simulstate, const TFCSTruthState *truth, const TFCSExtrapolationState *extrapol) const
bool	compare (const TFCSParametrizationBase &ref) const

Protected Attributes
std::vector< unsigned int >	m_bin_start
	Contains the index where the TFCSParametrizationBase* instances to run for a given bin start.

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
int	GetAlphaBinsForRBin (const TAxis *x, int ix, int yBinNum) const

Static Private Member Functions
static int	GetBinsInFours (double const bins)

Private Attributes
std::vector< int >	m_bin_ninit
TFCSGANEtaSlice *	m_slice = nullptr
TFCSGANXMLParameters	m_param
std::mutex	m_mutex
Chain_t	m_chain
std::vector< std::unique_ptr< TFCSParametrizationBase > >	m_writtenBases
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 21 of file TFCSEnergyAndHitGANV2.h.

Member Typedef Documentation

Chain_t

typedef std::vector<TFCSParametrizationBase *> TFCSParametrizationChain::Chain_t

inherited

Definition at line 40 of file TFCSParametrizationChain.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 TFCSEnergyAndHitGANV2::FCSEnergyInitializationStatusBits

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 45 of file TFCSEnergyAndHitGANV2.h.

                                         {
    kOnlyScaleEnergy =
        BIT(18) 
  };

FCSGANfreemem

enum TFCSEnergyAndHitGANV2::FCSGANfreemem

Status bit for chain persistency.

Enumerator
kGANfreemem	Set this bit in the TObject bit if the memory for m_input should be freed after reading in athena.

Definition at line 35 of file TFCSEnergyAndHitGANV2.h.

                     {
    kGANfreemem = BIT(17) 
  };

FCSSplitChainObjects

enum TFCSParametrizationChain::FCSSplitChainObjects

inherited

Status bit for chain persistency.

Enumerator
kSplitChainObjects	Set this bit in the TObject bit field if the TFCSParametrizationBase objects in the chain should be written as separate keys into the root file instead of directly writing the objects. This is needed if the sum of all objects in the chain use >1GB of memory, which can't be handeled by TBuffer. Drawback is that identical objects will get stored as multiple instances
kRetryChainFromStart

Definition at line 20 of file TFCSParametrizationChain.h.

                            {
    kSplitChainObjects =
        BIT(16), 
    kRetryChainFromStart = BIT(17)
  };

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

TFCSEnergyAndHitGANV2()

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

Definition at line 36 of file TFCSEnergyAndHitGANV2.cxx.

    : TFCSParametrizationBinnedChain(name, title) {
  set_GANfreemem();
}

~TFCSEnergyAndHitGANV2()

TFCSEnergyAndHitGANV2::~TFCSEnergyAndHitGANV2 ( )

virtual

Definition at line 42 of file TFCSEnergyAndHitGANV2.cxx.

                                              {
  if (m_slice != nullptr) {
    delete m_slice;
  }
}

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

chain() [1/2]

Chain_t & TFCSParametrizationChain::chain ( )

inlineinherited

Definition at line 54 of file TFCSParametrizationChain.h.

{ return m_chain; };

chain() [2/2]

const Chain_t & TFCSParametrizationChain::chain ( ) const

inlineinherited

Definition at line 53 of file TFCSParametrizationChain.h.

{ return m_chain; };

CleanAuxInfo()

virtual void TFCSParametrizationBase::CleanAuxInfo ( TFCSSimulationState & ) const

inlinevirtualinherited

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

Definition at line 130 of file TFCSParametrizationBase.h.

{};

clear()

void TFCSParametrization::clear ( )

inherited

Definition at line 18 of file TFCSParametrization.cxx.

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

clear_pdgid()

void TFCSParametrization::clear_pdgid ( )

virtualinherited

Reimplemented in TFCSEnergyBinParametrization.

Definition at line 37 of file TFCSParametrization.cxx.

{ m_pdgid.clear(); }

compare()

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

protectedinherited

Definition at line 78 of file TFCSParametrization.cxx.

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

Ekin_max()

double TFCSParametrization::Ekin_max ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 37 of file TFCSParametrization.h.

{ return m_Ekin_max; };

Ekin_min()

double TFCSParametrization::Ekin_min ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 36 of file TFCSParametrization.h.

{ return m_Ekin_min; };

Ekin_nominal()

double TFCSParametrization::Ekin_nominal ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 35 of file TFCSParametrization.h.

{ return m_Ekin_nominal; };

eta_max()

double TFCSParametrization::eta_max ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 40 of file TFCSParametrization.h.

{ return m_eta_max; };

eta_min()

double TFCSParametrization::eta_min ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 39 of file TFCSParametrization.h.

{ return m_eta_min; };

eta_nominal()

double TFCSParametrization::eta_nominal ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 38 of file TFCSParametrization.h.

{ return m_eta_nominal; };

fillEnergy()

bool TFCSEnergyAndHitGANV2::fillEnergy	(	TFCSSimulationState &	simulstate,
		const TFCSTruthState *	truth,
		const TFCSExtrapolationState *	extrapol ) const

Definition at line 114 of file TFCSEnergyAndHitGANV2.cxx.

                                                  {
  if (!truth) {
    ATH_MSG_ERROR("Invalid truth pointer");
    return false;
  }
 
  if (!extrapol) {
    ATH_MSG_ERROR("Invalid extrapolation pointer");
    return false;
  }
 
  const int pdgId = truth->pdgid();
  const float charge = MC::charge(pdgId);
 
  float Einit;
  const float Ekin = truth->Ekin();
  if (OnlyScaleEnergy())
    Einit = simulstate.E();
  else
    Einit = Ekin;
 
  ATH_MSG_VERBOSE("Momentum " << truth->P() << " pdgId " << truth->pdgid());
  // check that the network exists
  if (!m_slice->IsGanCorrectlyLoaded()) {
    ATH_MSG_WARNING("GAN not loaded correctly.");
    return false;
  }
  // This lock is an attempt to fix ATLASSIM-7031. remove if not necessary
  // Hold until NetworkOutputs goes out of scope
  std::scoped_lock lock(m_mutex);
  TFCSGANEtaSlice::NetworkOutputs outputs =
      m_slice->GetNetworkOutputs(truth, extrapol, simulstate);
  ATH_MSG_VERBOSE("network outputs size: " << outputs.size());
 
  const TFCSGANXMLParameters::Binning &binsInLayers = m_param.GetBinning();
  const auto ganVersion = m_param.GetGANVersion();
  const TFCSGANEtaSlice::FitResultsPerLayer &fitResults =
      m_slice->GetFitResults();  // used only if GAN version > 1
 
  ATH_MSG_DEBUG("energy voxels size = " << outputs.size());
 
  double totalEnergy = 0;
  for (const auto & output : outputs) {
    totalEnergy += output.second;
  }
  if (totalEnergy < 0) {
    ATH_MSG_WARNING("Energy from GAN is negative, skipping particle");
    return false;
  }
 
  ATH_MSG_VERBOSE("Get binning");
 
  simulstate.set_E(0);
 
  int vox = 0;
  for (const auto &[layer, h] : binsInLayers) {
    // attempt to debug intermittent ci issues described in
    // https://its.cern.ch/jira/browse/ATLASSIM-7031
    if (*reinterpret_cast<void*const*>(&h) == nullptr || h.IsZombie() || h.IsOnHeap() || dynamic_cast<const TH2D*>(&h) == nullptr) {
      ATH_MSG_ERROR("Histogram for layer " << layer << " at " << &h <<
                    " is broken; " <<
                    "See ATLASSIM-7031.");
 
      std::ostringstream ss;
      ss << "Node dump:\n";
      CxxUtils::safeHexdump (ss, (reinterpret_cast<const char*>(&h)) - 8*sizeof(void*), 4096);
      ss << "Map head:\n";
      CxxUtils::safeHexdump (ss, &binsInLayers, sizeof(binsInLayers));
      ATH_MSG_INFO(ss.str());
      ss.str("");
      ss << "First node:\n";
      auto it = binsInLayers.begin();
      CxxUtils::safeHexdump (ss, (reinterpret_cast<const char*>(&*it)) - 8*sizeof(void*), 4096);
      ATH_MSG_INFO(ss.str());
      ATH_MSG_INFO("zombie: " << h.IsZombie() <<
                   " on heap: " << h.IsOnHeap());
      ATH_MSG_INFO(" dynamic type: " << typeid(h).name());
 
      ATH_MSG_INFO("Got truth state: ");
      truth->Print();
 
      ATH_MSG_INFO("Got extrapolation state: ");
      extrapol->Print();
 
      ATH_MSG_INFO("Got simulation state: ");
      simulstate.Print();
 
      ATH_MSG_INFO("Got GAN XML parameters: ");
      m_param.Print();
 
      return false;
    }
 
    const int xBinNum = h.GetNbinsX();
    const int yBinNum = h.GetNbinsY();
    const TAxis *x = h.GetXaxis();
 
    // If only one bin in r means layer is empty, no value should be added
    if (xBinNum == 1) {
      ATH_MSG_VERBOSE(" Layer "
                      << layer
                      << " has only one bin in r, this means is it not used, "
                         "skipping (this is needed to keep correct "
                         "syncronisation of voxel and layers)");
      // delete h;
      continue;
    }
 
    ATH_MSG_VERBOSE(" Getting energy for Layer " << layer);
 
    // First fill energies
    for (int ix = 1; ix <= xBinNum; ++ix) {
      double binsInAlphaInRBin = GetAlphaBinsForRBin(x, ix, yBinNum);
      for (int iy = 1; iy <= binsInAlphaInRBin; ++iy) {
        const double energyInVoxel = outputs.at(std::to_string(vox));
        ATH_MSG_VERBOSE(" Vox " << vox << " energy " << energyInVoxel
                                << " binx " << ix << " biny " << iy);
 
        if (energyInVoxel <= 0) {
          vox++;
          continue;
        }
 
        simulstate.add_E(layer, Einit * energyInVoxel);
        vox++;
      }
    }
  }
 
  for (unsigned int ichain = m_bin_start.back(); ichain < size(); ++ichain) {
    ATH_MSG_DEBUG("now run for all bins: " << chain()[ichain]->GetName());
    if (simulate_and_retry(chain()[ichain], simulstate, truth, extrapol) !=
        FCSSuccess) {
      return FCSFatal;
    }
  }
 
  vox = 0;
  for (const auto &[layer, h] : binsInLayers) {
    const int xBinNum = h.GetNbinsX();
    const int yBinNum = h.GetNbinsY();
    const TAxis *x = h.GetXaxis();
    const TAxis *y = h.GetYaxis();
 
    simulstate.setAuxInfo<int>("GANlayer"_FCShash, layer);
    TFCSLateralShapeParametrizationHitBase::Hit hit;
 
    // If only one bin in r means layer is empty, no value should be added
    if (xBinNum == 1) {
      ATH_MSG_VERBOSE(" Layer "
                      << layer
                      << " has only one bin in r, this means is it not used, "
                         "skipping (this is needed to keep correct "
                         "syncronisation of voxel and layers)");
      // delete h;
      continue;
    }
 
    if (get_number_of_bins() > 0) {
      const int bin = get_bin(simulstate, truth, extrapol);
      if (bin >= 0 && bin < (int)get_number_of_bins()) {
        for (unsigned int ichain = m_bin_start[bin];
             ichain < TMath::Min(m_bin_start[bin] + get_nr_of_init(bin),
                                 m_bin_start[bin + 1]);
             ++ichain) {
          ATH_MSG_DEBUG("for " << get_variable_text(simulstate, truth, extrapol)
                               << " run init " << get_bin_text(bin) << ": "
                               << chain()[ichain]->GetName());
          if (chain()[ichain]->InheritsFrom(
                  TFCSLateralShapeParametrizationHitBase::Class())) {
            TFCSLateralShapeParametrizationHitBase *sim =
                static_cast<TFCSLateralShapeParametrizationHitBase *>(chain()[ichain]);
            if (sim->simulate_hit(hit, simulstate, truth, extrapol) !=
                FCSSuccess) {
              ATH_MSG_ERROR("error for "
                            << get_variable_text(simulstate, truth, extrapol)
                            << " run init " << get_bin_text(bin) << ": "
                            << chain()[ichain]->GetName());
              return false;
            }
          } else {
            ATH_MSG_ERROR("for "
                          << get_variable_text(simulstate, truth, extrapol)
                          << " run init " << get_bin_text(bin) << ": "
                          << chain()[ichain]->GetName()
                          << " does not inherit from "
                             "TFCSLateralShapeParametrizationHitBase");
            return false;
          }
        }
      } else {
        ATH_MSG_WARNING("nothing to init for "
                        << get_variable_text(simulstate, truth, extrapol)
                        << ": " << get_bin_text(bin));
      }
    }
 
    int binResolution = 5;
    if (layer == 1 || layer == 5) {
      binResolution = 1;
    }
 
    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 << " Extrap eta " << center_eta
                              << " phi " << center_phi << " R " << center_r);
 
    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)));
 
    int nHitsAlpha{};
    int nHitsR{};
 
    // Now create hits
    for (int ix = 1; ix <= xBinNum; ++ix) {
      const int binsInAlphaInRBin = GetAlphaBinsForRBin(x, ix, yBinNum);
 
      // Horrible work around for variable # of bins along alpha direction
      const int binsToMerge = yBinNum == 32 ? 32 / binsInAlphaInRBin : 1;
      for (int iy = 1; iy <= binsInAlphaInRBin; ++iy) {
        const double energyInVoxel = outputs.at(std::to_string(vox));
        const int lowEdgeIndex = (iy - 1) * binsToMerge + 1;
 
        ATH_MSG_VERBOSE(" Vox " << vox << " energy " << energyInVoxel
                                << " binx " << ix << " biny " << iy);
 
        if (energyInVoxel <= 0) {
          vox++;
          continue;
        }
 
        if (std::abs(pdgId) == 22 || std::abs(pdgId) == 11) {
          // maximum 10 MeV per hit, equaly distributed in alpha and r
          int maxHitsInVoxel = energyInVoxel * truth->Ekin() / 10;
          if (maxHitsInVoxel < 1)
            maxHitsInVoxel = 1;
          nHitsAlpha = std::sqrt(maxHitsInVoxel);
          nHitsR = std::sqrt(maxHitsInVoxel);
        } else {
          // One hit per mm along r
          nHitsR = x->GetBinUpEdge(ix) - x->GetBinLowEdge(ix);
          if (yBinNum == 1) {
            // nbins in alpha depend on circumference lenght
            const double r = x->GetBinUpEdge(ix);
            nHitsAlpha = ceil(2 * TMath::Pi() * r / binResolution);
          } else {
            // d = 2*r*sin (a/2r) this distance at the upper r must be 1mm for
            // layer 1 or 5, 5mm otherwise.
            const double angle = y->GetBinUpEdge(iy) - y->GetBinLowEdge(iy);
            const double r = x->GetBinUpEdge(ix);
            const double d = 2 * r * sin(angle / 2 * r);
            nHitsAlpha = ceil(d / binResolution);
          }
 
          if (layer != 1 && layer != 5) {
            // For layers that are not EMB1 or EMEC1 use a maximum of 10 hits
            // per direction, a higher granularity is needed for the other
            // layers
            const int maxNhits = 10;
            nHitsAlpha = std::min(maxNhits, std::max(1, nHitsAlpha));
            nHitsR = std::min(maxNhits, std::max(1, nHitsR));
          }
        }
 
        for (int ir = 0; ir < nHitsR; ++ir) {
          double r =
              x->GetBinLowEdge(ix) + x->GetBinWidth(ix) / (nHitsR + 1) * ir;
 
          for (int ialpha = 1; ialpha <= nHitsAlpha; ++ialpha) {
            if (ganVersion > 1) {
              if (fitResults.at(layer)[ix - 1] != 0) {
                int tries = 0;
                double a = CLHEP::RandFlat::shoot(simulstate.randomEngine(),
                                                  x->GetBinLowEdge(ix),
                                                  x->GetBinUpEdge(ix));
                double rand_r =
                    log((a - x->GetBinLowEdge(ix)) / (x->GetBinWidth(ix))) /
                    fitResults.at(layer)[ix - 1];
                while ((rand_r < x->GetBinLowEdge(ix) ||
                        rand_r > x->GetBinUpEdge(ix)) &&
                       tries < 100) {
                  a = CLHEP::RandFlat::shoot(simulstate.randomEngine(),
                                             x->GetBinLowEdge(ix),
                                             x->GetBinUpEdge(ix));
                  rand_r =
                      log((a - x->GetBinLowEdge(ix)) / (x->GetBinWidth(ix))) /
                      fitResults.at(layer)[ix - 1];
                  tries++;
                }
                if (tries >= 100) {
                  ATH_MSG_VERBOSE(" Too many tries for bin ["
                                  << x->GetBinLowEdge(ix) << "-"
                                  << x->GetBinUpEdge(ix) << "] having slope "
                                  << fitResults.at(layer)[ix - 1]
                                  << " will use grid (old method)");
                } else {
                  r = rand_r;
                }
              }
            }
 
            double alpha;
            if (binsInAlphaInRBin == 1) {
              alpha = CLHEP::RandFlat::shoot(simulstate.randomEngine(),
                                             -TMath::Pi(), TMath::Pi());
            } else {
              alpha =
                  y->GetBinLowEdge(lowEdgeIndex) +
                  y->GetBinWidth(iy) * binsToMerge / (nHitsAlpha + 1) * ialpha;
 
              if (m_param.IsSymmetrisedAlpha()) {
                if (CLHEP::RandFlat::shoot(simulstate.randomEngine(),
                                           -TMath::Pi(), TMath::Pi()) < 0) {
                  alpha = -alpha;
                }
              }
            }
 
            hit.reset();
            hit.E() = Einit * energyInVoxel / (nHitsAlpha * nHitsR);
 
            if (layer <= 20) {
              float delta_eta_mm = r * cos(alpha);
              float delta_phi_mm = r * sin(alpha);
 
              ATH_MSG_VERBOSE("delta_eta_mm " << delta_eta_mm
                                              << " delta_phi_mm "
                                              << delta_phi_mm);
 
              // 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);
            } else {  // FCAL is in (x,y,z)
              const float hit_r = r * cos(alpha) + center_r;
              float delta_phi = r * 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 * cos(hit_phi);
              hit.y() = hit_r * sin(hit_phi);
              hit.z() = center_z;
              ATH_MSG_VERBOSE(" Hit x " << hit.x() << " y " << hit.y()
                                        << " layer " << layer);
            }
 
            if (get_number_of_bins() > 0) {
              const int bin = get_bin(simulstate, truth, extrapol);
              if (bin >= 0 && bin < (int)get_number_of_bins()) {
                for (unsigned int ichain =
                         m_bin_start[bin] + get_nr_of_init(bin);
                     ichain < m_bin_start[bin + 1]; ++ichain) {
                  ATH_MSG_DEBUG(
                      "for " << get_variable_text(simulstate, truth, extrapol)
                             << " run " << get_bin_text(bin) << ": "
                             << chain()[ichain]->GetName());
                  if (chain()[ichain]->InheritsFrom(
                          TFCSLateralShapeParametrizationHitBase::Class())) {
                    TFCSLateralShapeParametrizationHitBase *sim =
                        static_cast<TFCSLateralShapeParametrizationHitBase
                             *>(chain()[ichain]);
                    if (sim->simulate_hit(hit, simulstate, truth, extrapol) !=
                        FCSSuccess) {
                      ATH_MSG_ERROR(
                          "error for "
                          << get_variable_text(simulstate, truth, extrapol)
                          << " run init " << get_bin_text(bin) << ": "
                          << chain()[ichain]->GetName());
                      return false;
                    }
                  } else {
                    ATH_MSG_ERROR(
                        "for " << get_variable_text(simulstate, truth, extrapol)
                               << " run init " << get_bin_text(bin) << ": "
                               << chain()[ichain]->GetName()
                               << " does not inherit from "
                                  "TFCSLateralShapeParametrizationHitBase");
                    return false;
                  }
                }
              } else {
                ATH_MSG_WARNING(
                    "nothing to do for "
                    << get_variable_text(simulstate, truth, extrapol) << ": "
                    << get_bin_text(bin));
              }
            } else {
              ATH_MSG_WARNING("no bins defined, is this intended?");
            }
          }
        }
        vox++;
      }
    }
 
    ATH_MSG_VERBOSE("Number of voxels " << vox);
    ATH_MSG_VERBOSE("Done layer " << layer);
  }
 
  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());
    }
  }
 
  ATH_MSG_VERBOSE("Done particle");
  return true;
}

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

GANfreemem()

bool TFCSEnergyAndHitGANV2::GANfreemem ( ) const

inline

Definition at line 40 of file TFCSEnergyAndHitGANV2.h.

{ return TestBit(kGANfreemem); };

get_bin()

virtual int TFCSEnergyAndHitGANV2::get_bin ( TFCSSimulationState & simulstate, const TFCSTruthState * , const TFCSExtrapolationState *  ) const

inlineoverridevirtual

use the layer to be done as binning of the GAN chain

Reimplemented from TFCSParametrizationBinnedChain.

Definition at line 56 of file TFCSEnergyAndHitGANV2.h.

                                                                     {
    return simulstate.getAuxInfo<int>("GANlayer"_FCShash);
  };

get_bin_text()

const std::string TFCSParametrizationBinnedChain::get_bin_text ( int bin ) const

virtualinherited

print the range of a bin; for bin -1, print the allowed range

Reimplemented in TFCSParametrizationAbsEtaSelectChain, TFCSParametrizationEkinSelectChain, and TFCSParametrizationEtaSelectChain.

Definition at line 60 of file TFCSParametrizationBinnedChain.cxx.

                                                                          {
  return std::string(Form("bin %d", bin));
}

get_Binning()

const TFCSGANXMLParameters::Binning & TFCSEnergyAndHitGANV2::get_Binning ( ) const

inline

Definition at line 67 of file TFCSEnergyAndHitGANV2.h.

                                                         {
    return m_param.GetBinning();
  };

get_ExtrapolationWeights()

const TFCSGANEtaSlice::ExtrapolatorWeights TFCSEnergyAndHitGANV2::get_ExtrapolationWeights ( ) const

inline

Definition at line 70 of file TFCSEnergyAndHitGANV2.h.

                                                                            {
    return m_slice->GetExtrapolatorWeights();
  };

get_nr_of_init()

unsigned int TFCSEnergyAndHitGANV2::get_nr_of_init

(

unsigned int

bin

)

const

Definition at line 56 of file TFCSEnergyAndHitGANV2.cxx.

                                                                         {
  if (bin >= m_bin_ninit.size())
    return 0;
  return m_bin_ninit[bin];
}

get_number_of_bins()

virtual unsigned int TFCSParametrizationBinnedChain::get_number_of_bins ( ) const

inlinevirtualinherited

Definition at line 22 of file TFCSParametrizationBinnedChain.h.

                                                  {
    return m_bin_start.size() - 1;
  };

get_variable_text()

const std::string TFCSEnergyAndHitGANV2::get_variable_text ( TFCSSimulationState & simulstate, const TFCSTruthState * , const TFCSExtrapolationState *  ) const

overridevirtual

Reimplemented from TFCSParametrizationBinnedChain.

Definition at line 107 of file TFCSEnergyAndHitGANV2.cxx.

                                          {
  return std::string(
      Form("layer=%d", simulstate.getAuxInfo<int>("GANlayer"_FCShash)));
}

GetAlphaBinsForRBin()

int TFCSEnergyAndHitGANV2::GetAlphaBinsForRBin ( const TAxis * x, int ix, int yBinNum ) const

private

Definition at line 729 of file TFCSEnergyAndHitGANV2.cxx.

                                                                  {
  double binsInAlphaInRBin = yBinNum;
  if (yBinNum == 32) {
    ATH_MSG_DEBUG("yBinNum is special value 32");
    const double widthX = x->GetBinWidth(ix);
    const double radious = x->GetBinCenter(ix);
    double circumference = radious * 2. * TMath::Pi();
    if (m_param.IsSymmetrisedAlpha()) {
      circumference = radious * TMath::Pi();
    }
 
    const double bins = circumference / widthX;
    binsInAlphaInRBin = GetBinsInFours(bins);
    ATH_MSG_DEBUG("Bin in alpha: " << binsInAlphaInRBin << " for r bin: " << ix
                                   << " (" << x->GetBinLowEdge(ix) << "-"
                                   << x->GetBinUpEdge(ix) << ")");
  }
  return binsInAlphaInRBin;
}

GetBinsInFours()

int TFCSEnergyAndHitGANV2::GetBinsInFours ( double const bins )

staticprivate

Definition at line 718 of file TFCSEnergyAndHitGANV2.cxx.

                                                           {
  if (bins < 4)
    return 4;
  else if (bins < 8)
    return 8;
  else if (bins < 16)
    return 16;
  else
    return 32;
}

initializeNetwork()

bool TFCSEnergyAndHitGANV2::initializeNetwork	(	int const &	pid,
		int const &	etaMin,
		const std::string &	FastCaloGANInputFolderName )

Definition at line 72 of file TFCSEnergyAndHitGANV2.cxx.

                                                 {
 
  // initialize all necessary constants
  // FIXME eventually all these could be stored in the .json file
 
  ATH_MSG_INFO(
      "Using FastCaloGANInputFolderName: " << FastCaloGANInputFolderName);
  // get neural net JSON file as an std::istream object
  const int etaMax = etaMin + 5;
 
  reset_match_all_pdgid();
  set_pdgid(pid);
  if (pid == 11)
    add_pdgid(-pid);
  if (pid == 211)
    add_pdgid(-pid);
  set_eta_min(etaMin / 100.0);
  set_eta_max(etaMax / 100.0);
  set_eta_nominal((etaMin + etaMax) / 200.0);
 
  int pidForXml = pid;
  if (pid != 22 && pid != 11) {
    pidForXml = 211;
  }
 
  const int etaMid = (etaMin + etaMax) / 2;
  m_param.InitialiseFromXML(pidForXml, etaMid, FastCaloGANInputFolderName);
  m_param.Print();
  m_slice = new TFCSGANEtaSlice(pid, etaMin, etaMax, m_param);
  m_slice->Print();
  return m_slice->LoadGAN();
}

is_match_all_calosample()

virtual bool TFCSEnergyAndHitGANV2::is_match_all_calosample ( ) const

inlineoverridevirtual

Reimplemented from TFCSParametrization.

Definition at line 32 of file TFCSEnergyAndHitGANV2.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 TFCSEnergyAndHitGANV2::is_match_all_Ekin_bin ( ) const

inlineoverridevirtual

Reimplemented from TFCSParametrization.

Definition at line 31 of file TFCSEnergyAndHitGANV2.h.

{ return true; };

is_match_all_eta()

virtual bool TFCSParametrization::is_match_all_eta ( ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 28 of file TFCSParametrization.h.

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

is_match_all_pdgid()

virtual bool TFCSParametrizationBase::is_match_all_pdgid ( ) const

inlinevirtualinherited

Definition at line 66 of file TFCSParametrizationBase.h.

{ return TestBit(kMatchAllPDGID); };

is_match_calosample()

bool TFCSEnergyAndHitGANV2::is_match_calosample ( int calosample ) const

overridevirtual

Reimplemented from TFCSParametrizationChain.

Definition at line 48 of file TFCSEnergyAndHitGANV2.cxx.

                                                                    {
  if (get_Binning().find(calosample) == get_Binning().cend())
    return false;
  if (get_Binning().at(calosample).GetNbinsX() == 1)
    return false;
  return true;
}

is_match_Ekin()

virtual bool TFCSParametrization::is_match_Ekin ( float Ekin ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 18 of file TFCSParametrization.h.

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

is_match_Ekin_bin()

virtual bool TFCSEnergyAndHitGANV2::is_match_Ekin_bin ( int ) const

inlineoverridevirtual

Reimplemented from TFCSParametrizationChain.

Definition at line 27 of file TFCSEnergyAndHitGANV2.h.

                                                       {
    return true;
  };

is_match_eta()

virtual bool TFCSParametrization::is_match_eta ( float eta ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 21 of file TFCSParametrization.h.

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

is_match_pdgid()

virtual bool TFCSParametrization::is_match_pdgid ( int id ) const

inlineoverridevirtualinherited

Reimplemented from TFCSParametrizationBase.

Definition at line 15 of file TFCSParametrization.h.

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

level()

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

inlineinherited

Retrieve output level.

Definition at line 201 of file MLogging.h.

{ return msg().level(); }

msg() [1/2]

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

inlineinherited

Return a stream for sending messages directly (no decoration).

Definition at line 231 of file MLogging.h.

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

msg() [2/2]

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

inlineinherited

Return a decorated starting stream for sending messages.

Definition at line 240 of file MLogging.h.

                                                        {
  return msg() << lvl;
}

msgLvl()

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

inlineinherited

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

Definition at line 222 of file MLogging.h.

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

OnlyScaleEnergy()

bool TFCSEnergyAndHitGANV2::OnlyScaleEnergy ( ) const

inline

Definition at line 51 of file TFCSEnergyAndHitGANV2.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 const TFCSParametrizationBase * TFCSParametrizationChain::operator[] ( unsigned int ) const

inlineoverridevirtualinherited

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 43 of file TFCSParametrizationChain.h.

                                              {
    return m_chain[ind];
  };

operator[]() [2/2]

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

inlineoverridevirtualinherited

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 46 of file TFCSParametrizationChain.h.

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

overridevirtual

Definition at line 572 of file TFCSEnergyAndHitGANV2.cxx.

                                                        {
  TFCSParametrization::Print(option);
  TString opt(option);
  const bool shortprint = opt.Index("short") >= 0;
  const bool longprint =
      msgLvl(MSG::DEBUG) || (msgLvl(MSG::INFO) && !shortprint);
  TString optprint = opt;
  optprint.ReplaceAll("short", "");
 
  TString prefix = "- ";
  for (unsigned int ichain = 0; ichain < size(); ++ichain) {
    if (ichain == 0 && ichain != m_bin_start.front()) {
      prefix = "> ";
      if (longprint)
        ATH_MSG_INFO(optprint << prefix << "Run for all bins");
    }
    for (unsigned int ibin = 0; ibin < get_number_of_bins(); ++ibin) {
      if (ichain == m_bin_start[ibin]) {
        if (ibin < get_number_of_bins() - 1)
          if (ichain == m_bin_start[ibin + 1])
            continue;
        prefix = Form("%-2d", ibin);
        if (longprint)
          ATH_MSG_INFO(optprint << prefix << "Run for " << get_bin_text(ibin));
      }
    }
    if (ichain == m_bin_start.back()) {
      prefix = "< ";
      if (longprint)
        ATH_MSG_INFO(optprint << prefix << "Run for all bins");
    }
    chain()[ichain]->Print(opt + prefix);
  }
}

push_back()

void TFCSParametrizationChain::push_back ( const Chain_t::value_type & param )

inlineinherited

Definition at line 55 of file TFCSParametrizationChain.h.

                                                 {
    m_chain.push_back(param);
    recalc();
  };

push_back_in_bin()

void TFCSParametrizationBinnedChain::push_back_in_bin ( TFCSParametrizationBase * param, unsigned int bin )

virtualinherited

Reimplemented in TFCSParametrizationEkinSelectChain, TFCSParametrizationEtaSelectChain, and TFCSParametrizationFloatSelectChain.

Definition at line 31 of file TFCSParametrizationBinnedChain.cxx.

                                                      {
  if (m_bin_start.size() <= bin + 1) {
    m_bin_start.resize(bin + 2, m_bin_start.back());
    m_bin_start.shrink_to_fit();
  }
  if (m_bin_start[bin + 1] == size()) {
    chain().push_back(param);
  } else {
    Chain_t::iterator it(&chain()[m_bin_start[bin + 1]]);
    chain().insert(it, param);
  }
  for (unsigned int ibin = bin + 1; ibin < m_bin_start.size(); ++ibin)
    m_bin_start[ibin]++;
  recalc();
}

push_before_first_bin()

void TFCSParametrizationBinnedChain::push_before_first_bin ( TFCSParametrizationBase * param )

virtualinherited

Definition at line 18 of file TFCSParametrizationBinnedChain.cxx.

                                    {
  if (m_bin_start[0] == size()) {
    chain().push_back(param);
  } else {
    Chain_t::iterator it(&chain()[m_bin_start[0]]);
    chain().insert(it, param);
  }
  for (auto &ele : m_bin_start)
    ele++;
  recalc();
}

recalc()

void TFCSParametrizationChain::recalc ( )

protectedvirtualinherited

Default is to call recalc_pdgid_intersect() and recalc_Ekin_eta_intersect().

Reimplemented in TFCSParametrizationEkinSelectChain, TFCSParametrizationEtaSelectChain, and TFCSParametrizationPDGIDSelectChain.

Definition at line 119 of file TFCSParametrizationChain.cxx.

                                      {
  clear();
  if (m_chain.empty())
    return;
 
  recalc_pdgid_intersect();
  recalc_Ekin_eta_intersect();
 
  m_chain.shrink_to_fit();
}

recalc_Ekin_eta_intersect()

void TFCSParametrizationChain::recalc_Ekin_eta_intersect ( )

protectedinherited

Definition at line 77 of file TFCSParametrizationChain.cxx.

                                                         {
  recalc_Ekin_intersect();
  recalc_eta_intersect();
}

recalc_Ekin_eta_union()

void TFCSParametrizationChain::recalc_Ekin_eta_union ( )

protectedinherited

Definition at line 114 of file TFCSParametrizationChain.cxx.

                                                     {
  recalc_Ekin_union();
  recalc_eta_union();
}

recalc_Ekin_intersect()

void TFCSParametrizationChain::recalc_Ekin_intersect ( )

protectedinherited

Definition at line 45 of file TFCSParametrizationChain.cxx.

                                                     {
  set_Ekin(*m_chain[0]);
 
  for (const auto &param : m_chain) {
    if (param->Ekin_min() > Ekin_min())
      set_Ekin_min(param->Ekin_min());
    if (param->Ekin_max() < Ekin_max())
      set_Ekin_max(param->Ekin_max());
    if (Ekin_nominal() < Ekin_min() || Ekin_nominal() > Ekin_max())
      set_Ekin_nominal(param->Ekin_nominal());
  }
 
  if (Ekin_nominal() < Ekin_min() || Ekin_nominal() > Ekin_max())
    set_Ekin_nominal(0.5 * (Ekin_min() + Ekin_max()));
}

recalc_Ekin_union()

void TFCSParametrizationChain::recalc_Ekin_union ( )

protectedinherited

Definition at line 82 of file TFCSParametrizationChain.cxx.

                                                 {
  set_Ekin(*m_chain[0]);
 
  for (const auto &param : m_chain) {
    if (param->Ekin_min() < Ekin_min())
      set_Ekin_min(param->Ekin_min());
    if (param->Ekin_max() > Ekin_max())
      set_Ekin_max(param->Ekin_max());
    if (Ekin_nominal() < Ekin_min() || Ekin_nominal() > Ekin_max())
      set_Ekin_nominal(param->Ekin_nominal());
  }
 
  if (Ekin_nominal() < Ekin_min() || Ekin_nominal() > Ekin_max())
    set_Ekin_nominal(0.5 * (Ekin_min() + Ekin_max()));
}

recalc_eta_intersect()

void TFCSParametrizationChain::recalc_eta_intersect ( )

protectedinherited

Definition at line 61 of file TFCSParametrizationChain.cxx.

                                                    {
  set_eta(*m_chain[0]);
 
  for (const auto &param : m_chain) {
    if (param->eta_min() > eta_min())
      set_eta_min(param->eta_min());
    if (param->eta_max() < eta_max())
      set_eta_max(param->eta_max());
    if (eta_nominal() < eta_min() || eta_nominal() > eta_max())
      set_eta_nominal(param->eta_nominal());
  }
 
  if (eta_nominal() < eta_min() || eta_nominal() > eta_max())
    set_eta_nominal(0.5 * (eta_min() + eta_max()));
}

recalc_eta_union()

void TFCSParametrizationChain::recalc_eta_union ( )

protectedinherited

Definition at line 98 of file TFCSParametrizationChain.cxx.

                                                {
  set_eta(*m_chain[0]);
 
  for (const auto &param : m_chain) {
    if (param->eta_min() < eta_min())
      set_eta_min(param->eta_min());
    if (param->eta_max() > eta_max())
      set_eta_max(param->eta_max());
    if (eta_nominal() < eta_min() || eta_nominal() > eta_max())
      set_eta_nominal(param->eta_nominal());
  }
 
  if (eta_nominal() < eta_min() || eta_nominal() > eta_max())
    set_eta_nominal(0.5 * (eta_min() + eta_max()));
}

recalc_pdgid_intersect()

void TFCSParametrizationChain::recalc_pdgid_intersect ( )

protectedinherited

Definition at line 20 of file TFCSParametrizationChain.cxx.

                                                      {
  set_pdgid(m_chain[0]->pdgid());
 
  for (const auto &param : m_chain) {
    std::set<int> tmp;
 
    std::set_intersection(pdgid().begin(), pdgid().end(),
                          param->pdgid().begin(), param->pdgid().end(),
                          std::inserter(tmp, tmp.begin()));
    set_pdgid(tmp);
  }
}

recalc_pdgid_union()

void TFCSParametrizationChain::recalc_pdgid_union ( )

protectedinherited

Definition at line 33 of file TFCSParametrizationChain.cxx.

                                                  {
  set_pdgid(chain()[0]->pdgid());
 
  for (const auto &param : chain()) {
    std::set<int> tmp;
 
    std::set_union(pdgid().begin(), pdgid().end(), param->pdgid().begin(),
                   param->pdgid().end(), std::inserter(tmp, tmp.begin()));
    set_pdgid(tmp);
  }
}

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

void TFCSEnergyAndHitGANV2::reset_GANfreemem ( )

inline

Definition at line 42 of file TFCSEnergyAndHitGANV2.h.

{ ResetBit(kGANfreemem); };

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 TFCSEnergyAndHitGANV2::reset_OnlyScaleEnergy ( )

inline

Definition at line 53 of file TFCSEnergyAndHitGANV2.h.

{ ResetBit(kOnlyScaleEnergy); };

reset_RetryChainFromStart()

void TFCSParametrizationChain::reset_RetryChainFromStart ( )

inlineinherited

Definition at line 38 of file TFCSParametrizationChain.h.

{ ResetBit(kRetryChainFromStart); };

reset_SplitChainObjects()

void TFCSParametrizationChain::reset_SplitChainObjects ( )

inlineinherited

Definition at line 34 of file TFCSParametrizationChain.h.

{ ResetBit(kSplitChainObjects); };

RetryChainFromStart()

bool TFCSParametrizationChain::RetryChainFromStart ( ) const

inlineinherited

Definition at line 36 of file TFCSParametrizationChain.h.

{ return TestBit(kRetryChainFromStart); };

set_daughter()

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

inlineoverridevirtualinherited

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 49 of file TFCSParametrizationChain.h.

                                                                     {
    m_chain.at(ind) = param;
  };

set_Ekin()

void TFCSParametrization::set_Ekin ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 55 of file TFCSParametrization.cxx.

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

set_Ekin_eta()

void TFCSParametrization::set_Ekin_eta ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 67 of file TFCSParametrization.cxx.

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

set_Ekin_max()

void TFCSParametrization::set_Ekin_max ( double max )

virtualinherited

Definition at line 45 of file TFCSParametrization.cxx.

{ m_Ekin_max = max; }

set_Ekin_min()

void TFCSParametrization::set_Ekin_min ( double min )

virtualinherited

Definition at line 43 of file TFCSParametrization.cxx.

{ m_Ekin_min = min; }

set_Ekin_nominal()

void TFCSParametrization::set_Ekin_nominal ( double min )

virtualinherited

Definition at line 39 of file TFCSParametrization.cxx.

                                                         {
  m_Ekin_nominal = nominal;
}

set_eta()

void TFCSParametrization::set_eta ( const TFCSParametrizationBase & ref )

virtualinherited

Definition at line 61 of file TFCSParametrization.cxx.

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

set_eta_max()

void TFCSParametrization::set_eta_max ( double max )

virtualinherited

Definition at line 53 of file TFCSParametrization.cxx.

{ m_eta_max = max; }

set_eta_min()

void TFCSParametrization::set_eta_min ( double min )

virtualinherited

Definition at line 51 of file TFCSParametrization.cxx.

{ m_eta_min = min; }

set_eta_nominal()

void TFCSParametrization::set_eta_nominal ( double min )

virtualinherited

Definition at line 47 of file TFCSParametrization.cxx.

                                                        {
  m_eta_nominal = nominal;
}

set_GANfreemem()

void TFCSEnergyAndHitGANV2::set_GANfreemem ( )

inline

Definition at line 41 of file TFCSEnergyAndHitGANV2.h.

{ SetBit(kGANfreemem); };

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 TFCSEnergyAndHitGANV2::set_nr_of_init	(	unsigned int	bin,
		unsigned int	ninit )

Definition at line 62 of file TFCSEnergyAndHitGANV2.cxx.

                                                               {
  if (bin >= m_bin_ninit.size()) {
    m_bin_ninit.resize(bin + 1, 0);
    m_bin_ninit.shrink_to_fit();
  }
  m_bin_ninit[bin] = ninit;
}

set_OnlyScaleEnergy()

void TFCSEnergyAndHitGANV2::set_OnlyScaleEnergy ( )

inline

Definition at line 52 of file TFCSEnergyAndHitGANV2.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_RetryChainFromStart()

void TFCSParametrizationChain::set_RetryChainFromStart ( )

inlineinherited

Definition at line 37 of file TFCSParametrizationChain.h.

{ SetBit(kRetryChainFromStart); };

set_SplitChainObjects()

void TFCSParametrizationChain::set_SplitChainObjects ( )

inlineinherited

Definition at line 33 of file TFCSParametrizationChain.h.

{ SetBit(kSplitChainObjects); };

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

SetRegionAndSliceFromXML()

void TFCSEnergyAndHitGANV2::SetRegionAndSliceFromXML ( int pid, int etaMax, std::string FastCaloGANInputFolderName )

protected

simulate()

FCSReturnCode TFCSEnergyAndHitGANV2::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 TFCSParametrizationBinnedChain.

Definition at line 551 of file TFCSEnergyAndHitGANV2.cxx.

                                                  {
  for (unsigned int ichain = 0; ichain < m_bin_start[0]; ++ichain) {
    ATH_MSG_DEBUG("now run for all bins: " << chain()[ichain]->GetName());
    if (simulate_and_retry(chain()[ichain], simulstate, truth, extrapol) !=
        FCSSuccess) {
      return FCSFatal;
    }
  }
 
  ATH_MSG_VERBOSE("Fill Energies");
  if (!fillEnergy(simulstate, truth, extrapol)) {
    ATH_MSG_WARNING("Could not fill energies ");
    // bail out but do not stop the job
    return FCSSuccess;
  }
 
  return FCSSuccess;
}

simulate_and_retry()

FCSReturnCode TFCSParametrizationChain::simulate_and_retry ( TFCSParametrizationBase * parametrization, TFCSSimulationState & simulstate, const TFCSTruthState * truth, const TFCSExtrapolationState * extrapol ) const

protectedinherited

size()

virtual unsigned int TFCSParametrizationChain::size ( ) const

inlineoverridevirtualinherited

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 41 of file TFCSParametrizationChain.h.

{ return m_chain.size(); };

SplitChainObjects()

bool TFCSParametrizationChain::SplitChainObjects ( ) const

inlineinherited

Definition at line 32 of file TFCSParametrizationChain.h.

{ return TestBit(kSplitChainObjects); };

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

test_path()

void TFCSEnergyAndHitGANV2::test_path ( const std::string & path, TFCSSimulationState * simulstate = nullptr, const TFCSTruthState * truth = nullptr, const TFCSExtrapolationState * extrapol = nullptr, const std::string & outputname = "unnamed", int pid = 211 )

static

Definition at line 625 of file TFCSEnergyAndHitGANV2.cxx.

                                                                            {
  ISF_FCS::MLogging logger;
  ATH_MSG_NOCLASS(logger, "Running test on " << path << std::endl);
  if (!simulstate) {
    simulstate = new TFCSSimulationState();
#if defined(__FastCaloSimStandAlone__)
    simulstate->setRandomEngine(new CLHEP::TRandomEngine());
#else
    simulstate->setRandomEngine(new CLHEP::RanluxEngine());
#endif
  }
  if (!truth) {
    ATH_MSG_NOCLASS(logger, "New particle");
    TFCSTruthState *t = new TFCSTruthState();
    t->SetPtEtaPhiM(65536, 0, 0, ParticleConstants::chargedPionMassInMeV);
    t->set_pdgid(pid);
    truth = t;
  }
  if (!extrapol) {
    TFCSExtrapolationState *e = new TFCSExtrapolationState();
    e->set_IDCaloBoundary_eta(truth->Eta());
    for (int i = 0; i < 24; ++i) {
      e->set_eta(i, TFCSExtrapolationState::SUBPOS_ENT, truth->Eta());
      e->set_eta(i, TFCSExtrapolationState::SUBPOS_EXT, truth->Eta());
      e->set_eta(i, TFCSExtrapolationState::SUBPOS_MID, truth->Eta());
      e->set_phi(i, TFCSExtrapolationState::SUBPOS_ENT, 0);
      e->set_phi(i, TFCSExtrapolationState::SUBPOS_EXT, 0);
      e->set_phi(i, TFCSExtrapolationState::SUBPOS_MID, 0);
      e->set_r(i, TFCSExtrapolationState::SUBPOS_ENT, 1500 + i * 10);
      e->set_r(i, TFCSExtrapolationState::SUBPOS_EXT, 1510 + i * 10);
      e->set_r(i, TFCSExtrapolationState::SUBPOS_MID, 1505 + i * 10);
      e->set_z(i, TFCSExtrapolationState::SUBPOS_ENT, 3500 + i * 10);
      e->set_z(i, TFCSExtrapolationState::SUBPOS_EXT, 3510 + i * 10);
      e->set_z(i, TFCSExtrapolationState::SUBPOS_MID, 3505 + i * 10);
    }
    extrapol = e;
  }
 
  TFCSEnergyAndHitGANV2 GAN("GAN", "GAN");
  GAN.setLevel(MSG::INFO);
  const int etaMin = 20;
  const int etaMax = etaMin + 5;
  ATH_MSG_NOCLASS(logger, "Initialize Networks");
  GAN.initializeNetwork(pid, etaMin, path);
  for (int i = 0; i < 24; ++i)
    if (GAN.is_match_calosample(i)) {
      TFCSCenterPositionCalculation *c = new TFCSCenterPositionCalculation(
          Form("center%d", i), Form("center layer %d", i));
      c->set_calosample(i);
      c->setExtrapWeight(0.5);
      c->setLevel(MSG::INFO);
      c->set_pdgid(pid);
      if (pid == 11)
        c->add_pdgid(-pid);
      if (pid == 211)
        c->add_pdgid(-pid);
      c->set_eta_min(etaMin / 100.0);
      c->set_eta_max(etaMax / 100.0);
      c->set_eta_nominal((etaMin + etaMax) / 200.0);
 
      GAN.push_back_in_bin(c, i);
      GAN.set_nr_of_init(i, 1);
    }
 
  GAN.Print();
 
  ATH_MSG_NOCLASS(logger, "Writing GAN to " << outputname);
  const std::string outname = "FCSGANtest_" + outputname + ".root";
  TFile *fGAN = TFile::Open(outname.c_str(), "recreate");
  fGAN->cd();
  // GAN.Write();
  fGAN->WriteObjectAny(&GAN, "TFCSEnergyAndHitGANV2", "GAN");
 
  fGAN->ls();
  fGAN->Close();
  delete fGAN;
 
  ATH_MSG_NOCLASS(logger, "Open " << outname);
  fGAN = TFile::Open(outname.c_str());
  TFCSEnergyAndHitGANV2 *GAN2 = (TFCSEnergyAndHitGANV2 *)(fGAN->Get("GAN"));
  GAN2->setLevel(MSG::INFO);
  GAN2->Print();
 
  ATH_MSG_NOCLASS(logger, "Before running GAN2->simulate()");
  GAN2->simulate(*simulstate, truth, extrapol);
  simulstate->Print();
  delete fGAN;
}

unit_test()

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

static

Definition at line 607 of file TFCSEnergyAndHitGANV2.cxx.

                                                                              {
  ISF_FCS::MLogging logger;
  ATH_MSG_NOCLASS(logger, "Start lwtnn test" << std::endl);
  std::string path =
      "/eos/atlas/atlascerngroupdisk/proj-simul/AF3_Run3/"
      "InputsToBigParamFiles/FastCaloGANWeightsVer02/";
  test_path(path, simulstate, truth, extrapol, "lwtnn");
 
  ATH_MSG_NOCLASS(logger, "Start onnx test" << std::endl);
  path =
      "/eos/atlas/atlascerngroupdisk/proj-simul/AF3_Run3/"
      "InputsToBigParamFiles/FastCaloGANWeightsONNXVer08/";
  test_path(path, simulstate, truth, extrapol, "onnx");
  ATH_MSG_NOCLASS(logger, "Finish all tests" << std::endl);
}

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_bin_ninit

std::vector<int> TFCSEnergyAndHitGANV2::m_bin_ninit

private

Definition at line 103 of file TFCSEnergyAndHitGANV2.h.

m_bin_start

std::vector<unsigned int> TFCSParametrizationBinnedChain::m_bin_start

protectedinherited

Contains the index where the TFCSParametrizationBase* instances to run for a given bin start.

The last entry of the vector correponds to the index from where on TFCSParametrizationBase* objects should be run again for all bins. This way one can loop over some instances for all bins, then only specific ones for one bin and at the end again over some for all bins

Definition at line 53 of file TFCSParametrizationBinnedChain.h.

m_chain

Chain_t TFCSParametrizationChain::m_chain

privateinherited

Definition at line 99 of file TFCSParametrizationChain.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_mutex

std::mutex TFCSEnergyAndHitGANV2::m_mutex

mutableprivate

Definition at line 112 of file TFCSEnergyAndHitGANV2.h.

m_nm

std::string ISF_FCS::MLogging::m_nm

privateinherited

Message source name.

Definition at line 211 of file MLogging.h.

m_param

TFCSGANXMLParameters TFCSEnergyAndHitGANV2::m_param

private

Definition at line 111 of file TFCSEnergyAndHitGANV2.h.

m_pdgid

std::set<int> TFCSParametrization::m_pdgid

privateinherited

Definition at line 63 of file TFCSParametrization.h.

m_slice

TFCSGANEtaSlice* TFCSEnergyAndHitGANV2::m_slice = nullptr

private

Definition at line 110 of file TFCSEnergyAndHitGANV2.h.

m_writtenBases

std::vector<std::unique_ptr<TFCSParametrizationBase> > TFCSParametrizationChain::m_writtenBases

privateinherited

Definition at line 102 of file TFCSParametrizationChain.h.

---

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

• TFCSEnergyAndHitGANV2.h
• TFCSEnergyAndHitGANV2.cxx