InDet::InDetBeamSpotFinder Class Reference

#include <InDetBeamSpotFinder.h>

Inheritance diagram for InDet::InDetBeamSpotFinder:

Classes
struct	beamSpotNtuple_struct
struct	vertexNtuple_struct

Public Member Functions
	InDetBeamSpotFinder (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode	initialize () override
virtual StatusCode	execute () override
virtual StatusCode	finalize () override
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
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	passEventSelection (const xAOD::EventInfo &)
BeamSpot::Event	readEvent (const xAOD::EventInfo &, const xAOD::VertexContainer &)
void	writeToVertexTree (BeamSpot::Event &, BeamSpot::VrtHolder &)
void	writeToBeamSpotTree (const IInDetBeamSpotTool *bs, std::vector< BeamSpot::Event > &, std::vector< BeamSpot::VrtHolder > &)
void	sortEvents ()
StatusCode	setupVertexTree ()
StatusCode	setupBeamSpotTree ()
StatusCode	performFits ()
bool	iequals (const std::string &, const std::string &)
int	min_lb (std::vector< BeamSpot::Event > &)
int	max_lb (std::vector< BeamSpot::Event > &)
int	min_run (std::vector< BeamSpot::Event > &)
int	max_run (std::vector< BeamSpot::Event > &)
void	convertVtxTypeNames ()
IInDetBeamSpotTool *	cloneTool (int)
bool	passVertexSelection (const xAOD::Vertex *)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ServiceHandle< IToolSvc >	m_toolSvc {this, "ToolSvc", "ToolSvc"}
ServiceHandle< ITHistSvc >	m_thistSvc {this, "THistSvc", "THistSvc"}
ToolHandleArray< IInDetBeamSpotTool >	m_beamSpotToolList
SG::ReadCondHandleKey< BunchCrossingCondData >	m_bcDataKey {this, "BunchCrossingCondDataKey", "BunchCrossingData" ,"SG Key of BunchCrossing CDO"}
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfo {this, "EvtInfo", "EventInfo", "EventInfo name"}
SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexContainer {this, "VertexContainer", "PrimaryVertices", "Vertex container name"}
unsigned int	m_maxRunsPerFit
unsigned int	m_maxEventsPerFit
unsigned int	m_maxLBsPerFit
std::vector< unsigned int >	m_BCIDsToAccept
bool	m_writeVertexNtuple
std::string	m_beamSpotNtupleName
TTree *	m_root_bs =0
TTree *	m_root_vrt =0
beamSpotNtuple_struct	m_beamSpotNtuple
vertexNtuple_struct	m_root_vtx
BeamSpotStatusCode	m_BeamStatusCode
std::string	m_vertexTreeName
std::vector< BeamSpot::Event >	m_eventList
std::map< BeamSpot::ID, std::vector< BeamSpot::Event > >	m_eventMap
std::vector< std::vector< BeamSpot::Event > >	m_sortedEventList
std::vector< std::vector< BeamSpot::VrtHolder > >	m_sortedVertices
int	m_minTrackNum {}
int	m_maxTrackNum {}
double	m_maxChi2Vertex {}
double	m_maxTransverseError {}
double	m_maxAbsCorrelXY {}
double	m_minVtxProb {}
unsigned int	m_minVertexNum {}
std::vector< std::string >	m_vertexTypeNames {}
std::vector< xAOD::VxType::VertexType >	m_vertexTypes {}
bool	m_useFilledBCIDsOnly {}
std::string	m_fitSortingKey {}
bool	m_writeAllVertices {}
unsigned long	m_secondsPerFit {}
DataObjIDColl	m_extendedExtraObjects
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

Definition at line 35 of file InDetBeamSpotFinder.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

InDetBeamSpotFinder()

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

Definition at line 30 of file InDetBeamSpotFinder.cxx.

                                                                                            :
  AthAlgorithm(name, pSvcLocator)
{
  declareProperty( "BeamSpotToolList"  , m_beamSpotToolList );
  declareProperty( "RunRange"     , m_maxRunsPerFit  = 0 );
  declareProperty( "LumiRange"      , m_maxLBsPerFit = 0 );
  declareProperty( "EventRange"   , m_maxEventsPerFit = 0);
  declareProperty( "UseBCID"     , m_BCIDsToAccept );
  declareProperty( "UseFilledBCIDsOnly", m_useFilledBCIDsOnly = true );
  declareProperty( "MinTracksPerVtx", m_minTrackNum = 5);
  declareProperty( "MaxTracksPerVtx", m_maxTrackNum = 1000000);
  declareProperty( "MinVtxNum"      , m_minVertexNum = 100);
  declareProperty( "MaxVtxChi2"     , m_maxChi2Vertex = 10);
  declareProperty( "MaxTransverseErr", m_maxTransverseError=1000000);
  declareProperty( "MaxAbsCorrelXY", m_maxAbsCorrelXY=0.8);
  declareProperty( "VertexTypes"       , m_vertexTypeNames);
  declareProperty( "MinVtxProb"     , m_minVtxProb=0.001);
  declareProperty( "GroupFitsBy"       , m_fitSortingKey = "none");
  declareProperty( "VertexNtuple" , m_writeVertexNtuple = true);
  declareProperty( "WriteAllVertices"  , m_writeAllVertices=false);
  declareProperty( "VertexTreeName"    , m_vertexTreeName = "Vertices");
  declareProperty( "SecondsPerFit",  m_secondsPerFit = 1);
}

Member Function Documentation

cloneTool()

InDet::IInDetBeamSpotTool * InDet::InDetBeamSpotFinder::cloneTool ( int i )

private

Definition at line 446 of file InDetBeamSpotFinder.cxx.

                                                                  {
  IInDetBeamSpotTool  * orig = &(*m_beamSpotToolList[i]);
  IInDetBeamSpotTool * temp = orig->Clone();
  return temp;
}

convertVtxTypeNames()

void InDet::InDetBeamSpotFinder::convertVtxTypeNames ( )

private

Definition at line 207 of file InDetBeamSpotFinder.cxx.

                                                  {
  if ( not m_vertexTypeNames.empty() ) {
    for ( std::vector<std::string>::const_iterator it = m_vertexTypeNames.begin();
    it != m_vertexTypeNames.end(); ++it) {
      if ((*it) == "NoVtx") ;
      else if ((*it) == "PriVtx") m_vertexTypes.push_back(xAOD::VxType::PriVtx);
      else if ((*it) == "SecVtx") m_vertexTypes.push_back(xAOD::VxType::SecVtx);
      else if ((*it) == "PileUp") m_vertexTypes.push_back(xAOD::VxType::PileUp);
      else if ((*it) == "ConvVtx") m_vertexTypes.push_back(xAOD::VxType::ConvVtx);
      else if ((*it) == "V0Vtx") m_vertexTypes.push_back(xAOD::VxType::V0Vtx);
      else if ((*it) == "KinkVtx") m_vertexTypes.push_back(xAOD::VxType::KinkVtx);
      else if ((*it) =="NotSpecified")m_vertexTypes.push_back(xAOD::VxType::NotSpecified) ;
    }
    ATH_MSG_INFO("Allowing " << m_vertexTypes.size() << " Vertex types" );
  }
  else {
    ATH_MSG_DEBUG( "No selection based on vertexType will be done" );
  }
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< 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< 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< 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< 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 InDet::InDetBeamSpotFinder::execute ( )

overridevirtual

Definition at line 76 of file InDetBeamSpotFinder.cxx.

                                            {
  SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInfo);
  SG::ReadHandle<xAOD::VertexContainer> vertexContainer(m_vertexContainer);
  if ( !passEventSelection( *eventInfo ) ) return StatusCode::SUCCESS;
  BeamSpot::Event currentEvent = readEvent(*eventInfo, *vertexContainer);
 
  m_eventList.push_back( currentEvent );
  if( m_writeVertexNtuple ){
    for( auto & thisVertex: currentEvent.vertices){
      if( thisVertex.passed || m_writeAllVertices ){
        writeToVertexTree( currentEvent, thisVertex );
      }
    }
  }
  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 InDet::InDetBeamSpotFinder::finalize ( )

overridevirtual

Definition at line 93 of file InDetBeamSpotFinder.cxx.

                                              {
  ATH_MSG_DEBUG( "in finalize()" );
  sortEvents();
  ATH_CHECK(performFits());
  return StatusCode::SUCCESS;
}

iequals()

bool InDet::InDetBeamSpotFinder::iequals ( const std::string & a, const std::string & b )

private

Definition at line 381 of file InDetBeamSpotFinder.cxx.

{
    unsigned int sz = a.size();
    if (b.size() != sz)
        return false;
    for (unsigned int i = 0; i < sz; ++i)
        if (tolower(a[i]) != tolower(b[i]))
            return false;
    return true;
}

initialize()

StatusCode InDet::InDetBeamSpotFinder::initialize ( )

overridevirtual

Definition at line 54 of file InDetBeamSpotFinder.cxx.

                                                {
  ATH_MSG_DEBUG( "in initialize()");
  if ( m_beamSpotToolList.empty() ){
    ATH_MSG_FATAL("FATAL ERROR: must provide at least one beamspot tool in beamSpotToolList");
    return StatusCode::FAILURE;
  }
 
  ATH_CHECK( m_thistSvc.retrieve() );
  ATH_CHECK( m_toolSvc.retrieve() );
 
  ATH_CHECK( m_bcDataKey.initialize(m_useFilledBCIDsOnly) );
 
  ATH_CHECK( m_eventInfo.initialize() );
  ATH_CHECK( m_vertexContainer.initialize() );
 
  for ( unsigned int i = 0; i < m_beamSpotToolList.size(); i++){ ATH_CHECK( m_beamSpotToolList[i].retrieve() );}
  if( m_writeVertexNtuple ){ ATH_CHECK(setupVertexTree()); }
  ATH_CHECK( setupBeamSpotTree() );
  convertVtxTypeNames();
  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.

max_lb()

int InDet::InDetBeamSpotFinder::max_lb ( std::vector< BeamSpot::Event > & eventList )

private

Definition at line 457 of file InDetBeamSpotFinder.cxx.

                                                                          {
  const auto largestLbEvent=std::max_element(eventList.begin(),eventList.end(), lesserLumiBlock);
  return largestLbEvent->lumiBlock;
}

max_run()

int InDet::InDetBeamSpotFinder::max_run ( std::vector< BeamSpot::Event > & eventList )

private

Definition at line 467 of file InDetBeamSpotFinder.cxx.

                                                                           {
  const auto largestRunEvent=std::max_element(eventList.begin(),eventList.end(), lesserRunNumber);
  return largestRunEvent->runNumber;
}

min_lb()

int InDet::InDetBeamSpotFinder::min_lb ( std::vector< BeamSpot::Event > & eventList )

private

Definition at line 452 of file InDetBeamSpotFinder.cxx.

                                                                          {
  const auto smallestLbEvent=std::min_element(eventList.begin(),eventList.end(), lesserLumiBlock);
  return smallestLbEvent->lumiBlock;
}

min_run()

int InDet::InDetBeamSpotFinder::min_run ( std::vector< BeamSpot::Event > & eventList )

private

Definition at line 462 of file InDetBeamSpotFinder.cxx.

                                                                           {
  const auto smallestRunEvent=std::min_element(eventList.begin(),eventList.end(), lesserRunNumber);
  return smallestRunEvent->runNumber;
}

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

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.

passEventSelection()

bool InDet::InDetBeamSpotFinder::passEventSelection ( const xAOD::EventInfo & eventInfo )

private

Definition at line 227 of file InDetBeamSpotFinder.cxx.

                                                                                {
  const int bcid = eventInfo.bcid();
  if ( m_useFilledBCIDsOnly) {
    SG::ReadCondHandle<BunchCrossingCondData> bcData(m_bcDataKey);
    if ( !bcData->isFilled(bcid) ) return false;
  }
  if( m_BCIDsToAccept.begin() !=  m_BCIDsToAccept.end() )
    return ( std::find(m_BCIDsToAccept.begin(), m_BCIDsToAccept.end(), bcid) != m_BCIDsToAccept.end());
  else 
    return true;
}

passVertexSelection()

bool InDet::InDetBeamSpotFinder::passVertexSelection ( const xAOD::Vertex * vtx )

private

Definition at line 239 of file InDetBeamSpotFinder.cxx.

                                                                          {
  if(!vtx)                                                     { return false; }
  if(m_vertexTypes.end() == std::find(  m_vertexTypes.begin(), m_vertexTypes.end(), vtx->vertexType() ) ) { return false; }
  if(vtx->chiSquared()/vtx->numberDoF() > m_maxChi2Vertex)     { return false; }
  if(static_cast<int>(vtx->nTrackParticles()) < m_minTrackNum ){ return false; }
  if(static_cast<int>(vtx->nTrackParticles()) > m_maxTrackNum ){ return false; }
  if(TMath::Prob(vtx->chiSquared(), vtx->numberDoF()) < m_minVtxProb ) { return false; }
  if(vtx->covariancePosition()(0,0) <= 0 || vtx->covariancePosition()(1,1) <= 0 || vtx->covariancePosition()(2,2) <= 0 ) { return false; }
  double maxTransverseError2 = m_maxTransverseError * m_maxTransverseError;
  if(vtx->covariancePosition()(0,0) > maxTransverseError2 || vtx->covariancePosition()(1,1) > maxTransverseError2) {return false;}
  return true;
}

performFits()

StatusCode InDet::InDetBeamSpotFinder::performFits ( )

private

Definition at line 260 of file InDetBeamSpotFinder.cxx.

                                                {
  IInDetBeamSpotTool::FitStatus bsFitStatus;
  std::vector<BeamSpot::VrtHolder> verticesToFit;
 
  for( auto & eventList: m_sortedEventList){
    verticesToFit.clear();
    if(eventList.size() > 0) {
      BeamSpot::Event ev = *eventList.begin();
      ATH_MSG_INFO( "Event list size/LB/bcid: " << eventList.size() << " " << ev.lumiBlock << " " << ev.bcid);
    }
    for( const auto & thisEvent: eventList){
      for( const auto & thisVertex: thisEvent.vertices){
        if( thisVertex.passed ) { 
      verticesToFit.push_back( thisVertex ); 
    }
      }
    }
 
    for( unsigned int j = 0; j < m_beamSpotToolList.size(); j++){
      IInDetBeamSpotTool * bs(nullptr);
      bs = cloneTool(j);
      if(!bs){ return StatusCode::FAILURE; }
      if(not verticesToFit.empty()) { bsFitStatus = bs->fit(verticesToFit); }
      else { bsFitStatus = IInDetBeamSpotTool::unsolved; }
      
      m_BeamStatusCode.clearWord();
      m_BeamStatusCode.setOnlineStatus(false);
      m_BeamStatusCode.setAlgType( bs->getFitID() );
      int fitStat;
      if ( bsFitStatus == IInDetBeamSpotTool::successful)
        fitStat = 3;
      else if ( bsFitStatus == IInDetBeamSpotTool::problems)
        fitStat = 1;
      else if ( bsFitStatus == IInDetBeamSpotTool::failed || bsFitStatus == IInDetBeamSpotTool::unsolved)
        fitStat = 0;  
      else 
        fitStat = 0;
      m_BeamStatusCode.setFitStatus(fitStat);
      if (bs->getParamMap()["sigmaX"] == 0 && bs->getParamMap()["sigmaY"] ==0 ) { m_BeamStatusCode.setFitWidth( false); }
      else { m_BeamStatusCode.setFitWidth(true); } 
 
      if(not eventList.empty()) writeToBeamSpotTree( bs, eventList, verticesToFit );
    }
  }
  return StatusCode::SUCCESS;
}

readEvent()

BeamSpot::Event InDet::InDetBeamSpotFinder::readEvent ( const xAOD::EventInfo & eventInfo, const xAOD::VertexContainer & vertexContainer )

private

Definition at line 100 of file InDetBeamSpotFinder.cxx.

                                                                                                                              {
  BeamSpot::Event event;
  BeamSpot::VrtHolder vertex;
  event.pileup = ceil(eventInfo.actualInteractionsPerCrossing());
  event.runNumber = eventInfo.runNumber();
  event.lumiBlock  = eventInfo.lumiBlock();
  event.bcid  = eventInfo.bcid();
  event.eventTime    = eventInfo.timeStamp();
  event.eventTime_NS = eventInfo.timeStampNSOffset();
  event.eventNumber  = eventInfo.eventNumber();
  const EventInfo* BSeventInfo; 
  //This is required for pseudo lumiblocks 
 
  if(!eventInfo.eventType(xAOD::EventInfo::EventType::IS_SIMULATION))
    {
      if( evtStore()->retrieve(BSeventInfo) != StatusCode::SUCCESS){
    ATH_MSG_ERROR("Cannot get event info.");
    return event;
      }
 
      // Remove the redundent logic, event.lumiBlock != BSeventInfo->event_ID()->lumi_block()
      event.lumiBlock = BSeventInfo->event_ID()->lumi_block();
    };
 
  for(const xAOD::Vertex* vtx:vertexContainer) {
    if (vtx->vertexType() == xAOD::VxType::NoVtx) continue; 
    vertex.x          = vtx->x();
    vertex.y          = vtx->y();
    vertex.z          = vtx->z();
    vertex.vxx        = vtx->covariancePosition()(0,0);
    vertex.vxy        = vtx->covariancePosition()(0,1);
    vertex.vyy        = vtx->covariancePosition()(1,1);
    vertex.vzz        = vtx->covariancePosition()(2,2);
    vertex.vertexType = vtx->vertexType();
    vertex.nTracks    = vtx->nTrackParticles();
    vertex.passed     = passVertexSelection( vtx );
    vertex.valid      = vertex.passed;    
    //Remove vertices with wrong x-y correlation
    if(vertex.vxy*vertex.vxy/vertex.vxx/vertex.vyy > 1.0 ||
       vertex.vzz < 0. ||
       vertex.vyy < 0. ||
       vertex.vxx < 0.) {
      ATH_MSG_DEBUG("Bad vertex: " << event.eventNumber << " " << vtx->vertexType() << " "
           << vertex.x << " " << vertex.y << " " << vertex.z << " "
           << vertex.vxx << " " << vertex.vyy << " " << vertex.vzz << " " 
           << vertex.nTracks << " " << vertex.vxy*vertex.vxy/vertex.vxx/vertex.vyy);
    }
    if(vertex.vxy*vertex.vxy/vertex.vxx/vertex.vyy > m_maxAbsCorrelXY*m_maxAbsCorrelXY ||
       vertex.vzz < 0. ||
       vertex.vyy < 0. ||
       vertex.vxx < 0.) continue;
    event.vertices.push_back( vertex );
  }
  return event;
}

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();
  }

setupBeamSpotTree()

StatusCode InDet::InDetBeamSpotFinder::setupBeamSpotTree ( )

private

Definition at line 307 of file InDetBeamSpotFinder.cxx.

                                                      {
  const std::string inRootID = "/INDETBEAMSPOTFINDER/";
  const std::string sbs       = "BeamSpotNt";//m_root_beamspotName;
  m_root_bs = new TTree(sbs.data(),"Beamspot Solutions");
  m_root_bs->Branch("bcid",          &m_beamSpotNtuple.bcid,"bcid/I");
  m_root_bs->Branch("pileup",        &m_beamSpotNtuple.pileup,"pileup/I");
  m_root_bs->Branch("defectWord",    &m_beamSpotNtuple.defectWord, "defectWord/I");
  m_root_bs->Branch("fill",          &m_beamSpotNtuple.fill,       "fill/I");
  m_root_bs->Branch("lbEnd",         &m_beamSpotNtuple.lbEnd,      "lbEnd/I");
  m_root_bs->Branch("lbStart",       &m_beamSpotNtuple.lbStart,    "lbStart/I");
  m_root_bs->Branch("nEvents",       &m_beamSpotNtuple.nEvents,    "nEvents/I");
  m_root_bs->Branch("nValid",        &m_beamSpotNtuple.nValid,     "nValid/I");
  m_root_bs->Branch("nVtxAll",       &m_beamSpotNtuple.nVtxAll,    "nVtxAll/I");
  m_root_bs->Branch("nVtxPrim",      &m_beamSpotNtuple.nVtxPrim,   "nVtxPrim/I");
  m_root_bs->Branch("separation",    &m_beamSpotNtuple.separation, "separation/I");
  m_root_bs->Branch("status",        &m_beamSpotNtuple.status,        "status/I");
  m_root_bs->Branch("timeEnd",       &m_beamSpotNtuple.timeEnd,       "timeEnd/I");
  m_root_bs->Branch("timeStart",     &m_beamSpotNtuple.timeStart,     "timeStart/I");
  m_root_bs->Branch("run",           &m_beamSpotNtuple.run,           "run/I");
  m_root_bs->Branch("runEnd",        &m_beamSpotNtuple.runEnd,        "runEnd/I");
  
  for( auto &tool : m_beamSpotToolList ){
    std::map<std::string,double> paramMap = tool->getParamMap();
    std::map<std::string,double> covMap   = tool->getCovMap();
 
    std::string slashD = "/D";
    std::string keySlashD; 
  
    //Loop over the parameters for a given fit tool, and create a branch for each if it doesn't exist
    for( std::map<std::string,double>::iterator iter = paramMap.begin(); iter != paramMap.end(); ++iter){
      std::string key =  iter->first;
      //double val =  iter->second;
      if( !(m_root_bs->GetBranch(key.c_str())) ){
        m_beamSpotNtuple.paramMap[key] = 0;
        keySlashD = key + slashD;
        m_root_bs->Branch(key.c_str(), &m_beamSpotNtuple.paramMap[key], keySlashD.c_str());
      }
    }
    //Loop over the covariance matrix for a given fit tool and create a branch for each element, if it doesn't already exist.
    for( std::map<std::string,double>::iterator iter = covMap.begin(); iter != covMap.end(); ++iter){
      const std::string & key =  iter->first;
      //double val  =  iter->second;
      if( !(m_root_bs->GetBranch(key.c_str())) ){
        m_beamSpotNtuple.covMap[key] = 0;
        keySlashD = key + slashD;
        m_root_bs->Branch( key.c_str(), &m_beamSpotNtuple.covMap[key], keySlashD.c_str());
      }
    }
  }
  return m_thistSvc->regTree(inRootID+sbs,m_root_bs);
}

setupVertexTree()

StatusCode InDet::InDetBeamSpotFinder::setupVertexTree ( )

private

Definition at line 252 of file InDetBeamSpotFinder.cxx.

                                                    {
  const std::string inRootID = "/INDETBEAMSPOTFINDER/";
  const std::string svrts       = m_vertexTreeName;
  m_root_vrt = new TTree(svrts.data(),"Vertices");
  m_root_vrt->Branch("vrt",&m_root_vtx,"x/D:y:z:vxx:vxy:vyy:vzz:vType/i:run:lb:bcid:pileup:nTracks:eventNumber/l:eventTime:eventTime_NS:passed/O:valid");
  return m_thistSvc->regTree(inRootID+svrts,m_root_vrt);
}

sortEvents()

void InDet::InDetBeamSpotFinder::sortEvents ( )

private

Definition at line 156 of file InDetBeamSpotFinder.cxx.

                                         {
  for( const auto & thisEvent: m_eventList){
    BeamSpot::ID id;
    id.runNumber( (m_maxRunsPerFit > 0) ? thisEvent.runNumber : 0 );
    id.lumiBlock( (m_maxLBsPerFit  > 0) ? thisEvent.lumiBlock : 0 );
    id.pileup   ( iequals(m_fitSortingKey,"pileup") ? thisEvent.pileup : 0 );
    id.bcid     ( iequals(m_fitSortingKey,"bcid"  ) ? thisEvent.bcid   : 0 );
    id.timeStamp( iequals(m_fitSortingKey,"time"  ) ? thisEvent.eventTime/m_secondsPerFit  : 0 );
    m_eventMap[id].push_back( thisEvent );
  }
  auto iter = m_eventMap.begin();
  BeamSpot::ID lastID = iter->first;
  BeamSpot::ID currentID = iter->first;
  unsigned int nRuns = 0;
  unsigned int nLBs = 0;
  unsigned int nFits = 1;
  m_sortedEventList.resize( nFits );
  
  for( iter = m_eventMap.begin(); iter != m_eventMap.end(); ++iter){
    currentID = iter->first;
    if( iter == m_eventMap.begin() || currentID.runNumber() != lastID.runNumber() ){ nRuns++; }
    if( iter == m_eventMap.begin() || currentID.lumiBlock() != lastID.lumiBlock() ){ nLBs++;  }
    if( currentID.timeStamp() != lastID.timeStamp() || 
        currentID.pileup() != lastID.pileup() 
        || currentID.bcid() != lastID.bcid() 
        || ( m_maxRunsPerFit > 0 && nRuns > m_maxRunsPerFit )
        || ( m_maxLBsPerFit  > 0 && nLBs  > m_maxLBsPerFit )){
      ATH_MSG_INFO( "New Fit " << currentID.timeStamp() << "<--"<<  lastID.timeStamp() << "\n\t" 
        << currentID.pileup() << " <-- " << lastID.pileup() << "\n\t" 
        << currentID.bcid() << " <-- " <<  lastID.bcid() << "\n\t"
        << nRuns << " " << m_maxRunsPerFit << "\n\t"
        << nLBs << "  " << m_maxLBsPerFit );
     
      nFits++;
      m_sortedEventList.resize(nFits);
      nRuns = 1; nLBs = 1;
    }
    for( unsigned int i = 0; i < iter->second.size(); i++){
      if( m_sortedEventList.at(nFits-1).size() == m_maxEventsPerFit && m_maxEventsPerFit > 0 ){
        nFits++;
        m_sortedEventList.resize(nFits);
        nRuns = 1; nLBs = 1;
      }
      
      m_sortedEventList.at(nFits-1).push_back( iter->second.at(i) );
    }
    lastID = iter->first;
  }
}

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, AthHistogramAlgorithm, and PyAthena::Alg.

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);
      }
    }
  }

