SiHitsTestTool Class Reference

#include <SiHitsTestTool.h>

Inheritance diagram for SiHitsTestTool:

Collaboration diagram for SiHitsTestTool:

Public Member Functions
	SiHitsTestTool (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
TH1 *	m_hits_edep
TH1 *	m_hits_log_edep
TH2 *	m_hits_xy
TH2 *	m_hits_zr
TH2 *	m_hits_edep_zr
TH1 *	m_hits_edep_z
TH1 *	m_hits_edep_r
TH1 *	m_hits_log_barcode
TH1 *	m_hits_step_length
TH1 *	m_hits_log_step_length

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

Detailed Description

Definition at line 11 of file SiHitsTestTool.h.

Constructor & Destructor Documentation

SiHitsTestTool()

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

Definition at line 16 of file SiHitsTestTool.cxx.

  : SimTestToolBase(type, name, parent),
    m_collection("PixelHits"),
    m_indetBarrel(0), m_indetLongView(0),
    m_hits_time(0), m_hits_edep(0), m_hits_log_edep(0),
    m_hits_xy(0), m_hits_zr(0),
    m_hits_edep_zr(0),
    m_hits_edep_z(0), m_hits_edep_r(0),
    m_hits_log_barcode(0),
    m_hits_step_length(0),
    m_hits_log_step_length(0)
{  
    declareProperty("CollectionName",  m_collection="PixelHits");
}

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

virtual

Reimplemented from SimTestToolBase.

Definition at line 31 of file SiHitsTestTool.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","x [mm]","y [mm]");
  // perhaps better make r-phi?
  _TH2D(m_indetLongView,"idet_zr",100,-3200.,3200.,100,0.,1200.);
  _SET_TITLE(m_indetLongView, "indet hit distribution","z [mm]","r [mm]");
 
  std::string detName("Pixel"); std::string detpathName("Pixel");
  // initialise pixel or SCT variables:
  if (m_collection=="PixelHits") {
    detName = "Pixel"; detpathName = "Pixel";
  }
  else if (m_collection=="PileupPixelHits") {
    detName = "PileupPixel"; detpathName = "Pixel";
  }
  else if (m_collection=="SCT_Hits") {
    detName = "SCT"; detpathName = "SCT";
  }
  else if (m_collection=="PileupSCT_Hits") {
    detName = "PileupSCT"; detpathName = "SCT";  
  }
  else {
    ATH_MSG_ERROR("collection: "<<m_collection<<", SiHitsTestTool for "<<name()<<" not supported !!!\n");
    return StatusCode::FAILURE;
  }
  m_path+=detpathName+"/";
  _TH2D(m_hits_xy,(detName+"_xy").c_str(),100,-1200.,1200.,100,-1200.,1200.);
  _SET_TITLE(m_hits_xy, "hit distribution","x [mm]","y [mm]");
  _TH2D(m_hits_zr,(detName+"_zr").c_str(),100,-3200.,3200.,100,0.,1200.);
  _SET_TITLE(m_hits_zr, "hit distribution","z [mm]","r [mm]");
  _TH1D_WEIGHTED(m_hits_time,(detName+"_edep_time").c_str(),100,0, 25);
  _SET_TITLE(m_hits_time, "energy weighted hit distribution","t [ns]","dE/dt");
  _TH1D(m_hits_edep,(detName+"_edep").c_str(),100,0.,1.);
  _SET_TITLE(m_hits_edep, "hit distribution","edep [MeV]","dN/dE [1/MeV]");
  _TH1D(m_hits_log_edep,(detName+"_log_edep").c_str(),100,-15.,2.);
  _SET_TITLE(m_hits_log_edep, "log hit distribution","log(edep [MeV])","dN/dlog(E) [1/log(MeV)]");
 
  _TH2D_WEIGHTED(m_hits_edep_zr,(detName+"_edep_zr").c_str(),100,-3200.,3200.,100,0.,1200.);
  _SET_TITLE(m_hits_edep_zr, "energy weighted indet hit distribution","z [mm]","r [mm]");
 
  _TH1D_WEIGHTED(m_hits_edep_z,(detName+"_edep_z").c_str(),100,-3200.,3200.);
  _SET_TITLE(m_hits_edep_z, "energy weighted indet hit distribution (barrel only)","z [mm]","dE/dz [MeV/mm]");
  _TH1D_WEIGHTED(m_hits_edep_r,(detName+"_edep_r").c_str(),100,0.,1200.);
  _SET_TITLE(m_hits_edep_r, "energy weighted indet hit distribution (endcap only)","r [mm]","dE/dr [MeV/mm]");
 
  _TH1D(m_hits_log_barcode,(detName+"_log_barcode").c_str(),50,0.,16.);
  _SET_TITLE(m_hits_log_barcode, "log(barcode)","log(barcode)","dN/dlog(barcode)");
 
  _TH1D(m_hits_step_length,(detName+"_step_length").c_str(),100,0.,2.);
  _SET_TITLE(m_hits_step_length, "Step Length","Step Length l","dN/dl");
  _TH1D(m_hits_log_step_length,(detName+"_log_step_length").c_str(), 100, -13, 2.);
  _SET_TITLE(m_hits_log_step_length, "Logarithm of Step Length","log(l [mm])","dN/dlog(l)");
 
  return StatusCode::SUCCESS;
}

