ParabolaCscClusterFitter Class Reference

#include <ParabolaCscClusterFitter.h>

Inheritance diagram for ParabolaCscClusterFitter:

Collaboration diagram for ParabolaCscClusterFitter:

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
	ParabolaCscClusterFitter (const std::string &type, const std::string &aname, const IInterface *parent)
	Constructor. More...
	~ParabolaCscClusterFitter ()=default
StatusCode	initialize ()
const DataNames &	dataNames () const
	data names for ntuple output in csc_cluster tree More...
Results	fit (const StripFitList &sfits) const
Results	fit (const StripFitList &sfits, double dposdz) const
	Main fit method. More...
double	getCorrectedError (const Muon::CscPrepData *pclu, double slope) const
	Correct the positon error for track angle. More...
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 double	ParabolaCorrection (CscPlane &plane, double &raw)
	Correction of raw parabola positions. More...
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
	Max. More...
double	m_error_tantheta
	error contribution in mm for the tan(theta) track angle correction More...
double	m_xtan_css_eta_offset
	constant to Calculate the angle of incidence. More...
double	m_xtan_css_eta_slope
	constant to Calculate the angle of incidence. More...
double	m_xtan_csl_eta_offset
	constant to Calculate the angle of incidence. More...
double	m_xtan_csl_eta_slope
	constant to Calculate the angle of incidence. More...
double	m_multi
	threshold multiplier for cluster peak finding More...
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager >	m_DetectorManagerKey
	retrieve MuonDetectorManager from the conditions store More...
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 28 of file ParabolaCscClusterFitter.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

ParabolaCscClusterFitter()

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

Constructor.

This class is derived from AthAlgTool(type, aname, parent). The 3 args are used to initialize the base class. declareProperty(...) is used to declare job options.

Definition at line 102 of file ParabolaCscClusterFitter.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(3);                                                // CSL phi
    declareProperty("max_width", m_max_width);                               // Maximum width (strips) for unspoiled clusters
    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 = -.0045);  // in mm
    declareProperty("xtan_csl_eta_slope", m_xtan_csl_eta_slope = 0.000131);  // in mm
    declareProperty("multi", m_multi = 3.1);                                 //  threshold multiplier  for cluster peak finding
}

~ParabolaCscClusterFitter()

ParabolaCscClusterFitter::~ParabolaCscClusterFitter ( )

default

Member Function Documentation

dataNames()

const DataNames & ParabolaCscClusterFitter::dataNames ( ) const

virtual

data names for ntuple output in csc_cluster tree

Reimplemented from ICscClusterFitter.