writeToBeamSpotTree()

void InDet::InDetBeamSpotFinder::writeToBeamSpotTree ( const IInDetBeamSpotTool * bs, std::vector< BeamSpot::Event > & eventList, std::vector< BeamSpot::VrtHolder > & vertexList )

private

leave these here; illustrates coverity defect 29456 for( std::map<std::string,double>::iterator iter = m_beamSpotNtuple.paramMap.begin(); iter != m_beamSpotNtuple.paramMap.end(); ++iter){
std::string key = iter->first; iter->second = ( bs->getParamMap().find(key) == bs->getParamMap().end() ) ? 0 : bs->getParamMap()[key]; }

leave these here; illustrates coverity defect 29456 for( std::map<std::string,double>::iterator iter = m_beamSpotNtuple.covMap.begin(); iter != m_beamSpotNtuple.covMap.end(); ++iter){
std::string key = iter->first; iter->second = ( bs->getCovMap().find(key) == bs->getCovMap().end() ) ? 0 : bs->getCovMap()[key]; }

Definition at line 392 of file InDetBeamSpotFinder.cxx.

                                                                                                                                                             {
  m_beamSpotNtuple.pileup = iequals(m_fitSortingKey,"pileup") ? eventList.at(0).pileup : 0;
  m_beamSpotNtuple.bcid   = iequals(m_fitSortingKey,"bcid")   ? eventList.at(0).bcid   : 0;
  m_beamSpotNtuple.defectWord = 0;
  m_beamSpotNtuple.fill = 0;
  m_beamSpotNtuple.lbEnd   = max_lb( eventList );
  m_beamSpotNtuple.lbStart = min_lb( eventList );
  m_beamSpotNtuple.nEvents = eventList.size();  
  m_beamSpotNtuple.nValid = vertexList.size();
  unsigned int nVtxAll = 0;
  unsigned int nVtxPrim = 0; 
  auto isPrimaryVertex=[](const BeamSpot::VrtHolder & vertex){ return (vertex.vertexType == xAOD::VxType::PriVtx); };
  for( const auto & thisEvent : eventList){
    const auto & theseVertices=thisEvent.vertices;
    nVtxAll += theseVertices.size();
    nVtxPrim+= std::count_if(theseVertices.begin(), theseVertices.end(), isPrimaryVertex);
  }
  m_beamSpotNtuple.nVtxAll = nVtxAll;
  m_beamSpotNtuple.nVtxPrim = nVtxPrim;
  m_beamSpotNtuple.run = min_run( eventList );
  m_beamSpotNtuple.status = m_BeamStatusCode.getWord();  
  m_beamSpotNtuple.separation = 0;
  m_beamSpotNtuple.timeEnd = iequals(m_fitSortingKey,"time")   ? eventList.back().eventTime   : 0;
  m_beamSpotNtuple.timeStart = iequals(m_fitSortingKey,"time")   ? eventList.front().eventTime   : 0;
  m_beamSpotNtuple.runEnd = max_run( eventList );
  const auto & bsToolParamMap = bs->getParamMap();
  const auto & bsToolCovMap= bs->getCovMap();
  for (auto & param:m_beamSpotNtuple.paramMap){
    const std::string & key = param.first;
    const auto & bsToolEquivalent =bsToolParamMap.find(key);
    param.second = ( bsToolEquivalent == bsToolParamMap.end() ) ? 0 : bsToolEquivalent->second;
  }
  for (auto & covariance:m_beamSpotNtuple.covMap){
    const std::string & key = covariance.first;
    const auto & bsToolEquivalent = bsToolCovMap.find(key);
    covariance.second = ( bsToolEquivalent == bsToolCovMap.end() ) ? 0 : bsToolEquivalent->second;
  }
  m_root_bs->Fill();
}

