#include <CaloSteppingAction.h>

Inheritance diagram for CaloSteppingAction:

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

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

	~CaloSteppingAction () override

Public Member Functions inherited from SimProducer
void	registerProducts (edm::ProducesCollector producesCollector)

	SimProducer ()

Public Member Functions inherited from SimWatcher
const SimWatcher &	operator= (const SimWatcher &)=delete

	SimWatcher ()

	SimWatcher (const SimWatcher &)=delete

virtual	~SimWatcher ()

Public Member Functions inherited from Observer< const BeginOfJob * >
	Observer ()

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

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 EndOfEvent * >
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const G4Step * >
	Observer ()

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

Private Types
typedef std::tuple< const G4LogicalVolume *, uint32_t, int, int, double, double, double, double, double, double, double >	PassiveData

Private Member Functions
double	curve_LY (double crystalLength, double crystalDepth) const

void	fillHit (uint32_t id, double dE, double time, int primID, uint16_t depth, double em, int flag)

void	fillHits (edm::PCaloHitContainer &cc, int type)

void	fillPassiveHits (edm::PassiveHitContainer &cc)

double	getBirkHC (double dE, double step, double chg, double dens) const

double	getBirkL3 (double dE, double step, double chg, double dens) const

uint16_t	getDepth (bool flag, double crystalDepth, double radl) const

uint32_t	getDetIDHC (int det, int lay, int depth, const math::XYZVectorD &pos) const

void	NaNTrap (const G4Step *) const

void	saveHits (int flag)

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

void	update (const BeginOfJob *job) 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 EndOfEvent *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
int	allSteps_

double	birkC1EC_

double	birkC1HC_

double	birkC2HC_

double	birkC3HC_

double	birkCutEC_

double	birkSlopeEC_

int	count_

std::unique_ptr< EcalBarrelNumberingScheme >	ebNumberingScheme_

std::unique_ptr< EcalEndcapNumberingScheme >	eeNumberingScheme_

int	eventID_

std::unique_ptr< HcalNumberingFromPS >	hcNumberingPS_

std::unique_ptr< HcalNumberingScheme >	hcNumberingScheme_

std::map< std::pair< int, CaloHitID >, CaloGVHit >	hitMap_ [nSD_]

std::map< const G4LogicalVolume *, std::string >	mapLV_

std::vector< std::string >	nameEBSD_

std::vector< std::string >	nameEESD_

std::vector< std::string >	nameHCSD_

std::vector< std::string >	nameHitC_

std::unique_ptr< CaloSlaveSD >	slave_ [nSD_]

double	slopeLY_

std::vector< PassiveData >	store_

double	timeSliceUnit_

std::vector< const G4LogicalVolume * >	volEBSD_

std::vector< const G4LogicalVolume * >	volEESD_

std::vector< const G4LogicalVolume * >	volHCSD_

std::map< const G4LogicalVolume *, double >	xtalMap_

Static Private Attributes
static const int	nSD_ = 3

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 52 of file CaloSteppingAction.h.

Member Typedef Documentation

◆ PassiveData

typedef std::tuple<const G4LogicalVolume *, uint32_t, int, int, double, double, double, double, double, double, double> CaloSteppingAction::PassiveData

private

Definition at line 104 of file CaloSteppingAction.h.

Constructor & Destructor Documentation

◆ CaloSteppingAction()

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

