Inheritance diagram for CaloSteppingAction:

Public Member Functions
void	beginRun (edm::EventSetup const &) override

	CaloSteppingAction (const edm::ParameterSet &p)

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

void	registerConsumes (edm::ConsumesCollector) override

	~CaloSteppingAction () override

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

void	registerProducts (edm::ProducesCollector producesCollector)

	SimProducer ()

	SimProducer (const SimProducer &)=delete

Public Member Functions inherited from SimWatcher
bool	isMT () const

const SimWatcher &	operator= (const SimWatcher &)=delete

	SimWatcher ()

	SimWatcher (const SimWatcher &)=delete

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

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

Private Attributes
const int	allSteps_

const double	birkC1EC_

const double	birkC1HC_

const double	birkC2HC_

const double	birkC3HC_

const double	birkCutEC_

const 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_]

const edm::ParameterSet	iC_

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

const std::vector< std::string >	nameEBSD_

const std::vector< std::string >	nameEESD_

const std::vector< std::string >	nameHCSD_

const std::vector< std::string >	nameHitC_

std::unique_ptr< CaloSlaveSD >	slave_ [nSD_]

const double	slopeLY_

std::vector< PassiveData >	store_

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

Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 67 of file CaloSteppingAction.cc.

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 122 of file CaloSteppingAction.cc.

Constructor & Destructor Documentation

◆ CaloSteppingAction()

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

Definition at line 126 of file CaloSteppingAction.cc.

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

     : 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)),
       count_(0) {
   edm::ParameterSet iC = p.getParameter<edm::ParameterSet>("CaloSteppingAction");
 
   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";
 }

◆ ~CaloSteppingAction()

CaloSteppingAction::~CaloSteppingAction ( )

override

Definition at line 178 of file CaloSteppingAction.cc.

References count_.

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

Member Function Documentation

◆ beginRun()

void CaloSteppingAction::beginRun ( edm::EventSetup const & es )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 233 of file CaloSteppingAction.cc.

References edm::EventSetup::getData().

                                                          {
   edm::LogVerbatim("Step") << "CaloSteppingAction:: Enter BeginOfRun";
 
 #ifdef HcalNumberingTest
   // Numbering From DDD
   const HcalDDDSimConstants* hcons_ = &es.getData(ddconsToken_);
   edm::LogVerbatim("Step") << "CaloSteppingAction:: Initialise "
                            << "HcalNumberingFromDDD";
   hcNumbering_ = std::make_unique<HcalNumberingFromDDD>(hcons_);
 #endif
 }

◆ curve_LY()

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

private

Definition at line 487 of file CaloSteppingAction.cc.

References slopeLY_, and mps_merge::weight.

Referenced by update().

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

◆ fillHit()

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

private

Definition at line 460 of file CaloSteppingAction.cc.

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

Referenced by update().

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

◆ fillHits()

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

private

Definition at line 204 of file CaloSteppingAction.cc.

References gpuPixelDoublets::cc, and slave_.

Referenced by produce().

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

◆ fillPassiveHits()

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

private

Definition at line 211 of file CaloSteppingAction.cc.

References gpuPixelDoublets::cc, mapLV_, and store_.

Referenced by produce().

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

◆ getBirkHC()

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

private

Definition at line 525 of file CaloSteppingAction.cc.

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

Referenced by update().

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

◆ getBirkL3()

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

private

Definition at line 505 of file CaloSteppingAction.cc.

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

Referenced by update().

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

◆ getDepth()

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

private

Definition at line 477 of file CaloSteppingAction.cc.

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

Referenced by update().

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

◆ getDetIDHC()

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

private

Definition at line 448 of file CaloSteppingAction.cc.

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

Referenced by update().

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

◆ NaNTrap()

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

private

Definition at line 434 of file CaloSteppingAction.cc.

References Exception, and edm::isNotFinite().

Referenced by update().

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

◆ produce()

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

overridevirtual

Implements SimProducer.

Definition at line 190 of file CaloSteppingAction.cc.

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

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

◆ registerConsumes()

void CaloSteppingAction::registerConsumes ( edm::ConsumesCollector cc )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 183 of file CaloSteppingAction.cc.

References edm::BeginRun, and gpuPixelDoublets::cc.

                                                                  {
 #ifdef HcalNumberingTest
   ddconsToken_ = cc.esConsumes<HcalDDDSimConstants, HcalSimNumberingRecord, edm::Transition::BeginRun>();
   edm::LogVerbatim("Step") << "CaloSteppingAction::Initialize ESGetToken for HcalDDDSimConstants";
 #endif
 }

◆ saveHits()

void CaloSteppingAction::saveHits ( int flag )

private

Definition at line 542 of file CaloSteppingAction.cc.

References hitMap_, findQualityFiles::size, and slave_.

Referenced by produce().

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

◆ update() [1/4]

void CaloSteppingAction::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 246 of file CaloSteppingAction.cc.

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

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

◆ update() [2/4]

void CaloSteppingAction::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 320 of file CaloSteppingAction.cc.

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

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

◆ update() [3/4]

void CaloSteppingAction::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 332 of file CaloSteppingAction.cc.

