HGCalTBMB Class Reference

Inheritance diagram for HGCalTBMB:

Public Member Functions
	HGCalTBMB (const edm::ParameterSet &)
	HGCalTBMB (const HGCalTBMB &)=delete
const HGCalTBMB &	operator= (const HGCalTBMB &)=delete
	~HGCalTBMB () override
Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)
bool	isMT () const
const SimWatcher &	operator= (const SimWatcher &)=delete
virtual void	registerConsumes (edm::ConsumesCollector)
	SimWatcher ()
	SimWatcher (const SimWatcher &)=delete
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
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
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 24 of file HGCalTBMB.cc.

Constructor & Destructor Documentation

HGCalTBMB::HGCalTBMB

(

const edm::ParameterSet &

p

)

Definition at line 50 of file HGCalTBMB.cc.

References edm::ParameterSet::getParameter(), gpuVertexFinder::hist, mps_fire::i, edm::Service< T >::isAvailable(), isotrackApplyRegressor::k, listNames_, TFileService::make(), me100_, me200_, me300_, mergeVDriftHistosByStation::name, nList_, stopName_, stopZ_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, and runGCPTkAlMap::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::LogVerbatim("HGCSim") << "HGCalTBMB initialized for " << nList_ << " volumes";
  for (unsigned int k = 0; k < nList_; ++k)
    edm::LogVerbatim("HGCSim") << " [" << k << "] " << listNames_[k];
  edm::LogVerbatim("HGCSim") << "Stop after " << stopZ_ << " or reaching volume " << stopName_;

  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, 100000, 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, 100000, 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, 100000, 0.0, 50000.0);
    hist->Sumw2(true);
    me300_.push_back(hist);
  }
  edm::LogVerbatim("HGCSim") << "HGCalTBMB: Booking user histos done ===";
}

HGCalTBMB::HGCalTBMB ( const HGCalTBMB &  )

delete

HGCalTBMB::~HGCalTBMB ( )

override

Definition at line 88 of file HGCalTBMB.cc.

{}

Member Function Documentation

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

private

Definition at line 164 of file HGCalTBMB.cc.

References cuy::ii, isotrackApplyRegressor::k, testEve_cfg::level, listNames_, mergeVDriftHistosByStation::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 &  )

delete

bool HGCalTBMB::stopAfter ( const G4Step *  aStep )

private

Definition at line 148 of file HGCalTBMB.cc.

References findVolume(), and stopZ_.

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
  edm::LogVerbatim("HGCSim") << " HGCalTBMB::Name " << touch->GetVolume(0)->GetName() << " z " << zz << " Flag" << flag;
#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 90 of file HGCalTBMB.cc.

References GeV, intLen_, nList_, radLen_, and stepLen_.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                              {
  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());
  edm::LogVerbatim("HGCSim") << "MaterialBudgetHcalHistos: Track " << aTrack->GetTrackID() << " Code " << theID
                             << " Energy " << theEnergy / CLHEP::GeV << " GeV; Momentum " << mom;
#endif
}

void HGCalTBMB::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 105 of file HGCalTBMB.cc.

References findVolume(), intLen_, nList_, radLen_, stepLen_, stopAfter(), and HLT_FULL_cff::track.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                          {
  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
  edm::LogVerbatim("HGCSim") << "HGCalTBMB::Step in " << touch->GetVolume(0)->GetLogicalVolume()->GetName() << " Index "
                             << indx << " Step " << step << " RadL " << step / radl << " IntL " << step / intl;
#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 133 of file HGCalTBMB.cc.

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

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                            {
  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];
    edm::LogVerbatim("HGCSim") << "HGCalTBMB::Volume[" << ii << "]: " << name << " == Step " << stepLen_[ii] << " RadL "
                               << radLen_[ii] << " IntL " << intLen_[ii];
#endif
  }
}

Member Data Documentation

std::vector<double> HGCalTBMB::intLen_

private

Definition at line 46 of file HGCalTBMB.cc.

Referenced by update().

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

private

Definition at line 42 of file HGCalTBMB.cc.

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

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

private

Definition at line 47 of file HGCalTBMB.cc.

Referenced by HGCalTBMB(), and update().

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

private

Definition at line 47 of file HGCalTBMB.cc.

Referenced by HGCalTBMB(), and update().

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

private

Definition at line 47 of file HGCalTBMB.cc.

Referenced by HGCalTBMB(), and update().

unsigned int HGCalTBMB::nList_

private

Definition at line 45 of file HGCalTBMB.cc.

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

std::vector<double> HGCalTBMB::radLen_

private

Definition at line 46 of file HGCalTBMB.cc.

Referenced by update().

std::vector<double> HGCalTBMB::stepLen_

private

Definition at line 46 of file HGCalTBMB.cc.

Referenced by update().

std::string HGCalTBMB::stopName_

private

Definition at line 43 of file HGCalTBMB.cc.

Referenced by findVolume(), and HGCalTBMB().

double HGCalTBMB::stopZ_

private

Definition at line 44 of file HGCalTBMB.cc.

Referenced by HGCalTBMB(), and stopAfter().