writeToVertexTree()

void InDet::InDetBeamSpotFinder::writeToVertexTree ( BeamSpot::Event & evt, BeamSpot::VrtHolder & vtx )

private

Definition at line 359 of file InDetBeamSpotFinder.cxx.

                                                                                           {
  m_root_vtx.x    = vtx.x;
  m_root_vtx.y    = vtx.y;
  m_root_vtx.z    = vtx.z;
  m_root_vtx.vxx  = vtx.vxx;
  m_root_vtx.vxy  = vtx.vxy;
  m_root_vtx.vyy  = vtx.vyy;
  m_root_vtx.vzz  = vtx.vzz;
  m_root_vtx.vType = vtx.vertexType;
  m_root_vtx.run = evt.runNumber;
  m_root_vtx.lb  = evt.lumiBlock;
  m_root_vtx.bcid = evt.bcid;
  m_root_vtx.pileup = evt.pileup;
  m_root_vtx.nTracks = vtx.nTracks;
  m_root_vtx.passed = vtx.passed;
  m_root_vtx.valid        = vtx.valid;
  m_root_vtx.eventNumber  = evt.eventNumber;
  m_root_vtx.eventTime    = evt.eventTime;
  m_root_vtx.eventTime_NS = evt.eventTime_NS;
  m_root_vrt->Fill();
}

