CscThresholdClusterBuilderTool Class Reference

#include <CscThresholdClusterBuilderTool.h>

Inheritance diagram for CscThresholdClusterBuilderTool:

Collaboration diagram for CscThresholdClusterBuilderTool:

Public Types
enum	NoiseOption { rms = 0, sigma, f001 }

Public Member Functions
	CscThresholdClusterBuilderTool (const std::string &type, const std::string &aname, const IInterface *)
	~CscThresholdClusterBuilderTool ()
StatusCode	initialize ()
	AlgTool InterfaceID. More...
StatusCode	getClusters (std::vector< IdentifierHash > &idVect, std::vector< IdentifierHash > &selectedIdVect, Muon::CscPrepDataContainer *object)
StatusCode	finalize ()
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
int	make_clusters (bool measphi, const std::vector< const Muon::CscStripPrepData * > &strips, Muon::CscPrepDataCollection *&collection)
StatusCode	getClusters (IdentifierHash idVect, std::vector< IdentifierHash > &selectedIdVect)
StatusCode	getClusters (IdentifierHash idVect, std::vector< IdentifierHash > &selectedIdVect, Muon::CscPrepDataContainer *pclusters)
StatusCode	getClusters (std::vector< IdentifierHash > &selectedIdVect, Muon::CscPrepDataContainer *pclusters)
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
double	m_threshold
float	m_kFactor
std::string	m_noiseOptionStr
NoiseOption	m_noiseOption
SG::ReadHandleKey< Muon::CscStripPrepDataContainer >	m_digit_key
ToolHandle< ICscCalibTool >	m_cscCalibTool
ToolHandle< ICscStripFitter >	m_pstrip_fitter
ToolHandle< ICscClusterFitter >	m_pfitter_def
ToolHandle< ICscClusterFitter >	m_pfitter_prec
ToolHandle< ICscClusterFitter >	m_pfitter_split
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager >	m_DetectorManagerKey
	retrieve MuonDetectorManager from the conditions store More...
bool	m_makeNarrowClusterThreeStrips
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 75 of file CscThresholdClusterBuilderTool.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

NoiseOption

enum CscThresholdClusterBuilderTool::NoiseOption

Enumerator
rms
sigma
f001

Definition at line 77 of file CscThresholdClusterBuilderTool.h.

{ rms = 0, sigma, f001 };

Constructor & Destructor Documentation

CscThresholdClusterBuilderTool()

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

Definition at line 69 of file CscThresholdClusterBuilderTool.cxx.

                                                                                         :
    AthAlgTool(type, aname, parent), m_noiseOption(rms), m_digit_key("CSC_Measurements") {
    declareInterface<ICscClusterBuilder>(this);
 
    declareProperty("threshold", m_threshold = 20000.0);
    declareProperty("kFactor", m_kFactor = 6.5);
    declareProperty("noiseOption", m_noiseOptionStr = "f001");
    declareProperty("digit_key", m_digit_key);
    declareProperty("makeNarrowClusterThreeStrips", m_makeNarrowClusterThreeStrips = true);
}

~CscThresholdClusterBuilderTool()

CscThresholdClusterBuilderTool::~CscThresholdClusterBuilderTool ( )

default

Member Function Documentation

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

finalize()

StatusCode CscThresholdClusterBuilderTool::finalize

(

)

Definition at line 347 of file CscThresholdClusterBuilderTool.cxx.

                                                    {
    ATH_MSG_VERBOSE("Finalizing " << name());
    return StatusCode::SUCCESS;
}

getClusters() [1/4]

StatusCode CscThresholdClusterBuilderTool::getClusters ( IdentifierHash  idVect, std::vector< IdentifierHash > &  selectedIdVect  )

private

getClusters() [2/4]

StatusCode CscThresholdClusterBuilderTool::getClusters ( IdentifierHash  idVect, std::vector< IdentifierHash > &  selectedIdVect, Muon::CscPrepDataContainer *  pclusters  )

private

