MuonHitTestToolBase Class Reference

#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.
double	m_BarrelEtaCut
	MDT barrel eta cut, applicable to the MDT 2D cross section plot.
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.
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 CSCHitsTestTool, MDTHitsTestTool, MMHitsTestTool, RPCHitsTestTool, sTGCHitsTestTool, and TGCHitsTestTool.

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.

{this, "EventInfo", "EventInfo", "EventInfo name"};

m_histSvc

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

protectedinherited

Definition at line 35 of file SimTestHisto.h.

{"THistSvc", "SimTestHisto"};

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.

{"/truth/"};

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:

• MuonHitTestToolBase.h
• MuonHitTestToolBase.cxx