Inheritance diagram for MaterialBudgetHcalProducer:

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

	MaterialBudgetHcalProducer (const edm::ParameterSet &)

	MaterialBudgetHcalProducer (const MaterialBudgetHcalProducer &)=delete

const MaterialBudgetHcalProducer &	operator= (const MaterialBudgetHcalProducer &)=delete

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

void	registerConsumes (edm::ConsumesCollector) override

	~MaterialBudgetHcalProducer () override=default

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

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

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

void	slotForUpdate (const BeginOfTrack * iT)

virtual	~Observer ()

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

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

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

void	slotForUpdate (const EndOfTrack * iT)

virtual	~Observer ()

Private Member Functions
std::vector< double >	getDDDArray (const std::string &str, const DDsvalues_type &sv)

std::vector< std::string >	getNames (DDFilteredView &fv)

std::vector< std::string >	getNames (cms::DDFilteredView &fv)

bool	isItEC (const std::string &)

bool	isItHF (const G4VTouchable *)

bool	isSensitive (const std::string &)

bool	stopAfter (const G4Step *)

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

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

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

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

Private Attributes
edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord >	cpvTokenDD4hep_

edm::ESGetToken< DDCompactView, IdealGeometryRecord >	cpvTokenDDD_

double	eta_

double	etaMaxP_

double	etaMinP_

bool	fromdd4hep_

std::vector< int >	hfLevels_

std::vector< std::string >	hfNames_

int	id_

double	intLen_

std::vector< double >	intLength_

double	intLens_ [maxSet_]

int	layer_

MaterialAccountingCaloCollection	matcoll_

std::vector< std::string >	matList_

std::string	name_

int	nlayHB_

int	nlayHE_

int	nlayHF_

int	nlayHO_

double	phi_

bool	printSum_

double	radLen_

std::vector< double >	radLength_

double	radLens_ [maxSet_]

double	rMax_

std::vector< std::string >	sensitiveEC_

std::vector< std::string >	sensitives_

double	stepLen_

std::vector< double >	stepLength_

double	stepLens_ [maxSet_]

int	steps_

double	zMax_

Static Private Attributes
static const int	maxSet2_ = 9

static const int	maxSet_ = 25

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 43 of file MaterialBudgetHcalProducer.cc.

Constructor & Destructor Documentation

◆ MaterialBudgetHcalProducer() [1/2]

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

Definition at line 93 of file MaterialBudgetHcalProducer.cc.

References edm::ParameterSet::getUntrackedParameter(), visualization-live-secondInstance_cfg::m, AlCaHLTBitMon_ParallelJobs::p, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                {
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("MaterialBudgetHcalProducer");
   rMax_ = m_p.getUntrackedParameter<double>("RMax", 4.5) * CLHEP::m;
   zMax_ = m_p.getUntrackedParameter<double>("ZMax", 13.0) * CLHEP::m;
   etaMinP_ = m_p.getUntrackedParameter<double>("EtaMinP", 5.2);
   etaMaxP_ = m_p.getUntrackedParameter<double>("EtaMaxP", 0.0);
   fromdd4hep_ = m_p.getUntrackedParameter<bool>("Fromdd4hep", false);
   printSum_ = m_p.getUntrackedParameter<bool>("PrintSummary", false);
   name_ = m_p.getUntrackedParameter<std::string>("Name", "Hcal");
   edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalProducer initialized with rMax " << rMax_ << " mm and zMax "
                                      << zMax_ << " mm printSummary is set to " << printSum_ << " and Fromdd4hep to "
                                      << fromdd4hep_;
 
   produces<MaterialAccountingCaloCollection>(Form("%sMatBCalo", name_.c_str()));
 }

◆ ~MaterialBudgetHcalProducer()

MaterialBudgetHcalProducer::~MaterialBudgetHcalProducer ( )

overridedefault

◆ MaterialBudgetHcalProducer() [2/2]

MaterialBudgetHcalProducer::MaterialBudgetHcalProducer ( const MaterialBudgetHcalProducer & )

delete

Member Function Documentation

◆ beginRun()

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

overridevirtual

Reimplemented from SimWatcher.

Definition at line 125 of file MaterialBudgetHcalProducer.cc.