Definition at line 164 of file CscThresholdClusterBuilderTool.cxx.

                                                                                            {
    // identifiers of collections already decoded and stored in the container will be skipped
    if (pclusters->indexFindPtr(givenHashId) != nullptr) {
        decodedIds.push_back(givenHashId);
        ATH_MSG_DEBUG("A collection already exists in the container for offline id hash. " << (int)givenHashId);
        return StatusCode::SUCCESS;
    }
 
    // Retrieve the CSC digits for this event.
    SG::ReadHandle<CscStripPrepDataContainer> pdigcon(m_digit_key);
    if (pdigcon.isValid()) {
        ATH_MSG_DEBUG("Retrieved strip container " << m_digit_key.key() << " with " << pdigcon->size() << " entries.");
    } else {
        ATH_MSG_WARNING("Failure to retrieve strip container " << m_digit_key.key());
        return StatusCode::SUCCESS;
    }
 
    //**********************************************
    // retrieve specific collection for the givenID
    const CscStripPrepDataCollection* col = pdigcon->indexFindPtr(givenHashId);
    if (nullptr == col) {
        unsigned int coll_hash = givenHashId;
        ATH_MSG_WARNING("Specific CSC Strip PrepData collection retrieving failed for collection hash = " << coll_hash);
        return StatusCode::SUCCESS;
    }
 
    ATH_MSG_DEBUG("Retrieved " << col->size() << " CSC Strip PrepDatas.");
 
    Identifier colid = col->identify();
    int istation = m_idHelperSvc->cscIdHelper().stationName(colid) - 49;
    int zsec = m_idHelperSvc->cscIdHelper().stationEta(colid);
    int phisec = m_idHelperSvc->cscIdHelper().stationPhi(colid);
 
    ATH_MSG_DEBUG("  Strip collection " << chamber(istation, zsec, phisec) << " has " << col->size() << " strips");
 
    // Create arrays to hold digits and cathode plane parameters.
    vector<const CscStripPrepData*> strips[8];
    int maxstrip[8] = {0, 0, 0, 0, 0, 0, 0, 0};
 
    // 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 StatusCode::FAILURE;
    }
 
    IdentifierHash hash;
    // Loop over digits and fill these arrays.
    for (CscStripPrepDataCollection::const_iterator idig = col->begin(); idig != col->end(); ++idig) {
        const CscStripPrepData& dig = **idig;
        Identifier did = dig.identify();
        hash = dig.collectionHash();
        const CscReadoutElement* pro = MuonDetMgr->getCscReadoutElement(did);
        int wlay = m_idHelperSvc->cscIdHelper().wireLayer(did);
        int measphi = m_idHelperSvc->cscIdHelper().measuresPhi(did);
        int idx = 2 * (wlay - 1) + measphi;
        // First entry for a cathode plane, initialize.
        if (maxstrip[idx] == 0) {
            maxstrip[idx] = pro->maxNumberOfStrips(measphi);
            for (int istrip = 0; istrip < maxstrip[idx]; ++istrip) strips[idx].push_back(nullptr);
        }
        int istrip = m_idHelperSvc->cscIdHelper().strip(did) - 1;
        if (istrip < 0 || istrip >= maxstrip[idx]) {
            ATH_MSG_WARNING("Invalid strip number");
            continue;
        }
        strips[idx][istrip] = &dig;
    }
 
    // Cluster.
    CscPrepDataCollection* newCollection = nullptr;
    for (int measphi = 0; measphi < 2; ++measphi) {
        for (int wlay = 1; wlay < 5; ++wlay) {
            int idx = 2 * (wlay - 1) + measphi;
            if (maxstrip[idx]) {
                make_clusters(measphi, strips[idx], newCollection);
                ATH_MSG_DEBUG("    " << wlay << "th layer ");
            }
        }
    }
    if (newCollection) {
        if (pclusters->addCollection(newCollection, hash).isFailure()) {
            ATH_MSG_ERROR("Couldn't add CscPrepdataCollection to container!");
            return StatusCode::RECOVERABLE;
        }
        decodedIds.push_back(hash);  // Record that this collection contains data
    }
 
    return StatusCode::SUCCESS;
}

getClusters() [3/4]

StatusCode CscThresholdClusterBuilderTool::getClusters ( std::vector< IdentifierHash > &  idVect, std::vector< IdentifierHash > &  selectedIdVect, Muon::CscPrepDataContainer *  object  )

virtual

Implements ICscClusterBuilder.

getClusters() [4/4]

StatusCode CscThresholdClusterBuilderTool::getClusters ( std::vector< IdentifierHash > &  selectedIdVect, Muon::CscPrepDataContainer *  pclusters  )

private

