HGCalTBPassive Class Reference

Inheritance diagram for HGCalTBPassive:

Public Member Functions
	HGCalTBPassive (const edm::ParameterSet &p)
	HGCalTBPassive (const HGCalTBPassive &)=delete
const HGCalTBPassive &	operator= (const HGCalTBPassive &)=delete
void	produce (edm::Event &, const edm::EventSetup &) override
	~HGCalTBPassive () 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 G4Step *>
	Observer ()
void	slotForUpdate (const G4Step * iT)
virtual	~Observer ()

Private Types
typedef std::map< G4LogicalVolume *, std::pair< unsigned int, std::string > >::iterator	volumeIterator

Private Member Functions
void	endOfEvent (edm::PassiveHitContainer &hgcPH, unsigned int k)
volumeIterator	findLV (G4LogicalVolume *plv)
G4VPhysicalVolume *	getTopPV ()
void	storeInfo (const volumeIterator itr, G4LogicalVolume *plv, unsigned int copy, double time, double energy, bool flag)
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 G4Step *step) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
const int	addlevel_
unsigned int	count_
bool	init_
const std::vector< std::string >	LVNames_
const edm::ParameterSet	m_Passive
std::map< G4LogicalVolume *, std::pair< unsigned int, std::string > >	mapLV_
const std::string	motherName_
std::map< std::pair< G4LogicalVolume *, unsigned int >, std::array< double, 3 > >	store_
G4LogicalVolume *	topLV_
G4VPhysicalVolume *	topPV_

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 HGCalTBPassive.cc.

Member Typedef Documentation

volumeIterator

typedef std::map<G4LogicalVolume *, std::pair<unsigned int, std::string> >::iterator HGCalTBPassive::volumeIterator

private

Definition at line 57 of file HGCalTBPassive.cc.

Constructor & Destructor Documentation

HGCalTBPassive() [1/2]

HGCalTBPassive::HGCalTBPassive

(

const edm::ParameterSet &

p

)

Definition at line 78 of file HGCalTBPassive.cc.

References addlevel_, isotrackApplyRegressor::k, LVNames_, motherName_, and mergeVDriftHistosByStation::name.

    : m_Passive(p.getParameter<edm::ParameterSet>("HGCalTBPassive")),
      LVNames_(m_Passive.getParameter<std::vector<std::string>>("LVNames")),
      motherName_(m_Passive.getParameter<std::string>("MotherName")),
      addlevel_((m_Passive.getParameter<bool>("IfDD4hep")) ? 1 : 0),
      topPV_(nullptr),
      topLV_(nullptr),
      count_(0),
      init_(false) {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCSim") << "Name of the mother volume " << motherName_ << " AddLevel " << addlevel_;
  unsigned k(0);
#endif
  for (const auto &name : LVNames_) {
    produces<edm::PassiveHitContainer>(Form("%sPassiveHits", name.c_str()));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCSim") << "Collection name[" << k << "] " << name;
    ++k;
#endif
  }
}

HGCalTBPassive() [2/2]

HGCalTBPassive::HGCalTBPassive ( const HGCalTBPassive &  )

delete

~HGCalTBPassive()

HGCalTBPassive::~HGCalTBPassive ( )

override

Definition at line 100 of file HGCalTBPassive.cc.

{}

Member Function Documentation

endOfEvent()

void HGCalTBPassive::endOfEvent ( edm::PassiveHitContainer &  hgcPH, unsigned int  k  )

private

Definition at line 214 of file HGCalTBPassive.cc.

References dqmdumpme::first, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, isotrackApplyRegressor::k, mapLV_, and store_.

Referenced by produce().

                                                                             {
#ifdef EDM_ML_DEBUG
  unsigned int kount(0);
#endif
  for (const auto &element : store_) {
    G4LogicalVolume *lv = (element.first).first;
    auto it = mapLV_.find(lv);
    if (it != mapLV_.end()) {
      if ((it->second).first == k) {
        PassiveHit hit(
            (it->second).second, (element.first).second, (element.second)[1], (element.second)[2], (element.second)[0]);
        hgcPH.push_back(hit);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCSim") << "HGCalTBPassive[" << k << "] Hit[" << kount << "] " << hit;
        ++kount;
#endif
      }
    }
  }
}

findLV()

HGCalTBPassive::volumeIterator HGCalTBPassive::findLV ( G4LogicalVolume *  plv )

private

Definition at line 239 of file HGCalTBPassive.cc.

References isotrackApplyRegressor::k, LVNames_, mapLV_, motherName_, mergeVDriftHistosByStation::name, DD4hep2DDDName::noNameSpace(), AlCaHLTBitMon_QueryRunRegistry::string, and topLV_.

