ALFA_LocRec Class Reference

#include <ALFA_LocRec.h>

Inheritance diagram for ALFA_LocRec:

Public Member Functions
	ALFA_LocRec (const std::string &name, ISvcLocator *pSvcLocator)
	~ALFA_LocRec ()
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
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	ReadGeometryDetCS ()
bool	StoreReconstructionGeometry (const eRPotName eRPName, const eFiberType eFType, const char *szDataDestination)
void	SaveGeometry ()
void	ClearGeometry ()
StatusCode	ALFACollectionReading (std::list< MDHIT > &ListMDHits, std::list< ODHIT > &ListODHits)
StatusCode	RecordCollection ()
StatusCode	RecordODCollection ()
StatusCode	ExecuteRecoMethod (const std::string &strAlgo, const eRPotName eRPName, const std::list< MDHIT > &ListMDHits, const std::list< ODHIT > &ListODHits)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
GEOMETRYCONFIGURATION	m_Config
ALFA_GeometryReader *	m_pGeometryReader
ALFA_LocRecEvCollection *	m_pLocRecEvCollection
ALFA_LocRecEvent *	m_pLocRecEvent
ALFA_LocRecODEvCollection *	m_pLocRecODEvCollection
ALFA_LocRecODEvent *	m_pLocRecODEvent
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfoKey {this, "EvtInfo", "EventInfo", "EventInfo name"}
std::list< eRPotName >	m_ListExistingRPots
UInt_t	m_eventNum
Int_t	m_iDataType
Int_t	m_iEvent
Int_t	m_iMultiplicityCutMD
Int_t	m_iNumLayerCutMD
Int_t	m_iUVCutMD
Int_t	m_iUVCutMDHalfReco
Int_t	m_iMultiplicityCutOD
Float_t	m_fOverlapCutMD
Float_t	m_fDistanceCutOD
Float_t	m_iLayerCutOD
Bool_t	m_bEdgeMethod_Opt_Sisters
Bool_t	m_bEdgeMethod_Opt_UseGaps
Float_t	m_faMD [RPOTSCNT][ALFALAYERSCNT *ALFAPLATESCNT][ALFAFIBERSCNT]
Float_t	m_fbMD [RPOTSCNT][ALFALAYERSCNT *ALFAPLATESCNT][ALFAFIBERSCNT]
Float_t	m_fzMD [RPOTSCNT][ALFALAYERSCNT *ALFAPLATESCNT][ALFAFIBERSCNT]
Float_t	m_faOD [RPOTSCNT][ODPLATESCNT][ODSIDESCNT][ODLAYERSCNT *ODFIBERSCNT]
Float_t	m_fbOD [RPOTSCNT][ODPLATESCNT][ODSIDESCNT][ODLAYERSCNT *ODFIBERSCNT]
Float_t	m_fzOD [RPOTSCNT][ODPLATESCNT][ODSIDESCNT][ODLAYERSCNT *ODFIBERSCNT]
std::string	m_strKeyGeometryForReco
std::vector< std::string >	m_vecListAlgoMD
std::vector< std::string >	m_vecListAlgoOD
std::string	m_strAlgoOD
std::string	m_strAlgoMD
std::string	m_strKeyLocRecEvCollection
std::string	m_strKeyLocRecODEvCollection
std::string	m_strCollectionName
std::string	m_strODCollectionName
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 43 of file ALFA_LocRec.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

ALFA_LocRec()

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

Definition at line 24 of file ALFA_LocRec.cxx.

                                                                        :
    AthAlgorithm(name, pSvcLocator)
{
    ATH_MSG_DEBUG("begin ALFA_LocRec::ALFA_LocRec");
 
    m_pGeometryReader = new ALFA_GeometryReader();
 
//  m_Config.eOwner=EOT_RECO;
    m_Config.clear();
 
    //begin of Geometry properties
    declareProperty("MetrologyType",m_Config.eRPMetrologyGeoType=EMT_NOMINAL);//EGST_FILE
    declareProperty("MetrologySource",m_Config.strRPMetrologyConnString=std::string(""));//"rpmetrology.dat"
    declareProperty("ShiftToX97Pos",m_Config.bShiftToX97Pos=false);
 
    declareProperty("B7L1U_PosType",m_Config.CfgRPosParams[0].eRPPosType=ERPPT_ACTIVE);
    declareProperty("B7L1U_MDGeometryType",m_Config.CfgRPosParams[0].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7L1U_ODGeometryType",m_Config.CfgRPosParams[0].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7L1U_MDSource",m_Config.CfgRPosParams[0].strMDConnString=std::string(""));
    declareProperty("B7L1U_ODSource",m_Config.CfgRPosParams[0].strODConnString=std::string(""));
//  declareProperty("B7L1U_CurrentLVDT",m_Config.CfgRPosParams[0].fCurrentLVDTmm=7.0*CLHEP::mm);
 
    declareProperty("B7L1L_PosType",m_Config.CfgRPosParams[1].eRPPosType=ERPPT_ACTIVE);
    declareProperty("B7L1L_MDGeometryType",m_Config.CfgRPosParams[1].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7L1L_ODGeometryType",m_Config.CfgRPosParams[1].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7L1L_MDSource",m_Config.CfgRPosParams[1].strMDConnString=std::string(""));
    declareProperty("B7L1L_ODSource",m_Config.CfgRPosParams[1].strODConnString=std::string(""));
//  declareProperty("B7L1L_CurrentLVDT",m_Config.CfgRPosParams[1].fCurrentLVDTmm=7.0*CLHEP::mm);
 
    declareProperty("A7L1U_PosType",m_Config.CfgRPosParams[2].eRPPosType=ERPPT_ACTIVE);
    declareProperty("A7L1U_MDGeometryType",m_Config.CfgRPosParams[2].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7L1U_ODGeometryType",m_Config.CfgRPosParams[2].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7L1U_MDSource",m_Config.CfgRPosParams[2].strMDConnString=std::string(""));
    declareProperty("A7L1U_ODSource",m_Config.CfgRPosParams[2].strODConnString=std::string(""));
//  declareProperty("A7L1U_CurrentLVDT",m_Config.CfgRPosParams[2].fCurrentLVDTmm=7.0*CLHEP::mm);
 
    declareProperty("A7L1L_PosType",m_Config.CfgRPosParams[3].eRPPosType=ERPPT_ACTIVE);
    declareProperty("A7L1L_MDGeometryType",m_Config.CfgRPosParams[3].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7L1L_ODGeometryType",m_Config.CfgRPosParams[3].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7L1L_MDSource",m_Config.CfgRPosParams[3].strMDConnString=std::string(""));
    declareProperty("A7L1L_ODSource",m_Config.CfgRPosParams[3].strODConnString=std::string(""));
//  declareProperty("A7L1L_CurrentLVDT",m_Config.CfgRPosParams[3].fCurrentLVDTmm=7.0*CLHEP::mm);
 
    declareProperty("A7R1U_PosType",m_Config.CfgRPosParams[4].eRPPosType=ERPPT_ACTIVE);
    declareProperty("A7R1U_MDGeometryType",m_Config.CfgRPosParams[4].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7R1U_ODGeometryType",m_Config.CfgRPosParams[4].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7R1U_MDSource",m_Config.CfgRPosParams[4].strMDConnString=std::string(""));
    declareProperty("A7R1U_ODSource",m_Config.CfgRPosParams[4].strODConnString=std::string(""));
//  declareProperty("A7R1U_CurrentLVDT",m_Config.CfgRPosParams[4].fCurrentLVDTmm=7.0*CLHEP::mm);
 
    declareProperty("A7R1L_PosType",m_Config.CfgRPosParams[5].eRPPosType=ERPPT_ACTIVE);
    declareProperty("A7R1L_MDGeometryType",m_Config.CfgRPosParams[5].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7R1L_ODGeometryType",m_Config.CfgRPosParams[5].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("A7R1L_MDSource",m_Config.CfgRPosParams[5].strMDConnString=std::string(""));
    declareProperty("A7R1L_ODSource",m_Config.CfgRPosParams[5].strODConnString=std::string(""));
//  declareProperty("A7L1U_CurrentLVDT",m_Config.CfgRPosParams[5].fCurrentLVDTmm=7.0*CLHEP::mm);
 
    declareProperty("B7R1U_PosType",m_Config.CfgRPosParams[6].eRPPosType=ERPPT_ACTIVE);
    declareProperty("B7R1U_MDGeometryType",m_Config.CfgRPosParams[6].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7R1U_ODGeometryType",m_Config.CfgRPosParams[6].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7R1U_MDSource",m_Config.CfgRPosParams[6].strMDConnString=std::string(""));
    declareProperty("B7R1U_ODSource",m_Config.CfgRPosParams[6].strODConnString=std::string(""));
//  declareProperty("B7R1U_CurrentLVDT",m_Config.CfgRPosParams[6].fCurrentLVDTmm=7.0*CLHEP::mm);
 
    declareProperty("B7R1L_PosType",m_Config.CfgRPosParams[7].eRPPosType=ERPPT_ACTIVE);
    declareProperty("B7R1L_MDGeometryType",m_Config.CfgRPosParams[7].eMDGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7R1L_ODGeometryType",m_Config.CfgRPosParams[7].eODGeoType=EGST_IDEALGEOMETRY);
    declareProperty("B7R1L_MDSource",m_Config.CfgRPosParams[7].strMDConnString=std::string(""));
    declareProperty("B7R1L_ODSource",m_Config.CfgRPosParams[7].strODConnString=std::string(""));
//  declareProperty("B7R1L_CurrentLVDT",m_Config.CfgRPosParams[7].fCurrentLVDTmm=7.0*CLHEP::mm);
    //end of Geometry properties
 
    // data type using in the local reconstruction
    // for the simulation data the value is 0, for the real data the value is 1. Unset value is -1
    declareProperty("DataType", m_iDataType=1, "data type using in the local reconstruction");
 
    //reconstruction methods properties - MDTracking, MDMultiple, ODTracking
    declareProperty("MultiplicityCutOD", m_iMultiplicityCutOD=30);
    declareProperty("DistanceCutOD", m_fDistanceCutOD=0.5);
    declareProperty("LayerCutOD", m_iLayerCutOD=3);
 
    declareProperty("OverlapCutMD", m_fOverlapCutMD=0.5);
    declareProperty("MultiplicityCutMD", m_iMultiplicityCutMD=5);
    declareProperty("NumLayerCutMD", m_iNumLayerCutMD=3);
    declareProperty("UVCutMD", m_iUVCutMD=3);
    declareProperty("UVCutMDHalfReco", m_iUVCutMDHalfReco=3);
 
    //reconstruction methods properties - EdgeMethod
    declareProperty("EdgeMethod_Opt_Sisters", m_bEdgeMethod_Opt_Sisters = false);
    declareProperty("EdgeMethod_Opt_UseGaps", m_bEdgeMethod_Opt_UseGaps = false);
 
    //reconstruction method selection for OD & MD
    declareProperty("AlgoOD", m_strAlgoOD="ODTracking");
    declareProperty("AlgoMD", m_strAlgoMD="MDTracking");
 
    m_vecListAlgoMD.clear();
    declareProperty("ListAlgoMD", m_vecListAlgoMD);
    m_vecListAlgoOD.clear();
    declareProperty("ListAlgoOD", m_vecListAlgoOD);
 
    m_ListExistingRPots.clear();
 
    m_strKeyGeometryForReco         = "ALFA_GeometryForReco";
    m_strCollectionName             = "ALFA_DigitCollection";
    m_strODCollectionName           = "ALFA_ODDigitCollection";
    m_strKeyLocRecEvCollection      = "ALFA_LocRecEvCollection";
    m_strKeyLocRecODEvCollection    = "ALFA_LocRecODEvCollection";
 
    m_pLocRecEvCollection = NULL;
    m_pLocRecEvent = NULL;
    m_pLocRecODEvCollection = NULL;
    m_pLocRecODEvent = NULL;
 
    m_eventNum = 0;
    m_iEvent   = 0;
 
    ATH_MSG_DEBUG("end ALFA_LocRec::ALFA_LocRec");
}

~ALFA_LocRec()

ALFA_LocRec::~ALFA_LocRec

(

)

Definition at line 143 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::~ALFA_LocRec");
 
//  if(m_pGeometryReader!=NULL)
//  {
//      delete m_pGeometryReader;
//      m_pGeometryReader = NULL;
//  }
 
    ATH_MSG_DEBUG("begin ALFA_LocRec::~ALFA_LocRec");
}

Member Function Documentation

ALFACollectionReading()

StatusCode ALFA_LocRec::ALFACollectionReading ( std::list< MDHIT > & ListMDHits, std::list< ODHIT > & ListODHits )

private

Definition at line 361 of file ALFA_LocRec.cxx.

{
    StatusCode sc;
 
    ATH_MSG_DEBUG("begin ALFA_LocRec::ALFACollectionReading()");
 
    ODHIT ODHit;
    MDHIT MDHit;
 
    ListMDHits.clear();
    ListODHits.clear();
 
    const ALFA_DigitCollection* mcGen = 0;
    const ALFA_ODDigitCollection* mcODGen = 0;
 
    sc = evtStore()->retrieve(mcGen,m_strCollectionName);
    if(sc.isFailure() || !mcGen)
    {
        return StatusCode::FAILURE;
    }
 
    ALFA_DigitCollection::const_iterator mcGenBeg = mcGen->begin();
    ALFA_DigitCollection::const_iterator mcGenEnd = mcGen->end();
 
    for(;mcGenBeg!=mcGenEnd;++mcGenBeg)
    {
        MDHit.iEvent = m_eventNum;
 
        MDHit.iRPot  = ((*mcGenBeg)->getStation());
        MDHit.iPlate = ((*mcGenBeg)->getPlate());
        MDHit.iFiber = ((*mcGenBeg)->getFiber());
 
        ListMDHits.push_back(MDHit);
    }
 
    sc = evtStore()->retrieve(mcODGen,m_strODCollectionName);
    if(sc.isFailure() || !mcODGen)
    {
        return StatusCode::FAILURE;
    }
 
    ALFA_ODDigitCollection::const_iterator mcODGenBeg = mcODGen->begin();
    ALFA_ODDigitCollection::const_iterator mcODGenEnd = mcODGen->end();
 
    for(;mcODGenBeg!=mcODGenEnd;++mcODGenBeg)
    {
        ODHit.iEvent = m_eventNum;
        ODHit.iRPot  = ((*mcODGenBeg)->getStation());
        ODHit.iPlate = ((*mcODGenBeg)->getPlate());
        ODHit.iSide  = ((*mcODGenBeg)->getSide());
        ODHit.iFiber = ((*mcODGenBeg)->getFiber());
 
        ListODHits.push_back(ODHit);
    }
 
    ATH_MSG_DEBUG("end ALFA_LocRec::ALFACollectionReading()");
 
    return StatusCode::SUCCESS;
}

ClearGeometry()

void ALFA_LocRec::ClearGeometry ( )

private

Definition at line 347 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::ClearGeometry()");
 
    memset(&m_faMD, 0, sizeof(m_faMD));
    memset(&m_fbMD, 0, sizeof(m_fbMD));
    memset(&m_fzMD, 0, sizeof(m_fzMD));
    memset(&m_faOD, 0, sizeof(m_faOD));
    memset(&m_fbOD, 0, sizeof(m_fbOD));
    memset(&m_fzOD, 0, sizeof(m_fzOD));
 
    ATH_MSG_DEBUG("end ALFA_LocRec::ClearGeometry()");
}

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 ALFA_LocRec::execute

