MuonG4Numbering.cc

Go to the documentation of this file.

#include "SimG4CMS/Muon/interface/MuonG4Numbering.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "Geometry/MuonNumbering/interface/MuonGeometryConstants.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "G4VPhysicalVolume.hh"
#include "G4VTouchable.hh"
#include "G4Step.hh"
 
#include <iostream>
 
//#define EDM_ML_DEBUG
 
MuonG4Numbering::MuonG4Numbering(const MuonGeometryConstants& muonConstants) {
  theLevelPart = muonConstants.getValue("level");
  theSuperPart = muonConstants.getValue("super");
  theBasePart = muonConstants.getValue("base");
  theStartCopyNo = muonConstants.getValue("xml_starts_with_copyno");
 
  // some consistency checks
 
  if (theBasePart != 1) {
    edm::LogVerbatim("MuonSim") << "MuonGeometryNumbering finds unusual base constant:" << theBasePart;
  }
  if (theSuperPart < 100) {
    edm::LogVerbatim("MuonSim") << "MuonGeometryNumbering finds unusual super constant:" << theSuperPart;
  }
  if (theLevelPart < 10 * theSuperPart) {
    edm::LogVerbatim("MuonSim") << "MuonGeometryNumbering finds unusual level constant:" << theLevelPart;
  }
  if ((theStartCopyNo != 0) && (theStartCopyNo != 1)) {
    std::cout << "MuonGeometryNumbering finds unusual start value for copy numbers:" << theStartCopyNo << std::endl;
  }
 
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("MuonSim") << "StartCopyNo = " << theStartCopyNo;
  edm::LogVerbatim("MuonSim") << "MuonG4Numbering configured with"
                              << "Level = " << theLevelPart << " Super = " << theSuperPart << " Base = " << theBasePart
                              << " StartCopyNo = " << theStartCopyNo;
#endif
}
 
MuonBaseNumber MuonG4Numbering::PhysicalVolumeToBaseNumber(const G4Step* aStep) {
  MuonBaseNumber num;
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
 
  for (int ii = 0; ii < touch->GetHistoryDepth(); ii++) {
    G4VPhysicalVolume* vol = touch->GetVolume(ii);
    int copyno = vol->GetCopyNo();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("MuonSim") << "MuonG4Numbering: " << vol->GetName() << " " << copyno << std::endl
                                << "Split " << copyNoRelevant(copyno) << ":" << theLevelPart << ":" << theSuperPart
                                << " ";
#endif
    if (copyNoRelevant(copyno)) {
      num.addBase(getCopyNoLevel(copyno), getCopyNoSuperNo(copyno), getCopyNoBaseNo(copyno) - theStartCopyNo);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("MuonSim") << " NoLevel " << getCopyNoLevel(copyno) << " Super " << getCopyNoSuperNo(copyno)
                                  << " Base " << getCopyNoBaseNo(copyno) << " Start " << theStartCopyNo;
#endif
    }
  }
 
  return num;
}
 
const int MuonG4Numbering::getCopyNoLevel(const int copyno) { return copyno / theLevelPart; }
 
const int MuonG4Numbering::getCopyNoSuperNo(const int copyno) { return (copyno % theLevelPart) / theSuperPart; }
 
const int MuonG4Numbering::getCopyNoBaseNo(const int copyno) { return copyno % theSuperPart; }
 
const bool MuonG4Numbering::copyNoRelevant(const int copyno) { return (copyno / theLevelPart) > 0; }