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)