Member Data Documentation

m_bcDataKey

SG::ReadCondHandleKey<BunchCrossingCondData> InDet::InDetBeamSpotFinder::m_bcDataKey {this, "BunchCrossingCondDataKey", "BunchCrossingData" ,"SG Key of BunchCrossing CDO"}

private

Definition at line 50 of file InDetBeamSpotFinder.h.

{this, "BunchCrossingCondDataKey", "BunchCrossingData" ,"SG Key of BunchCrossing CDO"};

m_BCIDsToAccept

std::vector<unsigned int> InDet::InDetBeamSpotFinder::m_BCIDsToAccept

private

Definition at line 63 of file InDetBeamSpotFinder.h.

m_beamSpotNtuple

beamSpotNtuple_struct InDet::InDetBeamSpotFinder::m_beamSpotNtuple

private

Definition at line 87 of file InDetBeamSpotFinder.h.

m_beamSpotNtupleName

std::string InDet::InDetBeamSpotFinder::m_beamSpotNtupleName

private

Definition at line 66 of file InDetBeamSpotFinder.h.

m_beamSpotToolList

ToolHandleArray<IInDetBeamSpotTool> InDet::InDetBeamSpotFinder::m_beamSpotToolList

private

Definition at line 48 of file InDetBeamSpotFinder.h.

