HGCalTBMB Class Reference

#include <HGCalTBMB.h>

Inheritance diagram for HGCalTBMB:

Public Member Functions
	HGCalTBMB (const edm::ParameterSet &)
	~HGCalTBMB () override
Public Member Functions inherited from SimWatcher
	SimWatcher ()
virtual	~SimWatcher ()
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
int	findVolume (const G4VTouchable *touch, bool stop) const
	HGCalTBMB (const HGCalTBMB &)=delete
const HGCalTBMB &	operator= (const HGCalTBMB &)=delete
bool	stopAfter (const G4Step *)
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
std::vector< double >	intLen_
std::vector< std::string >	listNames_
std::vector< TH1D * >	me100_
std::vector< TH1D * >	me200_
std::vector< TH1D * >	me300_
unsigned int	nList_
std::vector< double >	radLen_
std::vector< double >	stepLen_
std::string	stopName_
double	stopZ_

Additional Inherited Members

Detailed Description

Definition at line 17 of file HGCalTBMB.h.

Constructor & Destructor Documentation

HGCalTBMB::HGCalTBMB

(

const edm::ParameterSet &

p

)

Definition at line 19 of file HGCalTBMB.cc.

References edm::ParameterSet::getParameter(), create_public_lumi_plots::hist, mps_fire::i, edm::Service< T >::isAvailable(), gen::k, listNames_, TFileService::make(), me100_, me200_, me300_, dataset::name, nList_, stopName_, stopZ_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, and fftjetcommon_cfi::title.

                                             {
  
  edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("HGCalTBMB");
  listNames_   = m_p.getParameter<std::vector<std::string> >("DetectorNames");
  stopName_    = m_p.getParameter<std::string>("StopName");
  stopZ_       = m_p.getParameter<double>("MaximumZ");
  nList_       = listNames_.size();
  edm::LogInfo("HGCSim") << "HGCalTBMB initialized for " << nList_ <<" volumes\n";
  for (unsigned int k=0; k<nList_; ++k)
    edm::LogInfo("HGCSim") << " [" << k << "] " << listNames_[k] << std::endl;
  edm::LogInfo("HGCSim") << "Stop after " << stopZ_ << " or reaching volume "
             << stopName_ << std::endl;

  edm::Service<TFileService> tfile;
  if ( !tfile.isAvailable() )
    throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                      << "please add it to config file";
  char  name[20], title[80];
  TH1D* hist;
  for (unsigned int i=0; i<=nList_; i++) {
    std::string named = (i == nList_) ? "Total" : listNames_[i];
    sprintf(name, "RadL%d", i);
    sprintf(title, "MB(X0) for (%s)", named.c_str());
    hist = tfile->make<TH1D>(name,title,1000,0.0,100.0);
    hist->Sumw2(true);
    me100_.push_back(hist);
    sprintf(name, "IntL%d", i);
    sprintf(title, "MB(L0) for (%s)", named.c_str());
    hist = tfile->make<TH1D>(name,title,1000,0.0,10.0);
    hist->Sumw2(true);
    me200_.push_back(hist);
    sprintf(name, "StepL%d", i);
    sprintf(title, "MB(Step) for (%s)", named.c_str());
    hist = tfile->make<TH1D>(name,title,1000,0.0,50000.0);
    hist->Sumw2(true);
    me300_.push_back(hist);
  }
  edm::LogInfo("HGCSim") << "HGCalTBMB: Booking user histos done ===";
}

HGCalTBMB::~HGCalTBMB ( )

override

Definition at line 59 of file HGCalTBMB.cc.

{ }

HGCalTBMB::HGCalTBMB ( const HGCalTBMB &  )

privatedelete

Member Function Documentation

int HGCalTBMB::findVolume ( const G4VTouchable *  touch, bool  stop  ) const

private

Definition at line 143 of file HGCalTBMB.cc.

References cuy::ii, gen::k, hcalDigis_cfi::level, listNames_, dataset::name, nList_, stopName_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by stopAfter(), and update().

                                                                    {

  int ivol  =-1;
  int level = (touch->GetHistoryDepth())+1;
  for (int ii = 0; ii < level; ii++) {
    std::string name = touch->GetVolume(ii)->GetName();
    if (stop) {
      if (strcmp(name.c_str(),stopName_.c_str()) == 0) ivol = 0;
    } else {
      for (unsigned int k=0; k<nList_; ++k) {
    if (strcmp(name.c_str(),listNames_[k].c_str()) == 0) {
      ivol = k; break;
    }
      }
    }
    if (ivol >= 0) break;
  }
  return ivol;
}

const HGCalTBMB& HGCalTBMB::operator= ( const HGCalTBMB &  )

privatedelete

bool HGCalTBMB::stopAfter ( const G4Step *  aStep )

private

Definition at line 126 of file HGCalTBMB.cc.

References gather_cfg::cout, findVolume(), RemoveAddSevLevel::flag, stopZ_, and geometryCSVtoXML::zz.

Referenced by update().

                                             {

  bool   flag(false);
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  G4ThreeVector hitPoint    = aStep->GetPreStepPoint()->GetPosition();
  if (aStep->GetPostStepPoint() != nullptr) 
    hitPoint = aStep->GetPostStepPoint()->GetPosition();
  double zz    = hitPoint.z();

  if ((findVolume(touch,true) == 0) || (zz > stopZ_)) flag = true;
#ifdef EDM_ML_DEBUG
  std::cout << " HGCalTBMB::Name " << touch->GetVolume(0)->GetName() << " z "
        << zz << " Flag" << flag << std::endl;
#endif
  return flag;
}

