#include <MaterialAccountingGroup.h>

Classes
class	BoundingBox

Public Member Functions
bool	addDetector (const MaterialAccountingDetector &detector)
	buffer material from a detector, if the detector is inside the DetLayer bounds More...

MaterialAccountingStep	average (void) const
	return the average normalized material accounting informations More...

double	averageEnergyLoss (void) const
	return the average normalized energy loss density factor for Bethe-Bloch More...

double	averageLength (void) const
	return the average normalized layer thickness More...

double	averageRadiationLengths (void) const
	return the average normalized number of radiation lengths More...

const std::vector< GlobalPoint > &	elements (void) const

void	endOfTrack (void)
	commit the buffer and reset the "already hit by this track" flag More...

std::pair< double, double >	getBoundingR () const
	Return the bouding limit in R for the hosted Group. More...

std::pair< double, double >	getBoundingZ () const
	Return the bouding limit in Z for the hosted Group. More...

std::string	info (void) const
	get some infos More...

bool	inside (const MaterialAccountingDetector &detector) const
	check if detector is inside any part of this layer More...

	MaterialAccountingGroup (const std::string &name, const DDCompactView &geometry)
	explicit constructors More...

const std::string &	name (void) const
	get the layer name More...

void	savePlots (void)
	save the plots More...

double	sigmaEnergyLoss (void) const
	return the sigma of the normalized energy loss density factor for Bethe-Bloch More...

double	sigmaLength (void) const
	return the sigma of the normalized layer thickness More...

double	sigmaRadiationLengths (void) const
	return the sigma of the normalized number of radiation lengths More...

unsigned int	tracks (void) const
	return the number of tracks that hit this layer More...

	~MaterialAccountingGroup (void)
	destructor More...

Private Member Functions
	MaterialAccountingGroup (const MaterialAccountingGroup &layer)=delete
	stop default copy ctor More...

MaterialAccountingGroup &	operator= (const MaterialAccountingGroup &layer)=delete
	stop default assignment operator More...

void	savePlot (TH1F *plot, const std::string &name)

void	savePlot (TProfile *plot, float average, const std::string &name)

Private Attributes
MaterialAccountingStep	m_accounting

BoundingBox	m_boundingbox

MaterialAccountingStep	m_buffer

bool	m_counted

TH1F *	m_dedx_spectrum

TProfile *	m_dedx_vs_eta

TProfile *	m_dedx_vs_r

TProfile *	m_dedx_vs_z

std::vector< GlobalPoint >	m_elements

MaterialAccountingStep	m_errors

TFile *	m_file

std::string	m_name

TH1F *	m_radlen_spectrum

TProfile *	m_radlen_vs_eta

TProfile *	m_radlen_vs_r

TProfile *	m_radlen_vs_z

unsigned int	m_tracks

Static Private Attributes
static const double	s_tolerance

Detailed Description

Definition at line 19 of file MaterialAccountingGroup.h.

Constructor & Destructor Documentation

◆ MaterialAccountingGroup() [1/2]

MaterialAccountingGroup::MaterialAccountingGroup	(	const std::string &	name,
		const DDCompactView &	geometry
	)

explicit constructors

Definition at line 23 of file MaterialAccountingGroup.cc.

     : m_name(name),
       m_elements(),
       m_boundingbox(),
       m_accounting(),
       m_errors(),
       m_tracks(0),
       m_counted(false),
       m_file(nullptr) {
   // retrieve the elements from DDD
   DDValue namevalue;
   if (TString(name.c_str()).Contains("Tracker")) {
     namevalue = DDValue("TrackingMaterialGroup", name);
   } else if (TString(name.c_str()).Contains("HGCal")) {
     namevalue = DDValue("Volume", name);
   } else if (TString(name.c_str()).Contains("HFNose")) {
     namevalue = DDValue("Volume", name);
   } else {
     LogTrace("MaterialAccountingGroup") << name << std::endl;
     edm::LogError("MaterialAccountingGroup") << "Only Tracker , HGCal and HFNose are supported" << std::endl;
   }
  
   DDSpecificsMatchesValueFilter filter{namevalue};
   DDFilteredView fv(geometry, filter);
   LogTrace("MaterialAccountingGroup") << "Elements within: " << name << std::endl;
   while (fv.next()) {
     // DD3Vector and DDTranslation are the same type as math::XYZVector
     math::XYZVector position = fv.translation() / 10.;  // mm -> cm
     LogTrace("MaterialAccountingGroup") << "Adding element at(r,z): ("
                                         << GlobalPoint(position.x(), position.y(), position.z()).perp() << ", "
                                         << GlobalPoint(position.x(), position.y(), position.z()).z() << ") cm"
                                         << std::endl;
     LogTrace("MaterialAccountingGroup") << "Name of added element: " << fv.logicalPart().toString() << std::endl;
     m_elements.push_back(GlobalPoint(position.x(), position.y(), position.z()));
   }
  
   // grow the bounding box
   for (unsigned int i = 0; i < m_elements.size(); ++i) {
     m_boundingbox.grow(m_elements[i].perp(), m_elements[i].z());
   }
   m_boundingbox.grow(s_tolerance);
   LogTrace("MaterialAccountingGroup") << "Final BBox r_range: " << m_boundingbox.range_r().first << ", "
                                       << m_boundingbox.range_r().second << std::endl
                                       << "Final BBox z_range: " << m_boundingbox.range_z().first << ", "
                                       << m_boundingbox.range_z().second << std::endl;
  
   // initialize the histograms
   m_dedx_spectrum = new TH1F((m_name + "_dedx_spectrum").c_str(), "Energy loss spectrum", 1000, 0, 1);
   m_radlen_spectrum = new TH1F((m_name + "_radlen_spectrum").c_str(), "Radiation lengths spectrum", 1000, 0, 1);
   m_dedx_vs_eta = new TProfile((m_name + "_dedx_vs_eta").c_str(), "Energy loss vs. eta", 600, -3, 3);
   m_dedx_vs_z = new TProfile((m_name + "_dedx_vs_z").c_str(), "Energy loss vs. Z", 6000, -300, 300);
   m_dedx_vs_r = new TProfile((m_name + "_dedx_vs_r").c_str(), "Energy loss vs. R", 1200, 0, 120);
   m_radlen_vs_eta = new TProfile((m_name + "_radlen_vs_eta").c_str(), "Radiation lengths vs. eta", 600, -3, 3);
   m_radlen_vs_z = new TProfile((m_name + "_radlen_vs_z").c_str(), "Radiation lengths vs. Z", 6000, -300, 300);
   m_radlen_vs_r = new TProfile((m_name + "_radlen_vs_r").c_str(), "Radiation lengths vs. R", 1200, 0, 120);
   m_dedx_spectrum->SetDirectory(nullptr);
   m_radlen_spectrum->SetDirectory(nullptr);
   m_dedx_vs_eta->SetDirectory(nullptr);
   m_dedx_vs_z->SetDirectory(nullptr);
   m_dedx_vs_r->SetDirectory(nullptr);
   m_radlen_vs_eta->SetDirectory(nullptr);
   m_radlen_vs_z->SetDirectory(nullptr);
   m_radlen_vs_r->SetDirectory(nullptr);
 }

References ALCARECOTkAlBeamHalo_cff::filter, MaterialAccountingGroup::BoundingBox::grow(), mps_fire::i, DDFilteredView::logicalPart(), LogTrace, m_boundingbox, m_dedx_spectrum, m_dedx_vs_eta, m_dedx_vs_r, m_dedx_vs_z, m_elements, m_name, m_radlen_spectrum, m_radlen_vs_eta, m_radlen_vs_r, m_radlen_vs_z, name(), DDFilteredView::next(), perp(), position, MaterialAccountingGroup::BoundingBox::range_r(), MaterialAccountingGroup::BoundingBox::range_z(), s_tolerance, DDBase< N, C >::toString(), DDFilteredView::translation(), and z.

◆ ~MaterialAccountingGroup()