References funct::abs(), EcalBaseNumber::addLevel(), allSteps_, filterCSVwithJSON::copy, curve_LY(), hcalRecHitTable_cff::depth, ebNumberingScheme_, eeNumberingScheme_, hcalRecHitTable_cff::energy, fillHit(), spr::find(), RemoveAddSevLevel::flag, g, getBirkHC(), getBirkL3(), EcalBaseNumber::getCapacity(), getDepth(), getDetIDHC(), cuy::ii, G4TrackToParticleID::isGammaElectronPositron(), submitPVResolutionJobs::key, DTRecHitClients_cfi::local, mapLV_, NaNTrap(), EcalBaseNumber::setSize(), findQualityFiles::size, store_, hcalRecHitTable_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().

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

◆ update() [4/4]

void CaloSteppingAction::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 428 of file CaloSteppingAction.cc.

References count_.

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

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

Member Data Documentation

◆ allSteps_

const int CaloSteppingAction::allSteps_

private

Definition at line 115 of file CaloSteppingAction.cc.

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

◆ birkC1EC_

const double CaloSteppingAction::birkC1EC_

private

Definition at line 116 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkL3().

◆ birkC1HC_

const double CaloSteppingAction::birkC1HC_

private

Definition at line 117 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkHC().

◆ birkC2HC_

const double CaloSteppingAction::birkC2HC_

private

Definition at line 117 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkHC().

◆ birkC3HC_

const double CaloSteppingAction::birkC3HC_

private

Definition at line 118 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkHC().

◆ birkCutEC_

const double CaloSteppingAction::birkCutEC_

private

Definition at line 117 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkL3().

◆ birkSlopeEC_

const double CaloSteppingAction::birkSlopeEC_

private

Definition at line 116 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkL3().

◆ count_

int CaloSteppingAction::count_

private

Definition at line 119 of file CaloSteppingAction.cc.

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

◆ ebNumberingScheme_

std::unique_ptr<EcalBarrelNumberingScheme> CaloSteppingAction::ebNumberingScheme_

private

Definition at line 99 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

◆ eeNumberingScheme_

std::unique_ptr<EcalEndcapNumberingScheme> CaloSteppingAction::eeNumberingScheme_

private

Definition at line 100 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

◆ eventID_

int CaloSteppingAction::eventID_

private

Definition at line 119 of file CaloSteppingAction.cc.

Referenced by fillHit(), and update().

◆ hcNumberingPS_

std::unique_ptr<HcalNumberingFromPS> CaloSteppingAction::hcNumberingPS_

private

Definition at line 101 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getDetIDHC().

◆ hcNumberingScheme_

std::unique_ptr<HcalNumberingScheme> CaloSteppingAction::hcNumberingScheme_

private

Definition at line 106 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getDetIDHC().

◆ hitMap_

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

private

Definition at line 120 of file CaloSteppingAction.cc.

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

◆ iC_

const edm::ParameterSet CaloSteppingAction::iC_

private

Definition at line 109 of file CaloSteppingAction.cc.

◆ mapLV_

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

private

Definition at line 114 of file CaloSteppingAction.cc.

Referenced by fillPassiveHits(), and update().

◆ nameEBSD_

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

private

Definition at line 110 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

◆ nameEESD_

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

private

Definition at line 110 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

◆ nameHCSD_

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

private

Definition at line 110 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

◆ nameHitC_

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

private

Definition at line 111 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and produce().

◆ nSD_

const int CaloSteppingAction::nSD_ = 3

staticprivate

Definition at line 98 of file CaloSteppingAction.cc.

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

◆ slave_

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

private

Definition at line 107 of file CaloSteppingAction.cc.

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

◆ slopeLY_

const double CaloSteppingAction::slopeLY_

private

Definition at line 116 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and curve_LY().

◆ store_

std::vector<PassiveData> CaloSteppingAction::store_

private

Definition at line 123 of file CaloSteppingAction.cc.

Referenced by fillPassiveHits(), and update().

◆ timeSliceUnit_

const double CaloSteppingAction::timeSliceUnit_

private

Definition at line 118 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and fillHit().

◆ volEBSD_

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

private

Definition at line 112 of file CaloSteppingAction.cc.

Referenced by update().

◆ volEESD_

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

private

Definition at line 112 of file CaloSteppingAction.cc.

Referenced by update().

◆ volHCSD_

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

private

Definition at line 112 of file CaloSteppingAction.cc.

Referenced by update().

◆ xtalMap_

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

private

Definition at line 113 of file CaloSteppingAction.cc.

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

◆ beginRun()

◆ curve_LY()

◆ fillHit()

◆ fillHits()

◆ fillPassiveHits()

◆ getBirkHC()

◆ getBirkL3()

◆ getDepth()

◆ getDetIDHC()

◆ NaNTrap()

◆ produce()

◆ registerConsumes()

◆ saveHits()

◆ update() [1/4]

◆ update() [2/4]

◆ update() [3/4]

◆ update() [4/4]

Member Data Documentation

◆ allSteps_

◆ birkC1EC_

◆ birkC1HC_

◆ birkC2HC_

◆ birkC3HC_

◆ birkCutEC_

◆ birkSlopeEC_

◆ count_

◆ ebNumberingScheme_

◆ eeNumberingScheme_

◆ eventID_

◆ hcNumberingPS_

◆ hcNumberingScheme_

◆ hitMap_

◆ iC_

◆ mapLV_

◆ nameEBSD_

◆ nameEESD_

◆ nameHCSD_

◆ nameHitC_

◆ nSD_

◆ slave_

◆ slopeLY_

◆ store_

◆ timeSliceUnit_

◆ volEBSD_

◆ volEESD_

◆ volHCSD_

◆ xtalMap_