processEvent()

StatusCode SiHitsTestTool::processEvent

(

)

Definition at line 90 of file SiHitsTestTool.cxx.

                                        {
 
  const SiHitCollection* p_collection = nullptr;
  if (evtStore()->retrieve(p_collection,m_collection).isSuccess())  {
    for (const SiHit& hit : *p_collection) {
      GeoSiHit ghit(hit);
      HepGeom::Point3D<double> u = ghit.getGlobalPosition();
      if (m_hits_xy) {
    m_hits_xy->Fill(u.x(),u.y());
    m_hits_zr->Fill(u.z(),u.perp());
      }
      m_indetBarrel->Fill(u.x(),u.y());
      m_indetLongView->Fill(u.z(),u.perp());
 
      m_hits_time->Fill(hit.meanTime(),hit.energyLoss());
      m_hits_edep->Fill(hit.energyLoss());
      m_hits_log_edep->Fill( log(hit.energyLoss()) );
      m_hits_edep_zr->Fill(u.z(),u.perp(),hit.energyLoss());
      int barcode = HepMC::barcode(hit.particleLink()); // FIXME barcode-based
      m_hits_log_barcode->Fill( barcode > 0 ? log(barcode) : -1 );
      double step_length = ( hit.localStartPosition()-hit.localEndPosition() ).mag();
      m_hits_step_length->Fill(step_length);
      m_hits_log_step_length->Fill( log(step_length));
 
      if (hit.getBarrelEndcap()==0) {
    // barrel
    m_hits_edep_z->Fill(u.z(),hit.energyLoss());
      }
      else {
    // discs
    m_hits_edep_r->Fill(u.perp(),hit.energyLoss());
      }
 
      // debuging only:
      //hit.print();
      /*
          *** SCT Hit 
          Barrel/ EndCap Number -2
          Layer/Disk Number     5
          Eta Number            0
          Phi Number            43
          Side Number            1
       */
 
    }
  }
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_collection

std::string SiHitsTestTool::m_collection

private

Definition at line 23 of file SiHitsTestTool.h.

m_histSvc

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

protectedinherited

Definition at line 35 of file SimTestHisto.h.

{"THistSvc", "SimTestHisto"};

m_hits_edep

TH1 * SiHitsTestTool::m_hits_edep

private

Definition at line 27 of file SiHitsTestTool.h.

m_hits_edep_r

TH1 * SiHitsTestTool::m_hits_edep_r

private

Definition at line 30 of file SiHitsTestTool.h.

m_hits_edep_z

TH1* SiHitsTestTool::m_hits_edep_z

private

Definition at line 30 of file SiHitsTestTool.h.

m_hits_edep_zr

TH2* SiHitsTestTool::m_hits_edep_zr

private

Definition at line 29 of file SiHitsTestTool.h.

m_hits_log_barcode

TH1* SiHitsTestTool::m_hits_log_barcode

private

Definition at line 31 of file SiHitsTestTool.h.

m_hits_log_edep

TH1 * SiHitsTestTool::m_hits_log_edep

private

Definition at line 27 of file SiHitsTestTool.h.

m_hits_log_step_length

TH1 * SiHitsTestTool::m_hits_log_step_length

private

Definition at line 32 of file SiHitsTestTool.h.

m_hits_step_length

TH1* SiHitsTestTool::m_hits_step_length

private

Definition at line 32 of file SiHitsTestTool.h.

m_hits_time

TH1* SiHitsTestTool::m_hits_time

private

Definition at line 27 of file SiHitsTestTool.h.

m_hits_xy

TH2* SiHitsTestTool::m_hits_xy

private

Definition at line 28 of file SiHitsTestTool.h.

m_hits_zr

TH2 * SiHitsTestTool::m_hits_zr

private

Definition at line 28 of file SiHitsTestTool.h.

m_indetBarrel

TH2* SiHitsTestTool::m_indetBarrel

private

Definition at line 25 of file SiHitsTestTool.h.

m_indetLongView

TH2 * SiHitsTestTool::m_indetLongView

private

Definition at line 25 of file SiHitsTestTool.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.

{"/truth/"};

---

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

• SiHitsTestTool.h
• SiHitsTestTool.cxx