MaterialBudget Class Reference

#include <MaterialBudget.h>

Inheritance diagram for MaterialBudget:

Public Member Functions
	MaterialBudget (const edm::ParameterSet &)
	MaterialBudget (const MaterialBudget &)=delete
const MaterialBudget &	operator= (const MaterialBudget &)=delete
	~MaterialBudget () override
Public Member Functions inherited from SimWatcher
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 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
void	book (const edm::ParameterSet &)
bool	stopAfter (const G4Step *)
void	update (const BeginOfRun *) 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 EndOfTrack *) 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...

Private Attributes
std::vector< double >	boundaries
std::vector< int >	constituents
std::vector< int >	detLevels
std::vector< std::string >	detNames
std::vector< std::string >	detTypes
double	eta_
std::vector< double >	etaRegions
std::vector< double >	intLen
std::vector< G4LogicalVolume * >	logVolumes
std::vector< TProfile * >	me100
std::vector< TProfile * >	me200
std::vector< TProfile * >	me300
std::vector< TProfile * >	me400
std::vector< TProfile * >	me500
std::vector< TProfile * >	me600
double	phi_
std::vector< double >	radLen
std::vector< int >	regionTypes
std::vector< int >	stackOrder
std::vector< double >	stepLen
double	stepT

Additional Inherited Members

Detailed Description

Definition at line 24 of file MaterialBudget.h.

Constructor & Destructor Documentation

MaterialBudget() [1/2]

MaterialBudget::MaterialBudget

(

const edm::ParameterSet &

p

)

Definition at line 19 of file MaterialBudget.cc.

                                                       {
  edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("MaterialBudget");
  detTypes = m_p.getParameter<std::vector<std::string> >("DetectorTypes");
  constituents = m_p.getParameter<std::vector<int> >("Constituents");
  stackOrder = m_p.getParameter<std::vector<int> >("StackOrder");
  detNames = m_p.getParameter<std::vector<std::string> >("DetectorNames");
  detLevels = m_p.getParameter<std::vector<int> >("DetectorLevels");
  etaRegions = m_p.getParameter<std::vector<double> >("EtaBoundaries");
  regionTypes = m_p.getParameter<std::vector<int> >("RegionTypes");
  boundaries = m_p.getParameter<std::vector<double> >("Boundaries");
  edm::LogInfo("MaterialBudget") << "MaterialBudget initialized for " << detTypes.size() << " detector types";
  unsigned int nc = 0;
  for (unsigned int ii = 0; ii < detTypes.size(); ++ii) {
    edm::LogInfo("MaterialBudget") << "Type[" << ii << "] : " << detTypes[ii] << " with " << constituents[ii]
                                   << ", order " << stackOrder[ii] << " constituents --> ";
    for (int kk = 0; kk < constituents[ii]; ++kk) {
      std::string name = "Unknown";
      int level = -1;
      if (nc < detNames.size()) {
        name = detNames[nc];
        level = detLevels[nc];
        ++nc;
      }
      edm::LogInfo("MaterialBudget") << "    Constituent[" << kk << "]: " << name << " at level " << level;
    }
  }
  edm::LogInfo("MaterialBudget") << "MaterialBudget Stop condition for " << etaRegions.size() << " eta regions";
  for (unsigned int ii = 0; ii < etaRegions.size(); ++ii) {
    edm::LogInfo("MaterialBudget") << "Region[" << ii << "] : Eta < " << etaRegions[ii] << " boundary type "
                                   << regionTypes[ii] << " limit " << boundaries[ii];
  }
  book(m_p);
}

References book(), boundaries, constituents, detLevels, detNames, detTypes, etaRegions, edm::ParameterSet::getParameter(), cuy::ii, GetRecoTauVFromDQM_MC_cff::kk, personalPlayback::level, Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, regionTypes, stackOrder, and AlCaHLTBitMon_QueryRunRegistry::string.

MaterialBudget() [2/2]