Definition at line 137 of file ParabolaCscClusterFitter.cxx.

                                                           {
    static const DataNames dnames ={"qA","qB","qC"};
    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.

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

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

{ 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 ParabolaCscClusterFitter::fit ( const StripFitList &  sfits ) const

virtual

Implements ICscClusterFitter.

Definition at line 398 of file ParabolaCscClusterFitter.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 = std::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 ParabolaCscClusterFitter::fit ( const StripFitList &  sfits, double  tantheta  ) const

virtual

Main fit method.

The fit is really an interpolation, which uses a parabola through the 3 points given by (-1, qA), (0, qB), and (1, qC). This raw estimate is corrected to account for non-linearites and provides a position in the range of -0.5, 0.5. Finally, the pitch is used to convert the position in mm relative to the center of the chamber.

Parameters

sfits	The list of amplitudes in the StripFitList.
tantheta	The theta angle of the track, used for error estimate. @Note How do you get the track angle before this hit reconstruction?

Returns

The Results vector, usually containing only one result, that of the largest hit.

Todo:

Remove this useless test.

Todo:

check segment finding efficiency with Y pos = strip number.

Todo:

Base errors on error propagation in parabola.

Todo:

Check errors against residuals.

Todo:

Check for saturation and increase errors.

Todo:

Check for large amps and increase errors.

Todo:

Check for small amps and increase errors.

Todo:

Check cluster profile for overlapping hits.

Implements ICscClusterFitter.

Definition at line 155 of file ParabolaCscClusterFitter.cxx.

                                                                                      {
    ATH_MSG_VERBOSE("Parabola fit with tool " << name());
 
    Results results;  // Vector of fit results
 
    // Check that the input has at least three strips.
    unsigned int nstrip = sfits.size();  // number of strips in this cluster
    if (nstrip < 3) {
        ATH_MSG_VERBOSE(" CscStatusNarrow: Input has fewer than three strips.");
        results.emplace_back(1, Muon::CscStatusNarrow);
        return results;
    }
 
    // Check the input array for NULL pointers.
    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);
 
    //  const CscReadoutElement* pro = pstrip->detectorElement(); fixed by Woochun
    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;
    }
 
    // Count strips above threshold:
    ATH_MSG_VERBOSE("Parabola fitter input has " << nstrip << " strips:");
    //unsigned int nstrip_threshold = 0;
    for (unsigned int istrip = 0; istrip < nstrip; ++istrip) {
        Identifier id = sfits[istrip].strip->identify();
        //if (sfits[istrip].charge >= 20000) ++nstrip_threshold;
        ATH_MSG_VERBOSE(" index: " << istrip << " chn:" << m_idHelperSvc->cscIdHelper().strip(id)
                                   << " amp:" << (int)(sfits[istrip].charge / 1000) << " ke.");
    }
 
    // Find the highest peak and count all peaks above threshold.
    // Peaks have to be above threshold to avoid counting noise fluctuations
    unsigned int istrip_peak = 0;
    int numPeaks = 0;  // count peaks above threshold
    float qpeak = -1;  // charge of the peak strip
    // Loop over strips excluding the first and last strip
    double charge_clu = sfits[0].charge + sfits[nstrip - 1].charge;  // cluster sum
    for (unsigned int istrip = 1; istrip < nstrip - 1; ++istrip) {
        float qthis = sfits[istrip].charge;
        float qlast = sfits[istrip - 1].charge;
        float qnext = sfits[istrip + 1].charge;
        double thr = m_multi * sfits[istrip].dcharge / 10;  // correct noise*10
        Identifier id = sfits[istrip].strip->identify();
        ATH_MSG_VERBOSE(" index: " << istrip << " chn:" << m_idHelperSvc->cscIdHelper().strip(id)
                                   << " amp:" << (int)(sfits[istrip].charge / 1000) << " ke, thr: " << (int)(thr / 1000) << " ke "
                                   << ((qthis > thr) ? "signal" : "noise") << sfits[istrip].dcharge / 1000);
        charge_clu += qthis;
 
        // Peak if the adjacent strips have less charge.
        if ((qthis >= qlast) && (qthis >= qnext)) {  // There is a peak.
            if (qthis > qpeak) {                     // larger peak then before
                istrip_peak = istrip;                // record new peak location
                qpeak = qthis;
            }
            if (qthis > thr) numPeaks++;  // Only count peaks above a threshold
        }
    }  // next istrip
 
    ATH_MSG_VERBOSE(" Peak is at index " << istrip_peak << " amp = " << qpeak / 1000);
    if (numPeaks > 1) {  // Error if more than one peak above threshold was found.
        results.emplace_back(6, Muon::CscStatusMultiPeak);
        ATH_MSG_VERBOSE(" CscStatusMultiPeak: multiple peaks are found: " << numPeaks);
        return results;
    }
 
    if (istrip_peak == 0) {  // no peak (even below threshold) was found.
        results.emplace_back(11);
        ATH_MSG_VERBOSE(" No peak was found.");
        return results;
    }
 
    // Check that the peak is not on the chamber edge.
    // This cannot happen due to the prev. loop
    if (strip0 + istrip_peak <= 0 || strip0 + istrip_peak >= maxstrip - 1) {
        results.emplace_back(4, Muon::CscStatusEdge);
        ATH_MSG_VERBOSE(" CscStatusEdge: strip0+istrip_peak = " << strip0 + istrip_peak);
        return results;
    }
 
    // Cluster width cut, should not be done...
    /*****************************************************
    if ( nstrip_threshold > m_max_width[plane] ) {
      results.push_back(Result(5, Muon::CscStatusWide));
      ATH_MSG_VERBOSE(" CscStatusWide: nstrip_threshold = "
                      << nstrip_threshold);
      return results;
    }
    ********************************************************/
 
    // Cluster is spoiled if peak is not at the center.
    // should not be done
    /*******************************************************
    bool is_even = 2*(nstrip/2) == nstrip;
    bool atcenter = istrip_peak == nstrip/2 ||
                    (is_even && istrip_peak+1 == nstrip/2);
    if ( ! atcenter ) {
      results.push_back(Result(7, Muon::CscStatusSkewed));
      ATH_MSG_VERBOSE(" CscStatusSkewed: peak is not at the center.");
      return results;
    }
    *********************************************************/
 
    // reconstructed position in strip number units:
    double savg;
 
    // Charge amplitude per strip:
    double qA = sfits[istrip_peak - 1].charge;
    double qB = sfits[istrip_peak].charge;
    double qC = sfits[istrip_peak + 1].charge;
 
    // charge errors
    double dqA = sfits[istrip_peak - 1].dcharge;
    double dqB = sfits[istrip_peak].dcharge;
    double dqC = sfits[istrip_peak + 1].dcharge;
 
    double raw;                             // uncorrected parabola peak position
    double denominator = 2 * qB - qA - qC;  // denominator for parabola
    if (denominator <= 0) {                 // not a peak, should not happen
        ATH_MSG_WARNING("  Bad parabola denominator: " << denominator);
        results.emplace_back(9);
        return results;
    } else {
        raw = 0.5 * (qC - qA) / denominator;  // peak of parabola through (qA, qB, qC)
    }
 
    double pos = ParabolaCorrection(plane, raw);
    double cog = 1.2 * (qC - qA) / (qA + qB + qC);  // for comparism
    // pos = cog;
 
    // error calculation: S/N =  charge sum / ave charge error.
    double dpos = (dqA + dqB + dqC) / 3 / (qA + qB + qC) * pitch * std::sqrt(2.0);  // pos error
    dpos = 0.08;                                                                    // 80 micron error
    if (measphi) dpos = 2.5;                                                        // worse phi resolution of ~2.5 mm
    // dpos = 200;//*= 2; /** todo Tune this to correct chi 2 of segments */
    // Calculate strip position.
    savg = pos + strip0 + istrip_peak;  // in strip numbers
    //  savg = 0.1 + strip0 + istrip_peak; // in strip numbers
    ATH_MSG_VERBOSE(" Parabola correction: plane = '" << splane(plane) << "' qA=" << qA << " qB=" << qB << " qC=" << qC << " raw=" << raw
                                                      << " pos=" << pos << ", cog=" << cog);
 
    // Position error contribution from incident angle:
    double dpostht = m_error_tantheta * std::abs(tantheta);
 
    DataMap dmap;
    dmap["qA"] = qA;
    dmap["qB"] = qB;
    dmap["qC"] = qC;
 
    // Set return values.
    Result res(0, Muon::CscStatusUnspoiled);  // status 0???
    res.position = pitch * (savg + 0.5 - 0.5 * maxstrip);
    res.dposition = std::sqrt(dpos * dpos + dpostht * dpostht);  // pos. error estimate
    res.strip = istrip_peak;                                     // strip number in this cluster
    res.fstrip = 0;
    res.lstrip = nstrip - 1;           // cluster width
    res.charge = charge_clu;           // total cluster charge over all strips
    res.time = sfits[res.strip].time;  // peaking time at center strip
    res.qleft = qA;                    // left charge
    res.qpeak = qB;                    // center charge
    res.qright = qC;                   // right charge
    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 ParabolaCscClusterFitter::getCorrectedError ( const Muon::CscPrepData *  pclu, double  slope  ) const

virtual

Correct the positon error for track angle.

Note

This should be done in the segment, not here.

Todo:

Read, understand, and verify correctness. @pclu The CSC Prep Data for the cluster @slope The track angle

Returns

The new error estimate.

Implements ICscClusterFitter.

Definition at line 368 of file ParabolaCscClusterFitter.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 = std::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);
 
    return newError;
}