MaterialAccountingGroup::~MaterialAccountingGroup ( void )

destructor

Definition at line 88 of file MaterialAccountingGroup.cc.

                                                       {
   delete m_dedx_spectrum;
   delete m_dedx_vs_eta;
   delete m_dedx_vs_z;
   delete m_dedx_vs_r;
   delete m_radlen_spectrum;
   delete m_radlen_vs_eta;
   delete m_radlen_vs_z;
   delete m_radlen_vs_r;
 }

References m_dedx_spectrum, m_dedx_vs_eta, m_dedx_vs_r, m_dedx_vs_z, m_radlen_spectrum, m_radlen_vs_eta, m_radlen_vs_r, and m_radlen_vs_z.

◆ MaterialAccountingGroup() [2/2]

MaterialAccountingGroup::MaterialAccountingGroup ( const MaterialAccountingGroup & layer )

privatedelete

stop default copy ctor

Member Function Documentation

◆ addDetector()

bool MaterialAccountingGroup::addDetector ( const MaterialAccountingDetector & detector )

buffer material from a detector, if the detector is inside the DetLayer bounds

Definition at line 132 of file MaterialAccountingGroup.cc.

                                                                                     {
   if (not inside(detector))
     return false;
  
   // multiple hits in the same layer (overlaps, etc.) from a single track still count as one for averaging,
   // since the energy deposits from the track have been already split between the different detectors
   m_buffer += detector.material();
   m_counted = true;
  
   return true;
 }

References hgcalTestNeighbor_cfi::detector, inside(), m_buffer, and m_counted.

◆ average()

MaterialAccountingStep MaterialAccountingGroup::average ( void ) const

inline

return the average normalized material accounting informations

Definition at line 91 of file MaterialAccountingGroup.h.

91 { return m_tracks ? m_accounting / m_tracks : MaterialAccountingStep(); }

References m_accounting, and m_tracks.

Referenced by endOfTrack().

◆ averageEnergyLoss()

double MaterialAccountingGroup::averageEnergyLoss ( void ) const

inline

return the average normalized energy loss density factor for Bethe-Bloch

Definition at line 100 of file MaterialAccountingGroup.h.

100 { return m_tracks ? m_accounting.energyLoss() / m_tracks : 0.; }

References MaterialAccountingStep::energyLoss(), m_accounting, and m_tracks.

Referenced by savePlots(), and sigmaEnergyLoss().

◆ averageLength()

double MaterialAccountingGroup::averageLength ( void ) const

inline

return the average normalized layer thickness

Definition at line 94 of file MaterialAccountingGroup.h.

94 { return m_tracks ? m_accounting.length() / m_tracks : 0.; }

References MaterialAccountingStep::length(), m_accounting, and m_tracks.

Referenced by sigmaLength().

◆ averageRadiationLengths()

double MaterialAccountingGroup::averageRadiationLengths ( void ) const

inline

return the average normalized number of radiation lengths

Definition at line 97 of file MaterialAccountingGroup.h.

97 { return m_tracks ? m_accounting.radiationLengths() / m_tracks : 0.; }

References m_accounting, m_tracks, and MaterialAccountingStep::radiationLengths().

Referenced by savePlots(), and sigmaRadiationLengths().

◆ elements()

const std::vector<GlobalPoint>& MaterialAccountingGroup::elements ( void ) const

inline

Definition at line 135 of file MaterialAccountingGroup.h.

135 { return m_elements; }

References m_elements.

◆ endOfTrack()

void MaterialAccountingGroup::endOfTrack ( void )

commit the buffer and reset the "already hit by this track" flag

