#include <sTGCHitsTestTool.h>

Inheritance diagram for sTGCHitsTestTool:

Collaboration diagram for sTGCHitsTestTool:

Public Member Functions
	sTGCHitsTestTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual StatusCode	initialize () override final
virtual StatusCode	processEvent () override final

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
TH2 *	m_sTgc_TransverseView_QS1C_posZ = nullptr
TH2 *	m_sTgc_TransverseView_QS2C_posZ = nullptr
TH2 *	m_sTgc_TransverseView_QS3C_posZ = nullptr
TH2 *	m_sTgc_TransverseView_QL1P_posZ = nullptr
TH2 *	m_sTgc_TransverseView_QL2P_posZ = nullptr
TH2 *	m_sTgc_TransverseView_QL3P_posZ = nullptr
TH2 *	m_sTgc_rZview_positiveZ = nullptr
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"}

Private Attributes
bool	m_DosTGCTest = true

structors and AlgTool implementation
std::string	m_key
	The MC truth key.
HepMC::ConstGenParticlePtr	getPrimary ()

Detailed Description

Definition at line 11 of file sTGCHitsTestTool.h.

Constructor & Destructor Documentation

◆ sTGCHitsTestTool()

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

inline

Definition at line 16 of file sTGCHitsTestTool.h.

    : MuonHitTestToolBase(type, name, parent)
  {
    declareProperty("DosTGCTest", m_DosTGCTest);
    m_BarrelEtaCut=10.;
  }

Member Function Documentation

◆ 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

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

finaloverridevirtual

Reimplemented from MuonHitTestToolBase.

Definition at line 106 of file sTGCHitsTestTool.cxx.

                                        {
  CHECK(MuonHitTestToolBase::initialize());
  _TH2D( m_sTgc_TransverseView_QS1C_posZ,"sTGC_TransverseView_QS1C_posZ",1200,-6000.,6000.,1200,-6000.,6000.);
  _TH2D( m_sTgc_TransverseView_QS2C_posZ,"sTGC_TransverseView_QS2C_posZ",1200,-6000.,6000.,1200,-6000.,6000.);
  _TH2D( m_sTgc_TransverseView_QS3C_posZ,"sTGC_TransverseView_QS3C_posZ",1200,-6000.,6000.,1200,-6000.,6000.);
 
  _TH2D( m_sTgc_TransverseView_QL1P_posZ,"sTGC_TransverseView_QL1P_posZ",1200,-6000.,6000.,1200,-6000.,6000.);
  _TH2D( m_sTgc_TransverseView_QL2P_posZ,"sTGC_TransverseView_QL2P_posZ",1200,-6000.,6000.,1200,-6000.,6000.);
  _TH2D( m_sTgc_TransverseView_QL3P_posZ,"sTGC_TransverseView_QL3P_posZ",1200,-6000.,6000.,1200,-6000.,6000.);
 
  _TH2D(m_sTgc_rZview_positiveZ,"sTGC_rZView_posZ",2000,6500.,8500.,5000,0.,5000.);
  return StatusCode::SUCCESS;
}

◆ processEvent()

StatusCode sTGCHitsTestTool::processEvent ( )

finaloverridevirtual

