#include <QratCscClusterFitter.h>

Inheritance diagram for QratCscClusterFitter:

Collaboration diagram for QratCscClusterFitter:

Public Types
typedef std::vector< Muon::CscStripPrepData * >	StripList

typedef ICscStripFitter::Result	StripFit

typedef std::vector< StripFit >	StripFitList

typedef std::vector< double >	ChargeList

typedef std::vector< std::string >	DataNames

typedef std::map< std::string, double >	DataMap

typedef std::vector< Result >	Results

Public Member Functions
	QratCscClusterFitter (const std::string &, const std::string &, const IInterface *)

	~QratCscClusterFitter ()=default

StatusCode	initialize ()

const DataNames &	dataNames () const

Results	fit (const StripFitList &sfits) const

Results	fit (const StripFitList &sfits, double dposdz) const

double	getCorrectedError (const Muon::CscPrepData *pclu, double slope) const

virtual Results	fit (const StripFitList &sfits) const=0

virtual Results	fit (const StripFitList &sfits, double dposdz) const=0

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
std::vector< unsigned int >	m_max_width

std::string	m_posopt_eta

std::string	m_posopt_phi

std::string	m_erropt_eta

std::string	m_erropt_phi

double	m_error_eta

double	m_error_phi

double	m_precisionErrorScaler

double	m_qrat_maxdiff

double	m_qrat_maxsig

double	m_error_tantheta

double	m_xtan_css_eta_offset

double	m_xtan_css_eta_slope

double	m_xtan_csl_eta_offset

double	m_xtan_csl_eta_slope

double	m_qratmin_css_eta

double	m_qratmin_csl_eta

std::vector< double >	m_qratcor_css_eta

std::vector< double >	m_qratcor_csl_eta

double	m_atanh_a_css_eta

double	m_atanh_b_css_eta

double	m_atanh_c_css_eta

double	m_atanh_x0_css_eta

double	m_atanh_a_csl_eta

double	m_atanh_b_csl_eta

double	m_atanh_c_csl_eta

double	m_atanh_x0_csl_eta

double	m_dposmin

SG::ReadCondHandleKey< MuonGM::MuonDetectorManager >	m_DetectorManagerKey
	retrieve MuonDetectorManager from the conditions store More...

ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc

ToolHandle< ICscAlignmentTool >	m_alignmentTool

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 26 of file QratCscClusterFitter.h.

Member Typedef Documentation

◆ ChargeList

typedef std::vector<double> ICscClusterFitter::ChargeList

inherited

Definition at line 43 of file ICscClusterFitter.h.

◆ DataMap

typedef std::map<std::string, double> ICscClusterFitter::DataMap

inherited

Definition at line 47 of file ICscClusterFitter.h.

◆ DataNames

typedef std::vector<std::string> ICscClusterFitter::DataNames

inherited

Definition at line 46 of file ICscClusterFitter.h.

◆ Results

typedef std::vector<Result> ICscClusterFitter::Results

inherited

Definition at line 101 of file ICscClusterFitter.h.

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ StripFit

typedef ICscStripFitter::Result ICscClusterFitter::StripFit

inherited

Definition at line 39 of file ICscClusterFitter.h.

◆ StripFitList

typedef std::vector<StripFit> ICscClusterFitter::StripFitList

inherited

Definition at line 40 of file ICscClusterFitter.h.

◆ StripList

typedef std::vector<Muon::CscStripPrepData*> ICscClusterFitter::StripList

inherited

Definition at line 36 of file ICscClusterFitter.h.

Constructor & Destructor Documentation

◆ QratCscClusterFitter()

QratCscClusterFitter::QratCscClusterFitter	(	const std::string &	type,
		const std::string &	aname,
		const IInterface *	parent
	)