References ALCARECOTkAlBeamHalo_cff::filter, spr::find(), DDFilteredView::firstChild(), cms::DDFilteredView::get(), edm::EventSetup::getData(), MsgTools::getNames(), mps_fire::i, dqmdumpme::k, personalPlayback::level, DDFilteredView::mergedSpecifics(), Skims_PA_cff::name, names, AlCaHLTBitMon_QueryRunRegistry::string, pfDeepBoostedJetPreprocessParams_cfi::sv, and groupFilesInBlocks::temp.

                                                                  {
   //----- Check that selected volumes are indeed part of the geometry
   // Numbering From DDD
   if (fromdd4hep_) {
     const cms::DDCompactView cpv = es.getData(cpvTokenDD4hep_);
     constexpr int32_t addLevel = 1;
     std::string attribute = "ReadOutName";
     std::string value = "HcalHits";
     const cms::DDFilter filter1(attribute, value);
     cms::DDFilteredView fv1(cpv, filter1);
     std::vector<std::string> names = getNames(fv1);
     for (auto& name : names) {
       std::string namx = (name.find('_') == std::string::npos) ? name : name.substr(0, name.find('_'));
       if (std::find(sensitives_.begin(), sensitives_.end(), namx) == sensitives_.end())
         sensitives_.emplace_back(namx);
     }
     edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer: Names to be tested for " << attribute
                                            << " = " << value << " has " << sensitives_.size() << " elements";
     for (unsigned int i = 0; i < sensitives_.size(); i++)
       edm::LogVerbatim("MaterialBudgetFull")
           << "MaterialBudgetHcalProducer: sensitives[" << i << "] = " << sensitives_[i];
     attribute = "Volume";
     value = "HF";
     const cms::DDFilter filter2(attribute, value);
     cms::DDFilteredView fv2(cpv, filter2);
     std::vector<int> temp = fv2.get<std::vector<int> >("hf", "Levels");
     hfNames_ = getNames(fv2);
     edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer: Names to be tested for " << attribute
                                            << " = " << value << " has " << hfNames_.size() << " elements";
     for (unsigned int i = 0; i < hfNames_.size(); i++) {
       hfLevels_.push_back(temp[i] + addLevel);
       edm::LogVerbatim("MaterialBudgetFull")
           << "MaterialBudgetHcalProducer:  HF[" << i << "] = " << hfNames_[i] << " at level " << hfLevels_[i];
     }
 
     const std::string ecalRO[2] = {"EcalHitsEB", "EcalHitsEE"};
     attribute = "ReadOutName";
     for (int k = 0; k < 2; k++) {
       value = ecalRO[k];
       const cms::DDFilter filter(attribute, value);
       cms::DDFilteredView fv(cpv, filter);
       std::vector<std::string> senstmp = getNames(fv);
       edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer: Names to be tested for " << attribute
                                              << " = " << value << " has " << senstmp.size() << " elements";
       for (unsigned int i = 0; i < senstmp.size(); i++) {
         std::string name = senstmp[i].substr(0, 4);
         if (std::find(sensitiveEC_.begin(), sensitiveEC_.end(), name) == sensitiveEC_.end())
           sensitiveEC_.push_back(name);
       }
     }
     for (unsigned int i = 0; i < sensitiveEC_.size(); i++)
       edm::LogVerbatim("MaterialBudgetFull")
           << "MaterialBudgetHcalProducer:sensitiveEC[" << i << "] = " << sensitiveEC_[i];
   } else {  // if not from dd4hep --> ddd
     const DDCompactView& cpv = es.getData(cpvTokenDDD_);
     constexpr int32_t addLevel = 0;
     std::string attribute = "ReadOutName";
     std::string value = "HcalHits";
     DDSpecificsMatchesValueFilter filter1{DDValue(attribute, value, 0)};
     DDFilteredView fv1(cpv, filter1);
     std::vector<std::string> names = getNames(fv1);
     for (auto& name : names) {
       std::string namx = (name.find('_') == std::string::npos) ? name : name.substr(0, name.find('_'));
       if (std::find(sensitives_.begin(), sensitives_.end(), namx) == sensitives_.end())
         sensitives_.emplace_back(namx);
     }
     edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer: Names to be tested for " << attribute
                                            << " = " << value << " has " << sensitives_.size() << " elements";
     for (unsigned int i = 0; i < sensitives_.size(); i++)
       edm::LogVerbatim("MaterialBudgetFull")
           << "MaterialBudgetHcalProducer: sensitives[" << i << "] = " << sensitives_[i];
     attribute = "Volume";
     value = "HF";
     DDSpecificsMatchesValueFilter filter2{DDValue(attribute, value, 0)};
     DDFilteredView fv2(cpv, filter2);
     hfNames_ = getNames(fv2);
     fv2.firstChild();
     DDsvalues_type sv(fv2.mergedSpecifics());
     std::vector<double> temp = getDDDArray("Levels", sv);
     edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer: Names to be tested for " << attribute
                                            << " = " << value << " has " << hfNames_.size() << " elements";
     for (unsigned int i = 0; i < hfNames_.size(); i++) {
       int level = static_cast<int>(temp[i]);
       hfLevels_.push_back(level + addLevel);
       edm::LogVerbatim("MaterialBudgetFull")
           << "MaterialBudgetHcalProducer:  HF[" << i << "] = " << hfNames_[i] << " at level " << hfLevels_[i];
     }
 
     const std::string ecalRO[2] = {"EcalHitsEB", "EcalHitsEE"};
     attribute = "ReadOutName";
     for (int k = 0; k < 2; k++) {
       value = ecalRO[k];
       DDSpecificsMatchesValueFilter filter3{DDValue(attribute, value, 0)};
       DDFilteredView fv3(cpv, filter3);
       std::vector<std::string> senstmp = getNames(fv3);
       edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer: Names to be tested for " << attribute
                                              << " = " << value << " has " << senstmp.size() << " elements";
       for (unsigned int i = 0; i < senstmp.size(); i++) {
         std::string name = senstmp[i].substr(0, 4);
         if (std::find(sensitiveEC_.begin(), sensitiveEC_.end(), name) == sensitiveEC_.end())
           sensitiveEC_.push_back(name);
       }
     }
     for (unsigned int i = 0; i < sensitiveEC_.size(); i++)
       edm::LogVerbatim("MaterialBudgetFull")
           << "MaterialBudgetHcalProducer:sensitiveEC[" << i << "] = " << sensitiveEC_[i];
   }
 }

