CaloSteppingAction Class Reference

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

Constructor & Destructor Documentation

CaloSteppingAction()

CaloSteppingAction::CaloSteppingAction

(

const edm::ParameterSet &

p

)

Definition at line 127 of file CaloSteppingAction.cc.

References allSteps_, birkC1EC_, birkC1HC_, birkC2HC_, birkC3HC_, birkCutEC_, birkSlopeEC_, ebNumberingScheme_, eeNumberingScheme_, hcNumberingPS_, hcNumberingScheme_, isotrackApplyRegressor::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 179 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 234 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 488 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 461 of file CaloSteppingAction.cc.

References CaloGVHit::addEnergyDeposit(), hcalRecHitTable_cff::depth, mps_fire::end, eventID_, RemoveAddSevLevel::flag, hitMap_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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 205 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 212 of file CaloSteppingAction.cc.

References gpuPixelDoublets::cc, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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 526 of file CaloSteppingAction.cc.

References funct::abs(), birkC1HC_, birkC2HC_, birkC3HC_, DummyCfis::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 506 of file CaloSteppingAction.cc.

References birkC1EC_, birkCutEC_, birkSlopeEC_, ALCARECOTkAlJpsiMuMu_cff::charge, fastSimProducer_cff::density, CrabHelper::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 478 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 449 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 435 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 191 of file CaloSteppingAction.cc.

References allSteps_, MillePedeFileConverter_cfg::e, fillHits(), fillPassiveHits(), isotrackApplyRegressor::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 184 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 543 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 247 of file CaloSteppingAction.cc.

References allSteps_, PVValHelper::dz, isotrackApplyRegressor::k, mapLV_, mergeVDriftHistosByStation::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 321 of file CaloSteppingAction.cc.

References allSteps_, SplitLinear::begin, mps_fire::end, eventID_, hitMap_, isotrackApplyRegressor::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 333 of file CaloSteppingAction.cc.

References funct::abs(), EcalBaseNumber::addLevel(), allSteps_, filterCSVwithJSON::copy, curve_LY(), hcalRecHitTable_cff::depth, ebNumberingScheme_, eeNumberingScheme_, HBHEDarkening_cff::energy, fillHit(), spr::find(), RemoveAddSevLevel::flag, g, getBirkHC(), getBirkL3(), EcalBaseNumber::getCapacity(), getDepth(), getDetIDHC(), cuy::ii, G4TrackToParticleID::isGammaElectronPositron(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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 429 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 116 of file CaloSteppingAction.cc.

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

birkC1EC_

const double CaloSteppingAction::birkC1EC_

private

Definition at line 117 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkL3().

birkC1HC_

const double CaloSteppingAction::birkC1HC_

private

Definition at line 118 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkHC().

birkC2HC_

const double CaloSteppingAction::birkC2HC_

private

Definition at line 118 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkHC().

birkC3HC_

const double CaloSteppingAction::birkC3HC_

private

Definition at line 119 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkHC().

birkCutEC_

const double CaloSteppingAction::birkCutEC_

private

Definition at line 118 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkL3().

birkSlopeEC_

const double CaloSteppingAction::birkSlopeEC_

private

Definition at line 117 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getBirkL3().

count_

int CaloSteppingAction::count_

private

Definition at line 120 of file CaloSteppingAction.cc.

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

ebNumberingScheme_

std::unique_ptr<EcalBarrelNumberingScheme> CaloSteppingAction::ebNumberingScheme_

private

Definition at line 100 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

eeNumberingScheme_

std::unique_ptr<EcalEndcapNumberingScheme> CaloSteppingAction::eeNumberingScheme_

private

Definition at line 101 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

eventID_

int CaloSteppingAction::eventID_

private

Definition at line 120 of file CaloSteppingAction.cc.

Referenced by fillHit(), and update().

hcNumberingPS_

std::unique_ptr<HcalNumberingFromPS> CaloSteppingAction::hcNumberingPS_

private

Definition at line 102 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getDetIDHC().

hcNumberingScheme_

std::unique_ptr<HcalNumberingScheme> CaloSteppingAction::hcNumberingScheme_

private

Definition at line 107 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and getDetIDHC().

hitMap_

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

private

Definition at line 121 of file CaloSteppingAction.cc.

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

iC_

const edm::ParameterSet CaloSteppingAction::iC_

private

Definition at line 110 of file CaloSteppingAction.cc.

mapLV_

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

private

Definition at line 115 of file CaloSteppingAction.cc.

Referenced by fillPassiveHits(), and update().

nameEBSD_

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

private

Definition at line 111 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

nameEESD_

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

private

Definition at line 111 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

nameHCSD_

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

private

Definition at line 111 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and update().

nameHitC_

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

private

Definition at line 112 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and produce().

nSD_

const int CaloSteppingAction::nSD_ = 3

staticprivate

Definition at line 99 of file CaloSteppingAction.cc.

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

slave_

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

private

Definition at line 108 of file CaloSteppingAction.cc.

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

slopeLY_

const double CaloSteppingAction::slopeLY_

private

Definition at line 117 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and curve_LY().

store_

std::vector<PassiveData> CaloSteppingAction::store_

private

Definition at line 124 of file CaloSteppingAction.cc.

Referenced by fillPassiveHits(), and update().

timeSliceUnit_

const double CaloSteppingAction::timeSliceUnit_

private

Definition at line 119 of file CaloSteppingAction.cc.

Referenced by CaloSteppingAction(), and fillHit().

volEBSD_

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

private

Definition at line 113 of file CaloSteppingAction.cc.

Referenced by update().

volEESD_

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

private

Definition at line 113 of file CaloSteppingAction.cc.

Referenced by update().

volHCSD_

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

private

Definition at line 113 of file CaloSteppingAction.cc.

Referenced by update().

xtalMap_

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

private

Definition at line 114 of file CaloSteppingAction.cc.

Referenced by update().