#include <RPCHitsTestTool.h>

Inheritance diagram for RPCHitsTestTool:

Collaboration diagram for RPCHitsTestTool:

Public Member Functions
	RPCHitsTestTool (const std::string &type, const std::string &name, const IInterface *parent)

StatusCode	initialize ()

StatusCode	processEvent ()

Protected Member Functions
Identifier	getIdentifier (HitID rpchit)

StatusCode	checkIdentifier (Identifier offid)

StatusCode	executeCheckEventInfo ()

StatusCode	executeFillHistos (const Amg::Vector3D &)

StatusCode	executeFillHistosSectors_Wedge1 (const Amg::Vector3D &, std::string)

StatusCode	executeFillHistosSectors_Wedge2 (const Amg::Vector3D &, std::string)

StatusCode	executeFillHistos_sTGc (const Amg::Vector3D &, std::string)

Protected Attributes
std::string	m_detname

const MuonGM::MuonDetectorManager *	m_pMuonMgr

SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfoKey {this, "EventInfo", "EventInfo", "EventInfo name"}
	SG key for Event Info. More...

double	m_BarrelEtaCut
	MDT barrel eta cut, applicable to the MDT 2D cross section plot. More...

TH1 *	m_muonevnt

TH1 *	m_muonrun

TH1 *	m_muoneta

TH1 *	m_muontheta

TH1 *	m_muonphi

TH1 *	m_muonzResid

TH1 *	m_muonphiResid

TH2 *	m_muondetBarrel

TH2 *	m_muonlongView

TH1 *	m_eta

TH1 *	m_theta

TH1 *	m_phi

TH1 *	m_zResid

TH1 *	m_phiResid

TH2 *	m_detBarrel

TH2 *	m_longView

Amg::Vector3D	m_direction

std::string	m_path {"/truth/"}

ServiceHandle< ITHistSvc >	m_histSvc {"THistSvc", "SimTestHisto"}

Private Attributes
bool	m_DoRPCTest

const RpcIdHelper *	m_pRpcIdHelper

structors and AlgTool implementation
HepMC::ConstGenParticlePtr	getPrimary ()

std::string	m_key
	The MC truth key. More...

Detailed Description

Definition at line 14 of file RPCHitsTestTool.h.

Constructor & Destructor Documentation

◆ RPCHitsTestTool()

RPCHitsTestTool::RPCHitsTestTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

inline

Definition at line 19 of file RPCHitsTestTool.h.

     : MuonHitTestToolBase(type, name, parent),
       m_pRpcIdHelper(0)
   {
     declareProperty("DoRPCTest", m_DoRPCTest=true);
   }

Member Function Documentation

◆ checkIdentifier()

StatusCode RPCHitsTestTool::checkIdentifier ( Identifier offid )

protected

Definition at line 48 of file RPCHitsTestTool.cxx.

                                                             {
  
   const MuonGM::RpcReadoutElement* descriptor = m_pMuonMgr->getRpcReadoutElement(offid);
   if (!descriptor) {
     ATH_MSG_FATAL("RPC readout element not found for Id = " << m_pRpcIdHelper->show_to_string(offid));
     return StatusCode::FAILURE;
   }
  
   // Testing Sim to Offline ID conversion
   Amg::Vector3D rpc_gPos   = descriptor->stripPos(offid);
   double rpc_strip_phi   = rpc_gPos.phi();
   if (rpc_strip_phi < 0.) rpc_strip_phi += 2.*M_PI;
  
   return StatusCode::SUCCESS;
 }

◆ executeCheckEventInfo()

StatusCode MuonHitTestToolBase::executeCheckEventInfo ( )

protectedinherited

Definition at line 35 of file MuonHitTestToolBase.cxx.

 {
   SG::ReadHandle<xAOD::EventInfo> eventInfo (m_eventInfoKey,Gaudi::Hive::currentContext());
   ATH_CHECK(eventInfo.isValid());
   uint64_t evt = eventInfo->eventNumber();
   int numrun = eventInfo->runNumber();
   ATH_MSG_VERBOSE("Processing EventInfo event #"<< evt<< " run: " << numrun);
   m_muonevnt->Fill(evt);
   m_muonrun->Fill(numrun);
   const McEventCollection* mcEvent;
   CHECK(evtStore()->retrieve(mcEvent,m_key));
  
   // *AS* Why only if mcEvent ==1? when would there be more than one event?
   if (mcEvent->size()!=1) {
     m_direction=Amg::Vector3D(0.,0.,0.);
     return StatusCode::SUCCESS;
   }
  
   // *AS* Why this (double) loop, if only the last entry is preserved?
   //    changed it to take the gen particle
   DataVector<HepMC::GenEvent>::const_iterator e;
   for (e=mcEvent->begin();e!=mcEvent->end(); ++e) {
     for (auto p: (**e)) {
       if ( !HepMC::is_simulation_particle(p)) {
         Amg::Vector3D temp_momentum(p->momentum().px(),
                                     p->momentum().py(),
                                     p->momentum().pz());
         m_direction = temp_momentum.unit();
         break;
       }
     }
   }
   return StatusCode::SUCCESS;
 }