(

)

Definition at line 211 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::execute()");
 
    StatusCode sc;
    std::list<eRPotName>::const_iterator iterRPName;
    eRPotName eRPName;
 
    std::list<MDHIT> ListMDHits;
    std::list<ODHIT> ListODHits;
    ListMDHits.clear();
    ListODHits.clear();
 
    m_eventNum = 0;
    SG::ReadHandle<xAOD::EventInfo> eventInfo (m_eventInfoKey,getContext());
    if(!eventInfo.isValid()) {
        ATH_MSG_ERROR("ALFA_LocRec, Cannot get event info.");
//      return sc;
    }
    else
    {
        // current event number
        m_eventNum = eventInfo->eventNumber();
    }
 
    //************************************************************
//  if (m_eventNum!=73681035) return StatusCode::SUCCESS;
    //************************************************************
 
    sc = RecordCollection();
    if (sc.isFailure())
    {
        ATH_MSG_ERROR("RecordCollection() failed");
        return StatusCode::SUCCESS;
    }
 
    sc = RecordODCollection();
    if (sc.isFailure())
    {
        ATH_MSG_ERROR("RecordODCollection() failed");
        return StatusCode::SUCCESS;
    }
 
    sc = ALFACollectionReading(ListMDHits, ListODHits);
    if (sc.isFailure())
    {
        ATH_MSG_WARNING("ALFA_Collection_Reading Failed");
//      return StatusCode::SUCCESS;
    }
    else
    {
        std::string strAlgoMD;
        for(auto & i : m_vecListAlgoMD)
        {
            strAlgoMD = i;
            for(iterRPName=m_ListExistingRPots.begin();iterRPName!=m_ListExistingRPots.end();++iterRPName)
            {
                eRPName = *iterRPName;
 
//              cout << "eRPName = " << eRPName << endl;
                //execute MD algorithm
                sc = ExecuteRecoMethod(strAlgoMD, eRPName, ListMDHits, ListODHits);
                if (sc.isFailure())
                {
                    ATH_MSG_FATAL("Algorithm " << strAlgoMD << " for the pot " << m_pGeometryReader->GetRPotLabel(eRPName) << " failure!");
                    return sc;
                }
            }
        }
 
        std::string strAlgoOD;
        for(auto & i : m_vecListAlgoOD)
        {
            strAlgoOD = i;
            for(iterRPName=m_ListExistingRPots.begin();iterRPName!=m_ListExistingRPots.end();++iterRPName)
            {
                eRPName = *iterRPName;
 
    //          cout << "eRPName = " << eRPName << endl;
                //execute OD algorithm
                sc = ExecuteRecoMethod(strAlgoOD, eRPName, ListMDHits, ListODHits);
                if (sc.isFailure())
                {
                    ATH_MSG_FATAL("Algorithm " << strAlgoOD << " for the pot " << m_pGeometryReader->GetRPotLabel(eRPName) << " failure!");
                    return sc;
                }
            }
        }
 
 
 
 
//      // local reconstruction process for all existing pots, MD & OD
//      for(iterRPName=m_ListExistingRPots.begin();iterRPName!=m_ListExistingRPots.end();iterRPName++)
//      {
//          eRPName = *iterRPName;
    
//          //execute MD reconstruction method
//          if (!m_strAlgoMD.empty())
//          {
//              sc = ExecuteRecoMethod(m_strAlgoMD, eRPName, ListMDHits, ListODHits);
//              if (sc.isFailure())
//              {
//                  ATH_MSG_ERROR("Algorithm " << m_strAlgoMD << " for pot " << m_pGeometryReader->GetRPotLabel(eRPName) << " failure!");
//  //              return sc;
//              }
//          }
    
//          //execute OD reconstruction method
//          if (!m_strAlgoOD.empty())
//          {
//              sc = ExecuteRecoMethod(m_strAlgoOD, eRPName, ListMDHits, ListODHits);
//              if (sc.isFailure())
//              {
//                  ATH_MSG_ERROR("Algorithm " << m_strAlgoOD << " for pot " << m_pGeometryReader->GetRPotLabel(eRPName) << " failure!");
//  //              return sc;
//              }
//          }
//      }
    }
 
    m_iEvent++;
    ATH_MSG_DEBUG("end ALFA_LocRec::execute()");
    return StatusCode::SUCCESS;
}

ExecuteRecoMethod()

StatusCode ALFA_LocRec::ExecuteRecoMethod ( const std::string & strAlgo, const eRPotName eRPName, const std::list< MDHIT > & ListMDHits, const std::list< ODHIT > & ListODHits )

private

Definition at line 677 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::ExecuteRecoMethod()");
    if (eRPName < 1){
        ATH_MSG_ERROR("ALFA_LocRec::ExecuteRecoMethod called with invalid detector enum");
        return StatusCode::FAILURE;
    }
    StatusCode sc = StatusCode::SUCCESS;
    ODRESULT ODResults;
    std::list<ODRESULT> listODResults;
    listODResults.clear();
    Float_t fRecPosX = -9999.0;
    Float_t fRecPosY = -9999.0;
 
    std::map<string, int> mapRecoMethods;
    mapRecoMethods.clear();
    mapRecoMethods.insert(std::pair<string, int>("ODTracking", 1));
    mapRecoMethods.insert(std::pair<string, int>("MDTracking", 2));
    mapRecoMethods.insert(std::pair<string, int>("MDOverlap", 3));
    mapRecoMethods.insert(std::pair<string, int>("CenterGravity", 4));
    mapRecoMethods.insert(std::pair<string, int>("MDMultiple", 5));
    mapRecoMethods.insert(std::pair<string, int>("HalfReco", 6));
    mapRecoMethods.insert(std::pair<string, int>("MDGap", 7));
    mapRecoMethods.insert(std::pair<string, int>("EdgeMethod", 8));
 
    std::vector<int> vecFibSel;
    vecFibSel.clear();
 
    switch(mapRecoMethods[strAlgo])
    {
    case 1:
        {
            ALFA_ODTracking* pODTracking = new ALFA_ODTracking();
 
            sc = pODTracking->Initialize(m_iMultiplicityCutOD, m_fDistanceCutOD, m_iLayerCutOD, m_iDataType);
            sc = pODTracking->Execute(eRPName-1, ListODHits, m_faOD, m_fbOD);
            sc = pODTracking->Finalize(&listODResults);
 
            std::list<ODRESULT>::const_iterator iter;
            for(iter=listODResults.begin(); iter!=listODResults.end(); ++iter)
            {
                ODResults.iSide   = (*iter).iSide;
                ODResults.fRecPos = (*iter).fRecPos;
                ODResults.fOverY  = (*iter).fOverY;
                ODResults.iNumY   = (*iter).iNumY;
 
                if ((ODResults.fRecPos != -9999.0))
                {
                    Int_t iFibSel[ODSIDESCNT][ODPLATESCNT];
                    std::fill_n(&iFibSel[0][0], sizeof(iFibSel)/sizeof(Int_t), -9999);
 
                    pODTracking->GetData(iFibSel);
                    vecFibSel.clear();
                    for (int iPlate=0; iPlate<ODPLATESCNT; iPlate++)
                    {
                        vecFibSel.push_back(iFibSel[ODResults.iSide][iPlate]);
                    }
 
                    m_pLocRecODEvCollection->push_back(new ALFA_LocRecODEvent(1, eRPName-1, ODResults.iSide, ODResults.fRecPos, ODResults.fOverY, ODResults.iNumY, vecFibSel));
                }
            }
 
            delete pODTracking;
            break;
        }
    case 2:
        {
            ALFA_MDTracking* pMDTracking = new ALFA_MDTracking();
            sc = pMDTracking->Initialize(m_faMD, m_fbMD, m_iMultiplicityCutMD, m_iUVCutMD, m_fOverlapCutMD);
            sc = pMDTracking->Execute(eRPName-1, ListMDHits);
            sc = pMDTracking->Finalize(fRecPosX, fRecPosY);
 
            if (fRecPosX!=-9999.0 && fRecPosY!=-9999.0)
            {
                Int_t iFibSel[ALFALAYERSCNT*ALFAPLATESCNT];
                Int_t iNumU = 0, iNumV = 0;
                Float_t fOverlapU = -9999.0, fOverlapV = -9999.0;
                std::fill_n(&iFibSel[0], sizeof(iFibSel)/sizeof(Int_t), -9999);
 
                pMDTracking->GetData(iNumU, iNumV, fOverlapU, fOverlapV, iFibSel);
                vecFibSel.clear();
                for (int & iLayer : iFibSel)
                {
                    vecFibSel.push_back(iLayer);
                }
 
                m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(2, eRPName-1, fRecPosX, fRecPosY, fOverlapU, fOverlapV, iNumU, iNumV, std::move(vecFibSel)));
            }
 
            delete pMDTracking;
            break;
        }
    case 3:
        {
            ALFA_MDOverlap* pMDOverlap = new ALFA_MDOverlap();
 
            sc = pMDOverlap->Initialize(eRPName, ListMDHits, m_faMD, m_fbMD, m_fzMD);
            sc = pMDOverlap->Execute();
            sc = pMDOverlap->Finalize(fRecPosX, fRecPosY);
 
            Int_t iFibSel[ALFALAYERSCNT*ALFAPLATESCNT];
            Int_t iNumU = 0, iNumV = 0;
            Float_t fOverlapU = -9999.0, fOverlapV = -9999.0;
            std::fill_n(&iFibSel[0], sizeof(iFibSel)/sizeof(Int_t), -9999);
 
            pMDOverlap->GetData(iFibSel);
            if (fRecPosX!=-9999 && fRecPosY!=-9999)
            {
                vecFibSel.clear();
                for (int & iLayer : iFibSel)
                {
                    vecFibSel.push_back(iLayer);
                }
 
                m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(3, eRPName-1, fRecPosX, fRecPosY, fOverlapU, fOverlapV, iNumU, iNumV, std::move(vecFibSel)));
            }
 
            delete pMDOverlap;
            break;
        }
    case 4:
        {
            ALFA_CenterGravity* pCenterGravity = new ALFA_CenterGravity();
 
            sc = pCenterGravity->Initialize(eRPName, ListMDHits, m_faMD, m_fbMD, m_fzMD);
            sc = pCenterGravity->Execute();
            sc = pCenterGravity->Finalize(fRecPosX, fRecPosY);
 
            Int_t iFibSel[ALFALAYERSCNT*ALFAPLATESCNT];
            Int_t iNumU = 0, iNumV = 0;
            Float_t fOverlapU = -9999.0, fOverlapV = -9999.0;
            std::fill_n(&iFibSel[0], sizeof(iFibSel)/sizeof(Int_t), -9999);
 
            pCenterGravity->GetData(iFibSel);
            if (fRecPosX!=-9999 && fRecPosY!=-9999)
            {
                vecFibSel.clear();
                for (int & iLayer : iFibSel)
                {
                    vecFibSel.push_back(iLayer);
                }
 
                m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(4, eRPName-1, fRecPosX, fRecPosY, fOverlapU, fOverlapV, iNumU, iNumV, std::move(vecFibSel)));
            }
 
            delete pCenterGravity;
            break;
        }
    case 5:
        {
            // MDMultiple algorihm ------------------------------------------------------
            Float_t fRecPosX[NTRACK];
            Float_t fRecPosY[NTRACK];
            ALFA_MDMultiple* pMDMultiple = new ALFA_MDMultiple();
 
            sc = pMDMultiple->Initialize(eRPName-1, m_faMD, m_fbMD, m_iMultiplicityCutMD, m_iNumLayerCutMD, m_iUVCutMD, m_fOverlapCutMD);
            sc = pMDMultiple->Execute(ListMDHits);
            sc = pMDMultiple->Finalize(fRecPosX, fRecPosY);
 
            Int_t iNumU[NTRACK], iNumV[NTRACK], iFibSel[NTRACK][ALFALAYERSCNT*ALFAPLATESCNT];
            Float_t fOverU[NTRACK], fOverV[NTRACK];
            std::fill_n(&iFibSel[0][0], sizeof(iFibSel)/sizeof(Int_t), -9999);
            std::fill_n(&fOverU[0], sizeof(fOverU)/sizeof(Float_t), -9999.0);
            std::fill_n(&fOverV[0], sizeof(fOverV)/sizeof(Float_t), -9999.0);
            memset(&iNumU, 0, sizeof(iNumU));
            memset(&iNumV, 0, sizeof(iNumV));
 
            pMDMultiple->GetData(iNumU, iNumV, fOverU, fOverV, iFibSel);
            for (Int_t iTrack=0; iTrack<NTRACK; iTrack++)
            {
                if (fRecPosX[iTrack] != -9999.0 && fRecPosY[iTrack] != -9999.0)
                {
                    vecFibSel.clear();
                    for (Int_t iLayer=0; iLayer<ALFALAYERSCNT*ALFAPLATESCNT; iLayer++)
                    {
                        vecFibSel.push_back(iFibSel[iTrack][iLayer]);
                    }
 
                    m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(5, eRPName-1, fRecPosX[iTrack], fRecPosY[iTrack], fOverU[iTrack], fOverV[iTrack], iNumU[iTrack], iNumV[iTrack], vecFibSel));
                }
            }
 
            delete pMDMultiple;
            break;
        }
    case 6:
        {
        //TODO - cele prepracovat
            ALFA_HalfReco* pHalfReco = new ALFA_HalfReco();
 
            Int_t iNumU[2] = {0, 0};
            Int_t iNumV[2] = {0, 0};
            Int_t iHalfFirst=0;
            Int_t iHalfSecond=1;
            Float_t fOverlapU[2] = {-9999, -9999};
            Float_t fOverlapV[2] = {-9999, -9999};
            Float_t fRecPosX[2] = {-9999, -9999};
            Float_t fRecPosY[2] = {-9999, -9999};
 
            Int_t iFibSel[ALFALAYERSCNT*ALFAPLATESCNT];
            std::fill_n(&iFibSel[0], sizeof(iFibSel)/sizeof(Int_t), -9999);
 
            std::vector<int> vecFibSel0;
            std::vector<int> vecFibSel1;
            vecFibSel0.clear();
            vecFibSel1.clear();
 
            if (eRPName==1 || eRPName==3 || eRPName==6 || eRPName==8)
            {
                iHalfFirst  = 0;
                iHalfSecond = 1;
            }
            else
            {
                iHalfFirst  = 1;
                iHalfSecond = 0;
            }
 
            sc = pHalfReco->Initialize(m_faMD, m_fbMD, iHalfFirst, m_iMultiplicityCutMD, m_iUVCutMDHalfReco, m_fOverlapCutMD);
            sc = pHalfReco->Execute(eRPName-1, ListMDHits);
            sc = pHalfReco->Finalize(fRecPosX[0], fRecPosY[0]);
 
            pHalfReco->GetData(iNumU[0], iNumV[0], fOverlapU[0], fOverlapV[0], iFibSel);
            vecFibSel0.clear();
            for (int & iLayer : iFibSel)
            {
                vecFibSel0.push_back(iLayer);
            }
 
            sc = pHalfReco->Initialize(m_faMD, m_fbMD, iHalfSecond, m_iMultiplicityCutMD, m_iUVCutMDHalfReco, m_fOverlapCutMD);
            sc = pHalfReco->Execute(eRPName-1, ListMDHits);
            sc = pHalfReco->Finalize(fRecPosX[1], fRecPosY[1]);
 
            pHalfReco->GetData(iNumU[1], iNumV[1], fOverlapU[1], fOverlapV[1], iFibSel);
            vecFibSel1.clear();
            for (int & iLayer : iFibSel)
            {
                vecFibSel1.push_back(iLayer);
            }
 
            if (fRecPosX[0]!=-9999.0 && fRecPosY[0]!=-9999.0 && fRecPosX[1]!=-9999.0 && fRecPosY[1]!=-9999.0)
            {
                m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(6, eRPName-1, fRecPosX[0], fRecPosY[0], fOverlapU[0], fOverlapV[0], iNumU[0], iNumV[0], std::move(vecFibSel0)));
                m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(6, eRPName-1, fRecPosX[1], fRecPosY[1], fOverlapU[1], fOverlapV[1], iNumU[1], iNumV[1], std::move(vecFibSel1)));
            }
 
            delete pHalfReco;
            break;
        }
    case 7:
        {
            Float_t fRecPosX[NTRACK];
            Float_t fRecPosY[NTRACK];
            ALFA_MDGap* pMDGap = new ALFA_MDGap();
 
            sc = pMDGap->Initialize(eRPName-1, m_faMD, m_fbMD, m_iUVCutMD, m_fOverlapCutMD);
            sc = pMDGap->Execute(ListMDHits);
            sc = pMDGap->Finalize(fRecPosX, fRecPosY);
 
            Int_t iNumU[NTRACK], iNumV[NTRACK];
            Int_t iFibSel[ALFALAYERSCNT*ALFAPLATESCNT];
            Float_t fOverU[NTRACK], fOverV[NTRACK];
 
            std::fill_n(&iFibSel[0], sizeof(iFibSel)/sizeof(Int_t), -9999);
            std::fill_n(&fOverU[0], sizeof(fOverU)/sizeof(Float_t), -9999.0);
            std::fill_n(&fOverV[0], sizeof(fOverV)/sizeof(Float_t), -9999.0);
 
            pMDGap->GetData(iNumU, iNumV, fOverU, fOverV, iFibSel);
            for (Int_t iTrack=0; iTrack<NTRACK; iTrack++)
            {
                if (fRecPosX[iTrack] != -9999.0 && fRecPosY[iTrack] != -9999.0)
                {
                    vecFibSel.clear();
                    for (int & iLayer : iFibSel)
                    {
                        vecFibSel.push_back(iLayer);
                    }
 
                    m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(7, eRPName-1, fRecPosX[iTrack], fRecPosY[iTrack], fOverU[iTrack], fOverV[iTrack], iNumU[iTrack], iNumV[iTrack], vecFibSel));
                }
            }
 
            delete pMDGap;
            break;
        }
    case 8:
        {
            Float_t fRecPosX[NTRACK];
            Float_t fRecPosY[NTRACK];
            ALFA_EdgeMethod* pEdgeMethod = new ALFA_EdgeMethod(m_bEdgeMethod_Opt_Sisters, m_bEdgeMethod_Opt_UseGaps);
            pEdgeMethod->Initialize(eRPName-1, m_faMD, m_fbMD, ListMDHits);
 
            Int_t iNumU[NTRACK], iNumV[NTRACK];
            Int_t iFibSel[ALFALAYERSCNT*ALFAPLATESCNT];
            Float_t fOverU[NTRACK], fOverV[NTRACK];
 
            std::fill_n(&iFibSel[0], sizeof(iFibSel)/sizeof(Int_t), -9999);
            std::fill_n(&fOverU[0], sizeof(fOverU)/sizeof(Float_t), -9999.0);
            std::fill_n(&fOverV[0], sizeof(fOverV)/sizeof(Float_t), -9999.0);
 
 
            std::vector<ALFA_EdgeMethod::Track> tracks;
            pEdgeMethod->EdgeMethod(eRPName - 1, tracks);
 
            for (UInt_t iTrack=0; iTrack<tracks.size(); iTrack++)
            {
                pEdgeMethod->selectedFibers(eRPName - 1, tracks.at(iTrack), iFibSel);
 
                iNumU[iTrack]    = tracks.at(iTrack).first.second;
                iNumV[iTrack]    = tracks.at(iTrack).second.second;
                fOverU[iTrack]   = tracks.at(iTrack).first.first.second;
                fOverV[iTrack]   = tracks.at(iTrack).second.first.second;
                fRecPosX[iTrack] = 0.5*TMath::Sqrt2()*(-tracks.at(iTrack).first.first.first-tracks.at(iTrack).second.first.first);
                fRecPosY[iTrack] = 0.5*TMath::Sqrt2()*(-tracks.at(iTrack).first.first.first+tracks.at(iTrack).second.first.first);
 
                if (fRecPosX[iTrack] != -9999.0 && fRecPosY[iTrack] != -9999.0)
                {
                    vecFibSel.clear();
                    for (int & iLayer : iFibSel)
                    {
                        vecFibSel.push_back(iLayer);
                    }
 
                    m_pLocRecEvCollection->push_back(new ALFA_LocRecEvent(8, eRPName-1, fRecPosX[iTrack], fRecPosY[iTrack], fOverU[iTrack], fOverV[iTrack], iNumU[iTrack], iNumV[iTrack], vecFibSel));
                }
            }
 
            delete pEdgeMethod;
            sc = StatusCode::SUCCESS;
            break;
        }
    default:
        {
            ATH_MSG_ERROR("Unable to recognize selected algorithm");
            std::cout << "strAlgo: " << strAlgo << std::endl;
            return StatusCode::SUCCESS;
        }
    }
 
    ATH_MSG_DEBUG("end ALFA_LocRec::ExecuteRecoMethod()");
 
    return sc;
}

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 ALFA_LocRec::finalize