Definition at line 144 of file MaterialAccountingGroup.cc.

                                              {
   // add a detector
   if (m_counted) {
     m_accounting += m_buffer;
     m_errors += m_buffer * m_buffer;
     ++m_tracks;
  
     GlobalPoint average((m_buffer.in().x() + m_buffer.out().x()) / 2.,
                         (m_buffer.in().y() + m_buffer.out().y()) / 2.,
                         (m_buffer.in().z() + m_buffer.out().z()) / 2.);
     m_dedx_spectrum->Fill(m_buffer.energyLoss());
     m_radlen_spectrum->Fill(m_buffer.radiationLengths());
     m_dedx_vs_eta->Fill(average.eta(), m_buffer.energyLoss(), 1.);
     m_dedx_vs_z->Fill(average.z(), m_buffer.energyLoss(), 1.);
     m_dedx_vs_r->Fill(average.perp(), m_buffer.energyLoss(), 1.);
     m_radlen_vs_eta->Fill(average.eta(), m_buffer.radiationLengths(), 1.);
     m_radlen_vs_z->Fill(average.z(), m_buffer.radiationLengths(), 1.);
     m_radlen_vs_r->Fill(average.perp(), m_buffer.radiationLengths(), 1.);
   }
   m_counted = false;
   m_buffer = MaterialAccountingStep();
 }