Definition at line 259 of file CscThresholdClusterBuilderTool.cxx.

                                                                                                                                 {
    // Retrieve the CSC digits for this event.
    SG::ReadHandle<CscStripPrepDataContainer> pdigcon(m_digit_key);
    if (pdigcon.isValid()) {
        ATH_MSG_DEBUG("Retrieved strip container " << m_digit_key.key() << " with " << pdigcon->size() << " entries.");
    } else {
        ATH_MSG_WARNING("Failure to retrieve strip container " << m_digit_key.key());
        return StatusCode::SUCCESS;
    }
 
    // Loop over digit collections.
    // This a loop over chambers (each with 4 wire planes).
    const CscStripPrepDataContainer& con = *pdigcon;
    for (CscStripPrepDataContainer::const_iterator icol = con.begin(); icol != con.end(); ++icol) {
        const CscStripPrepDataCollection& col = **icol;
        // check if the collection is already used
        if (pclusters->indexFindPtr(col.identifyHash()) != nullptr) {
            // store the identifier hash and continue
            decodedIds.push_back(col.identifyHash());
            continue;
        }
        Identifier colid = col.identify();
        int istation = m_idHelperSvc->cscIdHelper().stationName(colid) - 49;
        int zsec = m_idHelperSvc->cscIdHelper().stationEta(colid);
        int phisec = m_idHelperSvc->cscIdHelper().stationPhi(colid);
        ATH_MSG_DEBUG("**Strip collection " << chamber(istation, zsec, phisec) << " sector " << m_idHelperSvc->cscIdHelper().sector(colid)
                                            << " has " << col.size() << " strips");
 
        // Create arrays to hold digits and cathode plane parameters.
        vector<const CscStripPrepData*> strips[8];
        int maxstrip[8] = {0, 0, 0, 0, 0, 0, 0, 0};
 
        // 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 StatusCode::FAILURE;
        }
 
        IdentifierHash hash;
        // Loop over digits and fill these arrays.
        for (CscStripPrepDataCollection::const_iterator idig = col.begin(); idig != col.end(); ++idig) {
            const CscStripPrepData& dig = **idig;
            Identifier did = dig.identify();
            hash = dig.collectionHash();
            const CscReadoutElement* pro = MuonDetMgr->getCscReadoutElement(did);
            int wlay = m_idHelperSvc->cscIdHelper().wireLayer(did);
            int measphi = m_idHelperSvc->cscIdHelper().measuresPhi(did);
            int idx = 2 * (wlay - 1) + measphi;
            // First entry for a cathode plane, initialize.
            if (maxstrip[idx] == 0) {
                maxstrip[idx] = pro->maxNumberOfStrips(measphi);
                for (int istrip = 0; istrip < maxstrip[idx]; ++istrip) strips[idx].push_back(nullptr);
            }
            int istrip = m_idHelperSvc->cscIdHelper().strip(did) - 1;
            if (istrip < 0 || istrip >= maxstrip[idx]) {
                ATH_MSG_WARNING("Invalid strip number");
                continue;
            }
            strips[idx][istrip] = &dig;
        }
 
        // Cluster.
        CscPrepDataCollection* newCollection = nullptr;
        for (int measphi = 0; measphi < 2; ++measphi) {
            for (int wlay = 1; wlay < 5; ++wlay) {
                int idx = 2 * (wlay - 1) + measphi;
                if (maxstrip[idx]) {
                    ATH_MSG_DEBUG("*** " << chamber(istation, zsec, phisec) << " sector " << m_idHelperSvc->cscIdHelper().sector(colid)
                                         << "  " << wlay << "th layer ");
                    make_clusters(measphi, strips[idx], newCollection);
                }
            }
        }
        if (newCollection) {
            if (pclusters->addCollection(newCollection, hash).isFailure()) {
                ATH_MSG_ERROR("Couldn't add CscPrepdataCollection to container!");
                return StatusCode::FAILURE;
            }
            decodedIds.push_back(hash);  // Record that this collection contains data
        }
    }  // end loop over chambers
 
    return StatusCode::SUCCESS;
}

initialize()

StatusCode CscThresholdClusterBuilderTool::initialize

(

)

AlgTool InterfaceID.

