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

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

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

private

Definition at line 103 of file CaloSteppingAction.h.

Constructor & Destructor Documentation

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

Definition at line 33 of file CaloSteppingAction.cc.

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_, slave_, slopeLY_, and timeSliceUnit_.

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

CaloSteppingAction::~CaloSteppingAction ( )

override

Definition at line 83 of file CaloSteppingAction.cc.

References count_.

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

Member Function Documentation

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

private

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

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

private

Definition at line 337 of file CaloSteppingAction.cc.

References CaloGVHit::addEnergyDeposit(), end, eventID_, RemoveAddSevLevel::flag, hitMap_, CaloGVHit::setEventID(), CaloGVHit::setID(), 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;
   }
 }

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

private

Definition at line 102 of file CaloSteppingAction.cc.

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

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

private

Definition at line 109 of file CaloSteppingAction.cc.

References 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));
       cc.emplace_back(hit);
     }
   }
 }

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

private

Definition at line 402 of file CaloSteppingAction.cc.

References funct::abs(), birkC1HC_, birkC2HC_, birkC3HC_, HltBtagPostValidation_cff::c, 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;
 }

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

private

Definition at line 382 of file CaloSteppingAction.cc.

References birkC1EC_, birkCutEC_, birkSlopeEC_, 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;
 }

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

private

Definition at line 354 of file CaloSteppingAction.cc.

References LEDCalibrationChannels::depth, 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;
 }

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

private

Definition at line 325 of file CaloSteppingAction.cc.

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

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

private

Definition at line 311 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";
   }
 }

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

overridevirtual

Implements SimProducer.

Definition at line 88 of file CaloSteppingAction.cc.

References allSteps_, fillHits(), fillPassiveHits(), dqmdumpme::k, eostools::move(), nameHitC_, nSD_, edm::Event::put(), 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");
   }
 }

void CaloSteppingAction::saveHits ( int flag )

private

Definition at line 419 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());
   }
 }

void CaloSteppingAction::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 128 of file CaloSteppingAction.cc.

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

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

void CaloSteppingAction::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 143 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_.

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

void CaloSteppingAction::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 202 of file CaloSteppingAction.cc.

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

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

void CaloSteppingAction::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 214 of file CaloSteppingAction.cc.

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(), AnalysisDataFormats_SUSYBSMObjects::hc, cuy::ii, G4TrackToParticleID::isGammaElectronPositron(), crabWrapper::key, DTRecHitClients_cfi::local, mapLV_, MeV, NaNTrap(), EcalBaseNumber::setSize(), findQualityFiles::size, store_, ntuplemaker::time, volEBSD_, volEESD_, volHCSD_, and xtalMap_.

                                                    {
   //  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());
   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;
       auto 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));
         }
       }
       auto 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);
     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);
 #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;
 #endif
     uint32_t copy = 0;
     if (((aStep->GetPostStepPoint() == nullptr) || (aStep->GetTrack()->GetNextVolume() == nullptr)) &&
         (aStep->IsLastStepInVolume())) {
       energy += (aStep->GetPreStepPoint()->GetKineticEnergy() / CLHEP::MeV);
     } else {
       time = aStep->GetPostStepPoint()->GetGlobalTime() / CLHEP::nanosecond;
       copy = (touch->GetHistoryDepth() < 1)
                  ? static_cast<uint32_t>(touch->GetReplicaNumber(0))
                  : static_cast<uint32_t>(touch->GetReplicaNumber(0) + 1000 * touch->GetReplicaNumber(1));
     }
     if (it != mapLV_.end()) {
       PassiveData key(std::make_tuple(lv, copy, trackId, pdg, time, energy, energy, stepl));
       store_.push_back(key);
     }
   }
 }

void CaloSteppingAction::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 305 of file CaloSteppingAction.cc.

References count_.

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