MaterialBudget::MaterialBudget ( const MaterialBudget &  )

delete

~MaterialBudget()

MaterialBudget::~MaterialBudget ( )

override

Definition at line 53 of file MaterialBudget.cc.

{}

Member Function Documentation

book()

void MaterialBudget::book ( const edm::ParameterSet &  m_p )

private

Definition at line 200 of file MaterialBudget.cc.

                                                    {
  // Book histograms
  edm::Service<TFileService> tfile;
 
  if (!tfile.isAvailable())
    throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                      << "please add it to config file";
 
  int binEta = m_p.getUntrackedParameter<int>("NBinEta", 320);
  int binPhi = m_p.getUntrackedParameter<int>("NBinPhi", 180);
  double minEta = m_p.getUntrackedParameter<double>("MinEta", -8.0);
  double maxEta = m_p.getUntrackedParameter<double>("MaxEta", 8.0);
  double maxPhi = CLHEP::pi;
  edm::LogInfo("MaterialBudget") << "MaterialBudget: Booking user "
                                 << "histos === with " << binEta << " bins "
                                 << "in eta from " << minEta << " to " << maxEta << " and " << binPhi << " bins "
                                 << "in phi from " << -maxPhi << " to " << maxPhi;
 
  char name[20], title[80];
  std::string named;
  for (unsigned int i = 0; i <= detTypes.size(); i++) {
    if (i >= detTypes.size())
      named = "Unknown";
    else
      named = detTypes[i];
    sprintf(name, "%d", i + 100);
    sprintf(title, "MB(X0) prof Eta in %s", named.c_str());
    me100.push_back(tfile->make<TProfile>(name, title, binEta, minEta, maxEta));
    sprintf(name, "%d", i + 200);
    sprintf(title, "MB(L0) prof Eta in %s", named.c_str());
    me200.push_back(tfile->make<TProfile>(name, title, binEta, minEta, maxEta));
    sprintf(name, "%d", i + 300);
    sprintf(title, "MB(Step) prof Eta in %s", named.c_str());
    me300.push_back(tfile->make<TProfile>(name, title, binEta, minEta, maxEta));
    sprintf(name, "%d", i + 400);
    sprintf(title, "MB(X0) prof Phi in %s", named.c_str());
    me400.push_back(tfile->make<TProfile>(name, title, binPhi, -maxPhi, maxPhi));
    sprintf(name, "%d", i + 500);
    sprintf(title, "MB(L0) prof Phi in %s", named.c_str());
    me500.push_back(tfile->make<TProfile>(name, title, binPhi, -maxPhi, maxPhi));
    sprintf(name, "%d", i + 600);
    sprintf(title, "MB(Step) prof Phi in %s", named.c_str());
    me600.push_back(tfile->make<TProfile>(name, title, binPhi, -maxPhi, maxPhi));
  }
 
  edm::LogInfo("MaterialBudget") << "MaterialBudget: Booking user "
                                 << "histos done ===";
}

References detTypes, edm::ParameterSet::getUntrackedParameter(), mps_fire::i, maxEta, HLT_FULL_cff::maxPhi, me100, me200, me300, me400, me500, me600, EgHLTOffEleSelection_cfi::minEta, Skims_PA_cff::name, pi, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, and runGCPTkAlMap::title.

Referenced by MaterialBudget().

operator=()

const MaterialBudget& MaterialBudget::operator= ( const MaterialBudget &  )

delete

stopAfter()

bool MaterialBudget::stopAfter ( const G4Step *  aStep )

private

