Inheritance diagram for SimG4HGCalValidation:

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

	SimG4HGCalValidation (const edm::ParameterSet &p)

	~SimG4HGCalValidation () override

Public Member Functions inherited from SimProducer
void	registerProducts (edm::ProducerBase &iProd)

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

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

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

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

Private Member Functions
void	clear ()

void	init ()

void	layerAnalysis (PHGCalValidInfo &)

const SimG4HGCalValidation &	operator= (const SimG4HGCalValidation &)

	SimG4HGCalValidation (const SimG4HGCalValidation &)

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

Private Attributes
unsigned int	count_

std::vector< unsigned int >	dets_

std::vector< int >	detTypes_

double	edepEE_

double	edepHEB_

double	edepHEF_

std::vector< HGCalNumberingScheme * >	hgcalNumbering_

std::vector< double >	hgcEEedep_

std::vector< double >	hgcHEBedep_

std::vector< double >	hgcHEFedep_

std::vector< unsigned int >	hgchitDets_

std::vector< unsigned int >	hgchitIndex_

std::vector< double >	hgchitX_

std::vector< double >	hgchitY_

std::vector< double >	hgchitZ_

std::vector< HGCNumberingScheme * >	hgcNumbering_

std::string	labelLayer_

int	levelT1_

int	levelT2_

std::vector< std::string >	names_

HcalNumberingFromDDD *	numberingFromDDD_

std::vector< int >	subdet_

std::vector< int >	types_

int	verbosity_

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 56 of file SimG4HGCalValidation.cc.

Constructor & Destructor Documentation

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

Definition at line 107 of file SimG4HGCalValidation.cc.

References detTypes_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gen::k, labelLayer_, names_, AlCaHLTBitMon_QueryRunRegistry::string, types_, and verbosity_.

                                                                   : 
   numberingFromDDD_(nullptr), levelT1_(999), levelT2_(999), count_(0) {
 
   edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("SimG4HGCalValidation");
   names_        = m_Anal.getParameter<std::vector<std::string> >("Names");
   types_        = m_Anal.getParameter<std::vector<int> >("Types");
   detTypes_     = m_Anal.getParameter<std::vector<int> >("DetTypes");
   labelLayer_   = m_Anal.getParameter<std::string>("LabelLayerInfo");
   verbosity_    = m_Anal.getUntrackedParameter<int>("Verbosity",0);
 
   produces<PHGCalValidInfo>(labelLayer_);
 
   edm::LogVerbatim("ValidHGCal") << "HGCalTestAnalysis:: Initialised as "
                  << "observer of begin events and of G4step "
                  << "with Parameter values: \n\tLabel : " 
                  << labelLayer_ << " and with " 
                  << names_.size() << " detectors";
   for (unsigned int k=0; k<names_.size(); ++k)
     edm::LogVerbatim("ValidHGCal") << " [" << k << "] " << names_[k] 
                    << " Type " << types_[k] << " DetType "
                    << detTypes_[k];
 } 

SimG4HGCalValidation::~SimG4HGCalValidation ( )

override

Definition at line 130 of file SimG4HGCalValidation.cc.

References hgcalNumbering_, hgcNumbering_, and numberingFromDDD_.

                                             {
   delete numberingFromDDD_;
   for (auto number : hgcNumbering_)   delete number;
   for (auto number : hgcalNumbering_) delete number;
 }

SimG4HGCalValidation::SimG4HGCalValidation ( const SimG4HGCalValidation & )

private

Member Function Documentation

void SimG4HGCalValidation::clear ( void )

private

Definition at line 376 of file SimG4HGCalValidation.cc.

References DEFINE_SIMWATCHER, hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, and hgchitZ_.

Referenced by update().

                                  {
 
   hgchitDets_.erase(hgchitDets_.begin(),hgchitDets_.end()); 
   hgchitIndex_.erase(hgchitIndex_.begin(),hgchitIndex_.end()); 
   hgchitX_.erase(hgchitX_.begin(),hgchitX_.end()); 
   hgchitY_.erase(hgchitY_.begin(),hgchitY_.end()); 
   hgchitZ_.erase(hgchitZ_.begin(),hgchitZ_.end());      
 }

void SimG4HGCalValidation::init ( )

private

void SimG4HGCalValidation::layerAnalysis ( PHGCalValidInfo & product )

private

Definition at line 360 of file SimG4HGCalValidation.cc.

References edepEE_, edepHEB_, edepHEF_, PHGCalValidInfo::fillhgcHits(), PHGCalValidInfo::fillhgcLayers(), hgcEEedep_, hgcHEBedep_, hgcHEFedep_, hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, hgchitZ_, MeV, and verbosity_.

