#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 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
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

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::HGCPassive ( const edm::ParameterSet & p )

Definition at line 22 of file HGCPassive.cc.

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

                                                : 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("ValidHGCal") << "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("ValidHGCal") << "Collection name[" << k << "] " << name;
     ++k;
 #endif
   }
 }

HGCPassive::~HGCPassive ( )

override

Definition at line 40 of file HGCPassive.cc.

40 {}

HGCPassive::HGCPassive ( const HGCPassive & )

privatedelete

Member Function Documentation

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

private

Definition at line 151 of file HGCPassive.cc.

References dqmdumpme::first, dqmdumpme::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("ValidHGCal") << "HGCPassive[" << k << "] Hit[" << kount << "] " << hit;
         ++kount;
 #endif
       }
     }
   }
 }

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

private

Definition at line 176 of file HGCPassive.cc.

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

Referenced by update().

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

G4VPhysicalVolume * HGCPassive::getTopPV ( )

private

Definition at line 172 of file HGCPassive.cc.

Referenced by update().

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

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

privatedelete

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

overridevirtual

Implements SimProducer.

Definition at line 42 of file HGCPassive.cc.

References endOfEvent(), dqmdumpme::k, LVNames_, eostools::move(), and edm::Event::put().

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

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.

References DEFINE_SIMWATCHER, HCALHighEnergyHPDFilter_cfi::energy, dqmdumpme::first, RemoveAddSevLevel::flag, crabWrapper::key, store_, and ntuplemaker::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("ValidHGCal") << "HGCPassive: Element " << (it->second).first << ":" << (it->second).second << ":"
                                  << copy << " T " << (itr->second)[0] << " E " << (itr->second)[1] << ":"
                                  << (itr->second)[2];
 #endif
 }

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.

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

                                              {
   topPV_ = getTopPV();
   if (topPV_ == nullptr) {
     edm::LogWarning("HGCPassive") << "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("ValidHGCal") << "HGCPassive::Finds " << mapLV_.size() << " logical volumes";
     unsigned int k(0);
     for (const auto& lvs : mapLV_) {
       edm::LogVerbatim("ValidHGCal") << "Entry[" << k << "] " << lvs.first << ": (" << (lvs.second).first << ", "
                                      << (lvs.second).second << ")";
       ++k;
     }
 #endif
   }
 }

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.

References count_, and store_.

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

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.

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

                                            {
   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("ValidHGCal") << 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("ValidHGCal") << "Level: " << ii << ":" << 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

Member Data Documentation

unsigned int HGCPassive::count_

private

Definition at line 73 of file HGCPassive.h.

Referenced by update().

bool HGCPassive::init_

private

Definition at line 74 of file HGCPassive.h.

Referenced by update().

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

private

Definition at line 66 of file HGCPassive.h.

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

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().

std::string HGCPassive::motherName_

private

Definition at line 70 of file HGCPassive.h.

Referenced by findLV(), and HGCPassive().

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().

G4LogicalVolume* HGCPassive::topLV_

private

Definition at line 68 of file HGCPassive.h.

Referenced by findLV(), and update().

G4VPhysicalVolume* HGCPassive::topPV_

private

Definition at line 67 of file HGCPassive.h.

Referenced by update().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation