HGCalTBMBProducer Class Reference

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_
const std::vector< std::string >	listNames_
const edm::ParameterSet	m_p
MaterialAccountingCaloCollection	matcoll_
const unsigned int	nList_
std::vector< double >	radLen_
std::vector< double >	stepLen_
const std::string	stopName_
const 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() [1/2]

HGCalTBMBProducer::HGCalTBMBProducer

(

const edm::ParameterSet &

p

)

Definition at line 52 of file HGCalTBMBProducer.cc.

References dqmdumpme::k, listNames_, nList_, stopName_, and stopZ_.

    : 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() [2/2]

HGCalTBMBProducer::HGCalTBMBProducer ( const HGCalTBMBProducer &  )

delete

~HGCalTBMBProducer()

HGCalTBMBProducer::~HGCalTBMBProducer ( )

overridedefault

Member Function Documentation

findVolume()

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

private

Definition at line 154 of file HGCalTBMBProducer.cc.

References cuy::ii, dqmdumpme::k, personalPlayback::level, listNames_, Skims_PA_cff::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;
}

operator=()

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

delete

produce()

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

overridevirtual

Implements SimProducer.

Definition at line 66 of file HGCalTBMBProducer.cc.

References MillePedeFileConverter_cfg::e, caloTruthProducer_cfi::hgc, matcoll_, and eostools::move().

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

stopAfter()

bool HGCalTBMBProducer::stopAfter ( const G4Step *  aStep )

private

Definition at line 137 of file HGCalTBMBProducer.cc.

References findVolume(), RemoveAddSevLevel::flag, stopZ_, and geometryCSVtoXML::zz.

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

update() [1/3]

void HGCalTBMBProducer::update ( const BeginOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 74 of file HGCalTBMBProducer.cc.

References createfilelist::int, 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
}

update() [2/3]

void HGCalTBMBProducer::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 89 of file HGCalTBMBProducer.cc.

References findVolume(), intLen_, nList_, radLen_, stepLen_, stopAfter(), and HLT_2022v11_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);
  }
}

update() [3/3]

void HGCalTBMBProducer::update ( const EndOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack *>.

Definition at line 118 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_, Skims_PA_cff::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

intLen_

std::vector<double> HGCalTBMBProducer::intLen_

private

Definition at line 49 of file HGCalTBMBProducer.cc.

Referenced by update().

listNames_

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

private

Definition at line 44 of file HGCalTBMBProducer.cc.

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

m_p

const edm::ParameterSet HGCalTBMBProducer::m_p

private

Definition at line 43 of file HGCalTBMBProducer.cc.

matcoll_

MaterialAccountingCaloCollection HGCalTBMBProducer::matcoll_

private

Definition at line 48 of file HGCalTBMBProducer.cc.

Referenced by produce(), and update().

nList_

const unsigned int HGCalTBMBProducer::nList_

private

Definition at line 47 of file HGCalTBMBProducer.cc.

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

radLen_

std::vector<double> HGCalTBMBProducer::radLen_

private

Definition at line 49 of file HGCalTBMBProducer.cc.

Referenced by update().

stepLen_

std::vector<double> HGCalTBMBProducer::stepLen_

private

Definition at line 49 of file HGCalTBMBProducer.cc.

Referenced by update().

stopName_

const std::string HGCalTBMBProducer::stopName_

private

Definition at line 45 of file HGCalTBMBProducer.cc.

Referenced by findVolume(), and HGCalTBMBProducer().

stopZ_

const double HGCalTBMBProducer::stopZ_

private

Definition at line 46 of file HGCalTBMBProducer.cc.

Referenced by HGCalTBMBProducer(), and stopAfter().