ATLAS Offline Software
Loading...
Searching...
No Matches
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
const InDetDD::TRT_DetectorManagerm_detMgr {nullptr}

structors and AlgTool implementation

std::string m_key
 The MC truth key.
HepMC::ConstGenParticlePtr getPrimary ()

Detailed Description

Definition at line 14 of file TrtHitsTestTool.h.

Constructor & Destructor Documentation

◆ TrtHitsTestTool()

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

Definition at line 14 of file TrtHitsTestTool.cxx.

15 : SimTestToolBase(type, name, parent),
16 m_collection("TRTUncompressedHits"),
20 m_hits_xy(0), m_hits_zr(0),
23{
24 declareProperty("CollectionName", m_collection="TRTUncompressedHits");
25}
SimTestToolBase(const std::string &type, const std::string &name, const IInterface *parent)
std::string m_collection
TH2 * m_hits_edep_zr_nonphotons

Member Function Documentation

◆ getPrimary()

HepMC::ConstGenParticlePtr SimTestToolBase::getPrimary ( )
protectedinherited

Definition at line 20 of file SimTestToolBase.cxx.

21{
22 const McEventCollection* mcCollection;
23 if (evtStore()->retrieve(mcCollection,m_key).isSuccess()) {
25 for (e=mcCollection->begin();e!=mcCollection->end(); ++e) {
26 for (auto p : (**e)) {
28 return p;
29 }
30 }
31 }
32 }
33 return 0;
34}
DataModel_detail::const_iterator< DataVector > const_iterator
Standard const_iterator.
Definition DataVector.h:838
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
std::string m_key
The MC truth key.
bool is_simulation_particle(const T &p)
Method to establish if a particle (or barcode) was created during the simulation (TODO update to be s...
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

◆ initialize()

StatusCode TrtHitsTestTool::initialize ( )
virtual

Reimplemented from SimTestToolBase.

Definition at line 27 of file TrtHitsTestTool.cxx.

28{
29 // initialise general (if not done yet)
30 m_path+="InDet/";
31 _TH2D(m_indetBarrel,"idet_xy",100,-1200.,1200.,100,-1200.,1200.);
32 _SET_TITLE(m_indetBarrel, "indet hit distribution [All]","x [mm]","y [mm]");
33 _TH2D(m_indetLongView,"idet_zr",100,-3200.,3200.,100,0.,1200.);
34 _SET_TITLE(m_indetLongView, "indet hit distribution [All]","z [mm]","r [mm]");
35
36 if (m_collection=="TRTUncompressedHits") m_path+="TRT/";
37 else if (m_collection=="PileupTRTUncompressedHits") m_path+="PileupTRT/";
38 else { std::cerr<<"unsupported CollectionName for "<< this->name()<<std::endl; }
39
40 _TH2D(m_hits_xy,"trt_xy",100,-1200.,1200.,100,-1200.,1200.);
41 _SET_TITLE(m_hits_xy, "hit distribution [All]","x [mm]","y [mm]");
42 _TH2D(m_hits_zr,"trt_zr",100,-3200.,3200.,100,0.,1200.);
43 _SET_TITLE(m_hits_zr, "hit distribution [All]","z [mm]","r [mm]");
44 _TH1D(m_hits_time_photons,"trt_edep_time_photons",100,0,25);
45 _SET_TITLE(m_hits_time_photons, "energy weighted hit distribution [Photons]","t [ns]","dE/dt");
46 _TH1D(m_hits_edep_photons,"trt_edep_photons",100,0.,40.);
47 _SET_TITLE(m_hits_edep_photons, "TRT hit distribution [Photons]","edep [MeV]","dN/dE [1/MeV]");
48 _TH2D(m_hits_edep_zr_photons,"trt_edep_zr_photons",100,-3200.,3200.,100,0.,1200.);
49 _SET_TITLE(m_hits_edep_zr_photons, "TRT energy weighted indet hit distribution [Photons]","z [mm]","r [mm]");
50 _TH1D(m_hits_time_nonphotons,"trt_edep_time_nonphotons",100,0,25);
51 _SET_TITLE(m_hits_time_nonphotons, "energy weighted hit distribution [Non-Photons]","t [ns]","dE/dt");
52 _TH1D(m_hits_edep_nonphotons,"trt_edep_nonphotons",100,0.,5.);
53 _SET_TITLE(m_hits_edep_nonphotons, "TRT hit distribution [Non-Photons]","edep [MeV]","dN/dE [1/MeV]");
54 _TH2D(m_hits_edep_zr_nonphotons,"trt_edep_zr_nonphotons",100,-3200.,3200.,100,0.,1200.);
55 _SET_TITLE(m_hits_edep_zr_nonphotons, "TRT energy weighted indet hit distribution [Non-Photons]","z [mm]","r [mm]");
56
57 _TH1D(m_hits_log_barcode,"trt_log_barcode",50,0.,16.);
58 _SET_TITLE(m_hits_log_barcode, "log(barcode)","log(barcode)","dN/dlog(barcode)");
59
61
62 return StatusCode::SUCCESS;
63}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define _TH2D(var, name, nbinx, xmin, xmax, nbiny, ymin, ymax)
#define _TH1D(var, name, nbin, xmin, xmax)
#define _SET_TITLE(var, title, xaxis, yaxis)
std::string m_path
const InDetDD::TRT_DetectorManager * m_detMgr

◆ processEvent()

StatusCode TrtHitsTestTool::processEvent ( )

Definition at line 66 of file TrtHitsTestTool.cxx.

66 {
67
68 const TRTUncompressedHitCollection* p_collection = nullptr;
69 if (evtStore()->retrieve(p_collection,m_collection).isSuccess()) {
70 for (const TRTUncompressedHit& hit : *p_collection) {
71 GeoTRTUncompressedHit ghit(hit);
72 HepGeom::Point3D<double> u = ghit.getGlobalPosition(m_detMgr);
73 m_indetBarrel->Fill(u.x(),u.y());
74 m_indetLongView->Fill(u.z(),u.perp());
75 m_hits_xy->Fill(u.x(),u.y());
76 m_hits_zr->Fill(u.z(),u.perp());
77 int barcode = HepMC::barcode(hit.particleLink()); // FIXME barcode-based
78 m_hits_log_barcode->Fill( barcode > 0 ? log(barcode) : -1);
79 int particleId(hit.GetParticleEncoding());
80 if (particleId == 22 || static_cast<int>(abs(particleId)/100000) == 41 || static_cast<int>(abs(particleId)/10000000) == 1)
81 {
82 //criteria set to match TRT_HitCollectionCnv_p3.cxx
83 //only photon energy is persistified properly.
84 m_hits_time_photons->Fill(hit.GetGlobalTime(),hit.GetEnergyDeposit());
85 m_hits_edep_photons->Fill(hit.GetEnergyDeposit());
86 m_hits_edep_zr_photons->Fill(u.z(),u.perp(),hit.GetEnergyDeposit());
87 }
88 else
89 {
90 m_hits_time_nonphotons->Fill(hit.GetGlobalTime(),hit.GetEnergyDeposit());
91 m_hits_edep_nonphotons->Fill(hit.GetEnergyDeposit());
92 m_hits_edep_zr_nonphotons->Fill(u.z(),u.perp(),hit.GetEnergyDeposit());
93 }
94 }
95 }
96 else {
97 ATH_MSG_WARNING("TRT Hit collection not found!");
98 }
99
100 return StatusCode::SUCCESS;
101}
#define ATH_MSG_WARNING(x)
AtlasHitsVector< TRTUncompressedHit > TRTUncompressedHitCollection
int barcode(const T *p)
Definition Barcode.h:16
@ u
Enums for curvilinear frames.
Definition ParamDefs.h:77

Member Data Documentation

◆ m_collection

std::string TrtHitsTestTool::m_collection
private

Definition at line 25 of file TrtHitsTestTool.h.

◆ m_detMgr

const InDetDD::TRT_DetectorManager* TrtHitsTestTool::m_detMgr {nullptr}
private

Definition at line 34 of file TrtHitsTestTool.h.

34{nullptr};

◆ m_histSvc

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

Definition at line 35 of file SimTestHisto.h.

35{"THistSvc", "SimTestHisto"};

◆ m_hits_edep_nonphotons

TH1 * TrtHitsTestTool::m_hits_edep_nonphotons
private

Definition at line 30 of file TrtHitsTestTool.h.

◆ m_hits_edep_photons

TH1 * TrtHitsTestTool::m_hits_edep_photons
private

Definition at line 29 of file TrtHitsTestTool.h.

◆ m_hits_edep_zr_nonphotons

TH2 * TrtHitsTestTool::m_hits_edep_zr_nonphotons
private

Definition at line 32 of file TrtHitsTestTool.h.

◆ m_hits_edep_zr_photons

TH2* TrtHitsTestTool::m_hits_edep_zr_photons
private

Definition at line 32 of file TrtHitsTestTool.h.

◆ m_hits_log_barcode

TH1* TrtHitsTestTool::m_hits_log_barcode
private

Definition at line 33 of file TrtHitsTestTool.h.

◆ m_hits_time_nonphotons

TH1* TrtHitsTestTool::m_hits_time_nonphotons
private

Definition at line 30 of file TrtHitsTestTool.h.

◆ m_hits_time_photons

TH1* TrtHitsTestTool::m_hits_time_photons
private

Definition at line 29 of file TrtHitsTestTool.h.

◆ m_hits_xy

TH2* TrtHitsTestTool::m_hits_xy
private

Definition at line 31 of file TrtHitsTestTool.h.

◆ m_hits_zr

TH2 * TrtHitsTestTool::m_hits_zr
private

Definition at line 31 of file TrtHitsTestTool.h.

◆ m_indetBarrel

TH2* TrtHitsTestTool::m_indetBarrel
private

Definition at line 27 of file TrtHitsTestTool.h.

◆ m_indetLongView

TH2 * TrtHitsTestTool::m_indetLongView
private

Definition at line 27 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.

34{"/truth/"};

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