Inheritance diagram for HGCalTBMBProducer:

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

	HGCalTBMBProducer (const HGCalTBMBProducer &)=delete

const HGCalTBMBProducer &	operator= (const HGCalTBMBProducer &)=delete

void	produce (edm::Event &, const edm::EventSetup &) override

	~HGCalTBMBProducer () override=default

Public Member Functions inherited from SimProducer
const SimProducer &	operator= (const SimProducer &)=delete

void	registerProducts (edm::ProducesCollector producesCollector)

	SimProducer ()

	SimProducer (const SimProducer &)=delete

Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)

bool	isMT () const

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

virtual void	registerConsumes (edm::ConsumesCollector)

	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
int	findVolume (const G4VTouchable *touch, bool stop) const

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
std::vector< double >	intLen_

std::vector< std::string >	listNames_

MaterialAccountingCaloCollection	matcoll_

unsigned int	nList_

std::vector< double >	radLen_

std::vector< double >	stepLen_

std::string	stopName_

double	stopZ_

Additional Inherited Members
Protected Member Functions inherited from SimProducer
template<class T >
void	produces ()

template<class T >
void	produces (const std::string &instanceName)

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

Detailed Description

Definition at line 23 of file HGCalTBMBProducer.cc.

Constructor & Destructor Documentation

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

Definition at line 51 of file HGCalTBMBProducer.cc.

References edm::ParameterSet::getParameter(), isotrackApplyRegressor::k, listNames_, nList_, stopName_, stopZ_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                              {
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("HGCalTBMB");
   listNames_ = m_p.getParameter<std::vector<std::string> >("DetectorNames");
   stopName_ = m_p.getParameter<std::string>("StopName");
   stopZ_ = m_p.getParameter<double>("MaximumZ");
   nList_ = listNames_.size();
   edm::LogVerbatim("HGCSim") << "HGCalTBMBProducer initialized for " << nList_ << " volumes";
   for (unsigned int k = 0; k < nList_; ++k)
     edm::LogVerbatim("HGCSim") << " [" << k << "] " << listNames_[k];
   edm::LogVerbatim("HGCSim") << "Stop after " << stopZ_ << " or reaching volume " << stopName_;
 
   produces<MaterialAccountingCaloCollection>("HGCalTBMB");
 }

HGCalTBMBProducer::HGCalTBMBProducer ( const HGCalTBMBProducer & )

delete

HGCalTBMBProducer::~HGCalTBMBProducer ( )

overridedefault

Member Function Documentation

int HGCalTBMBProducer::findVolume	(	const G4VTouchable *	touch,
		bool	stop
	)		const

private

Definition at line 153 of file HGCalTBMBProducer.cc.

References cuy::ii, isotrackApplyRegressor::k, testEve_cfg::level, listNames_, mergeVDriftHistosByStation::name, nList_, stopName_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by stopAfter(), and update().

                                                                             {
   int ivol = -1;
   int level = (touch->GetHistoryDepth()) + 1;
   for (int ii = 0; ii < level; ii++) {
     std::string name = touch->GetVolume(ii)->GetName();
     if (stop) {
       if (strcmp(name.c_str(), stopName_.c_str()) == 0)
         ivol = 0;
     } else {
       for (unsigned int k = 0; k < nList_; ++k) {
         if (strcmp(name.c_str(), listNames_[k].c_str()) == 0) {
           ivol = k;
           break;
         }
       }
     }
     if (ivol >= 0)
       break;
   }
   return ivol;
 }

const HGCalTBMBProducer& HGCalTBMBProducer::operator= ( const HGCalTBMBProducer & )

delete

void HGCalTBMBProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &
	)

overridevirtual

Implements SimProducer.

Definition at line 65 of file HGCalTBMBProducer.cc.

References matcoll_, eostools::move(), and edm::Event::put().

                                                                  {
   std::unique_ptr<MaterialAccountingCaloCollection> hgc(new MaterialAccountingCaloCollection);
   for (auto const& mbc : matcoll_) {
     hgc->emplace_back(mbc);
   }
   e.put(std::move(hgc), "HGCalTBMB");
 }

bool HGCalTBMBProducer::stopAfter ( const G4Step * aStep )

private

Definition at line 136 of file HGCalTBMBProducer.cc.

References findVolume(), and stopZ_.

Referenced by update().

                                                      {
   bool flag(false);
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
   if (aStep->GetPostStepPoint() != nullptr)
     hitPoint = aStep->GetPostStepPoint()->GetPosition();
   double zz = hitPoint.z();
 
   if ((findVolume(touch, true) == 0) || (zz > stopZ_))
     flag = true;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << " HGCalTBMBProducer::Name " << touch->GetVolume(0)->GetName() << " z " << zz << " Flag"
                              << flag;
 #endif
   return flag;
 }

