---

CheckOverlap.cc

Go to the documentation of this file.

#include "Validation/CheckOverlap/interface/CheckOverlap.h"

#include "SimG4Core/Notification/interface/BeginOfJob.h"
#include "SimG4Core/Notification/interface/BeginOfRun.h"

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"

#include "G4Run.hh"
#include "G4PhysicalVolumeStore.hh"
#include "G4PVPlacement.hh"
#include "G4PVParameterised.hh"
#include "G4LogicalVolume.hh"
#include "G4Material.hh"
#include "G4TransportationManager.hh"

#include <set>

CheckOverlap::CheckOverlap(const edm::ParameterSet &p) : topLV(0) {
  nodeName = p.getUntrackedParameter<std::string>("NodeName", "");
  nPoints  = p.getUntrackedParameter<int>("Resolution", 1000);
  edm::LogInfo("G4cout") << "CheckOverlap:: initialised with Node Name "
                         << " " << nodeName << " and Resolution " 
                         << nPoints;
}
 
CheckOverlap::~CheckOverlap() {}
  
void CheckOverlap::update(const BeginOfRun * run) {
  
  edm::LogInfo("G4cout") << "Node Name " << nodeName;
  if (nodeName != "") {
    const G4PhysicalVolumeStore * pvs = G4PhysicalVolumeStore::GetInstance();
    std::vector<G4VPhysicalVolume *>::const_iterator pvcite;
    int i = 0;
    for (pvcite = pvs->begin(); pvcite != pvs->end(); pvcite++) {
      edm::LogInfo("G4cout") << "Name of node " << (++i) << " : " 
                             << (*pvcite)->GetName();
      if ((*pvcite)->GetName() == (G4String)(nodeName)) {
        topLV = (*pvcite)->GetLogicalVolume();
        break;
      }
    }
  } else {
    G4VPhysicalVolume * theTopPV = getTopPV();
    topLV = theTopPV->GetLogicalVolume();
  }
  if (topLV != 0) edm::LogInfo("G4cout") << "Top LV Name " << topLV->GetName();
  else            edm::LogInfo("G4cout") << "No Top LV Found";
  //---------- Check all PV's
  if (topLV) 
    checkHierarchyLeafPVLV(topLV, 0);
}

void CheckOverlap::checkHierarchyLeafPVLV(G4LogicalVolume * lv, 
                                          uint leafDepth) {

  //----- Get LV daughters from list of PV daughters
  mmlvpv lvpvDaughters;
  std::set< G4LogicalVolume * > lvDaughters;
  int NoDaughters = lv->GetNoDaughters();
  while ((NoDaughters--)>0) {
    G4VPhysicalVolume * pvD = lv->GetDaughter(NoDaughters);
    lvpvDaughters.insert(mmlvpv::value_type(pvD->GetLogicalVolume(), pvD));
    lvDaughters.insert(pvD->GetLogicalVolume());
  }
 
  std::set< G4LogicalVolume * >::const_iterator scite;
  mmlvpv::const_iterator mmcite;

  //----- Check daughters of LV
  for (scite = lvDaughters.begin(); scite != lvDaughters.end(); scite++) {
    std::pair< mmlvpv::iterator, mmlvpv::iterator > mmER = lvpvDaughters.equal_range(*scite);    
    //----- Check daughters PV of this LV
    for (mmcite = mmER.first ; mmcite != mmER.second; mmcite++) 
      checkPV((*mmcite).second, leafDepth+1);
    //----- Check daughters LV
    checkHierarchyLeafPVLV(*scite, leafDepth+1);
  }
}

void CheckOverlap::checkPV(G4VPhysicalVolume * pv, uint leafDepth) {

  //----- PV info
#ifndef G4V7
  std::string mother = "World";
  if (pv->GetMotherLogical()) mother = pv->GetMotherLogical()->GetName();
  if (!pv->IsReplicated()) {
    G4PVPlacement* pvplace = dynamic_cast<G4PVPlacement* >(pv);
    G4bool ok = pvplace->CheckOverlaps(nPoints);
    edm::LogInfo("G4cout") << "Placed PV " << pvplace->GetName() 
                           << " Number " << pvplace->GetCopyNo() 
                           << " in mother " << mother << " at depth " 
                           << leafDepth << " Status " << ok;
    if (ok) {
      if(pv->GetRotation() == 0) {
        edm::LogInfo("G4cout") << "Translation " << pv->GetTranslation()
                               << " and with no rotation";
      } else {
        edm::LogInfo("G4cout") << "Translation " << pv->GetTranslation()
                               << " and with rotation "<< *(pv->GetRotation());
      }
      G4LogicalVolume* lv = pv->GetLogicalVolume();
      dumpLV(lv, "Self");
      if (pv->GetMotherLogical()) {
        lv = pv->GetMotherLogical();
        dumpLV (lv, "Mother");
      }
    }
  } else {
    if (pv->GetParameterisation() != 0) {
      G4PVParameterised* pvparam = dynamic_cast<G4PVParameterised* >(pv);
      G4bool ok = pvparam->CheckOverlaps(nPoints);
      edm::LogInfo("G4cout") << "Parametrized PV " << pvparam->GetName()
                             << " in mother " << mother << " at depth "
                             << leafDepth << " Status " << ok;
    }
  }
#endif
}

G4VPhysicalVolume * CheckOverlap::getTopPV() {
  return G4TransportationManager::GetTransportationManager()
    ->GetNavigatorForTracking()->GetWorldVolume();
}

void CheckOverlap::dumpLV(G4LogicalVolume* lv, std::string str) {
  edm::LogInfo("G4cout") << "Dump of " << str << " Logical Volume " 
                         << lv->GetName() << "  Solid: " 
                         << lv->GetSolid()->GetName() << "  Material: "
                         << lv->GetMaterial()->GetName();
  edm::LogInfo("G4cout") << *(lv->GetSolid());
}

---