Definition at line 89 of file CscThresholdClusterBuilderTool.cxx.

                                                      {
    // Display algorithm properties.
    ATH_MSG_DEBUG("Properties for " << name() << ":");
    ATH_MSG_DEBUG("  Strip threshold is Max( " << m_threshold << ", " << m_kFactor << "*stripNoise ) where stripNoise is from "
                                               << m_noiseOptionStr);
    ATH_CHECK(m_digit_key.initialize());
    if (m_noiseOptionStr != "rms" && m_noiseOptionStr != "sigma" && m_noiseOptionStr != "f001") {
        ATH_MSG_DEBUG(" noiseOption is not among rms/sigma/f001. rms is used for default!!");
        m_noiseOptionStr = "rms";
    }
    if (m_noiseOptionStr == "rms")
        m_noiseOption = rms;
    else if (m_noiseOptionStr == "sigma")
        m_noiseOption = sigma;
    else if (m_noiseOptionStr == "f001")
        m_noiseOption = f001;
 
    ATH_MSG_DEBUG("  Strip fitter is " << m_pstrip_fitter.typeAndName());
    ATH_MSG_DEBUG("  Default cluster fitter is " << m_pfitter_def.typeAndName());
    ATH_MSG_DEBUG("  Precision cluster fitter is " << m_pfitter_prec.typeAndName());
    ATH_MSG_DEBUG("  Split cluster fitter is " << m_pfitter_split.typeAndName());
    ATH_MSG_DEBUG("  Input digit key is " << m_digit_key.key());
 
    // CSC calibratin tool for the Condtiions Data base access //
    ATH_CHECK_RECOVERABLE(m_cscCalibTool.retrieve());
 
    // Retrieve the strip fitting tool.
    ATH_CHECK_RECOVERABLE(m_pstrip_fitter.retrieve());
    ATH_MSG_DEBUG("Retrieved strip fitting tool " << m_pstrip_fitter);
 
    // Retrieve the default cluster fitting tool.
    ATH_CHECK_RECOVERABLE(m_pfitter_def.retrieve());
    ATH_MSG_DEBUG("Retrieved CSC default cluster fitting tool");
 
    // Retrieve the precision cluster fitting tool.
    ATH_CHECK_RECOVERABLE(m_pfitter_prec.retrieve());
    ATH_MSG_DEBUG("Retrieved CSC precision cluster fitting tool");
 
    // Retrieve the split cluster fitting tool.
    ATH_CHECK_RECOVERABLE(m_pfitter_split.retrieve());
    ATH_MSG_DEBUG("Retrieved CSC split cluster fitting tool");
 
    // retrieve MuonDetectorManager from the conditions store
    ATH_CHECK(m_DetectorManagerKey.initialize());
    ATH_CHECK(m_idHelperSvc.retrieve());
 
    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& ICscClusterBuilder::interfaceID ( )

inlinestaticinherited

Must declare this, with name of interface

Definition at line 30 of file ICscClusterBuilder.h.

                                            {
        static const InterfaceID IID_ICscClusterBuilder("ICscClusterBuilder", 1, 0);
        return IID_ICscClusterBuilder;
    }

make_clusters()

int CscThresholdClusterBuilderTool::make_clusters ( bool  measphi, const std::vector< const Muon::CscStripPrepData * > &  strips, Muon::CscPrepDataCollection *&  collection  )

private

Make every cluster have three strips ///////////////////////////////////////////

Definition at line 364 of file CscThresholdClusterBuilderTool.cxx.

                                                                                         {
    // Loop over channels.
    unsigned int maxstrip = strips.size();
 
    ATH_MSG_DEBUG("    Clustering for " << setaphi(measphi) << " plane with " << maxstrip << " strips");
 
    // Phase I:
    //
    // Loop over strips and fetch the charge and time for each.
    // Also set flag indicating if this strip has pointer and charge is above threshold(active)
    ICscClusterFitter::StripFitList allStripfits;
    std::vector<bool> astrip;  // check active strip
    std::vector<bool> bstrip;  // check bad strip
    IdentifierHash cscHashId;
 
    // Always [0, 191] or [0, 47]
    for (unsigned int istrip = 0; istrip < strips.size(); ++istrip) {
        const CscStripPrepData* pstrip = strips[istrip];
        ICscClusterFitter::StripFit res;
        bool active = false;
        bool isBadChannel = false;
        if (pstrip) {
            if (!newCollection) {
                Identifier elementId = m_idHelperSvc->cscIdHelper().elementID(pstrip->identify());
                cscHashId = pstrip->collectionHash();
                newCollection = new CscPrepDataCollection(cscHashId);
                newCollection->setIdentifier(elementId);
            }
            res = m_pstrip_fitter->fit(*pstrip);
 
            IdentifierHash stripHash;
            Identifier stripId = pstrip->identify();
            if (m_idHelperSvc->cscIdHelper().get_channel_hash(stripId, stripHash)) {
                ATH_MSG_WARNING("Unable to get CSC striphash id "
                                << " the identifier is ");
                stripId.show();
            }
 
            if (res.stripStatus == Muon::CscStrStatHot || res.stripStatus == Muon::CscStrStatDead) isBadChannel = true;
 
            float stripNoise = 0;
            if (m_noiseOption == rms) {
                stripNoise = m_cscCalibTool->stripRMS(stripHash);
            } else if (m_noiseOption == sigma) {
                stripNoise = m_cscCalibTool->stripNoise(stripHash);
            } else if (m_noiseOption == f001) {  // f001 is rawADC +1
                stripNoise = m_cscCalibTool->stripF001(stripHash) - m_cscCalibTool->stripPedestal(stripHash);
                stripNoise /= 3.251;
            }
 
            active = res.charge >= m_threshold && res.charge >= m_kFactor * stripNoise;
            if (isBadChannel) active = false;  // Let's remove Bad channel First...
 
            if (msgLvl(MSG::DEBUG)) {
                // Log message.
                ostringstream strlog;
                strlog << "      Strip " << setw(3) << istrip + 1 << ": charge= " << setw(7) << int(res.charge) << " dcharge= " << setw(7)
                       << int(res.dcharge);
                if (std::fabs(res.time) < 1e8)
                    strlog << " time=" << setw(3) << int(res.time + 0.5);
                else
                    strlog << " time=OVERFLOW";
                if (active)
                    strlog << " *";
                else if (isBadChannel)
                    strlog << " b";
                else
                    strlog << " .";
                if (res.status)
                    strlog << " x";
                else
                    strlog << " o";
                strlog << "   :" << res.status;
                ATH_MSG_DEBUG(strlog.str());
            }
        }
        allStripfits.push_back(res);
        astrip.push_back(active);
        bstrip.push_back(isBadChannel);
    }
 
    // Phase II:
    //
    // Bad Channel recovery in case of strip above strip being nearby
 
    // 1. identify strips to recover
    std::vector<bool> rstrip;  // check recover strip
    bool IsAnyStripRecovered = false;
    for (unsigned int istrip = 0; istrip < strips.size(); ++istrip) {
        bool adjacentActive = false;
        if (bstrip[istrip]) {
            if (istrip > 0 && astrip[istrip - 1]) adjacentActive = true;
            if (istrip + 1 < strips.size() && astrip[istrip + 1]) adjacentActive = true;
            if (adjacentActive) IsAnyStripRecovered = true;
        }
        rstrip.push_back(adjacentActive);
    }
 
    // 2. make it active if strip to recover is not active
    if (IsAnyStripRecovered) {  // This loop is needed if there is any bad strip recovered because of adjacent active strip
 
        if (msgLvl(MSG::DEBUG)) {
            ostringstream checklog1;
            ostringstream checklog2;
 
            for (unsigned int istrip = 0; istrip < strips.size(); ++istrip) {
                if (istrip % 24 == 0) checklog1 << "\n";
                checklog1 << int(astrip[istrip]) << " ";
 
                if (!astrip[istrip] && rstrip[istrip]) {  // not active but bad strip with adjacent strip active
                    ATH_MSG_DEBUG("**** Strip " << istrip << " is recovered!!");
                }
                if (istrip % 24 == 0) checklog2 << "\n";
                checklog2 << int(astrip[istrip]) << " ";
            }
            ATH_MSG_DEBUG("Strip active map before and after");
            ATH_MSG_DEBUG(checklog1.str());
            ATH_MSG_DEBUG(checklog2.str());
        }
 
        for (unsigned int istrip = 0; istrip < strips.size(); ++istrip) {
            if (!astrip[istrip] && rstrip[istrip]) {  // not active but bad strip with adjacent strip active
                astrip[istrip] = rstrip[istrip];      // insert true
            }
        }
    }
    // Phase III:
    //
    // Find clusters : save first strip and nstrip
    vector<unsigned int> strip0s;
    vector<unsigned int> nstrips;
 
    // Loop over strips and create clusters.
    int nstrip = 0;
    int first_strip = 0;  // First strip in the cluster.
    bool incluster = false;
    for (unsigned int istrip = 0; istrip < strips.size(); ++istrip) {
        // If the current strip is not active, skip it.
        if (!astrip[istrip]) continue;
        assert(strips[istrip] != 0);  //  CscStripPrepData* pstrip = strips[istrip];
 
        if (!incluster) {
            incluster = true;
            nstrip = 0;
            first_strip = istrip;
        }
        ++nstrip;
 
        // If this is not the last strip in the plane, and the next
        // strip is active, add the next strip to the cluster.
        if (istrip != maxstrip - 1 && astrip[istrip + 1]) continue;
 
        // Recover narrow cluster
        if (!measphi && nstrip < 3) {
            bool collectLeftStrip = false;
            bool collectRightStrip = false;
 
            if (nstrip == 1) {
                if (int(istrip) >= nstrip     // left adjacent strip should be inside of CSC chamber >0
                    && istrip + 1 < maxstrip  // the other side strip should be available < 192
                    && (allStripfits[istrip - 1].charge > 0.1 && allStripfits[istrip + 1].charge > 0.1)  // charge should be positive
                    && strips[istrip - 1] && strips[istrip + 1]) {  // both adjacent strip identifier should exist
                    collectLeftStrip = true;
                    collectRightStrip = true;
                }
            } else if (nstrip == 2) {
                if (allStripfits[istrip - 1].charge > allStripfits[istrip].charge) {  // In case of left strip not fired
                    if (int(istrip) >= nstrip                                         // nstrip 2
                        && allStripfits[istrip - 2].charge > 0.1                      // charge should be positive
                        && strips[istrip - 2])                                        // left strip Identifier should exist
                        collectLeftStrip = true;
                } else {  // In case of right strip not fired
                    if (istrip + 1 < maxstrip &&
                        allStripfits[istrip + 1].charge > 0.1  // charge should be positive if 0, then 0.341E-134 will enter
                        && strips[istrip + 1])                 // right strip Identifier should exist
                        collectRightStrip = true;
                }
            }
 
            if (collectLeftStrip) {
                first_strip = first_strip - 1;
                nstrip += 1;
            }
            if (collectRightStrip) { nstrip += 1; }
 
            if (msgLvl(MSG::DEBUG)) {
                // Log message.
                ostringstream narrowlog;
                narrowlog << "  ** narrow Clusters   " << first_strip + 1 << "  " << nstrip << "  L:R  " << collectLeftStrip << " "
                          << collectRightStrip;
                for (int i = 0; i < nstrip; ++i) { narrowlog << "  " << allStripfits[first_strip + i].charge; }
                for (int i = 0; i < nstrip; ++i) { narrowlog << "  " << strips[first_strip + i]; }
                ATH_MSG_DEBUG(narrowlog.str());
            }
        }  // Only for eta plane nstrip <3
 
        strip0s.push_back(first_strip);
        nstrips.push_back(nstrip);
 
        // Reset incluster.
        incluster = false;
    }
 
    // Phase IV:
    //
    // Merge narrow cluster into adjacent cluster if any exists.
    // Only for eta strips...
    vector<unsigned int> newStrip0s;
    vector<unsigned int> newNstrips;
 
    int nMerged = 0;  // the difference b/w old Nclu and new Nclu
    for (unsigned int icl = 0; icl < nstrips.size(); ++icl) {
        unsigned int nstrip = nstrips[icl];
        unsigned int strip0 = strip0s[icl];
 
        ATH_MSG_VERBOSE("          " << icl << "th cluster merger " << strip0 << " " << nstrip);
 
        //#### if you find narrow cluster
        if (!measphi) {
            if (nstrip < 3) {
                // at least one cluster before to check left cluster and continuous
                if (icl > 0 && (strip0 == strip0s[icl - 1] + nstrips[icl - 1])) {
                    unsigned int newStrip0 = strip0s[icl - 1];
                    unsigned int newNstrip = nstrips[icl - 1] + nstrip;
 
                    ATH_MSG_DEBUG("   " << icl << "   ** narrow Cluster merger Type I" << newStrip0 << "  " << newNstrip);
 
                    newStrip0s[icl - 1 - nMerged] = newStrip0;
                    newNstrips[icl - 1 - nMerged] = newNstrip;
                    ++nMerged;
                    continue;
                }
                // at least one cluster after to check right cluster and continuous
                if (icl + 1 < nstrips.size() && (strip0 + nstrip == strip0s[icl + 1])) {
                    unsigned int newStrip0 = strip0;
                    unsigned int newNstrip = nstrip + nstrips[icl + 1];
 
                    ATH_MSG_DEBUG("   " << icl << "   ** narrow Cluster merger Type II" << newStrip0 << "  " << newNstrip);
 
                    newStrip0s.push_back(newStrip0);
                    newNstrips.push_back(newNstrip);
 
                    icl += 1;
                    ++nMerged;
                    continue;
                }
            }
        }  // !measphi
        // if nstrip >2 OR
        // still narrow strip then just keep it...
        newStrip0s.push_back(strip0);
        newNstrips.push_back(nstrip);
    }  // for
 
    if (strip0s.size() != newStrip0s.size()) {
        ATH_MSG_DEBUG(" Phase II -> III Merged " << strip0s.size() << ":" << nstrips.size() << "  " << newStrip0s.size() << ":"
                                                 << newNstrips.size());
        for (unsigned int icl = 0; icl < nstrips.size(); ++icl)
            ATH_MSG_DEBUG("   *** " << icl << "   [" << strip0s[icl] << "," << strip0s[icl] + nstrips[icl] - 1 << "]  " << nstrips[icl]);
        for (unsigned int icl = 0; icl < newNstrips.size(); ++icl)
            ATH_MSG_DEBUG("   ****** " << icl << "   [" << newStrip0s[icl] << "," << newStrip0s[icl] + newNstrips[icl] - 1 << "]  "
                                       << newNstrips[icl]);
    }
 
    // Phase V:
    //
    // Using strip0 and nstrip fill up collection
 
    ICscClusterFitter::StripFitList sfits;
    std::vector<const CscStripPrepData*> clusterStrips;
    clusterStrips.reserve(50);
    std::vector<Identifier> prd_digit_ids;
    prd_digit_ids.reserve(50);
 
    for (unsigned int icl = 0; icl < newNstrips.size(); ++icl) {  // for each cluster
 
        ATH_MSG_VERBOSE("      Creating " << icl << "th cluster");
 
        unsigned int nstrip = newNstrips[icl];  // only used here
        unsigned int strip0 = newStrip0s[icl];  // only used here
 
        sfits.clear();
        clusterStrips.clear();
        prd_digit_ids.clear();
 
        for (unsigned int ist = strip0; ist < strip0 + nstrip; ++ist) {
            const CscStripPrepData* pstrip = strips[ist];
            const ICscClusterFitter::StripFit& sfit = allStripfits[ist];
 
            sfits.push_back(sfit);
            clusterStrips.push_back(pstrip);
            prd_digit_ids.push_back(pstrip->identify());
        }
 
        ATH_MSG_VERBOSE("     ++++++++++++++  nstrip +++++   " << nstrip);
        if (nstrip < 3 && m_makeNarrowClusterThreeStrips) {
 
            bool leftToFill = false;
            bool rightToFill = false;
            if (nstrip == 1) {
                leftToFill = true;
                rightToFill = true;
            } else {
                if (sfits[0].charge > sfits[1].charge) {
                    leftToFill = true;
                } else if (sfits[0].charge < sfits[1].charge) {
                    rightToFill = true;
                } else {
                    ATH_MSG_DEBUG("  It should be CHECKED!!! ");
                    if (strip0 > 0) {
                        if (strips[strip0 - 1]) {
                            leftToFill = true;
                        } else if (strips[strip0 + 2]) {
                            rightToFill = true;
                        }
                    } else if (strips[strip0 + 2]) {
                        rightToFill = true;
                    }
                }
            }
 
            ATH_MSG_VERBOSE(" strip0  nstrip  filling left or right " << strip0 << "  " << nstrip << " " << leftToFill << " "
                                                                      << rightToFill);
            ATH_MSG_VERBOSE("  sfits[0] " << sfits[0].charge);
            if (nstrip == 2) ATH_MSG_VERBOSE("  sfits[1] " << sfits[1].charge);
 
            for (unsigned int i = 0; i < allStripfits.size(); ++i) { ATH_MSG_VERBOSE("index " << i << "  " << allStripfits[i].charge); }
 
            if (leftToFill) {
                //        ATH_MSG_DEBUG( " Left to fill " << allStripfits[strip0-1].charge);
                bool fillTheOtherSide = false;
                if (strip0 == 0) {
                    fillTheOtherSide = true;
                } else {
                    if (strips[strip0 - 1] == nullptr) fillTheOtherSide = true;
                }
 
                if (strip0 + nstrip >= allStripfits.size()) { fillTheOtherSide = false; }
 
                if (!fillTheOtherSide) {
                    if (strips[strip0 - 1]) {
                        sfits.insert(sfits.begin(), allStripfits[strip0 - 1]);
                        clusterStrips.insert(clusterStrips.begin(), strips[strip0 - 1]);
                        prd_digit_ids.insert(prd_digit_ids.begin(), strips[strip0 - 1]->identify());
                        strip0--;
                        nstrip = prd_digit_ids.size();
                    }
                } else {
                    if (strips[strip0 + nstrip]) {                       // for edge this can happen
                        sfits.push_back(allStripfits[strip0 + nstrip]);  // This is the case for example
                        // 12799.6 39183.9 39698
                        clusterStrips.push_back(strips[strip0 + nstrip]);
                        prd_digit_ids.push_back(strips[strip0 + nstrip]->identify());
                        nstrip = prd_digit_ids.size();
                    }
                }
            }
 
            if (rightToFill) {
                bool fillTheOtherSide = false;
                if (strip0 + nstrip >= allStripfits.size()) {
                    fillTheOtherSide = true;
                } else {
                    if (strips[strip0 + nstrip] == nullptr) fillTheOtherSide = true;
                }
 
                if (strip0 == 0) { fillTheOtherSide = false; }
 
                if (!fillTheOtherSide) {
                    if (strips[strip0 + nstrip]) {
                        sfits.push_back(allStripfits[strip0 + nstrip]);
                        clusterStrips.push_back(strips[strip0 + nstrip]);
                        prd_digit_ids.push_back(strips[strip0 + nstrip]->identify());
                        nstrip = prd_digit_ids.size();
                    }
                } else {
                    if (strips[strip0 - 1]) {  // for edge this can happen
                        sfits.insert(sfits.begin(), allStripfits[strip0 - 1]);
                        clusterStrips.insert(clusterStrips.begin(), strips[strip0 - 1]);
                        prd_digit_ids.insert(prd_digit_ids.begin(), strips[strip0 - 1]->identify());
                        strip0--;
                        nstrip = prd_digit_ids.size();
                    }
                }
            }
        }
 
        int fitresult = 99;
        std::vector<ICscClusterFitter::Result> results;
 
        // Precision fit.
        if (!measphi) {
            results = m_pfitter_prec->fit(sfits);
            fitresult = results[0].fitStatus;
            ATH_MSG_VERBOSE("        Performing precision fit " << m_pfitter_prec << " result return=" << fitresult);
 
            // in case of multipeak cluster
            if (fitresult == 6) {
                results = m_pfitter_split->fit(sfits);
                fitresult = results[0].fitStatus;
                for (unsigned int i = 0; i < results.size(); ++i)
                    ATH_MSG_VERBOSE("    Performing split fit with " << m_pfitter_split << " result return=" << results[i].fitStatus);
            }
        }
 
        bool precisionFitFailed = fitresult > 0 && fitresult < 20;  // splitclusterFit fail => 19
        // Default fit for phi and eta failed
        if (measphi || precisionFitFailed) {
            ICscClusterFitter::Result res;
            CscClusterStatus oldclustatus;
            if (!measphi) {
                res = results[0];
                oldclustatus = res.clusterStatus;
            } else {
                oldclustatus = Muon::CscStatusSimple;
            }
            results = m_pfitter_def->fit(sfits);
            if (!results.empty()) {
                res = results[0];
                fitresult = results[0].fitStatus;
                if (msgLvl(MSG::VERBOSE)) {
                    ostringstream deflog;
                    deflog << "      Performing default fit with " << m_pfitter_def;
                    if (fitresult) {
                        deflog << " failed: return=" << fitresult;
                    } else {
                        deflog << " succeeded";
                    }
                    ATH_MSG_VERBOSE(deflog.str());
                }
                // Keep the status from the first fit if it is defined.
                if (oldclustatus != Muon::CscStatusUndefined) {
                    res.clusterStatus = oldclustatus;
                    // we want to keep oldcluster status
                    results[0] = res;
                }
            }
        }
 
        //
        // Phase V. For multiple results, fill up collection
        //
        unsigned int nresults = results.size();
        for (unsigned int ire = 0; ire < nresults; ++ire) {
            CscClusterStatus clustatus = results[ire].clusterStatus;
            Muon::CscTimeStatus timeStatus = results[ire].timeStatus;
            double pos = results[ire].position;
            double err = results[ire].dposition;
            unsigned int id_strip = results[ire].strip;  // return peak strip index (unsigned integer)
            double cluster_charge = results[ire].charge;
            double cluster_time = results[ire].time;
            if (clustatus == Muon::CscStatusUndefined) ATH_MSG_DEBUG("      Csc Cluster Status is not defined.");
 
            if (id_strip >= sfits.size()) {
                ATH_MSG_WARNING("        Fit size check failed: ");
                continue;
            }
            // Fetch the strip used to identify this cluster.
            const CscStripPrepData* pstrip_id = nullptr;
            if (id_strip < clusterStrips.size()) pstrip_id = clusterStrips[id_strip];
            if (!pstrip_id) {
                ATH_MSG_WARNING("        Fit ID check failed: ");
                continue;
            }
 
            // Create ATLAS CSC cluster.
            Identifier cluster_id = pstrip_id->identify();
            IdentifierHash cluster_hash = pstrip_id->collectionHash();
            int zsec = m_idHelperSvc->cscIdHelper().stationEta(cluster_id);
            int wlay = m_idHelperSvc->cscIdHelper().wireLayer(cluster_id);
            // This local position is in the muon (not tracking) coordinate system.
            // 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 0;
            }
            const CscReadoutElement* pro = MuonDetMgr->getCscReadoutElement(cluster_id);
            // Amg::Vector3D local_pos = pro->localClusterPos(zsec, wlay, measphi, pos);
            Amg::Vector3D localTrk_pos = pro->nominalLocalClusterPos(zsec, wlay, measphi, pos);
 
            auto  cov = Amg::MatrixX(1, 1);
            (cov)(0, 0) = err * err;
            Amg::Vector2D plpos(measphi ? localTrk_pos.y() : localTrk_pos.z(), measphi ? localTrk_pos.z() : localTrk_pos.y());
            if (msgLvl(MSG::DEBUG)) {
                ATH_MSG_DEBUG("        Cluster parameters: " << nresults);
                ATH_MSG_DEBUG("                ID strip: " << first_strip + id_strip << "(" << first_strip << ":" << id_strip << ")");
                ATH_MSG_DEBUG("          local position: " << plpos.x() << " " << plpos.y() << " error: " << Amg::toString(cov));
                ATH_MSG_DEBUG("                  charge: " << cluster_charge);
                ATH_MSG_DEBUG("                    time: " << cluster_time);
                ATH_MSG_DEBUG("                  status: " << Muon::toString(clustatus));
            }
            unsigned int fstrip = results[ire].fstrip;
            unsigned int lstrip = results[ire].lstrip;
            std::vector<Identifier> prd_digit_ids_submit;
            for (unsigned int ids_index = fstrip; ids_index < lstrip + 1; ++ids_index) {
                if (ids_index >= prd_digit_ids.size()) {
                    ATH_MSG_WARNING("index out of range " << ids_index << " size " << prd_digit_ids.size());
                    continue;
                }
                prd_digit_ids_submit.push_back(prd_digit_ids[ids_index]);
            }
            unsigned int nstrip = prd_digit_ids_submit.size();
            ATH_MSG_DEBUG("                    size: " << nstrip << " " << sfits.size());
            ATH_MSG_DEBUG("   all              size: " << strips.size() << " " << allStripfits.size());
 
            //      allStripfits.push_back(res);
 
            CscPrepData* pclus = new CscPrepData(cluster_id,
                                                 cluster_hash,
                                                 plpos,
                                                 prd_digit_ids_submit,
                                                 cov,
                                                 pro,
                                                 int(cluster_charge + 0.5),
                                                 cluster_time,
                                                 clustatus,
                                                 timeStatus);
            pclus->setHashAndIndex(newCollection->identifyHash(),
                                   newCollection->size());
 
            newCollection->push_back(pclus);
        }
    }  // end loop over clusters
 
    return 0;
}

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_cscCalibTool

