/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/SimG4CMS/Muon/src/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/MuonDDDConstants.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "G4VPhysicalVolume.hh"
#include "G4VTouchable.hh"
#include "G4Step.hh"

#include <iostream>

MuonG4Numbering::MuonG4Numbering(const DDCompactView& cpv){
  MuonDDDConstants muonConstants(cpv);
  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) {
    std::cout << "MuonDDDNumbering finds unusual base constant:"
         <<theBasePart<<std::endl;
  }
  if (theSuperPart<100) {
    std::cout << "MuonDDDNumbering finds unusual super constant:"
         <<theSuperPart<<std::endl;
  }
  if (theLevelPart<10*theSuperPart) {
    std::cout << "MuonDDDNumbering finds unusual level constant:"
         <<theLevelPart<<std::endl;
  }
  if ((theStartCopyNo!=0)&&(theStartCopyNo!=1)) {
    std::cout << "MuonDDDNumbering finds unusual start value for copy numbers:"
         <<theStartCopyNo<<std::endl;
  }

  LogDebug("MuonSimDebug") << "MuonG4Numbering configured with"<<std::endl;
  LogDebug("MuonSimDebug") << "Level = "<<theLevelPart<<" ";
  LogDebug("MuonSimDebug") << "Super = "<<theSuperPart<<" ";
  LogDebug("MuonSimDebug") << "Base = "<<theBasePart<<" ";
  LogDebug("MuonSimDebug") << "StartCopyNo = "<<theStartCopyNo<<std::endl;

}

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();
    LogDebug("MuonSimDebug") << "MuonG4Numbering: " << vol->GetName()<<" "<<copyno<<std::endl;
    if (copyNoRelevant(copyno)) {
      num.addBase(getCopyNoLevel(copyno),
                  getCopyNoSuperNo(copyno),
                  getCopyNoBaseNo(copyno)-theStartCopyNo);
    }
  }

  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;
}