#include <MuonHitTestToolBase.h>

Inheritance diagram for MuonHitTestToolBase:

Collaboration diagram for MuonHitTestToolBase:

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

virtual StatusCode	initialize () override

Protected Member Functions
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"}

structors and AlgTool implementation
std::string	m_key
	The MC truth key. More...

HepMC::ConstGenParticlePtr	getPrimary ()

Detailed Description

Definition at line 21 of file MuonHitTestToolBase.h.

Constructor & Destructor Documentation

◆ MuonHitTestToolBase()

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

Definition at line 17 of file MuonHitTestToolBase.cxx.

   : SimTestToolBase(type, name, parent),
     m_detname("MDT"),
     m_pMuonMgr(0),
     m_muonevnt(0),m_muonrun(0),
     m_muoneta(0), m_muontheta(0), m_muonphi(0),
     m_muonzResid(0), m_muonphiResid(0),
     m_muondetBarrel(0), m_muonlongView(0),
     m_eta(0), m_theta(0), m_phi(0),
     m_zResid(0),m_phiResid(0),
     m_detBarrel(0), m_longView(0)
  
 {
   declareProperty("BarrelEtaCut", m_BarrelEtaCut=99.);
   declareProperty("DetectorName", m_detname="MDT");
  
 }

Member Function Documentation

◆ executeCheckEventInfo()

StatusCode MuonHitTestToolBase::executeCheckEventInfo ( )

protected

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 )

protected

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
	)

protected

◆ executeFillHistosSectors_Wedge1()

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

protected

◆ executeFillHistosSectors_Wedge2()

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

protected

◆ 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 MuonHitTestToolBase::initialize ( )

overridevirtual

Generic Muon part, these histograms are filled from all subalgorithms (MDT,RPC,CSC,TGC) For detailed description look in the .cxx code

specific part, these histograms are filled from subalgorithm For detailed description look, e.g, in the CSCHitsTestAlg.cxx code

Reimplemented from SimTestToolBase.

Reimplemented in MMHitsTestTool, sTGCHitsTestTool, MDTHitsTestTool, TGCHitsTestTool, RPCHitsTestTool, and CSCHitsTestTool.

Definition at line 111 of file MuonHitTestToolBase.cxx.

                                            {
  
   CHECK(detStore()->retrieve(m_pMuonMgr));
   ATH_CHECK (m_eventInfoKey.initialize());
  
   //MuonSpectrometer
   m_path+="Muon/";
   // will only create new if not already registered (uses m_path)
   _TH1D(m_muonevnt,"event_num",100,0.,1000.);
   _TH1D(m_muonrun,"run_num",100,-300.,300.);
  
   /*
     TH1D *etamuon=new TH1D("muonhitpos_eta","muonhitpos_eta",50,-5.,5.);
     registerHistogram("/truth/muonhitpos_eta",etamuon);
  
     TH1D *theta=new TH1D("muonhitpos_theta","muonhitpos_theta",50,-10.,10.);
     registerHistogram("/truth/muonhitpos_theta",theta);
  
     TH1D *phimuon=new TH1D("muonhitpos_phi","muonhitpos_phi",50,-5.,5.);
     registerHistogram("/truth/muonhitpos_phi",phimuon);
   */
  
   _TH1D(m_muonzResid,"muonhitpos_zResid",50,-300.,300.);
   _TH1D(m_muonphiResid,"muonhitpos_phiResid",50,-300.,300.);
  
   _TH2D(m_muondetBarrel,"muondet_barrel",200,-11000.,11000.,200,-11000.,11000.);
   _TH2D(m_muonlongView,"muonlong_view",200,-24000.,24000.,200,0.,13000.);
   //  ================================================================================
  
   m_path+=m_detname+"/";
   // book Histograms
   /*
     TH1D *eta=new TH1D("hitpos_eta","hitpos_eta",50,-5.,5.);
     registerHistogram("/truth/hitpos_eta",eta);
  
     TH1D *theta=new TH1D("hitpos_theta","hitpos_theta",50,-10.,10.);
     registerHistogram("/truth/hitpos_theta",theta);
  
     TH1D *phi=new TH1D("hitpos_phi","hitpos_phi",50,-5.,5.);
     registerHistogram("/truth/hitpos_phi",phi);
   */
   _TH1D(m_zResid,(m_detname+"_hitpos_zResid").c_str(),50,-300.,300.);
   _TH1D(m_phiResid,(m_detname+"_hitpos_phiResid").c_str(),50,-300.,300.);
  
   _TH2D(m_detBarrel,(m_detname+"_det_barrel").c_str(),200,-11000.,11000.,200,-11000.,11000.);
   _TH2D(m_longView,(m_detname+"_long_view").c_str(),200,-24000.,24000.,200,0.,14000.);
  
   return StatusCode::SUCCESS;
  
 }

