TrtHitsTestTool Class Reference

#include <TrtHitsTestTool.h>

Inheritance diagram for TrtHitsTestTool:

Collaboration diagram for TrtHitsTestTool:

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

Protected Attributes
std::string	m_path {"/truth/"}
ServiceHandle< ITHistSvc >	m_histSvc {"THistSvc", "SimTestHisto"}

Private Attributes
std::string	m_collection
TH2 *	m_indetBarrel
TH2 *	m_indetLongView
TH1 *	m_hits_time_photons
TH1 *	m_hits_edep_photons
TH1 *	m_hits_time_nonphotons
TH1 *	m_hits_edep_nonphotons
TH2 *	m_hits_xy
TH2 *	m_hits_zr
TH2 *	m_hits_edep_zr_photons
TH2 *	m_hits_edep_zr_nonphotons
TH1 *	m_hits_log_barcode

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

Detailed Description

Definition at line 11 of file TrtHitsTestTool.h.

Constructor & Destructor Documentation

TrtHitsTestTool()

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

Definition at line 13 of file TrtHitsTestTool.cxx.

  : SimTestToolBase(type, name, parent),
    m_collection("TRTUncompressedHits"),
    m_indetBarrel(0), m_indetLongView(0),
    m_hits_time_photons(0), m_hits_edep_photons(0), 
    m_hits_time_nonphotons(0), m_hits_edep_nonphotons(0),
    m_hits_xy(0), m_hits_zr(0),
    m_hits_edep_zr_photons(0), m_hits_edep_zr_nonphotons(0),
    m_hits_log_barcode(0)
{  
    declareProperty("CollectionName",  m_collection="TRTUncompressedHits");
}

Member Function Documentation

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

virtual

Reimplemented from SimTestToolBase.

Definition at line 26 of file TrtHitsTestTool.cxx.

{
  // initialise general (if not done yet)
  m_path+="InDet/";
  _TH2D(m_indetBarrel,"idet_xy",100,-1200.,1200.,100,-1200.,1200.);
  _SET_TITLE(m_indetBarrel, "indet hit distribution [All]","x [mm]","y [mm]");
  _TH2D(m_indetLongView,"idet_zr",100,-3200.,3200.,100,0.,1200.);
  _SET_TITLE(m_indetLongView, "indet hit distribution [All]","z [mm]","r [mm]");
 
  if (m_collection=="TRTUncompressedHits") m_path+="TRT/";
  else if (m_collection=="PileupTRTUncompressedHits") m_path+="PileupTRT/";
  else { std::cerr<<"unsupported CollectionName for "<< this->name()<<std::endl; }
 
  _TH2D(m_hits_xy,"trt_xy",100,-1200.,1200.,100,-1200.,1200.);
  _SET_TITLE(m_hits_xy, "hit distribution [All]","x [mm]","y [mm]");
  _TH2D(m_hits_zr,"trt_zr",100,-3200.,3200.,100,0.,1200.);
  _SET_TITLE(m_hits_zr, "hit distribution [All]","z [mm]","r [mm]");
  _TH1D(m_hits_time_photons,"trt_edep_time_photons",100,0,25);
  _SET_TITLE(m_hits_time_photons, "energy weighted hit distribution [Photons]","t [ns]","dE/dt");
  _TH1D(m_hits_edep_photons,"trt_edep_photons",100,0.,40.);
  _SET_TITLE(m_hits_edep_photons, "TRT hit distribution [Photons]","edep [MeV]","dN/dE [1/MeV]");
  _TH2D(m_hits_edep_zr_photons,"trt_edep_zr_photons",100,-3200.,3200.,100,0.,1200.);
  _SET_TITLE(m_hits_edep_zr_photons, "TRT energy weighted indet hit distribution [Photons]","z [mm]","r [mm]");
  _TH1D(m_hits_time_nonphotons,"trt_edep_time_nonphotons",100,0,25);
  _SET_TITLE(m_hits_time_nonphotons, "energy weighted hit distribution [Non-Photons]","t [ns]","dE/dt");
  _TH1D(m_hits_edep_nonphotons,"trt_edep_nonphotons",100,0.,5.);
  _SET_TITLE(m_hits_edep_nonphotons, "TRT hit distribution [Non-Photons]","edep [MeV]","dN/dE [1/MeV]");
  _TH2D(m_hits_edep_zr_nonphotons,"trt_edep_zr_nonphotons",100,-3200.,3200.,100,0.,1200.);
  _SET_TITLE(m_hits_edep_zr_nonphotons, "TRT energy weighted indet hit distribution [Non-Photons]","z [mm]","r [mm]");
 
  _TH1D(m_hits_log_barcode,"trt_log_barcode",50,0.,16.);
  _SET_TITLE(m_hits_log_barcode, "log(barcode)","log(barcode)","dN/dlog(barcode)");
 
  return StatusCode::SUCCESS;
}

