#include <MaterialBudgetHcal.h>

Inheritance diagram for MaterialBudgetHcal:

Public Member Functions
	MaterialBudgetHcal (const edm::ParameterSet &)

	MaterialBudgetHcal (const MaterialBudgetHcal &)=delete

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

Public Member Functions inherited from SimWatcher
	SimWatcher ()

virtual	~SimWatcher ()

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

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

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 BeginOfJob *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfTrack *) 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...

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

Private Attributes
bool	fromdd4hep_

double	rMax_

std::unique_ptr< MaterialBudgetCastorHistos >	theHistoCastor_

std::unique_ptr< MaterialBudgetHcalHistos >	theHistoHcal_

double	zMax_

Additional Inherited Members

Detailed Description

Definition at line 18 of file MaterialBudgetHcal.h.

Constructor & Destructor Documentation

◆ MaterialBudgetHcal() [1/2]

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

Definition at line 22 of file MaterialBudgetHcal.cc.

                                                                {
   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);
   }
 }

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

◆ MaterialBudgetHcal() [2/2]

MaterialBudgetHcal::MaterialBudgetHcal ( const MaterialBudgetHcal & )

delete

Member Function Documentation

◆ operator=()

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

delete

◆ stopAfter()

bool MaterialBudgetHcal::stopAfter ( const G4Step * aStep )

private

Definition at line 84 of file MaterialBudgetHcal.cc.

                                                       {
   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;
   }
 }

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

Referenced by update().

◆ update() [1/4]

void MaterialBudgetHcal::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 39 of file MaterialBudgetHcal.cc.

                                                      {
   //----- Check that selected volumes are indeed part of the geometry
   // Numbering From DDD
   if (fromdd4hep_) {
     edm::ESTransientHandle<cms::DDCompactView> pDD;
     (*job)()->get<IdealGeometryRecord>().get(pDD);
     if (theHistoHcal_)
       theHistoHcal_->fillBeginJob((*pDD));
   } else {
     edm::ESTransientHandle<DDCompactView> pDD;
     (*job)()->get<IdealGeometryRecord>().get(pDD);
     if (theHistoHcal_)
       theHistoHcal_->fillBeginJob((*pDD));
   }
 }

References fromdd4hep_, and theHistoHcal_.

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

◆ update() [2/4]

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.

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

References theHistoCastor_, and theHistoHcal_.

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

◆ update() [3/4]

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.

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

References theHistoCastor_, and theHistoHcal_.

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

◆ update() [4/4]

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.

                                                    {
   //---------- 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);
   }
 }