#include <TGCHitsTestTool.h>

Inheritance diagram for TGCHitsTestTool:

Collaboration diagram for TGCHitsTestTool:

Public Member Functions
	TGCHitsTestTool (const std::string &type, const std::string &name, const IInterface *parent)
StatusCode	initialize ()
StatusCode	processEvent ()

Protected Member Functions
Identifier	getIdentifier (HitID)
StatusCode	checkIdentifier (Identifier)
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"}

Private Attributes
bool	m_DoTGCTest
const TgcIdHelper *	m_pTgcIdHelper

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

Detailed Description

Definition at line 15 of file TGCHitsTestTool.h.

Constructor & Destructor Documentation

◆ TGCHitsTestTool()

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

inline

Definition at line 20 of file TGCHitsTestTool.h.

    : MuonHitTestToolBase(type, name, parent),
    m_pTgcIdHelper(0)
  {
    declareProperty("DoTGCTest", m_DoTGCTest=true);
  }

Member Function Documentation

◆ checkIdentifier()

StatusCode TGCHitsTestTool::checkIdentifier ( Identifier offid )

protected

Definition at line 45 of file TGCHitsTestTool.cxx.

                                                            {
  const MuonGM::TgcReadoutElement* descriptor = m_pMuonMgr->getTgcReadoutElement(offid);
  if (!descriptor) {
    ATH_MSG_FATAL("TGC readout element not found for Id = " << m_pTgcIdHelper->show_to_string(offid));
    return StatusCode::FAILURE;
  }
 
  // -------- Testing Sim to Offline ID conversion ---------
  Amg::Vector3D tgc_gPos = descriptor->channelPos(offid);
  double tgc_strip_phi = tgc_gPos.phi();
  if (tgc_strip_phi < 0.) tgc_strip_phi += 2.*M_PI;
  /*
    double tgc_strip_perp  = tgc_gPos.perp();
    double tgc_strip_z   = tgc_gPos.z();
    double tgc_strip_cot   = 1./tan(tgc_gPos.theta());
  */
 
  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 TGCHitsTestTool::getIdentifier ( HitID tgchit )

protected

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.

Definition at line 23 of file TGCHitsTestTool.cxx.

                                                      {
  const TgcHitIdHelper* tgchelper = TgcHitIdHelper::GetHelper();
  std::string tgc_stname = tgchelper->GetStationName(tgchit);
  int         tgc_steta  = tgchelper->GetStationEta(tgchit);
  int         tgc_stphi  = tgchelper->GetStationPhi(tgchit);
  int         tgc_gg     = tgchelper->GetGasGap(tgchit);
  int         tgc_strip  = 1;//tgchelper->GetIsStrip(tgchit);
  int         tgc_channel = 1;//tgchelper->GetChannel(tgchit);
 
  Identifier offid = m_pTgcIdHelper->channelID(tgc_stname, tgc_steta, tgc_stphi, tgc_gg, tgc_strip, tgc_channel);
  if (offid == 0) {
    ATH_MSG_FATAL("TGC: Cannot build a valid Identifier; skip ");
  }
  ATH_MSG_VERBOSE(" TGC: Offline Id "<<m_pTgcIdHelper->show_to_string(offid));
 
  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 TGCHitsTestTool::initialize ( )

virtual

Reimplemented from MuonHitTestToolBase.

Definition at line 116 of file TGCHitsTestTool.cxx.

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

◆ processEvent()

StatusCode TGCHitsTestTool::processEvent ( )

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

Check the Hits identifiers, access the functions that give: Station name, station eta, station phi, doublet Z, doublet phi, doublet R, GasGap, Measures Phi. The values of these variables are written out to the AANtuple (variable content and range explained in the code section where AANTuple variables are registered)

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)

Definition at line 66 of file TGCHitsTestTool.cxx.

                                         {
  CHECK(executeCheckEventInfo());
 
  //MuonSpectrometer
 
  // TGC
  // Get the TGC histograms
  // Short description:
  //   MuonSpectrometer TGC cross section
  //   MuonSpectrometer TGC longitudinal
  //   Eta, Theta, Phi, z residual (...), phi residual (...) of hits in TGCs
 
  // Enter the main algorithm loop
  if (m_DoTGCTest) {
 
    const TGCSimHitCollection* p_collection = nullptr;
    if (evtStore()->retrieve(p_collection,"TGC_Hits") == StatusCode::SUCCESS) {
      for (const TGCSimHit& hit : *p_collection) {
        HitID tgchit= hit.TGCid();
        Identifier offid= getIdentifier(tgchit);
        CHECK(checkIdentifier(offid));
 
         //Check Hits
        GeoTGCHit 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_DoTGCTest

bool TGCHitsTestTool::m_DoTGCTest

private

Definition at line 38 of file TGCHitsTestTool.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_pTgcIdHelper

const TgcIdHelper* TGCHitsTestTool::m_pTgcIdHelper

private

Definition at line 40 of file TGCHitsTestTool.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

◆ TGCHitsTestTool()

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_DoTGCTest

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

◆ m_theta

◆ m_zResid