Geant4DetectorSensitivesConstruction.cpp

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation 
//--------------------------------------------------------------------------
// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
// All rights reserved.
//
// For the licensing terms see $DD4hepINSTALL/LICENSE.
// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
//
//  \author Markus Frank
//  \date   2015-11-09
//
//==========================================================================
 
// Framework include files
#include <DDG4/Geant4DetectorConstruction.h>
 
namespace dd4hep {
 
  namespace sim {
 
 
    class Geant4DetectorSensitivesConstruction : public Geant4DetectorConstruction   {
    public:
      Geant4DetectorSensitivesConstruction(Geant4Context* ctxt, const std::string& nam);
      virtual ~Geant4DetectorSensitivesConstruction();
      void constructSensitives(Geant4DetectorConstructionContext* ctxt);
    };
  }    // End namespace sim
}      // End namespace dd4hep
 
 
// Framework include files
#include <DD4hep/InstanceCount.h>
#include <DD4hep/Printout.h>
#include <DD4hep/Plugins.h>
#include <DD4hep/Detector.h>
 
#include <DDG4/Geant4Mapping.h>
#include <DDG4/Geant4Kernel.h>
#include <DDG4/Factories.h>
 
// ROOT include files
#include <TGeoManager.h>
// Geant4 include files
#include <G4SDManager.hh>
#include <G4PVPlacement.hh>
#include <G4VSensitiveDetector.hh>
 
using namespace dd4hep::sim;
 
DECLARE_GEANT4ACTION(Geant4DetectorSensitivesConstruction)
 
Geant4DetectorSensitivesConstruction::Geant4DetectorSensitivesConstruction(Geant4Context* ctxt, const std::string& nam)
: Geant4DetectorConstruction(ctxt,nam)
{
  InstanceCount::increment(this);
}
 
Geant4DetectorSensitivesConstruction::~Geant4DetectorSensitivesConstruction() {
  InstanceCount::decrement(this);
}
 
void Geant4DetectorSensitivesConstruction::constructSensitives(Geant4DetectorConstructionContext* ctxt)   {
  Geant4GeometryInfo* p = Geant4Mapping::instance().ptr();
  const Geant4Kernel& kernel = context()->kernel();
  const auto&         types  = kernel.sensitiveDetectorTypes();
  const std::string&  dflt   = kernel.defaultSensitiveDetectorType();
  for( const auto& iv : p->sensitives )  {
    SensitiveDetector sd = iv.first;
    std::string nam  = sd.name();
    auto        iter = types.find(nam);
    std::string typ  = (iter != types.end()) ? (*iter).second : dflt;
    G4VSensitiveDetector* g4sd = this->createSensitiveDetector(typ, nam);
    for( const TGeoVolume* vol : iv.second )  {
      G4LogicalVolume* g4v = p->g4Volumes[vol];
      if( !g4v )  {
        except("ConstructSDandField: Failed to access G4LogicalVolume for SD %s of type %s.",
               nam.c_str(), typ.c_str());
      }
      ctxt->setSensitiveDetector(g4v, g4sd);
    }
  }
  print("+++ Handled %ld sensitive detectors.",p->sensitives.size());
}