m_BeamStatusCode

BeamSpotStatusCode InDet::InDetBeamSpotFinder::m_BeamStatusCode

private

Definition at line 89 of file InDetBeamSpotFinder.h.

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_eventInfo

SG::ReadHandleKey<xAOD::EventInfo> InDet::InDetBeamSpotFinder::m_eventInfo {this, "EvtInfo", "EventInfo", "EventInfo name"}

private

Definition at line 53 of file InDetBeamSpotFinder.h.

{this, "EvtInfo", "EventInfo", "EventInfo name"};

m_eventList

std::vector<BeamSpot::Event> InDet::InDetBeamSpotFinder::m_eventList

private

Definition at line 92 of file InDetBeamSpotFinder.h.

m_eventMap

std::map< BeamSpot::ID, std::vector<BeamSpot::Event> > InDet::InDetBeamSpotFinder::m_eventMap

private

Definition at line 93 of file InDetBeamSpotFinder.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_fitSortingKey

std::string InDet::InDetBeamSpotFinder::m_fitSortingKey {}

private

Definition at line 131 of file InDetBeamSpotFinder.h.

{};

m_maxAbsCorrelXY

double InDet::InDetBeamSpotFinder::m_maxAbsCorrelXY {}

private

Definition at line 118 of file InDetBeamSpotFinder.h.