processEvent()

StatusCode TrtHitsTestTool::processEvent

(

)

Definition at line 63 of file TrtHitsTestTool.cxx.

                                         {
 
  const TRTUncompressedHitCollection* p_collection = nullptr;
  if (evtStore()->retrieve(p_collection,m_collection).isSuccess()) {
    for (const TRTUncompressedHit& hit : *p_collection) {
      GeoTRTUncompressedHit ghit(hit);
      HepGeom::Point3D<double> u = ghit.getGlobalPosition();
      m_indetBarrel->Fill(u.x(),u.y());
      m_indetLongView->Fill(u.z(),u.perp());
      m_hits_xy->Fill(u.x(),u.y());
      m_hits_zr->Fill(u.z(),u.perp());
      int barcode = HepMC::barcode(hit.particleLink()); // FIXME barcode-based
      m_hits_log_barcode->Fill( barcode > 0 ? log(barcode) : -1);
      int particleId(hit.GetParticleEncoding());
      if (particleId == 22 || static_cast<int>(abs(particleId)/100000) == 41 ||  static_cast<int>(abs(particleId)/10000000) == 1)
    {
      //criteria set to match TRT_HitCollectionCnv_p3.cxx
      //only photon energy is persistified properly.
      m_hits_time_photons->Fill(hit.GetGlobalTime(),hit.GetEnergyDeposit());
      m_hits_edep_photons->Fill(hit.GetEnergyDeposit());
      m_hits_edep_zr_photons->Fill(u.z(),u.perp(),hit.GetEnergyDeposit());
    }
      else
    {
      m_hits_time_nonphotons->Fill(hit.GetGlobalTime(),hit.GetEnergyDeposit());
      m_hits_edep_nonphotons->Fill(hit.GetEnergyDeposit());
      m_hits_edep_zr_nonphotons->Fill(u.z(),u.perp(),hit.GetEnergyDeposit());
    }
    }
  }
  else {
    ATH_MSG_WARNING("TRT Hit collection not found!");
  }
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_collection

std::string TrtHitsTestTool::m_collection

private

Definition at line 22 of file TrtHitsTestTool.h.

m_histSvc

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

protectedinherited

Definition at line 35 of file SimTestHisto.h.

m_hits_edep_nonphotons

TH1 * TrtHitsTestTool::m_hits_edep_nonphotons

private

Definition at line 27 of file TrtHitsTestTool.h.

m_hits_edep_photons

TH1 * TrtHitsTestTool::m_hits_edep_photons

private

Definition at line 26 of file TrtHitsTestTool.h.

m_hits_edep_zr_nonphotons

TH2 * TrtHitsTestTool::m_hits_edep_zr_nonphotons

private

Definition at line 29 of file TrtHitsTestTool.h.

m_hits_edep_zr_photons

TH2* TrtHitsTestTool::m_hits_edep_zr_photons

private

Definition at line 29 of file TrtHitsTestTool.h.

m_hits_log_barcode

TH1* TrtHitsTestTool::m_hits_log_barcode

private

Definition at line 30 of file TrtHitsTestTool.h.

m_hits_time_nonphotons

TH1* TrtHitsTestTool::m_hits_time_nonphotons

private

Definition at line 27 of file TrtHitsTestTool.h.

m_hits_time_photons

TH1* TrtHitsTestTool::m_hits_time_photons

private

Definition at line 26 of file TrtHitsTestTool.h.

m_hits_xy

TH2* TrtHitsTestTool::m_hits_xy

private

Definition at line 28 of file TrtHitsTestTool.h.

m_hits_zr

TH2 * TrtHitsTestTool::m_hits_zr

private

Definition at line 28 of file TrtHitsTestTool.h.

m_indetBarrel

TH2* TrtHitsTestTool::m_indetBarrel

private

Definition at line 24 of file TrtHitsTestTool.h.

m_indetLongView

TH2 * TrtHitsTestTool::m_indetLongView

private

Definition at line 24 of file TrtHitsTestTool.h.

m_key

std::string SimTestToolBase::m_key

protectedinherited

The MC truth key.

Definition at line 34 of file SimTestToolBase.h.

m_path

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

protectedinherited

Definition at line 34 of file SimTestHisto.h.

---

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

• TrtHitsTestTool.h
• TrtHitsTestTool.cxx