#include <MaterialBudgetHcal.h>

Inheritance diagram for MaterialBudgetHcal:

Public Member Functions
void	beginRun (edm::EventSetup const &) override

	MaterialBudgetHcal (const edm::ParameterSet &)

	MaterialBudgetHcal (const MaterialBudgetHcal &)=delete

const MaterialBudgetHcal &	operator= (const MaterialBudgetHcal &)=delete

void	registerConsumes (edm::ConsumesCollector) override

Public Member Functions inherited from SimWatcher
bool	isMT () const

const SimWatcher &	operator= (const SimWatcher &)=delete

	SimWatcher ()

	SimWatcher (const SimWatcher &)=delete

virtual	~SimWatcher ()

Public Member Functions inherited from Observer< const BeginOfTrack *>
	Observer ()

void	slotForUpdate (const BeginOfTrack * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const G4Step *>
	Observer ()

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfTrack *>
	Observer ()

void	slotForUpdate (const EndOfTrack * iT)

virtual	~Observer ()

Private Member Functions
bool	stopAfter (const G4Step *)

void	update (const BeginOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const G4Step *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord >	cpvTokenDD4hep_

edm::ESGetToken< DDCompactView, IdealGeometryRecord >	cpvTokenDDD_

bool	fromdd4hep_

double	rMax_

std::unique_ptr< MaterialBudgetCastorHistos >	theHistoCastor_

std::unique_ptr< MaterialBudgetHcalHistos >	theHistoHcal_

double	zMax_

Additional Inherited Members
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 21 of file MaterialBudgetHcal.h.

Constructor & Destructor Documentation

◆ MaterialBudgetHcal() [1/2]

MaterialBudgetHcal::MaterialBudgetHcal ( const edm::ParameterSet & p )

Definition at line 16 of file MaterialBudgetHcal.cc.

References fromdd4hep_, edm::ParameterSet::getUntrackedParameter(), visualization-live-secondInstance_cfg::m, AlCaHLTBitMon_ParallelJobs::p, rMax_, theHistoCastor_, theHistoHcal_, and zMax_.

                                                                {
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("MaterialBudgetHcal");
   rMax_ = m_p.getUntrackedParameter<double>("RMax", 4.5) * CLHEP::m;
   zMax_ = m_p.getUntrackedParameter<double>("ZMax", 13.0) * CLHEP::m;
   fromdd4hep_ = m_p.getUntrackedParameter<bool>("Fromdd4hep", false);
   bool doHcal = m_p.getUntrackedParameter<bool>("DoHCAL", true);
   edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcal initialized with rMax " << rMax_ << " mm and zMax " << zMax_
                                      << " mm doHcal is set to " << doHcal << " and Fromdd4hep to " << fromdd4hep_;
   if (doHcal) {
     theHistoHcal_ = std::make_unique<MaterialBudgetHcalHistos>(m_p);
     theHistoCastor_.reset(nullptr);
   } else {
     theHistoHcal_.reset(nullptr);
     theHistoCastor_ = std::make_unique<MaterialBudgetCastorHistos>(m_p);
   }
 }

◆ MaterialBudgetHcal() [2/2]

MaterialBudgetHcal::MaterialBudgetHcal ( const MaterialBudgetHcal & )

delete

Member Function Documentation

◆ beginRun()

void MaterialBudgetHcal::beginRun ( edm::EventSetup const & es )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 41 of file MaterialBudgetHcal.cc.

References cpvTokenDD4hep_, cpvTokenDDD_, fromdd4hep_, edm::EventSetup::getHandle(), and theHistoHcal_.

                                                          {
   //----- Check that selected volumes are indeed part of the geometry
   // Numbering From DDD
   if (fromdd4hep_) {
     const auto& pDD = es.getHandle(cpvTokenDD4hep_);
     if (theHistoHcal_)
       theHistoHcal_->fillBeginJob((*pDD));
   } else {
     const auto& pDD = es.getHandle(cpvTokenDDD_);
     if (theHistoHcal_)
       theHistoHcal_->fillBeginJob((*pDD));
   }
 }

◆ operator=()

const MaterialBudgetHcal& MaterialBudgetHcal::operator= ( const MaterialBudgetHcal & )

delete

◆ registerConsumes()

void MaterialBudgetHcal::registerConsumes ( edm::ConsumesCollector cc )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 33 of file MaterialBudgetHcal.cc.

References edm::BeginRun, cpvTokenDD4hep_, cpvTokenDDD_, edm::ConsumesCollector::esConsumes(), and fromdd4hep_.

                                                                  {
   if (fromdd4hep_)
     cpvTokenDD4hep_ = cc.esConsumes<edm::Transition::BeginRun>();
   else
     cpvTokenDDD_ = cc.esConsumes<edm::Transition::BeginRun>();
   edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcal: Initialize the token for CompactView";
 }

◆ stopAfter()

bool MaterialBudgetHcal::stopAfter ( const G4Step * aStep )

private

Definition at line 84 of file MaterialBudgetHcal.cc.

References funct::abs(), rMax_, findQualityFiles::rr, zMax_, and geometryCSVtoXML::zz.

Referenced by update().

                                                       {
   G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
   double rr = hitPoint.perp();
   double zz = std::abs(hitPoint.z());
 
   if (rr > rMax_ || zz > zMax_) {
     edm::LogVerbatim("MaterialBudget") << " MaterialBudgetHcal::StopAfter R = " << rr << " and Z = " << zz;
     return true;
   } else {
     return false;
   }
 }

◆ update() [1/3]

void MaterialBudgetHcal::update ( const BeginOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 55 of file MaterialBudgetHcal.cc.

References theHistoCastor_, and theHistoHcal_.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                        {
   const G4Track* aTrack = (*trk)();  // recover G4 pointer if wanted
   if (theHistoHcal_)
     theHistoHcal_->fillStartTrack(aTrack);
   if (theHistoCastor_)
     theHistoCastor_->fillStartTrack(aTrack);
 }

◆ update() [2/3]

void MaterialBudgetHcal::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 63 of file MaterialBudgetHcal.cc.

References stopAfter(), theHistoCastor_, theHistoHcal_, and HLT_2022v15_cff::track.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                    {
   //---------- each step
   if (theHistoHcal_)
     theHistoHcal_->fillPerStep(aStep);
   if (theHistoCastor_)
     theHistoCastor_->fillPerStep(aStep);
 
   //----- Stop tracking after selected position
   if (stopAfter(aStep)) {
     G4Track* track = aStep->GetTrack();
     track->SetTrackStatus(fStopAndKill);
   }
 }

◆ update() [3/3]

void MaterialBudgetHcal::update ( const EndOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack *>.

Definition at line 77 of file MaterialBudgetHcal.cc.

References theHistoCastor_, and theHistoHcal_.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                      {
   if (theHistoHcal_)
     theHistoHcal_->fillEndTrack();
   if (theHistoCastor_)
     theHistoCastor_->fillEndTrack();
 }