PrintMTDSens Class Reference

Inheritance diagram for PrintMTDSens:

Public Member Functions
	PrintMTDSens (edm::ParameterSet const &p)
	~PrintMTDSens () 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 BeginOfRun *>
	Observer ()
void	slotForUpdate (const BeginOfRun * iT)
virtual	~Observer ()

Private Member Functions
int	dumpTouch (G4VPhysicalVolume *pv, unsigned int leafDepth, int ns)
bool	getBaseNumber ()
G4VPhysicalVolume *	getTopPV ()
void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
BTLNumberingScheme	btlNS_
ETLNumberingScheme	etlNS_
G4NavigationHistory	fHistory
std::vector< std::string >	names_
std::vector< size_t >	nsens_
MTDBaseNumber	thisN_

Additional Inherited Members
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 28 of file PrintMTDSens.cc.

Constructor & Destructor Documentation

PrintMTDSens()

PrintMTDSens::PrintMTDSens

(

edm::ParameterSet const &

p

)

Definition at line 49 of file PrintMTDSens.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::G4cout, names_, nsens_, and AlCaHLTBitMon_ParallelJobs::p.

                                                   : thisN_(), btlNS_(), etlNS_() {
  names_ = p.getUntrackedParameter<std::vector<std::string>>("Name");
  nsens_.resize(names_.size());
  for (size_t index = 0; index < nsens_.size(); index++) {
    nsens_[index] = 0;
  }
  G4cout << "PrintMTDSens:: Print position of MTD Sensitive Touchables: " << G4endl;
  for (const auto &thisName : names_) {
    G4cout << " for name " << thisName << "\n";
  }
  G4cout << " Total of " << names_.size() << " sensitive volume types" << G4endl;
}

~PrintMTDSens()

PrintMTDSens::~PrintMTDSens ( )

override

Definition at line 62 of file PrintMTDSens.cc.

{}

Member Function Documentation

dumpTouch()

int PrintMTDSens::dumpTouch ( G4VPhysicalVolume *  pv, unsigned int  leafDepth, int  ns  )

private

Definition at line 75 of file PrintMTDSens.cc.

References btlNS_, edm::contains(), etlNS_, fHistory, alignBH_cfg::fixed, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, getBaseNumber(), ETLNumberingScheme::getUnitID(), BTLNumberingScheme::getUnitID(), recoMuon::in, PixelPluginsPhase0_cfi::isBarrel, names_, nsens_, groupFilesInBlocks::ntotal, ALPAKA_ACCELERATOR_NAMESPACE::riemannFit::printIt(), DetId::rawId(), cms_rounding::roundIfNear0(), contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, and thisN_.

Referenced by update().

                                                                                 {
  if (leafDepth == 0)
    fHistory.SetFirstEntry(pv);
  else
    fHistory.NewLevel(pv, kNormal, pv->GetCopyNo());

  int ntotal(ns);
  G4ThreeVector globalpoint = fHistory.GetTopTransform().Inverse().TransformPoint(G4ThreeVector(0, 0, 0));
  G4LogicalVolume *lv = pv->GetLogicalVolume();

  std::string mother = "World";
  bool printIt(false);
  if (pv->GetMotherLogical()) {
    mother = pv->GetMotherLogical()->GetName();
  }
  size_t index(0);
  for (const auto &thisName : names_) {
    G4String g4name(thisName);
    if (G4StrUtil::contains(lv->GetName(), g4name)) {
      printIt = true;
      break;
    }
    index++;
  }

  if (lv->GetSensitiveDetector() && printIt) {
    std::stringstream sunitt;
    ++ntotal;
    ++nsens_[index];
    bool isBarrel = getBaseNumber();
    if (isBarrel) {
      BTLDetId theId(btlNS_.getUnitID(thisN_));
      sunitt << theId.rawId();
#ifdef EDM_ML_DEBUG
      G4cout << theId << G4endl;
#endif
    } else {
      ETLDetId theId(etlNS_.getUnitID(thisN_));
      sunitt << theId.rawId();
#ifdef EDM_ML_DEBUG
      G4cout << theId << G4endl;
#endif
    }

    auto fround = [&](double in) {
      std::stringstream ss;
      ss << std::fixed << std::setw(14) << roundIfNear0(in);
      return ss.str();
    };

    G4Box *thisSens = static_cast<G4Box *>(lv->GetSolid());
    G4ThreeVector cn1Global = fHistory.GetTopTransform().Inverse().TransformPoint(
        G4ThreeVector(thisSens->GetXHalfLength(), thisSens->GetYHalfLength(), thisSens->GetZHalfLength()));

    sunitt << fround(globalpoint.x()) << fround(globalpoint.y()) << fround(globalpoint.z()) << fround(cn1Global.x())
           << fround(cn1Global.y()) << fround(cn1Global.z());
    edm::LogVerbatim("MTDG4sensUnitTest") << sunitt.str();
  }

  int NoDaughters = lv->GetNoDaughters();
  while ((NoDaughters--) > 0) {
    G4VPhysicalVolume *pvD = lv->GetDaughter(NoDaughters);
    if (!pvD->IsReplicated())
      ntotal = dumpTouch(pvD, leafDepth + 1, ntotal);
  }

  if (leafDepth > 0)
    fHistory.BackLevel();
  return ntotal;
}

