#include <ALFA_GloRec.h>

Inheritance diagram for ALFA_GloRec:

Collaboration diagram for ALFA_GloRec:

Public Member Functions
	ALFA_GloRec (const std::string &name, ISvcLocator *pSvcLocator)
	~ALFA_GloRec ()
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
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
StatusCode	Truth_info ()
void	InitHistos ()
void	FillTrackCandHistos (AlfaTrackCand *trkcand)
void	WriteHistos ()
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ALFA_GloRecEvCollection *	m_pGloRecEvCollection {}
int	m_iDataType = 0
std::string	m_strTrackPathPatterns
TObjArray	m_TObjArrTrackPathPatterns {}
std::string	m_strGloRecAnalysisFile
std::string	m_strGloRecCollectionName
std::string	m_strLocRecCorrCollectionName
std::string	m_TruthCollectionName
double	m_px_g_pos {}
double	m_py_g_pos {}
double	m_pz_g_pos {}
double	m_x_g_pos {}
double	m_y_g_pos {}
double	m_z_g_pos {}
double	m_px_g_neg {}
double	m_py_g_neg {}
double	m_pz_g_neg {}
double	m_x_g_neg {}
double	m_y_g_neg {}
double	m_z_g_neg {}
TH1F *	m_th1_x_g {}
TH1F *	m_th1_y_g {}
TH1F *	m_th1_xslope_g {}
TH1F *	m_th1_yslope_g {}
TH1F *	m_th1_xnearuppotresiduals {}
TH1F *	m_th1_ynearuppotresiduals {}
TH1F *	m_th1_xfaruppotresiduals {}
TH1F *	m_th1_yfaruppotresiduals {}
TH1F *	m_th1_xnearlwpotresiduals {}
TH1F *	m_th1_ynearlwpotresiduals {}
TH1F *	m_th1_xfarlwpotresiduals {}
TH1F *	m_th1_yfarlwpotresiduals {}
TH2F *	m_th2_truexvsrecx {}
TH2F *	m_th2_trueyvsrecy {}
TH2F *	m_th2_truexslopevsrecxslope {}
TH2F *	m_th2_trueyslopevsrecyslope {}
TH1F *	m_th1_recxovertruex {}
TH1F *	m_th1_recyovertruey {}
TH1F *	m_th1_recxslopeovertruexslope {}
TH1F *	m_th1_recyslopeovertrueyslope {}
TH1F *	m_th1_recxminustruex {}
TH1F *	m_th1_recyminustruey {}
TH1F *	m_th1_recxslopeminustruexslope {}
TH1F *	m_th1_recyslopeminustrueyslope {}
TH2F *	m_th2_extrapxvsrecxnearpot {}
TH2F *	m_th2_extrapyvsrecynearpot {}
TH1F *	m_th1_recxoverextrapxnearpot {}
TH1F *	m_th1_recyoverextrapynearpot {}
TH1F *	m_th1_recxminusextrapxnearpot {}
TH1F *	m_th1_recyminusextrapynearpot {}
TH2F *	m_th2_extrapxvsrecxfarpot {}
TH2F *	m_th2_extrapyvsrecyfarpot {}
TH1F *	m_th1_recxoverextrapxfarpot {}
TH1F *	m_th1_recyoverextrapyfarpot {}
TH1F *	m_th1_recxminusextrapxfarpot {}
TH1F *	m_th1_recyminusextrapyfarpot {}
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 28 of file ALFA_GloRec.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_GloRec()

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

Definition at line 84 of file ALFA_GloRec.cxx.

                                                                        :
AthAlgorithm(name, pSvcLocator)
{
 
  MsgStream LogStream(Athena::getMessageSvc(), "ALFA_GloRec::ALFA_GloRec");
    LogStream << MSG::DEBUG << "begin ALFA_GloRec::ALFA_GloRec" << endmsg;
    // ... steerable job options 
  declareProperty("DataType", m_iDataType);
  declareProperty("TrackPathPatterns", m_strTrackPathPatterns);
  declareProperty("AnalysisFile", m_strGloRecAnalysisFile);
  
    declareProperty("OutputGloRecCollectionName",m_strGloRecCollectionName);
    declareProperty("InputLocRecCollectionName",m_strLocRecCorrCollectionName);
    declareProperty("TruthEventCollectionName",m_TruthCollectionName);   
 
}

◆ ~ALFA_GloRec()

ALFA_GloRec::~ALFA_GloRec ( )

Definition at line 109 of file ALFA_GloRec.cxx.

                          {
    MsgStream LogStream(Athena::getMessageSvc(), "ALFA_GloRec::~ALFA_GloRec");
    LogStream << MSG::DEBUG << "begin ALFA_GloRec::~ALFA_GloRec" << endmsg;
    LogStream << MSG::DEBUG << "end ALFA_GloRec::~ALFA_GloRec" << endmsg;
}

Member Function Documentation

◆ 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.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ 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.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ execute()

StatusCode ALFA_GloRec::execute ( )

private