Definition at line 210 of file QratCscClusterFitter.cxx.

                                                                                                                   :
     AthAlgTool(type, aname, parent) {
     declareInterface<ICscClusterFitter>(this);
     m_max_width.push_back(5);                   // CSS eta
     m_max_width.push_back(5);                   // CSL eta
     m_max_width.push_back(3);                   // CSS phi
     m_max_width.push_back(2);                   // CSL phi
     declareProperty("max_width", m_max_width);  // Maximum width (strips) for unspoiled clusters
     declareProperty("position_option_eta", m_posopt_eta = "ATANH");
     declareProperty("position_option_phi", m_posopt_phi = "NONE");
     declareProperty("error_option_eta", m_erropt_eta = "CHARGE");
     declareProperty("error_option_phi", m_erropt_phi = "NONE");
     declareProperty("error_eta", m_error_eta = 0.050);  // in mm this is for fixed error
     declareProperty("error_phi", m_error_phi = 0.140);  // in mm this is for fixed error
  
     declareProperty("precisionErrorScaler", m_precisionErrorScaler = 1.0);  // in mm this is for fixed error
  
     declareProperty("qrat_maxdiff", m_qrat_maxdiff = 0.15);                   // in strips
     declareProperty("qrat_maxsig", m_qrat_maxsig = 6.0);                      // in strips
     declareProperty("error_tantheta", m_error_tantheta = 0.57);               // in mm
     declareProperty("xtan_css_eta_offset", m_xtan_css_eta_offset = 0.0015);   // in mm
     declareProperty("xtan_css_eta_slope", m_xtan_css_eta_slope = 0.000137);   // in mm
     declareProperty("xtan_csl_eta_offset", m_xtan_csl_eta_offset = -0.0045);  // in mm
     declareProperty("xtan_csl_eta_slope", m_xtan_csl_eta_slope = 0.000131);   // in mm
     declareProperty("qratmin_css_eta", m_qratmin_css_eta = 0.0940459);
     declareProperty("qratmin_csl_eta", m_qratmin_csl_eta = 0.108975);
     declareProperty("qratcor_css_eta", m_qratcor_css_eta);  // in strips
     declareProperty("qratcor_csl_eta", m_qratcor_csl_eta);  // in strips
  
     declareProperty("atanh_a_css_eta", m_atanh_a_css_eta = 1.5);
     declareProperty("atanh_b_css_eta", m_atanh_b_css_eta = 1.411);
     declareProperty("atanh_c_css_eta", m_atanh_c_css_eta = 2.329);
     declareProperty("atanh_x0_css_eta", m_atanh_x0_css_eta = 0.6601);
     declareProperty("atanh_a_csl_eta", m_atanh_a_csl_eta = 1.516);
     declareProperty("atanh_b_csl_eta", m_atanh_b_csl_eta = 1.427);
     declareProperty("atanh_c_csl_eta", m_atanh_c_csl_eta = 2.35);
     declareProperty("atanh_x0_csl_eta", m_atanh_x0_csl_eta = 0.6615);
  
     declareProperty("dposmin", m_dposmin = 0.082);
 }

◆ ~QratCscClusterFitter()

QratCscClusterFitter::~QratCscClusterFitter ( )

default

Member Function Documentation

◆ dataNames()

const DataNames & QratCscClusterFitter::dataNames ( ) const

virtual

Reimplemented from ICscClusterFitter.

Definition at line 302 of file QratCscClusterFitter.cxx.

                                                        {
     auto init = [&]() {
         DataNames dnames;
         bool dofixed = false;
         bool docharge = false;
         if (m_posopt_phi == "POLYNOMIAL" || m_posopt_phi == "TABLE" || m_posopt_phi == "ATANH") {
             if (m_erropt_phi == "FIXED") dofixed = true;
             if (m_erropt_phi == "CHARGE") docharge = true;
         }
         if (m_posopt_eta == "POLYNOMIAL" || m_posopt_eta == "TABLE" || m_posopt_eta == "ATANH") {
             if (m_erropt_eta == "FIXED") dofixed = true;
             if (m_erropt_eta == "CHARGE") docharge = true;
         }
         if (dofixed || docharge) {
             dnames.emplace_back("scor1");
             dnames.emplace_back("scor2");
             dnames.emplace_back("scor");
         }
         if (docharge) {
             dnames.emplace_back("dscor1");
             dnames.emplace_back("dscor2");
             dnames.emplace_back("dscor");
             dnames.emplace_back("scordiff");
             dnames.emplace_back("dscordiff");
         }
         return dnames;
     };
     static const DataNames dnames = init();
     return dnames;
 }

◆ declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

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

inlineprivateinherited

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

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 245 of file AthCommonDataStore.h.

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

