#include <TrigEgammaForwardPrecisionCaloHypoTool.h>

Inheritance diagram for TrigEgammaForwardPrecisionCaloHypoTool:

Collaboration diagram for TrigEgammaForwardPrecisionCaloHypoTool:

Public Member Functions
	TrigEgammaForwardPrecisionCaloHypoTool (const std::string &type, const std::string &name, const IInterface *parent)

virtual StatusCode	initialize () override

virtual StatusCode	decide (std::vector< ITrigEgammaForwardPrecisionCaloHypoTool::ClusterInfo > &input) const override

virtual bool	decide (const ITrigEgammaForwardPrecisionCaloHypoTool::ClusterInfo &i) const override

Private Member Functions
int	findCutIndex (float eta) const

Private Attributes
HLT::Identifier	m_decisionId

Gaudi::Property< std::vector< float > >	m_etabin { this, "EtaBins", {} , "Bins of eta" }
	selection variable for PRECISION calo selection:eta bins More...

Gaudi::Property< std::vector< float > >	m_eTthr { this, "ETthr", {}, "ET Threshold" }

Gaudi::Property< std::vector< float > >	m_eT2thr { this, "ET2thr", {}, "Second layer ET threshold" }

Gaudi::Property< float >	m_detacluster { this, "dETACLUSTERthr", 0. , "" }

Gaudi::Property< float >	m_dphicluster { this, "dPHICLUSTERthr", 0. , "" }

ToolHandle< GenericMonitoringTool >	m_monTool { this, "MonTool", "", "Monitoring tool" }

Detailed Description

Definition at line 18 of file TrigEgammaForwardPrecisionCaloHypoTool.h.

Constructor & Destructor Documentation

◆ TrigEgammaForwardPrecisionCaloHypoTool()

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

Definition at line 16 of file TrigEgammaForwardPrecisionCaloHypoTool.cxx.

19 : base_class( type, name, parent ),

20 m_decisionId( HLT::Identifier::fromToolName( name ) ) {}

Member Function Documentation

◆ decide() [1/2]

bool TrigEgammaForwardPrecisionCaloHypoTool::decide ( const ITrigEgammaForwardPrecisionCaloHypoTool::ClusterInfo & i ) const

overridevirtual

Definition at line 55 of file TrigEgammaForwardPrecisionCaloHypoTool.cxx.

                                                                                                                            {
  
   bool pass = false;
  
   auto dEta         = Monitored::Scalar( "dEta", -1. ); 
   auto dPhi         = Monitored::Scalar( "dPhi", -1. );
   auto eT_Cluster   = Monitored::Scalar( "Et_em"   , -1.0 );
   auto etaBin       = Monitored::Scalar( "EtaBin", -1. );
   auto monEta       = Monitored::Scalar( "Eta", -99. ); 
   auto monPhi       = Monitored::Scalar( "Phi", -99. );
   auto PassedCuts   = Monitored::Scalar<int>( "CutCounter", -1 );  
   auto monitorIt    = Monitored::Group( m_monTool,dEta, dPhi, eT_Cluster,
                                         etaBin, monEta,monPhi,PassedCuts );
  // when leaving scope it will ship data to monTool
   PassedCuts = PassedCuts + 1; //got called (data in place)
  
   auto roiDescriptor = input.roi;
  
  
   if ( std::abs( roiDescriptor->eta() ) > 2.6 ) {
       ATH_MSG_DEBUG( "REJECT The cluster had eta coordinates beyond the EM fiducial volume : " << roiDescriptor->eta() << "; stop the chain now" );
       pass=false; // special case       
       return pass;
   } 
  
   ATH_MSG_INFO( "; RoI ID = " << roiDescriptor->roiId()
       << ": Eta = " << roiDescriptor->eta()
       << ", Phi = " << roiDescriptor->phi() );
  
   // fill local variables for RoI reference position
   double etaRef = roiDescriptor->eta();
   double phiRef = roiDescriptor->phi();
   // correct phi the to right range ( probably not needed anymore )   
   if ( std::abs( phiRef ) > M_PI ) phiRef -= 2*M_PI; // correct phi if outside range
  
  
   auto pClus = input.cluster;
   
   float absEta = std::abs( pClus->eta() );
   const int cutIndex = findCutIndex( absEta );
   
  
   
   dEta =  pClus->eta() - etaRef;
   //  Deal with angle diferences greater than Pi
   dPhi =  std::abs( pClus->phi() - phiRef );
   dPhi = ( dPhi < M_PI ? dPhi : 2*M_PI - dPhi ); // TB why only <
   eT_Cluster  = pClus->et();
   // apply cuts: DeltaEta( clus-ROI )
   ATH_MSG_INFO( "CaloCluster: eta="  << pClus->eta()
         << " roi eta=" << etaRef << " DeltaEta=" << dEta
         << " cut: <"   << m_detacluster          );
   
   if ( std::abs( pClus->eta() - etaRef ) > m_detacluster ) {
     ATH_MSG_DEBUG("REJECT Cluster dEta cut failed");
     return pass;
   }
   PassedCuts = PassedCuts + 1; //Deta
   
   // DeltaPhi( clus-ROI )
   ATH_MSG_DEBUG( ": phi="  << pClus->phi()
        << " roi phi="<< phiRef    << " DeltaPhi="<< dPhi
        << " cut: <"  << m_dphicluster );
   
   if( dPhi > m_dphicluster ) {
     ATH_MSG_DEBUG("REJECT Clsuter dPhi cut failed");
     return pass;
   }
   PassedCuts = PassedCuts + 1; //DPhi
  
   // eta range
   if ( cutIndex == -1 ) {  // VD
     ATH_MSG_DEBUG( "Cluster eta: " << absEta << " outside eta range " << m_etabin[m_etabin.size()-1] );
     return pass;
   } else { 
     ATH_MSG_DEBUG( "eta bin used for cuts " << cutIndex );
   }
   PassedCuts = PassedCuts + 1; // passed eta cut
   
   // ET_em
   ATH_MSG_DEBUG( "CaloCluster: ET_em=" << eT_Cluster << " cut: >"  << m_eTthr[cutIndex] );
   if ( eT_Cluster < m_eTthr[cutIndex] ) {
     ATH_MSG_DEBUG("REJECT et cut failed");
     return pass;
   }
   PassedCuts = PassedCuts + 1; // ET_em
   
   // got this far => passed!
   pass = true;
  
   // Reach this point successfully  
   ATH_MSG_DEBUG( "pass = " << pass );
  
   return pass;
  
 }