(

)

Definition at line 337 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::finalize()");
 
 
 
    ATH_MSG_DEBUG("end ALFA_LocRec::finalize()");
    return StatusCode::SUCCESS;
}

initialize()

StatusCode ALFA_LocRec::initialize

(

)

Definition at line 156 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::initialize()");
 
    StatusCode sc;
    ClearGeometry();
 
//  sc = service("StoreGateSvc",m_storeGate);
//  if(sc.isFailure())
//  {
//      ATH_MSG_ERROR("reconstruction: unable to retrieve pointer to StoreGateSvc");
//      return sc;
//  }
    
    //read geometry
    if(ReadGeometryDetCS())
    {
        ATH_MSG_DEBUG("Geometry loaded successfully");
    }
    else
    {
        ATH_MSG_FATAL("Could not load geometry");
        return StatusCode::FAILURE;
    }
 
    //write ALFA_GeometryReader to StoreGate
//  if (StatusCode::SUCCESS!=service("DetectorStore",m_pDetStore))
//  {
//      ATH_MSG_ERROR("Detector store not found");
//      return StatusCode::FAILURE;
//  }
    sc = detStore()->record(m_pGeometryReader, m_strKeyGeometryForReco);
    if(sc.isFailure())
    {
        ATH_MSG_ERROR("m_pGeometryReader: unable to record to StoreGate");
        return sc;
    }
 
