TrigEgammaFastElectronReAlgo Class Reference

TrigEgammaFastElectronReAlgo is a Trigger Fex Algorithm that retrieves the TrigEMCluster container and the TrackCollection containers and then creates a TrigElectron Container with a subset of calorimeter-ID selection variables that are calculated (eg. More...

#include <TrigEgammaFastElectronReAlgo.h>

Inheritance diagram for TrigEgammaFastElectronReAlgo:

Collaboration diagram for TrigEgammaFastElectronReAlgo:

Public Member Functions
	TrigEgammaFastElectronReAlgo (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) const override
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
bool	extrapolate (const EventContext &ctx, const CaloDetDescrManager &caloDD, const xAOD::TrigEMCluster *, const xAOD::TrackParticle *, double &, double &) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Static Private Member Functions
static double	getCaloPt (const xAOD::TrigElectron *el)
static double	getTkPt (const xAOD::TrigElectron *el)

Private Attributes
Gaudi::Property< bool >	m_acceptAll {this, "AcceptAll", false, "Build electrons for all tracks"}
Gaudi::Property< float >	m_clusEtthr {this, "ClusEt", 20.0*Gaudi::Units::GeV , " lower limit on cluster Et"}
Gaudi::Property< float >	m_trackPtthr {this, "TrackPt", 5.0*Gaudi::Units::GeV , "lower limit on TrackPt cut" }
Gaudi::Property< float >	m_calotrkdeta_noextrap {this, "CaloTrackdEtaNoExtrap", 0.5, "Upper limit on DEta between Calo cluster and Track for track preselection before extrapolation"}
Gaudi::Property< float >	m_trackPtthr_highet {this, "TrackPtHighEt", 2.0*Gaudi::Units::GeV , "lower limit on TrackPt cut High Et Cluster (20GeV)"}
Gaudi::Property< float >	m_calotrkdeta_noextrap_highet {this, "CaloTrackdEtaNoExtrapHighEt", 0, "upper limit on DEta between Calo cluster and Track for track preselection before extrapolation for High Et cluster (20GeV)"}
Gaudi::Property< float >	m_calotrackdeta {this, "CaloTrackdETA", 0, "Upper limit on DEta between Calo cluster and Track"}
Gaudi::Property< float >	m_calotrackdphi {this, "CaloTrackdPHI", 0, "Upper limit on DPhi between Calo cluster and Track"}
Gaudi::Property< float >	m_calotrackdeoverp_low {this, "CaloTrackdEoverPLow", 0, "lower limit on E(calo)/p(track)"}
Gaudi::Property< float >	m_calotrackdeoverp_high {this, "CaloTrackdEoverPHigh", 0, "upper limit on track E(calo)/p(track)"}
Gaudi::Property< float >	m_RCAL {this, "RCalBarrelFace", 1470.0*Gaudi::Units::mm , "Radius of inner face of the barrel calorimeter"}
Gaudi::Property< float >	m_ZCAL {this, "ZCalEndcapFace", 3800.0*Gaudi::Units::mm, "z of the inner face of endcap calorimeter"}
ToolHandle< Trk::IParticleCaloExtensionTool >	m_caloExtensionTool
SG::ReadHandleKey< TrigRoiDescriptorCollection >	m_roiCollectionKey
SG::ReadHandleKey< xAOD::TrigEMClusterContainer >	m_TrigEMClusterContainerKey
SG::ReadHandleKey< xAOD::TrackParticleContainer >	m_TrackParticleContainerKey
SG::ReadCondHandleKey< CaloDetDescrManager >	m_caloDetDescrMgrKey
SG::WriteHandleKey< xAOD::TrigElectronContainer >	m_outputElectronsKey
SG::WriteHandleKey< xAOD::TrigElectronContainer >	m_outputDummyElectronsKey
ToolHandle< GenericMonitoringTool >	m_monTool { this, "MonTool", "", "Monitoring tool" }
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default).
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

TrigEgammaFastElectronReAlgo is a Trigger Fex Algorithm that retrieves the TrigEMCluster container and the TrackCollection containers and then creates a TrigElectron Container with a subset of calorimeter-ID selection variables that are calculated (eg.

E/p) The algo will apply some very loose selection cuts to the TrigElectronContainer created which is of TrigParticle type TrigElectron The TrigElectron conatiner will then be retrieved by the hypothesis algorithm TrigEgammaFastElectronHypo that will perform the corresponding fast electron selection