Definition at line 154 of file ALFA_GloRec.cxx.

                                {
 
    MsgStream LogStream(Athena::getMessageSvc(), "ALFA_GloRec::execute()");
    LogStream << MSG::DEBUG << "begin ALFA_GloRec::execute()" << endmsg;
    StatusCode sc = StatusCode::SUCCESS;
 
    
 
 
  //  get truth information
  if(m_iDataType == 0){
      LogStream << MSG::DEBUG << "truth info()" << endmsg;        
      sc = Truth_info();  
      if  (sc.isFailure()) {
        LogStream << MSG::WARNING << "ALFA digi, Truth_info failed" << endmsg;  
        return StatusCode::SUCCESS; 
      }
    }
    // ... record output collection 
    m_pGloRecEvCollection = new ALFA_GloRecEvCollection();
    sc = evtStore()->record(m_pGloRecEvCollection, "ALFA_GloRecEvCollection");
    if (sc.isFailure()) {
        msg(MSG::ERROR) << "Could not create the empty LocRecEv collection in StoreGate" << endmsg;
    }
 
    // ... get LocRecCorr data
    const ALFA_LocRecCorrEvCollection* pLocRecCorrCol = nullptr;
    sc = evtStore()->retrieve(pLocRecCorrCol, m_strLocRecCorrCollectionName);
    if(sc.isFailure() || !pLocRecCorrCol)
    {
        msg(MSG::ERROR) << "Container "<< m_strLocRecCorrCollectionName <<" NOT FOUND !!!!!!!" << endmsg;
        return sc;
    }
 
    // ... collect the local hits and contain them in an auxiliary hit collection
    ALFA_LocRecCorrEvCollection::const_iterator mcColBeg = pLocRecCorrCol->begin();
    ALFA_LocRecCorrEvCollection::const_iterator mcColEnd = pLocRecCorrCol->end();
    AlfaLocalHits hits;
    for(; mcColBeg!=mcColEnd; ++mcColBeg) {
            HepGeom::Point3D<double> PointInAtlasCS;
        const ALFA_LocRecCorrEvent* pLocRecCorr_tmp = static_cast<const ALFA_LocRecCorrEvent*> (*mcColBeg);
        AlfaLocalHit l_hit(pLocRecCorr_tmp);
        hits.AddHit(l_hit);
        //PointInAtlasCS = m_pGeometryReader->GetDetPointInAtlas((eRPotName)(iRPot+1), Point3D<double>(fRecPosX, fRecPosY, 10.0));
    }
  //hits.Dump();
   
   
  // ... apply path pattern
    for(int i = 0; i < m_TObjArrTrackPathPatterns.GetEntries(); i++){
        TString l_pathpattern(((TObjString*) m_TObjArrTrackPathPatterns.At(i))->GetString());
        hits.ApplyPathPattern(l_pathpattern.Data());
        //hits.DumpPaths();
        for(int i = 0 ; i < hits.GetNpaths() ; i++){
            AlfaTrackCand trkcand(hits.GetPathHits(i));
            if(trkcand.IsValidTrack()){
                m_pGloRecEvCollection->push_back(new ALFA_GloRecEvent(trkcand.GetArmNum(),trkcand.GetX(),trkcand.GetY(),trkcand.GetXslope(),trkcand.GetXslope(),trkcand.GetNearLocRecCorr(),trkcand.GetFarLocRecCorr()));  
                FillTrackCandHistos(&trkcand);
            }
      }
      hits.ResetPaths();
  }
 
 
 
 
    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();
}

◆ FillTrackCandHistos()

void ALFA_GloRec::FillTrackCandHistos ( AlfaTrackCand * trkcand )

private

