Muon::RPC_TimingTool Class Reference

#include <RPC_TimingTool.h>

Inheritance diagram for Muon::RPC_TimingTool:

Collaboration diagram for Muon::RPC_TimingTool:

Public Member Functions
	RPC_TimingTool (const std::string &, const std::string &, const IInterface *)
virtual	~RPC_TimingTool ()=default
virtual StatusCode	initialize () override final
virtual TimingResult	calculateTimingResult (const std::vector< const MuonClusterOnTrack * > &hits) const override final
	Calculate the time offset of a given set of hits wrt to the current bunch. More...

Private Member Functions
double	getError (const Muon::MuonClusterOnTrack &) const
	calculate error on the RPC time More...

Private Attributes
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

Detailed Description

Definition at line 24 of file RPC_TimingTool.h.

Constructor & Destructor Documentation

RPC_TimingTool()

RPC_TimingTool::RPC_TimingTool	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 11 of file RPC_TimingTool.cxx.

                                                                                            :
      base_class(t,n,p) {}

~RPC_TimingTool()

virtual Muon::RPC_TimingTool::~RPC_TimingTool ( )

virtualdefault

Member Function Documentation

calculateTimingResult()

IMuonHitTimingTool::TimingResult RPC_TimingTool::calculateTimingResult ( const std::vector< const MuonClusterOnTrack * > &  hits ) const

finaloverridevirtual

Calculate the time offset of a given set of hits wrt to the current bunch.

Definition at line 19 of file RPC_TimingTool.cxx.

                                                                                                                               {
    
    // sanity check
    if( hits.empty() )     return {};
    
    ATH_MSG_DEBUG(" calculating timing for hits " << hits.size() );
 
    std::vector<int> histogram(100,0);
    float minTime = -100.;
    float binwidth = 2*std::abs(minTime)/histogram.size();
    float invbinwidth = 1./binwidth;
    
    
    for( const auto *hit : hits ){
      const RpcClusterOnTrack* rpc = dynamic_cast<const RpcClusterOnTrack*>(hit);
      if( !rpc ) continue;
      int bin = invbinwidth*(rpc->time()-minTime);
      ATH_MSG_DEBUG(m_idHelperSvc->toString(hit->identify()) << " time " << rpc->time() << " bin " << bin );
      if( bin < 0 || bin >= (int)histogram.size() ) continue;
      ++histogram[bin];
    }
    int maxbin = -1;
    int maxval = 0;
    for( unsigned int i=0; i<histogram.size();++i ){
      int val = histogram[i];
      if( val > maxval ){
        maxbin = i;
        maxval = val;
      }
    }
    if( maxbin != -1 ){      
      double time = 0.;
      unsigned int nhits = 0;
      double error = getError(*hits.front());
      for( const auto *hit : hits ){
        const RpcClusterOnTrack* rpc = dynamic_cast<const RpcClusterOnTrack*>(hit);
        if( !rpc ) continue;
        int bin = invbinwidth*(rpc->time()-minTime);
        ATH_MSG_DEBUG(m_idHelperSvc->toString(hit->identify()) << " time " << rpc->time() << " bin " << bin );
        // select hits in the max bin or the two neighbouring ones
        if( bin > maxbin+1 || bin < maxbin-1 ) continue;
        time += rpc->time();
        ++nhits;
      }
      if( nhits == 0 ) return {};
      time /= nhits;
      ATH_MSG_DEBUG( " final time " << time << " error " << error );
      return {true,time,error};
    }
    return {};
  }

getError()

double RPC_TimingTool::getError ( const Muon::MuonClusterOnTrack &  rpcRIO ) const

private

calculate error on the RPC time

Definition at line 71 of file RPC_TimingTool.cxx.

                                                                        { 
 
    // S. Spagnolo: Aug. 5th 2011 - implementation for rel 17 reprocessing (based on slide 8 and 9 of 
    // https://indico.cern.ch/getFile.py/access?contribId=6&resId=0&materialId=slides&confId=147791
    // ... no calibration offline yet, no correction for signal propagation time)
    Identifier id = rpcRIO.identify();
    int measPhi = m_idHelperSvc->measuresPhi(id);
    MuonStationIndex::StIndex stIndex = m_idHelperSvc->stationIndex(id);
    if( stIndex == MuonStationIndex::BO) {
      // outer layer station => high pt 
      if (measPhi==1) return 5.10;
      return 4.84;     
    }else{
      // pivot plane of Middle station dbR=2 or low pT confirm dbR=1 but names like "BMS", "BML", etc
      if (measPhi==1) return 5.04;
      else return 4.18;      
    }
  }

initialize()

StatusCode RPC_TimingTool::initialize ( )

finaloverridevirtual

Definition at line 14 of file RPC_TimingTool.cxx.

                                        {
    ATH_CHECK(m_idHelperSvc.retrieve());
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> Muon::RPC_TimingTool::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 40 of file RPC_TimingTool.h.

---

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

• RPC_TimingTool.h
• RPC_TimingTool.cxx