Referenced by produce().

                                                                  {
 
   if (verbosity_ > 0)
     edm::LogVerbatim("ValidHGCal") << "\n ===>>> SimG4HGCalValidation: Energy "
                    << "deposit\n at EE : " << std::setw(6) 
                    << edepEE_/CLHEP::MeV << "\n at HEF: " 
                    << std::setw(6) << edepHEF_/CLHEP::MeV 
                    << "\n at HEB: " << std::setw(6)
                    << edepHEB_/CLHEP::MeV;
   
   //Fill HGC Variables
   product.fillhgcHits(hgchitDets_,hgchitIndex_, hgchitX_, hgchitY_, hgchitZ_);
   product.fillhgcLayers(edepEE_,edepHEF_,edepHEB_,hgcEEedep_, hgcHEFedep_,hgcHEBedep_);
 }

const SimG4HGCalValidation& SimG4HGCalValidation::operator= ( const SimG4HGCalValidation & )

private

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

overridevirtual

Implements SimProducer.

Definition at line 136 of file SimG4HGCalValidation.cc.

References labelLayer_, layerAnalysis(), eostools::move(), and edm::Event::put().

                                                                     {
 
   std::unique_ptr<PHGCalValidInfo> productLayer(new PHGCalValidInfo);
   layerAnalysis(*productLayer);
   e.put(std::move(productLayer),labelLayer_);
 }

void SimG4HGCalValidation::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 143 of file SimG4HGCalValidation.cc.

References dets_, detTypes_, Exception, DetId::Forward, ForwardEmpty, HGCalDDDConstants::geomMode(), edm::EventSetup::get(), DetId::Hcal, HcalEndcap, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, hgcalNumbering_, HGCEE, hgcEEedep_, HGCHEB, hgcHEBedep_, HGCHEF, hgcHEFedep_, hgcNumbering_, mps_fire::i, edm::ESHandleBase::isValid(), LayerTriplets::layers(), HGCalDDDConstants::layers(), levelT1_, levelT2_, HGCalDDDConstants::levelTop(), numberingFromDDD_, edm::ESHandle< T >::product(), subdet_, HGCalGeometryMode::Trapezoid, and types_.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                         {
 
   const edm::EventSetup* es = (*job)();
   for (unsigned int type=0; type<types_.size(); ++type) {
     int      layers(0);
     int      detType = detTypes_[type];
     G4String nameX = "HGCal";
     if (types_[type] <= 1) {
       if (types_[type] == 0) {
     subdet_.emplace_back(ForwardSubdetector::ForwardEmpty);
     if      (detType == 0) dets_.emplace_back((int)(DetId::HGCalEE));
     else if (detType == 1) dets_.emplace_back((int)(DetId::HGCalHSi));
     else                   dets_.emplace_back((int)(DetId::HGCalHSc));
       } else {
     dets_.push_back((unsigned int)(DetId::Forward));
     if      (detType == 0) subdet_.emplace_back((int)(ForwardSubdetector::HGCEE));
     else if (detType == 1) subdet_.emplace_back((int)(ForwardSubdetector::HGCHEF));
     else                   subdet_.emplace_back((int)(ForwardSubdetector::HGCHEB));
       }
       if        (detType == 0) nameX = "HGCalEESensitive";
       else if   (detType == 1) nameX = "HGCalHESiliconSensitive";
       else                     nameX = "HGCalHEScintillatorSensitive";
       edm::ESHandle<HGCalDDDConstants>    hdc;
       es->get<IdealGeometryRecord>().get(nameX,hdc);
       if (hdc.isValid()) {
     HGCalGeometryMode::GeometryMode m_mode = hdc->geomMode();
     levelT1_ = hdc->levelTop(0); 
     levelT2_ = hdc->levelTop(1); 
     if (m_mode == HGCalGeometryMode::Trapezoid) {
       types_[type] =-1;
       hgcalNumbering_.emplace_back(new HGCalNumberingScheme(*hdc,(DetId::Detector)(dets_[type]),nameX));
     } else if ((m_mode == HGCalGeometryMode::Hexagon) ||
            (m_mode == HGCalGeometryMode::HexagonFull)) {
       types_[type] = 1;
       hgcNumbering_.push_back(new HGCNumberingScheme(*hdc,nameX));
     } else {
       types_[type] = 0;
       hgcalNumbering_.emplace_back(new HGCalNumberingScheme(*hdc,(DetId::Detector)(dets_[type]),nameX));
     }
     layers = hdc->layers(false);
       } else {
     edm::LogError("ValidHGCal") << "Cannot find HGCalDDDConstants for "
                     << nameX;
     throw cms::Exception("Unknown", "ValidHGCal") 
       << "Cannot find HGCalDDDConstants for " << nameX << "\n";
       }
     } else {
       nameX = "HcalEndcap";
       dets_.push_back((unsigned int)(DetId::Hcal));
       subdet_.push_back((int)(HcalSubdetector::HcalEndcap));
       edm::ESHandle<HcalDDDSimConstants>    hdc;
       es->get<HcalSimNumberingRecord>().get(hdc);
       if (hdc.isValid()) {
     numberingFromDDD_ = new HcalNumberingFromDDD(hdc.product());
     layers = 18;
       } else {
     edm::LogError("ValidHGCal") << "Cannot find HcalDDDSimConstant";
     throw cms::Exception("Unknown", "ValidHGCal") 
       << "Cannot find HcalDDDSimConstant\n";
       }
     }
     if        (detType == 0) {
       for (int i=0; i<layers; ++i) hgcEEedep_.push_back(0);
     } else if (detType == 1) {
       for (int i=0; i<layers; ++i) hgcHEFedep_.push_back(0);
     } else                   {
       for (int i=0; i<layers; ++i) hgcHEBedep_.push_back(0);
     }
     edm::LogVerbatim("ValidHGCal") << "[" << type << "]: " << nameX << " det "
                    << dets_[type] << " subdet " 
                    << subdet_[type] << " with " << layers 
                    << " layers";
   }
 }

