#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::ProducerBase &iProd)

	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 61 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(), gen::k, LVNames_, motherName_, dataset::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 42 of file HGCPassive.cc.

42 { }

HGCPassive::HGCPassive ( const HGCPassive & )

privatedelete

Member Function Documentation

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

private

Definition at line 167 of file HGCPassive.cc.

References plotBeamSpotDB::first, gen::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 194 of file HGCPassive.cc.

References gen::k, LVNames_, mapLV_, motherName_, dataset::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 190 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 44 of file HGCPassive.cc.

References endOfEvent(), gen::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 212 of file HGCPassive.cc.

References DEFINE_SIMWATCHER, plotBeamSpotDB::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 53 of file HGCPassive.cc.

References findLV(), plotBeamSpotDB::first, getTopPV(), init_, gen::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 80 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 91 of file HGCPassive.cc.

References popcon2dropbox::copy(), findLV(), GeV, mps_fire::i, cuy::ii, init_, createfilelist::int, hcalDigis_cfi::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 76 of file HGCPassive.h.

Referenced by update().

bool HGCPassive::init_

private

Definition at line 77 of file HGCPassive.h.

Referenced by update().

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

private

Definition at line 69 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 72 of file HGCPassive.h.

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

std::string HGCPassive::motherName_

private

Definition at line 73 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 78 of file HGCPassive.h.

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

G4LogicalVolume* HGCPassive::topLV_

private

Definition at line 71 of file HGCPassive.h.

Referenced by findLV(), and update().

G4VPhysicalVolume* HGCPassive::topPV_

private

Definition at line 70 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