void HGCalTBMB::update ( const BeginOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 62 of file HGCalTBMB.cc.

References gather_cfg::cout, GeV, createfilelist::int, intLen_, nList_, radLen_, and stepLen_.

Referenced by progressbar.ProgressBar::__next__(), 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(), progressbar.ProgressBar::finish(), 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().

                                              {

  radLen_  = std::vector<double>(nList_+1,0);
  intLen_  = std::vector<double>(nList_+1,0);
  stepLen_ = std::vector<double>(nList_+1,0);

#ifdef EDM_ML_DEBUG
  const G4Track *   aTrack = (*trk)(); // recover G4 pointer if wanted
  const G4ThreeVector& mom = aTrack->GetMomentum() ;
  double         theEnergy = aTrack->GetTotalEnergy();
  int            theID     = (int)(aTrack->GetDefinition()->GetPDGEncoding());
  std::cout << "MaterialBudgetHcalHistos: Track " << aTrack->GetTrackID()
        << " Code " << theID << " Energy " <<theEnergy/CLHEP::GeV
        << " GeV; Momentum " << mom << std::endl;
#endif
}

void HGCalTBMB::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 79 of file HGCalTBMB.cc.

References gather_cfg::cout, findVolume(), intLen_, nList_, radLen_, stepLen_, stopAfter(), and HiIsolationCommonParameters_cff::track.

Referenced by progressbar.ProgressBar::__next__(), 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(), progressbar.ProgressBar::finish(), 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().

                                          {

  G4Material * material = aStep->GetPreStepPoint()->GetMaterial();
  double step    = aStep->GetStepLength();
  double radl    = material->GetRadlen();
  double intl    = material->GetNuclearInterLength();

  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  int indx = findVolume(touch,false);

  if (indx >= 0) {
    stepLen_[indx] += step;
    radLen_[indx]  += (step/radl);
    intLen_[indx]  += (step/intl);
  }
  stepLen_[nList_] += step;
  radLen_[nList_]  += (step/radl);
  intLen_[nList_]  += (step/intl);
#ifdef EDM_ML_DEBUG
  std::cout << "HGCalTBMB::Step in "
        << touch->GetVolume(0)->GetLogicalVolume()->GetName()
        << " Index " << indx <<" Step " << step << " RadL " << step/radl
        << " IntL " << step/intl << std::endl;
#endif

  if (stopAfter(aStep)) {
    G4Track* track = aStep->GetTrack();
    track->SetTrackStatus(fStopAndKill);
  }
}

void HGCalTBMB::update ( const EndOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 110 of file HGCalTBMB.cc.

References gather_cfg::cout, cuy::ii, intLen_, listNames_, me100_, me200_, me300_, dataset::name, nList_, radLen_, stepLen_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by progressbar.ProgressBar::__next__(), 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(), progressbar.ProgressBar::finish(), 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().

                                            {

  for (unsigned int ii=0; ii<=nList_; ++ii) {
    me100_[ii]->Fill(radLen_[ii]);
    me200_[ii]->Fill(intLen_[ii]);
    me300_[ii]->Fill(stepLen_[ii]);
#ifdef EDM_ML_DEBUG
    std::string name("Total");
    if (ii < nList_) name = listNames_[ii];
    std::cout << "HGCalTBMB::Volume[" << ii << "]: " << name  << " == Step "
          << stepLen_[ii] << " RadL " << radLen_[ii] << " IntL "
          << intLen_[ii] << std::endl;
#endif
  }
}

Member Data Documentation

std::vector<double> HGCalTBMB::intLen_

private

Definition at line 43 of file HGCalTBMB.h.

Referenced by update().

std::vector<std::string> HGCalTBMB::listNames_

private

Definition at line 39 of file HGCalTBMB.h.

Referenced by findVolume(), HGCalTBMB(), and update().

std::vector<TH1D*> HGCalTBMB::me100_

private

Definition at line 44 of file HGCalTBMB.h.

Referenced by HGCalTBMB(), and update().

std::vector<TH1D*> HGCalTBMB::me200_

private

Definition at line 44 of file HGCalTBMB.h.

Referenced by HGCalTBMB(), and update().

std::vector<TH1D*> HGCalTBMB::me300_

private

Definition at line 44 of file HGCalTBMB.h.

Referenced by HGCalTBMB(), and update().

unsigned int HGCalTBMB::nList_

private

Definition at line 42 of file HGCalTBMB.h.

Referenced by findVolume(), HGCalTBMB(), and update().

std::vector<double> HGCalTBMB::radLen_

private

Definition at line 43 of file HGCalTBMB.h.

Referenced by update().

std::vector<double> HGCalTBMB::stepLen_

private

Definition at line 43 of file HGCalTBMB.h.

Referenced by update().

std::string HGCalTBMB::stopName_

private

Definition at line 40 of file HGCalTBMB.h.

Referenced by findVolume(), and HGCalTBMB().

double HGCalTBMB::stopZ_

private

Definition at line 41 of file HGCalTBMB.h.

Referenced by HGCalTBMB(), and stopAfter().