◆ decide() [2/2]

virtual StatusCode TrigEgammaForwardPrecisionCaloHypoTool::decide ( std::vector< ITrigEgammaForwardPrecisionCaloHypoTool::ClusterInfo > & input ) const

overridevirtual

◆ findCutIndex()

int TrigEgammaForwardPrecisionCaloHypoTool::findCutIndex ( float eta ) const

private

Definition at line 152 of file TrigEgammaForwardPrecisionCaloHypoTool.cxx.

                                                                           {
   const float absEta = std::abs(eta);
   
   auto binIterator = std::adjacent_find( m_etabin.begin(), m_etabin.end(), [=](float left, float right){ return left < absEta and absEta < right; }  );
   if ( binIterator == m_etabin.end() ) {
     return -1;
   }
   return  binIterator - m_etabin.begin();
 }

◆ initialize()

StatusCode TrigEgammaForwardPrecisionCaloHypoTool::initialize ( )

overridevirtual

Definition at line 22 of file TrigEgammaForwardPrecisionCaloHypoTool.cxx.

 {
   ATH_MSG_DEBUG( "Initialization completed successfully"   );    
   ATH_MSG_DEBUG( "EtaBins        = " << m_etabin      );
   ATH_MSG_DEBUG( "ETthr          = " << m_eTthr    << "( lo )/" << m_eT2thr    << "( hi )" );  
   ATH_MSG_DEBUG( "dPHICLUSTERthr = " << m_dphicluster );
   ATH_MSG_DEBUG( "dETACLUSTERthr = " << m_detacluster );
   
   if ( m_etabin.size() == 0 ) {
     ATH_MSG_ERROR(  " There are no cuts set (EtaBins property is an empty list)" );
     return StatusCode::FAILURE;
   }
  
   unsigned int nEtaBin = m_etabin.size();
   
 #define CHECK_SIZE( __n) if ( m_##__n.size() !=  (nEtaBin - 1) )        \
     { ATH_MSG_DEBUG(" __n size is " << m_##__n.size() << " but needs to be " << (nEtaBin - 1) ); return StatusCode::FAILURE; }
  
   CHECK_SIZE( eTthr );
   CHECK_SIZE( eT2thr );
 #undef CHECK_SIZE
  
  
  
   ATH_MSG_DEBUG( "Tool configured for chain/id: " << m_decisionId );
  
   if ( not m_monTool.name().empty() ) 
     CHECK( m_monTool.retrieve() );
  
   return StatusCode::SUCCESS;
 }