◆ getDDDArray()

std::vector< double > MaterialBudgetHcalProducer::getDDDArray	(	const std::string &	str,
		const DDsvalues_type &	sv
	)

private

Definition at line 457 of file MaterialBudgetHcalProducer.cc.

References DDfetch(), Exception, str, pfDeepBoostedJetPreprocessParams_cfi::sv, and relativeConstraints::value.

                                                                                                         {
   edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer:getDDDArray called for " << str;
   DDValue value(str);
   if (DDfetch(&sv, value)) {
     edm::LogVerbatim("MaterialBudgetFull") << value;
     const std::vector<double>& fvec = value.doubles();
     int nval = fvec.size();
     if (nval < 1) {
       throw cms::Exception("MaterialBudgetHcalProducer") << "nval = " << nval << " < 1 for array " << str << "\n";
     }
 
     return fvec;
   } else {
     throw cms::Exception("MaterialBudgetHcalProducer") << "cannot get array " << str << "\n";
   }
 }

◆ getNames() [1/2]

std::vector< std::string > MaterialBudgetHcalProducer::getNames ( DDFilteredView & fv )

private

Definition at line 428 of file MaterialBudgetHcalProducer.cc.

References spr::find(), DDFilteredView::firstChild(), dqm-mbProfile::log, DDFilteredView::logicalPart(), DDFilteredView::next(), AlCaHLTBitMon_QueryRunRegistry::string, and createJobs::tmp.

                                                                             {
   std::vector<std::string> tmp;
   bool dodet = fv.firstChild();
   while (dodet) {
     const DDLogicalPart& log = fv.logicalPart();
     std::string namx = log.name().name();
     if (std::find(tmp.begin(), tmp.end(), namx) == tmp.end())
       tmp.push_back(namx);
     dodet = fv.next();
   }
   return tmp;
 }

◆ getNames() [2/2]

std::vector< std::string > MaterialBudgetHcalProducer::getNames ( cms::DDFilteredView & fv )

private

Definition at line 441 of file MaterialBudgetHcalProducer.cc.

