TBPlaneTrackingAlgo Class Reference

#include <TBPlaneTrackingAlgo.h>

Inheritance diagram for TBPlaneTrackingAlgo:

Collaboration diagram for TBPlaneTrackingAlgo:

Public Member Functions
	TBPlaneTrackingAlgo (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~TBPlaneTrackingAlgo ()
virtual StatusCode	initialize () override
virtual StatusCode	execute () override
virtual StatusCode	finalize () 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
void	FillRandomHit ()
StatusCode	getnewcalib ()
bool	fitPlane (const TBHitPlaneCont *hitPlaneCont, double &a1, double &a2, double &chi2, std::vector< double > &residual)
	Fit data to the function u = a1 + a2*w. More...
bool	fitHits (const std::vector< double > &u, const std::vector< double > &w, const std::vector< double > &eu, double &a1, double &a2)
	Fit data to the function u = a1 + a2*w. More...
double	getResidual (const double &u, const double &w, const double &a1, const double &a2)
	Calculates the residual-> r = (u_i - a1 - a2*w_i) More...
double	getChi2 (const std::vector< double > &v_u, const std::vector< double > &v_w, const std::vector< double > &v_eu, const double &a1, const double &a2)
	Calculates the chi2 += ((u_i - a1 - a2*w_i)/eu)^2. More...
StatusCode	buildHits ()
	Build HitPlaneCont from BPC. More...
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
bool	m_testAlgo
TBHitPlaneCont	m_hitPlaneCont_u
TBHitPlaneCont	m_hitPlaneCont_v
std::vector< float >	m_bpc_posZX
std::vector< float >	m_bpc_posZY
std::vector< float >	m_bpc_posX
std::vector< float >	m_bpc_posY
std::vector< float >	m_bpc_errposZX
std::vector< float >	m_bpc_errposZY
std::vector< float >	m_bpc_errposX
std::vector< float >	m_bpc_errposY
std::vector< float >	m_bpc_errmeasX
std::vector< float >	m_bpc_errmeasY
std::vector< std::string >	m_bpc_names
std::string	m_calib_filename
unsigned int	m_runnumber
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 20 of file TBPlaneTrackingAlgo.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

TBPlaneTrackingAlgo()

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

Definition at line 21 of file TBPlaneTrackingAlgo.cxx.

  : AthAlgorithm(name, pSvcLocator),
    m_runnumber(0)
{
  declareProperty("TestAlgo",m_testAlgo=false);
  declareProperty("CalibFileName",  m_calib_filename="H6HitCalib.txt");
  declareProperty("BPCnames",  m_bpc_names);
  declareProperty("BPCposX",  m_bpc_posX);
  declareProperty("BPCposY",  m_bpc_posY);
  declareProperty("BPCposZX",  m_bpc_posZX);
  declareProperty("BPCposZY",  m_bpc_posZY);
  declareProperty("BPCerrposX",  m_bpc_errposX);
  declareProperty("BPCerrposY",  m_bpc_errposY);
  declareProperty("BPCerrposZX",  m_bpc_errposZX);
  declareProperty("BPCerrposZY",  m_bpc_errposZY);
  declareProperty("BPCerrmeasX",  m_bpc_errmeasX);
  declareProperty("BPCerrmeasY",  m_bpc_errmeasY);
  
}

~TBPlaneTrackingAlgo()

virtual TBPlaneTrackingAlgo::~TBPlaneTrackingAlgo ( )

inlinevirtual

Definition at line 27 of file TBPlaneTrackingAlgo.h.

{};

Member Function Documentation

buildHits()

StatusCode TBPlaneTrackingAlgo::buildHits ( )

private

Build HitPlaneCont from BPC.

Definition at line 355 of file TBPlaneTrackingAlgo.cxx.

{
  ATH_MSG_DEBUG (" In buildHits() ");
  m_hitPlaneCont_u.clear();
  m_hitPlaneCont_v.clear();
 
  TBBPCCont * bpcCont = nullptr;
  ATH_CHECK( evtStore()->retrieve(bpcCont, "BPCCont") );
 
  // Loop over BPC
  for (const TBBPC* bpc : *bpcCont) {
    std::string name = bpc->getDetectorName();
    ATH_MSG_DEBUG (" Hits in BPC "<< name);
    // Find calibration index for this BPC
    unsigned int ind=0;
    while(ind<m_bpc_names.size()) 
    {
      if(name==m_bpc_names[ind]) break; 
      else ind++;
    }
    if(ind==m_bpc_names.size()){
      ATH_MSG_ERROR ("No calibrations for BPC" << name);
      continue;
    }
    
    ATH_MSG_DEBUG (" BPC number "<< ind);
   
    
    if(!(bpc->isXPosOverflow()||bpc->isXPulseOverflow())){
      TBHitPlane * hitu = new TBHitPlane();
      float w = m_bpc_posZX[ind];
      //      float u =m_bpc_posX[ind] + m_bpc_invX[ind]*bpc->getXPos();
      float u =m_bpc_posX[ind] + bpc->getXPos();
      ATH_MSG_DEBUG ("BPC" << ind << "X =" << u);
      hitu->setPosu(u);  hitu->setPosw(w);
      float uerr=sqrt(m_bpc_errposX[ind]*m_bpc_errposX[ind]+m_bpc_errmeasX[ind]*m_bpc_errmeasX[ind]);
      hitu->setErroru(uerr);  hitu->setErrorw(m_bpc_errposZX[ind]);
      hitu->setValidFlag(bpc->isXPosOverflow()&&bpc->isYPosOverflow());
      m_hitPlaneCont_u.push_back(hitu);
    }
 
    if(!(bpc->isYPosOverflow()||bpc->isYPulseOverflow())){
      TBHitPlane * hitv = new TBHitPlane();
      float w = m_bpc_posZY[ind];
      //      float v =m_bpc_posY[ind] + m_bpc_invY[ind]*bpc->getYPos();
      float v =m_bpc_posY[ind] + bpc->getYPos();
      ATH_MSG_DEBUG ("BPC" << ind << "Y =" << v);
      hitv->setPosu(v);  hitv->setPosw(w);
      float verr=sqrt(m_bpc_errposY[ind]*m_bpc_errposY[ind]+m_bpc_errmeasY[ind]*m_bpc_errmeasY[ind]);
      hitv->setErroru(verr);  hitv->setErrorw(m_bpc_errposZY[ind]);
      hitv->setValidFlag(bpc->isXPosOverflow()&&bpc->isYPosOverflow());
      m_hitPlaneCont_v.push_back(hitv);
    }
  }
 
  return StatusCode::SUCCESS;
}

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 TBPlaneTrackingAlgo::execute ( )

overridevirtual

Definition at line 49 of file TBPlaneTrackingAlgo.cxx.

{
  ATH_MSG_DEBUG ("Executing TBPlaneTracking algorithm");
  const EventContext& ctx = Gaudi::Hive::currentContext();
  
  if(m_testAlgo == true){
    ATH_MSG_WARNING ("TBPlaneTrackingAlgo: " <<
                     "Using FillRandomHit method to get data");
    FillRandomHit();
  }
 
  // Check if this is not a physics event and exit OK
  // retrieve Event Info
  const TBEventInfo* theEventInfo;
  StatusCode checkOut = evtStore()->retrieve(theEventInfo,"TBEventInfo");
  if ( checkOut.isFailure() )
    {
      ATH_MSG_ERROR ("cannot retrieve TBEventInfo from StoreGate");
      return StatusCode::FAILURE;
    }
  int evtType = theEventInfo->getEventType();
  ATH_MSG_DEBUG ("Event Type found " << evtType);
  if (evtType != 1)  return StatusCode::SUCCESS;
 
  StatusCode sc1;
 
  // Get run number and get new calib constants -----------------------------
  unsigned int thisrun=ctx.eventID().run_number();
 
  if(thisrun != m_runnumber)
    {
      m_runnumber= thisrun;
      ATH_CHECK(getnewcalib());
    }
  // ------------------------------------------------------------------------
 
 
 
  sc1 = buildHits();
  
//   sc1 = m_StoreGate->retrieve(m_hitPlaneCont_u,"HitPlaneCont_U");
//   sc2 = m_StoreGate->retrieve(m_hitPlaneCont_v,"HitPlaneCont_V");
  
  if(sc1.isFailure()){
    ATH_MSG_WARNING ("TBPlaneTrackingAlgo: Retrieval of HitPlane failed");
    setFilterPassed(false);
    return StatusCode::SUCCESS;
  } else {
 
    int hitPlaneNumU, hitPlaneNumV;
    hitPlaneNumU = m_hitPlaneCont_u.size();
    hitPlaneNumV = m_hitPlaneCont_v.size();
    
    if( (hitPlaneNumU < 2) || (hitPlaneNumV < 2) ){
      ATH_MSG_WARNING ("TBPlaneTrackingAlgo: "
                       << "Not enough hits in one or both planes, " 
                       << "Cannot make track.");
      
      setFilterPassed(false);
      return StatusCode::SUCCESS;
    }
 
    // Setting the chi2's for the track //
    double chi2_u = 0;
    double chi2_v = 0;
    std::vector<double> residual_u, residual_v;
    double a1_u = 0;
    double a2_u = 0;
    double a1_v = 0;
    double a2_v = 0;
    bool check = true;
 
    check = fitPlane(&m_hitPlaneCont_u, a1_u, a2_u, chi2_u, residual_u);
    if(check == false){
      ATH_MSG_ERROR ("TBPlaneTrackingAlgo: " << "Fit failure.");
      setFilterPassed(false);
      return StatusCode::SUCCESS;
    }
    
    check = fitPlane(&m_hitPlaneCont_v, a1_v, a2_v, chi2_v, residual_v);
    if(check == false){
      ATH_MSG_ERROR ("TBPlaneTrackingAlgo: " << "Fit failure.");
      setFilterPassed(false);
      return StatusCode::SUCCESS;
    } 
 
    // Setting the slopes & intercepts for each track //
    ATH_MSG_DEBUG ("Setting fit parameters of track.");
   
    TBTrack *track = new TBTrack(hitPlaneNumU, hitPlaneNumV);
    track->setUintercept(a1_u);
    track->setVintercept(a1_v);
    track->setUslope(a2_u);
    track->setVslope(a2_v);
 
    // Setting the residual for plane u //
    for(int i = 0; i < hitPlaneNumU; i++){
      track->setResidualu(i, residual_u[i]);
    }
    
    // Setting the residual for plane v //
    for(int i = 0; i < hitPlaneNumV; i++){
      track->setResidualv(i, residual_v[i]);
    }
    
    ATH_MSG_DEBUG ("chi2 in u: " << chi2_u);
    ATH_MSG_DEBUG ("chi2 in v: " << chi2_v);
    ATH_MSG_DEBUG ("Setting chi2s of track.");
 
    track->setChi2_u(chi2_u);
    track->setChi2_v(chi2_v);
   
    // Record track info. into storeGate //
    StatusCode sc = evtStore()->record(track,"Track");
    
    if ( sc.isFailure( ) ) {
      ATH_MSG_FATAL ("Cannot record Track");
    }
  }
  
  setFilterPassed(true);
  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();
}

FillRandomHit()

void TBPlaneTrackingAlgo::FillRandomHit ( )

private

Definition at line 314 of file TBPlaneTrackingAlgo.cxx.

{
 
  ATH_MSG_DEBUG ("Starting FillRandom");
 
  TBHitPlaneCont * hitPlaneCont_U = new TBHitPlaneCont(); 
  TBHitPlaneCont * hitPlaneCont_V = new TBHitPlaneCont();
  
  TBHitPlane *hitPaneU1 = new TBHitPlane();
  hitPaneU1->setPosu(1);  hitPaneU1->setPosw(1);
  hitPlaneCont_U->push_back(hitPaneU1);
  
 
  TBHitPlane *hitPaneU2 = new TBHitPlane();
  hitPaneU2->setPosu(2);  hitPaneU2->setPosw(2);
  hitPlaneCont_U->push_back(hitPaneU2);
 
  StatusCode sc1 = evtStore()->record(hitPlaneCont_U,"HitPlaneCont_U");
  if ( sc1.isFailure( ) ) {
    ATH_MSG_FATAL ("Cannot record HitPlaneCont_U");
  }
 
  TBHitPlane *hitPaneV1 = new TBHitPlane();
  hitPaneV1->setPosu(1);  hitPaneV1->setPosw(1);
  hitPlaneCont_V->push_back(hitPaneV1);
 
  TBHitPlane *hitPaneV2 = new TBHitPlane();
  hitPaneV2->setPosu(2);  hitPaneV2->setPosw(2);
  hitPlaneCont_V->push_back(hitPaneV2);
 
  StatusCode sc2 = evtStore()->record(hitPlaneCont_V,"HitPlaneCont_V");
  if ( sc2.isFailure( ) ) {
    ATH_MSG_FATAL ("Cannot record HitPlaneCont_V");
  }
 
  ATH_MSG_DEBUG ("FillRandom Done");
 
}

finalize()

StatusCode TBPlaneTrackingAlgo::finalize ( )

overridevirtual

Definition at line 175 of file TBPlaneTrackingAlgo.cxx.

{
  ATH_MSG_INFO ("Finalizing TBTracking algorithm");
  return StatusCode::SUCCESS; 
}

fitHits()

bool TBPlaneTrackingAlgo::fitHits ( const std::vector< double > &  u, const std::vector< double > &  w, const std::vector< double > &  eu, double &  a1, double &  a2  )

private

Fit data to the function u = a1 + a2*w.

Definition at line 236 of file TBPlaneTrackingAlgo.cxx.

{
  // v_u = vector of u data points
  // v_w = vector of w data points
  // v_eu = vector of errors in u data points
  // v_ew = vector of errors in w data points
  // a1 and a2 = fit parameters: u = a1 + a2*w
  
  // *** note that the fit algorithm used  (given in 'Numerical Methods 
  // in C' section 15.2) assumes that the w data points are known exactly 
  // i.e. that v_ew[i]=0 for all i
 
  // 'Numerical Methods' notes that the task of fitting a straight 
  // line model with errors in both coordinates is "considerably harder"
  // (section 15.3) - but clearly it could be done 
  
  double m_s = 0;
  double m_su = 0;
  double m_sww = 0;
  double m_sw = 0;
  double m_suw = 0;
 
  int numHits = v_u.size();
  for(int i = 0; i < numHits; i++){
    m_s += 1 / (v_eu[i]*v_eu[i]);
    m_su += v_u[i] / (v_eu[i]*v_eu[i]); 
    m_sww += v_w[i]*v_w[i] / (v_eu[i]*v_eu[i]);
    m_sw += v_w[i] / (v_eu[i]*v_eu[i]);
    m_suw += v_u[i]*v_w[i] / (v_eu[i]*v_eu[i]);
  }
 
  const double denom = (m_s*m_sww-m_sw*m_sw);
  if(denom == 0){
    return false;
  }
 
  const double inv_denom = 1. / denom;
  a1 = (m_su*m_sww - m_sw*m_suw) * inv_denom;
  a2 = (m_s*m_suw - m_su*m_sw) * inv_denom;
  return true;
}

fitPlane()

bool TBPlaneTrackingAlgo::fitPlane ( const TBHitPlaneCont *  hitPlaneCont, double &  a1, double &  a2, double &  chi2, std::vector< double > &  residual  )

private

Fit data to the function u = a1 + a2*w.

and determines intercept, slope, residual for each BPC, and chi2 on fit

Definition at line 183 of file TBPlaneTrackingAlgo.cxx.

{
  // Get number of hits in plane //
  int hitPlaneNum = hitPlaneCont->size();
  
  ATH_MSG_DEBUG ("The hit plane container size is: " << hitPlaneNum);
  
  std::vector<double> vec_u;
  std::vector<double> vec_w;
  std::vector<double> err_vec_u;
  std::vector<double> err_vec_w;
 
  for (const TBHitPlane* hp : *hitPlaneCont) {
    ATH_MSG_DEBUG ("Position in u: " << hp->getPosu());
    ATH_MSG_DEBUG ("Position in w: " << hp->getPosw());
    ATH_MSG_DEBUG ("Error in u: " << hp->getErroru());
    ATH_MSG_DEBUG ("Error in w: " << hp->getErrorw());
    vec_u.push_back(hp->getPosu());
    vec_w.push_back(hp->getPosw());
    err_vec_u.push_back(hp->getErroru());
    err_vec_w.push_back(hp->getErrorw());
  }
  
  if(vec_u.size() != vec_w.size()){
    ATH_MSG_ERROR ("TBPlaneTrackingAlgo: Invalid hits");
    return false;
  }
 
  bool check;
//  check = fitHits(vec_u, vec_w, err_vec_u, err_vec_w, a1, a2);
  check = fitHits(vec_u, vec_w, err_vec_u, a1, a2);
  if(check == false){
    ATH_MSG_ERROR ("TBPlaneTrackingAlgo: Invalid denumerator");
    return false;
  }
 
//  chi2 = getChi2(vec_u, vec_w, err_vec_u, err_vec_w, a1, a2);
  chi2 = getChi2(vec_u, vec_w, err_vec_u, a1, a2);
 
  ATH_MSG_DEBUG ("Fit results:" 
                 << " intercept = " << a1 
                 << " and slope = " << a2);
  
  // Must be wrong //
  for (int i = 0; i < hitPlaneNum; i++) {
    residual.push_back(getResidual(vec_u[i],vec_w[i],a1, a2));
  }
  return true;
}

getChi2()

double TBPlaneTrackingAlgo::getChi2 ( const std::vector< double > &  v_u, const std::vector< double > &  v_w, const std::vector< double > &  v_eu, const double &  a1, const double &  a2  )

private

Calculates the chi2 += ((u_i - a1 - a2*w_i)/eu)^2.

Definition at line 291 of file TBPlaneTrackingAlgo.cxx.

{ 
  // *** as in fitHits() above, the chi squared calculation below
  // (taken from 'Numerical Methods in C' section 15.2) 
  // assumes that the w data points are known exactly i.e. that 
  // v_ew[i]=0 for all i
  
  int numX = v_u.size();
  double chi2 = 0;
  
  for(int i = 0; i < numX; i++){
    chi2 += (v_u[i] - a1 - a2*v_w[i])*(v_u[i] - a1 - a2*v_w[i])/(v_eu[i]*v_eu[i]);
  }  
  return chi2;
}

getnewcalib()

StatusCode TBPlaneTrackingAlgo::getnewcalib ( )

private

Definition at line 415 of file TBPlaneTrackingAlgo.cxx.

{
  // Get calib constant from an ASCII file with the following structure :
  // 
  // runnumber
  // bpcnumber1 coeff1 coeff2 ... coeff8
  // bpcnumber2 coeff1 coeff2 ... coeff8
  // ...
  // bpcnumber6 coeff1 coeff2 ... coeff8
  // runnumber
  // ...
  //
  // coeff must have the following order :
  // bpcnumber posX posY posZX posZY errposX errposY errposZX errposZY errmeasX errmeasY
 
  ATH_MSG_DEBUG ("Get new calibs for run " << m_runnumber);
  
  int bpcnumber= m_bpc_names.size();
 
  m_bpc_posX.clear(); m_bpc_posX.resize(bpcnumber);
  m_bpc_posY.clear(); m_bpc_posY.resize(bpcnumber);
  m_bpc_posZX.clear(); m_bpc_posZX.resize(bpcnumber);
  m_bpc_posZY.clear(); m_bpc_posZY.resize(bpcnumber);
  m_bpc_errposX.clear(); m_bpc_errposX.resize(bpcnumber);
  m_bpc_errposY.clear(); m_bpc_errposY.resize(bpcnumber);
  m_bpc_errposZX.clear(); m_bpc_errposZX.resize(bpcnumber);
  m_bpc_errposZY.clear(); m_bpc_errposZY.resize(bpcnumber);
  m_bpc_errmeasX.clear(); m_bpc_errmeasX.resize(bpcnumber);
  m_bpc_errmeasY.clear(); m_bpc_errmeasY.resize(bpcnumber);
 
  int pos;
 
  std::ifstream calibfile;
  std::string filename = PathResolver::find_file (m_calib_filename, "DATAPATH");
  calibfile.open(filename.c_str());
  if(!calibfile.good()){
    ATH_MSG_WARNING (" Problem with file named "<< m_calib_filename << " in $DATAPATH");
    return StatusCode::FAILURE;
  } else {
    ATH_MSG_DEBUG (" File: "<< filename << " opened");
  }
  unsigned int runnumber;
  calibfile >> runnumber;
  pos = calibfile.tellg();
  ATH_MSG_DEBUG (" Run number "<< runnumber);
  while((runnumber<m_runnumber)&&(!calibfile.eof()))
    {
      runnumber=0;
      pos = calibfile.tellg();
      // discard next lines
      for(int j=0;j<bpcnumber+1;j++) calibfile.ignore(5000,'\n');
      // check next runnumber
      calibfile >> runnumber;
      if(runnumber==0) {  // reached an empty line : exit.
        ATH_MSG_DEBUG ("empty line");
        calibfile.clear();
        break;
      }
      ATH_MSG_DEBUG (" Run number "<< runnumber);
    }
  
  // Now we found the good set of constant (the ones following pos)
  if(runnumber==m_runnumber)  pos = calibfile.tellg();
  ATH_MSG_DEBUG (" Pos = "<< pos);
  calibfile.seekg(pos);
  ATH_MSG_DEBUG (" Will use the following constants :");
  for(int j=0;j<bpcnumber;j++) 
    {
      int bpcn;
      calibfile >> bpcn; 
      calibfile >> m_bpc_posX[j];
      calibfile >> m_bpc_posY[j];
      calibfile >> m_bpc_posZX[j];
      calibfile >> m_bpc_posZY[j];
      calibfile >> m_bpc_errposX[j];
      calibfile >> m_bpc_errposY[j];
      calibfile >> m_bpc_errposZX[j];
      calibfile >> m_bpc_errposZY[j];
      calibfile >> m_bpc_errmeasX[j];
      calibfile >> m_bpc_errmeasY[j];
 
      ATH_MSG_DEBUG (bpcn << " "<<m_bpc_posX[j]
                     << " "<< m_bpc_posY[j]    
                     << " "<< m_bpc_posZX[j]   
                     << " "<< m_bpc_posZY[j]   
                     << " "<< m_bpc_errposX[j] 
                     << " "<< m_bpc_errposY[j] 
                     << " "<< m_bpc_errposZX[j]  
                     << " "<< m_bpc_errposZY[j]
                     << " "<< m_bpc_errmeasX[j]
                     << " "<< m_bpc_errmeasY[j]);
    }
  
  calibfile.close();
 
  return StatusCode::SUCCESS;  
  
}

getResidual()

double TBPlaneTrackingAlgo::getResidual ( const double &  u, const double &  w, const double &  a1, const double &  a2  )

private

Calculates the residual-> r = (u_i - a1 - a2*w_i)

Definition at line 283 of file TBPlaneTrackingAlgo.cxx.

{
  return (u - a1 - a2*w);
}

initialize()

StatusCode TBPlaneTrackingAlgo::initialize ( )

overridevirtual

Definition at line 42 of file TBPlaneTrackingAlgo.cxx.

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

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_bpc_errmeasX

std::vector<float> TBPlaneTrackingAlgo::m_bpc_errmeasX

private

Definition at line 83 of file TBPlaneTrackingAlgo.h.

m_bpc_errmeasY

std::vector<float> TBPlaneTrackingAlgo::m_bpc_errmeasY

private

Definition at line 83 of file TBPlaneTrackingAlgo.h.

m_bpc_errposX

std::vector<float> TBPlaneTrackingAlgo::m_bpc_errposX

private

Definition at line 82 of file TBPlaneTrackingAlgo.h.

m_bpc_errposY

std::vector<float> TBPlaneTrackingAlgo::m_bpc_errposY

private

Definition at line 82 of file TBPlaneTrackingAlgo.h.

m_bpc_errposZX

std::vector<float> TBPlaneTrackingAlgo::m_bpc_errposZX

private

Definition at line 82 of file TBPlaneTrackingAlgo.h.

m_bpc_errposZY

std::vector<float> TBPlaneTrackingAlgo::m_bpc_errposZY

private

Definition at line 82 of file TBPlaneTrackingAlgo.h.

m_bpc_names

std::vector<std::string> TBPlaneTrackingAlgo::m_bpc_names

private

Definition at line 85 of file TBPlaneTrackingAlgo.h.

m_bpc_posX

std::vector<float> TBPlaneTrackingAlgo::m_bpc_posX

private

Definition at line 81 of file TBPlaneTrackingAlgo.h.

m_bpc_posY

std::vector<float> TBPlaneTrackingAlgo::m_bpc_posY

private

Definition at line 81 of file TBPlaneTrackingAlgo.h.

m_bpc_posZX

std::vector<float> TBPlaneTrackingAlgo::m_bpc_posZX

private

Definition at line 81 of file TBPlaneTrackingAlgo.h.

m_bpc_posZY

std::vector<float> TBPlaneTrackingAlgo::m_bpc_posZY

private

Definition at line 81 of file TBPlaneTrackingAlgo.h.

m_calib_filename

std::string TBPlaneTrackingAlgo::m_calib_filename

private

Definition at line 88 of file TBPlaneTrackingAlgo.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_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_hitPlaneCont_u

TBHitPlaneCont TBPlaneTrackingAlgo::m_hitPlaneCont_u

private

Definition at line 77 of file TBPlaneTrackingAlgo.h.

m_hitPlaneCont_v

TBHitPlaneCont TBPlaneTrackingAlgo::m_hitPlaneCont_v

private

Definition at line 78 of file TBPlaneTrackingAlgo.h.

m_runnumber

unsigned int TBPlaneTrackingAlgo::m_runnumber

private

Definition at line 89 of file TBPlaneTrackingAlgo.h.

m_testAlgo

bool TBPlaneTrackingAlgo::m_testAlgo

private

Definition at line 76 of file TBPlaneTrackingAlgo.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.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:

• TBPlaneTrackingAlgo.h
• TBPlaneTrackingAlgo.cxx