Definition at line 389 of file ALFA_GloRec.cxx.

                                                            {
    
    if(m_iDataType==0){
        double x_g, y_g, z_g, px_g, py_g, pz_g;
        if(trkcand->IsLeftTrack()){ // left is positive
            x_g = m_x_g_pos;
            y_g = m_y_g_pos;
            z_g = m_z_g_pos;
            px_g = m_px_g_pos;
            py_g = m_py_g_pos;
            pz_g = m_pz_g_pos;
        }else if(trkcand->IsRightTrack()){ // right is negative
            x_g = m_x_g_neg;
            y_g = m_y_g_neg;
            z_g = m_z_g_neg;
            px_g = m_px_g_neg;
            py_g = m_py_g_neg;
            pz_g = m_pz_g_neg;
        }else{
            ATH_MSG_ERROR("ALFA_GloRec::FillTrackCandHistos : Track candidate neither left nor right !!!");
            return;
        }
        m_th1_xnearuppotresiduals -> Fill(trkcand->GetNearImpactPoint()->x()-trkcand->GetNearLocalHit()->GetX(),1.);
        m_th1_ynearuppotresiduals -> Fill(trkcand->GetNearImpactPoint()->y()-trkcand->GetNearLocalHit()->GetY(),1.);
        m_th1_xfaruppotresiduals -> Fill(trkcand->GetFarImpactPoint()->x()-trkcand->GetFarLocalHit()->GetX(),1.);
        m_th1_yfaruppotresiduals -> Fill(trkcand->GetFarImpactPoint()->y()-trkcand->GetFarLocalHit()->GetY(),1.);
        m_th2_truexvsrecx -> Fill(x_g, trkcand->GetX(), 1.);
        m_th2_trueyvsrecy -> Fill(y_g, trkcand->GetY(), 1.);
        m_th1_recxovertruex -> Fill(trkcand->GetX()/x_g, 1.);
        m_th1_recyovertruey -> Fill(trkcand->GetY()/y_g, 1.);
        //double p_gen = TMath::Sqrt(px_g*px_g+px_g*px_g+py_g*pz_g);
        double px_norm = px_g / pz_g;//p_gen;
        double py_norm = py_g / pz_g;//p_gen;
        m_th2_truexslopevsrecxslope -> Fill(px_norm, trkcand->GetXslope(), 1.);
        m_th2_trueyslopevsrecyslope -> Fill(py_norm, trkcand->GetYslope(), 1.);
        m_th1_recxovertruex -> Fill(trkcand->GetX()/x_g, 1.);
        m_th1_recyovertruey -> Fill(trkcand->GetY()/y_g, 1.);
        m_th1_recxslopeovertruexslope -> Fill(trkcand->GetXslope()/px_norm, 1.);
        m_th1_recyslopeovertrueyslope -> Fill(trkcand->GetYslope()/py_norm, 1.);
        m_th1_recxminustruex -> Fill(trkcand->GetX()-x_g, 1.);
        m_th1_recyminustruey -> Fill(trkcand->GetY()-y_g, 1.);
        m_th1_recxslopeminustruexslope -> Fill(trkcand->GetXslope()-px_norm, 1.);
        m_th1_recyslopeminustrueyslope -> Fill(trkcand->GetYslope()-py_norm, 1.);
        // ... extrapolate x any of gen track to near pot
        double nearhitx = trkcand->GetNearLocalHit()->GetX();
        double nearhity = trkcand->GetNearLocalHit()->GetY();
        double nearhitz = trkcand->GetNearLocalHit()->GetZ();
        double genxextrtonearpot = x_g + px_norm * (nearhitz - z_g);
        double genyextrtonearpot = y_g + py_norm * (nearhitz - z_g);
        //double genzextrtonearpot = nearhitz;
        m_th2_extrapxvsrecxnearpot -> Fill(genxextrtonearpot, nearhitx,1.);
        m_th2_extrapyvsrecynearpot -> Fill(genyextrtonearpot, nearhity,1.);
        m_th1_recxoverextrapxnearpot -> Fill(nearhitx / genxextrtonearpot, 1.);
        m_th1_recyoverextrapynearpot -> Fill(nearhity / genyextrtonearpot, 1.);
        m_th1_recxminusextrapxnearpot -> Fill(nearhitx - genxextrtonearpot, 1.);
        m_th1_recyminusextrapynearpot -> Fill(nearhity - genyextrtonearpot, 1.);
        // ... extrapolate x any of gen track to far pot
        double farhitx = trkcand->GetFarLocalHit()->GetX();
        double farhity = trkcand->GetFarLocalHit()->GetY();
        double farhitz = trkcand->GetFarLocalHit()->GetZ();
        double genxextrtofarpot = x_g + px_norm * (farhitz - z_g);
        double genyextrtofarpot = y_g + py_norm * (farhitz - z_g);
        //double genzextrtofarpot = farhitz;
        m_th2_extrapxvsrecxfarpot -> Fill(genxextrtofarpot, farhitx,1.);
        m_th2_extrapyvsrecyfarpot -> Fill(genyextrtofarpot, farhity,1.);
        m_th1_recxoverextrapxfarpot -> Fill(farhitx / genxextrtofarpot, 1.);
        m_th1_recyoverextrapyfarpot -> Fill(farhity / genyextrtofarpot, 1.);
        m_th1_recxminusextrapxfarpot -> Fill(farhitx - genxextrtofarpot, 1.);
        m_th1_recyminusextrapyfarpot -> Fill(farhity - genyextrtofarpot, 1.);
    }
}

◆ finalize()

StatusCode ALFA_GloRec::finalize ( )

private

Definition at line 246 of file ALFA_GloRec.cxx.

                                 {
 
    MsgStream LogStream(Athena::getMessageSvc(), "ALFA_GloRec::finalize()");
    LogStream << MSG::DEBUG << "begin ALFA_GloRec::finalize()" << endmsg;
 
    WriteHistos();
 
    return StatusCode::SUCCESS;
}

◆ InitHistos()

void ALFA_GloRec::InitHistos ( )

private