{}; //max absolute XY correlation of a vertex to be added 

m_maxChi2Vertex

double InDet::InDetBeamSpotFinder::m_maxChi2Vertex {}

private

Definition at line 116 of file InDetBeamSpotFinder.h.

{};    //maxChi2 of vertex

m_maxEventsPerFit

unsigned int InDet::InDetBeamSpotFinder::m_maxEventsPerFit

private

Definition at line 60 of file InDetBeamSpotFinder.h.

m_maxLBsPerFit

unsigned int InDet::InDetBeamSpotFinder::m_maxLBsPerFit

private

Definition at line 61 of file InDetBeamSpotFinder.h.

m_maxRunsPerFit

unsigned int InDet::InDetBeamSpotFinder::m_maxRunsPerFit

private

Definition at line 59 of file InDetBeamSpotFinder.h.

m_maxTrackNum

int InDet::InDetBeamSpotFinder::m_maxTrackNum {}

private

Definition at line 115 of file InDetBeamSpotFinder.h.

{};         //max Tracks in a vertex to be added

m_maxTransverseError

double InDet::InDetBeamSpotFinder::m_maxTransverseError {}

private

Definition at line 117 of file InDetBeamSpotFinder.h.

{}; // max transverse vertex resolution

m_minTrackNum

int InDet::InDetBeamSpotFinder::m_minTrackNum {}