//  if (!detStore()->contains<ALFA_GeometryReader>(m_strKeyGeometryForReco))
//  {
//      ATH_MSG_ERROR("m_pGeometryReader is not in StoreGate");
//  }
//  else
//  {
//      ATH_MSG_DEBUG("m_pGeometryReader is in StoreGate");
//  }
 
    ATH_CHECK(m_eventInfoKey.initialize());
 
    m_iEvent = 0;
 
    ATH_MSG_DEBUG("end ALFA_LocRec::initialize()");
    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()

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.

ReadGeometryDetCS()

bool ALFA_LocRec::ReadGeometryDetCS ( )

private

Definition at line 421 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::ReadGeometryDetCS()");
 
    Int_t nPlateID, nFiberID;
    FIBERPARAMS FiberParams;
    eRPotName eRPName;
 
    bool bRes=false;
    std::list<eRPotName>::const_iterator iterRPName;
 
    if((bRes=m_pGeometryReader->Initialize(&m_Config,EFCS_ATLAS))==true)
    {
        ATH_MSG_DEBUG("Geometry successfully initialized");
        m_pGeometryReader->GetListOfExistingRPotIDs(&m_ListExistingRPots);
    }
    else
    {
        ATH_MSG_FATAL("Cannot initialize geometry");
        return bRes;
    }
 
    for(iterRPName=m_ListExistingRPots.begin();iterRPName!=m_ListExistingRPots.end();++iterRPName)
    {
        eRPName = *iterRPName;
 
        //for MD-fibers
        for (nPlateID = 0; nPlateID < ALFAPLATESCNT; nPlateID++)
        {
            for (nFiberID = 0; nFiberID < ALFAFIBERSCNT; nFiberID++)
            {
                //for V-fibers
                if (m_pGeometryReader->GetVFiberParams(&FiberParams, eRPName, nPlateID+1, nFiberID+1))
                {
                    m_faMD[eRPName-1][2*nPlateID+1][nFiberID] = FiberParams.fSlope;
                    m_fbMD[eRPName-1][2*nPlateID+1][nFiberID] = FiberParams.fOffset;
                    m_fzMD[eRPName-1][2*nPlateID+1][nFiberID] = FiberParams.fZPos;
                }
                else
                {
                    ATH_MSG_ERROR("Unable to get FiberV parameters");
                    return false;
                }
 
                //for U-fibers
                if (m_pGeometryReader->GetUFiberParams(&FiberParams, eRPName, nPlateID+1, nFiberID+1))
                {
                    m_faMD[eRPName-1][2*nPlateID][nFiberID] = FiberParams.fSlope;
                    m_fbMD[eRPName-1][2*nPlateID][nFiberID] = FiberParams.fOffset;
                    m_fzMD[eRPName-1][2*nPlateID][nFiberID] = FiberParams.fZPos;
                }
                else
                {
                    ATH_MSG_ERROR("Unable to get FiberU parameters");
                    return false;
                }
            }
        }
 
        //m_ODfiber = (m_ODside==0)? m_ODfiber : 29-m_ODfiber;
        //for OD-fibers
        for (nPlateID = 0; nPlateID < ODPLATESCNT; nPlateID++)
        {
            for (nFiberID = 0; nFiberID < ODFIBERSCNT; nFiberID++)
            {
                //(+16 because U0-nFiberID is indexed from 16 to 30 in Geometry package)
                if (m_pGeometryReader->GetODFiberParams(&FiberParams, EFT_ODFIBERU0, eRPName, nPlateID+1, nFiberID+16))
                {
                    m_faOD[eRPName-1][nPlateID][0][nFiberID+15] = FiberParams.fSlope;
                    m_fbOD[eRPName-1][nPlateID][0][nFiberID+15] = FiberParams.fOffset;
                    m_fzOD[eRPName-1][nPlateID][0][nFiberID+15] = FiberParams.fZPos;
                }
                else
                {
                    ATH_MSG_ERROR("Unable to get ODFiberU0 parameters");
                    return false;
                }
 
                if (m_pGeometryReader->GetODFiberParams(&FiberParams, EFT_ODFIBERV0, eRPName, nPlateID+1, nFiberID+1))
                {
                    m_faOD[eRPName-1][nPlateID][0][nFiberID] = FiberParams.fSlope;
                    m_fbOD[eRPName-1][nPlateID][0][nFiberID] = FiberParams.fOffset;
                    m_fzOD[eRPName-1][nPlateID][0][nFiberID] = FiberParams.fZPos;
                }
                else
                {
                    ATH_MSG_ERROR("Unable to get ODFiberV0 parameters");
                    return false;
                }
 
                if (m_pGeometryReader->GetODFiberParams(&FiberParams, EFT_ODFIBERU1, eRPName, nPlateID+1, nFiberID+16))
                {
                    if (m_iDataType==1) //for the tunnel and testbeam data
                    {
                        m_faOD[eRPName-1][nPlateID][1][14-nFiberID] = FiberParams.fSlope;
                        m_fbOD[eRPName-1][nPlateID][1][14-nFiberID] = FiberParams.fOffset;
                        m_fzOD[eRPName-1][nPlateID][1][14-nFiberID] = FiberParams.fZPos;
                    }
                    else                //simulation
                    {
                        m_faOD[eRPName-1][nPlateID][1][nFiberID+15] = FiberParams.fSlope;
                        m_fbOD[eRPName-1][nPlateID][1][nFiberID+15] = FiberParams.fOffset;
                        m_fzOD[eRPName-1][nPlateID][1][nFiberID+15] = FiberParams.fZPos;
                    }
                }
                else
                {
                    ATH_MSG_ERROR("Unable to get ODFiberU1 parameters");
                    return false;
                }
                //(+16 because V1-nFiberID is indexed from 16 to 30 in Geometry package)
                if (m_pGeometryReader->GetODFiberParams(&FiberParams, EFT_ODFIBERV1, eRPName, nPlateID+1, nFiberID+1))
                {
                    if (m_iDataType==1) //for the tunnel and testbeam data
                    {
                        m_faOD[eRPName-1][nPlateID][1][29-nFiberID] = FiberParams.fSlope;
                        m_fbOD[eRPName-1][nPlateID][1][29-nFiberID] = FiberParams.fOffset;
                        m_fzOD[eRPName-1][nPlateID][1][29-nFiberID] = FiberParams.fZPos;
                    }
                    else                //simulation
                    {
                        m_faOD[eRPName-1][nPlateID][1][nFiberID] = FiberParams.fSlope;
                        m_fbOD[eRPName-1][nPlateID][1][nFiberID] = FiberParams.fOffset;
                        m_fzOD[eRPName-1][nPlateID][1][nFiberID] = FiberParams.fZPos;
                    }
                }
                else
                {
                    ATH_MSG_ERROR("Unable to get ODFiberV1 parameters");
                    return false;
                }
            }
        }
    }
 
