JetHitAssociation Class Reference

#include <JetHitAssociation.h>

Inheritance diagram for JetHitAssociation:

Collaboration diagram for JetHitAssociation:

Public Member Functions
	JetHitAssociation (const std::string &name, ISvcLocator *svcloc)
	Constructor with parameters: More...
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual StatusCode	execute () override
virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
virtual StatusCode	sysStart () override
	Handle START transition. More...
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...
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. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	saveHits (const std::vector< const xAOD::Jet * > &jets, const xAOD::Vertex *const vertex, SG::ReadHandle< xAOD::TrackMeasurementValidationContainer > &hits, const SG::WriteHandle< xAOD::TrackMeasurementValidationContainer > &writeHandle, unsigned long long int &nStoredHits, unsigned long long int &nTotalHits)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
SG::ReadHandleKey< xAOD::JetContainer >	m_jetCollectionName
SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexCollectionName
SG::ReadHandleKey< xAOD::TrackMeasurementValidationContainer >	m_inputPixHitCollectionName
SG::ReadHandleKey< xAOD::TrackMeasurementValidationContainer >	m_inputSCTHitCollectionName
SG::WriteHandleKey< xAOD::TrackMeasurementValidationContainer >	m_outputPixHitCollectionName
SG::WriteHandleKey< xAOD::TrackMeasurementValidationContainer >	m_outputSCTHitCollectionName
float	m_jetPtThreshold
float	m_dRmatchHitToJet
unsigned long long int	m_nStoredPixHits
unsigned long long int	m_nStoredSCTHits
unsigned long long int	m_nTotalPixHits
unsigned long long int	m_nTotalSCTHits
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 24 of file JetHitAssociation.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

JetHitAssociation()

JetHitAssociation::JetHitAssociation	(	const std::string &	name,
		ISvcLocator *	svcloc
	)

Constructor with parameters:

Definition at line 11 of file JetHitAssociation.cxx.

  : AthAlgorithm(name, svcloc),
    m_nStoredPixHits(0),
    m_nStoredSCTHits(0),
    m_nTotalPixHits(0),
    m_nTotalSCTHits(0) {
    
  declareProperty("jetCollectionName", m_jetCollectionName = "AntiKt4EMPFlowJets", 
              "Jet collection that will be used for matching the Hits");
  declareProperty("vertexCollectionName", m_vertexCollectionName = "PrimaryVertices", 
              "Name of PV collection. Default is PrimaryVertices");
  declareProperty("outputPixHitCollectionName", m_outputPixHitCollectionName = "JetAssociatedPixelClusters", 
              "Name of output pixel hit collection. JetAssociatedPixelClusters is default.");
  declareProperty("outputSCTHitCollectionName", m_outputSCTHitCollectionName = "JetAssociatedSCTClusters", 
              "Name of output SCT hit collection. JetAssociatedSCTClusters is default.");
  declareProperty("inputPixHitCollectionName", m_inputPixHitCollectionName = "PixelClusters", 
              "Name of input pixel hit collection. PixelClusters is default");
  declareProperty("inputSCTHitCollectionName", m_inputSCTHitCollectionName = "SCT_Clusters", 
              "Name of input SCT hit collection. SCT_Clusters is default");
  declareProperty("jetPtThreshold", m_jetPtThreshold = 250000.0, 
              "Hits are saved only if they match to jets with pT > jetPtThreshold [in MeV]");
  declareProperty("dRmatchHitToJet", m_dRmatchHitToJet = 0.4, 
              "The radius used for matching hits to jets");
          
  declare(m_outputPixHitCollectionName);
  declare(m_outputSCTHitCollectionName);      
}

Member Function Documentation

declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  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());
 
  }

declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  t, const SG::NotHandleType &    )

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(t);
  }

declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleBase &  hndl, const std::string &  doc, const SG::VarHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKey &  hndl, const std::string &  doc, const SG::VarHandleKeyType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKeyArray &  hndArr, const std::string &  doc, const SG::VarHandleKeyArrayType &    )

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc, const SG::NotHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(name, property, doc);
  }

declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc = "none"  )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

  {
    typedef typename SG::HandleClassifier<T>::type htype;
    return declareProperty (name, property, doc, htype());
  }

declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T > &  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< 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() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< 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; }

evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode JetHitAssociation::execute ( )

overridevirtual