private

Definition at line 114 of file InDetBeamSpotFinder.h.

{};         //min Tracks in a vertex to be added

m_minVertexNum

unsigned int InDet::InDetBeamSpotFinder::m_minVertexNum {}

private

Definition at line 120 of file InDetBeamSpotFinder.h.

{};        //min vertex count for solution

m_minVtxProb

double InDet::InDetBeamSpotFinder::m_minVtxProb {}

private

Definition at line 119 of file InDetBeamSpotFinder.h.

{}; // probability cut on chi2/ndf

m_root_bs

TTree* InDet::InDetBeamSpotFinder::m_root_bs =0

private

Definition at line 68 of file InDetBeamSpotFinder.h.

m_root_vrt

TTree* InDet::InDetBeamSpotFinder::m_root_vrt =0

private

Definition at line 69 of file InDetBeamSpotFinder.h.

m_root_vtx

vertexNtuple_struct InDet::InDetBeamSpotFinder::m_root_vtx

private

Definition at line 88 of file InDetBeamSpotFinder.h.

m_secondsPerFit

unsigned long InDet::InDetBeamSpotFinder::m_secondsPerFit {}

private

Definition at line 137 of file InDetBeamSpotFinder.h.

{}; 

m_sortedEventList

std::vector<std::vector<BeamSpot::Event> > InDet::InDetBeamSpotFinder::m_sortedEventList