//  SaveGeometry();
 
    ATH_MSG_DEBUG("end ALFA_LocRec::ReadGeometryDetCS()");
 
    return bRes;
}

RecordCollection()

StatusCode ALFA_LocRec::RecordCollection ( )

private

Definition at line 1021 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::RecordCollection()");
 
    StatusCode sc = StatusCode::SUCCESS;
 
    m_pLocRecEvCollection = new ALFA_LocRecEvCollection();
    sc = evtStore()->record(m_pLocRecEvCollection, m_strKeyLocRecEvCollection);
 
    if (sc.isFailure())
    {
        ATH_MSG_ERROR(m_strAlgoMD << "MD - Could not record the empty LocRecEv collection in StoreGate");
        return sc;
    }
    else
    {
        ATH_MSG_DEBUG("MD - LocRecEv collection was recorded in StoreGate");
    }
 
    ATH_MSG_DEBUG("end ALFA_LocRec::RecordCollection()");
 
    return sc;
}

RecordODCollection()

StatusCode ALFA_LocRec::RecordODCollection ( )

private

Definition at line 1045 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::RecordODCollection()");
 
    StatusCode sc = StatusCode::SUCCESS;
 
    m_pLocRecODEvCollection = new ALFA_LocRecODEvCollection();
    sc = evtStore()->record(m_pLocRecODEvCollection, m_strKeyLocRecODEvCollection);
 
    if (sc.isFailure())
    {
        ATH_MSG_ERROR(m_strAlgoOD << "OD - Could not record the empty LocRecEv collection in StoreGate");
 
        return sc;
    }
    else
    {
         ATH_MSG_DEBUG("OD - LocRecEv collection is recorded in StoreGate");
    }
 
    ATH_MSG_DEBUG("end ALFA_LocRec::RecordODCollection()");
 
    return sc;
}

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

SaveGeometry()

void ALFA_LocRec::SaveGeometry ( )

private

Definition at line 614 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::SaveGeometry()");
 
 
    eGeoSourceType eGeoType;
    char szFilename[64];
    std::list<eRPotName>::const_iterator iterRPName;
 
    for(iterRPName=m_ListExistingRPots.begin();iterRPName!=m_ListExistingRPots.end();++iterRPName)
    {
        eGeoType=m_pGeometryReader->GetRPGeometryType(*iterRPName, EFT_FIBERMD);
        switch(eGeoType)
        {
            case EGST_IDEALGEOMETRY:
                sprintf(szFilename,"ALFA_LocRec_Idealgeometry_RP-%s.txt",m_pGeometryReader->GetRPotLabel(*iterRPName));
                StoreReconstructionGeometry(*iterRPName, EFT_FIBERMD, szFilename);
                ATH_MSG_DEBUG("The ALFA Fiber geometry was stored in the "<<szFilename<<" file");
                break;
            case EGST_FILE:
                sprintf(szFilename,"ALFA_LocRec_Realgeometry_RP-%s.txt",m_pGeometryReader->GetRPotLabel(*iterRPName));
                StoreReconstructionGeometry(*iterRPName, EFT_FIBERMD, szFilename);
                ATH_MSG_DEBUG("The ALFA Fiber geometry was stored in the "<<szFilename<<" file");
                break;
            case EGST_DATABASE:
                sprintf(szFilename,"ALFA_LocRec_Realdbgeometry_RP-%s.txt",m_pGeometryReader->GetRPotLabel(*iterRPName));
                StoreReconstructionGeometry(*iterRPName, EFT_FIBERMD, szFilename);
                ATH_MSG_DEBUG("The ALFA Fiber geometry was stored in the "<<szFilename<<" file");
                break;
            default:
                ATH_MSG_DEBUG("Unrecognized MD geometry!");
                return;
                break;
        }
 
        eGeoType=m_pGeometryReader->GetRPGeometryType(*iterRPName, EFT_FIBEROD);
        switch(eGeoType)
        {
            case EGST_IDEALGEOMETRY:
                sprintf(szFilename,"ALFA_LocRec_Idealgeometry_OD_RP-%s.txt",m_pGeometryReader->GetRPotLabel(*iterRPName));
                StoreReconstructionGeometry(*iterRPName, EFT_FIBEROD, szFilename);
                ATH_MSG_DEBUG("The ODFiber geometry was stored in the "<<szFilename<<" file");
                break;
            case EGST_FILE:
                sprintf(szFilename,"ALFA_LocRec_Realgeometry_OD_RP-%s.txt",m_pGeometryReader->GetRPotLabel(*iterRPName));
                StoreReconstructionGeometry(*iterRPName, EFT_FIBEROD, szFilename);
                ATH_MSG_DEBUG("The ODFiber geometry was stored in the "<<szFilename<<" file");
                break;
            case EGST_DATABASE:
                sprintf(szFilename,"ALFA_LocRec_Realdbgeometry_OD_RP-%s.txt",m_pGeometryReader->GetRPotLabel(*iterRPName));
                StoreReconstructionGeometry(*iterRPName, EFT_FIBEROD, szFilename);
                ATH_MSG_DEBUG("The ODFiber geometry was stored in the "<<szFilename<<" file");
                break;
            default:
                ATH_MSG_DEBUG("Unrecognized OD geometry!");
                return;
                break;
        }
    }
 
    ATH_MSG_DEBUG("end ALFA_LocRec::SaveGeometry()");
}

StoreReconstructionGeometry()

bool ALFA_LocRec::StoreReconstructionGeometry ( const eRPotName eRPName, const eFiberType eFType, const char * szDataDestination )

private

Definition at line 563 of file ALFA_LocRec.cxx.

{
    ATH_MSG_DEBUG("begin ALFA_LocRec::StoreReconstructionGeometry()");
 
    Int_t iLayer, iFiber, iPlate, iSide;
    Float_t fParamB, fZ, fSlope;
    FILE * pFile;
 
    pFile = fopen(szDataDestination, "w");
    if(pFile==NULL) return false;
    fprintf(pFile, "xxxxxxxxxxxxxxxxxxx\n");
    if (eFType==EFT_FIBERMD)
    {
        fprintf(pFile, "20\n");
        for(iLayer=0; iLayer<2*ALFAPLATESCNT; iLayer++)
        {
            for(iFiber=0; iFiber<ALFAFIBERSCNT; iFiber++)
            {
                fSlope=m_faMD[eRPName-1][iLayer][iFiber];
                fParamB=m_fbMD[eRPName-1][iLayer][iFiber];
                fZ=m_fzMD[eRPName-1][iLayer][iFiber];
 
                fprintf(pFile, " %2d  %2d  %7.5f  %7.4f  %7.3f\n", iLayer+1, iFiber+1, fSlope, fParamB, fZ);
            }
        }
    }
    else if (eFType==EFT_FIBEROD)
    {
        fprintf(pFile, "6\n");
        for(iPlate=0; iPlate<ODPLATESCNT; iPlate++)
        {
            for (iSide=0; iSide<ODSIDESCNT; iSide++)
            {
                for(iFiber=0; iFiber<ODLAYERSCNT*ODFIBERSCNT; iFiber++)
                {
                    fSlope=m_faOD[eRPName-1][iPlate][iSide][iFiber];
                    fParamB=m_fbOD[eRPName-1][iPlate][iSide][iFiber];
                    fZ=m_fzOD[eRPName-1][iPlate][iSide][iFiber];
 
                    fprintf(pFile, " %2d  %2d  %7.5f  %7.4f  %7.3f\n", (2*iPlate+iSide+1), iFiber+1, fSlope, fParamB, fZ);
                }
            }
        }
    }
    fclose(pFile);
 
    ATH_MSG_DEBUG("end ALFA_LocRec::StoreReconstructionGeometry()");
 
    return true;
}

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

