MaterialBudgetVolume Class Reference

Inheritance diagram for MaterialBudgetVolume:

Classes
struct	MatInfo

Public Member Functions
	MaterialBudgetVolume (const edm::ParameterSet &p)
	MaterialBudgetVolume (const MaterialBudgetVolume &)=delete
const MaterialBudgetVolume &	operator= (const MaterialBudgetVolume &)=delete
void	produce (edm::Event &, const edm::EventSetup &) override
	~MaterialBudgetVolume () override
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 BeginOfRun *>
	Observer ()
void	slotForUpdate (const BeginOfRun * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()
void	slotForUpdate (const BeginOfEvent * 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
void	endOfEvent (edm::MaterialInformationContainer &matbg)
int	findLV (const G4VTouchable *)
bool	loadLV ()
bool	stopAfter (const G4Step *)
void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfEvent *evt) 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 *step) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
int	iaddLevel_
bool	init_
std::vector< MatInfo >	lengths_
std::vector< int >	lvLevel_
std::vector< std::string >	lvNames_
std::map< int, std::pair< G4LogicalVolume *, int > >	mapLV_
double	rMax_
std::vector< MaterialInformation >	store_
double	zMax_

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 36 of file MaterialBudgetVolume.cc.

Constructor & Destructor Documentation

MaterialBudgetVolume() [1/2]

MaterialBudgetVolume::MaterialBudgetVolume

(

const edm::ParameterSet &

p

)

Definition at line 79 of file MaterialBudgetVolume.cc.

References edm::ParameterSet::getParameter(), iaddLevel_, isotrackApplyRegressor::k, lvLevel_, lvNames_, visualization-live-secondInstance_cfg::m, AlCaHLTBitMon_ParallelJobs::p, rMax_, and zMax_.

                                                                   : init_(false) {
  edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("MaterialBudgetVolume");

  lvNames_ = m_p.getParameter<std::vector<std::string> >("lvNames");
  lvLevel_ = m_p.getParameter<std::vector<int> >("lvLevels");
  iaddLevel_ = (m_p.getParameter<bool>("useDD4hep")) ? 1 : 0;
  rMax_ = m_p.getParameter<double>("rMax") * CLHEP::m;
  zMax_ = m_p.getParameter<double>("zMax") * CLHEP::m;

  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume: Studies Material budget for " << lvNames_.size()
                                     << " volumes with addLevel " << iaddLevel_ << " and with rMax " << rMax_
                                     << " mm; zMax " << zMax_ << " mm";
  std::ostringstream st1;
  for (unsigned int k = 0; k < lvNames_.size(); ++k)
    st1 << " [" << k << "] " << lvNames_[k] << " at " << lvLevel_[k];
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume: Volumes" << st1.str();

  produces<edm::MaterialInformationContainer>("MaterialInformation");
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume: will produce MaterialInformationContainer";
#endif
}

MaterialBudgetVolume() [2/2]

MaterialBudgetVolume::MaterialBudgetVolume ( const MaterialBudgetVolume &  )

delete

~MaterialBudgetVolume()

MaterialBudgetVolume::~MaterialBudgetVolume ( )

inlineoverride

Definition at line 46 of file MaterialBudgetVolume.cc.

{}

Member Function Documentation

endOfEvent()

void MaterialBudgetVolume::endOfEvent ( edm::MaterialInformationContainer &  matbg )

private

Definition at line 181 of file MaterialBudgetVolume.cc.

References info(), and store_.

Referenced by produce().

                                                                            {
#ifdef EDM_ML_DEBUG
  unsigned int kount(0);
#endif
  for (const auto& element : store_) {
    MaterialInformation info(element.vname(),
                             element.id(),
                             element.trackEta(),
                             element.trackPhi(),
                             element.stepLength(),
                             element.radiationLength(),
                             element.interactionLength());
    matbg.push_back(info);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume:: Info[" << kount << "] " << info;
    ++kount;
#endif
  }
}

findLV()

int MaterialBudgetVolume::findLV ( const G4VTouchable *  touch )

private

Definition at line 233 of file MaterialBudgetVolume.cc.

References mps_fire::i, cuy::ii, gammaJetAnalysis_CHSJECs_cff::level, GCP_tree_cfg::levels, mapLV_, and edm::second().