Definition at line 340 of file ALFA_GloRec.cxx.

                            {
    m_th1_x_g = new TH1F("x_g","x^{vertex}_{gen}", 100, -5,5);
    m_th1_y_g = new TH1F("y_g","y^{vertex}_{gen}", 100, -5,5);
    m_th1_xslope_g = new TH1F("xslope_g","x'^{vertex}_{gen}", 100, -1.e-3,1.e3);
    m_th1_yslope_g = new TH1F("yslope_g","y'^{vertex}_{gen}", 100, -1.e-3,1.e3);
 
  m_th1_xnearuppotresiduals = new TH1F("xnearuppotresiduals","trk x extrap. minus near IP U loc hit",100, -10., 10.);
  m_th1_ynearuppotresiduals = new TH1F("ynearuppotresiduals","trk y extrap. minus near IP U loc hit",100, -10., 10.);
  m_th1_xfaruppotresiduals  = new TH1F("xfaruppotresiduals","trk x extrap. minus far IP U loc hit",100, -10., 10.);
  m_th1_yfaruppotresiduals  = new TH1F("yfaruppotresiduals","trk y extrap. minus far IP U loc hit",100, -10., 10.);
  m_th1_xnearlwpotresiduals = new TH1F("xnearlwpotresiduals","trk x extrap. minus near IP L loc hit",100, -10., 10.);
  m_th1_ynearlwpotresiduals = new TH1F("ynearlwpotresiduals","trk y extrap. minus near IP L loc hit",100, -10., 10.);
  m_th1_xfarlwpotresiduals  = new TH1F("xfarlwpotresiduals","trk x extrap. minus far IP L loc hit",100, -10., 10.);
  m_th1_yfarlwpotresiduals  = new TH1F("yfarlwpotresiduals","trk y extrap. minus far IP L loc hit",100, -10., 10.);
      
  m_th2_truexvsrecx = new TH2F("truexvsrecx","gen x vs rec x", 100, -50, 50, 100, -50, 50);
  m_th2_trueyvsrecy = new TH2F("trueyvsrecy","gen y vs rec y", 100, -50, 50, 100, -50, 50);
  m_th2_truexslopevsrecxslope = new TH2F("truexslopevsrecxslope","gen xslope vs rec xslope", 100,-0.002,0.002,100,-0.002,0.002);
  m_th2_trueyslopevsrecyslope = new TH2F("trueyslopevsrecyslope","gen yslope vs rec yslope", 100,-0.002,0.002,100,-0.002,0.002);
  m_th1_recxovertruex = new TH1F("recxovertruex","rec x over true x", 100, -2, 2);
  m_th1_recyovertruey = new TH1F("recyovertruey","rec y over true y", 100, -2, 2);
  m_th1_recxslopeovertruexslope = new TH1F("recxslopeovertruexslope","rec xslope over true xslope", 100,-2,2);
  m_th1_recyslopeovertrueyslope = new TH1F("recyslopeovertrueyslope","rec yslope over true yslope", 100,-2,2);
  m_th1_recxminustruex = new TH1F("recxminustruex","rec x minus true x", 100, -10, 10);
  m_th1_recyminustruey = new TH1F("recyminustruey","rec y minus true y", 100, -10, 10);
  m_th1_recxslopeminustruexslope = new TH1F("recxslopeminustruexslope","rec xslope minus true xslope", 100,-0.001,0.001);
  m_th1_recyslopeminustrueyslope = new TH1F("recyslopeminustrueyslope","rec yslope minus true yslope", 100,-0.001,0.001);
 
  m_th2_extrapxvsrecxnearpot = new TH2F("extrapxvsrecxnearpot","extrapolated gen x vs rec x nearpot", 100, -50, 50, 100, -50, 50);
  m_th2_extrapyvsrecynearpot = new TH2F("extrapyvsrecynearpot","extrapolated gen y vs rec y nearpot", 100, -50, 50, 100, -50, 50);
  m_th1_recxoverextrapxnearpot = new TH1F("recxoverextrapxnearpot","rec x over extrap x nearpot", 100, 1.-0.2, 1.+0.2);
  m_th1_recyoverextrapynearpot = new TH1F("recyoverextrapynearpot","rec y over extrap y nearpot", 100, 1.-0.2, 1.+0.2);
  m_th1_recxminusextrapxnearpot = new TH1F("recxminusextrapxnearpot","rec x minus extrap x nearpot", 100, -0.5, 0.5);
  m_th1_recyminusextrapynearpot = new TH1F("recyminusextrapynearpot","rec y minus extrap y nearpot", 100, -0.5, 0.5);
 
  m_th2_extrapxvsrecxfarpot = new TH2F("extrapxvsrecxfarpot","extrapolated gen x vs rec x farpot", 100, -50, 50, 100, -50, 50);
  m_th2_extrapyvsrecyfarpot = new TH2F("extrapyvsrecyfarpot","extrapolated gen y vs rec y farpot", 100, -50, 50, 100, -50, 50);
  m_th1_recxoverextrapxfarpot = new TH1F("recxoverextrapxfarpot","rec x over extrap x farpot", 100, 1.-0.2, 1.+0.2);
  m_th1_recyoverextrapyfarpot = new TH1F("recyoverextrapyfarpot","rec y over extrap y farpot", 100, 1.-0.2, 1.+0.2);
  m_th1_recxminusextrapxfarpot = new TH1F("recxminusextrapxfarpot","rec x minus extrap x farpot", 100, -0.5, 0.5);
  m_th1_recyminusextrapyfarpot = new TH1F("recyminusextrapyfarpot","rec y minus extrap y farpot", 100, -0.5, 0.5);
}