getBaseNumber()

bool PrintMTDSens::getBaseNumber ( )

private

Definition at line 150 of file PrintMTDSens.cc.

References MTDBaseNumber::addLevel(), edm::contains(), fHistory, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, MTDBaseNumber::getCapacity(), cuy::ii, MTDBaseNumber::reset(), MTDBaseNumber::setSize(), and thisN_.

Referenced by dumpTouch().

                                 {
  bool isBTL(false);
  thisN_.reset();
  thisN_.setSize(fHistory.GetMaxDepth());
  int theSize = fHistory.GetDepth() + 1;
  if (thisN_.getCapacity() < theSize)
    thisN_.setSize(theSize);
  //Get name and copy numbers
  if (theSize > 1) {
#ifdef EDM_ML_DEBUG
    G4cout << "Building MTD basenumber:" << G4endl;
#endif
    for (int ii = theSize; ii-- > 0;) {
      thisN_.addLevel(fHistory.GetVolume(ii)->GetName(), fHistory.GetReplicaNo(ii));
#ifdef EDM_ML_DEBUG
      G4cout << "PrintMTDSens::getBaseNumber(): Adding level " << theSize - 1 - ii << ": "
             << fHistory.GetVolume(ii)->GetName() << "[" << fHistory.GetReplicaNo(ii) << "]" << G4endl;
#endif
      if (!isBTL) {
        isBTL = G4StrUtil::contains(fHistory.GetVolume(ii)->GetName(), "BarrelTimingLayer");
      }
    }
  }
  return isBTL;
}

getTopPV()

G4VPhysicalVolume * PrintMTDSens::getTopPV ( )

private

Definition at line 146 of file PrintMTDSens.cc.

Referenced by update().

                                          {
  return G4TransportationManager::GetTransportationManager()->GetNavigatorForTracking()->GetWorldVolume();
}

update()

void PrintMTDSens::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 64 of file PrintMTDSens.cc.

References dumpTouch(), ecalTB2006H4_GenSimDigiReco_cfg::G4cout, getTopPV(), names_, nsens_, and groupFilesInBlocks::ntotal.

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

                                               {
  G4VPhysicalVolume *theTopPV = getTopPV();
  int ntotal = dumpTouch(theTopPV, 0, 0);
  G4cout << "\nTotal number of sensitive detector volumes for MTD is " << ntotal << G4endl;
  for (size_t index = 0; index < nsens_.size(); index++) {
    G4cout << "Sensitive volume " << names_[index] << " # copies " << nsens_[index] << G4endl;
  }
}

Member Data Documentation

btlNS_

BTLNumberingScheme PrintMTDSens::btlNS_

private

Definition at line 45 of file PrintMTDSens.cc.

Referenced by dumpTouch().

etlNS_

ETLNumberingScheme PrintMTDSens::etlNS_

private

Definition at line 46 of file PrintMTDSens.cc.

Referenced by dumpTouch().

fHistory

G4NavigationHistory PrintMTDSens::fHistory

private

Definition at line 42 of file PrintMTDSens.cc.

Referenced by dumpTouch(), and getBaseNumber().

names_

std::vector<std::string> PrintMTDSens::names_

private

Definition at line 40 of file PrintMTDSens.cc.

Referenced by dumpTouch(), PrintMTDSens(), and update().

nsens_

std::vector<size_t> PrintMTDSens::nsens_

private

Definition at line 41 of file PrintMTDSens.cc.

Referenced by dumpTouch(), PrintMTDSens(), and update().

thisN_

MTDBaseNumber PrintMTDSens::thisN_

private

Definition at line 44 of file PrintMTDSens.cc.

Referenced by dumpTouch(), and getBaseNumber().