Definition at line 33 of file CaloSteppingAction.cc.

                                                                : count_(0) {
   edm::ParameterSet iC = p.getParameter<edm::ParameterSet>("CaloSteppingAction");
   nameEBSD_ = iC.getParameter<std::vector<std::string> >("EBSDNames");
   nameEESD_ = iC.getParameter<std::vector<std::string> >("EESDNames");
   nameHCSD_ = iC.getParameter<std::vector<std::string> >("HCSDNames");
   nameHitC_ = iC.getParameter<std::vector<std::string> >("HitCollNames");
   allSteps_ = iC.getParameter<int>("AllSteps");
   slopeLY_ = iC.getParameter<double>("SlopeLightYield");
   birkC1EC_ = iC.getParameter<double>("BirkC1EC");
   birkSlopeEC_ = iC.getParameter<double>("BirkSlopeEC");
   birkCutEC_ = iC.getParameter<double>("BirkCutEC");
   birkC1HC_ = iC.getParameter<double>("BirkC1HC");
   birkC2HC_ = iC.getParameter<double>("BirkC2HC");
   birkC3HC_ = iC.getParameter<double>("BirkC3HC");
   timeSliceUnit_ = iC.getUntrackedParameter<double>("TimeSliceUnit", 1.0);
  
   edm::LogVerbatim("Step") << "CaloSteppingAction:: " << nameEBSD_.size() << " names for EB SD's";
   for (unsigned int k = 0; k < nameEBSD_.size(); ++k)
     edm::LogVerbatim("Step") << "[" << k << "] " << nameEBSD_[k];
   edm::LogVerbatim("Step") << "CaloSteppingAction:: " << nameEESD_.size() << " names for EE SD's";
   for (unsigned int k = 0; k < nameEESD_.size(); ++k)
     edm::LogVerbatim("Step") << "[" << k << "] " << nameEESD_[k];
   edm::LogVerbatim("Step") << "CaloSteppingAction:: " << nameHCSD_.size() << " names for HC SD's";
   for (unsigned int k = 0; k < nameHCSD_.size(); ++k)
     edm::LogVerbatim("Step") << "[" << k << "] " << nameHCSD_[k];
   edm::LogVerbatim("Step") << "CaloSteppingAction::Constants for ECAL: slope " << slopeLY_ << " Birk constants "
                            << birkC1EC_ << ":" << birkSlopeEC_ << ":" << birkCutEC_;
   edm::LogVerbatim("Step") << "CaloSteppingAction::Constants for HCAL: Birk "
                            << "constants " << birkC1HC_ << ":" << birkC2HC_ << ":" << birkC3HC_;
   edm::LogVerbatim("Step") << "CaloSteppingAction::Constant for time slice " << timeSliceUnit_;
   edm::LogVerbatim("Step") << "CaloSteppingAction:: " << nameHitC_.size() << " hit collections";
   for (unsigned int k = 0; k < nameHitC_.size(); ++k)
     edm::LogVerbatim("Step") << "[" << k << "] " << nameHitC_[k];
  
   ebNumberingScheme_ = std::make_unique<EcalBarrelNumberingScheme>();
   eeNumberingScheme_ = std::make_unique<EcalEndcapNumberingScheme>();
   hcNumberingPS_ = std::make_unique<HcalNumberingFromPS>(iC);
   hcNumberingScheme_ = std::make_unique<HcalNumberingScheme>();
 #ifdef HcalNumberingTest
   hcNumbering_.reset(nullptr);
 #endif
   for (int k = 0; k < CaloSteppingAction::nSD_; ++k) {
     slave_[k] = std::make_unique<CaloSlaveSD>(nameHitC_[k]);
     produces<edm::PCaloHitContainer>(nameHitC_[k]);
   }
   if (allSteps_ != 0)
     produces<edm::PassiveHitContainer>("AllPassiveHits");
   edm::LogVerbatim("Step") << "CaloSteppingAction:: All Steps Flag " << allSteps_ << " for passive hits";
 }

References allSteps_, birkC1EC_, birkC1HC_, birkC2HC_, birkC3HC_, birkCutEC_, birkSlopeEC_, ebNumberingScheme_, eeNumberingScheme_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcNumberingPS_, hcNumberingScheme_, dqmdumpme::k, nameEBSD_, nameEESD_, nameHCSD_, nameHitC_, nSD_, AlCaHLTBitMon_ParallelJobs::p, slave_, slopeLY_, and timeSliceUnit_.

◆ ~CaloSteppingAction()

CaloSteppingAction::~CaloSteppingAction ( )

override

Definition at line 83 of file CaloSteppingAction.cc.

                                         {
   edm::LogVerbatim("Step") << "CaloSteppingAction: -------->  Total number of "
                            << "selected entries : " << count_;
 }

References count_.

Member Function Documentation

◆ curve_LY()

double CaloSteppingAction::curve_LY	(	double	crystalLength,
		double	crystalDepth
	)		const

private

Definition at line 387 of file CaloSteppingAction.cc.

                                                                                    {
   double weight = 1.;
   double dapd = crystalLength - crystalDepth;
   if (dapd >= -0.1 || dapd <= crystalLength + 0.1) {
     if (dapd <= 100.)
       weight = 1.0 + slopeLY_ - dapd * 0.01 * slopeLY_;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("Step") << "CaloSteppingAction::curve_LY " << crystalDepth << ":" << crystalLength << ":" << dapd
                              << ":" << weight;
 #endif
   } else {
     edm::LogWarning("Step") << "CaloSteppingAction: light coll curve : wrong "
                             << "distance to APD " << dapd << " crlength = " << crystalLength
                             << " crystal Depth = " << crystalDepth << " weight = " << weight;
   }
   return weight;
 }

References slopeLY_, and mps_merge::weight.

Referenced by update().

◆ fillHit()

void CaloSteppingAction::fillHit	(	uint32_t	id,
		double	dE,
		double	time,
		int	primID,
		uint16_t	depth,
		double	em,
		int	flag
	)

private