void SimG4HGCalValidation::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 219 of file SimG4HGCalValidation.cc.

References clear(), count_, edepEE_, edepHEB_, edepHEF_, hgcEEedep_, hgcHEBedep_, hgcHEFedep_, mps_fire::i, and verbosity_.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                           {
  
   int iev = (*evt)()->GetEventID();
   if (verbosity_ > 0) 
     edm::LogVerbatim("ValidHGCal") << "SimG4HGCalValidation: =====> Begin "
                    << "event = " << iev;
   
   ++count_;
   edepEE_ = edepHEF_ = edepHEB_ = 0.;
 
   //HGCal variables 
   for (unsigned int i = 0; i<hgcEEedep_.size();  i++) hgcEEedep_[i]  = 0.;
   for (unsigned int i = 0; i<hgcHEFedep_.size(); i++) hgcHEFedep_[i] = 0.; 
   for (unsigned int i = 0; i<hgcHEBedep_.size(); i++) hgcHEBedep_[i] = 0.; 
 
   //Cache reset  
   clear();
 }

void SimG4HGCalValidation::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 239 of file SimG4HGCalValidation.cc.

References TauDecayModes::dec, HcalNumberingFromDDD::HcalID::depth, egammaForCoreTracking_cff::depth, dets_, detTypes_, edepEE_, edepHEB_, edepHEF_, HcalNumberingFromDDD::HcalID::etaR, spr::find(), hgcalNumbering_, hgcEEedep_, hgcHEBedep_, hgcHEFedep_, hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, hgchitZ_, hgcNumbering_, gen::k, HcalNumberingFromDDD::HcalID::lay, levelT1_, levelT2_, dataset::name, names_, numberingFromDDD_, HcalTestNumbering::packHcalIndex(), HcalNumberingFromDDD::HcalID::phis, HcalNumberingFromDDD::HcalID::subdet, subdet_, tmp, types_, HcalNumberingFromDDD::unitID(), verbosity_, mps_merge::weight, and HcalNumberingFromDDD::HcalID::zside.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                       {
 
   if (aStep != nullptr) {
     G4VPhysicalVolume* curPV  = aStep->GetPreStepPoint()->GetPhysicalVolume();
     const G4String& name = curPV->GetName();
     G4VSensitiveDetector* curSD = aStep->GetPreStepPoint()->GetSensitiveDetector();
     if (verbosity_ > 1)
       edm::LogVerbatim("ValidHGCal") << "ValidHGCal::Step in " << name <<" at "
                      << aStep->GetPreStepPoint()->GetPosition();
 
     // Only for Sensitive detector
     if (curSD != nullptr) {
       int type(-1);
       for (unsigned int k=0; k<names_.size(); ++k) {
     if (name.find(names_[k].c_str()) != std::string::npos) {
       type = k; break;
     }
       }
       int detType = (type >= 0) ? detTypes_[type] : -1;
       if (verbosity_ > 0)
     edm::LogVerbatim("ValidHGCal") << "ValidHGCal::In SD " << name 
                        << " type " << type << ":" << detType;
 
       // Right type of SD
       if (type >= 0) {
     //Get the 32-bit index of the hit cell
     const G4TouchableHistory* touchable = static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
     unsigned int index(0);
     int          layer(0);
     G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
     if (types_[type] <= 1) {
       // HGCal
       G4ThreeVector localpos = touchable->GetHistory()->GetTopTransform().TransformPoint(hitPoint);   
       float globalZ = touchable->GetTranslation(0).z();
       int iz(globalZ>0 ? 1 : -1);
       int module(-1), cell(-1);
       if (types_[type] == 1) {
         if (touchable->GetHistoryDepth() == levelT1_) {
           layer  = touchable->GetReplicaNumber(0);
         } else {
           layer  = touchable->GetReplicaNumber(2);
           module = touchable->GetReplicaNumber(1);
           cell   = touchable->GetReplicaNumber(0);
         }
         index = hgcNumbering_[type]->getUnitID((ForwardSubdetector)(subdet_[type]), layer, module, cell, iz, localpos);
       } else {
         if ((touchable->GetHistoryDepth() == levelT1_) || 
         (touchable->GetHistoryDepth() == levelT2_)) {
           layer  = touchable->GetReplicaNumber(0);
         } else {
           layer  = touchable->GetReplicaNumber(3);
           module = touchable->GetReplicaNumber(2);
           cell   = touchable->GetReplicaNumber(1);
         }
         double weight(0);
         index = hgcalNumbering_[type]->getUnitID(layer, module, cell, iz,
                              hitPoint, weight);
       }
       if (verbosity_ > 1)
         edm::LogVerbatim("ValidHGCal") << "HGCal: " << name << " Layer "
                        << layer << " Module " << module
                        << " Cell " << cell;
     } else {
       // Hcal
       int depth = (touchable->GetReplicaNumber(0))%10 + 1;
       int lay   = (touchable->GetReplicaNumber(0)/10)%100 + 1;
       int det   = (touchable->GetReplicaNumber(1))/1000;
       HcalNumberingFromDDD::HcalID tmp = numberingFromDDD_->unitID(det, hitPoint, depth, lay);
       index = HcalTestNumbering::packHcalIndex(tmp.subdet,tmp.zside,tmp.depth,tmp.etaR,tmp.phis,tmp.lay);
       layer = tmp.lay;
       if (verbosity_ > 1)
         edm::LogVerbatim("ValidHGCal") << "HCAL: " << det << ":" << depth 
                        << ":" << lay << " o/p " 
                        << tmp.subdet << ":" << tmp.zside 
                        << ":" << tmp.depth << ":" 
                        << tmp.etaR << ":" << tmp.phis 
                        << ":" << tmp.lay << " point "
                        << hitPoint << " " << hitPoint.rho()
                        << ":"  << hitPoint.eta() << ":" 
                        << hitPoint.phi();
     }
 
     double edeposit = aStep->GetTotalEnergyDeposit();
     if (verbosity_ > 0)
       edm::LogVerbatim("ValidHGCal") << "Layer " << layer << " Index " 
                      << std::hex << index << std::dec
                      << " Edep " << edeposit << " hit at "
                      << hitPoint;
     if        (detType == 0) {
       edepEE_  += edeposit;
       if (layer < (int)(hgcEEedep_.size()))
         hgcEEedep_[layer]  += edeposit;
     } else if (detType == 1) {
       edepHEF_ += edeposit;
       if (layer < (int)(hgcHEFedep_.size())) 
         hgcHEFedep_[layer] += edeposit;
     } else                   {
       edepHEB_ += edeposit;
       if (layer < (int)(hgcHEBedep_.size())) 
         hgcHEBedep_[layer] += edeposit;
     }
     G4String nextVolume("XXX");
     if (aStep->GetTrack()->GetNextVolume()!=nullptr)
       nextVolume = aStep->GetTrack()->GetNextVolume()->GetName();
 
     if (nextVolume.c_str()!=name.c_str()) { //save hit when it exits cell
       if (std::find(hgchitIndex_.begin(),hgchitIndex_.end(),index) == hgchitIndex_.end()) {
         hgchitDets_.push_back(dets_[type]);
         hgchitIndex_.push_back(index);
         hgchitX_.push_back(hitPoint.x());
         hgchitY_.push_back(hitPoint.y()); 
         hgchitZ_.push_back(hitPoint.z());
       }
     }
       } // it is right type of SD
     } // it is in a SD
   }//if aStep!=NULL
 }//end update aStep