Referenced by update().

                                                          {
  int level(-1);
  int levels = ((touch->GetHistoryDepth()) + 1);
  for (const auto& lvs : mapLV_) {
    if ((lvs.second).second <= levels) {
      int ii = levels - (lvs.second).second;
      if ((touch->GetVolume(ii)->GetLogicalVolume()) == (lvs.second).first) {
        level = lvs.first;
        break;
      }
    }
  }
#ifdef EDM_ML_DEBUG
  if (level < 0) {
    edm::LogVerbatim("MaterialBudget") << "findLV: Gets " << level << " from " << levels << " levels in touchables";
    for (int i = 0; i < levels; ++i)
      edm::LogVerbatim("MaterialBudget") << "[" << (levels - i) << "] " << touch->GetVolume(i)->GetLogicalVolume()
                                         << " : " << touch->GetVolume(i)->GetLogicalVolume()->GetName();
  }
#endif
  return level;
}

loadLV()

bool MaterialBudgetVolume::loadLV ( )

private

Definition at line 201 of file MaterialBudgetVolume.cc.

References dqmdumpme::first, RemoveAddSevLevel::flag, mps_fire::i, iaddLevel_, isotrackApplyRegressor::k, lvLevel_, lvNames_, mapLV_, mergeVDriftHistosByStation::name, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by update().

                                  {
  const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
  bool flag(false);
  if (lvs != nullptr) {
    std::vector<G4LogicalVolume*>::const_iterator lvcite;
    for (unsigned int i = 0; i < lvNames_.size(); i++) {
      G4LogicalVolume* lv = nullptr;
      std::string name(lvNames_[i]);
      for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
        std::string namx(dd4hep::dd::noNamespace((*lvcite)->GetName()));
        if (namx == name) {
          lv = (*lvcite);
          break;
        }
      }
      if (lv != nullptr)
        mapLV_[i] = std::make_pair(lv, (lvLevel_[i] + iaddLevel_));
    }
    flag = true;
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume::Finds " << mapLV_.size() << " logical volumes";
    unsigned int k(0);
    for (const auto& lvs : mapLV_) {
      edm::LogVerbatim("MaterialBudget") << "Entry[" << k << "] " << lvs.first << ": (" << (lvs.second).first << ", "
                                         << (lvs.second).second << ") : " << lvNames_[lvs.first];
      ++k;
    }
#endif
  }
  return flag;
}

operator=()

const MaterialBudgetVolume& MaterialBudgetVolume::operator= ( const MaterialBudgetVolume &  )

delete

produce()

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

overridevirtual

Implements SimProducer.

Definition at line 102 of file MaterialBudgetVolume.cc.

References MillePedeFileConverter_cfg::e, endOfEvent(), and eostools::move().

                                                                    {
  std::unique_ptr<edm::MaterialInformationContainer> matbg(new edm::MaterialInformationContainer);
  endOfEvent(*matbg);
  e.put(std::move(matbg), "MaterialInformation");
}

stopAfter()

bool MaterialBudgetVolume::stopAfter ( const G4Step *  aStep )

private

Definition at line 256 of file MaterialBudgetVolume.cc.

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

Referenced by update().

                                                        {
  bool flag(false);
  if ((rMax_ > 0.0) && (zMax_ > 0.0)) {
    G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
    double rr = hitPoint.perp();
    double zz = std::abs(hitPoint.z());

    if (rr > rMax_ || zz > zMax_) {
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("MaterialBudget") << " MaterialBudgetHcal::StopAfter R = " << rr << " and Z = " << zz;
#endif
      flag = true;
    }
  }
  return flag;
}

update() [1/5]

void MaterialBudgetVolume::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 108 of file MaterialBudgetVolume.cc.

References init_, loadLV(), and mapLV_.

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

                                                       {
  init_ = loadLV();

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume::Finds " << mapLV_.size()
                                     << " logical volumes with return flag " << init_;
#endif
}

update() [2/5]

void MaterialBudgetVolume::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 117 of file MaterialBudgetVolume.cc.