Definition at line 249 of file MaterialBudget.cc.

                                                  {
  G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
  if (aStep->GetPostStepPoint() != nullptr)
    hitPoint = aStep->GetPostStepPoint()->GetPosition();
  double rr = hitPoint.perp();
  double zz = std::abs(hitPoint.z());
 
  bool flag(false);
  for (unsigned int ii = 0; ii < etaRegions.size(); ++ii) {
#ifdef EDM_ML_DEBUG
    edm::LogInfo("MaterialBudget") << " MaterialBudget::Eta " << eta_ << " in Region[" << ii << "] with "
                                   << etaRegions[ii] << " type " << regionTypes[ii] << "|" << boundaries[ii];
#endif
    if (fabs(eta_) < etaRegions[ii]) {
      if (regionTypes[ii] == 0) {
        if (rr >= boundaries[ii] - 0.001)
          flag = true;
      } else {
        if (zz >= boundaries[ii] - 0.001)
          flag = true;
      }
#ifdef EDM_ML_DEBUG
      if (flag)
        edm::LogInfo("MaterialBudget") << " MaterialBudget::Stop after R = " << rr << " and Z = " << zz;
#endif
      break;
    }
  }
#ifdef EDM_ML_DEBUG
  edm::LogInfo("MaterialBudget") << " MaterialBudget:: R = " << rr << " and Z = " << zz << " Flag " << flag;
#endif
  return flag;
}

References funct::abs(), boundaries, eta_, etaRegions, RemoveAddSevLevel::flag, cuy::ii, regionTypes, findQualityFiles::rr, and geometryCSVtoXML::zz.

Referenced by update().

update() [1/4]

void MaterialBudget::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 55 of file MaterialBudget.cc.

                                             {
  const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
  std::vector<G4LogicalVolume*>::const_iterator lvcite;
 
  unsigned int kount = detNames.size();
  for (unsigned int ii = 0; ii < kount; ++ii)
    logVolumes.push_back(nullptr);
 
  for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
    for (unsigned int ii = 0; ii < detNames.size(); ++ii) {
      if (strcmp(detNames[ii].c_str(),
                 static_cast<std::string>(dd4hep::dd::noNamespace((*lvcite)->GetName())).c_str()) == 0) {
        logVolumes[ii] = (*lvcite);
        kount--;
        break;
      }
    }
    if (kount <= 0)
      break;
  }
  edm::LogInfo("MaterialBudget") << "MaterialBudget::Finds " << detNames.size() - kount << " out of " << detNames.size()
                                 << " LV addresses";
  for (unsigned int ii = 0; ii < detNames.size(); ++ii) {
    std::string name("Unknown");
    if (logVolumes[ii] != nullptr)
      name = dd4hep::dd::noNamespace(logVolumes[ii]->GetName());
    edm::LogInfo("MaterialBudget") << "LV[" << ii << "] : " << detNames[ii] << " Address " << logVolumes[ii] << " | "
                                   << name;
  }
 
  for (unsigned int ii = 0; ii < (detTypes.size() + 1); ++ii) {
    stepLen.push_back(0);
    radLen.push_back(0);
    intLen.push_back(0);
  }
  stackOrder.push_back(0);
}

References detNames, detTypes, cuy::ii, intLen, logVolumes, Skims_PA_cff::name, radLen, stackOrder, stepLen, and AlCaHLTBitMon_QueryRunRegistry::string.

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

update() [2/4]