Referenced by update().

                                                                        {
  auto itr = mapLV_.find(plv);
  if (itr == mapLV_.end()) {
    std::string name = DD4hep2DDDName::noNameSpace(static_cast<std::string>(plv->GetName()));
    for (unsigned int k = 0; k < LVNames_.size(); ++k) {
      if (name.find(LVNames_[k]) != std::string::npos) {
        mapLV_[plv] = std::pair<unsigned int, std::string>(k, name);
        itr = mapLV_.find(plv);
        break;
      }
    }
  }
  if (topLV_ == nullptr) {
    if (DD4hep2DDDName::noNameSpace(static_cast<std::string>(plv->GetName())) == motherName_)
      topLV_ = plv;
  }
  return itr;
}

getTopPV()

G4VPhysicalVolume * HGCalTBPassive::getTopPV ( )

private

Definition at line 235 of file HGCalTBPassive.cc.

Referenced by update().

                                            {
  return G4TransportationManager::GetTransportationManager()->GetNavigatorForTracking()->GetWorldVolume();
}

operator=()

const HGCalTBPassive& HGCalTBPassive::operator= ( const HGCalTBPassive &  )

delete

produce()

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

overridevirtual

Implements SimProducer.

Definition at line 102 of file HGCalTBPassive.cc.

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

                                                               {
  for (unsigned int k = 0; k < LVNames_.size(); ++k) {
    std::unique_ptr<edm::PassiveHitContainer> hgcPH(new edm::PassiveHitContainer);
    endOfEvent(*hgcPH, k);
    e.put(std::move(hgcPH), Form("%sPassiveHits", LVNames_[k].c_str()));
  }
}

storeInfo()

void HGCalTBPassive::storeInfo ( const volumeIterator  itr, G4LogicalVolume *  plv, unsigned int  copy, double  time, double  energy, bool  flag  )

private

Definition at line 258 of file HGCalTBPassive.cc.

References filterCSVwithJSON::copy, HBHEDarkening_cff::energy, dqmdumpme::first, RemoveAddSevLevel::flag, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, submitPVResolutionJobs::key, store_, and hcalRecHitTable_cff::time.

Referenced by update().

                                          {
  std::pair<G4LogicalVolume *, unsigned int> key(plv, copy);
  auto itr = store_.find(key);
  double ee = (flag) ? energy : 0;
  if (itr == store_.end()) {
    store_[key] = {{time, energy, energy}};
  } else {
    (itr->second)[1] += ee;
    (itr->second)[2] += energy;
  }
#ifdef EDM_ML_DEBUG
  itr = store_.find(key);
  edm::LogVerbatim("HGCSim") << "HGCalTBPassive: Element " << (it->second).first << ":" << (it->second).second << ":"
                             << copy << " T " << (itr->second)[0] << " E " << (itr->second)[1] << ":"
                             << (itr->second)[2];
#endif
}

update() [1/3]

void HGCalTBPassive::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 110 of file HGCalTBPassive.cc.

References findLV(), dqmdumpme::first, getTopPV(), init_, isotrackApplyRegressor::k, mapLV_, and topPV_.

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

                                                 {
  topPV_ = getTopPV();
  if (topPV_ == nullptr) {
    edm::LogWarning("HGCSim") << "Cannot find top level volume\n";
  } else {
    init_ = true;
    const G4LogicalVolumeStore *lvs = G4LogicalVolumeStore::GetInstance();
    for (auto lvcite : *lvs) {
      findLV(lvcite);
    }

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCSim") << "HGCalTBPassive::Finds " << mapLV_.size() << " logical volumes";
    unsigned int k(0);
    for (const auto &lvs : mapLV_) {
      edm::LogVerbatim("HGCSim") << "Entry[" << k << "] " << lvs.first << ": (" << (lvs.second).first << ", "
                                 << (lvs.second).second << ")";
      ++k;
    }
#endif
  }
}

update() [2/3]

void HGCalTBPassive::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 134 of file HGCalTBPassive.cc.

References count_, caHitNtupletGeneratorKernels::iev, and store_.

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

                                                   {
  int iev = (*evt)()->GetEventID();
  edm::LogVerbatim("HGCSim") << "HGCalTBPassive: =====> Begin event = " << iev << std::endl;

  ++count_;
  store_.clear();
}

update() [3/3]

void HGCalTBPassive::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 144 of file HGCalTBPassive.cc.

References filterCSVwithJSON::copy, HBHEDarkening_cff::energy, findLV(), mps_fire::i, init_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, gammaJetAnalysis_CHSJECs_cff::level, mapLV_, DD4hep2DDDName::noNameSpace(), storeInfo(), hcalRecHitTable_cff::time, and topLV_.

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

                                               {
  if (aStep != nullptr) {
    G4VSensitiveDetector *curSD = aStep->GetPreStepPoint()->GetSensitiveDetector();
    const G4VTouchable *touchable = aStep->GetPreStepPoint()->GetTouchable();

    int level = (touchable->GetHistoryDepth());
    if (curSD == nullptr) {
      G4LogicalVolume *plv = touchable->GetVolume()->GetLogicalVolume();
      auto it = (init_) ? mapLV_.find(plv) : findLV(plv);
      double time = aStep->GetTrack()->GetGlobalTime();
      double energy = (aStep->GetTotalEnergyDeposit()) / CLHEP::GeV;

      unsigned int copy(0);
      if (((aStep->GetPostStepPoint() == nullptr) || (aStep->GetTrack()->GetNextVolume() == nullptr)) &&
          (aStep->IsLastStepInVolume())) {
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCSim") << DD4hep2DDDName::noNameSpace(static_cast<std::string>(plv->GetName()))
                                   << " F|L Step " << aStep->IsFirstStepInVolume() << ":" << aStep->IsLastStepInVolume()
                                   << " Position" << aStep->GetPreStepPoint()->GetPosition() << " Track "
                                   << aStep->GetTrack()->GetDefinition()->GetParticleName() << " at"
                                   << aStep->GetTrack()->GetPosition() << " Volume " << aStep->GetTrack()->GetVolume()
                                   << ":" << aStep->GetTrack()->GetNextVolume() << " Status "
                                   << aStep->GetTrack()->GetTrackStatus() << " KE "
                                   << aStep->GetTrack()->GetKineticEnergy() << " Deposit "
                                   << aStep->GetTotalEnergyDeposit() << " Map " << (it != mapLV_.end());
#endif
        energy += (aStep->GetPreStepPoint()->GetKineticEnergy() / CLHEP::GeV);
      } else {
        time = (aStep->GetPostStepPoint()->GetGlobalTime());
        copy = (level < 2)
                   ? 0
                   : static_cast<unsigned int>(touchable->GetReplicaNumber(0) + 1000 * touchable->GetReplicaNumber(1));
      }
      if (it != mapLV_.end()) {
        storeInfo(it, plv, copy, time, energy, true);
      } else if (topLV_ != nullptr) {
        auto itr = findLV(topLV_);
        if (itr != mapLV_.end()) {
          storeInfo(itr, topLV_, copy, time, energy, true);
        }
      }
    }  // if (curSD==NULL)

    // Now for the mother volumes
    if (level > 0) {
      double energy = (aStep->GetTotalEnergyDeposit()) / CLHEP::GeV;
      double time = (aStep->GetTrack()->GetGlobalTime());

      for (int i = level; i > 0; --i) {
        G4LogicalVolume *plv = touchable->GetVolume(i)->GetLogicalVolume();
        auto it = (init_) ? mapLV_.find(plv) : findLV(plv);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCSim") << "Level: " << level << ":" << i << " "
                                   << DD4hep2DDDName::noNameSpace(static_cast<std::string>(plv->GetName()))
                                   << " flag in the List " << (it != mapLV_.end());
#endif
        if (it != mapLV_.end()) {
          unsigned int copy =
              (i == level) ? 0
                           : (unsigned int)(touchable->GetReplicaNumber(i) + 1000 * touchable->GetReplicaNumber(i + 1));
          storeInfo(it, plv, copy, time, energy, false);
        }
      }
    }
  }  // if (aStep != NULL)

}  // end update aStep

Member Data Documentation

addlevel_

const int HGCalTBPassive::addlevel_

private

Definition at line 67 of file HGCalTBPassive.cc.

Referenced by HGCalTBPassive().

count_

unsigned int HGCalTBPassive::count_

private

Definition at line 73 of file HGCalTBPassive.cc.

Referenced by update().

init_

bool HGCalTBPassive::init_

private

Definition at line 74 of file HGCalTBPassive.cc.

Referenced by update().

LVNames_

const std::vector<std::string> HGCalTBPassive::LVNames_

private

Definition at line 65 of file HGCalTBPassive.cc.

Referenced by findLV(), HGCalTBPassive(), and produce().

m_Passive

const edm::ParameterSet HGCalTBPassive::m_Passive

private

Definition at line 64 of file HGCalTBPassive.cc.

mapLV_

std::map<G4LogicalVolume *, std::pair<unsigned int, std::string> > HGCalTBPassive::mapLV_

private

Definition at line 70 of file HGCalTBPassive.cc.

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

motherName_

const std::string HGCalTBPassive::motherName_

private

Definition at line 66 of file HGCalTBPassive.cc.

Referenced by findLV(), and HGCalTBPassive().

store_

std::map<std::pair<G4LogicalVolume *, unsigned int>, std::array<double, 3> > HGCalTBPassive::store_

private

Definition at line 75 of file HGCalTBPassive.cc.

Referenced by endOfEvent(), storeInfo(), and update().

topLV_

G4LogicalVolume* HGCalTBPassive::topLV_

private

Definition at line 69 of file HGCalTBPassive.cc.

Referenced by findLV(), and update().

topPV_

G4VPhysicalVolume* HGCalTBPassive::topPV_

private

Definition at line 68 of file HGCalTBPassive.cc.

Referenced by update().