GlobalSim::Egamma1BDTAlgTool Class Reference

#include <Egamma1BDTAlgTool.h>

Inheritance diagram for GlobalSim::Egamma1BDTAlgTool:

Collaboration diagram for GlobalSim::Egamma1BDTAlgTool:

Public Member Functions
	Egamma1BDTAlgTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~Egamma1BDTAlgTool ()=default
StatusCode	initialize () override
virtual StatusCode	run (const EventContext &ctx) const override
virtual std::string	toString () const override

Private Member Functions
std::vector< double >	combine_phi (const LArStripNeighborhood *) const
std::vector< ap_int< 10 > >	digitize (const std::vector< double > &) const

Private Attributes
Gaudi::Property< bool >	m_enableDump
SG::ReadHandleKey< LArStripNeighborhoodContainer >	m_nbhdContainerReadKey

Static Private Attributes
static constexpr int	s_required_phi_len = 17
static constexpr int	s_combination_len = 18

Detailed Description

Definition at line 26 of file Egamma1BDTAlgTool.h.

Constructor & Destructor Documentation

Egamma1BDTAlgTool()

GlobalSim::Egamma1BDTAlgTool::Egamma1BDTAlgTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 14 of file Egamma1BDTAlgTool.cxx.

                                                 :
    base_class(type, name, parent){
  }

~Egamma1BDTAlgTool()

virtual GlobalSim::Egamma1BDTAlgTool::~Egamma1BDTAlgTool ( )

virtualdefault

Member Function Documentation

combine_phi()

std::vector< double > GlobalSim::Egamma1BDTAlgTool::combine_phi ( const LArStripNeighborhood *  nbhd ) const

private

Definition at line 76 of file Egamma1BDTAlgTool.cxx.

                                                                        {
    auto result = std::vector<double>();
 
    const auto& phi_low = nbhd->phi_low();
    if (phi_low.size() != s_required_phi_len) {return result;}
 
    const auto& phi_center = nbhd->phi_center();
    if (phi_center.size() != s_required_phi_len) {return result;}
 
    
    const auto& phi_high = nbhd->phi_high();
    if (phi_high.size() != s_required_phi_len) {return result;}
 
    result.resize(s_combination_len);
 
    constexpr int c{8};
 
    result.at(0) = phi_center.at(c).m_e;
 
    result.at(1) = std::max(phi_low.at(c).m_e, phi_high.at(c).m_e);
 
    int ri{2};
    for (int diff = 1; diff != 9; ++diff) { 
      result.at(ri) =
    std::max({phi_center.at(c-diff).m_e,
        phi_center.at(c+diff).m_e});
 
      result.at(ri+1) =
    std::max({phi_low.at(c-diff).m_e,
      phi_low.at(c+diff).m_e,
      phi_high.at(c-diff).m_e,
      phi_high.at(c+diff).m_e});
      
      ri += 2;
    }
    
    return result;
  }

digitize()

std::vector< ap_int< 10 > > GlobalSim::Egamma1BDTAlgTool::digitize ( const std::vector< double > &  v ) const

private

Definition at line 117 of file Egamma1BDTAlgTool.cxx.

                                                              {
    auto sf = [](double v) {
      if (v < 0) {return 0.;}
      if (v < 8000) {return v / 31.25;}
      if (v < 40000) {return 192. + v / 125;}
      if (v < 168000) {return 432. + v / 500;}
      if (v < 678000) {return 686. + v / 2000;}
      return 1023.;
    };
 
    auto result = std::vector<ap_int<10>>();
    result.reserve(s_combination_len);
      
    std::transform(std::cbegin(v),
           std::cend(v),
           std::back_inserter(result),
           sf);
           
    return result;
  }

initialize()

StatusCode GlobalSim::Egamma1BDTAlgTool::initialize ( )

override

Definition at line 20 of file Egamma1BDTAlgTool.cxx.

                                           {
       
    CHECK(m_nbhdContainerReadKey.initialize());
 
    return StatusCode::SUCCESS;
  }

run()

StatusCode GlobalSim::Egamma1BDTAlgTool::run ( const EventContext &  ctx ) const

overridevirtual

Definition at line 28 of file Egamma1BDTAlgTool.cxx.

                                                      {
    ATH_MSG_DEBUG("run()");
 
  
    // read in LArStrip neighborhoods from the event store
    // there is one neighborhood per EFex RoI
    auto in =
      SG::ReadHandle<LArStripNeighborhoodContainer>(m_nbhdContainerReadKey,
                            ctx);
    CHECK(in.isValid());
 
    ATH_MSG_DEBUG("read in " << (*in).size() << " neighborhoods");
 
    
    for (const auto& nbhd : *in) {
      auto c_phi = combine_phi(nbhd);
      if (c_phi.empty()) {continue;}  // corner case: not all phi have len 17
      auto input = digitize(c_phi);
 
      assert(input.size() == n_features);
      ap_int<10>* c_input = &input[0];  // vector->array
 
      score_t scores[GlobalSim::BDT::fn_classes(n_classes)];
 
      // the bdt variable is already set up 
      bdt.decision_function(c_input, scores);
      {
    std::stringstream ss;
    ss << "BDT input: ";
    for (const auto& i : input) {ss << i << ' ';}
    ATH_MSG_DEBUG(ss.str());
      }
 
      {
    std::stringstream ss;
    ss << "C BDT output: ";
    for (const auto& i : scores) {ss << i << ' ';}
    ATH_MSG_DEBUG(ss.str());
      }
 
    }
 
   
    return StatusCode::SUCCESS;
  }

toString()

std::string GlobalSim::Egamma1BDTAlgTool::toString ( ) const

overridevirtual

Definition at line 139 of file Egamma1BDTAlgTool.cxx.

                                              {
 
    std::stringstream ss;
    ss << "Egamma1BDTAlgTool. name: " << name() << '\n'
       << m_nbhdContainerReadKey << '\n'
       << '\n';
    return ss.str();
  }

Member Data Documentation

m_enableDump

Gaudi::Property<bool> GlobalSim::Egamma1BDTAlgTool::m_enableDump

private

Initial value:

{this,
         "enableDump",
         {false},
         "flag to enable dumps"}

Definition at line 45 of file Egamma1BDTAlgTool.h.

m_nbhdContainerReadKey

SG::ReadHandleKey<LArStripNeighborhoodContainer> GlobalSim::Egamma1BDTAlgTool::m_nbhdContainerReadKey

private

Initial value:

{
      this,
      "LArNeighborhoodContainerReadKey",
      "stripNeighborhoodContainer",
      "key to read inLArNeighborhoodReadKeys"}

Definition at line 52 of file Egamma1BDTAlgTool.h.

s_combination_len

constexpr int GlobalSim::Egamma1BDTAlgTool::s_combination_len = 18

inlinestaticconstexprprivate

Definition at line 68 of file Egamma1BDTAlgTool.h.

s_required_phi_len

constexpr int GlobalSim::Egamma1BDTAlgTool::s_required_phi_len = 17

inlinestaticconstexprprivate

Definition at line 64 of file Egamma1BDTAlgTool.h.

---

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

• Egamma1BDTAlgTool.h
• Egamma1BDTAlgTool.cxx