HGCPassive Class Reference

#include <HGCPassive.h>

Inheritance diagram for HGCPassive:

Public Member Functions
	HGCPassive (const edm::ParameterSet &p)
void	produce (edm::Event &, const edm::EventSetup &) override
	~HGCPassive () override
Public Member Functions inherited from SimProducer
void	registerProducts (edm::ProducesCollector producesCollector)
	SimProducer ()
Public Member Functions inherited from SimWatcher
	SimWatcher ()
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 ()
	HGCPassive (const HGCPassive &)=delete
const HGCPassive &	operator= (const HGCPassive &)=delete
void	storeInfo (const volumeIterator itr, G4LogicalVolume *plv, unsigned int copy, double time, double energy, bool flag)
void	update (const BeginOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfRun *run) 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
unsigned int	count_
bool	init_
std::vector< std::string >	LVNames_
std::map< G4LogicalVolume *, std::pair< unsigned int, std::string > >	mapLV_
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)

Detailed Description

Definition at line 37 of file HGCPassive.h.

Member Typedef Documentation

volumeIterator

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

private

Definition at line 59 of file HGCPassive.h.

Constructor & Destructor Documentation

HGCPassive() [1/2]

HGCPassive::HGCPassive

(

const edm::ParameterSet &

p

)

Definition at line 22 of file HGCPassive.cc.

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

References edm::ParameterSet::getParameter(), dqmdumpme::k, LVNames_, motherName_, Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, and AlCaHLTBitMon_QueryRunRegistry::string.

~HGCPassive()

HGCPassive::~HGCPassive ( )

override

Definition at line 40 of file HGCPassive.cc.

{}

HGCPassive() [2/2]

HGCPassive::HGCPassive ( const HGCPassive &  )

privatedelete

Member Function Documentation

endOfEvent()

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

private

Definition at line 151 of file HGCPassive.cc.

                                                                         {
#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") << "HGCPassive[" << k << "] Hit[" << kount << "] " << hit;
        ++kount;
#endif
      }
    }
  }
}

References dqmdumpme::first, dqmdumpme::k, mapLV_, and store_.

Referenced by produce().

findLV()

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

private

Definition at line 176 of file HGCPassive.cc.

                                                                {
  auto itr = mapLV_.find(plv);
  if (itr == mapLV_.end()) {
    std::string name = 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 (std::string(plv->GetName()) == motherName_)
      topLV_ = plv;
  }
  return itr;
}

References dqmdumpme::k, LVNames_, mapLV_, motherName_, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, and topLV_.

Referenced by update().

getTopPV()

G4VPhysicalVolume * HGCPassive::getTopPV ( )

private

Definition at line 172 of file HGCPassive.cc.

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

Referenced by update().

operator=()

const HGCPassive& HGCPassive::operator= ( const HGCPassive &  )

privatedelete

produce()

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

overridevirtual

Implements SimProducer.

Definition at line 42 of file HGCPassive.cc.

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

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

storeInfo()

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

private

Definition at line 195 of file HGCPassive.cc.

                                      {
  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") << "HGCPassive: Element " << (it->second).first << ":" << (it->second).second << ":"
                             << copy << " T " << (itr->second)[0] << " E " << (itr->second)[1] << ":"
                             << (itr->second)[2];
#endif
}

References filterCSVwithJSON::copy, HCALHighEnergyHPDFilter_cfi::energy, dqmdumpme::first, RemoveAddSevLevel::flag, crabWrapper::key, store_, and ntuplemaker::time.

Referenced by update().

update() [1/3]

void HGCPassive::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 74 of file HGCPassive.cc.

                                               {
  int iev = (*evt)()->GetEventID();
  edm::LogVerbatim("HGCSim") << "HGCPassive: =====> Begin event = " << iev << std::endl;
 
  ++count_;
  store_.clear();
}

References count_, and store_.

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

update() [2/3]

void HGCPassive::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 50 of file HGCPassive.cc.

                                             {
  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") << "HGCPassive::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
  }
}

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

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

update() [3/3]

void HGCPassive::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 84 of file HGCPassive.cc.

                                           {
  if (aStep != nullptr) {
    G4VSensitiveDetector* curSD = aStep->GetPreStepPoint()->GetSensitiveDetector();
    const G4VTouchable* touchable = aStep->GetPreStepPoint()->GetTouchable();
 
    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") << 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 = (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
    int level = (touchable->GetHistoryDepth());
    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 << " " << 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

References filterCSVwithJSON::copy, HCALHighEnergyHPDFilter_cfi::energy, findLV(), GeV, mps_fire::i, init_, createfilelist::int, personalPlayback::level, mapLV_, storeInfo(), ntuplemaker::time, and topLV_.

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

Member Data Documentation

count_

unsigned int HGCPassive::count_

private

Definition at line 73 of file HGCPassive.h.

Referenced by update().

init_

bool HGCPassive::init_

private

Definition at line 74 of file HGCPassive.h.

Referenced by update().

LVNames_

std::vector<std::string> HGCPassive::LVNames_

private

Definition at line 66 of file HGCPassive.h.

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

mapLV_

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

private

Definition at line 69 of file HGCPassive.h.

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

motherName_

std::string HGCPassive::motherName_

private

Definition at line 70 of file HGCPassive.h.

Referenced by findLV(), and HGCPassive().

store_

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

private

Definition at line 75 of file HGCPassive.h.

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

topLV_

G4LogicalVolume* HGCPassive::topLV_

private

Definition at line 68 of file HGCPassive.h.

Referenced by findLV(), and update().

topPV_

G4VPhysicalVolume* HGCPassive::topPV_

private

Definition at line 67 of file HGCPassive.h.

Referenced by update().