void HGCalTBMBProducer::update ( const BeginOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 73 of file HGCalTBMBProducer.cc.

References GeV, intLen_, nList_, radLen_, and stepLen_.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                       {
   radLen_ = std::vector<double>(nList_ + 1, 0);
   intLen_ = std::vector<double>(nList_ + 1, 0);
   stepLen_ = std::vector<double>(nList_ + 1, 0);
 
 #ifdef EDM_ML_DEBUG
   const G4Track* aTrack = (*trk)();  // recover G4 pointer if wanted
   const G4ThreeVector& mom = aTrack->GetMomentum();
   double theEnergy = aTrack->GetTotalEnergy();
   int theID = (int)(aTrack->GetDefinition()->GetPDGEncoding());
   edm::LogVerbatim("HGCSim") << "HGCalTBMBProducer: Track " << aTrack->GetTrackID() << " Code " << theID << " Energy "
                              << theEnergy / CLHEP::GeV << " GeV; Momentum " << mom;
 #endif
 }

void HGCalTBMBProducer::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 88 of file HGCalTBMBProducer.cc.

References findVolume(), intLen_, nList_, radLen_, stepLen_, stopAfter(), and HLT_FULL_cff::track.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                   {
   G4Material* material = aStep->GetPreStepPoint()->GetMaterial();
   double step = aStep->GetStepLength();
   double radl = material->GetRadlen();
   double intl = material->GetNuclearInterLength();
 
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   int indx = findVolume(touch, false);
 
   if (indx >= 0) {
     stepLen_[indx] += step;
     radLen_[indx] += (step / radl);
     intLen_[indx] += (step / intl);
   }
   stepLen_[nList_] += step;
   radLen_[nList_] += (step / radl);
   intLen_[nList_] += (step / intl);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "HGCalTBMBProducer::Step in " << touch->GetVolume(0)->GetLogicalVolume()->GetName()
                              << " Index " << indx << " Step " << step << " RadL " << step / radl << " IntL "
                              << step / intl;
 #endif
 
   if (stopAfter(aStep)) {
     G4Track* track = aStep->GetTrack();
     track->SetTrackStatus(fStopAndKill);
   }
 }

void HGCalTBMBProducer::update ( const EndOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 117 of file HGCalTBMBProducer.cc.

References cuy::ii, intLen_, listNames_, MaterialAccountingCalo::m_eta, MaterialAccountingCalo::m_intLen, MaterialAccountingCalo::m_phi, MaterialAccountingCalo::m_radLen, MaterialAccountingCalo::m_stepLen, matcoll_, mergeVDriftHistosByStation::name, nList_, radLen_, stepLen_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                     {
   MaterialAccountingCalo matCalo;
   matCalo.m_eta = 0;
   matCalo.m_phi = 0;
   for (unsigned int ii = 0; ii <= nList_; ++ii) {
     matCalo.m_stepLen.emplace_back(stepLen_[ii]);
     matCalo.m_radLen.emplace_back(radLen_[ii]);
     matCalo.m_intLen.emplace_back(intLen_[ii]);
 #ifdef EDM_ML_DEBUG
     std::string name("Total");
     if (ii < nList_)
       name = listNames_[ii];
     edm::LogVerbatim("HGCSim") << "HGCalTBMBProducer::Volume[" << ii << "]: " << name << " == Step " << stepLen_[ii]
                                << " RadL " << radLen_[ii] << " IntL " << intLen_[ii];
 #endif
   }
   matcoll_.emplace_back(matCalo);
 }

Member Data Documentation

std::vector<double> HGCalTBMBProducer::intLen_

private

Definition at line 48 of file HGCalTBMBProducer.cc.

Referenced by update().

std::vector<std::string> HGCalTBMBProducer::listNames_

private

Definition at line 43 of file HGCalTBMBProducer.cc.

Referenced by findVolume(), HGCalTBMBProducer(), and update().

MaterialAccountingCaloCollection HGCalTBMBProducer::matcoll_

private

Definition at line 47 of file HGCalTBMBProducer.cc.

Referenced by produce(), and update().

unsigned int HGCalTBMBProducer::nList_

private

Definition at line 46 of file HGCalTBMBProducer.cc.

Referenced by findVolume(), HGCalTBMBProducer(), and update().

std::vector<double> HGCalTBMBProducer::radLen_

private

Definition at line 48 of file HGCalTBMBProducer.cc.

Referenced by update().

std::vector<double> HGCalTBMBProducer::stepLen_

private

Definition at line 48 of file HGCalTBMBProducer.cc.

Referenced by update().

std::string HGCalTBMBProducer::stopName_

private

Definition at line 44 of file HGCalTBMBProducer.cc.

Referenced by findVolume(), and HGCalTBMBProducer().

double HGCalTBMBProducer::stopZ_

private

Definition at line 45 of file HGCalTBMBProducer.cc.

Referenced by HGCalTBMBProducer(), and stopAfter().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation