LucidHitsTestTool Class Reference

#include <LucidHitsTestTool.h>

Inheritance diagram for LucidHitsTestTool:

Collaboration diagram for LucidHitsTestTool:

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

Protected Attributes
TH2 *	m_x_v_y_edep
TH2 *	m_x_v_y_hits
TH2 *	m_R_v_z_edep
TH2 *	m_R_v_z_hits
TH2 *	m_x_v_y_post
TH2 *	m_R_v_z_post
TH2 *	m_time_edep
TH1 *	m_pdgid
TH1 *	m_tubeid
TH1 *	m_gen_volume
TH1 *	m_wavelength
TH1 *	m_hit_edep
TH1 *	m_total_e
std::string	m_path {"/truth/"}
ServiceHandle< ITHistSvc >	m_histSvc {"THistSvc", "SimTestHisto"}

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

Detailed Description

Definition at line 10 of file LucidHitsTestTool.h.

Constructor & Destructor Documentation

LucidHitsTestTool()

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

Definition at line 15 of file LucidHitsTestTool.cxx.

  : SimTestToolBase(type, name, parent)
  , m_x_v_y_edep(0), m_x_v_y_hits(0), m_R_v_z_edep(0), m_R_v_z_hits(0)
  , m_x_v_y_post(0), m_R_v_z_post(0), m_time_edep(0)
  , m_pdgid(0), m_tubeid(0), m_gen_volume(0), m_wavelength(0), m_hit_edep(0), m_total_e(0)
{;}

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

virtual

Reimplemented from SimTestToolBase.

Definition at line 22 of file LucidHitsTestTool.cxx.

{
  m_path+="Lucid/";
  // ----------------------------------------
  _TH2D_WEIGHTED(m_x_v_y_edep,"x_v_y_edep",100,-500.,500.,100,-500.,500.);
  _SET_TITLE(m_x_v_y_edep, "Energy distribution [MeV]","x [mm]","y [mm]");
  _TH2D(m_x_v_y_hits,"x_v_y_hits",100,-500.,500.,100,-500.,500.);
  _SET_TITLE(m_x_v_y_hits, "Hit distribution","x [mm]","y [mm]");
 
  _TH2D_WEIGHTED(m_R_v_z_edep,"R_v_z_edep",100,0.,500.,100,5500.,8500.);
  _SET_TITLE(m_R_v_z_edep, "Energy distribution [MeV]","R [mm]","z [mm]");
  _TH2D(m_R_v_z_hits,"R_v_z_hits",100,0.,500.,100,5500.,8500.);
  _SET_TITLE(m_R_v_z_hits, "Hit distribution","R [mm]","z [mm]");
 
  _TH2D(m_x_v_y_post,"x_v_y_post",100,-500.,500.,100,-500.,500.);
  _SET_TITLE(m_x_v_y_post, "Post hit distribution","x [mm]","y [mm]");
  _TH2D(m_R_v_z_post,"R_v_z_post",100,0.,500.,100,5500.,8500.);
  _SET_TITLE(m_R_v_z_post, "Post hit distribution","R [mm]","z [mm]");
 
  _TH1D(m_pdgid,"pdg_id",1000,-3000,3000.);
  _SET_TITLE(m_pdgid, "PDG ID","PDG ID","Hits");
 
  _TH2D_WEIGHTED(m_time_edep,"Time_edep",100,-10.,50.,100,-10.,50.);
  _SET_TITLE(m_time_edep, "Energy distribution [MeV]","Time 1 [ns]","Time 2 [ns]");
 
  _TH1D(m_tubeid,"tube_id",1000,0,10000.);
  _SET_TITLE(m_tubeid, "Tube ID","Tube ID","Hits");
 
  _TH1D(m_gen_volume,"gen_volume",1000,0,100000.);
  _SET_TITLE(m_gen_volume, "Volume","Volume","Hits");
 
  _TH1D(m_wavelength,"wavelength",1000,0,3000.);
  _SET_TITLE(m_wavelength, "Wavelength","Wavelength","Hits");
 
  _TH1D(m_hit_edep,"hit_edep",100,0,10000.);
  _SET_TITLE(m_hit_edep, "Energy deposited per hit","Energy [MeV]","Hits");
 
  _TH1D(m_total_e,"total_e",1000,0,100000.);
  _SET_TITLE(m_total_e, "Total energy deposited","Energy [MeV]","Events");
 
  return StatusCode::SUCCESS;
}