References average(), MaterialAccountingStep::energyLoss(), MaterialAccountingStep::in(), m_accounting, m_buffer, m_counted, m_dedx_spectrum, m_dedx_vs_eta, m_dedx_vs_r, m_dedx_vs_z, m_errors, m_radlen_spectrum, m_radlen_vs_eta, m_radlen_vs_r, m_radlen_vs_z, m_tracks, MaterialAccountingStep::out(), MaterialAccountingStep::radiationLengths(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ getBoundingR()

std::pair<double, double> MaterialAccountingGroup::getBoundingR ( ) const

inline

Return the bouding limit in R for the hosted Group.

Definition at line 85 of file MaterialAccountingGroup.h.

85 { return m_boundingbox.range_r(); };

References m_boundingbox, and MaterialAccountingGroup::BoundingBox::range_r().

◆ getBoundingZ()

std::pair<double, double> MaterialAccountingGroup::getBoundingZ ( ) const

inline

Return the bouding limit in Z for the hosted Group.

Definition at line 88 of file MaterialAccountingGroup.h.

88 { return m_boundingbox.range_z(); };

References m_boundingbox, and MaterialAccountingGroup::BoundingBox::range_z().

◆ info()

std::string MaterialAccountingGroup::info ( void ) const

get some infos

Definition at line 208 of file MaterialAccountingGroup.cc.

                                                   {
   std::stringstream out;
   out << std::setw(48) << std::left << m_name << std::right << std::fixed;
   ;
   out << "BBox: " << std::setprecision(1) << std::setw(6) << m_boundingbox.range_z().first << " < Z < "
       << std::setprecision(1) << std::setw(6) << m_boundingbox.range_z().second;
   out << ", " << std::setprecision(1) << std::setw(5) << m_boundingbox.range_r().first << " < R < "
       << std::setprecision(1) << std::setw(5) << m_boundingbox.range_r().second;
   out << "   Elements: " << std::setw(6) << m_elements.size();
   return out.str();
 }

References alignBH_cfg::fixed, m_boundingbox, m_elements, m_name, MillePedeFileConverter_cfg::out, MaterialAccountingGroup::BoundingBox::range_r(), and MaterialAccountingGroup::BoundingBox::range_z().

◆ inside()

bool MaterialAccountingGroup::inside ( const MaterialAccountingDetector & detector ) const

check if detector is inside any part of this layer

Definition at line 103 of file MaterialAccountingGroup.cc.

                                                                                      {
   const GlobalPoint& position = detector.position();
   // first check to see if the point is inside the bounding box
   LogTrace("MaterialAccountingGroup") << "Testing position: (x, y, z, r) = " << position.x() << ", " << position.y()
                                       << ", " << position.z() << ", " << position.perp() << std::endl;
   if (not m_boundingbox.inside(position.perp(), position.z())) {
     LogTrace("MaterialAccountingGroup") << "r outside of: (" << m_boundingbox.range_r().first << ", "
                                         << m_boundingbox.range_r().second << "), Z ouside of: ("
                                         << m_boundingbox.range_z().first << ", " << m_boundingbox.range_z().second
                                         << ")" << std::endl;
     return false;
   } else {
     // now check if the point is actually close enough to any element
     LogTrace("MaterialAccountingGroup") << "r within: (" << m_boundingbox.range_r().first << ", "
                                         << m_boundingbox.range_r().second << "), Z within: ("
                                         << m_boundingbox.range_z().first << ", " << m_boundingbox.range_z().second
                                         << ")" << std::endl;
     for (unsigned int i = 0; i < m_elements.size(); ++i) {
       LogTrace("MaterialAccountingGroup")
           << "Closest testing agains(x, y, z, r): (" << m_elements[i].x() << ", " << m_elements[i].y() << ", "
           << m_elements[i].z() << ", " << m_elements[i].perp() << ") --> " << (position - m_elements[i]).mag()
           << " vs tolerance: " << s_tolerance << std::endl;
       if ((position - m_elements[i]).mag2() < (s_tolerance * s_tolerance))
         return true;
     }
     return false;
   }
 }

References hgcalTestNeighbor_cfi::detector, mps_fire::i, MaterialAccountingGroup::BoundingBox::inside(), LogTrace, m_boundingbox, m_elements, mag(), position, MaterialAccountingGroup::BoundingBox::range_r(), MaterialAccountingGroup::BoundingBox::range_z(), and s_tolerance.

Referenced by addDetector().

◆ name()

const std::string& MaterialAccountingGroup::name ( void ) const

inline

get the layer name

Definition at line 123 of file MaterialAccountingGroup.h.

123 { return m_name; }

References m_name.

Referenced by config.CFG::__str__(), validation.Sample::digest(), MaterialAccountingGroup(), and savePlot().

◆ operator=()

MaterialAccountingGroup& MaterialAccountingGroup::operator= ( const MaterialAccountingGroup & layer )

privatedelete

stop default assignment operator

◆ savePlot() [1/2]

void MaterialAccountingGroup::savePlot	(	TH1F *	plot,
		const std::string &	name
	)

private

Definition at line 167 of file MaterialAccountingGroup.cc.

                                                                         {
   TCanvas canvas(name.c_str(), plot->GetTitle(), 1280, 1024);
   plot->SetFillColor(15);  // grey
   plot->SetLineColor(1);   // black
   plot->Draw("c e");
   canvas.GetFrame()->SetFillColor(kWhite);
   canvas.Draw();
   canvas.SaveAs((name + ".png").c_str(), "");
  
   // store te plot into m_file
   plot->SetDirectory(m_file);
 }

References svgfig::canvas(), m_file, name(), and plotFactory::plot.

Referenced by savePlots().

◆ savePlot() [2/2]

void MaterialAccountingGroup::savePlot	(	TProfile *	plot,
		float	average,
		const std::string &	name
	)

private

Definition at line 180 of file MaterialAccountingGroup.cc.

                                                                                            {
   // Nota Bene:
   // these "line" plots are not deleted explicitly since
   //   - deleting them before saving them to a TFile will not save them
   //   - deleting them after the TFile they're stored into results in a SEGV
   // ROOT is probably "taking care" (read: messing things up) somehow...
   TH1F* line =
       new TH1F((name + "_par").c_str(), "Parametrization", 1, plot->GetXaxis()->GetXmin(), plot->GetXaxis()->GetXmax());
   line->SetBinContent(1, average);
  
   TCanvas canvas(name.c_str(), plot->GetTitle(), 1280, 1024);
   plot->SetFillColor(15);  // grey
   plot->SetLineColor(1);   // black
   plot->SetLineWidth(2);
   plot->Draw("c e6");
   line->SetLineColor(2);  // red
   line->SetLineWidth(2);
   line->Draw("same");
   canvas.GetFrame()->SetFillColor(kWhite);
   canvas.Draw();
   canvas.SaveAs((name + ".png").c_str(), "");
  
   // store te plots into m_file
   plot->SetDirectory(m_file);
   line->SetDirectory(m_file);
 }

References svgfig::canvas(), mps_splice::line, m_file, name(), and plotFactory::plot.

◆ savePlots()

void MaterialAccountingGroup::savePlots ( void )

save the plots

Definition at line 220 of file MaterialAccountingGroup.cc.

                                             {
   m_file = new TFile((m_name + ".root").c_str(), "RECREATE");
   savePlot(m_dedx_spectrum, m_name + "_dedx_spectrum");
   savePlot(m_radlen_spectrum, m_name + "_radlen_spectrum");
   savePlot(m_dedx_vs_eta, averageEnergyLoss(), m_name + "_dedx_vs_eta");
   savePlot(m_dedx_vs_z, averageEnergyLoss(), m_name + "_dedx_vs_z");
   savePlot(m_dedx_vs_r, averageEnergyLoss(), m_name + "_dedx_vs_r");
   savePlot(m_radlen_vs_eta, averageRadiationLengths(), m_name + "_radlen_vs_eta");
   savePlot(m_radlen_vs_z, averageRadiationLengths(), m_name + "_radlen_vs_z");
   savePlot(m_radlen_vs_r, averageRadiationLengths(), m_name + "_radlen_vs_r");
   m_file->Write();
   m_file->Close();
  
   delete m_file;
 }

References averageEnergyLoss(), averageRadiationLengths(), m_dedx_spectrum, m_dedx_vs_eta, m_dedx_vs_r, m_dedx_vs_z, m_file, m_name, m_radlen_spectrum, m_radlen_vs_eta, m_radlen_vs_r, m_radlen_vs_z, and savePlot().

◆ sigmaEnergyLoss()

double MaterialAccountingGroup::sigmaEnergyLoss ( void ) const

inline

return the sigma of the normalized energy loss density factor for Bethe-Bloch

Definition at line 115 of file MaterialAccountingGroup.h.

                                      {
     return m_tracks ? std::sqrt(m_errors.energyLoss() / m_tracks - averageEnergyLoss() * averageEnergyLoss()) : 0.;
   }

References averageEnergyLoss(), MaterialAccountingStep::energyLoss(), m_errors, m_tracks, and mathSSE::sqrt().

◆ sigmaLength()

double MaterialAccountingGroup::sigmaLength ( void ) const

inline

return the sigma of the normalized layer thickness

Definition at line 103 of file MaterialAccountingGroup.h.

                                  {
     return m_tracks ? std::sqrt(m_errors.length() / m_tracks - averageLength() * averageLength()) : 0.;
   }

References averageLength(), MaterialAccountingStep::length(), m_errors, m_tracks, and mathSSE::sqrt().

◆ sigmaRadiationLengths()

double MaterialAccountingGroup::sigmaRadiationLengths ( void ) const

inline

return the sigma of the normalized number of radiation lengths

Definition at line 108 of file MaterialAccountingGroup.h.

                                            {
     return m_tracks ? std::sqrt(m_errors.radiationLengths() / m_tracks -
                                 averageRadiationLengths() * averageRadiationLengths())
                     : 0.;
   }

References averageRadiationLengths(), m_errors, m_tracks, MaterialAccountingStep::radiationLengths(), and mathSSE::sqrt().

◆ tracks()

unsigned int MaterialAccountingGroup::tracks ( void ) const

inline

return the number of tracks that hit this layer

Definition at line 120 of file MaterialAccountingGroup.h.

120 { return m_tracks; }

References m_tracks.

Member Data Documentation

◆ m_accounting

MaterialAccountingStep MaterialAccountingGroup::m_accounting

private

Definition at line 144 of file MaterialAccountingGroup.h.

Referenced by average(), averageEnergyLoss(), averageLength(), averageRadiationLengths(), and endOfTrack().

◆ m_boundingbox

BoundingBox MaterialAccountingGroup::m_boundingbox

private

Definition at line 143 of file MaterialAccountingGroup.h.

Referenced by getBoundingR(), getBoundingZ(), info(), inside(), and MaterialAccountingGroup().

◆ m_buffer

MaterialAccountingStep MaterialAccountingGroup::m_buffer

private

Definition at line 148 of file MaterialAccountingGroup.h.

Referenced by addDetector(), and endOfTrack().

◆ m_counted

bool MaterialAccountingGroup::m_counted

private

Definition at line 147 of file MaterialAccountingGroup.h.

Referenced by addDetector(), and endOfTrack().

◆ m_dedx_spectrum

TH1F* MaterialAccountingGroup::m_dedx_spectrum

private

Definition at line 151 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_dedx_vs_eta

TProfile* MaterialAccountingGroup::m_dedx_vs_eta

private

Definition at line 154 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_dedx_vs_r

TProfile* MaterialAccountingGroup::m_dedx_vs_r

private

Definition at line 156 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_dedx_vs_z

TProfile* MaterialAccountingGroup::m_dedx_vs_z

private

Definition at line 155 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_elements

std::vector<GlobalPoint> MaterialAccountingGroup::m_elements

private

Definition at line 142 of file MaterialAccountingGroup.h.

Referenced by elements(), info(), inside(), and MaterialAccountingGroup().

◆ m_errors

MaterialAccountingStep MaterialAccountingGroup::m_errors

private

Definition at line 145 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), sigmaEnergyLoss(), sigmaLength(), and sigmaRadiationLengths().

◆ m_file

TFile* MaterialAccountingGroup::m_file

mutableprivate

Definition at line 162 of file MaterialAccountingGroup.h.

Referenced by savePlot(), and savePlots().

◆ m_name

std::string MaterialAccountingGroup::m_name

private

Definition at line 141 of file MaterialAccountingGroup.h.

Referenced by info(), MaterialAccountingGroup(), name(), and savePlots().

◆ m_radlen_spectrum

TH1F* MaterialAccountingGroup::m_radlen_spectrum

private

Definition at line 152 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_radlen_vs_eta

TProfile* MaterialAccountingGroup::m_radlen_vs_eta

private

Definition at line 157 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_radlen_vs_r

TProfile* MaterialAccountingGroup::m_radlen_vs_r

private

Definition at line 159 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_radlen_vs_z

TProfile* MaterialAccountingGroup::m_radlen_vs_z

private

Definition at line 158 of file MaterialAccountingGroup.h.

Referenced by endOfTrack(), MaterialAccountingGroup(), savePlots(), and ~MaterialAccountingGroup().

◆ m_tracks

unsigned int MaterialAccountingGroup::m_tracks

private

Definition at line 146 of file MaterialAccountingGroup.h.

Referenced by average(), averageEnergyLoss(), averageLength(), averageRadiationLengths(), endOfTrack(), sigmaEnergyLoss(), sigmaLength(), sigmaRadiationLengths(), and tracks().

◆ s_tolerance

const double MaterialAccountingGroup::s_tolerance

staticprivate

Initial value:

=

0.01

Definition at line 165 of file MaterialAccountingGroup.h.

Referenced by inside(), and MaterialAccountingGroup().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MaterialAccountingGroup() [1/2]

◆ ~MaterialAccountingGroup()

◆ MaterialAccountingGroup() [2/2]

Member Function Documentation

◆ addDetector()

◆ average()

◆ averageEnergyLoss()

◆ averageLength()

◆ averageRadiationLengths()

◆ elements()

◆ endOfTrack()

◆ getBoundingR()

◆ getBoundingZ()

◆ info()

◆ inside()

◆ name()

◆ operator=()

◆ savePlot() [1/2]

◆ savePlot() [2/2]

◆ savePlots()

◆ sigmaEnergyLoss()

◆ sigmaLength()

◆ sigmaRadiationLengths()

◆ tracks()

Member Data Documentation

◆ m_accounting

◆ m_boundingbox

◆ m_buffer

◆ m_counted

◆ m_dedx_spectrum

◆ m_dedx_vs_eta

◆ m_dedx_vs_r

◆ m_dedx_vs_z

◆ m_elements

◆ m_errors

◆ m_file

◆ m_name

◆ m_radlen_spectrum

◆ m_radlen_vs_eta

◆ m_radlen_vs_r

◆ m_radlen_vs_z

◆ m_tracks

◆ s_tolerance