Definition at line 360 of file CaloSteppingAction.cc.

                                                                                                                      {
   CaloHitID currentID(id, time, primID, depth, timeSliceUnit_);
   double edepEM = (em ? dE : 0);
   double edepHAD = (em ? 0 : dE);
   std::pair<int, CaloHitID> evID = std::make_pair(eventID_, currentID);
   auto it = hitMap_[flag].find(evID);
   if (it != hitMap_[flag].end()) {
     (it->second).addEnergyDeposit(edepEM, edepHAD);
   } else {
     CaloGVHit aHit;
     aHit.setEventID(eventID_);
     aHit.setID(currentID);
     aHit.addEnergyDeposit(edepEM, edepHAD);
     hitMap_[flag][evID] = aHit;
   }
 }

References CaloGVHit::addEnergyDeposit(), LEDCalibrationChannels::depth, mps_fire::end, eventID_, RemoveAddSevLevel::flag, hitMap_, CaloGVHit::setEventID(), CaloGVHit::setID(), protons_cff::time, and timeSliceUnit_.

Referenced by update().

◆ fillHits()

void CaloSteppingAction::fillHits	(	edm::PCaloHitContainer &	cc,
		int	type
	)

private

Definition at line 102 of file CaloSteppingAction.cc.

                                                                     {
   edm::LogVerbatim("Step") << "CaloSteppingAction::fillHits for type " << type << " with "
                            << slave_[type].get()->hits().size() << " hits";
   cc = slave_[type].get()->hits();
   slave_[type].get()->Clean();
 }

References slave_.

Referenced by produce().

◆ fillPassiveHits()

void CaloSteppingAction::fillPassiveHits ( edm::PassiveHitContainer & cc )

private

Definition at line 109 of file CaloSteppingAction.cc.

                                                                    {
   edm::LogVerbatim("Step") << "CaloSteppingAction::fillPassiveHits with " << store_.size() << " hits";
   for (const auto& element : store_) {
     auto lv = std::get<0>(element);
     auto it = mapLV_.find(lv);
     if (it != mapLV_.end()) {
       PassiveHit hit(it->second,
                      std::get<1>(element),
                      std::get<5>(element),
                      std::get<6>(element),
                      std::get<4>(element),
                      std::get<2>(element),
                      std::get<3>(element),
                      std::get<7>(element),
                      std::get<8>(element),
                      std::get<9>(element),
                      std::get<10>(element));
       cc.emplace_back(hit);
     }
   }
 }

References mapLV_, and store_.

Referenced by produce().

◆ getBirkHC()

double CaloSteppingAction::getBirkHC	(	double	dE,
		double	step,
		double	chg,
		double	dens
	)		const

private

Definition at line 425 of file CaloSteppingAction.cc.

                                                                                                     {
   double weight = 1.;
   if (charge != 0. && step > 0.) {
     double dedx = dEStep / step;
     double rkb = birkC1HC_ / density;
     double c = birkC2HC_ * rkb * rkb;
     if (std::abs(charge) >= 2.)
       rkb /= birkC3HC_;
     weight = 1. / (1. + rkb * dedx + c * dedx * dedx);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("Step") << "CaloSteppingAction::getBirkHC Charge " << charge << " dE/dx " << dedx << " Birk Const "
                              << rkb << ", " << c << " Weight = " << weight << " dE " << dEStep;
 #endif
   }
   return weight;
 }

References funct::abs(), birkC1HC_, birkC2HC_, birkC3HC_, c, ALCARECOTkAlJpsiMuMu_cff::charge, fastSimProducer_cff::density, and mps_merge::weight.

Referenced by update().

◆ getBirkL3()

double CaloSteppingAction::getBirkL3	(	double	dE,
		double	step,
		double	chg,
		double	dens
	)		const

private