References caHitNtupletGeneratorKernels::iev, init_, loadLV(), and store_.

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

                                                         {
#ifdef EDM_ML_DEBUG
  int iev = (*evt)()->GetEventID();
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume: =====> Begin event = " << iev << std::endl;
#endif
  if (!init_)
    init_ = loadLV();

  store_.clear();
}

update() [3/5]

void MaterialBudgetVolume::update ( const BeginOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 128 of file MaterialBudgetVolume.cc.

References lengths_, and mapLV_.

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

                                                         {
  lengths_ = std::vector<MatInfo>(mapLV_.size());

#ifdef EDM_ML_DEBUG
  const G4Track* aTrack = (*trk)();
  const G4ThreeVector& mom = aTrack->GetMomentum();
  double theEnergy = aTrack->GetTotalEnergy();
  int theID = static_cast<int>(aTrack->GetDefinition()->GetPDGEncoding());
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolumme: Track " << aTrack->GetTrackID() << " Code " << theID
                                     << " Energy " << theEnergy / CLHEP::GeV << " GeV; Momentum " << mom / CLHEP::GeV
                                     << " GeV";
#endif
}

update() [4/5]

void MaterialBudgetVolume::update ( const EndOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack *>.

Definition at line 167 of file MaterialBudgetVolume.cc.

References PVValHelper::eta, info(), isotrackApplyRegressor::k, lengths_, lvNames_, phi, and store_.

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

                                                       {
  const G4Track* aTrack = (*trk)();
  int id = aTrack->GetTrackID();
  double eta = aTrack->GetMomentumDirection().eta();
  double phi = aTrack->GetMomentumDirection().phi();
  for (unsigned int k = 0; k < lengths_.size(); ++k) {
    MaterialInformation info(lvNames_[k], id, eta, phi, lengths_[k].stepL, lengths_[k].radL, lengths_[k].intL);
    store_.emplace_back(info);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume::Volume[" << k << "]: " << info;
#endif
  }
}

update() [5/5]

void MaterialBudgetVolume::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 142 of file MaterialBudgetVolume.cc.

References findLV(), lengths_, stopAfter(), and HLT_2024v14_cff::track.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), 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 index = findLV(touch);
  if (index >= 0) {
    lengths_[index].stepL += step;
    lengths_[index].radL += (step / radl);
    lengths_[index].intL += (step / intl);
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetVolume::Step in "
                                     << touch->GetVolume(0)->GetLogicalVolume()->GetName() << " Index " << index
                                     << " Step " << step << " RadL " << step / radl << " IntL " << step / intl;
#endif

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

Member Data Documentation

iaddLevel_

int MaterialBudgetVolume::iaddLevel_

private

Definition at line 71 of file MaterialBudgetVolume.cc.

Referenced by loadLV(), and MaterialBudgetVolume().

init_

bool MaterialBudgetVolume::init_

private

Definition at line 73 of file MaterialBudgetVolume.cc.

Referenced by update().

lengths_

std::vector<MatInfo> MaterialBudgetVolume::lengths_

private

Definition at line 75 of file MaterialBudgetVolume.cc.

Referenced by update().

lvLevel_

std::vector<int> MaterialBudgetVolume::lvLevel_

private

Definition at line 70 of file MaterialBudgetVolume.cc.

Referenced by loadLV(), and MaterialBudgetVolume().

lvNames_

std::vector<std::string> MaterialBudgetVolume::lvNames_

private

Definition at line 69 of file MaterialBudgetVolume.cc.

Referenced by loadLV(), MaterialBudgetVolume(), and update().

mapLV_

std::map<int, std::pair<G4LogicalVolume*, int> > MaterialBudgetVolume::mapLV_

private

Definition at line 74 of file MaterialBudgetVolume.cc.

Referenced by findLV(), loadLV(), and update().

rMax_

double MaterialBudgetVolume::rMax_

private

Definition at line 72 of file MaterialBudgetVolume.cc.

Referenced by MaterialBudgetVolume(), and stopAfter().

store_

std::vector<MaterialInformation> MaterialBudgetVolume::store_

private

Definition at line 76 of file MaterialBudgetVolume.cc.

Referenced by endOfEvent(), and update().

zMax_

double MaterialBudgetVolume::zMax_

private

Definition at line 72 of file MaterialBudgetVolume.cc.

Referenced by MaterialBudgetVolume(), and stopAfter().