Definition at line 66 of file JetHitAssociation.cxx.

                                      {
  
  const EventContext& ctx = Gaudi::Hive::currentContext();
  
  // All jets
  SG::ReadHandle<xAOD::JetContainer> jetReadHandle(m_jetCollectionName, ctx);
  if ( !jetReadHandle.isValid() ) {
    ATH_MSG_ERROR("Failed to retrieve jet container with key " << m_jetCollectionName.key() );
    return StatusCode::FAILURE;
  }
  
  // Keep only jets above threashold
  std::vector<const xAOD::Jet*> jets;
  jets.reserve( jetReadHandle->size() );
  for (const xAOD::Jet *jet : *jetReadHandle) {
    if (jet->pt() > m_jetPtThreshold) jets.push_back(jet);
  }
  
// Create write handles
  SG::WriteHandle<xAOD::TrackMeasurementValidationContainer> outputPixHits(m_outputPixHitCollectionName, ctx);
  SG::WriteHandle<xAOD::TrackMeasurementValidationContainer> outputSCTHits(m_outputSCTHitCollectionName, ctx);
 
  
  // Vertex collection
  SG::ReadHandle<xAOD::VertexContainer> vertexReadHandle(m_vertexCollectionName, ctx);
  if ( !vertexReadHandle.isValid() ) {
    ATH_MSG_ERROR("Failed to retrieve PrimaryVertices container" );
    return StatusCode::FAILURE;
  }
  
  // Find primary vertex
  const xAOD::Vertex *primVtx = nullptr;
  for (const xAOD::Vertex *vertex : *vertexReadHandle) {
    if (vertex->vertexType() == xAOD::VxType::PriVtx) { 
      primVtx = vertex;
      break;
    }
  }
  
  // Read pixel hits
  SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> pixHitReadHandle(m_inputPixHitCollectionName, ctx);
  if ( !pixHitReadHandle.isValid() ) {
    ATH_MSG_ERROR("Failed to retrieve pixel hit container with key " << m_inputPixHitCollectionName);
    return StatusCode::FAILURE;
  }
  
  // Read SCT hits
  SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> sctHitReadHandle(m_inputSCTHitCollectionName, ctx);
  if ( !sctHitReadHandle.isValid() ) {
    ATH_MSG_ERROR("Failed to retrieve SCT hit container with key " << m_inputSCTHitCollectionName);
    return StatusCode::FAILURE;
  }
  
  // These subroutines check for viable jets and write out hits if any exist. 
  // A JetAssociatedPixelCluster or JetAssociatedSCTCluster container is created for each event regardless
  ATH_CHECK( saveHits(jets, primVtx, pixHitReadHandle, outputPixHits, m_nStoredPixHits, m_nTotalPixHits) );
  ATH_CHECK( saveHits(jets, primVtx, sctHitReadHandle, outputSCTHits, m_nStoredSCTHits, m_nTotalSCTHits) );
  
  return StatusCode::SUCCESS;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< 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 & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.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 Algorithm::extraOutputDeps();
}

finalize()

StatusCode JetHitAssociation::finalize ( )

overridevirtual

Definition at line 190 of file JetHitAssociation.cxx.

                                       {
 
  ATH_MSG_DEBUG("JetHitAssociation: Total number of Pix Hits tested = " << m_nTotalPixHits);
  ATH_MSG_DEBUG("JetHitAssociation: Total number of Pix Hits Stored = " << m_nStoredPixHits);
  ATH_MSG_DEBUG("JetHitAssociation: Total number of SCT Hits tested = " << m_nTotalSCTHits);
  ATH_MSG_DEBUG("JetHitAssociation: Total number of SCT Hits Stored = " << m_nStoredSCTHits);
  
  return StatusCode::SUCCESS;
}

initialize()

StatusCode JetHitAssociation::initialize ( )

overridevirtual

Definition at line 39 of file JetHitAssociation.cxx.

                                         {
 
  ATH_CHECK( AthAlgorithm::initialize() );
  
  ATH_MSG_DEBUG("jetCollectionName           = " << m_jetCollectionName);
  ATH_MSG_DEBUG("vertexCollectionName        = " << m_vertexCollectionName);
  ATH_MSG_DEBUG("inputPixHitCollectionName   = " << m_inputPixHitCollectionName);
  ATH_MSG_DEBUG("inputSCTHitCollectionName   = " << m_inputSCTHitCollectionName);
  ATH_MSG_DEBUG("outputPixHitCollectionName  = " << m_outputPixHitCollectionName);
  ATH_MSG_DEBUG("outputSCTHitCollectionName  = " << m_outputSCTHitCollectionName);
  ATH_MSG_DEBUG("jetPtThreshold [MeV]         = " << m_jetPtThreshold);
  ATH_MSG_DEBUG("dRmatchHitToJet             = " << m_dRmatchHitToJet);
  
  ATH_CHECK(m_inputPixHitCollectionName.initialize());
  ATH_CHECK(m_inputSCTHitCollectionName.initialize());
  ATH_CHECK(m_jetCollectionName.initialize());
  ATH_CHECK(m_vertexCollectionName.initialize());
  ATH_CHECK(m_outputPixHitCollectionName.initialize());
  ATH_CHECK(m_outputSCTHitCollectionName.initialize());
 
  ATH_MSG_DEBUG("JetHitAssociation: Initialized");
  
  return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< 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.

msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

bool AthCommonMsg< 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< 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< Algorithm > >::renounce ( T &  h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

saveHits()

StatusCode JetHitAssociation::saveHits ( const std::vector< const xAOD::Jet * > &  jets, const xAOD::Vertex *const  vertex, SG::ReadHandle< xAOD::TrackMeasurementValidationContainer > &  hits, const SG::WriteHandle< xAOD::TrackMeasurementValidationContainer > &  writeHandle, unsigned long long int &  nStoredHits, unsigned long long int &  nTotalHits  )

private

Definition at line 128 of file JetHitAssociation.cxx.

                                                           {
 
  // Collection where output hits are stored
  auto outputHitCollection = std::make_unique<xAOD::TrackMeasurementValidationContainer>();
  auto outputHitCollectionAux = std::make_unique<xAOD::TrackMeasurementValidationAuxContainer>();
  outputHitCollection->setStore(outputHitCollectionAux.get());
    
    
  // if there are no jets you need to skip this loop but you still need to write out the 
  // if the vertex is a null pointer (no primary vertex in the event) do the same
  // JetAssociatedSCTClusters and JetAssociatedPixelClusters
  if (!jets.empty() && vertex!=nullptr) {
 
    // Get the x,y,z of the primary vertex
    TVector3 PVposition(vertex->x(), vertex->y(), vertex->z());   
    ATH_MSG_DEBUG("JetHitAssociation: PrimaryVertex Z Position = " << PVposition.Z());
 
    // Loop over hits
    for (const xAOD::TrackMeasurementValidation *hit : *hits) {
      // Get pixel globalX,Y,Z
      float x = hit->globalX();
      float y = hit->globalY();
      float z = hit->globalZ();
 
      // Put the x,y,z in a TVector3
      TVector3 hitPosition(x,y,z);
 
      // Correct for the PV position
      hitPosition = hitPosition - PVposition;
 
      // Loop over jets
      for (const xAOD::Jet *jet : jets) {
    // Calculate dR(hit,jet)
    float dR = hitPosition.DeltaR(jet->p4().Vect());
      
    // if the dR is smaller than the dR association threshold save the hit
    if (dR < m_dRmatchHitToJet) {
      outputHitCollection->push_back(std::make_unique<xAOD::TrackMeasurementValidation>(*hit));
      *outputHitCollection->back() = *hit; //Without this line. entries are filled with zero instead of the correct values
      ++nStoredHits;
      break;
    }
      }   
    }
  }  
  // count the total number of hits run
  nTotalHits += hits->size(); 
    
  // Write output container
  if ( ! outputHits.put( std::move(outputHitCollection), std::move(outputHitCollectionAux) ) ) {
    ATH_MSG_ERROR("Could not write output hit containers");
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
}

sysInitialize()

StatusCode AthAlgorithm::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< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

                                       {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::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< 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< 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_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_dRmatchHitToJet

float JetHitAssociation::m_dRmatchHitToJet

private

Definition at line 56 of file JetHitAssociation.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

m_inputPixHitCollectionName

SG::ReadHandleKey<xAOD::TrackMeasurementValidationContainer> JetHitAssociation::m_inputPixHitCollectionName

private

Definition at line 49 of file JetHitAssociation.h.

m_inputSCTHitCollectionName

SG::ReadHandleKey<xAOD::TrackMeasurementValidationContainer> JetHitAssociation::m_inputSCTHitCollectionName

private

Definition at line 50 of file JetHitAssociation.h.

m_jetCollectionName

SG::ReadHandleKey<xAOD::JetContainer> JetHitAssociation::m_jetCollectionName

private

Definition at line 47 of file JetHitAssociation.h.

m_jetPtThreshold

float JetHitAssociation::m_jetPtThreshold

private

Definition at line 55 of file JetHitAssociation.h.

m_nStoredPixHits

unsigned long long int JetHitAssociation::m_nStoredPixHits

private

Definition at line 57 of file JetHitAssociation.h.

m_nStoredSCTHits

unsigned long long int JetHitAssociation::m_nStoredSCTHits

private

Definition at line 58 of file JetHitAssociation.h.

m_nTotalPixHits

unsigned long long int JetHitAssociation::m_nTotalPixHits

private

Definition at line 59 of file JetHitAssociation.h.

m_nTotalSCTHits

unsigned long long int JetHitAssociation::m_nTotalSCTHits

private

Definition at line 60 of file JetHitAssociation.h.

m_outputPixHitCollectionName

SG::WriteHandleKey<xAOD::TrackMeasurementValidationContainer> JetHitAssociation::m_outputPixHitCollectionName

private

Definition at line 51 of file JetHitAssociation.h.

m_outputSCTHitCollectionName

SG::WriteHandleKey<xAOD::TrackMeasurementValidationContainer> JetHitAssociation::m_outputSCTHitCollectionName

private

Definition at line 52 of file JetHitAssociation.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vertexCollectionName

SG::ReadHandleKey<xAOD::VertexContainer> JetHitAssociation::m_vertexCollectionName

private

Definition at line 48 of file JetHitAssociation.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• JetHitAssociation.h
• JetHitAssociation.cxx