#include <HGCPassive.h>

Inheritance diagram for HGCPassive:

Public Member Functions
	HGCPassive (const edm::ParameterSet &p)

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

virtual	~HGCPassive ()

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 Member Functions
void	endOfEvent (edm::PassiveHitContainer &hgcPH, unsigned int k)

std::map< G4LogicalVolume *, std::pair< unsigned int, std::string > >::iterator	findLV (G4LogicalVolume *plv)

G4VPhysicalVolume *	getTopPV ()

	HGCPassive (const HGCPassive &)

const HGCPassive &	operator= (const HGCPassive &)

void	update (const BeginOfRun *run)
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfEvent *evt)
	This routine will be called when the appropriate signal arrives. More...

void	update (const G4Step *step)
	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::map< std::pair< G4LogicalVolume *, unsigned int >, std::pair< double, double > >	store_

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 34 of file HGCPassive.h.

Constructor & Destructor Documentation

HGCPassive::HGCPassive ( const edm::ParameterSet & p )

Definition at line 24 of file HGCPassive.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gen::k, LVNames_, and dataset::name.

                                                : count_(0), init_(false) {
 
   edm::ParameterSet m_Passive = p.getParameter<edm::ParameterSet>("HGCPassive");
   LVNames_  = m_Passive.getUntrackedParameter<std::vector<std::string> >("LVNames");
 
 #ifdef EDM_ML_DEBUG
   unsigned int k(0);
 #endif
   for (const auto& name : LVNames_) {
     produces<edm::PassiveHitContainer>(Form("%sPassiveHits",name.c_str()));
 #ifdef EDM_ML_DEBUG
     std::cout << "Collection name[" << k << "] " << name << std::endl;
     ++k;
 #endif
   }
 } 

HGCPassive::~HGCPassive ( )

virtual

Definition at line 41 of file HGCPassive.cc.

41 {

42 }

HGCPassive::HGCPassive ( const HGCPassive & )

private

Member Function Documentation

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

private

Definition at line 127 of file HGCPassive.cc.

References gather_cfg::cout, 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).second,(element.second).first);
     hgcPH.push_back(hit);
 #ifdef EDM_ML_DEBUG
     std::cout << "HGCPassive[" << k << "] Hit[" << kount << "] " << hit
           << std::endl;
     ++kount;
 #endif
       }
     }
   }
 }

std::map< G4LogicalVolume *, std::pair< unsigned int, std::string > >::iterator HGCPassive::findLV ( G4LogicalVolume * plv )

private

Definition at line 153 of file HGCPassive.cc.

References DEFINE_SIMWATCHER, gen::k, LVNames_, mapLV_, dataset::name, and AlCaHLTBitMon_QueryRunRegistry::string.

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;
       }
     }
   }
   return itr;
 }

G4VPhysicalVolume * HGCPassive::getTopPV ( )

private

Definition at line 149 of file HGCPassive.cc.

Referenced by update().

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

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

private

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

virtual

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::update ( const BeginOfRun * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 53 of file HGCPassive.cc.

References gather_cfg::cout, findLV(), plotBeamSpotDB::first, getTopPV(), init_, gen::k, mapLV_, and topPV_.

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

void HGCPassive::update ( const BeginOfEvent * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 78 of file HGCPassive.cc.

References count_, and store_.

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

void HGCPassive::update ( const G4Step * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 89 of file HGCPassive.cc.

References popcon2dropbox::copy(), gather_cfg::cout, findLV(), plotBeamSpotDB::first, init_, createfilelist::int, crabWrapper::key, mapLV_, NULL, edm::second(), store_, and ntuplemaker::time.

                                             {
 
   if (aStep != NULL) {
 
     G4VSensitiveDetector* curSD = aStep->GetPreStepPoint()->GetSensitiveDetector();
     if (curSD==NULL) {
       
       G4TouchableHistory* touchable = (G4TouchableHistory*)aStep->GetPreStepPoint()->GetTouchable();
       G4LogicalVolume* plv = (G4LogicalVolume*)touchable->GetVolume()->GetLogicalVolume();
       auto it = (init_) ? mapLV_.find(plv) : findLV(plv);
       if (it != mapLV_.end()) {
     unsigned int copy = (unsigned int)(touchable->GetReplicaNumber(0) + 
                        1000*touchable->GetReplicaNumber(1));
     std::pair<G4LogicalVolume*,unsigned int> key(plv,copy);
     auto itr = store_.find(key);
     double time = (aStep->GetPostStepPoint()->GetGlobalTime());
     if (itr == store_.end()) {
       store_[key] = std::pair<double,double>(time,0.0);
       itr         = store_.find(key);
     }
     double edeposit = aStep->GetTotalEnergyDeposit();
     (itr->second).second += edeposit;
 #ifdef EDM_ML_DEBUG
     std::cout << "HGCPassive: Element " << (it->second).first << ":" 
           << (it->second).second << ":" << copy << " T " 
           << (itr->second).first << " E " << (itr->second).second 
           << std::endl;
 #endif
       }//if( it != map.end() )
     }//if (curSD==NULL)
   }//if (aStep != NULL)
 
     
 }//end update aStep

Member Data Documentation

unsigned int HGCPassive::count_

private

Definition at line 67 of file HGCPassive.h.

Referenced by update().

bool HGCPassive::init_

private

Definition at line 68 of file HGCPassive.h.

Referenced by update().

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

private

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

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

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

private

Definition at line 69 of file HGCPassive.h.

Referenced by endOfEvent(), and update().

G4VPhysicalVolume* HGCPassive::topPV_

private

Definition at line 63 of file HGCPassive.h.

Referenced by update().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation