PhotoMultiplierSDAction.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     : M.Frank
//
//==========================================================================
 
// Framework include files
#include "PMTHit.h"
#include "DDG4/Geant4SensDetAction.inl"
#include "DDG4/Factories.h"
 
//using namespace CLHEP;
 
namespace npdet {
 
  class PhotoMultiplierSD {
  public:
    using Hit            = PMTHit;
    int    nPhotons      = 0;
    double averageEnergy = 0;
  };
}
 
namespace dd4hep {
 
  namespace sim   {
 
    using namespace npdet;
    
    // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    //               Geant4SensitiveAction<PhotoMultiplierSD>
    // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    template <> void Geant4SensitiveAction<PhotoMultiplierSD>::defineCollections()    {
      m_collectionID = declareReadoutFilteredCollection<PhotoMultiplierSD::Hit>();
    }
 
    template <> bool Geant4SensitiveAction<PhotoMultiplierSD>::process(G4Step* step,G4TouchableHistory* /*hist*/ ) {
      Geant4StepHandler h(step);
      Position prePos    = h.prePos();
      Position postPos   = h.postPos();
      Position direction = postPos - prePos;
      Position position  = mean_direction(prePos,postPos);
      double   hit_len   = direction.R();
 
      // Somehow extract here the physics you want
      PhotoMultiplierSD::Hit* hit = new PhotoMultiplierSD::Hit(h.trkID(), h.trkPdgID(), h.deposit(), h.track->GetGlobalTime());
      Geant4HitData::MonteCarloContrib contrib = Geant4HitData::extractContribution(step);
      hit->cellID        = cellID(step);
      hit->energy        = contrib.deposit;
      hit->position      = position;
      hit->momentum      = 0.5*(h. preMom() + h.postMom());
      hit->length        = hit_len;
      collection(m_collectionID)->add(hit);
      mark(h.track);
      if ( 0 == hit->cellID )  {
        hit->cellID        = volumeID(step);
        except("+++ Invalid CELL ID for hit!");
      }
      printP1("Hit with deposit:%f  Pos:%f %f %f ID=%016X",
            step->GetTotalEnergyDeposit(),position.X(),position.Y(),position.Z(),
            (void*)hit->cellID);
      Geant4TouchableHandler handler(step);
      print("    Geant4 path:%s",handler.path().c_str());
 
      // Do something with my personal data (can be also something more clever ;-):
      m_userData.averageEnergy += contrib.deposit;
      ++m_userData.nPhotons;
      return true;
    }
 
  }
}
 
//--- Factory declaration
namespace dd4hep { namespace sim {
    //typedef Geant4SensitiveAction<PhotoMultiplierSD> PhotoMultiplierSDAction;
    using PhotoMultiplierSDAction = Geant4SensitiveAction<PhotoMultiplierSD>;
  }}
DECLARE_GEANT4SENSITIVE(PhotoMultiplierSDAction)