private

Definition at line 94 of file InDetBeamSpotFinder.h.

m_sortedVertices

std::vector<std::vector<BeamSpot::VrtHolder> > InDet::InDetBeamSpotFinder::m_sortedVertices

private

Definition at line 95 of file InDetBeamSpotFinder.h.

m_thistSvc

ServiceHandle<ITHistSvc> InDet::InDetBeamSpotFinder::m_thistSvc {this, "THistSvc", "THistSvc"}

private

Definition at line 47 of file InDetBeamSpotFinder.h.

{this, "THistSvc", "THistSvc"};

m_toolSvc

ServiceHandle<IToolSvc> InDet::InDetBeamSpotFinder::m_toolSvc {this, "ToolSvc", "ToolSvc"}

private

Definition at line 46 of file InDetBeamSpotFinder.h.

{this, "ToolSvc", "ToolSvc"};

m_useFilledBCIDsOnly

bool InDet::InDetBeamSpotFinder::m_useFilledBCIDsOnly {}

private

Definition at line 130 of file InDetBeamSpotFinder.h.

{}; // Only use filled BCIDs not empty ones (data)

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vertexContainer

SG::ReadHandleKey<xAOD::VertexContainer> InDet::InDetBeamSpotFinder::m_vertexContainer {this, "VertexContainer", "PrimaryVertices", "Vertex container name"}

private

Definition at line 55 of file InDetBeamSpotFinder.h.

{this, "VertexContainer", "PrimaryVertices", "Vertex container name"};

m_vertexTreeName

std::string InDet::InDetBeamSpotFinder::m_vertexTreeName

private

Definition at line 91 of file InDetBeamSpotFinder.h.

m_vertexTypeNames

std::vector<std::string> InDet::InDetBeamSpotFinder::m_vertexTypeNames {}

private

Definition at line 122 of file InDetBeamSpotFinder.h.

{}; //names of vertexTypes

m_vertexTypes

std::vector<xAOD::VxType::VertexType> InDet::InDetBeamSpotFinder::m_vertexTypes {}

private

Definition at line 123 of file InDetBeamSpotFinder.h.

{};

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_writeAllVertices

bool InDet::InDetBeamSpotFinder::m_writeAllVertices {}

private

Definition at line 133 of file InDetBeamSpotFinder.h.

{};

m_writeVertexNtuple

bool InDet::InDetBeamSpotFinder::m_writeVertexNtuple

private

Definition at line 65 of file InDetBeamSpotFinder.h.

---

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

• InDetBeamSpotFinder.h
• InDetBeamSpotFinder.cxx