◆ executeFillHistos()

StatusCode MuonHitTestToolBase::executeFillHistos ( const Amg::Vector3D & u )

protectedinherited

for MDTs that have barrel + endcap section, take only the barrel when plotting the xy-2d plot of the detector

Definition at line 70 of file MuonHitTestToolBase.cxx.

                                                                         {
  
   if (m_direction.perp() > 0 && fabs(m_direction.eta())<m_BarrelEtaCut){
     //mdtdet->Fill(u.x(),u.y());
     m_muondetBarrel->Fill(u.x(),u.y());
     m_detBarrel->Fill(u.x(),u.y());
   }
  
   double rad=sqrt(u.x()*u.x()+u.y()*u.y());
   m_muonlongView->Fill(u.z(),rad);
   m_longView->Fill(u.z(),rad);
  
   // //m_direction vector is filled with truth above, so here it is wrong (no eta, theta, phi of the hit!!!)
   // *AS* why not use "u"?
   // theta->Fill(m_direction.theta());
   // theta->Fill(m_direction.theta());
   // eta->Fill(m_direction.eta());
   // eta->Fill(m_direction.eta());
   // phi->Fill(m_direction.phi());
   // phi->Fill(m_direction.phi());
  
   if (m_direction.perp() > 0) {
     m_muonzResid->Fill(u.cross(m_direction).dot(m_direction.cross(Amg::Vector3D(0,0,1)).unit()));
     m_muonphiResid->Fill(u.cross(m_direction).z());
  
     m_zResid->Fill(u.cross(m_direction).dot(m_direction.cross(Amg::Vector3D(0,0,1)).unit()));
     m_phiResid->Fill(u.cross(m_direction).z());
   }
   else {
     m_muonzResid->Fill(0);
     m_muonphiResid->Fill(0);
     m_zResid->Fill(0);
     m_phiResid->Fill(0);
   }
  
   return StatusCode::SUCCESS;
 }

◆ executeFillHistos_sTGc()

StatusCode MuonHitTestToolBase::executeFillHistos_sTGc	(	const Amg::Vector3D &	,
		std::string
	)

protectedinherited

◆ executeFillHistosSectors_Wedge1()

StatusCode MuonHitTestToolBase::executeFillHistosSectors_Wedge1	(	const Amg::Vector3D &	,
		std::string
	)

protectedinherited

◆ executeFillHistosSectors_Wedge2()

StatusCode MuonHitTestToolBase::executeFillHistosSectors_Wedge2	(	const Amg::Vector3D &	,
		std::string
	)

protectedinherited

◆ getIdentifier()

Identifier RPCHitsTestTool::getIdentifier ( HitID rpchit )

protected

Definition at line 25 of file RPCHitsTestTool.cxx.

                                                       {
   const RpcHitIdHelper* rpchelper = RpcHitIdHelper::GetHelper(m_pRpcIdHelper->gasGapMax());
   std::string rpc_stname = rpchelper->GetStationName(rpchit);
   int         rpc_steta  = rpchelper->GetZSector(rpchit);
   int         rpc_stphi  = rpchelper->GetPhiSector(rpchit);
   int         rpc_dbr    = rpchelper->GetDoubletR(rpchit);
   int         rpc_dbz    = rpchelper->GetDoubletZ(rpchit);
   int         rpc_dbp    = rpchelper->GetDoubletPhi(rpchit);
   int         rpc_gg     = rpchelper->GetGasGapLayer(rpchit);
   int         rpc_mfi    = rpchelper->GetMeasuresPhi(rpchit);
   int         rpc_strip  = 1;//rpchelper->GetStrip(rpchit);
  
   // Create the offline identifiers, fill them with hit info. Ready to be
   // accessed and retrieve info. Currently not used in this code except for
   // some checks.
   Identifier offid = m_pRpcIdHelper->channelID(rpc_stname, rpc_steta, rpc_stphi, rpc_dbr, rpc_dbz, rpc_dbp, rpc_gg, rpc_mfi, rpc_strip);
   if (offid == 0) {
     ATH_MSG_FATAL("RPC : Cannot build a valid Identifier; skip ");
   }
  
   return offid;
 }