◆ declareProperty() [2/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 221 of file AthCommonDataStore.h.

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

◆ declareProperty() [3/6]

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

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

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

◆ declareProperty() [4/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 333 of file AthCommonDataStore.h.

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

◆ declareProperty() [5/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

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

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

Definition at line 352 of file AthCommonDataStore.h.

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

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                  {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

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

◆ evtStore() [1/2]

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

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::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

◆ fit() [1/4]

Results QratCscClusterFitter::fit ( const StripFitList & sfits ) const

virtual

Implements ICscClusterFitter.

Definition at line 738 of file QratCscClusterFitter.cxx.

                                                                  {
     Results results = fit(sfits, 0.0);
     Results new_results;
     for (unsigned int iresult = 0; iresult < results.size(); ++iresult) {
         Result res = results[iresult];
         if (res.fitStatus) {
             new_results.push_back(res);
             continue;
         }
         // Fetch the chamber type.
         const CscStripPrepData* pstrip = sfits[0].strip;
         Identifier idStrip0 = pstrip->identify();
         int station = m_idHelperSvc->cscIdHelper().stationName(idStrip0) - 49;  // 1=CSS, 2=CSL
         // Calculate the angle of incidence.
         double tantht = 0.0;
         double pos = res.position;
         if (station == 1) {
             tantht = m_xtan_css_eta_offset + m_xtan_css_eta_slope * pos;
         } else {
             tantht = m_xtan_csl_eta_offset + m_xtan_csl_eta_slope * pos;
         }
         // Correct the error using this angle.
         double dpostht = m_error_tantheta * std::abs(tantht);
         double dpos = res.dposition;
  
         res.dposition = sqrt(dpos * dpos + dpostht * dpostht);
  
         // Return the updated result.
         new_results.push_back(res);
     }
  
     return new_results;
 }

◆ fit() [2/4]

virtual Results ICscClusterFitter::fit

◆ fit() [3/4]

Results QratCscClusterFitter::fit	(	const StripFitList &	sfits,
		double	dposdz
	)		const

virtual

Implements ICscClusterFitter.

Definition at line 335 of file QratCscClusterFitter.cxx.

                                                                                   {
     ATH_MSG_VERBOSE("QRAT fit with tool " << name());
  
     Results results;
  
     // Check input has at least three strips.
     unsigned int nstrip = sfits.size();
     if (nstrip < 3) {
         ATH_MSG_VERBOSE("  Input has fewer than three strips.");
         if (nstrip == 2) {
             Muon::CscTimeStatus tstatus = (sfits[0].charge > sfits[1].charge) ? sfits[0].timeStatus : sfits[1].timeStatus;
             results.emplace_back(1, Muon::CscStatusNarrow, tstatus);
         } else if (nstrip == 1) {
             Muon::CscTimeStatus tstatus = sfits[0].timeStatus;
             results.emplace_back(1, Muon::CscStatusNarrow, tstatus);
         }
         return results;
     }
  
     // Fetch the number of strips and check the input arrays.
     for (unsigned int istrip = 0; istrip < nstrip; ++istrip) {
         if (sfits[istrip].strip == nullptr) {
             ATH_MSG_WARNING("Strip pointer is null.");
             results.emplace_back(2);
             return results;
         }
     }
  
     // Use the first strip to extract the layer parameters.
     const CscStripPrepData* pstrip = sfits[0].strip;
     Identifier idStrip0 = pstrip->identify();
  
     // retrieve MuonDetectorManager from the conditions store
     SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
     const MuonGM::MuonDetectorManager* MuonDetMgr = DetectorManagerHandle.cptr();
     if (MuonDetMgr == nullptr) {
         ATH_MSG_ERROR("Null pointer to the MuonDetectorManager conditions object");
         return results;
     }
     const CscReadoutElement* pro = MuonDetMgr->getCscReadoutElement(idStrip0);
  
     bool measphi = m_idHelperSvc->cscIdHelper().CscIdHelper::measuresPhi(idStrip0);
     double pitch = pro->cathodeReadoutPitch(0, measphi);
     unsigned int maxstrip = pro->maxNumberOfStrips(measphi);
     unsigned int strip0 = m_idHelperSvc->cscIdHelper().strip(idStrip0) - 1;
     int station = m_idHelperSvc->cscIdHelper().stationName(idStrip0) - 49;  // 1=CSS, 2=CSL
  
     CscPlane plane = findPlane(station, measphi);
     if (plane == UNKNOWN_PLANE) {
         ATH_MSG_WARNING("Invalid CSC plane: station=" << station << "; measphi=" << measphi);
         results.emplace_back(3);
         return results;
     }
  
     // Display input strips.
     ATH_MSG_VERBOSE("QRAT fittter input has " << nstrip << " strips");
     for (unsigned int istrip = 0; istrip < nstrip; ++istrip) {
         Identifier id = sfits[istrip].strip->identify();
         ATH_MSG_VERBOSE("  " << station << " : " << measphi << "  " << m_idHelperSvc->cscIdHelper().wireLayer(id) << "  " << istrip << " "
                              << m_idHelperSvc->cscIdHelper().strip(id) << " " << sfits[istrip].charge);
     }
  
     // Find the peak strip and check the shape.
     unsigned int istrip_peak = 0;  // strip number within cluster
     // Loop over strips excluding the edges.
     for (unsigned int istrip = 1; istrip < nstrip - 1; ++istrip) {
         // Tell clang to optimize assuming that FP operations may trap.
         CXXUTILS_TRAPPING_FP;
         StripFit sfit = sfits[istrip];
         float qthis = sfit.charge;
         float qlast = sfits[istrip - 1].charge;
         float qnext = sfits[istrip + 1].charge;
         // Peak if the adjacent strips have less charge.
         bool ispeak = qthis > qlast && qthis > qnext;
         // Special case: next strip has the same charge.
         // Require the previous strip has less charge and the next following
         // strip be absent or have less charge.
         if (!ispeak) {
             if (qthis == qnext) { ispeak = (qthis > qlast) && (istrip + 2 == nstrip || sfits[istrip + 2].charge < qthis); }
         }
         // Special case: first and second strips have the same charge.
         // Require the third strip has less charge.
         if (!ispeak) {
             if (istrip == 1) {
                 if (qthis == qlast) {
                     ispeak = qthis > qnext;  // bug found 10/13/07
                 }
             }
         }
         // Record if peak.
         if (ispeak) {
             if (istrip_peak) {                                           // Error if multiple peaks are found.
                 results.emplace_back(6, Muon::CscStatusMultiPeak);  // Time status should be defined in SimpleClusterFit...
                 return results;
             }
             istrip_peak = istrip;
         }
     }
     ATH_MSG_VERBOSE("  Peak is at " << istrip_peak << "th strip in cluster");
  
     // Check we are not on the edge.
     if (strip0 <= 0 || strip0 + nstrip > maxstrip) {
         results.emplace_back(4, Muon::CscStatusEdge, sfits[istrip_peak].timeStatus);
         return results;
     }
  
     // MS: remove this check since width is amplitude dependent.
     // use saturation check instead
     /**************************************
     if ( nstrip > m_max_width[plane] ) {
       //  if ( nstrip_threshold > m_max_width[plane] ) {
       results.push_back(Result(5, Muon::CscStatusWide, sfits[istrip_peak].timeStatus));
       return results;
     }
     ***************************************/
  
     // Cluster is spoiled if peak is not at the center.
     bool is_even = 2 * (nstrip / 2) == nstrip;
     bool atcenter = istrip_peak == nstrip / 2 || (is_even && istrip_peak + 1 == nstrip / 2);
     if (!atcenter) {
         results.emplace_back(7, Muon::CscStatusSkewed, sfits[istrip_peak].timeStatus);
         return results;
     }
  
     if (sfits[istrip_peak].stripStatus == Muon::CscStrStatSaturated || sfits[istrip_peak - 1].stripStatus == Muon::CscStrStatSaturated ||
         sfits[istrip_peak + 1].stripStatus == Muon::CscStrStatSaturated) {
         results.emplace_back(15, Muon::CscStatusSaturated, sfits[istrip_peak].timeStatus);
         return results;
     }
  
     // left/peak/right strip should have good time information....
     if (sfits[istrip_peak].stripStatus != Muon::CscStrStatSuccess
         //       || sfits[istrip_peak-1].stripStatus != Muon::CscStrStatSuccess
         //       || sfits[istrip_peak+1].stripStatus != Muon::CscStrStatSuccess ) {
     ) {
         results.emplace_back(14, Muon::CscStatusStripFitFailed, sfits[istrip_peak].timeStatus);
         return results;
     } else if (sfits[istrip_peak - 1].stripStatus == Muon::CscStrStatHot || sfits[istrip_peak - 1].stripStatus == Muon::CscStrStatDead ||
                sfits[istrip_peak + 1].stripStatus == Muon::CscStrStatHot || sfits[istrip_peak + 1].stripStatus == Muon::CscStrStatDead) {
         results.emplace_back(14, Muon::CscStatusStripFitFailed, sfits[istrip_peak].timeStatus);
         return results;
     }
  
     // Set initial strip position.
     double savg = istrip_peak;
  
     // Calculate QRAT correction to strip position.
     std::string posopt = m_posopt_eta;
     std::string erropt = m_erropt_eta;
     double dpos = 0.0;
     if (measphi) {
         posopt = m_posopt_phi;
         erropt = m_erropt_phi;
     }
     double q1 = sfits[istrip_peak - 1].charge;
     double q0 = sfits[istrip_peak].charge;
     double q2 = sfits[istrip_peak + 1].charge;
     double qrat1 = q1 / q0;
     double qrat2 = q2 / q0;
     double dq1 = sfits[istrip_peak - 1].dcharge;
     double dq0 = sfits[istrip_peak].dcharge;
     double dq2 = sfits[istrip_peak + 1].dcharge;
  
     ATH_MSG_VERBOSE("  QRAT charge ratios: " << qrat1 << " " << qrat2);
     double scor1 = 0.;  // left side correction.
     double dscordqrat1 = 0.;
     double scor2 = 0;  // right side correction.
     double dscordqrat2 = 0.;
     int stat1 = 0;
     int stat2 = 0;
     if (posopt == "POLYNOMIAL") {
         stat1 = qrat_correction(plane, qrat1, scor1, dscordqrat1);
         stat2 = qrat_correction(plane, qrat2, scor2, dscordqrat2);
     } else if (posopt == "TABLE") {
         double qrmin = 0.0;
         const std::vector<double>* pcor = nullptr;
         if (plane == CSS_ETA) {
             qrmin = m_qratmin_css_eta;
             pcor = &m_qratcor_css_eta;
         } else if (plane == CSL_ETA) {
             qrmin = m_qratmin_csl_eta;
             pcor = &m_qratcor_csl_eta;
         } else {
             ATH_MSG_WARNING("    Invalid QRAT plane: " << splane(plane));
             results.emplace_back(8);
             return results;
         }
         if (pcor) {
             stat1 = qrat_interpolation(qrmin, *pcor, qrat1, scor1, dscordqrat1);
             stat2 = qrat_interpolation(qrmin, *pcor, qrat2, scor2, dscordqrat2);
         } else {
             stat1 = 99;
         }
     } else if (posopt == "ATANH") {  // MS: atanh parametrization
         double a, b, c, x0;          // parameters
         if (plane == CSS_ETA) {
             a = m_atanh_a_css_eta;
             b = m_atanh_b_css_eta;
             c = m_atanh_c_css_eta;
             x0 = m_atanh_x0_css_eta;
         } else if (plane == CSL_ETA) {
             a = m_atanh_a_csl_eta;
             b = m_atanh_b_csl_eta;
             c = m_atanh_c_csl_eta;
             x0 = m_atanh_x0_csl_eta;
         } else {
             ATH_MSG_WARNING("    Invalid QRAT plane: " << splane(plane));
             results.emplace_back(8);
             return results;
         }
         stat1 = qrat_atanh(a, b, c, x0, qrat1, scor1, dscordqrat1);
         stat2 = qrat_atanh(a, b, c, x0, qrat2, scor2, dscordqrat2);
  
     } else {
         ATH_MSG_WARNING("    Invalid position option: " << posopt);
         results.emplace_back(9);
         return results;
     }
     if (stat1 || stat2) {
         ATH_MSG_VERBOSE("    QRAT correction failed: SPOILED");
         results.emplace_back(10);
         return results;
     }
     ATH_MSG_VERBOSE("  QRAT strip corrs: " << scor1 << " " << scor2);
     ATH_MSG_VERBOSE("  QRAT derivatives: " << dscordqrat1 << " " << dscordqrat2);
  
     // Compare left and right.
     // Flip sign of the left side correction.
     scor1 = -scor1;
     dscordqrat1 = -dscordqrat1;
     double scor = 0.0;
     DataMap dmap;
     // Calculation weighting corrections by their derivatives, i.e. assuming the
     // two charge ratios have the same error.
     if (erropt == "FIXED") {
         double w1 = 1 / (dscordqrat1 * dscordqrat1);
         double w2 = 1 / (dscordqrat2 * dscordqrat2);
         scor = (w1 * scor1 + w2 * scor2) / (w1 + w2);
         double scor_diff = std::abs(scor2 - scor1);
         ATH_MSG_VERBOSE("  Combined corr: " << scor);
         dpos = measphi ? m_error_phi : m_error_eta;
         // Fill data map.
         dmap["scor1"] = scor1;
         dmap["scor2"] = scor2;
         dmap["scor"] = scor;
         // Exit if measurements are inconsistent.
         if (scor_diff > m_qrat_maxdiff) {
             ATH_MSG_VERBOSE("  SPOILED (scor_diff=" << scor_diff << ")");
             results.emplace_back(11, Muon::CscStatusQratInconsistent);
             return results;
         }
         // Calculation using the (independent) errors in the three charges.
     } else if (erropt == "CHARGE") {
         // Calculate intermediate variables.
         double x1 = dscordqrat1 * qrat1;
         double x2 = dscordqrat2 * qrat2;
         double dqq0 = dq0 / q0;
         double dqq1 = dq1 / q1;
         double dqq2 = dq2 / q2;
         double rnum = (-x1 * x1 * dqq1 * dqq1 + x2 * x2 * dqq2 * dqq2 + (x2 * x2 - x1 * x1) * dqq0 * dqq0);
         double rden = (x1 * x1 * dqq1 * dqq1 + x2 * x2 * dqq2 * dqq2 + (x2 - x1) * (x2 - x1) * dqq0 * dqq0);
         double rfac = rnum / rden;
         // Calculate the significance of the difference between the measurements.
         double dscor_diff = sqrt(rden);
         double scor_diff = scor2 - scor1;
         double scor_sig = std::abs(scor_diff) / dscor_diff;
         // Calculate the weighted average of the corrections.
         double w1 = 0.5 * (1.0 + rfac);
         double w2 = 0.5 * (1.0 - rfac);
         scor = w1 * scor1 + w2 * scor2;
         // Calculate the error in this average.
         double ddscor1 = w1 * x1 * dqq1;
         double ddscor2 = w2 * x2 * dqq2;
         double ddscor0 = (w1 * x1 + w2 * x2) * dqq0;
         double dscor = sqrt(ddscor1 * ddscor1 + ddscor2 * ddscor2 + ddscor0 * ddscor0);
         dpos = pitch * dscor;
         // add minimum error (in mm) in quadrature:
         dpos = sqrt(dpos * dpos + m_dposmin * m_dposmin);
         // Fill data map.
         double dscor1 = std::abs(x1) * sqrt(dqq1 * dqq1 + dqq0 * dqq0);
         double dscor2 = std::abs(x2) * sqrt(dqq2 * dqq2 + dqq0 * dqq0);
         dmap["scor1"] = scor1;
         dmap["dscor1"] = dscor1;
         dmap["scor2"] = scor2;
         dmap["dscor2"] = dscor2;
         dmap["scordiff"] = scor_diff;
         dmap["dscordiff"] = dscor_diff;
         dmap["scor"] = scor;
         dmap["dscor"] = dscor;
         // Debugging.
         ATH_MSG_VERBOSE("QRAT CHARGE calculation");
         ATH_MSG_VERBOSE("       q1, q0, q2: " << q1 << " " << q0 << " " << q2);
         ATH_MSG_VERBOSE("    dq1, dq0, dq2: " << dq1 << " " << dq0 << " " << dq2);
         ATH_MSG_VERBOSE("  dscordqrat1, x1: " << dscordqrat1 << " " << x1);
         ATH_MSG_VERBOSE("  dscordqrat2, x2: " << dscordqrat2 << " " << x2);
         ATH_MSG_VERBOSE("     scor1 = " << scor1 << " +/- " << dscor1);
         ATH_MSG_VERBOSE("     scor2 = " << scor2 << " +/- " << dscor2);
         ATH_MSG_VERBOSE("      scor = " << scor << " +/- " << dscor);
         ATH_MSG_VERBOSE("  scordiff = " << scor_diff << " +/- " << dscor_diff);
         ATH_MSG_VERBOSE("  scor_sig = " << scor_sig);
  
         // Exit if measurements are inconsistent.
         if (scor_sig > m_qrat_maxsig) {
             ATH_MSG_VERBOSE("  SPOILED (scor_sig=" << scor_sig << ")");
             results.emplace_back(12, Muon::CscStatusQratInconsistent, sfits[istrip_peak].timeStatus);
             return results;
         }
     } else {
         ATH_MSG_WARNING("    Invalid error option: " << erropt);
         results.emplace_back(13, Muon::CscStatusUndefined, sfits[istrip_peak].timeStatus);
         return results;
     }
  
     // Correct strip position.
     savg += scor + strip0;
     ATH_MSG_VERBOSE("    QRAT corr " << splane(plane) << " nstrip=" << nstrip << " savg=" << savg << " qrat1=" << qrat1
                                      << " qrat2=" << qrat2 << " scor1=" << scor1 << " scor2=" << scor2);
  
     // Error due to incident angle.
     double dpostht = m_error_tantheta * std::abs(tantheta);
  
     // Scale error on dpos
     dpos = dpos * m_precisionErrorScaler;
  
     // Set return values.
     Result res(0, Muon::CscStatusUnspoiled);
     res.strip = istrip_peak;
     res.position = pitch * (savg + 0.5 - 0.5 * maxstrip);
  
     // internal alignment ...
     Identifier id = sfits[res.strip].strip->identify();
     double offset = m_alignmentTool->getAlignmentOffset(id);
     res.position -= offset;
  
     res.dposition = sqrt(dpos * dpos + dpostht * dpostht);
  
     res.fstrip = 0;
     res.lstrip = nstrip - 1;
     res.time = sfits[istrip_peak].time;
     res.time_beforeT0Corr = sfits[istrip_peak].time_beforeT0Corr;
     res.time_beforeBPCorr = sfits[istrip_peak].time_beforeBPCorr;
     res.timeStatus = sfits[istrip_peak].timeStatus;
  
     res.qleft = q1;
     res.qpeak = q0;
     res.qright = q2;
  
     //  res.diff = dmap["scordiff"];
     //  res.sig  = dmap["scordiff"]/dmap["dscordiff"];
  
     // cluster charge should be qsum over three strip... 3/21/2011
     res.charge = res.qleft + res.qpeak + res.qright;
     res.charge_beforeBPCorr =
         sfits[istrip_peak].charge_beforeBPCorr + sfits[istrip_peak - 1].charge_beforeBPCorr + sfits[istrip_peak + 1].charge_beforeBPCorr;
  
     res.dataMap = dmap;
  
     ATH_MSG_VERBOSE("   Position :: pos=" << res.position << " dpos:dtht=" << dpos << ":" << dpostht << " ==>" << res.dposition
                                           << "  at tanth=" << tantheta);
  
     results.push_back(res);
     return results;
 }

◆ fit() [4/4]

virtual Results ICscClusterFitter::fit

◆ getCorrectedError()

double QratCscClusterFitter::getCorrectedError	(	const Muon::CscPrepData *	pclu,
		double	slope
	)		const

virtual

Implements ICscClusterFitter.

Definition at line 705 of file QratCscClusterFitter.cxx.

                                                                                           {
     // Cluster position.
     Trk::ParamDefs icor = Trk::loc1;
     Trk::ParamDefs ierr = Trk::loc1;
     double pos = pclu->localPosition()[icor];
     double dpos = Amg::error(pclu->localCovariance(), ierr);
  
     Identifier idStrip0 = pclu->identify();
     int station = m_idHelperSvc->cscIdHelper().stationName(idStrip0) - 49;  // 1=CSS, 2=CSL
     // Calculate the angle of incidence.
     double tantht = 0.0;
     if (station == 1) {
         tantht = m_xtan_css_eta_offset + m_xtan_css_eta_slope * pos;
     } else {
         tantht = m_xtan_csl_eta_offset + m_xtan_csl_eta_slope * pos;
     }
     // Correct the error using this angle.
     double old_dpostht = m_error_tantheta * std::abs(tantht);
  
     double new_dpostht = m_error_tantheta * std::abs(slope);
  
     double newError = sqrt(dpos * dpos - old_dpostht * old_dpostht + new_dpostht * new_dpostht);
  
     ATH_MSG_VERBOSE("   Position :: pos=" << pos << " dpos:newdpos=" << dpos << " : " << newError << "  " << old_dpostht << "  "
                                           << new_dpostht);
  
     if (slope < -990) { newError = sqrt(dpos * dpos + old_dpostht * old_dpostht); }
  
     return newError;
 }

◆ initialize()

StatusCode QratCscClusterFitter::initialize ( )

Definition at line 253 of file QratCscClusterFitter.cxx.

                                             {
     ATH_MSG_VERBOSE("Initalizing " << name());
  
     // retrieve MuonDetectorManager from the conditions store
     ATH_CHECK(m_DetectorManagerKey.initialize());
  
     ATH_CHECK(m_idHelperSvc.retrieve());
  
     if (m_alignmentTool.retrieve().isFailure()) {
         ATH_MSG_WARNING(name() << ": unable to retrieve cluster fitter " << m_alignmentTool);
     } else {
         ATH_MSG_DEBUG(name() << ": retrieved " << m_alignmentTool);
     }
  
     ATH_MSG_DEBUG("Properties for " << name() << ":");
     ATH_MSG_DEBUG("       Eta position option: " << m_posopt_eta);
     ATH_MSG_DEBUG("       Phi position option: " << m_posopt_phi);
     ATH_MSG_DEBUG("          Eta error option: " << m_erropt_eta);
     ATH_MSG_DEBUG("          Phi error option: " << m_erropt_phi);
     ATH_MSG_DEBUG("        Eta assigned error: " << m_error_eta);
     ATH_MSG_DEBUG("        Phi assigned error: " << m_error_phi);
     ATH_MSG_DEBUG("        Max strip pos diff: " << m_qrat_maxdiff);
     ATH_MSG_DEBUG("         Max strip pos sig: " << m_qrat_maxsig);
     ATH_MSG_DEBUG("    Non-normal error coeff: " << m_error_tantheta);
     ATH_MSG_DEBUG("  CSS eta pos-slope offset: " << m_xtan_css_eta_offset);
     ATH_MSG_DEBUG("   CSS eta pos-slope slope: " << m_xtan_css_eta_slope);
     ATH_MSG_DEBUG("  CSL eta pos-slope offset: " << m_xtan_csl_eta_offset);
     ATH_MSG_DEBUG("   CSL eta pos-slope slope: " << m_xtan_csl_eta_slope);
     ATH_MSG_DEBUG("      CSS eta table offset: " << m_qratmin_css_eta);
     ATH_MSG_DEBUG("      CSL eta table offset: " << m_qratmin_csl_eta);
     ATH_MSG_DEBUG("        CSS eta table size: " << m_qratcor_css_eta.size());
     ATH_MSG_DEBUG("        CSL eta table size: " << m_qratcor_csl_eta.size());
  
     ATH_MSG_DEBUG("atanh_a_css_eta:            " << m_atanh_a_css_eta);
     ATH_MSG_DEBUG("atanh_b_css_eta:            " << m_atanh_b_css_eta);
     ATH_MSG_DEBUG("atanh_c_css_eta:            " << m_atanh_c_css_eta);
     ATH_MSG_DEBUG("atanh_x0_css_eta:           " << m_atanh_x0_css_eta);
     ATH_MSG_DEBUG("atanh_a_csl_eta:            " << m_atanh_a_csl_eta);
     ATH_MSG_DEBUG("atanh_b_csl_eta:            " << m_atanh_b_csl_eta);
     ATH_MSG_DEBUG("atanh_c_csl_eta:            " << m_atanh_c_csl_eta);
     ATH_MSG_DEBUG("atanh_x0_csl_eta:           " << m_atanh_x0_csl_eta);
  
     ATH_MSG_DEBUG("Minimum pos error:          " << m_dposmin);
  
     return StatusCode::SUCCESS;
 }

◆ inputHandles()

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

◆ interfaceID()

static const InterfaceID& ICscClusterFitter::interfaceID ( )

inlinestaticinherited

Must declare this, with name of interface

Definition at line 89 of file ICscClusterFitter.h.

                                             { 
         static const InterfaceID IID_ICscClusterFitter("ICscClusterFitter", 1, 0);
         return IID_ICscClusterFitter;
     }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::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< AlgTool > >::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_alignmentTool

ToolHandle<ICscAlignmentTool> QratCscClusterFitter::m_alignmentTool

private

Initial value:

{
        this,
        "CscAlignmentTool",
        "CscAlignmentTool/CscAlignmentTool",
    }

Definition at line 88 of file QratCscClusterFitter.h.

◆ m_atanh_a_csl_eta

double QratCscClusterFitter::m_atanh_a_csl_eta

private

Definition at line 67 of file QratCscClusterFitter.h.

◆ m_atanh_a_css_eta

double QratCscClusterFitter::m_atanh_a_css_eta

private

Definition at line 63 of file QratCscClusterFitter.h.

◆ m_atanh_b_csl_eta

double QratCscClusterFitter::m_atanh_b_csl_eta

private

Definition at line 68 of file QratCscClusterFitter.h.

◆ m_atanh_b_css_eta

double QratCscClusterFitter::m_atanh_b_css_eta

private

Definition at line 64 of file QratCscClusterFitter.h.

◆ m_atanh_c_csl_eta

double QratCscClusterFitter::m_atanh_c_csl_eta

private

Definition at line 69 of file QratCscClusterFitter.h.

◆ m_atanh_c_css_eta

double QratCscClusterFitter::m_atanh_c_css_eta

private

Definition at line 65 of file QratCscClusterFitter.h.

◆ m_atanh_x0_csl_eta

double QratCscClusterFitter::m_atanh_x0_csl_eta

private

Definition at line 70 of file QratCscClusterFitter.h.

◆ m_atanh_x0_css_eta

double QratCscClusterFitter::m_atanh_x0_css_eta

private

Definition at line 66 of file QratCscClusterFitter.h.

◆ m_DetectorManagerKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> QratCscClusterFitter::m_DetectorManagerKey

private

Initial value:

{
        this,
        "DetectorManagerKey",
        "MuonDetectorManager",
        "Key of input MuonDetectorManager condition data",
    }

retrieve MuonDetectorManager from the conditions store

Definition at line 75 of file QratCscClusterFitter.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_dposmin

double QratCscClusterFitter::m_dposmin

private

Definition at line 72 of file QratCscClusterFitter.h.

◆ m_erropt_eta

std::string QratCscClusterFitter::m_erropt_eta

private

Definition at line 46 of file QratCscClusterFitter.h.

◆ m_erropt_phi

std::string QratCscClusterFitter::m_erropt_phi

private

Definition at line 47 of file QratCscClusterFitter.h.

◆ m_error_eta

double QratCscClusterFitter::m_error_eta

private

Definition at line 48 of file QratCscClusterFitter.h.

◆ m_error_phi

double QratCscClusterFitter::m_error_phi

private

Definition at line 49 of file QratCscClusterFitter.h.

◆ m_error_tantheta

double QratCscClusterFitter::m_error_tantheta

private

Definition at line 53 of file QratCscClusterFitter.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> QratCscClusterFitter::m_idHelperSvc

private

Initial value:

{
        this,
        "MuonIdHelperSvc",
        "Muon::MuonIdHelperSvc/MuonIdHelperSvc",
    }

Definition at line 82 of file QratCscClusterFitter.h.

◆ m_max_width

std::vector<unsigned int> QratCscClusterFitter::m_max_width

private

Definition at line 43 of file QratCscClusterFitter.h.

◆ m_posopt_eta

std::string QratCscClusterFitter::m_posopt_eta

private

Definition at line 44 of file QratCscClusterFitter.h.

◆ m_posopt_phi

std::string QratCscClusterFitter::m_posopt_phi

private

Definition at line 45 of file QratCscClusterFitter.h.

◆ m_precisionErrorScaler

double QratCscClusterFitter::m_precisionErrorScaler

private

Definition at line 50 of file QratCscClusterFitter.h.

◆ m_qrat_maxdiff

double QratCscClusterFitter::m_qrat_maxdiff

private

Definition at line 51 of file QratCscClusterFitter.h.

◆ m_qrat_maxsig

double QratCscClusterFitter::m_qrat_maxsig

private

Definition at line 52 of file QratCscClusterFitter.h.

◆ m_qratcor_csl_eta

std::vector<double> QratCscClusterFitter::m_qratcor_csl_eta

private

Definition at line 61 of file QratCscClusterFitter.h.

◆ m_qratcor_css_eta

std::vector<double> QratCscClusterFitter::m_qratcor_css_eta

private

Definition at line 60 of file QratCscClusterFitter.h.

◆ m_qratmin_csl_eta

double QratCscClusterFitter::m_qratmin_csl_eta

private

Definition at line 59 of file QratCscClusterFitter.h.

◆ m_qratmin_css_eta

double QratCscClusterFitter::m_qratmin_css_eta

private

Definition at line 58 of file QratCscClusterFitter.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_xtan_csl_eta_offset

double QratCscClusterFitter::m_xtan_csl_eta_offset

private

Definition at line 56 of file QratCscClusterFitter.h.

◆ m_xtan_csl_eta_slope

double QratCscClusterFitter::m_xtan_csl_eta_slope

private

Definition at line 57 of file QratCscClusterFitter.h.

◆ m_xtan_css_eta_offset

double QratCscClusterFitter::m_xtan_css_eta_offset

private

Definition at line 54 of file QratCscClusterFitter.h.

◆ m_xtan_css_eta_slope

double QratCscClusterFitter::m_xtan_css_eta_slope

private

Definition at line 55 of file QratCscClusterFitter.h.

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

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ ChargeList

◆ DataMap

◆ DataNames

◆ Results

◆ StoreGateSvc_t

◆ StripFit

◆ StripFitList

◆ StripList

Constructor & Destructor Documentation

◆ QratCscClusterFitter()

◆ ~QratCscClusterFitter()

Member Function Documentation

◆ dataNames()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ fit() [1/4]

◆ fit() [2/4]

◆ fit() [3/4]

◆ fit() [4/4]

◆ getCorrectedError()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_alignmentTool

◆ m_atanh_a_csl_eta

◆ m_atanh_a_css_eta

◆ m_atanh_b_csl_eta

◆ m_atanh_b_css_eta

◆ m_atanh_c_csl_eta

◆ m_atanh_c_css_eta

◆ m_atanh_x0_csl_eta

◆ m_atanh_x0_css_eta

◆ m_DetectorManagerKey

◆ m_detStore

◆ m_dposmin

◆ m_erropt_eta

◆ m_erropt_phi

◆ m_error_eta

◆ m_error_phi

◆ m_error_tantheta

◆ m_evtStore

◆ m_idHelperSvc

◆ m_max_width

◆ m_posopt_eta

◆ m_posopt_phi

◆ m_precisionErrorScaler

◆ m_qrat_maxdiff

◆ m_qrat_maxsig

◆ m_qratcor_csl_eta