processEvent()

StatusCode LucidHitsTestTool::processEvent

(

)

Definition at line 65 of file LucidHitsTestTool.cxx.

                                           {
 
  double etot = 0;
  const LUCID_SimHitCollection* iter = nullptr;
  CHECK(evtStore()->retrieve(iter));
  for (const LUCID_SimHit& hit : *iter) {
    m_x_v_y_edep->Fill( hit.GetX() , hit.GetY() , hit.GetEnergy() );
    m_x_v_y_hits->Fill( hit.GetX() , hit.GetY() );
    m_R_v_z_edep->Fill( std::sqrt( hit.GetX()*hit.GetX() + hit.GetY()*hit.GetY() ) , std::abs( hit.GetZ() ) , hit.GetEnergy() );
    m_R_v_z_hits->Fill( std::sqrt( hit.GetX()*hit.GetX() + hit.GetY()*hit.GetY() ) , std::abs( hit.GetZ() ) );
    m_x_v_y_post->Fill( hit.GetEPX() , hit.GetEPY() );
    m_R_v_z_post->Fill( std::sqrt( hit.GetEPX()*hit.GetEPX() + hit.GetEPY()*hit.GetEPY() ) , std::abs( hit.GetEPZ() ) );
 
    m_pdgid->Fill( hit.GetPdgCode() );
    m_time_edep->Fill( hit.GetPreStepTime() , hit.GetPostStepTime() , hit.GetEnergy() );
    m_tubeid->Fill( hit.GetTubeID() );
    m_gen_volume->Fill( hit.GetGenVolume() );
    m_wavelength->Fill( hit.GetWavelength() );
    m_hit_edep->Fill( hit.GetEnergy() );
 
    etot+=hit.GetEnergy();
  }
  m_total_e->Fill( etot );
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_gen_volume

TH1 * LucidHitsTestTool::m_gen_volume

protected

Definition at line 21 of file LucidHitsTestTool.h.

m_histSvc

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

protectedinherited

Definition at line 35 of file SimTestHisto.h.

{"THistSvc", "SimTestHisto"};

m_hit_edep

TH1 * LucidHitsTestTool::m_hit_edep

protected

Definition at line 21 of file LucidHitsTestTool.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/"};

m_pdgid

TH1* LucidHitsTestTool::m_pdgid

protected

Definition at line 21 of file LucidHitsTestTool.h.

m_R_v_z_edep

TH2 * LucidHitsTestTool::m_R_v_z_edep

protected

Definition at line 19 of file LucidHitsTestTool.h.

m_R_v_z_hits

TH2 * LucidHitsTestTool::m_R_v_z_hits

protected

Definition at line 19 of file LucidHitsTestTool.h.

m_R_v_z_post

TH2 * LucidHitsTestTool::m_R_v_z_post

protected

Definition at line 20 of file LucidHitsTestTool.h.

m_time_edep

TH2 * LucidHitsTestTool::m_time_edep

protected

Definition at line 20 of file LucidHitsTestTool.h.

m_total_e

TH1 * LucidHitsTestTool::m_total_e

protected

Definition at line 21 of file LucidHitsTestTool.h.

m_tubeid

TH1 * LucidHitsTestTool::m_tubeid

protected

Definition at line 21 of file LucidHitsTestTool.h.

m_wavelength

TH1 * LucidHitsTestTool::m_wavelength

protected

Definition at line 21 of file LucidHitsTestTool.h.

m_x_v_y_edep

TH2* LucidHitsTestTool::m_x_v_y_edep

protected

Definition at line 19 of file LucidHitsTestTool.h.

m_x_v_y_hits

TH2 * LucidHitsTestTool::m_x_v_y_hits

protected

Definition at line 19 of file LucidHitsTestTool.h.

m_x_v_y_post

TH2* LucidHitsTestTool::m_x_v_y_post

protected

Definition at line 20 of file LucidHitsTestTool.h.

---

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

• LucidHitsTestTool.h
• LucidHitsTestTool.cxx