initialize()

StatusCode ParabolaCscClusterFitter::initialize

(

)

Definition at line 119 of file ParabolaCscClusterFitter.cxx.

                                                {
    ATH_MSG_VERBOSE("Initalizing " << name());
 
    ATH_CHECK(m_idHelperSvc.retrieve());
 
    ATH_MSG_DEBUG("Properties for " << name() << ":");
    ATH_MSG_DEBUG("    tan(theta) 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);
    // retrieve MuonDetectorManager from the conditions store
    ATH_CHECK(m_DetectorManagerKey.initialize());
 
    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.

ParabolaCorrection()

double ParabolaCscClusterFitter::ParabolaCorrection ( CscPlane &  plane, double &  raw  )

static

Correction of raw parabola positions.

Correct the raw parabola interpolation, based on CSCPlane.

Parameters

plane	The csc plane enum for small or large chamber, X or Y plane
raw	The raw value of the parabola peak positon

Returns

The corrected value

The correction function takes the output of the parabola method, the peak position, and returns a position that is closer to the predicted position. The correction depends on the strip pitch and was found experimentally during various beam tests and cosmics tests.

Parameters

plane	The plane of the detector: Large or Small, X or Y.
raw	The uncorrected value

Returns

The corrected value, often in range from -0.5 to 0.5.

Definition at line 65 of file ParabolaCscClusterFitter.cxx.

                                                                                {
    double a, b, c;  // correction values
    switch (plane) {
        case CSS_ETA: {       // small Chamber, X strips
            a = 2.97043e-01;  // from segment186049.root
            b = 8.73895e+00;
            c = 2.08465e-01;
            break;
        }
        case CSS_PHI: {  // small chamber, Y strips
            a = 0.2662;  // 2rt test with small chambers, longy.root
            b = 62.19;
            c = 0.2191;
            break;
        }
        case CSL_ETA: {       // large chamber, X strips
            a = 2.99744e-01;  // from segment186049.root
            b = 8.30484e+00;
            c = 2.07538e-01;
            break;
        }
        case CSL_PHI: {  // large chamber, Y strips
            a = 0.2766;  // Run 1782 in bat 184
            b = 50.08;
            c = 0.0301;
            break;
        }
        default: {  // invalid plane, should not happen, check beforehand.
            return -1.0;
        }
    }  // end switch
 
    return a * std::atan(b * raw) + c * raw;
}

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_DetectorManagerKey

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

private

Initial value:

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

retrieve MuonDetectorManager from the conditions store

Definition at line 77 of file ParabolaCscClusterFitter.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_error_tantheta

double ParabolaCscClusterFitter::m_error_tantheta

private

error contribution in mm for the tan(theta) track angle correction

Definition at line 61 of file ParabolaCscClusterFitter.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> ParabolaCscClusterFitter::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 74 of file ParabolaCscClusterFitter.h.

m_max_width

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

private

Max.

width of a cluster to be unspoiled. Layer type in array[4]

Definition at line 58 of file ParabolaCscClusterFitter.h.

m_multi

double ParabolaCscClusterFitter::m_multi

private

threshold multiplier for cluster peak finding

Definition at line 72 of file ParabolaCscClusterFitter.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 ParabolaCscClusterFitter::m_xtan_csl_eta_offset

private

constant to Calculate the angle of incidence.

Definition at line 67 of file ParabolaCscClusterFitter.h.

m_xtan_csl_eta_slope

double ParabolaCscClusterFitter::m_xtan_csl_eta_slope

private

constant to Calculate the angle of incidence.

Definition at line 69 of file ParabolaCscClusterFitter.h.

m_xtan_css_eta_offset

double ParabolaCscClusterFitter::m_xtan_css_eta_offset

private

constant to Calculate the angle of incidence.

Definition at line 63 of file ParabolaCscClusterFitter.h.

m_xtan_css_eta_slope

double ParabolaCscClusterFitter::m_xtan_css_eta_slope

private

constant to Calculate the angle of incidence.

Definition at line 65 of file ParabolaCscClusterFitter.h.

---

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

• ParabolaCscClusterFitter.h
• ParabolaCscClusterFitter.cxx