Member Data Documentation

m_bEdgeMethod_Opt_Sisters

Bool_t ALFA_LocRec::m_bEdgeMethod_Opt_Sisters

private

Definition at line 75 of file ALFA_LocRec.h.

m_bEdgeMethod_Opt_UseGaps

Bool_t ALFA_LocRec::m_bEdgeMethod_Opt_UseGaps

private

Definition at line 76 of file ALFA_LocRec.h.

m_Config

GEOMETRYCONFIGURATION ALFA_LocRec::m_Config

private

Definition at line 50 of file ALFA_LocRec.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_eventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> ALFA_LocRec::m_eventInfoKey {this, "EvtInfo", "EventInfo", "EventInfo name"}

private

Definition at line 59 of file ALFA_LocRec.h.

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

m_eventNum

UInt_t ALFA_LocRec::m_eventNum

private

Definition at line 64 of file ALFA_LocRec.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_faMD

Float_t ALFA_LocRec::m_faMD[RPOTSCNT][ALFALAYERSCNT *ALFAPLATESCNT][ALFAFIBERSCNT]

private

Definition at line 79 of file ALFA_LocRec.h.

m_faOD

Float_t ALFA_LocRec::m_faOD[RPOTSCNT][ODPLATESCNT][ODSIDESCNT][ODLAYERSCNT *ODFIBERSCNT]

private

Definition at line 84 of file ALFA_LocRec.h.

m_fbMD

Float_t ALFA_LocRec::m_fbMD[RPOTSCNT][ALFALAYERSCNT *ALFAPLATESCNT][ALFAFIBERSCNT]

private

Definition at line 80 of file ALFA_LocRec.h.

m_fbOD

Float_t ALFA_LocRec::m_fbOD[RPOTSCNT][ODPLATESCNT][ODSIDESCNT][ODLAYERSCNT *ODFIBERSCNT]

private

Definition at line 85 of file ALFA_LocRec.h.

m_fDistanceCutOD

Float_t ALFA_LocRec::m_fDistanceCutOD

private

Definition at line 73 of file ALFA_LocRec.h.

m_fOverlapCutMD

Float_t ALFA_LocRec::m_fOverlapCutMD

private

Definition at line 72 of file ALFA_LocRec.h.

m_fzMD

Float_t ALFA_LocRec::m_fzMD[RPOTSCNT][ALFALAYERSCNT *ALFAPLATESCNT][ALFAFIBERSCNT]

private

Definition at line 81 of file ALFA_LocRec.h.

m_fzOD

Float_t ALFA_LocRec::m_fzOD[RPOTSCNT][ODPLATESCNT][ODSIDESCNT][ODLAYERSCNT *ODFIBERSCNT]

private

Definition at line 86 of file ALFA_LocRec.h.

m_iDataType

Int_t ALFA_LocRec::m_iDataType

private

Definition at line 65 of file ALFA_LocRec.h.

m_iEvent

Int_t ALFA_LocRec::m_iEvent

private

Definition at line 66 of file ALFA_LocRec.h.

m_iLayerCutOD

Float_t ALFA_LocRec::m_iLayerCutOD

private

Definition at line 74 of file ALFA_LocRec.h.

m_iMultiplicityCutMD

Int_t ALFA_LocRec::m_iMultiplicityCutMD

private

Definition at line 67 of file ALFA_LocRec.h.

m_iMultiplicityCutOD

Int_t ALFA_LocRec::m_iMultiplicityCutOD

private

Definition at line 71 of file ALFA_LocRec.h.

m_iNumLayerCutMD

Int_t ALFA_LocRec::m_iNumLayerCutMD

private

Definition at line 68 of file ALFA_LocRec.h.

m_iUVCutMD

Int_t ALFA_LocRec::m_iUVCutMD

private

Definition at line 69 of file ALFA_LocRec.h.

m_iUVCutMDHalfReco

Int_t ALFA_LocRec::m_iUVCutMDHalfReco

private

Definition at line 70 of file ALFA_LocRec.h.

m_ListExistingRPots

std::list<eRPotName> ALFA_LocRec::m_ListExistingRPots

private

Definition at line 62 of file ALFA_LocRec.h.

m_pGeometryReader

ALFA_GeometryReader* ALFA_LocRec::m_pGeometryReader

private

Definition at line 51 of file ALFA_LocRec.h.

m_pLocRecEvCollection

ALFA_LocRecEvCollection* ALFA_LocRec::m_pLocRecEvCollection

private

Definition at line 53 of file ALFA_LocRec.h.

m_pLocRecEvent

ALFA_LocRecEvent* ALFA_LocRec::m_pLocRecEvent

private

Definition at line 54 of file ALFA_LocRec.h.

m_pLocRecODEvCollection

ALFA_LocRecODEvCollection* ALFA_LocRec::m_pLocRecODEvCollection

private

Definition at line 56 of file ALFA_LocRec.h.

m_pLocRecODEvent

ALFA_LocRecODEvent* ALFA_LocRec::m_pLocRecODEvent

private

Definition at line 57 of file ALFA_LocRec.h.

m_strAlgoMD

std::string ALFA_LocRec::m_strAlgoMD

private

Definition at line 92 of file ALFA_LocRec.h.

m_strAlgoOD

std::string ALFA_LocRec::m_strAlgoOD

private

Definition at line 91 of file ALFA_LocRec.h.

m_strCollectionName

std::string ALFA_LocRec::m_strCollectionName

private

Definition at line 95 of file ALFA_LocRec.h.

m_strKeyGeometryForReco

std::string ALFA_LocRec::m_strKeyGeometryForReco

private

Definition at line 88 of file ALFA_LocRec.h.

m_strKeyLocRecEvCollection

std::string ALFA_LocRec::m_strKeyLocRecEvCollection

private

Definition at line 93 of file ALFA_LocRec.h.

m_strKeyLocRecODEvCollection

std::string ALFA_LocRec::m_strKeyLocRecODEvCollection

private

Definition at line 94 of file ALFA_LocRec.h.

m_strODCollectionName

std::string ALFA_LocRec::m_strODCollectionName

private

Definition at line 96 of file ALFA_LocRec.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vecListAlgoMD

std::vector<std::string> ALFA_LocRec::m_vecListAlgoMD

private

Definition at line 89 of file ALFA_LocRec.h.

m_vecListAlgoOD

std::vector<std::string> ALFA_LocRec::m_vecListAlgoOD

private

Definition at line 90 of file ALFA_LocRec.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:

• ALFA_LocRec.h
• ALFA_LocRec.cxx