Definition at line 49 of file TrigEgammaFastElectronReAlgo.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

TrigEgammaFastElectronReAlgo()

TrigEgammaFastElectronReAlgo::TrigEgammaFastElectronReAlgo	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 30 of file TrigEgammaFastElectronReAlgo.cxx.

    : AthReentrantAlgorithm(name, pSvcLocator)
{}

Member Function Documentation

cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode TrigEgammaFastElectronReAlgo::execute ( const EventContext & ctx ) const

overridevirtual

< monitor preselection between track eta and cluster before extrapolation

Create a TrigElectron corresponding to this candidate assume cluster quantities give better estimate of transverse energy (probably a safe assumption for large pT) and that track parameters at perigee give better estimates of angular quantities

Definition at line 76 of file TrigEgammaFastElectronReAlgo.cxx.

{
  using namespace xAOD;   
 
  ATH_MSG_DEBUG( "Executing " <<name());
 
  auto trigElecColl =   SG::makeHandle (m_outputElectronsKey, ctx);  
  ATH_CHECK( trigElecColl.record (std::make_unique<xAOD::TrigElectronContainer>(),
                           std::make_unique<xAOD::TrigEMClusterAuxContainer>()) );
  
  ATH_MSG_DEBUG( "Made WriteHandle " << m_outputElectronsKey );
 
  auto trigDummyElecColl =   SG::makeHandle (m_outputDummyElectronsKey, ctx);
 
  ATH_MSG_DEBUG( "Made Dummy WriteHandle " << m_outputDummyElectronsKey );
  ATH_CHECK( trigDummyElecColl.record (std::make_unique<xAOD::TrigElectronContainer>(),
                           std::make_unique<xAOD::TrigEMClusterAuxContainer>()) );  
 
  auto roiCollection = SG::makeHandle(m_roiCollectionKey, ctx);
  ATH_MSG_DEBUG( "Made handle " << m_roiCollectionKey  );
  
 
  if (roiCollection->size()==0) {
    ATH_MSG_DEBUG(" RoI collection size = 0");
    return StatusCode::SUCCESS;
  }
 
  const TrigRoiDescriptor* roiDescriptor = *(roiCollection->begin());
 
  ATH_MSG_DEBUG(" RoI ID = "   << (roiDescriptor)->roiId() << ": Eta = "      << (roiDescriptor)->eta()
                << ", Phi = "      << (roiDescriptor)->phi());
  
  float calo_eta(999), calo_phi(999), calo_et(-1);
 
  auto clusterCont = SG::makeHandle (m_TrigEMClusterContainerKey, ctx);
  ATH_MSG_DEBUG( "Made handle " << m_TrigEMClusterContainerKey  );
  
  //JTB Should only be 1 cluster in each RoI 
 
  const xAOD::TrigEMCluster *emCluster=(*clusterCont->begin());
  ElementLink<xAOD::TrigEMClusterContainer> clusEL=ElementLink<xAOD::TrigEMClusterContainer> (*clusterCont,0);
 
  // copy relevant quantities to local variables
  calo_eta = emCluster->eta();
  calo_phi = emCluster->phi();
  calo_et  = emCluster->et();
 
 
  // Transverse em energy
  ATH_MSG_DEBUG("Cluster: ET=" << calo_et);
  ATH_MSG_DEBUG("searching a matching track: loop over tracks");
 
  SG::ReadHandle<xAOD::TrackParticleContainer> tracks(m_TrackParticleContainerKey, ctx);
  ATH_MSG_DEBUG( "Made handle " << m_TrackParticleContainerKey  );
 
 
  if (tracks->size() == 0){
      return StatusCode::SUCCESS; // Exit early if there are no tracks
  }
 
  size_t coll_size = tracks->size();
  trigElecColl->reserve(coll_size);
 
  // monitoring
  std::vector<float> calotrkdeta_noextrap_mon; 
 
  auto caloPtMon = Monitored::Collection("PtCalo",   *trigElecColl,  getCaloPt );
  auto trackPtMon = Monitored::Collection("PtTrack", *trigElecColl, getTkPt );
  auto caloTrackDEtaMon = Monitored::Collection("CaloTrackdEta", *trigElecColl, &xAOD::TrigElectron::trkClusDeta );
  auto caloTrackDPhiMon = Monitored::Collection("CaloTrackdPhi", *trigElecColl, &xAOD::TrigElectron::trkClusDphi );
  auto etOverPtMon = Monitored::Collection("CaloTrackEoverP", *trigElecColl, &xAOD::TrigElectron::etOverPt );
  auto caloTrackDEtaNoExtrapMon = Monitored::Collection("CaloTrackdEtaNoExtrapMon", calotrkdeta_noextrap_mon );
 
  auto mon = Monitored::Group(m_monTool,  caloPtMon, trackPtMon, caloTrackDEtaMon, caloTrackDPhiMon, etOverPtMon, caloTrackDEtaNoExtrapMon );
 
  // Make Dummy Electron
  xAOD::TrigElectron* trigDummyElec = new xAOD::TrigElectron();
 
  ElementLink<xAOD::TrackParticleContainer> trackDummyEL = ElementLink<xAOD::TrackParticleContainer> (*tracks, 0);
     
  trigDummyElecColl->push_back(trigDummyElec);
  trigDummyElec->init( 0, 0, 0,  0, clusEL, trackDummyEL);
  
  //Pick the Calo Det Descr 
  SG::ReadCondHandle<CaloDetDescrManager> caloDetDescrMgrHandle{
      m_caloDetDescrMgrKey, ctx};
  ATH_CHECK(caloDetDescrMgrHandle.isValid());
  const CaloDetDescrManager* caloDD = *caloDetDescrMgrHandle;
 
  // loop over tracks
  for(unsigned int track_index = 0; track_index < tracks->size(); track_index++)
  {
    const xAOD::TrackParticle_v1* trkIter = (*tracks)[track_index];
    ATH_MSG_DEBUG("Track loop starts");
    ATH_MSG_VERBOSE("AlgoId = " << (trkIter)->patternRecoInfo());
    ATH_MSG_VERBOSE("At perigee:");
    ATH_MSG_DEBUG(" Pt  = " << std::abs((trkIter)->pt())); 
    ATH_MSG_DEBUG(" phi = " << std::abs((trkIter)->phi0()));
    ATH_MSG_DEBUG(" eta = " << std::abs((trkIter)->eta())); 
    ATH_MSG_DEBUG(" z0  = " << std::abs((trkIter)->z0()));
    ATH_MSG_DEBUG(" d0  = " << std::abs((trkIter)->d0()));  
    // ============================================= //
    // Pt cut
    float trkPt = std::abs((trkIter)->pt());
    float etoverpt = std::abs(calo_et/trkPt);
    float calotrkdeta_noextrap = (trkIter)->eta() - calo_eta;
    
    double etaAtCalo=999.;
    double phiAtCalo=999.;
    if(m_acceptAll){
      if(!extrapolate(ctx,*caloDD,emCluster,trkIter,etaAtCalo,phiAtCalo)){
        ATH_MSG_VERBOSE("extrapolator failed");
        continue; 
      }
      else{
        ATH_MSG_VERBOSE("REGTEST: TrigElectron: cluster index = " << clusEL.index() <<
               " eta = " << etaAtCalo << " phi = " << phiAtCalo);
 
        xAOD::TrigElectron* trigElec = new xAOD::TrigElectron();
        trigElecColl->push_back(trigElec);
        ElementLink<xAOD::TrackParticleContainer> trackEL = ElementLink<xAOD::TrackParticleContainer> (*tracks, track_index);
        trigElec->init(  roiDescriptor->roiWord(),
                    etaAtCalo, phiAtCalo,  etoverpt,        
                    clusEL,
                    trackEL);
        calotrkdeta_noextrap_mon.push_back(calotrkdeta_noextrap);
      }
    }else {
 
      ATH_MSG_DEBUG("Apply cuts on track with index: "<<track_index);
      if(trkPt < m_trackPtthr){
        ATH_MSG_DEBUG("Failed track pt cut " << trkPt);
        continue;
      }
      if(std::abs(calotrkdeta_noextrap) > m_calotrkdeta_noextrap){
        ATH_MSG_DEBUG("Failed pre extrapolation calo track deta " << calotrkdeta_noextrap);
        continue;
      }
      if(calo_et > m_clusEtthr){
        if(trkPt < m_trackPtthr_highet){
          ATH_MSG_DEBUG("Failed track pt cut for high et cluster");
          continue;
        }
        if(calotrkdeta_noextrap > m_calotrkdeta_noextrap_highet){
          ATH_MSG_DEBUG("Failed pre extrapolation calo track deta for high et");
          continue;
        }
      }
      if (etoverpt < m_calotrackdeoverp_low){ 
        ATH_MSG_DEBUG("failed low cut on ET/PT");
        continue;
      }
      if (etoverpt > m_calotrackdeoverp_high){ 
        ATH_MSG_DEBUG("failed high cut on ET/PT");
        continue;
      }
      if(!extrapolate(ctx,*caloDD,emCluster,trkIter,etaAtCalo,phiAtCalo)){
        ATH_MSG_DEBUG("extrapolator failed 1");
        continue; 
      }
      // all ok: do track-matching cuts
      ATH_MSG_DEBUG("extrapolated eta/phi=" << etaAtCalo << "/" << phiAtCalo);
      // match in eta
      float dEtaCalo = std::abs(etaAtCalo - calo_eta);
      ATH_MSG_DEBUG("deta = " << dEtaCalo); 
      if ( dEtaCalo > m_calotrackdeta){ 
        ATH_MSG_DEBUG("failed eta match cut " << dEtaCalo);
        continue;
      }
      // match in phi: deal with differences larger than Pi
      float dPhiCalo =  std::abs(phiAtCalo - calo_phi);
      dPhiCalo = ( dPhiCalo < M_PI ? dPhiCalo : 2*M_PI - dPhiCalo );
      ATH_MSG_DEBUG("dphi = " << dPhiCalo); 
      if ( dPhiCalo > m_calotrackdphi) {
        ATH_MSG_DEBUG("failed phi match cut " << dPhiCalo);
        continue;
      }
      // all cuts passed
      ElementLink<xAOD::TrackParticleContainer> trackEL = ElementLink<xAOD::TrackParticleContainer> (*tracks, track_index);
      ATH_MSG_DEBUG("pT of trackEL: "<<(*trackEL)->pt());     
 
      ATH_MSG_DEBUG("REGTEST: TrigElectron: cluster index = " << clusEL.index() <<
                " track index= "     << track_index << " eta = " <<  etaAtCalo << " phi = " << phiAtCalo << 
                " deta = " << dEtaCalo << "dphi = " << dPhiCalo);
      xAOD::TrigElectron* trigElec = new xAOD::TrigElectron();
      trigElecColl->push_back(trigElec);
      trigElec->init( roiDescriptor->roiWord(),
                etaAtCalo, phiAtCalo,  etoverpt,        
                clusEL,
                trackEL);
      
      ATH_MSG_DEBUG(" deta = " << dEtaCalo << " deta = " << trigElec->trkClusDeta() 
                << " dphi = " << dPhiCalo << " dphi = " << trigElec->trkClusDphi()
                << " caloEta = " << calo_eta << " caloEta = " << trigElec->caloEta()
                << " caloPhi = " << calo_phi << " calophi = " << trigElec->caloPhi()
                << " etaAtCalo = " << etaAtCalo << " etaAtCalo = " << trigElec->trkEtaAtCalo()
                << " phiAtCalo = " << phiAtCalo << " phiAtCalo = " << trigElec->trkPhiAtCalo()
                );
      calotrkdeta_noextrap_mon.push_back(calotrkdeta_noextrap);
     
    }
    
  }
 
 
  ATH_MSG_DEBUG("REGTEST:  returning an xAOD::TrigElectronContainer with size "<< trigElecColl->size() << ".");
  
  return StatusCode::SUCCESS;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

extrapolate()

bool TrigEgammaFastElectronReAlgo::extrapolate ( const EventContext & ctx, const CaloDetDescrManager & caloDD, const xAOD::TrigEMCluster * clus, const xAOD::TrackParticle * trk, double & etaAtCalo, double & phiAtCalo ) const

private

Definition at line 290 of file TrigEgammaFastElectronReAlgo.cxx.

                                                {
 
  // use the provided EM cluster to "guide" the extrapolation
  // 1st figure which layer we  want to shoot at
  // in this case we chose between EMB2 or EME2
  // given the cluster.
  std::vector<CaloSampling::CaloSample> samples;
  if (clus->energy(CaloSampling::CaloSample::EME2) >
      clus->energy(CaloSampling::CaloSample::EMB2)) {
    samples.push_back(CaloSampling::CaloSample::EME2);
  } else {
    samples.push_back(CaloSampling::CaloSample::EMB2);
  }
 
  // create the surfaces we want to reach.
  // Aka either a cylinder for EME2 or a disc for EMB2
  std::vector<std::unique_ptr<Trk::Surface>> caloSurfaces =
      m_caloExtensionTool->caloSurfacesFromLayers(samples, clus->eta(), caloDD);
 
  // And then try to reach them
  const auto extension = m_caloExtensionTool->surfaceCaloExtension(
      ctx, trk->perigeeParameters(), samples, caloSurfaces,
      Trk::nonInteracting);
 
  if (extension.empty()) {
    ATH_MSG_VERBOSE("extrapolator failed 1");
    return false;
  }
  // We target exactly one EMB2 or EME2  (the vector has 1 entry if not empty)
  etaAtCalo = extension[0].second->position().eta();
  phiAtCalo = extension[0].second->position().phi();
 
  ATH_MSG_VERBOSE("Hit sampling :" << extension.at(0).first << " at eta : "
                                   << etaAtCalo << " at phi : " << phiAtCalo);
 
  return true;
}

filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

getCaloPt()

double TrigEgammaFastElectronReAlgo::getCaloPt ( const xAOD::TrigElectron * el )

inlinestaticprivate

Definition at line 63 of file TrigEgammaFastElectronReAlgo.h.

                                                               {
      if(!el) return -999.0;
      return el->pt();
    }

getTkPt()

double TrigEgammaFastElectronReAlgo::getTkPt ( const xAOD::TrigElectron * el )

inlinestaticprivate

Definition at line 68 of file TrigEgammaFastElectronReAlgo.h.

                                                             {
      if(!el) return -999.0;
      return (el->trackParticle() ? el->trackParticle()->pt()   : -999);
    }

initialize()

StatusCode TrigEgammaFastElectronReAlgo::initialize ( )

overridevirtual

Definition at line 35 of file TrigEgammaFastElectronReAlgo.cxx.

{
  ATH_MSG_DEBUG("Initialization:");
 
  if (!m_monTool.empty()) {
    ATH_MSG_DEBUG("Retrieving monTool");
    CHECK(m_monTool.retrieve());
  } else {
    ATH_MSG_DEBUG("No monTool configured => NO MONITOING");
  }
 
  if ( m_caloExtensionTool.retrieve().isFailure() ) {  
        ATH_MSG_ERROR("Unable to locate TrackExtrapolator tool ");
        return StatusCode::FAILURE; 
  }  
 
  // print out settings
  ATH_MSG_DEBUG("Initialization completed successfully:"); 
  ATH_MSG_DEBUG( "AcceptAll            = "<< (m_acceptAll==true ? "True" : "False")); 
  ATH_MSG_DEBUG("TrackPt              = " << m_trackPtthr);          
  ATH_MSG_DEBUG("TrackPtHighEt        = " << m_trackPtthr_highet);          
  ATH_MSG_DEBUG("ClusEt               = " << m_clusEtthr);          
  ATH_MSG_DEBUG("CaloTrackdEtaNoExtrap= " << m_calotrkdeta_noextrap);        
  ATH_MSG_DEBUG("CaloTrackdEtaNoExtrapHighEt= " << m_calotrkdeta_noextrap_highet);        
  ATH_MSG_DEBUG("CaloTrackdETA        = " << m_calotrackdeta);        
  ATH_MSG_DEBUG("CaloTrackdPHI        = " << m_calotrackdphi);        
  ATH_MSG_DEBUG("CaloTrackdEoverPLow  = " << m_calotrackdeoverp_low); 
  ATH_MSG_DEBUG("CaloTrackdEoverPHigh = " << m_calotrackdeoverp_high);
 
  ATH_CHECK( m_roiCollectionKey.initialize() );
  ATH_CHECK( m_TrigEMClusterContainerKey.initialize() );
  ATH_CHECK( m_TrackParticleContainerKey.initialize() );
  ATH_CHECK( m_caloDetDescrMgrKey.initialize() );
  ATH_CHECK( m_outputElectronsKey.initialize() );
  ATH_CHECK( m_outputDummyElectronsKey.initialize() );
 
  return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

                                                               {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_acceptAll

Gaudi::Property<bool> TrigEgammaFastElectronReAlgo::m_acceptAll {this, "AcceptAll", false, "Build electrons for all tracks"}

private

Definition at line 80 of file TrigEgammaFastElectronReAlgo.h.

{this, "AcceptAll", false, "Build electrons for all tracks"};

m_caloDetDescrMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> TrigEgammaFastElectronReAlgo::m_caloDetDescrMgrKey

private

Initial value:

{
        this, 
        "CaloDetDescrManager", 
        "CaloDetDescrManager",
        "SG Key for CaloDetDescrManager in the Condition Store"}

Definition at line 114 of file TrigEgammaFastElectronReAlgo.h.

                                                                 {
        this, 
        "CaloDetDescrManager", 
        "CaloDetDescrManager",
        "SG Key for CaloDetDescrManager in the Condition Store"};

m_caloExtensionTool

ToolHandle<Trk::IParticleCaloExtensionTool> TrigEgammaFastElectronReAlgo::m_caloExtensionTool

private

Initial value:

{
        this, "ParticleCaloExtensionTool",
        "Trk::ParticleCaloExtensionTool/ParticleCaloExtensionTool",
        "Tool to extrapolate Track to Calo inner surface"}

Definition at line 94 of file TrigEgammaFastElectronReAlgo.h.

                                                                 {
        this, "ParticleCaloExtensionTool",
        "Trk::ParticleCaloExtensionTool/ParticleCaloExtensionTool",
        "Tool to extrapolate Track to Calo inner surface"};

m_calotrackdeoverp_high

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_calotrackdeoverp_high {this, "CaloTrackdEoverPHigh", 0, "upper limit on track E(calo)/p(track)"}

private

Definition at line 89 of file TrigEgammaFastElectronReAlgo.h.

{this,  "CaloTrackdEoverPHigh", 0, "upper limit on track E(calo)/p(track)"};

m_calotrackdeoverp_low

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_calotrackdeoverp_low {this, "CaloTrackdEoverPLow", 0, "lower limit on E(calo)/p(track)"}

private

Definition at line 88 of file TrigEgammaFastElectronReAlgo.h.

{this,  "CaloTrackdEoverPLow",  0, "lower limit on E(calo)/p(track)"};

m_calotrackdeta

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_calotrackdeta {this, "CaloTrackdETA", 0, "Upper limit on DEta between Calo cluster and Track"}

private

Definition at line 86 of file TrigEgammaFastElectronReAlgo.h.

{this,  "CaloTrackdETA",    0, "Upper limit on DEta between Calo cluster and Track"};

m_calotrackdphi

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_calotrackdphi {this, "CaloTrackdPHI", 0, "Upper limit on DPhi between Calo cluster and Track"}

private

Definition at line 87 of file TrigEgammaFastElectronReAlgo.h.

{this,  "CaloTrackdPHI",     0, "Upper limit on DPhi between Calo cluster and Track"}; 

m_calotrkdeta_noextrap

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_calotrkdeta_noextrap {this, "CaloTrackdEtaNoExtrap", 0.5, "Upper limit on DEta between Calo cluster and Track for track preselection before extrapolation"}

private

Definition at line 83 of file TrigEgammaFastElectronReAlgo.h.

{this,  "CaloTrackdEtaNoExtrap",   0.5, "Upper limit on DEta between Calo cluster and Track for track preselection before extrapolation"};

m_calotrkdeta_noextrap_highet

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_calotrkdeta_noextrap_highet {this, "CaloTrackdEtaNoExtrapHighEt", 0, "upper limit on DEta between Calo cluster and Track for track preselection before extrapolation for High Et cluster (20GeV)"}

private

Definition at line 85 of file TrigEgammaFastElectronReAlgo.h.

{this,  "CaloTrackdEtaNoExtrapHighEt",  0, "upper limit on DEta between Calo cluster and Track for track preselection before extrapolation for High Et cluster (20GeV)"};

m_clusEtthr

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_clusEtthr {this, "ClusEt", 20.0*Gaudi::Units::GeV , " lower limit on cluster Et"}

private

Definition at line 81 of file TrigEgammaFastElectronReAlgo.h.

{this,  "ClusEt",  20.0*Gaudi::Units::GeV , " lower limit on cluster Et"};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_monTool

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

private

Definition at line 131 of file TrigEgammaFastElectronReAlgo.h.

{ this, "MonTool", "", "Monitoring tool" };

m_outputDummyElectronsKey

SG::WriteHandleKey<xAOD::TrigElectronContainer> TrigEgammaFastElectronReAlgo::m_outputDummyElectronsKey

private

Initial value:

{ this,
        "DummyElectronsName",                  
        "Electrons",                           
        "output Dummy Electron container name "}

Definition at line 126 of file TrigEgammaFastElectronReAlgo.h.

                                                                         { this,
        "DummyElectronsName",                  // property name
        "Electrons",                           // default value of StoreGate key
        "output Dummy Electron container name "};

m_outputElectronsKey

SG::WriteHandleKey<xAOD::TrigElectronContainer> TrigEgammaFastElectronReAlgo::m_outputElectronsKey

private

Initial value:

{
        this,
        "ElectronsName",  
        "Electrons",      
        "output Electron container name "}

Definition at line 120 of file TrigEgammaFastElectronReAlgo.h.

                                                                    {
        this,
        "ElectronsName",  // property name
        "Electrons",      // default value of StoreGate key
        "output Electron container name "};

m_RCAL

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_RCAL {this, "RCalBarrelFace", 1470.0*Gaudi::Units::mm , "Radius of inner face of the barrel calorimeter"}

private

Definition at line 90 of file TrigEgammaFastElectronReAlgo.h.

{this,  "RCalBarrelFace",  1470.0*Gaudi::Units::mm , "Radius of inner face of the barrel calorimeter"};

m_roiCollectionKey

SG::ReadHandleKey<TrigRoiDescriptorCollection> TrigEgammaFastElectronReAlgo::m_roiCollectionKey

private

Initial value:

{ this, 
        "RoIs",                             
        "rois",                             
        "input RoI Collection name"}

Definition at line 99 of file TrigEgammaFastElectronReAlgo.h.

                                                                    { this, 
        "RoIs",                             // property name
        "rois",                             // default value of StoreGate key
        "input RoI Collection name"};

m_TrackParticleContainerKey

SG::ReadHandleKey<xAOD::TrackParticleContainer> TrigEgammaFastElectronReAlgo::m_TrackParticleContainerKey

private

Initial value:

{ this,
        "TrackParticlesName",         
        "Tracks",                     
        "input TrackParticle container name"}

Definition at line 109 of file TrigEgammaFastElectronReAlgo.h.

                                                                           { this,
        "TrackParticlesName",         // property name
        "Tracks",                     // default value of StoreGate key
        "input TrackParticle container name"};

m_trackPtthr

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_trackPtthr {this, "TrackPt", 5.0*Gaudi::Units::GeV , "lower limit on TrackPt cut" }

private

Definition at line 82 of file TrigEgammaFastElectronReAlgo.h.

{this,  "TrackPt",  5.0*Gaudi::Units::GeV , "lower limit on TrackPt cut" };

m_trackPtthr_highet

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_trackPtthr_highet {this, "TrackPtHighEt", 2.0*Gaudi::Units::GeV , "lower limit on TrackPt cut High Et Cluster (20GeV)"}

private

Definition at line 84 of file TrigEgammaFastElectronReAlgo.h.

{this,  "TrackPtHighEt",  2.0*Gaudi::Units::GeV , "lower limit on TrackPt cut High Et Cluster (20GeV)"};

m_TrigEMClusterContainerKey

SG::ReadHandleKey<xAOD::TrigEMClusterContainer> TrigEgammaFastElectronReAlgo::m_TrigEMClusterContainerKey

private

Initial value:

{ this,
        "TrigEMClusterName",       
        "clusters",                
        "input TrigEMCluster Container name"}

Definition at line 104 of file TrigEgammaFastElectronReAlgo.h.

                                                                           { this,
        "TrigEMClusterName",       // property name
        "clusters",                // default value of StoreGate key
        "input TrigEMCluster Container name"};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_ZCAL

Gaudi::Property<float> TrigEgammaFastElectronReAlgo::m_ZCAL {this, "ZCalEndcapFace", 3800.0*Gaudi::Units::mm, "z of the inner face of endcap calorimeter"}

private

Definition at line 91 of file TrigEgammaFastElectronReAlgo.h.

{this,  "ZCalEndcapFace",     3800.0*Gaudi::Units::mm, "z of the inner face of endcap calorimeter"};

---

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

• TrigEgammaFastElectronReAlgo.h
• TrigEgammaFastElectronReAlgo.cxx