References spr::find(), cms::DDFilteredView::firstChild(), dqmiodumpmetadata::n, cms::DDFilteredView::name(), AlCaHLTBitMon_QueryRunRegistry::string, and createJobs::tmp.

                                                                                {
   std::vector<std::string> tmp;
   const std::vector<std::string> notIn = {
       "CALO", "HCal", "MBBTL", "MBBTR", "MBBTC", "MBAT", "MBBT_R1M", "MBBT_R1P", "MBBT_R1MX", "MBBT_R1PX", "VCAL"};
   while (fv.firstChild()) {
     const std::string n{fv.name().data(), fv.name().size()};
     if (std::find(notIn.begin(), notIn.end(), n) == notIn.end()) {
       std::string::size_type pos = n.find(':');
       const std::string namx = (pos == std::string::npos) ? n : std::string(n, pos + 1, n.size() - 1);
       if (std::find(tmp.begin(), tmp.end(), namx) == tmp.end())
         tmp.push_back(namx);
     }
   }
   return tmp;
 }

◆ isItEC()

bool MaterialBudgetHcalProducer::isItEC ( const std::string & name )

private

Definition at line 499 of file MaterialBudgetHcalProducer.cc.

References Skims_PA_cff::name.

                                                              {
   std::vector<std::string>::const_iterator it = sensitiveEC_.begin();
   std::vector<std::string>::const_iterator itEnd = sensitiveEC_.end();
   for (; it != itEnd; ++it)
     if (name.substr(0, 4) == *it)
       return true;
   return false;
 }

◆ isItHF()

bool MaterialBudgetHcalProducer::isItHF ( const G4VTouchable * touch )

private

Definition at line 484 of file MaterialBudgetHcalProducer.cc.

References jetsAK4_CHS_cff::levels, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                  {
   int levels = ((touch->GetHistoryDepth()) + 1);
   for (unsigned int it = 0; it < hfNames_.size(); it++) {
     if (levels >= hfLevels_[it]) {
       std::string name =
           (static_cast<std::string>(dd4hep::dd::noNamespace(touch->GetVolume(levels - hfLevels_[it])->GetName())))
               .substr(0, 4);
       if (name == hfNames_[it]) {
         return true;
       }
     }
   }
   return false;
 }

◆ isSensitive()

bool MaterialBudgetHcalProducer::isSensitive ( const std::string & name )

private

Definition at line 474 of file MaterialBudgetHcalProducer.cc.

References Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                   {
   std::vector<std::string>::const_iterator it = sensitives_.begin();
   std::vector<std::string>::const_iterator itEnd = sensitives_.end();
   std::string namx = (name.find('_') == std::string::npos) ? name : name.substr(0, name.find('_'));
   for (; it != itEnd; ++it)
     if (namx == *it)
       return true;
   return false;
 }

◆ operator=()

const MaterialBudgetHcalProducer& MaterialBudgetHcalProducer::operator= ( const MaterialBudgetHcalProducer & )

delete

◆ produce()

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

overridevirtual

Implements SimProducer.

Definition at line 117 of file MaterialBudgetHcalProducer.cc.

References MillePedeFileConverter_cfg::e, caloTruthProducer_cfi::hgc, and eostools::move().

                                                                           {
   std::unique_ptr<MaterialAccountingCaloCollection> hgc(new MaterialAccountingCaloCollection);
   for (auto const& mbc : matcoll_) {
     hgc->emplace_back(mbc);
   }
   e.put(std::move(hgc), Form("%sMatBCalo", name_.c_str()));
 }

◆ registerConsumes()

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

overridevirtual

Reimplemented from SimWatcher.

Definition at line 109 of file MaterialBudgetHcalProducer.cc.

References edm::BeginRun, and edm::ConsumesCollector::esConsumes().

                                                                          {
   if (fromdd4hep_)
     cpvTokenDD4hep_ = cc.esConsumes<edm::Transition::BeginRun>();
   else
     cpvTokenDDD_ = cc.esConsumes<edm::Transition::BeginRun>();
   edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalProducer: Initialize the token for CompactView";
 }

◆ stopAfter()

bool MaterialBudgetHcalProducer::stopAfter ( const G4Step * aStep )

private

Definition at line 415 of file MaterialBudgetHcalProducer.cc.

References funct::abs(), findQualityFiles::rr, and geometryCSVtoXML::zz.

                                                               {
   G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
   double rr = hitPoint.perp();
   double zz = std::abs(hitPoint.z());
 
   if (rr > rMax_ || zz > zMax_) {
     edm::LogVerbatim("MaterialBudget") << " MaterialBudgetHcalProducer::StopAfter R = " << rr << " and Z = " << zz;
     return true;
   } else {
     return false;
   }
 }

◆ update() [1/4]

void MaterialBudgetHcalProducer::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 234 of file MaterialBudgetHcalProducer.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

234 { matcoll_.clear(); }

MaterialBudgetHcalProducer::matcoll_

MaterialAccountingCaloCollection matcoll_

Definition: MaterialBudgetHcalProducer.cc:83

◆ update() [2/4]

void MaterialBudgetHcalProducer::update ( const BeginOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 236 of file MaterialBudgetHcalProducer.cc.

References DeadROC_duringRun::dir, mps_fire::i, and sistrip::layer_.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                                {
   const G4Track* aTrack = (*trk)();  // recover G4 pointer if wanted
 
   id_ = layer_ = steps_ = 0;
   radLen_ = intLen_ = stepLen_ = 0;
   nlayHB_ = nlayHE_ = nlayHF_ = nlayHO_ = 0;
   for (int i = 0; i < maxSet_; ++i)
     stepLens_[i] = radLens_[i] = intLens_[i] = 0;
 
   const G4ThreeVector& dir = aTrack->GetMomentum();
   if (dir.theta() != 0) {
     eta_ = dir.eta();
   } else {
     eta_ = -99;
   }
   phi_ = dir.phi();
   edm::LogVerbatim("MaterialBudget") << "Start track with (eta, phi) = (" << eta_ << ", " << phi_ << " Energy "
                                      << aTrack->GetTotalEnergy() << " and ID "
                                      << (aTrack->GetDefinition()->GetPDGEncoding());
 
   if (printSum_) {
     matList_.clear();
     stepLength_.clear();
     radLength_.clear();
     intLength_.clear();
   }
 }

◆ update() [3/4]

void MaterialBudgetHcalProducer::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 264 of file MaterialBudgetHcalProducer.cc.

References funct::abs(), angle_units::operators::convertRadToDeg(), geant_units::operators::convertUnitsTo(), fastSimProducer_cff::density, newFWLiteAna::found, cuy::ii, sistrip::layer_, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, and HLT_2022v12_cff::track.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                            {
   //---------- each step
   G4Material* material = aStep->GetPreStepPoint()->GetMaterial();
   double step = aStep->GetStepLength();
   double radl = material->GetRadlen();
   double intl = material->GetNuclearInterLength();
   double density = convertUnitsTo(1._g_per_cm3, material->GetDensity());
 
   int idOld = id_;
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   std::string name = (static_cast<std::string>(dd4hep::dd::noNamespace(touch->GetVolume(0)->GetName())));
   std::string matName = (static_cast<std::string>(dd4hep::dd::noNamespace(material->GetName())));
   if (printSum_) {
     bool found = false;
     for (unsigned int ii = 0; ii < matList_.size(); ii++) {
       if (matList_[ii] == matName) {
         stepLength_[ii] += step;
         radLength_[ii] += (step / radl);
         intLength_[ii] += (step / intl);
         found = true;
         break;
       }
     }
     if (!found) {
       matList_.push_back(matName);
       stepLength_.push_back(step);
       radLength_.push_back(step / radl);
       intLength_.push_back(step / intl);
     }
     if ((std::abs(eta_) >= etaMinP_) && (std::abs(eta_) <= etaMaxP_))
       edm::LogVerbatim("MaterialBudget") << "Volume " << name << " id " << id_ << ":" << idOld << " Step " << step
                                          << " Material " << matName << " Old Length " << stepLen_ << " X0 "
                                          << step / radl << ":" << radLen_ << " Lambda " << step / intl << ":"
                                          << intLen_;
   } else {
     if ((std::abs(eta_) >= etaMinP_) && (std::abs(eta_) <= etaMaxP_))
       edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalProducer: Step at " << name << " id " << id_ << ":"
                                          << idOld << " Length " << step << " in " << matName << " of density "
                                          << density << " g/cc; Radiation Length " << radl << " mm; Interaction Length "
                                          << intl << " mm\n                          Position "
                                          << aStep->GetPreStepPoint()->GetPosition() << " Cylindrical R "
                                          << aStep->GetPreStepPoint()->GetPosition().perp() << " Length (so far) "
                                          << stepLen_ << " L/X0 " << step / radl << "/" << radLen_ << " L/Lambda "
                                          << step / intl << "/" << intLen_;
   }
 
   int det = 0, lay = 0;
   double abseta = std::abs(eta_);
   edm::LogVerbatim("MaterialBudgetFull") << "Volume " << name << ":" << matName << " EC:Sensitive:HF " << isItEC(name)
                                          << ":" << isSensitive(name) << ":" << isItHF(touch) << " Eta " << abseta
                                          << " HC " << ((touch->GetReplicaNumber(1)) / 1000) << ":"
                                          << ((touch->GetReplicaNumber(0) / 10) % 100 + 3) << " X0 "
                                          << (radLen_ + (step / radl)) << " Lambda " << (intLen_ + (step / intl));
 
   if (isItEC(name)) {
     det = 1;
     lay = 1;
   } else {
     if (isSensitive(name)) {
       if (isItHF(touch)) {
         det = 5;
         lay = 21;
         if (lay != layer_)
           ++nlayHF_;
       } else {
         det = (touch->GetReplicaNumber(1)) / 1000;
         lay = (touch->GetReplicaNumber(0) / 10) % 100 + 3;
         if (det == 4) {
           if (abseta < 1.479)
             lay = layer_ + 1;
           else
             lay--;
           if (lay < 3)
             lay = 3;
           if (lay == layer_)
             lay++;
           if (lay > 20)
             lay = 20;
           if (lay != layer_)
             ++nlayHE_;
         } else if (lay != layer_) {
           if (lay >= 20)
             ++nlayHO_;
           else
             ++nlayHB_;
         }
       }
       edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalProducer: Det " << det << " Layer " << lay << " Eta "
                                              << eta_ << " Phi " << convertRadToDeg(phi_);
     } else if (layer_ == 1) {
       det = -1;
       lay = 2;
     }
   }
   if (det != 0) {
     if (lay != layer_) {
       id_ = lay;
       layer_ = lay;
     }
   }
 
   if (id_ > idOld) {
     if ((abseta >= etaMinP_) && (abseta <= etaMaxP_))
       edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalProducer: Step at " << name << " calls filHisto with "
                                          << (id_ - 1);
     stepLens_[id_ - 1] = stepLen_;
     radLens_[id_ - 1] = radLen_;
     intLens_[id_ - 1] = intLen_;
   }
 
   stepLen_ += step;
   radLen_ += (step / radl);
   intLen_ += (step / intl);
   if (id_ == 21) {
     if (!isItHF(aStep->GetPostStepPoint()->GetTouchable())) {
       if ((abseta >= etaMinP_) && (abseta <= etaMaxP_))
         edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalProducer: After HF in " << name << ":"
                                            << static_cast<std::string>(dd4hep::dd::noNamespace(
                                                   aStep->GetPostStepPoint()->GetTouchable()->GetVolume(0)->GetName()))
                                            << " calls fillHisto with " << id_;
       stepLens_[idOld] = stepLen_;
       radLens_[idOld] = radLen_;
       intLens_[idOld] = intLen_;
       ++id_;
       layer_ = 0;
     }
   }
 
   //----- Stop tracking after selected position
   if (stopAfter(aStep)) {
     G4Track* track = aStep->GetTrack();
     track->SetTrackStatus(fStopAndKill);
   }
 }

◆ update() [4/4]

void MaterialBudgetHcalProducer::update ( const EndOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack *>.

Definition at line 399 of file MaterialBudgetHcalProducer.cc.

References dqmdumpme::k, MaterialAccountingCalo::m_eta, MaterialAccountingCalo::m_intLen, MaterialAccountingCalo::m_layers, MaterialAccountingCalo::m_phi, MaterialAccountingCalo::m_radLen, and MaterialAccountingCalo::m_stepLen.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                              {
   MaterialAccountingCalo matCalo;
   matCalo.m_eta = eta_;
   matCalo.m_phi = phi_;
   for (int k = 0; k < maxSet_; ++k) {
     matCalo.m_stepLen.emplace_back(stepLens_[k]);
     matCalo.m_radLen.emplace_back(radLens_[k]);
     matCalo.m_intLen.emplace_back(intLens_[k]);
   }
   matCalo.m_layers.emplace_back(nlayHB_);
   matCalo.m_layers.emplace_back(nlayHE_);
   matCalo.m_layers.emplace_back(nlayHO_);
   matCalo.m_layers.emplace_back(nlayHF_);
   matcoll_.emplace_back(matCalo);
 }