◆ initialize()

StatusCode ALFA_GloRec::initialize ( )

private

Definition at line 120 of file ALFA_GloRec.cxx.

                                   {
    MsgStream LogStream(Athena::getMessageSvc(), "ALFA_GloRec::initialize()");
    LogStream << MSG::DEBUG << "begin ALFA_GloRec::initialize()" << endmsg;
 
    msg(MSG::DEBUG) << "======================================" << endmsg;
    msg(MSG::DEBUG) << "ALFA_GloRec::initialize  - data type  " << endmsg;
    msg(MSG::DEBUG) << "======================================" << endmsg;
    msg(MSG::DEBUG) << "DataType = " << m_iDataType << endmsg;
 
    TString l_strtmp(m_strTrackPathPatterns);
    m_TObjArrTrackPathPatterns.Clear();
    msg(MSG::DEBUG) << "==========================================================" << endmsg;
    msg(MSG::DEBUG) << "ALFA_GloRec::initialize  - Path patterns from job options " << endmsg;
    msg(MSG::DEBUG) << "==========================================================" << endmsg;
    m_TObjArrTrackPathPatterns = *(l_strtmp.Tokenize(","));
    for(int i = 0; i < m_TObjArrTrackPathPatterns.GetEntries(); i++){
            msg(MSG::DEBUG) << "pattern[" << i << "] = " << ((TObjString*) m_TObjArrTrackPathPatterns.At(i))->GetString().Data() << endmsg;
    }
 
  
  // ... declare histograms for some analysis
  InitHistos();
  
  
  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.

◆ 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, 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.

◆ Truth_info()

StatusCode ALFA_GloRec::Truth_info ( )

private

Definition at line 269 of file ALFA_GloRec.cxx.

                                  {
 
  const McEventCollection* mcTru = nullptr;      
  StatusCode sc = evtStore()->retrieve(mcTru,m_TruthCollectionName);
  if(sc.isFailure() || !mcTru){
    msg(MSG::DEBUG) << "Container "<< m_TruthCollectionName <<" NOT FOUND !!!!!!!" << endmsg;
    return StatusCode::FAILURE;
  }
 
  
  //int VertexFoundInFrontOfALFA = 0;
 
 
  McEventCollection::const_iterator mcTruBeg = mcTru->begin();
  McEventCollection::const_iterator mcTruEnd = mcTru->end();  
 
  //loop over collection (container) with hits (i.e. over 1 event)
  for(;mcTruBeg!=mcTruEnd;++mcTruBeg){
   
    //loop over one event
    m_z_g_pos = -1.e10;
    m_z_g_neg =  1.e10;
    for(const HepMC::ConstGenParticlePtr& begGen:  (**mcTruBeg)){
            if(((std::abs(begGen->pdg_id())==2212)||(std::abs(begGen->pdg_id())==211))){
                double px, py, pz;
                double x, y, z;
                px = begGen->momentum().px();   
                py = begGen->momentum().py();
                pz = begGen->momentum().pz();
                                HepMC::ConstGenVertexPtr l_prodvert = begGen->production_vertex();
                x = l_prodvert->position().x();
                y = l_prodvert->position().y();
                z = l_prodvert->position().z();
                if(z > 236000.){
                    m_px_g_pos = px;
                    m_py_g_pos = py;
                    m_pz_g_pos = pz;
                    m_x_g_pos = x;
                    m_y_g_pos = y;
                    m_z_g_pos = z;                                        
                    continue;
                }else if(z < -236000.){
                    m_px_g_neg = px;
                    m_py_g_neg = py;
                    m_pz_g_neg = pz;
                    m_x_g_neg = x;
                    m_y_g_neg = y;
                    m_z_g_neg = z;            
                    continue;                           
                }
        
              // ... break the loop if both suitable vertices were found
              if (m_z_g_pos > 0. && m_z_g_neg < 0.){
                msg(MSG::DEBUG) << "gen px, py, pz = " << m_px_g_pos <<  "  " << m_py_g_pos << "  " << m_pz_g_pos << endmsg;
                    msg(MSG::DEBUG) << "prod x, y, z = " << m_x_g_pos <<  "  " << m_y_g_pos << "  " << m_z_g_pos << endmsg;
                    msg(MSG::DEBUG) << "gen px, py, pz = " << m_px_g_neg <<  "  " << m_py_g_neg << "  " << m_pz_g_neg << endmsg;
                    msg(MSG::DEBUG) << "prod x, y, z = " << m_x_g_neg <<  "  " << m_y_g_neg << "  " << m_z_g_neg << endmsg;
                    break;
                }
      }
    }
  }  
  return StatusCode::SUCCESS;
}

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

◆ WriteHistos()

void ALFA_GloRec::WriteHistos ( )

private

Definition at line 469 of file ALFA_GloRec.cxx.

                             {
  // ... write the histograms into a root file
    TFile fw_glorecanalysis(m_strGloRecAnalysisFile.data(), "RECREATE", m_strGloRecAnalysisFile.data(), 0);
  m_th1_xnearuppotresiduals -> Write();
  m_th1_ynearuppotresiduals -> Write();
  m_th1_xfaruppotresiduals -> Write();
  m_th1_yfaruppotresiduals -> Write();
  m_th1_xnearlwpotresiduals -> Write();
  m_th1_ynearlwpotresiduals -> Write();
  m_th1_xfarlwpotresiduals -> Write();
  m_th1_yfarlwpotresiduals -> Write();
  m_th2_truexvsrecx -> Write();    
    m_th2_trueyvsrecy -> Write();    
    m_th2_truexslopevsrecxslope -> Write();    
    m_th2_trueyslopevsrecyslope -> Write();    
    m_th1_recxovertruex -> Write();    
    m_th1_recyovertruey -> Write();    
    m_th1_recxslopeovertruexslope -> Write();    
    m_th1_recyslopeovertrueyslope -> Write();    
    m_th1_recxminustruex -> Write();    
    m_th1_recyminustruey -> Write();    
    m_th1_recxslopeminustruexslope -> Write();    
    m_th1_recyslopeminustrueyslope -> Write();    
  m_th2_extrapxvsrecxnearpot -> Write();   
  m_th2_extrapyvsrecynearpot -> Write();
  m_th1_recxoverextrapxnearpot -> Write();
  m_th1_recyoverextrapynearpot -> Write();
  m_th1_recxminusextrapxnearpot -> Write();
  m_th1_recyminusextrapynearpot -> Write();
  m_th2_extrapxvsrecxfarpot -> Write();   
  m_th2_extrapyvsrecyfarpot -> Write();
  m_th1_recxoverextrapxfarpot -> Write();
  m_th1_recyoverextrapyfarpot -> Write();
  m_th1_recxminusextrapxfarpot -> Write();
  m_th1_recyminusextrapyfarpot -> Write();
 
  fw_glorecanalysis.Write();
  fw_glorecanalysis.Close();
}

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_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_iDataType

int ALFA_GloRec::m_iDataType = 0

private

Definition at line 47 of file ALFA_GloRec.h.

◆ m_pGloRecEvCollection

ALFA_GloRecEvCollection* ALFA_GloRec::m_pGloRecEvCollection {}

private

Definition at line 41 of file ALFA_GloRec.h.

41{};

◆ m_px_g_neg

double ALFA_GloRec::m_px_g_neg {}

private

Definition at line 78 of file ALFA_GloRec.h.

78{};

◆ m_px_g_pos

double ALFA_GloRec::m_px_g_pos {}

private

Definition at line 71 of file ALFA_GloRec.h.

71{};

◆ m_py_g_neg

double ALFA_GloRec::m_py_g_neg {}

private

Definition at line 79 of file ALFA_GloRec.h.

79{};

◆ m_py_g_pos

double ALFA_GloRec::m_py_g_pos {}

private

Definition at line 72 of file ALFA_GloRec.h.

72{};

◆ m_pz_g_neg

double ALFA_GloRec::m_pz_g_neg {}

private

Definition at line 80 of file ALFA_GloRec.h.

80{};

◆ m_pz_g_pos

double ALFA_GloRec::m_pz_g_pos {}

private

Definition at line 73 of file ALFA_GloRec.h.

73{};

◆ m_strGloRecAnalysisFile

std::string ALFA_GloRec::m_strGloRecAnalysisFile

private

Definition at line 50 of file ALFA_GloRec.h.

◆ m_strGloRecCollectionName

std::string ALFA_GloRec::m_strGloRecCollectionName

private

Definition at line 51 of file ALFA_GloRec.h.

◆ m_strLocRecCorrCollectionName

std::string ALFA_GloRec::m_strLocRecCorrCollectionName

private

Definition at line 52 of file ALFA_GloRec.h.

◆ m_strTrackPathPatterns

std::string ALFA_GloRec::m_strTrackPathPatterns

private

Definition at line 48 of file ALFA_GloRec.h.

◆ m_th1_recxminusextrapxfarpot

TH1F* ALFA_GloRec::m_th1_recxminusextrapxfarpot {}

private

Definition at line 126 of file ALFA_GloRec.h.

126{};

◆ m_th1_recxminusextrapxnearpot

TH1F* ALFA_GloRec::m_th1_recxminusextrapxnearpot {}

private

Definition at line 119 of file ALFA_GloRec.h.

119{};

◆ m_th1_recxminustruex

TH1F* ALFA_GloRec::m_th1_recxminustruex {}

private

Definition at line 110 of file ALFA_GloRec.h.

110{};

◆ m_th1_recxoverextrapxfarpot

TH1F* ALFA_GloRec::m_th1_recxoverextrapxfarpot {}

private

Definition at line 124 of file ALFA_GloRec.h.

124{};

◆ m_th1_recxoverextrapxnearpot

TH1F* ALFA_GloRec::m_th1_recxoverextrapxnearpot {}

private

Definition at line 117 of file ALFA_GloRec.h.

117{};

◆ m_th1_recxovertruex

TH1F* ALFA_GloRec::m_th1_recxovertruex {}

private

Definition at line 106 of file ALFA_GloRec.h.

106{};

◆ m_th1_recxslopeminustruexslope

TH1F* ALFA_GloRec::m_th1_recxslopeminustruexslope {}

private

Definition at line 112 of file ALFA_GloRec.h.

112{};

◆ m_th1_recxslopeovertruexslope

TH1F* ALFA_GloRec::m_th1_recxslopeovertruexslope {}

private

Definition at line 108 of file ALFA_GloRec.h.

108{};

◆ m_th1_recyminusextrapyfarpot

TH1F* ALFA_GloRec::m_th1_recyminusextrapyfarpot {}

private

Definition at line 127 of file ALFA_GloRec.h.

127{};

◆ m_th1_recyminusextrapynearpot

TH1F* ALFA_GloRec::m_th1_recyminusextrapynearpot {}

private

Definition at line 120 of file ALFA_GloRec.h.

120{};

◆ m_th1_recyminustruey

TH1F* ALFA_GloRec::m_th1_recyminustruey {}

private

Definition at line 111 of file ALFA_GloRec.h.

111{};

◆ m_th1_recyoverextrapyfarpot

TH1F* ALFA_GloRec::m_th1_recyoverextrapyfarpot {}

private

Definition at line 125 of file ALFA_GloRec.h.

125{};

◆ m_th1_recyoverextrapynearpot

TH1F* ALFA_GloRec::m_th1_recyoverextrapynearpot {}

private

Definition at line 118 of file ALFA_GloRec.h.

118{};

◆ m_th1_recyovertruey

TH1F* ALFA_GloRec::m_th1_recyovertruey {}

private

Definition at line 107 of file ALFA_GloRec.h.

107{};

◆ m_th1_recyslopeminustrueyslope

TH1F* ALFA_GloRec::m_th1_recyslopeminustrueyslope {}

private

Definition at line 113 of file ALFA_GloRec.h.

113{};

◆ m_th1_recyslopeovertrueyslope

TH1F* ALFA_GloRec::m_th1_recyslopeovertrueyslope {}

private

Definition at line 109 of file ALFA_GloRec.h.

109{};

◆ m_th1_x_g

TH1F* ALFA_GloRec::m_th1_x_g {}

private

Definition at line 88 of file ALFA_GloRec.h.

88{};

◆ m_th1_xfarlwpotresiduals

TH1F* ALFA_GloRec::m_th1_xfarlwpotresiduals {}

private

Definition at line 99 of file ALFA_GloRec.h.

99{};

◆ m_th1_xfaruppotresiduals

TH1F* ALFA_GloRec::m_th1_xfaruppotresiduals {}

private

Definition at line 95 of file ALFA_GloRec.h.

95{};

◆ m_th1_xnearlwpotresiduals

TH1F* ALFA_GloRec::m_th1_xnearlwpotresiduals {}

private

Definition at line 97 of file ALFA_GloRec.h.

97{};

◆ m_th1_xnearuppotresiduals

TH1F* ALFA_GloRec::m_th1_xnearuppotresiduals {}

private

Definition at line 93 of file ALFA_GloRec.h.

93{};

◆ m_th1_xslope_g

TH1F* ALFA_GloRec::m_th1_xslope_g {}

private

Definition at line 90 of file ALFA_GloRec.h.

90{};

◆ m_th1_y_g

TH1F* ALFA_GloRec::m_th1_y_g {}

private

Definition at line 89 of file ALFA_GloRec.h.

89{};

◆ m_th1_yfarlwpotresiduals

TH1F* ALFA_GloRec::m_th1_yfarlwpotresiduals {}

private

Definition at line 100 of file ALFA_GloRec.h.

100{};

◆ m_th1_yfaruppotresiduals

TH1F* ALFA_GloRec::m_th1_yfaruppotresiduals {}

private

Definition at line 96 of file ALFA_GloRec.h.

96{};

◆ m_th1_ynearlwpotresiduals

TH1F* ALFA_GloRec::m_th1_ynearlwpotresiduals {}

private

Definition at line 98 of file ALFA_GloRec.h.

98{};

◆ m_th1_ynearuppotresiduals

TH1F* ALFA_GloRec::m_th1_ynearuppotresiduals {}

private

Definition at line 94 of file ALFA_GloRec.h.

94{};

◆ m_th1_yslope_g

TH1F* ALFA_GloRec::m_th1_yslope_g {}

private

Definition at line 91 of file ALFA_GloRec.h.

91{};

◆ m_th2_extrapxvsrecxfarpot

TH2F* ALFA_GloRec::m_th2_extrapxvsrecxfarpot {}

private

Definition at line 122 of file ALFA_GloRec.h.

122{};

◆ m_th2_extrapxvsrecxnearpot

TH2F* ALFA_GloRec::m_th2_extrapxvsrecxnearpot {}

private

Definition at line 115 of file ALFA_GloRec.h.

115{};

◆ m_th2_extrapyvsrecyfarpot

TH2F* ALFA_GloRec::m_th2_extrapyvsrecyfarpot {}

private

Definition at line 123 of file ALFA_GloRec.h.

123{};

◆ m_th2_extrapyvsrecynearpot

TH2F* ALFA_GloRec::m_th2_extrapyvsrecynearpot {}

private

Definition at line 116 of file ALFA_GloRec.h.

116{};

◆ m_th2_truexslopevsrecxslope

TH2F* ALFA_GloRec::m_th2_truexslopevsrecxslope {}

private

Definition at line 104 of file ALFA_GloRec.h.

104{};

◆ m_th2_truexvsrecx

TH2F* ALFA_GloRec::m_th2_truexvsrecx {}

private

Definition at line 102 of file ALFA_GloRec.h.

102{};

◆ m_th2_trueyslopevsrecyslope

TH2F* ALFA_GloRec::m_th2_trueyslopevsrecyslope {}

private

Definition at line 105 of file ALFA_GloRec.h.

105{};

◆ m_th2_trueyvsrecy

TH2F* ALFA_GloRec::m_th2_trueyvsrecy {}

private

Definition at line 103 of file ALFA_GloRec.h.

103{};

◆ m_TObjArrTrackPathPatterns

TObjArray ALFA_GloRec::m_TObjArrTrackPathPatterns {}

private

Definition at line 49 of file ALFA_GloRec.h.

49{};

◆ m_TruthCollectionName

std::string ALFA_GloRec::m_TruthCollectionName

private

Definition at line 53 of file ALFA_GloRec.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.

◆ m_x_g_neg

double ALFA_GloRec::m_x_g_neg {}

private

Definition at line 81 of file ALFA_GloRec.h.

81{};

◆ m_x_g_pos

double ALFA_GloRec::m_x_g_pos {}

private

Definition at line 74 of file ALFA_GloRec.h.

74{};

◆ m_y_g_neg

double ALFA_GloRec::m_y_g_neg {}

private

Definition at line 82 of file ALFA_GloRec.h.

82{};

◆ m_y_g_pos

double ALFA_GloRec::m_y_g_pos {}

private

Definition at line 75 of file ALFA_GloRec.h.

75{};

◆ m_z_g_neg

double ALFA_GloRec::m_z_g_neg {}

private

Definition at line 83 of file ALFA_GloRec.h.

83{};

◆ m_z_g_pos

double ALFA_GloRec::m_z_g_pos {}

private

Definition at line 76 of file ALFA_GloRec.h.

76{};

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

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ ALFA_GloRec()

◆ ~ALFA_GloRec()

Member Function Documentation

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ FillTrackCandHistos()

◆ finalize()

◆ InitHistos()

◆ initialize()

◆ inputHandles()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ Truth_info()

◆ updateVHKA()

◆ WriteHistos()

Member Data Documentation

◆ m_detStore

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_iDataType

◆ m_pGloRecEvCollection

◆ m_px_g_neg

◆ m_px_g_pos

◆ m_py_g_neg

◆ m_py_g_pos

◆ m_pz_g_neg

◆ m_pz_g_pos

◆ m_strGloRecAnalysisFile

◆ m_strGloRecCollectionName

◆ m_strLocRecCorrCollectionName

◆ m_strTrackPathPatterns

◆ m_th1_recxminusextrapxfarpot

◆ m_th1_recxminusextrapxnearpot

◆ m_th1_recxminustruex

◆ m_th1_recxoverextrapxfarpot

◆ m_th1_recxoverextrapxnearpot

◆ m_th1_recxovertruex

◆ m_th1_recxslopeminustruexslope

◆ m_th1_recxslopeovertruexslope

◆ m_th1_recyminusextrapyfarpot

◆ m_th1_recyminusextrapynearpot

◆ m_th1_recyminustruey

◆ m_th1_recyoverextrapyfarpot

◆ m_th1_recyoverextrapynearpot

◆ m_th1_recyovertruey

◆ m_th1_recyslopeminustrueyslope

◆ m_th1_recyslopeovertrueyslope

◆ m_th1_x_g

◆ m_th1_xfarlwpotresiduals

◆ m_th1_xfaruppotresiduals

◆ m_th1_xnearlwpotresiduals

◆ m_th1_xnearuppotresiduals

◆ m_th1_xslope_g

◆ m_th1_y_g

◆ m_th1_yfarlwpotresiduals

◆ m_th1_yfaruppotresiduals

◆ m_th1_ynearlwpotresiduals

◆ m_th1_ynearuppotresiduals

◆ m_th1_yslope_g

◆ m_th2_extrapxvsrecxfarpot

◆ m_th2_extrapxvsrecxnearpot