HI::AddEventShape Struct Reference

#include <HIEventShapeSummaryUtils.h>

Collaboration diagram for HI::AddEventShape:

Public Member Functions
	AddEventShape (float ww)
constexpr	AddEventShape ()
void	operator() (xAOD::HIEventShape *slice, const xAOD::HIEventShape *in_slice) const

Public Attributes
bool	match_num_harmonics

Private Attributes
float	m_weight

Detailed Description

Definition at line 29 of file HIEventShapeSummaryUtils.h.

Constructor & Destructor Documentation

AddEventShape() [1/2]

HI::AddEventShape::AddEventShape ( float  ww )

inline

Definition at line 34 of file HIEventShapeSummaryUtils.h.

:  match_num_harmonics(true), m_weight(ww){};

AddEventShape() [2/2]

constexpr HI::AddEventShape::AddEventShape ( )

inlineconstexpr

Definition at line 35 of file HIEventShapeSummaryUtils.h.

: match_num_harmonics(true), m_weight(1){};

Member Function Documentation

operator()()

void HI::AddEventShape::operator()	(	xAOD::HIEventShape *	slice,
		const xAOD::HIEventShape *	in_slice
	)		const

Definition at line 10 of file HIEventShapeSummaryUtils.cxx.

  {
    //update members
    slice->setNCells(slice->nCells()+m_weight*in_slice->nCells());
    slice->setEt(slice->et()+m_weight*in_slice->et());
    slice->setArea(slice->area() + m_weight*in_slice->area());
    slice->setRho(slice->rho() + m_weight*in_slice->rho());
 
    unsigned int nord=slice->etCos().size();
    if(nord!=slice->etSin().size()) throw std::domain_error("Input HIEventShape has unequal n-harmonics for Q_x and Q_y");
 
    //if output shape is empty, copy directly from input
    if(nord==0)
    {
      slice->etCos().assign(in_slice->etCos().begin(),in_slice->etCos().end());
      slice->etSin().assign(in_slice->etSin().begin(),in_slice->etSin().end());
      return;
    }
 
    unsigned int in_nord=in_slice->etCos().size();
    //can only do addition for harmonics in input slice, print warning?
    if(in_nord < nord) nord=in_nord;
    else if(in_nord > nord && match_num_harmonics)
    {
      //print a warning?
      auto itr=in_slice->etCos().begin();
      std::advance(itr,nord);
      slice->etCos().insert(slice->etCos().end(),itr,in_slice->etCos().end());
      itr=in_slice->etSin().begin();
      std::advance(itr,nord);
      slice->etSin().insert(slice->etSin().end(),itr,in_slice->etSin().end());
    }
    //only need to sum to nord, if prev condition was met i>nord set explicitly by insert
    for(unsigned int i=0; i<nord; i++)
    {
      float tmp_cos = slice->etCos().at(i);
      slice->etCos()[i] = tmp_cos + m_weight*in_slice->etCos().at(i);
      float tmp_sin = slice->etSin().at(i);
      slice->etSin()[i] = tmp_sin + m_weight*in_slice->etSin().at(i);
    }
  }

Member Data Documentation

m_weight

float HI::AddEventShape::m_weight

private

Definition at line 39 of file HIEventShapeSummaryUtils.h.

match_num_harmonics

bool HI::AddEventShape::match_num_harmonics

Definition at line 32 of file HIEventShapeSummaryUtils.h.

---

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

• HIEventShapeSummaryUtils.h
• HIEventShapeSummaryUtils.cxx