void MaterialBudget::update ( const BeginOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 93 of file MaterialBudget.cc.

                                                   {
  for (unsigned int ii = 0; ii < (detTypes.size() + 1); ++ii) {
    stepLen[ii] = 0;
    radLen[ii] = 0;
    intLen[ii] = 0;
  }
 
  const G4Track* aTrack = (*trk)();  // recover G4 pointer if wanted
  const G4ThreeVector& mom = aTrack->GetMomentum();
  if (mom.theta() != 0) {
    eta_ = mom.eta();
  } else {
    eta_ = -99;
  }
  phi_ = mom.phi();
  stepT = 0;
 
#ifdef EDM_ML_DEBUG
  double theEnergy = aTrack->GetTotalEnergy();
  int theID = (int)(aTrack->GetDefinition()->GetPDGEncoding());
  edm::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: Track " << aTrack->GetTrackID() << " Code " << theID
                                 << " Energy " << theEnergy / CLHEP::GeV << " GeV; Eta " << eta_ << " Phi "
                                 << phi_ / CLHEP::deg << " PT " << mom.perp() / CLHEP::GeV << " GeV *****";
#endif
}

References detTypes, eta_, GeV, cuy::ii, createfilelist::int, intLen, phi_, radLen, stepLen, and stepT.

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

update() [3/4]

void MaterialBudget::update ( const EndOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 161 of file MaterialBudget.cc.

                                                 {
  for (unsigned int ii = 0; ii < detTypes.size(); ++ii) {
    for (unsigned int jj = 0; jj <= detTypes.size(); ++jj) {
      if (stackOrder[jj] == (int)(ii + 1)) {
        for (unsigned int kk = 0; kk <= detTypes.size(); ++kk) {
          if (stackOrder[kk] == (int)(ii)) {
            radLen[jj] += radLen[kk];
            intLen[jj] += intLen[kk];
#ifdef EDM_ML_DEBUG
            edm::LogInfo("MaterialBudget")
                << "MaterialBudget::Add " << kk << ":" << stackOrder[kk] << " to " << jj << ":" << stackOrder[jj]
                << " RadL " << radLen[kk] << " : " << radLen[jj] << " IntL " << intLen[kk] << " : " << intLen[jj]
                << " StepL " << stepLen[kk] << " : " << stepLen[jj];
#endif
            break;
          }
        }
        break;
      }
    }
  }
 
  for (unsigned int ii = 0; ii <= detTypes.size(); ++ii) {
    me100[ii]->Fill(eta_, radLen[ii]);
    me200[ii]->Fill(eta_, intLen[ii]);
    me300[ii]->Fill(eta_, stepLen[ii]);
    me400[ii]->Fill(phi_, radLen[ii]);
    me500[ii]->Fill(phi_, intLen[ii]);
    me600[ii]->Fill(phi_, stepLen[ii]);
#ifdef EDM_ML_DEBUG
    std::string name("Unknown");
    if (ii < detTypes.size())
      name = detTypes[ii];
    edm::LogInfo("MaterialBudget") << "MaterialBudget::Volume[" << ii << "]: " << name << " eta " << eta_ << " == Step "
                                   << stepLen[ii] << " RadL " << radLen[ii] << " IntL " << intLen[ii];
#endif
  }
}

References detTypes, eta_, cuy::ii, intLen, findQualityFiles::jj, GetRecoTauVFromDQM_MC_cff::kk, me100, me200, me300, me400, me500, me600, Skims_PA_cff::name, phi_, radLen, stackOrder, stepLen, and AlCaHLTBitMon_QueryRunRegistry::string.

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

update() [4/4]

void MaterialBudget::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 119 of file MaterialBudget.cc.

                                               {
  //---------- each step
  G4Material* material = aStep->GetPreStepPoint()->GetMaterial();
  double step = aStep->GetStepLength();
  double radl = material->GetRadlen();
  double intl = material->GetNuclearInterLength();
 
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  unsigned int indx = detTypes.size();
  unsigned int nc = 0;
  for (unsigned int ii = 0; ii < detTypes.size(); ++ii) {
    for (int kk = 0; kk < constituents[ii]; ++kk) {
      if (detLevels[nc + kk] <= (int)((touch->GetHistoryDepth()) + 1)) {
        int jj = (int)((touch->GetHistoryDepth()) + 1) - detLevels[nc + kk];
        if ((touch->GetVolume(jj)->GetLogicalVolume()) == logVolumes[nc + kk]) {
          indx = ii;
          break;
        }
      }
    }
    nc += (unsigned int)(constituents[ii]);
    if (indx == ii)
      break;
  }
 
  stepT += step;
  stepLen[indx] = stepT;
  radLen[indx] += (step / radl);
  intLen[indx] += (step / intl);
#ifdef EDM_ML_DEBUG
  edm::LogInfo("MaterialBudget") << "MaterialBudget::Step in "
                                 << dd4hep::dd::noNamespace(touch->GetVolume(0)->GetLogicalVolume()->GetName())
                                 << " Index " << indx << " Step " << step << " RadL " << step / radl << " IntL "
                                 << step / intl;
#endif
  //----- Stop tracking after selected position
  if (stopAfter(aStep)) {
    G4Track* track = aStep->GetTrack();
    track->SetTrackStatus(fStopAndKill);
  }
}

References constituents, detLevels, detTypes, cuy::ii, createfilelist::int, intLen, findQualityFiles::jj, GetRecoTauVFromDQM_MC_cff::kk, logVolumes, radLen, stepLen, stepT, stopAfter(), and HLT_FULL_cff::track.

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

Member Data Documentation

boundaries

std::vector<double> MaterialBudget::boundaries

private

Definition at line 47 of file MaterialBudget.h.

Referenced by MaterialBudget(), and stopAfter().

constituents

std::vector<int> MaterialBudget::constituents

private

Definition at line 46 of file MaterialBudget.h.

Referenced by MaterialBudget(), and update().

detLevels

std::vector<int> MaterialBudget::detLevels

private

Definition at line 46 of file MaterialBudget.h.

Referenced by MaterialBudget(), and update().

detNames

std::vector<std::string> MaterialBudget::detNames

private

Definition at line 45 of file MaterialBudget.h.

Referenced by MaterialBudget(), and update().

detTypes

std::vector<std::string> MaterialBudget::detTypes

private

Definition at line 45 of file MaterialBudget.h.

Referenced by book(), MaterialBudget(), and update().

eta_

double MaterialBudget::eta_

private

Definition at line 51 of file MaterialBudget.h.

Referenced by stopAfter(), and update().

etaRegions

std::vector<double> MaterialBudget::etaRegions

private

Definition at line 47 of file MaterialBudget.h.

Referenced by MaterialBudget(), and stopAfter().

intLen

std::vector<double> MaterialBudget::intLen

private

Definition at line 50 of file MaterialBudget.h.

Referenced by update().

logVolumes

std::vector<G4LogicalVolume*> MaterialBudget::logVolumes

private

Definition at line 48 of file MaterialBudget.h.

Referenced by update().

me100

std::vector<TProfile*> MaterialBudget::me100

private

Definition at line 49 of file MaterialBudget.h.

Referenced by book(), and update().

me200

std::vector<TProfile*> MaterialBudget::me200

private

Definition at line 49 of file MaterialBudget.h.

Referenced by book(), and update().

me300

std::vector<TProfile*> MaterialBudget::me300

private

Definition at line 49 of file MaterialBudget.h.

Referenced by book(), and update().

me400

std::vector<TProfile*> MaterialBudget::me400

private

Definition at line 49 of file MaterialBudget.h.

Referenced by book(), and update().

me500

std::vector<TProfile*> MaterialBudget::me500

private

Definition at line 49 of file MaterialBudget.h.

Referenced by book(), and update().

me600

std::vector<TProfile*> MaterialBudget::me600

private

Definition at line 49 of file MaterialBudget.h.

Referenced by book(), and update().

phi_

double MaterialBudget::phi_

private

Definition at line 51 of file MaterialBudget.h.

Referenced by update().

radLen

std::vector<double> MaterialBudget::radLen

private

Definition at line 50 of file MaterialBudget.h.

Referenced by update().

regionTypes

std::vector<int> MaterialBudget::regionTypes

private

Definition at line 46 of file MaterialBudget.h.

Referenced by MaterialBudget(), and stopAfter().

stackOrder

std::vector<int> MaterialBudget::stackOrder

private

Definition at line 46 of file MaterialBudget.h.

Referenced by MaterialBudget(), and update().

stepLen

std::vector<double> MaterialBudget::stepLen

private

Definition at line 50 of file MaterialBudget.h.

Referenced by update().

stepT

double MaterialBudget::stepT

private

Definition at line 51 of file MaterialBudget.h.

Referenced by update().