ToolHandle<ICscCalibTool> CscThresholdClusterBuilderTool::m_cscCalibTool

private

Initial value:

{
        this,
        "cscCalibTool",
        "CscCalibTool/CscCalibTool",
    }

Definition at line 118 of file CscThresholdClusterBuilderTool.h.

m_DetectorManagerKey

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

private

Initial value:

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

retrieve MuonDetectorManager from the conditions store

Definition at line 155 of file CscThresholdClusterBuilderTool.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_digit_key

SG::ReadHandleKey<Muon::CscStripPrepDataContainer> CscThresholdClusterBuilderTool::m_digit_key

private

Definition at line 115 of file CscThresholdClusterBuilderTool.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> CscThresholdClusterBuilderTool::m_idHelperSvc

private

Initial value:

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

Definition at line 148 of file CscThresholdClusterBuilderTool.h.

m_kFactor

float CscThresholdClusterBuilderTool::m_kFactor

private

Definition at line 112 of file CscThresholdClusterBuilderTool.h.

m_makeNarrowClusterThreeStrips

bool CscThresholdClusterBuilderTool::m_makeNarrowClusterThreeStrips

private

Definition at line 162 of file CscThresholdClusterBuilderTool.h.

m_noiseOption

NoiseOption CscThresholdClusterBuilderTool::m_noiseOption

