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

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::pair< double, double > >	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 36 of file HGCPassive.h.

Member Typedef Documentation

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

private

Definition at line 60 of file HGCPassive.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file HGCPassive.cc.

References gather_cfg::cout, 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
   std::cout << "Name of the mother volume " << motherName_ << std::endl;
   unsigned k(0);
 #endif
   for (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 ( )

override

Definition at line 42 of file HGCPassive.cc.

42 {

43 }

HGCPassive::HGCPassive ( const HGCPassive & )

privatedelete

Member Function Documentation

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

private

Definition at line 150 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
       }
     }
   }
 }

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

private

Definition at line 176 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 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 45 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
	)

private

Definition at line 194 of file HGCPassive.cc.

References gather_cfg::cout, DEFINE_SIMWATCHER, plotBeamSpotDB::first, crabWrapper::key, edm::second(), store_, and ntuplemaker::time.

Referenced by update().

                                            {
 
   std::pair<G4LogicalVolume*,unsigned int> key(plv,copy);
   auto itr = store_.find(key);
   if (itr == store_.end()) {
     store_[key] = std::pair<double,double>(time,energy);
   } else {
     (itr->second).second += energy;
   }
 #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
 }

void HGCPassive::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 54 of file HGCPassive.cc.

References gather_cfg::cout, 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
     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 * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 79 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 * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 90 of file HGCPassive.cc.

References popcon2dropbox::copy(), gather_cfg::cout, findLV(), GeV, init_, createfilelist::int, mapLV_, storeInfo(), ntuplemaker::time, and topLV_.

                                             {
 
   if (aStep != nullptr) {
 
     G4VSensitiveDetector* curSD = aStep->GetPreStepPoint()->GetSensitiveDetector();
     if (curSD==nullptr) {
      
       G4TouchableHistory* touchable = (G4TouchableHistory*)aStep->GetPreStepPoint()->GetTouchable();
       G4LogicalVolume* plv = (G4LogicalVolume*)touchable->GetVolume()->GetLogicalVolume();
       auto it = (init_) ? mapLV_.find(plv) : findLV(plv);
       if (((aStep->GetPostStepPoint() == nullptr) || 
        (aStep->GetTrack()->GetNextVolume() == nullptr)) &&
       (aStep->IsLastStepInVolume())) {
 #ifdef EDM_ML_DEBUG
     std::cout << 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()) << std::endl;
 #endif
     double time   = aStep->GetTrack()->GetGlobalTime();
     double energy = (aStep->GetPreStepPoint()->GetKineticEnergy() +
              aStep->GetTotalEnergyDeposit())/CLHEP::GeV;
     if (it != mapLV_.end()) {
       storeInfo(it, plv, 0, time, energy);
     } else if (topLV_ != nullptr) {
       auto itr = (init_) ? mapLV_.find(topLV_) : findLV(topLV_);
       if (itr != mapLV_.end()) {
         storeInfo(itr, topLV_, 0, time, energy);
       }
     }
       } else if (it != mapLV_.end()) {
     unsigned int copy = (unsigned int)(touchable->GetReplicaNumber(0) + 
                        1000*touchable->GetReplicaNumber(1));
     double time = (aStep->GetPostStepPoint()->GetGlobalTime());
     double edeposit = (aStep->GetTotalEnergyDeposit())/CLHEP::GeV;
     storeInfo(it, plv, copy, time, edeposit);
       } else if (topLV_ != nullptr) {
     auto itr = findLV(topLV_);
     if (itr != mapLV_.end()) {
       double time = (aStep->GetPostStepPoint()->GetGlobalTime());
       double edeposit = (aStep->GetTotalEnergyDeposit())/CLHEP::GeV;
       storeInfo(itr, topLV_, 100, time, edeposit);
     }
       }
     }//if (curSD==NULL)
   }//if (aStep != NULL)
 
     
 }//end update aStep

Member Data Documentation

unsigned int HGCPassive::count_

private

Definition at line 75 of file HGCPassive.h.

Referenced by update().

bool HGCPassive::init_

private

Definition at line 76 of file HGCPassive.h.

Referenced by update().

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

private

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

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

std::string HGCPassive::motherName_

private

Definition at line 72 of file HGCPassive.h.

Referenced by findLV(), and HGCPassive().

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

private

Definition at line 77 of file HGCPassive.h.

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

G4LogicalVolume* HGCPassive::topLV_

private

Definition at line 70 of file HGCPassive.h.

Referenced by findLV(), and update().

G4VPhysicalVolume* HGCPassive::topPV_

private

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