Geant4OutputAction.cpp

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation for LCD
//--------------------------------------------------------------------------
// 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 "DD4hep/Printout.h"
#include "DD4hep/InstanceCount.h"
#include "DDG4/Geant4Particle.h"
#include "DDG4/Geant4RunAction.h"
#include "DDG4/Geant4OutputAction.h"

// Geant 4 includes
#include "G4HCofThisEvent.hh"
#include "G4Event.hh"

using namespace DD4hep::Simulation;
using namespace DD4hep;
using namespace std;

Geant4OutputAction::Geant4OutputAction(Geant4Context* ctxt, const string& nam)
  : Geant4EventAction(ctxt, nam), m_truth(0)
{
  InstanceCount::increment(this);
  declareProperty("Output", m_output);
  declareProperty("HandleErrorsAsFatal", m_errorFatal=true);
  // Need to instantiate run action to configure fibers
  ctxt->runAction();
}

Geant4OutputAction::~Geant4OutputAction() {
  InstanceCount::decrement(this);
}

void Geant4OutputAction::configureFiber(Geant4Context* thread_ctxt)  {
  Geant4EventAction::configureFiber(thread_ctxt);
  thread_ctxt->runAction().callAtBegin(this, &Geant4OutputAction::beginRun);
  thread_ctxt->runAction().callAtEnd(this, &Geant4OutputAction::endRun);
}

void Geant4OutputAction::begin(const G4Event* /* event */) {
}

void Geant4OutputAction::end(const G4Event* evt) {
  OutputContext < G4Event > ctxt(evt);
  G4HCofThisEvent* hce = evt->GetHCofThisEvent();
  if ( hce )  {
    int nCol = hce->GetNumberOfCollections();
    try  {
      m_truth = context()->event().extension<Geant4ParticleMap>(false);
      if ( m_truth && !m_truth->isValid() )  {
        m_truth = 0;
        printout(WARNING,name(),"+++ [Event:%d] No valid MC truth info present. "
                 "Is a Particle handler installed ?",evt->GetEventID());
      }
      try  {
        saveEvent(ctxt);
        for (int i = 0; i < nCol; ++i) {
          G4VHitsCollection* hc = hce->GetHC(i);
          saveCollection(ctxt, hc);
        }
      }
      catch(const exception& e)   {
        printout(ERROR,name(),"+++ [Event:%d] Exception while saving event:%s",
                 evt->GetEventID(),e.what());
        if ( m_errorFatal ) throw;
      }
      catch(...)   {
        printout(ERROR,name(),"+++ [Event:%d] UNKNWON Exception while saving event",
                 evt->GetEventID());
        if ( m_errorFatal ) throw;
      }
      commit(ctxt);
    }
    catch(const exception& e)   {
      printout(ERROR,name(),"+++ [Event:%d] Exception while saving event:%s",
               evt->GetEventID(),e.what());
      if ( m_errorFatal ) throw;
    }
    catch(...)   {
      printout(ERROR,name(),"+++ [Event:%d] UNKNWON Exception while saving event",
               evt->GetEventID());
      if ( m_errorFatal ) throw;
    }
    m_truth = 0;
    return;
  }
  printout(WARNING,"Geant4OutputAction",
           "+++ The value of G4HCofThisEvent is NULL. No collections saved!");
}

void Geant4OutputAction::commit(OutputContext<G4Event>& /* ctxt */) {
}

void Geant4OutputAction::beginRun(const G4Run* /* run */) {
}

void Geant4OutputAction::endRun(const G4Run* /* run */) {
}

void Geant4OutputAction::saveRun(const G4Run* /* run */) {
}

void Geant4OutputAction::saveEvent(OutputContext<G4Event>& /* ctxt */) {
}

void Geant4OutputAction::saveCollection(OutputContext<G4Event>& /* ctxt */, G4VHitsCollection* /* collection */) {
}