private

Definition at line 114 of file CscThresholdClusterBuilderTool.h.

m_noiseOptionStr

std::string CscThresholdClusterBuilderTool::m_noiseOptionStr

private

Definition at line 113 of file CscThresholdClusterBuilderTool.h.

m_pfitter_def

ToolHandle<ICscClusterFitter> CscThresholdClusterBuilderTool::m_pfitter_def

private

Initial value:

{
        this,
        "default_fitter",
        "SimpleCscClusterFitter/SimpleCscClusterFitter",
    }

Definition at line 132 of file CscThresholdClusterBuilderTool.h.

m_pfitter_prec

ToolHandle<ICscClusterFitter> CscThresholdClusterBuilderTool::m_pfitter_prec

private

Initial value:

{
        this,
        "precision_fitter",
        "QratCscClusterFitter/QratCscClusterFitter",
    }

Definition at line 137 of file CscThresholdClusterBuilderTool.h.

m_pfitter_split

ToolHandle<ICscClusterFitter> CscThresholdClusterBuilderTool::m_pfitter_split

private

Initial value:

{
        this,
        "split_fitter",
        "CscSplitClusterFitter/CscSplitClusterFitter",
    }

Definition at line 142 of file CscThresholdClusterBuilderTool.h.

m_pstrip_fitter

ToolHandle<ICscStripFitter> CscThresholdClusterBuilderTool::m_pstrip_fitter

private

Initial value:

{
        this,
        "strip_fitter",
        "CalibCscStripFitter/CalibCscStripFitter",
    }

Definition at line 125 of file CscThresholdClusterBuilderTool.h.

m_threshold

double CscThresholdClusterBuilderTool::m_threshold

private

Definition at line 111 of file CscThresholdClusterBuilderTool.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.

---

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

• CscThresholdClusterBuilderTool.h
• CscThresholdClusterBuilderTool.cxx