TCS::cTauMultiplicity Class Reference

#include <cTauMultiplicity.h>

Inheritance diagram for TCS::cTauMultiplicity:

Collaboration diagram for TCS::cTauMultiplicity:

Public Member Functions
	cTauMultiplicity (const std::string &name)
virtual	~cTauMultiplicity ()=default
virtual StatusCode	initialize () override
virtual StatusCode	processBitCorrect (const TCS::InputTOBArray &input, Count &count) override final
virtual StatusCode	process (const TCS::InputTOBArray &input, Count &count) override final
unsigned int	numberOutputBits () const
void	setNumberOutputBits (unsigned int numberOutputBits)
void	setThreshold (const TrigConf::L1Threshold &thr)
const TrigConf::L1Threshold *	getThreshold ()
void	setClassName (const std::string &className)
void	setLegacyMode (bool isLegacyTopo)
const std::string &	name () const
const std::string &	className () const
std::string	fullname () const
unsigned int	algoId () const
StatusCode	reset ()
bool	isSortingAlg () const
bool	isDecisionAlg () const
bool	isCountingAlg () const
bool	isLegacyTopo () const
unsigned int	calcDeltaPhiBW (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcDeltaEtaBW (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcInvMassBW (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcTMassBW (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcDeltaR2BW (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned long	quadraticSumBW (int i1, int i2)
unsigned int	calcDeltaPhi (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcDeltaEta (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcInvMass (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcTMass (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
unsigned int	calcDeltaR2 (const TCS::GenericTOB *tob1, const TCS::GenericTOB *tob2)
const Parameter &	parameter (const std::string &parameterName) const
const Parameter &	parameter (const std::string &parameterName, unsigned int selection) const
const ParameterSpace &	parameters () const
void	setParameters (const ParameterSpace &)
void	setAlgoId (unsigned int algoId)
void	setL1TopoHistSvc (std::shared_ptr< IL1TopoHistSvc >)
	sets the external hist service
void	bookHistMult (std::vector< std::string > &regName, const std::string &name, const std::string &title, const std::string &xtitle, const int binx, const int xmin, const int xmax)
void	bookHistMult (std::vector< std::string > &regName, const std::string &name, const std::string &title, const std::string &xtitle, const std::string &ytitle, const int binx, const int xmin, const int xmax, const int biny, const int ymin, const int ymax)
void	bookHist (std::vector< std::string > &regName, const std::string &name, const std::string &title, const int binx, const int xmin, const int xmax)
void	bookHist (std::vector< std::string > &regName, const std::string &name, const std::string &title, const int binx, const int xmin, const int xmax, const int biny, const int ymin, const int ymax)
bool	msgLvl (const MSGTC::Level lvl) const
	Test the output level.
MsgStreamTC &	msg () const
	The standard message stream.
MsgStreamTC &	msg (const MSGTC::Level lvl) const
	The standard message stream.
const std::string &	getName () const
	name accessor

Static Public Member Functions
static size_t	cTauMatching (const xAOD::eFexTauRoI &eTau, const xAOD::jFexTauRoIContainer &jTauRoIs)
static bool	cTauMatching (const xAOD::eFexTauRoI &eTau, const xAOD::jFexTauRoI &jTau)
static bool	checkIsolationWP (const xAOD::eFexTauRoI &eTau, const xAOD::jFexTauRoI &jTau, const TrigConf::L1Threshold_cTAU &thr)
static bool	checkeTAUWP (const xAOD::eFexTauRoI &eTau, const TrigConf::L1Threshold_cTAU &thr)

Protected Types
enum	AlgType { NONE , SORT , DECISION , COUNT }

Protected Member Functions
void	defineParameter (const std::string &name, TCS::parType_t value)
void	defineParameter (const std::string &name, TCS::parType_t value, unsigned int selection)
void	registerHist (TH1 *)
void	registerHist (TH2 *)
void	fillHist1D (const std::string &histName, double x)
void	fillHist2D (const std::string &histName, double x, double y)
bool	isocut (const std::string &threshold, const unsigned int bit) const
bool	isocut (const unsigned int threshold, const unsigned int bit) const

Protected Attributes
std::vector< std::string >	m_histAccept

Private Types
using	WP = TrigConf::Selection::WP

Private Member Functions
bool	checkIsolationWP (const TCS::cTauTOB *etauCand, const TCS::cTauTOB *jtauCand) const
bool	checkeTAUWP (const TCS::cTauTOB *etauCand) const
bool	cTauMatching (const TCS::cTauTOB *etauCand, const TCS::cTauTOB *jtauCand) const
virtual StatusCode	doReset ()
void	defineParameter (const Parameter &)
std::string	ToString (const int val)

Private Attributes
const TrigConf::L1Threshold_cTAU *	m_threshold {nullptr}
std::shared_ptr< TrigConf::L1ThrExtraInfo_cTAU >	m_extraInfo
std::vector< std::string >	m_histcTauEt
std::vector< std::string >	m_histcTauPhiEta
std::vector< std::string >	m_histcTauEtEta
std::vector< std::string >	m_histcTauIsoFraction
std::vector< std::string >	m_histcTauIsoMatchedPass
unsigned int	m_numberOutputBits {1}
std::unique_ptr< ConfigurableAlgImpl >	m_impl
std::string	m_name {""}
std::string	m_className {""}
unsigned int	m_algoId {0}
ParameterSpace	m_parameters {""}
AlgType	m_algType
bool	m_isLegacyTopo
boost::thread_specific_ptr< MsgStreamTC >	m_msg_tls
	MsgStreamTC instance (a std::cout like with print-out levels)

Detailed Description

Definition at line 28 of file cTauMultiplicity.h.

Member Typedef Documentation

WP

using TCS::cTauMultiplicity::WP = TrigConf::Selection::WP

private

Definition at line 68 of file cTauMultiplicity.h.

Member Enumeration Documentation

AlgType

enum TCS::ConfigurableAlg::AlgType

protectedinherited

Enumerator
NONE
SORT
DECISION
COUNT

Definition at line 32 of file ConfigurableAlg.h.

{ NONE, SORT, DECISION, COUNT };

Constructor & Destructor Documentation

cTauMultiplicity()

TCS::cTauMultiplicity::cTauMultiplicity

(

const std::string &

name

)

Definition at line 23 of file cTauMultiplicity.cxx.

                                                            : CountingAlg(name) {
  setNumberOutputBits(12); //To-Do: Make this flexible to adapt to the menu. Each counting requires more than one bit       
}

~cTauMultiplicity()

virtual TCS::cTauMultiplicity::~cTauMultiplicity ( )

virtualdefault

Member Function Documentation

algoId()

unsigned int TCS::ConfigurableAlg::algoId ( ) const

inlineinherited

Definition at line 53 of file ConfigurableAlg.h.

{ return m_algoId; }

bookHist() [1/2]

void ConfigurableAlg::bookHist ( std::vector< std::string > & regName, const std::string & name, const std::string & title, const int binx, const int xmin, const int xmax )

inherited

Definition at line 270 of file ConfigurableAlg.cxx.

                                                                                                                                                           {
  std::string xmin_str = ToString(xmin);
  std::string xmax_str = ToString(xmax);
  std::string newTitle = title;
  std::string newName = name;
 
  newTitle = xmin_str+title+xmax_str;
  newName = name+"_"+xmin_str+title+xmax_str;
  std::replace( newName.begin(), newName.end(), '-', 'n');
  std::replace( newName.begin(), newName.end(), ' ', '_');
 
  regName.push_back(m_name+"/"+newName);
 
  // Add units to axis labels
  std::string xTitle = title;
  if (m_isLegacyTopo) {
    if (title == "ETA" || title == "DETA" || title == "PHI" || title == "DPHI" || title == "DR") { xTitle = title+"#times10"; }
    if (title == "PT" || title == "ET" || title == "HT" || title == "INVM" || title == "MT")     { xTitle = title+" [GeV]"; } 
  } 
  else {
    if (title == "ETA" || title == "DETA" || title == "DR")                                  { xTitle = title+"#times40"; }
    if (title == "PHI" || title == "DPHI")                                                   { xTitle = title+"#times20"; }
    if (title == "PT" || title == "ET" || title == "HT" || title == "INVM" || title == "MT") { xTitle = title+" [100 MeV]"; } 
  }
 
  int xmin_new,xmax_new,binx_new=binx;
  if ( xmin > 0.0)
    { xmin_new=0.0; }
  else
    { xmin_new=1.5*xmin; }
 
  // if the maximum inv. mass cut defined by the menu over 4000 MeV,
  // set x-axis range maximum to 4000 MeV. This is only for legacy triggers,
  // phase1 does not exceed this limit- asonay 08/2022
  if (xmax > 4000 && m_isLegacyTopo) {
    xmax_new = 4000;
    binx_new = 200;
  }
  else {xmax_new = 1.5*xmax;}
  
  // if angular kinematics, use fixed range
  int eta_max = 50;
  int phi_max = 64;
  int dr_max = 30;
  if (not m_isLegacyTopo) {
    eta_max *= 4;
    phi_max *= 2;
    dr_max *= 4;
  }
 
  if ( title.find("ETA") != std::string::npos ){
    xmin_new=-eta_max;
    xmax_new=eta_max;
  }
  if ( title.find("PHI") != std::string::npos || title.find("DPHI") != std::string::npos ){
    xmin_new=0;
    xmax_new=phi_max;
  }
  if ( title.find("DETA") != std::string::npos || title.find("DR") != std::string::npos ){
    xmin_new=0;
    xmax_new=dr_max;
  }
 
  TH1 *h = new TH1F(newName.c_str(), newTitle.c_str(), binx_new, xmin_new, xmax_new);
  h->GetXaxis()->SetTitle(xTitle.c_str());
  m_impl->registerHist(h);
}

bookHist() [2/2]

void ConfigurableAlg::bookHist ( std::vector< std::string > & regName, const std::string & name, const std::string & title, const int binx, const int xmin, const int xmax, const int biny, const int ymin, const int ymax )

inherited

Definition at line 347 of file ConfigurableAlg.cxx.

                                                                                                                                                                                                           {
  auto usPos = title.find(" vs ");
  std::string xName = title.substr(0,usPos);
  std::string yName = title.substr(usPos+4);
  std::string xmin_str = ToString(xmin);
  std::string xmax_str = ToString(xmax);
  std::string ymin_str = ToString(ymin);
  std::string ymax_str = ToString(ymax);
  std::string newTitle = title;
  std::string newName = name;
 
  newTitle = xmin_str+xName+xmax_str+" vs "+ymin_str+yName+ymax_str;
  newName = name+"_"+xmin_str+xName+xmax_str+"_"+ymin_str+yName+ymax_str;
  std::replace( newName.begin(), newName.end(), '-', 'n');
  std::replace( newName.begin(), newName.end(), ' ', '_');
 
  regName.push_back(m_name+"/"+newName);
 
  // Add units to axis labels
  std::string xTitle = xName;
  std::string yTitle = yName;
  if (m_isLegacyTopo) {
    if (xName == "ETA" || xName == "DETA" || xName == "PHI" || xName == "DPHI" || xName == "DR") { xTitle = xName+"#times10"; }
    if (xName == "PT" || xName == "ET" || xName == "HT" || xName == "INVM" || xName == "MT")     { xTitle = xName+" [GeV]"; } 
 
    if (yName == "ETA" || yName == "DETA" || yName == "PHI" || yName == "DPHI" || yName == "DR") { yTitle = yName+"#times10"; }
    if (yName == "PT" || yName == "ET" || yName == "HT" || yName == "INVM" || yName == "MT")     { yTitle = yName+" [GeV]"; } 
  } 
  else {
    if (xName == "ETA" || xName == "DETA" || xName == "DR")                                  { xTitle = xName+"#times40"; }
    if (xName == "PHI" || xName == "DPHI")                                                   { xTitle = xName+"#times20"; }
    if (xName == "PT" || xName == "ET" || xName == "HT" || xName == "INVM" || xName == "MT") { xTitle = xName+" [100 MeV]"; } 
 
    if (yName == "ETA" || yName == "DETA" || yName == "DR")                                  { yTitle = yName+"#times40"; }
    if (yName == "PHI" || yName == "DPHI")                                                   { yTitle = yName+"#times20"; }
    if (yName == "PT" || yName == "ET" || yName == "HT" || yName == "INVM" || yName == "MT") { yTitle = yName+" [100 MeV]"; } 
  }
 
  int xmin_new,xmax_new,binx_new=binx;
  if ( xmin > 0.0)
    { xmin_new=0.0; }
  else
    { xmin_new=1.5*xmin; }
  
  // if the maximum inv. mass cut defined by the menu over 4000 MeV,
  // set x-axis range maximum to 4000 MeV. This is only for legacy triggers,
  // phase1 does not exceed this limit- asonay 08/2022
  if (xmax > 4000 && m_isLegacyTopo) {
    xmax_new = 4000;
    binx_new = 200;
  }
  else {xmax_new = 1.5*xmax;}
 
  int ymin_new,ymax_new,biny_new=biny;
  if ( ymin > 0.0)
    { ymin_new=0.0; }
  else
    { ymin_new=1.5*ymin; }
  
  // if the maximum inv. mass cut defined by the menu over 4000 MeV,
  // set y-axis range maximum to 4000 MeV. This is only for legacy triggers,
  // phase1 does not exceed this limit- asonay 08/2022
  if (ymax > 4000 && m_isLegacyTopo) {
    ymax_new = 4000;
    biny_new = 200;
  }
  else {ymax_new = 1.5*ymax;}
 
  
  // if angular kinematics, use fixed range
  int eta_max = 50;
  int phi_max = 64;
  int dr_max = 30;
  if (not m_isLegacyTopo) {
    eta_max *= 4;
    phi_max *= 2;
    dr_max *= 4;
  }
 
  if ( xName.find("ETA") != std::string::npos ){
    xmin_new=-eta_max;
    xmax_new=eta_max;
  }
  if ( yName.find("ETA") != std::string::npos ){
    ymin_new=-eta_max;
    ymax_new=eta_max;
  }
  if ( xName.find("PHI") != std::string::npos || xName.find("DPHI") != std::string::npos ){
    xmin_new=0;
    xmax_new=phi_max;
  }
  if ( yName.find("PHI") != std::string::npos || yName.find("DPHI") != std::string::npos ){
    ymin_new=0;
    ymax_new=phi_max;
  }
  if ( xName.find("DETA") != std::string::npos || xName.find("DR") != std::string::npos ){
    xmin_new=0;
    xmax_new=dr_max;
  }
  if ( yName.find("DETA") != std::string::npos || yName.find("DR") != std::string::npos ){
    ymin_new=0;
    ymax_new=dr_max;
  }
 
  TH2 *h = new TH2F(newName.c_str(), newTitle.c_str(), binx_new, xmin_new, xmax_new, biny_new, ymin_new, ymax_new);
  h->GetXaxis()->SetTitle(xTitle.c_str());
  h->GetYaxis()->SetTitle(yTitle.c_str());
  m_impl->registerHist(h);
}

bookHistMult() [1/2]

void ConfigurableAlg::bookHistMult ( std::vector< std::string > & regName, const std::string & name, const std::string & title, const std::string & xtitle, const int binx, const int xmin, const int xmax )

inherited

Definition at line 338 of file ConfigurableAlg.cxx.

                                                                                                                                                                                        {
 
  regName.push_back(m_name+"/"+name);
 
  TH1 *h = new TH1F(name.c_str(), title.c_str(), binx, xmin, xmax);
  h->GetXaxis()->SetTitle(xtitle.c_str());
  m_impl->registerHist(h);
}

bookHistMult() [2/2]

void ConfigurableAlg::bookHistMult ( std::vector< std::string > & regName, const std::string & name, const std::string & title, const std::string & xtitle, const std::string & ytitle, const int binx, const int xmin, const int xmax, const int biny, const int ymin, const int ymax )

inherited

Definition at line 457 of file ConfigurableAlg.cxx.

                                                                                                                                                                                                                                                                   {
 
  regName.push_back(m_name+"/"+name);
 
  TH2 *h = new TH2F(name.c_str(), title.c_str(), binx, xmin, xmax, biny, ymin, ymax);
  h->GetXaxis()->SetTitle(xtitle.c_str());
  h->GetYaxis()->SetTitle(ytitle.c_str());
  m_impl->registerHist(h);
}

calcDeltaEta()

unsigned int ConfigurableAlg::calcDeltaEta ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 173 of file ConfigurableAlg.cxx.

                                                                                  {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcDeltaEtaLegacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcDeltaEta(tob1,tob2);}
}

calcDeltaEtaBW()

unsigned int ConfigurableAlg::calcDeltaEtaBW ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 128 of file ConfigurableAlg.cxx.

                                                                                    {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcDeltaEtaBWLegacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcDeltaEtaBW(tob1,tob2);}
}

calcDeltaPhi()

unsigned int ConfigurableAlg::calcDeltaPhi ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 165 of file ConfigurableAlg.cxx.

                                                                                  {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcDeltaPhiLegacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcDeltaPhi(tob1,tob2);}
}

calcDeltaPhiBW()

unsigned int ConfigurableAlg::calcDeltaPhiBW ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 120 of file ConfigurableAlg.cxx.

                                                                                    {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcDeltaPhiBWLegacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcDeltaPhiBW(tob1,tob2);}
}

calcDeltaR2()

unsigned int ConfigurableAlg::calcDeltaR2 ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 191 of file ConfigurableAlg.cxx.

                                                                                 {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcDeltaR2Legacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcDeltaR2(tob1,tob2);}
}

calcDeltaR2BW()

unsigned int ConfigurableAlg::calcDeltaR2BW ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 152 of file ConfigurableAlg.cxx.

                                                                                   {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcDeltaR2BWLegacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcDeltaR2BW(tob1,tob2);}
}

calcInvMass()

unsigned int ConfigurableAlg::calcInvMass ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 181 of file ConfigurableAlg.cxx.

                                                                                 {
  return TSU::Kinematics::calcInvMass(tob1,tob2);
}

calcInvMassBW()

unsigned int ConfigurableAlg::calcInvMassBW ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 136 of file ConfigurableAlg.cxx.

                                                                                   {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcInvMassBWLegacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcInvMassBW(tob1,tob2);}
}

calcTMass()

unsigned int ConfigurableAlg::calcTMass ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 186 of file ConfigurableAlg.cxx.

                                                                               {
  return TSU::Kinematics::calcTMass(tob1,tob2);
}

calcTMassBW()

unsigned int ConfigurableAlg::calcTMassBW ( const TCS::GenericTOB * tob1, const TCS::GenericTOB * tob2 )

inherited

Definition at line 144 of file ConfigurableAlg.cxx.

                                                                                 {
  if (m_isLegacyTopo)
    {return TSU::Kinematics::calcTMassBWLegacy(tob1,tob2);}
  else
    {return TSU::Kinematics::calcTMassBW(tob1,tob2);}
}

checkeTAUWP() [1/2]

bool TCS::cTauMultiplicity::checkeTAUWP ( const TCS::cTauTOB * etauCand ) const

private

Definition at line 148 of file cTauMultiplicity.cxx.

                                                                      {
  if(m_threshold->isolation() == WP::NONE) return true;
  auto iso_wp = m_extraInfo->isolation(m_threshold->isolation(), etauCand->etaDouble());
  return etauCand->RCore() >= iso_wp.eTAU_rCoreMin_WP_fw() && etauCand->RHad() >= iso_wp.eTAU_rHadMin_WP_fw();
}

checkeTAUWP() [2/2]

bool TCS::cTauMultiplicity::checkeTAUWP ( const xAOD::eFexTauRoI & eTau, const TrigConf::L1Threshold_cTAU & thr )

static

Definition at line 249 of file cTauMultiplicity.cxx.

                                                                                                     {
  if(thr.isolation() == WP::NONE) return true;
  auto iso_wp = thr.getExtraInfo()->isolation(thr.isolation(), eTau.eta());
  return eTau.tauOneThresholds() >= iso_wp.eTAU_rCoreMin_WP_fw() && eTau.tauTwoThresholds() >= iso_wp.eTAU_rHadMin_WP_fw();
}

checkIsolationWP() [1/2]

bool TCS::cTauMultiplicity::checkIsolationWP ( const TCS::cTauTOB * etauCand, const TCS::cTauTOB * jtauCand ) const

private

Definition at line 141 of file cTauMultiplicity.cxx.

                                                                                                       {
  if(m_threshold->isolation() == WP::NONE) return true;
  auto iso_wp = m_extraInfo->isolation(m_threshold->isolation(), etauCand->etaDouble());
  return jtauCand->EtIso()*1024 + jtauCand->Et()*iso_wp.isolation_jTAUCoreScale_fw() < etauCand->Et()*iso_wp.isolation_fw();
}

checkIsolationWP() [2/2]

bool TCS::cTauMultiplicity::checkIsolationWP ( const xAOD::eFexTauRoI & eTau, const xAOD::jFexTauRoI & jTau, const TrigConf::L1Threshold_cTAU & thr )

static

Definition at line 242 of file cTauMultiplicity.cxx.

                                                                                                                                      {
  if(thr.isolation() == WP::NONE) return true;
  auto iso_wp = thr.getExtraInfo()->isolation(thr.isolation(), eTau.eta());
  return jTau.tobIso()*2*1024 + jTau.tobEt()*2*iso_wp.isolation_jTAUCoreScale_fw() < eTau.etTOB()*iso_wp.isolation_fw();
}

className()

const std::string & TCS::ConfigurableAlg::className ( ) const

inlineinherited

Definition at line 49 of file ConfigurableAlg.h.

{ return m_className; }

cTauMatching() [1/3]

bool TCS::cTauMultiplicity::cTauMatching ( const TCS::cTauTOB * etauCand, const TCS::cTauTOB * jtauCand ) const

private

Definition at line 155 of file cTauMultiplicity.cxx.

                                                                                                   {
 
  bool matching  = false; 
 
  // Matching is done comparing eta_tower and phi_tower (granularity = 0.1)
  // These coordinates represent the lower edge of the towers (both for eFEX and jFEX)
 
  // eTau eta = 4*eta_tower + seed, eta from -25 to 24
  int eTauEtaTower;
  if(etauCand->eta()%4 >= 0 ) eTauEtaTower = etauCand->eta() - etauCand->eta()%4;
  else                        eTauEtaTower = etauCand->eta() - etauCand->eta()%4 - 4;
  int jTauEtaTower;
  if (jtauCand->eta()%4 >= 0 ) jTauEtaTower = jtauCand->eta() - jtauCand->eta()%4;
  else                         jTauEtaTower = jtauCand->eta() - jtauCand->eta()%4 - 4;
 
  //int jTauEtaTower = jtauCand->eta();              // jTau eta = 4*eta_tower
  unsigned int eTauPhiTower = etauCand->phi() >> 1;     // eTau phi = 2*phi_tower 
  unsigned int jTauPhiTower = jtauCand->phi() >> 1;     // jTau phi = 2*phi_tower + 1 (jTau coordinates are at center of tower)
 
  matching = (eTauEtaTower == jTauEtaTower) && (eTauPhiTower == jTauPhiTower);
 
  return matching;
 
}

cTauMatching() [2/3]

bool TCS::cTauMultiplicity::cTauMatching ( const xAOD::eFexTauRoI & eTau, const xAOD::jFexTauRoI & jTau )

static

Definition at line 219 of file cTauMultiplicity.cxx.

                                                                                             {
 
  // eFEX: etaTower = iEta, phiTower = iPhi
  // jFEX: etaTower = globalEta, phiTower = globalPhi
  
  int eTauEtaTower;
  if(eTau.iEtaTopo()%4 >= 0 ) eTauEtaTower = eTau.iEtaTopo() - eTau.iEtaTopo()%4;
  else                        eTauEtaTower = eTau.iEtaTopo() - eTau.iEtaTopo()%4 - 4;
  
  int jTauEtaTopo = TSU::toTopoEta(jTau.eta());
  int jTauEtaTower;
  if(jTauEtaTopo%4 >= 0 ) jTauEtaTower = jTauEtaTopo - jTauEtaTopo%4;
  else                    jTauEtaTower = jTauEtaTopo - jTauEtaTopo%4 - 4;
 
  unsigned int jTauPhiTower = TSU::toTopoPhi(jTau.phi()) >> 1; //ignore lowest bit as jTau coordinates are taken at tower center
  unsigned int eTauPhiTower = static_cast<unsigned int>(eTau.iPhiTopo()) >> 1; //shift eTau location in the same way to stay consistent
    
  bool matching = ( eTauEtaTower == jTauEtaTower ) && ( eTauPhiTower == jTauPhiTower );
  return matching;
 
}

cTauMatching() [3/3]

size_t TCS::cTauMultiplicity::cTauMatching ( const xAOD::eFexTauRoI & eTau, const xAOD::jFexTauRoIContainer & jTauRoIs )

static

Definition at line 183 of file cTauMultiplicity.cxx.

                                                                                                            {
 
  // Return the index of the matched jTau if existent (otherwise return std::numeric_limits<size_t>::max())
  size_t i_matched{std::numeric_limits<size_t>::max()};
  size_t i_jTau{0};
  
  int eTauEtaTower;
  if(eTau.iEtaTopo()%4 >= 0 ) eTauEtaTower = eTau.iEtaTopo() - eTau.iEtaTopo()%4;
  else                        eTauEtaTower = eTau.iEtaTopo() - eTau.iEtaTopo()%4 - 4;
  
  for(const xAOD::jFexTauRoI* jTau : jTauRoIs) {
 
    // eFEX: etaTower = iEta, phiTower = iPhi
    // jFEX: etaTower = globalEta, phiTower = globalPhi
    
    int jTauEtaTopo = TSU::toTopoEta(jTau->eta());
    int jTauEtaTower;
    if(jTauEtaTopo%4 >= 0 ) jTauEtaTower = jTauEtaTopo - jTauEtaTopo%4;
    else                    jTauEtaTower = jTauEtaTopo - jTauEtaTopo%4 - 4;
  
    unsigned int jTauPhiTower = TSU::toTopoPhi(jTau->phi()) >> 1; //ignore lowest bit as jTau coordinates are taken at tower center
    unsigned int eTauPhiTower = static_cast<unsigned int>(eTau.iPhiTopo()) >> 1; //shift eTau location in the same way to stay consistent
    bool matching = ( eTauEtaTower == jTauEtaTower ) && ( eTauPhiTower == jTauPhiTower );
 
    if(matching) {
      i_matched = i_jTau;
      break; // Break the loop when a match is found
    }
    ++i_jTau;
  }
 
  return i_matched;
 
}

defineParameter() [1/3]

void ConfigurableAlg::defineParameter ( const Parameter & parameter )

privateinherited

Definition at line 213 of file ConfigurableAlg.cxx.

                                                            {
   if( parameters().contains( parameter.name(), parameter.selection() ) ) {
      TCS_EXCEPTION("Duplicate parameter definition '" << parameter.name() << "' for algorithm '" << name() << "'");
   }
   m_parameters.addParameter(parameter);
}

defineParameter() [2/3]

void ConfigurableAlg::defineParameter ( const std::string & name, TCS::parType_t value )

protectedinherited

Definition at line 201 of file ConfigurableAlg.cxx.

                                                                           {
   defineParameter( Parameter(name, value) );
}

defineParameter() [3/3]

void ConfigurableAlg::defineParameter ( const std::string & name, TCS::parType_t value, unsigned int selection )

protectedinherited

Definition at line 207 of file ConfigurableAlg.cxx.

                                                                                                   {
   defineParameter( Parameter(name, value, selection) );
}

doReset()

TCS::StatusCode TCS::CountingAlg::doReset ( )

privatevirtualinherited

Implements TCS::ConfigurableAlg.

Definition at line 11 of file CountingAlg.cxx.

                        {
   return StatusCode::SUCCESS;
}

fillHist1D()

void ConfigurableAlg::fillHist1D ( const std::string & histName, double x )

protectedinherited

Definition at line 467 of file ConfigurableAlg.cxx.

                                                                     {
  m_impl->fillHist1D(histName,x);
}

fillHist2D()

void ConfigurableAlg::fillHist2D ( const std::string & histName, double x, double y )

protectedinherited

Definition at line 471 of file ConfigurableAlg.cxx.

                                                                               {
  m_impl->fillHist2D(histName,x,y);
}

fullname()

std::string TCS::ConfigurableAlg::fullname ( ) const

inlineinherited

Definition at line 50 of file ConfigurableAlg.h.

{ return m_className + "/" + m_name; }

getName()

const std::string & TrigConf::TrigConfMessaging::getName ( ) const

inlineinherited

name accessor

Returns

the name

Definition at line 101 of file TrigConfMessaging.h.

                                                           {
    return m_name;
  }

getThreshold()

const TrigConf::L1Threshold * TCS::CountingAlg::getThreshold ( )

inlineinherited

Definition at line 47 of file CountingAlg.h.

{ return m_threshold; }

initialize()

TCS::StatusCode TCS::cTauMultiplicity::initialize ( )

overridevirtual

Reimplemented from TCS::ConfigurableAlg.

Definition at line 28 of file cTauMultiplicity.cxx.

                                              {
  m_threshold = dynamic_cast<const TrigConf::L1Threshold_cTAU*>(getThreshold());
  if (not m_threshold){
    TRG_MSG_ERROR("Dynamic cast failed in TCS::cTauMultiplicity::initialize");
    return StatusCode::FAILURE;
  }
 
  m_extraInfo = m_threshold->getExtraInfo();
 
  // Book monitoring histograms
  bookHistMult(m_histAccept, "cTauMultiplicity_accept_EtaPt_"+m_threshold->name(), "Mult_"+m_threshold->name(), "#eta#times40", "E_{t} [GeV]", 200, -200, 200, 100, 0, 100);
  bookHistMult(m_histAccept, "cTauMultiplicity_accept_counts_"+m_threshold->name(), "Mult_"+m_threshold->name(), "counts", 15, 0, 15);
 
  // cTau TOB monitoring histograms
  bookHistMult(m_histcTauEt, "cTauTOBEt", "Matched cTau TOB Et", "E_{t} [GeV]", 200, 0, 400);
  bookHistMult(m_histcTauPhiEta, "cTauTOBPhiEta", "Matched cTau TOB location", "#eta#times40", "#phi#times20", 200, -200, 200, 128, 0, 128);
  bookHistMult(m_histcTauEtEta, "cTauTOBEtEta", "Matched cTau TOB Et vs eta", "#eta#times40", "E_{t} [GeV]", 200, -200, 200, 200, 0, 400);
  
  const std::string wp_name = m_threshold->isolation() != WP::NONE ? TrigConf::Selection::wpToString(m_threshold->isolation()) : "NoIso";
  bookHistMult(m_histcTauIsoFraction, "cTauTOBIsoFraction", std::string("Matched cTAU ") + wp_name + " isolation fraction", wp_name + " isolation fraction", 200, 0, 10);
 
  bookHistMult(m_histcTauIsoMatchedPass, "cTauTOBIsoMatchedPass", "Matched cTau isolation pass", "isolation pass", 2, 0, 2);
 
  return StatusCode::SUCCESS;
}

isCountingAlg()

bool TCS::ConfigurableAlg::isCountingAlg ( ) const

inlineinherited

Definition at line 64 of file ConfigurableAlg.h.

{ return m_algType == COUNT; }

isDecisionAlg()

bool TCS::ConfigurableAlg::isDecisionAlg ( ) const

inlineinherited

Definition at line 62 of file ConfigurableAlg.h.

{ return m_algType == DECISION; }

isLegacyTopo()

bool TCS::ConfigurableAlg::isLegacyTopo ( ) const

inlineinherited

Definition at line 66 of file ConfigurableAlg.h.

{ return m_isLegacyTopo; }

isocut() [1/2]

bool ConfigurableAlg::isocut ( const std::string & threshold, const unsigned int bit ) const

protectedinherited

Definition at line 484 of file ConfigurableAlg.cxx.

                                                                                {
  unsigned int value = 0;
  if (threshold == "None") {value = 0;}
  else if (threshold == "Loose") {value = 1;}
  else if (threshold == "Medium") {value = 2;}
  else if (threshold == "HadMedium") {value = 2;}
  else if (threshold == "Tight") {value = 3;}
  else {
    TRG_MSG_WARNING("No isolation defined as " << threshold);
  }
  
  if (bit >= value) {return true;}
  else {return false;}
}

isocut() [2/2]

bool ConfigurableAlg::isocut ( const unsigned int threshold, const unsigned int bit ) const

protectedinherited

Definition at line 500 of file ConfigurableAlg.cxx.

                                                                                  {
  if (bit >= threshold) {return true;}
  else {return false;}
}

isSortingAlg()

bool TCS::ConfigurableAlg::isSortingAlg ( ) const

inlineinherited

Definition at line 60 of file ConfigurableAlg.h.

{ return m_algType == SORT; }

msg() [1/2]

MsgStreamTC & TrigConf::TrigConfMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 86 of file TrigConfMessaging.h.

  {
    MsgStreamTC* ms = m_msg_tls.get();
    if (!ms) {
      ms = new MsgStreamTC(m_name);
      m_msg_tls.reset(ms);
    }
    return *ms;
  }

msg() [2/2]

MsgStreamTC & TrigConf::TrigConfMessaging::msg ( const MSGTC::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 96 of file TrigConfMessaging.h.

  {
    return msg() << lvl;
  }

msgLvl()

bool TrigConf::TrigConfMessaging::msgLvl ( const MSGTC::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicting if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 75 of file TrigConfMessaging.h.

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

name()

const std::string & TCS::ConfigurableAlg::name ( ) const

inlineinherited

Definition at line 48 of file ConfigurableAlg.h.

{ return m_name; }

numberOutputBits()

unsigned int TCS::CountingAlg::numberOutputBits ( ) const

inlineinherited

Definition at line 40 of file CountingAlg.h.

{ return m_numberOutputBits; }

parameter() [1/2]

const TCS::Parameter & ConfigurableAlg::parameter ( const std::string & parameterName ) const

inherited

Definition at line 243 of file ConfigurableAlg.cxx.

                                                                {
   return parameters().parameter(parameterName);
}

parameter() [2/2]

const TCS::Parameter & ConfigurableAlg::parameter ( const std::string & parameterName, unsigned int selection ) const

inherited

Definition at line 249 of file ConfigurableAlg.cxx.

                                                                                        {
   return parameters().parameter(parameterName, selection);
}

parameters()

const ParameterSpace & TCS::ConfigurableAlg::parameters ( ) const

inlineinherited

Definition at line 87 of file ConfigurableAlg.h.

{ return m_parameters; }

process()

TCS::StatusCode TCS::cTauMultiplicity::process ( const TCS::InputTOBArray & input, Count & count )

finaloverridevirtual

Implements TCS::CountingAlg.

Definition at line 62 of file cTauMultiplicity.cxx.

{
  const cTauTOBArray& cTaus = dynamic_cast<const cTauTOBArray&>(input);
 
  int counting = 0;
  // Loop over eTau candidates
  for(cTauTOBArray::const_iterator etauCand = cTaus.begin(); etauCand != cTaus.end(); ++etauCand) {
    if((*etauCand)->tobType() != TCS::ETAU) continue;
    
    bool accept = false; // accept = (isMatched==true && isIsolated==true) || (isMatched==false)
    bool isMatched  = false; // Is the eTau matched to a jTau?
    bool isIsolated = false; // If matched: does the resulting cTau pass the isolation cut?
    float isolation_fraction = 0; // cTAU isolation fraction
 
    // Loop over jTau candidates
    for(cTauTOBArray::const_iterator jtauCand = cTaus.begin(); jtauCand != cTaus.end(); ++jtauCand) {
      if((*jtauCand)->tobType() != TCS::JTAU) continue;
 
      isMatched = cTauMatching(*etauCand, *jtauCand);
 
      if(isMatched) {
        float etauCand_et = static_cast<float>((*etauCand)->Et());
        float etauCand_eta = static_cast<float>((*etauCand)->etaDouble());
        float jtauCand_et = static_cast<float>((*jtauCand)->Et());
        float jtauCand_etIso = static_cast<float>((*jtauCand)->EtIso());
 
        // Updated isolation condition, WP-dependent (ATR-28641)
        // I = (E_T^{jTAU Iso} + jTAUCoreScale * (E_T^{jTAU Core} - E_T^{eTAU})) / E_T^{eTAU}
        if(m_threshold->isolation() != WP::NONE) {
          isolation_fraction = (jtauCand_etIso + m_extraInfo->isolation(m_threshold->isolation(), etauCand_eta).isolation_jTAUCoreScale_d() * (jtauCand_et - etauCand_et)) / etauCand_et;
        }
        // Old isolation condition coded as in firmware: https://indico.cern.ch/event/1079697/contributions/4541419/attachments/2315137/3940824/cTAU_FirmwareAlgoProposal.pdf page 8
 
        // Check the isolation WP
        isIsolated = checkIsolationWP(*etauCand, *jtauCand);
 
        break; // Break loop when a match is found
      }
 
    } // End of jTau loop
 
    // Fill cTau TOB histograms before threshold cuts (matched cTaus only)
    if(isMatched) {
      fillHist1D(m_histcTauEt[0], (*etauCand)->EtDouble());
      fillHist2D(m_histcTauPhiEta[0], (*etauCand)->eta(), (*etauCand)->phi());
      fillHist2D(m_histcTauEtEta[0], (*etauCand)->eta(), (*etauCand)->EtDouble());
      fillHist1D(m_histcTauIsoFraction[0], isolation_fraction);
      fillHist1D(m_histcTauIsoMatchedPass[0], isMatched && isIsolated);
    }
 
    if(isMatched && isIsolated) accept = true; // This is a good matched cTau
    if(!isMatched) accept = true; // This is a non-matched eTau
 
    if(!checkeTAUWP(*etauCand)) accept = false; // Check eTAU rCore/BDT and rHad WP
 
    // Menu threshold uses 0.1 eta granularity but eFex objects have 0.025 eta granularity
    // eFex eta is calculated as 4*eta_tower (0.1 gran.) + seed (0.025 gran.), eta from -25 to 24
    int eta_thr;
    if((*etauCand)->eta()%4 >= 0) eta_thr = (*etauCand)->eta() - (*etauCand)->eta()%4;
    else                          eta_thr = (*etauCand)->eta() - (*etauCand)->eta()%4 - 4;
 
    accept = accept && (*etauCand)->Et() > m_threshold->thrValue100MeV(eta_thr/4); // Convert eta_thr to units of 0.1 to pass as an argument
 
    if(accept) {
      counting++;
      fillHist2D(m_histAccept[0], (*etauCand)->eta(), (*etauCand)->EtDouble());
    }
    
  } // End of eTau loop
 
  fillHist1D(m_histAccept[1], counting);
  
  // Pass counting to TCS::Count object - output bits are composed there
  count.setSizeCount(counting);
 
  return TCS::StatusCode::SUCCESS;
}

processBitCorrect()

TCS::StatusCode TCS::cTauMultiplicity::processBitCorrect ( const TCS::InputTOBArray & input, Count & count )

finaloverridevirtual

Implements TCS::CountingAlg.

Definition at line 56 of file cTauMultiplicity.cxx.

{
  return process(input, count);
}

quadraticSumBW()

unsigned long ConfigurableAlg::quadraticSumBW ( int i1, int i2 )

inherited

Definition at line 160 of file ConfigurableAlg.cxx.

                                              {
  return TSU::Kinematics::quadraticSumBW(i1, i2);
}

registerHist() [1/2]

void ConfigurableAlg::registerHist ( TH1 * h )

protectedinherited

Definition at line 262 of file ConfigurableAlg.cxx.

                                          {
   m_impl->registerHist(h);
}

registerHist() [2/2]

void ConfigurableAlg::registerHist ( TH2 * h )

protectedinherited

Definition at line 266 of file ConfigurableAlg.cxx.

                                          {
   m_impl->registerHist(h);
}

reset()

StatusCode TCS::ConfigurableAlg::reset ( )

inlineinherited

Definition at line 58 of file ConfigurableAlg.h.

{ return doReset(); }

setAlgoId()

void TCS::ConfigurableAlg::setAlgoId ( unsigned int algoId )

inlineinherited

Definition at line 93 of file ConfigurableAlg.h.

{ m_algoId = algoId; }

setClassName()

void TCS::ConfigurableAlg::setClassName ( const std::string & className )

inlineinherited

Definition at line 43 of file ConfigurableAlg.h.

{ m_className = className; }

setL1TopoHistSvc()

void ConfigurableAlg::setL1TopoHistSvc ( std::shared_ptr< IL1TopoHistSvc > histSvc )

inherited

sets the external hist service

forwarding public interface

Definition at line 258 of file ConfigurableAlg.cxx.

                                                                            {
   m_impl->setL1TopoHistSvc(std::move(histSvc));
}

setLegacyMode()

void TCS::ConfigurableAlg::setLegacyMode ( bool isLegacyTopo )

inlineinherited

Definition at line 45 of file ConfigurableAlg.h.

{m_isLegacyTopo=isLegacyTopo;}

setNumberOutputBits()

void TCS::CountingAlg::setNumberOutputBits ( unsigned int numberOutputBits )

inlineinherited

Definition at line 41 of file CountingAlg.h.

{ m_numberOutputBits = numberOutputBits; }

setParameters()

void ConfigurableAlg::setParameters ( const ParameterSpace & parameters )

inherited

Definition at line 222 of file ConfigurableAlg.cxx.

                                                                {
   
   // check if parameterspace is for this algorithm
   if(name() != parameters.algName()) {
      TCS_EXCEPTION("Name of parameter set (" << parameters.algName() << ") do not match configurable name (" << name() << ")" );
   }
 
   // copy parameters
   for(const Parameter & p : parameters) {
      try {
         m_parameters.setParameter(p);
      }
      catch(const TCS::Exception& e) {
         TCS_EXCEPTION("Algorithm configuration failed for " << name() << "! Caught exception: " << e.what());
      }
   }
}

setThreshold()

void TCS::CountingAlg::setThreshold ( const TrigConf::L1Threshold & thr )

inlineinherited

Definition at line 44 of file CountingAlg.h.

{ m_threshold = &thr; }

ToString()

std::string ConfigurableAlg::ToString ( const int val )

privateinherited

Definition at line 475 of file ConfigurableAlg.cxx.

{
  const int val_int = static_cast<int>(val);
  std::ostringstream temp;
  temp << val_int;
  return temp.str();
}

Member Data Documentation

m_algoId

unsigned int TCS::ConfigurableAlg::m_algoId {0}

privateinherited

Definition at line 138 of file ConfigurableAlg.h.

{0};

m_algType

AlgType TCS::ConfigurableAlg::m_algType

privateinherited

Definition at line 142 of file ConfigurableAlg.h.

m_className

std::string TCS::ConfigurableAlg::m_className {""}

privateinherited

Definition at line 135 of file ConfigurableAlg.h.

{""};

m_extraInfo

std::shared_ptr<TrigConf::L1ThrExtraInfo_cTAU> TCS::cTauMultiplicity::m_extraInfo

private

Definition at line 53 of file cTauMultiplicity.h.

m_histAccept

std::vector<std::string> TCS::CountingAlg::m_histAccept

protectedinherited

Definition at line 51 of file CountingAlg.h.

m_histcTauEt

std::vector<std::string> TCS::cTauMultiplicity::m_histcTauEt

private

Definition at line 62 of file cTauMultiplicity.h.

m_histcTauEtEta

std::vector<std::string> TCS::cTauMultiplicity::m_histcTauEtEta

private

Definition at line 64 of file cTauMultiplicity.h.

m_histcTauIsoFraction

std::vector<std::string> TCS::cTauMultiplicity::m_histcTauIsoFraction

private

Definition at line 65 of file cTauMultiplicity.h.

m_histcTauIsoMatchedPass

std::vector<std::string> TCS::cTauMultiplicity::m_histcTauIsoMatchedPass

private

Definition at line 66 of file cTauMultiplicity.h.

m_histcTauPhiEta

std::vector<std::string> TCS::cTauMultiplicity::m_histcTauPhiEta

private

Definition at line 63 of file cTauMultiplicity.h.

m_impl

std::unique_ptr<ConfigurableAlgImpl> TCS::ConfigurableAlg::m_impl

privateinherited

Definition at line 129 of file ConfigurableAlg.h.

m_isLegacyTopo

bool TCS::ConfigurableAlg::m_isLegacyTopo

privateinherited

Definition at line 144 of file ConfigurableAlg.h.

m_msg_tls

boost::thread_specific_ptr<MsgStreamTC> TrigConf::TrigConfMessaging::m_msg_tls

mutableprivateinherited

MsgStreamTC instance (a std::cout like with print-out levels)

Definition at line 71 of file TrigConfMessaging.h.

m_name

std::string TCS::ConfigurableAlg::m_name {""}

privateinherited

Definition at line 134 of file ConfigurableAlg.h.

{""};

m_numberOutputBits

unsigned int TCS::CountingAlg::m_numberOutputBits {1}

privateinherited

Definition at line 60 of file CountingAlg.h.

{1};

m_parameters

ParameterSpace TCS::ConfigurableAlg::m_parameters {""}

privateinherited

Definition at line 140 of file ConfigurableAlg.h.

{""};

m_threshold

const TrigConf::L1Threshold_cTAU* TCS::cTauMultiplicity::m_threshold {nullptr}

private

Definition at line 52 of file cTauMultiplicity.h.

{nullptr};

---

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

• cTauMultiplicity.h
• cTauMultiplicity.cxx