Member Data Documentation

◆ m_BarrelEtaCut

double MuonHitTestToolBase::m_BarrelEtaCut

protected

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

protected

Definition at line 54 of file MuonHitTestToolBase.h.

◆ m_detname

std::string MuonHitTestToolBase::m_detname

protected

Definition at line 37 of file MuonHitTestToolBase.h.

◆ m_direction

Amg::Vector3D MuonHitTestToolBase::m_direction

protected

Definition at line 58 of file MuonHitTestToolBase.h.

◆ m_eta

TH1* MuonHitTestToolBase::m_eta

protected

Definition at line 52 of file MuonHitTestToolBase.h.

◆ m_eventInfoKey

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

protected

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

protected

Definition at line 54 of file MuonHitTestToolBase.h.

◆ m_muondetBarrel

TH2* MuonHitTestToolBase::m_muondetBarrel

protected

Definition at line 50 of file MuonHitTestToolBase.h.

◆ m_muoneta

TH1* MuonHitTestToolBase::m_muoneta

protected

Definition at line 48 of file MuonHitTestToolBase.h.

◆ m_muonevnt

TH1* MuonHitTestToolBase::m_muonevnt

protected

Definition at line 47 of file MuonHitTestToolBase.h.

◆ m_muonlongView

TH2 * MuonHitTestToolBase::m_muonlongView

protected

Definition at line 50 of file MuonHitTestToolBase.h.

◆ m_muonphi

TH1 * MuonHitTestToolBase::m_muonphi

protected

Definition at line 48 of file MuonHitTestToolBase.h.

◆ m_muonphiResid

TH1 * MuonHitTestToolBase::m_muonphiResid

protected

Definition at line 49 of file MuonHitTestToolBase.h.

◆ m_muonrun

TH1 * MuonHitTestToolBase::m_muonrun

protected

Definition at line 47 of file MuonHitTestToolBase.h.

◆ m_muontheta

TH1 * MuonHitTestToolBase::m_muontheta

protected

Definition at line 48 of file MuonHitTestToolBase.h.

◆ m_muonzResid

TH1* MuonHitTestToolBase::m_muonzResid

protected

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

protected

Definition at line 52 of file MuonHitTestToolBase.h.

◆ m_phiResid

TH1 * MuonHitTestToolBase::m_phiResid

protected

Definition at line 53 of file MuonHitTestToolBase.h.

◆ m_pMuonMgr

const MuonGM::MuonDetectorManager* MuonHitTestToolBase::m_pMuonMgr

protected

Definition at line 38 of file MuonHitTestToolBase.h.

◆ m_theta

TH1 * MuonHitTestToolBase::m_theta

protected

Definition at line 52 of file MuonHitTestToolBase.h.

◆ m_zResid

TH1* MuonHitTestToolBase::m_zResid

protected

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

structors and AlgTool implementation

Detailed Description

Constructor & Destructor Documentation

◆ MuonHitTestToolBase()

Member Function Documentation

◆ executeCheckEventInfo()

◆ executeFillHistos()

◆ executeFillHistos_sTGc()

◆ executeFillHistosSectors_Wedge1()

◆ executeFillHistosSectors_Wedge2()

◆ getPrimary()

◆ initialize()

Member Data Documentation

◆ m_BarrelEtaCut

◆ m_detBarrel

◆ m_detname

◆ m_direction

◆ 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_theta

◆ m_zResid