Definition at line 405 of file CaloSteppingAction.cc.

                                                                                                     {
   double weight = 1.;
   if (charge != 0. && step > 0.) {
     double dedx = dEStep / step;
     double rkb = birkC1EC_ / density;
     if (dedx > 0) {
       weight = 1. - birkSlopeEC_ * log(rkb * dedx);
       if (weight < birkCutEC_)
         weight = birkCutEC_;
       else if (weight > 1.)
         weight = 1.;
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("Step") << "CaloSteppingAction::getBirkL3 Charge " << charge << " dE/dx " << dedx << " Birk Const "
                              << rkb << " Weight = " << weight << " dE " << dEStep << " step " << step;
 #endif
   }
   return weight;
 }

References birkC1EC_, birkCutEC_, birkSlopeEC_, ALCARECOTkAlJpsiMuMu_cff::charge, fastSimProducer_cff::density, dqm-mbProfile::log, and mps_merge::weight.

Referenced by update().

◆ getDepth()

uint16_t CaloSteppingAction::getDepth	(	bool	flag,
		double	crystalDepth,
		double	radl
	)		const

private

Definition at line 377 of file CaloSteppingAction.cc.

                                                                                        {
   uint16_t depth1 = (flag ? 0 : PCaloHit::kEcalDepthRefz);
   uint16_t depth2 = (uint16_t)floor(crystalDepth / radl);
   uint16_t depth = (((depth2 & PCaloHit::kEcalDepthMask) << PCaloHit::kEcalDepthOffset) | depth1);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("Step") << "CaloSteppingAction::getDepth radl " << radl << ":" << crystalDepth << " depth " << depth;
 #endif
   return depth;
 }

References LEDCalibrationChannels::depth, RemoveAddSevLevel::flag, PCaloHit::kEcalDepthMask, PCaloHit::kEcalDepthOffset, and PCaloHit::kEcalDepthRefz.

Referenced by update().

◆ getDetIDHC()

uint32_t CaloSteppingAction::getDetIDHC	(	int	det,
		int	lay,
		int	depth,
		const math::XYZVectorD &	pos
	)		const

private

Definition at line 348 of file CaloSteppingAction.cc.

                                                                                                     {
   HcalNumberingFromDDD::HcalID tmp = hcNumberingPS_.get()->unitID(det, lay, depth, pos);
 #ifdef HcalNumberingTest
   auto id0 = hcNumberingScheme_.get()->getUnitID(tmp);
   HcalNumberingFromDDD::HcalID tmpO = hcNumbering_.get()->unitID(det, pos, depth, lay);
   auto idO = hcNumberingScheme_.get()->getUnitID(tmpO);
   std::string error = (id0 == idO) ? " ** OK **" : " ** ERROR **";
   edm::LogVerbatim("Step") << "Det ID " << HcalDetId(id0) << " Original " << HcalDetId(idO) << error;
 #endif
   return (hcNumberingScheme_.get()->getUnitID(tmp));
 }

References LEDCalibrationChannels::depth, relativeConstraints::error, hcNumberingPS_, hcNumberingScheme_, AlCaHLTBitMon_QueryRunRegistry::string, and createJobs::tmp.

Referenced by update().

◆ NaNTrap()

void CaloSteppingAction::NaNTrap ( const G4Step * aStep ) const

private

Definition at line 334 of file CaloSteppingAction.cc.

                                                           {
   auto currentPos = aStep->GetTrack()->GetPosition();
   double xyz = currentPos.x() + currentPos.y() + currentPos.z();
   auto currentMom = aStep->GetTrack()->GetMomentum();
   xyz += currentMom.x() + currentMom.y() + currentMom.z();
  
   if (edm::isNotFinite(xyz)) {
     auto pCurrentVol = aStep->GetPreStepPoint()->GetPhysicalVolume();
     auto& nameOfVol = pCurrentVol->GetName();
     throw cms::Exception("Unknown", "CaloSteppingAction")
         << " Corrupted Event - NaN detected in volume " << nameOfVol << "\n";
   }
 }

References Exception, and edm::isNotFinite().

Referenced by update().

◆ produce()

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

overridevirtual

Implements SimProducer.

Definition at line 88 of file CaloSteppingAction.cc.

                                                                   {
   for (int k = 0; k < CaloSteppingAction::nSD_; ++k) {
     saveHits(k);
     auto product = std::make_unique<edm::PCaloHitContainer>();
     fillHits(*product, k);
     e.put(std::move(product), nameHitC_[k]);
   }
   if (allSteps_ != 0) {
     std::unique_ptr<edm::PassiveHitContainer> hgcPH(new edm::PassiveHitContainer);
     fillPassiveHits(*hgcPH);
     e.put(std::move(hgcPH), "AllPassiveHits");
   }
 }

References allSteps_, MillePedeFileConverter_cfg::e, fillHits(), fillPassiveHits(), dqmdumpme::k, eostools::move(), nameHitC_, nSD_, and saveHits().

◆ saveHits()

void CaloSteppingAction::saveHits ( int flag )

private

Definition at line 442 of file CaloSteppingAction.cc.

                                           {
   edm::LogVerbatim("Step") << "CaloSteppingAction:: saveHits for type " << type << " with " << hitMap_[type].size()
                            << " hits";
   slave_[type].get()->ReserveMemory(hitMap_[type].size());
   for (auto const& hit : hitMap_[type]) {
     slave_[type].get()->processHits(hit.second.getUnitID(),
                                     0.001 * hit.second.getEM(),
                                     0.001 * hit.second.getHadr(),
                                     hit.second.getTimeSlice(),
                                     hit.second.getTrackID(),
                                     hit.second.getDepth());
   }
 }

References hitMap_, findQualityFiles::size, and slave_.

Referenced by produce().

◆ update() [1/5]

void CaloSteppingAction::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 220 of file CaloSteppingAction.cc.

                                                        {
   eventID_ = (*evt)()->GetEventID();
   edm::LogVerbatim("Step") << "CaloSteppingAction: Begin of event = " << eventID_;
   for (int k = 0; k < CaloSteppingAction::nSD_; ++k) {
     hitMap_[k].erase(hitMap_[k].begin(), hitMap_[k].end());
     slave_[k].get()->Initialize();
   }
   if (allSteps_ != 0)
     store_.clear();
 }

References allSteps_, mps_fire::end, eventID_, hitMap_, dqmdumpme::k, nSD_, slave_, and store_.

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

◆ update() [2/5]

void CaloSteppingAction::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 131 of file CaloSteppingAction.cc.

                                                      {
   edm::LogVerbatim("Step") << "CaloSteppingAction:: Enter BeginOfJob";
  
 #ifdef HcalNumberingTest
   // Numbering From DDD
   edm::ESHandle<HcalDDDSimConstants> hdc;
   (*job)()->get<HcalSimNumberingRecord>().get(hdc);
   const HcalDDDSimConstants* hcons_ = hdc.product();
   edm::LogVerbatim("Step") << "CaloSteppingAction:: Initialise "
                            << "HcalNumberingFromDDD";
   hcNumbering_ = std::make_unique<HcalNumberingFromDDD>(hcons_);
 #endif
 }

References edm::ESHandle< T >::product().

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

◆ update() [3/5]

void CaloSteppingAction::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 146 of file CaloSteppingAction.cc.

                                                      {
   int irun = (*run)()->GetRunID();
   edm::LogVerbatim("Step") << "CaloSteppingAction:: Begin of Run = " << irun;
  
   const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
   if (lvs) {
     std::map<const std::string, const G4LogicalVolume*> nameMap;
     std::map<const std::string, const G4LogicalVolume*>::const_iterator itr;
     for (auto lvi = lvs->begin(), lve = lvs->end(); lvi != lve; ++lvi) {
       nameMap.emplace((*lvi)->GetName(), *lvi);
       if (allSteps_ < 0)
         mapLV_[*lvi] = (*lvi)->GetName();
     }
  
     for (auto const& name : nameEBSD_) {
       for (itr = nameMap.begin(); itr != nameMap.end(); ++itr) {
         const std::string& lvname = itr->first;
         if (lvname.find(name) != std::string::npos) {
           volEBSD_.emplace_back(itr->second);
           int type = (lvname.find("refl") == std::string::npos) ? -1 : 1;
           G4Trap* solid = static_cast<G4Trap*>(itr->second->GetSolid());
           double dz = 2 * solid->GetZHalfLength() / CLHEP::mm;
           xtalMap_.insert(std::pair<const G4LogicalVolume*, double>(itr->second, dz * type));
           if ((allSteps_ > 0) && ((allSteps_ % 10) > 0))
             mapLV_[itr->second] = itr->first;
         }
       }
     }
     for (auto const& name : nameEESD_) {
       for (itr = nameMap.begin(); itr != nameMap.end(); ++itr) {
         const std::string& lvname = itr->first;
         if (lvname.find(name) != std::string::npos) {
           volEESD_.emplace_back(itr->second);
           int type = (lvname.find("refl") == std::string::npos) ? 1 : -1;
           G4Trap* solid = static_cast<G4Trap*>(itr->second->GetSolid());
           double dz = 2 * solid->GetZHalfLength() / CLHEP::mm;
           xtalMap_.insert(std::pair<const G4LogicalVolume*, double>(itr->second, dz * type));
           if ((allSteps_ > 0) && (((allSteps_ / 10) % 10) > 0))
             mapLV_[itr->second] = itr->first;
         }
       }
     }
  
     for (auto const& name : nameHCSD_) {
       for (itr = nameMap.begin(); itr != nameMap.end(); ++itr) {
         const std::string& lvname = itr->first;
         if (lvname.find(name) != std::string::npos) {
           volHCSD_.emplace_back(itr->second);
           if ((allSteps_ > 0) && (((allSteps_ / 100) % 10) > 0))
             mapLV_[itr->second] = itr->first;
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("Step") << volEBSD_.size() << " logical volumes for EB SD";
   for (unsigned int k = 0; k < volEBSD_.size(); ++k)
     edm::LogVerbatim("Step") << "[" << k << "] " << volEBSD_[k];
   edm::LogVerbatim("Step") << volEESD_.size() << " logical volumes for EE SD";
   for (unsigned int k = 0; k < volEESD_.size(); ++k)
     edm::LogVerbatim("Step") << "[" << k << "] " << volEESD_[k];
   edm::LogVerbatim("Step") << volHCSD_.size() << " logical volumes for HC SD";
   for (unsigned int k = 0; k < volHCSD_.size(); ++k)
     edm::LogVerbatim("Step") << "[" << k << "] " << volHCSD_[k];
   edm::LogVerbatim("Step") << mapLV_.size() << " logical volumes for Passive hits";
   unsigned int k(0);
   for (auto itr = mapLV_.begin(); itr != mapLV_.end(); ++itr) {
     edm::LogVerbatim("Step") << "[" << k << "] " << itr->second << ":" << itr->first;
     ++k;
   }
 #endif
 }

References allSteps_, PVValHelper::dz, dqmdumpme::k, mapLV_, Skims_PA_cff::name, nameEBSD_, nameEESD_, nameHCSD_, AlCaHLTBitMon_QueryRunRegistry::string, volEBSD_, volEESD_, volHCSD_, and xtalMap_.

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

◆ update() [4/5]

void CaloSteppingAction::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 328 of file CaloSteppingAction.cc.

                                                      {
   ++count_;
   // Fill event input
   edm::LogVerbatim("Step") << "CaloSteppingAction: EndOfEvent " << (*evt)()->GetEventID();
 }

References count_.

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

◆ update() [5/5]

void CaloSteppingAction::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 232 of file CaloSteppingAction.cc.

                                                    {
   //  edm::LogVerbatim("Step") <<"CaloSteppingAction: At each Step";
   NaNTrap(aStep);
   auto lv = aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
   bool hc = (std::find(volHCSD_.begin(), volHCSD_.end(), lv) != volHCSD_.end());
   bool eb = (std::find(volEBSD_.begin(), volEBSD_.end(), lv) != volEBSD_.end());
   bool ee = (std::find(volEESD_.begin(), volEESD_.end(), lv) != volEESD_.end());
   uint32_t unitID(0);
   if (hc || eb || ee) {
     double dEStep = aStep->GetTotalEnergyDeposit() / CLHEP::MeV;
     auto const theTrack = aStep->GetTrack();
     double time = theTrack->GetGlobalTime() / CLHEP::nanosecond;
     int primID = theTrack->GetTrackID();
     bool em = G4TrackToParticleID::isGammaElectronPositron(theTrack);
     auto const touch = aStep->GetPreStepPoint()->GetTouchable();
     auto const& hitPoint = aStep->GetPreStepPoint()->GetPosition();
     if (hc) {
       int depth = (touch->GetReplicaNumber(0)) % 10 + 1;
       int lay = (touch->GetReplicaNumber(0) / 10) % 100 + 1;
       int det = (touch->GetReplicaNumber(1)) / 1000;
       unitID = getDetIDHC(det, lay, depth, math::XYZVectorD(hitPoint.x(), hitPoint.y(), hitPoint.z()));
       if (unitID > 0 && dEStep > 0.0) {
         dEStep *= getBirkHC(dEStep,
                             (aStep->GetStepLength() / CLHEP::cm),
                             aStep->GetPreStepPoint()->GetCharge(),
                             (aStep->GetPreStepPoint()->GetMaterial()->GetDensity() / (CLHEP::g / CLHEP::cm3)));
         fillHit(unitID, dEStep, time, primID, 0, em, 2);
       }
     } else {
       EcalBaseNumber theBaseNumber;
       int size = touch->GetHistoryDepth() + 1;
       if (theBaseNumber.getCapacity() < size)
         theBaseNumber.setSize(size);
       //Get name and copy numbers
       if (size > 1) {
         for (int ii = 0; ii < size; ii++) {
           theBaseNumber.addLevel(touch->GetVolume(ii)->GetName(), touch->GetReplicaNumber(ii));
         }
       }
       unitID = (eb ? (ebNumberingScheme_->getUnitID(theBaseNumber)) : (eeNumberingScheme_->getUnitID(theBaseNumber)));
       if (unitID > 0 && dEStep > 0.0) {
         auto local = touch->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
         auto ite = xtalMap_.find(lv);
         double crystalLength = ((ite == xtalMap_.end()) ? 230.0 : std::abs(ite->second));
         double crystalDepth =
             ((ite == xtalMap_.end()) ? 0.0 : (std::abs(0.5 * (ite->second) + (local.z() / CLHEP::mm))));
         double radl = aStep->GetPreStepPoint()->GetMaterial()->GetRadlen() / CLHEP::mm;
         bool flag = ((ite == xtalMap_.end()) ? true : (((ite->second) >= 0) ? true : false));
         auto depth = getDepth(flag, crystalDepth, radl);
         dEStep *= (getBirkL3(dEStep,
                              (aStep->GetStepLength() / CLHEP::cm),
                              aStep->GetPreStepPoint()->GetCharge(),
                              (aStep->GetPreStepPoint()->GetMaterial()->GetDensity() / (CLHEP::g / CLHEP::cm3))) *
                    curve_LY(crystalLength, crystalDepth));
         fillHit(unitID, dEStep, time, primID, depth, em, (eb ? 0 : 1));
       }
     }
   }
  
   if (allSteps_ != 0) {
     auto it = mapLV_.find(lv);
     if (it != mapLV_.end()) {
       double energy = aStep->GetTotalEnergyDeposit() / CLHEP::MeV;
       auto const touch = aStep->GetPreStepPoint()->GetTouchable();
       double time = aStep->GetTrack()->GetGlobalTime() / CLHEP::nanosecond;
       int trackId = aStep->GetTrack()->GetTrackID();
       int pdg = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
       double stepl = (aStep->GetStepLength() / CLHEP::cm);
       double xp = aStep->GetPreStepPoint()->GetPosition().x() / CLHEP::cm;
       double yp = aStep->GetPreStepPoint()->GetPosition().y() / CLHEP::cm;
       double zp = aStep->GetPreStepPoint()->GetPosition().z() / CLHEP::cm;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("Step") << "CaloSteppingAction: Volume " << lv->GetName() << " History "
                                << touch->GetHistoryDepth() << " Pointers " << aStep->GetPostStepPoint() << ":"
                                << aStep->GetTrack()->GetNextVolume() << ":" << aStep->IsLastStepInVolume() << " E "
                                << energy << " T " << time << " PDG " << pdg << " step " << stepl << " Position (" << xp
                                << ", " << yp << ", " << zp << ")";
 #endif
       uint32_t copy = (allSteps_ < 0) ? 0 : unitID;
       if (((aStep->GetPostStepPoint() == nullptr) || (aStep->GetTrack()->GetNextVolume() == nullptr)) &&
           (aStep->IsLastStepInVolume())) {
         energy += (aStep->GetPreStepPoint()->GetKineticEnergy() / CLHEP::MeV);
       } else {
         time = aStep->GetPostStepPoint()->GetGlobalTime() / CLHEP::nanosecond;
         if (copy == 0)
           copy = (touch->GetHistoryDepth() < 1)
                      ? static_cast<uint32_t>(touch->GetReplicaNumber(0))
                      : static_cast<uint32_t>(touch->GetReplicaNumber(0) + 1000 * touch->GetReplicaNumber(1));
       }
       PassiveData key(std::make_tuple(lv, copy, trackId, pdg, time, energy, energy, stepl, xp, yp, zp));
       store_.push_back(key);
     }
   }
 }

References funct::abs(), EcalBaseNumber::addLevel(), allSteps_, filterCSVwithJSON::copy, curve_LY(), LEDCalibrationChannels::depth, ebNumberingScheme_, eeNumberingScheme_, HCALHighEnergyHPDFilter_cfi::energy, fillHit(), spr::find(), RemoveAddSevLevel::flag, g, getBirkHC(), getBirkL3(), EcalBaseNumber::getCapacity(), getDepth(), getDetIDHC(), cuy::ii, G4TrackToParticleID::isGammaElectronPositron(), crabWrapper::key, DTRecHitClients_cfi::local, mapLV_, MeV, NaNTrap(), EcalBaseNumber::setSize(), findQualityFiles::size, store_, protons_cff::time, funct::true, volEBSD_, volEESD_, volHCSD_, and xtalMap_.

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

Member Data Documentation

◆ allSteps_

int CaloSteppingAction::allSteps_

private

Definition at line 98 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), produce(), and update().

◆ birkC1EC_

double CaloSteppingAction::birkC1EC_

private

Definition at line 99 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getBirkL3().

◆ birkC1HC_

double CaloSteppingAction::birkC1HC_

private

Definition at line 100 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getBirkHC().

◆ birkC2HC_

double CaloSteppingAction::birkC2HC_

private

Definition at line 100 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getBirkHC().

◆ birkC3HC_

double CaloSteppingAction::birkC3HC_

private

Definition at line 101 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getBirkHC().

◆ birkCutEC_

double CaloSteppingAction::birkCutEC_

private

Definition at line 100 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getBirkL3().

◆ birkSlopeEC_

double CaloSteppingAction::birkSlopeEC_

private

Definition at line 99 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getBirkL3().

◆ count_

int CaloSteppingAction::count_

private

Definition at line 98 of file CaloSteppingAction.h.

Referenced by update(), and ~CaloSteppingAction().

◆ ebNumberingScheme_

std::unique_ptr<EcalBarrelNumberingScheme> CaloSteppingAction::ebNumberingScheme_

private

Definition at line 84 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and update().

◆ eeNumberingScheme_

std::unique_ptr<EcalEndcapNumberingScheme> CaloSteppingAction::eeNumberingScheme_

private

Definition at line 85 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and update().

◆ eventID_

int CaloSteppingAction::eventID_

private

Definition at line 98 of file CaloSteppingAction.h.

Referenced by fillHit(), and update().

◆ hcNumberingPS_

std::unique_ptr<HcalNumberingFromPS> CaloSteppingAction::hcNumberingPS_

private

Definition at line 86 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getDetIDHC().

◆ hcNumberingScheme_

std::unique_ptr<HcalNumberingScheme> CaloSteppingAction::hcNumberingScheme_

private

Definition at line 90 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and getDetIDHC().

◆ hitMap_

std::map<std::pair<int, CaloHitID>, CaloGVHit> CaloSteppingAction::hitMap_[nSD_]

private

Definition at line 102 of file CaloSteppingAction.h.

Referenced by fillHit(), saveHits(), and update().

◆ mapLV_

std::map<const G4LogicalVolume *, std::string> CaloSteppingAction::mapLV_

private

Definition at line 97 of file CaloSteppingAction.h.

Referenced by fillPassiveHits(), and update().

◆ nameEBSD_

std::vector<std::string> CaloSteppingAction::nameEBSD_

private

Definition at line 93 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and update().

◆ nameEESD_

std::vector<std::string> CaloSteppingAction::nameEESD_

private

Definition at line 93 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and update().

◆ nameHCSD_

std::vector<std::string> CaloSteppingAction::nameHCSD_

private

Definition at line 93 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and update().

◆ nameHitC_

std::vector<std::string> CaloSteppingAction::nameHitC_

private

Definition at line 94 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and produce().

◆ nSD_

const int CaloSteppingAction::nSD_ = 3

staticprivate

Definition at line 83 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), produce(), and update().

◆ slave_

std::unique_ptr<CaloSlaveSD> CaloSteppingAction::slave_[nSD_]

private

Definition at line 91 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), fillHits(), saveHits(), and update().

◆ slopeLY_

double CaloSteppingAction::slopeLY_

private

Definition at line 99 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and curve_LY().

◆ store_

std::vector<PassiveData> CaloSteppingAction::store_

private

Definition at line 105 of file CaloSteppingAction.h.

Referenced by fillPassiveHits(), and update().

◆ timeSliceUnit_

double CaloSteppingAction::timeSliceUnit_

private

Definition at line 101 of file CaloSteppingAction.h.

Referenced by CaloSteppingAction(), and fillHit().

◆ volEBSD_

std::vector<const G4LogicalVolume *> CaloSteppingAction::volEBSD_

private

Definition at line 95 of file CaloSteppingAction.h.

Referenced by update().

◆ volEESD_

std::vector<const G4LogicalVolume *> CaloSteppingAction::volEESD_

private

Definition at line 95 of file CaloSteppingAction.h.

Referenced by update().

◆ volHCSD_

std::vector<const G4LogicalVolume *> CaloSteppingAction::volHCSD_

private

Definition at line 95 of file CaloSteppingAction.h.

Referenced by update().

◆ xtalMap_

std::map<const G4LogicalVolume *, double> CaloSteppingAction::xtalMap_

private

Definition at line 96 of file CaloSteppingAction.h.

Referenced by update().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ PassiveData

Constructor & Destructor Documentation

◆ CaloSteppingAction()

◆ ~CaloSteppingAction()

Member Function Documentation

◆ curve_LY()

◆ fillHit()

◆ fillHits()

◆ fillPassiveHits()

◆ getBirkHC()

◆ getBirkL3()

◆ getDepth()

◆ getDetIDHC()

◆ NaNTrap()

◆ produce()

◆ saveHits()

◆ update() [1/5]

◆ update() [2/5]

◆ update() [3/5]

◆ update() [4/5]

◆ update() [5/5]

Member Data Documentation

◆ allSteps_

◆ birkC1EC_

◆ birkC1HC_

◆ birkC2HC_

◆ birkC3HC_

◆ birkCutEC_

◆ birkSlopeEC_

◆ count_

◆ ebNumberingScheme_

◆ eeNumberingScheme_

◆ eventID_

◆ hcNumberingPS_

◆ hcNumberingScheme_

◆ hitMap_

◆ mapLV_

◆ nameEBSD_

◆ nameEESD_

◆ nameHCSD_

◆ nameHitC_

◆ nSD_

◆ slave_

◆ slopeLY_

◆ store_

◆ timeSliceUnit_

◆ volEBSD_

◆ volEESD_

◆ volHCSD_

◆ xtalMap_