◆ getPrimary()

HepMC::ConstGenParticlePtr SimTestToolBase::getPrimary ( )

protectedinherited

Definition at line 20 of file SimTestToolBase.cxx.

 {
   const McEventCollection* mcCollection;
   if (evtStore()->retrieve(mcCollection,m_key).isSuccess()) {
     DataVector<HepMC::GenEvent>::const_iterator e;
     for (e=mcCollection->begin();e!=mcCollection->end(); ++e) {
       for (auto p : (**e)) {
     if (!HepMC::is_simulation_particle(p)) {
       return p;
     }
       }
     }
   }
   return 0;
 }

◆ initialize()

StatusCode RPCHitsTestTool::initialize ( )

virtual

Reimplemented from MuonHitTestToolBase.

Definition at line 104 of file RPCHitsTestTool.cxx.

                                        {
   CHECK(MuonHitTestToolBase::initialize());
   m_pRpcIdHelper = m_pMuonMgr->rpcIdHelper();
   return StatusCode::SUCCESS;
 }

◆ processEvent()

StatusCode RPCHitsTestTool::processEvent ( )

Definition at line 65 of file RPCHitsTestTool.cxx.

                                          {
   CHECK(executeCheckEventInfo());
  
   // MuonSpectrometer
   // Get the generic MuonSpectrometer (MDT,RPC,CSC,TGC) histograms
   // Short description:
   //   MuonSpectrometer cross section
   //   MuonSpectrometer longitudinal
   //   Number of events processed from truth
   //   Event specific run number
   //   Eta, Theta, Phi, z residual (...), phi residual (...) of hits in MuonSpectrometer
  
   if (m_DoRPCTest) {
  
     const RPCSimHitCollection* p_collection = nullptr;
     if (evtStore()->retrieve(p_collection,"RPC_Hits") == StatusCode::SUCCESS) {
       for (const RPCSimHit& hit : *p_collection) {
         // Check the Hits identifiers, access the functions that give:
         // Station name, station eta, station phi, doublet Z, doublet phi, doublet R, GasGap, Measures Phi.
         HitID rpchit= (hit).RPCid();
         Identifier offid= getIdentifier(rpchit);
         CHECK(checkIdentifier(offid));
        
         // Check Hits
         // For every hit within the event, get the global position Amg::Vector3D u and then retrieve all releveant info
         // either from the Amg::Vector3D or from the MC vector (direction)
         GeoRPCHit ghit(hit);
  
         if (!ghit) continue;
         Amg::Vector3D u = ghit.getGlobalPosition();
         CHECK(executeFillHistos(u));
       }
     }
   }
  
   return StatusCode::SUCCESS;
 }

Member Data Documentation

◆ m_BarrelEtaCut

double MuonHitTestToolBase::m_BarrelEtaCut

protectedinherited

MDT barrel eta cut, applicable to the MDT 2D cross section plot.

Definition at line 44 of file MuonHitTestToolBase.h.

◆ m_detBarrel

TH2* MuonHitTestToolBase::m_detBarrel

protectedinherited

Definition at line 54 of file MuonHitTestToolBase.h.

◆ m_detname

std::string MuonHitTestToolBase::m_detname

protectedinherited

Definition at line 37 of file MuonHitTestToolBase.h.

◆ m_direction

Amg::Vector3D MuonHitTestToolBase::m_direction

protectedinherited

Definition at line 58 of file MuonHitTestToolBase.h.

◆ m_DoRPCTest

bool RPCHitsTestTool::m_DoRPCTest

private

Definition at line 36 of file RPCHitsTestTool.h.

◆ m_eta

TH1* MuonHitTestToolBase::m_eta

protectedinherited

Definition at line 52 of file MuonHitTestToolBase.h.

◆ m_eventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> MuonHitTestToolBase::m_eventInfoKey {this, "EventInfo", "EventInfo", "EventInfo name"}

protectedinherited

SG key for Event Info.

Definition at line 41 of file MuonHitTestToolBase.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> SimTestHisto::m_histSvc {"THistSvc", "SimTestHisto"}

protectedinherited

Definition at line 35 of file SimTestHisto.h.

◆ m_key

std::string SimTestToolBase::m_key

protectedinherited

The MC truth key.

Definition at line 34 of file SimTestToolBase.h.

◆ m_longView

TH2 * MuonHitTestToolBase::m_longView

protectedinherited

Definition at line 54 of file MuonHitTestToolBase.h.

◆ m_muondetBarrel

TH2* MuonHitTestToolBase::m_muondetBarrel

protectedinherited

Definition at line 50 of file MuonHitTestToolBase.h.

◆ m_muoneta

TH1* MuonHitTestToolBase::m_muoneta

protectedinherited

Definition at line 48 of file MuonHitTestToolBase.h.

◆ m_muonevnt

TH1* MuonHitTestToolBase::m_muonevnt

protectedinherited

Definition at line 47 of file MuonHitTestToolBase.h.

◆ m_muonlongView

TH2 * MuonHitTestToolBase::m_muonlongView

protectedinherited

Definition at line 50 of file MuonHitTestToolBase.h.

◆ m_muonphi

TH1 * MuonHitTestToolBase::m_muonphi

protectedinherited

Definition at line 48 of file MuonHitTestToolBase.h.

◆ m_muonphiResid

TH1 * MuonHitTestToolBase::m_muonphiResid

protectedinherited

Definition at line 49 of file MuonHitTestToolBase.h.

◆ m_muonrun

TH1 * MuonHitTestToolBase::m_muonrun

protectedinherited

Definition at line 47 of file MuonHitTestToolBase.h.

◆ m_muontheta

TH1 * MuonHitTestToolBase::m_muontheta

protectedinherited

Definition at line 48 of file MuonHitTestToolBase.h.

◆ m_muonzResid

TH1* MuonHitTestToolBase::m_muonzResid

protectedinherited

Definition at line 49 of file MuonHitTestToolBase.h.

◆ m_path

std::string SimTestHisto::m_path {"/truth/"}

protectedinherited

Definition at line 34 of file SimTestHisto.h.

◆ m_phi

TH1 * MuonHitTestToolBase::m_phi

protectedinherited

Definition at line 52 of file MuonHitTestToolBase.h.

◆ m_phiResid

TH1 * MuonHitTestToolBase::m_phiResid

protectedinherited

Definition at line 53 of file MuonHitTestToolBase.h.

◆ m_pMuonMgr

const MuonGM::MuonDetectorManager* MuonHitTestToolBase::m_pMuonMgr

protectedinherited

Definition at line 38 of file MuonHitTestToolBase.h.

◆ m_pRpcIdHelper

const RpcIdHelper* RPCHitsTestTool::m_pRpcIdHelper

private

Definition at line 38 of file RPCHitsTestTool.h.

◆ m_theta

TH1 * MuonHitTestToolBase::m_theta

protectedinherited

Definition at line 52 of file MuonHitTestToolBase.h.

◆ m_zResid

TH1* MuonHitTestToolBase::m_zResid

protectedinherited

Definition at line 53 of file MuonHitTestToolBase.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

structors and AlgTool implementation

Detailed Description

Constructor & Destructor Documentation

◆ RPCHitsTestTool()

Member Function Documentation

◆ checkIdentifier()

◆ executeCheckEventInfo()

◆ executeFillHistos()

◆ executeFillHistos_sTGc()

◆ executeFillHistosSectors_Wedge1()

◆ executeFillHistosSectors_Wedge2()

◆ getIdentifier()

◆ getPrimary()

◆ initialize()

◆ processEvent()

Member Data Documentation

◆ m_BarrelEtaCut

◆ m_detBarrel

◆ m_detname

◆ m_direction

◆ m_DoRPCTest

◆ m_eta

◆ m_eventInfoKey

◆ m_histSvc

◆ m_key

◆ m_longView

◆ m_muondetBarrel

◆ m_muoneta

◆ m_muonevnt

◆ m_muonlongView

◆ m_muonphi

◆ m_muonphiResid

◆ m_muonrun

◆ m_muontheta

◆ m_muonzResid

◆ m_path

◆ m_phi

◆ m_phiResid

◆ m_pMuonMgr

◆ m_pRpcIdHelper

◆ m_theta

◆ m_zResid