Definition at line 39 of file sTGCHitsTestTool.cxx.

                                          {
  CHECK(executeCheckEventInfo());
 
   if (m_DosTGCTest) {
    const sTGCSimHitCollection* p_collection = nullptr;
    CHECK(evtStore()->retrieve(p_collection,"sTGC_Hits"));
    for (const sTGCSimHit& hit : *p_collection) {
      Amg::Vector3D u = hit.globalPosition();
      CHECK(executeFillHistos(u));
      
      
      const sTgcHitIdHelper* hitHelper = sTgcHitIdHelper::GetHelper();
      int simId = hit.sTGCId();
      std::string sim_stationName = hitHelper->GetStationName(simId);
 
      static const std::string QS1C("QS1C");
      static const std::string QS2C("QS2C");
      static const std::string QS3C("QS3C");
      static const std::string QL1P("QL1P");
      static const std::string QL2P("QL2P");
      static const std::string QL3P("QL3P");
 
      if (sim_stationName==QS1C && u.z()>0){
        m_sTgc_TransverseView_QS1C_posZ->Fill(u.x(),u.y());
      }
 
 
      if (sim_stationName==QS2C && u.z()>0){
        m_sTgc_TransverseView_QS2C_posZ->Fill(u.x(),u.y());
      }
 
 
      if (sim_stationName==QS3C && u.z()>0){
        m_sTgc_TransverseView_QS3C_posZ->Fill(u.x(),u.y());
      }
 
      if (sim_stationName==QL1P && u.z()>0){
        m_sTgc_TransverseView_QL1P_posZ->Fill(u.x(),u.y());
      }
 
 
      if (sim_stationName==QL2P && u.z()>0){
        m_sTgc_TransverseView_QL2P_posZ->Fill(u.x(),u.y());
      }
 
      if (sim_stationName==QL3P && u.z()>0){
        m_sTgc_TransverseView_QL3P_posZ->Fill(u.x(),u.y());
      }
 
      double r_sTGc = sqrt(u.x()*u.x()+u.y()*u.y());
 
      if (u.z() > 0){
        m_sTgc_rZview_positiveZ->Fill(u.z(), r_sTGc);
      }
 
 
      // GeoMMHit 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_DosTGCTest

bool sTGCHitsTestTool::m_DosTGCTest = true

private

Definition at line 39 of file sTGCHitsTestTool.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.

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

◆ m_histSvc

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

protectedinherited

Definition at line 35 of file SimTestHisto.h.

35{"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

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.

34{"/truth/"};

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

TH2* sTGCHitsTestTool::m_sTgc_rZview_positiveZ = nullptr

protected

Definition at line 35 of file sTGCHitsTestTool.h.

◆ m_sTgc_TransverseView_QL1P_posZ

TH2* sTGCHitsTestTool::m_sTgc_TransverseView_QL1P_posZ = nullptr

protected

Definition at line 32 of file sTGCHitsTestTool.h.

◆ m_sTgc_TransverseView_QL2P_posZ

TH2* sTGCHitsTestTool::m_sTgc_TransverseView_QL2P_posZ = nullptr

protected

Definition at line 33 of file sTGCHitsTestTool.h.

◆ m_sTgc_TransverseView_QL3P_posZ

TH2* sTGCHitsTestTool::m_sTgc_TransverseView_QL3P_posZ = nullptr

protected

Definition at line 34 of file sTGCHitsTestTool.h.

◆ m_sTgc_TransverseView_QS1C_posZ

TH2* sTGCHitsTestTool::m_sTgc_TransverseView_QS1C_posZ = nullptr

protected

Definition at line 29 of file sTGCHitsTestTool.h.

◆ m_sTgc_TransverseView_QS2C_posZ

TH2* sTGCHitsTestTool::m_sTgc_TransverseView_QS2C_posZ = nullptr

protected

Definition at line 30 of file sTGCHitsTestTool.h.

◆ m_sTgc_TransverseView_QS3C_posZ

TH2* sTGCHitsTestTool::m_sTgc_TransverseView_QS3C_posZ = nullptr

protected

Definition at line 31 of file sTGCHitsTestTool.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

◆ sTGCHitsTestTool()

Member Function Documentation

◆ executeCheckEventInfo()

◆ executeFillHistos()

◆ executeFillHistos_sTGc()

◆ executeFillHistosSectors_Wedge1()

◆ executeFillHistosSectors_Wedge2()

◆ getPrimary()

◆ initialize()

◆ processEvent()

Member Data Documentation

◆ m_BarrelEtaCut

◆ m_detBarrel

◆ m_detname

◆ m_direction

◆ m_DosTGCTest

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

◆ m_sTgc_TransverseView_QL1P_posZ

◆ m_sTgc_TransverseView_QL2P_posZ

◆ m_sTgc_TransverseView_QL3P_posZ

◆ m_sTgc_TransverseView_QS1C_posZ

◆ m_sTgc_TransverseView_QS2C_posZ

◆ m_sTgc_TransverseView_QS3C_posZ

◆ m_theta

◆ m_zResid