Deprecated/UpgradeASCIITests/G4TankSimulatorASCII/G4TankEventAction.cc

Go to the documentation of this file.

#include "G4TankSimulator.h"
#include "G4TankEventAction.h"
#include "G4TankTrackingAction.h"
#include "G4TankSteppingAction.h"
#include "ParticleCases.h"
using namespace G4TankSimulatorASCII;

#include <G4Event.hh>
#include <G4Timer.hh>

#include <fwk/CentralConfig.h>
#include <fwk/RunController.h>
#include <fwk/VModule.h>
#include <fwk/RandomEngineRegistry.h>
#include <CLHEP/Random/RandFlat.h>
#include <CLHEP/Random/RandPoisson.h>
#include <CLHEP/Random/RandExponential.h>
using namespace fwk;

#include <sevt/PMT.h>
#include <sevt/PMTSimData.h>
#include <sevt/SEvent.h>
#include <sevt/Station.h>
#include <sevt/StationSimData.h>
using namespace sevt;

#include <utl/ErrorLogger.h>
#include <utl/Particle.h>
#include <utl/Trace.h>
#include <utl/PhysicalConstants.h>
using namespace utl;

#include <iostream>
#include <sstream>
using namespace std;

using CLHEP::RandFlat;
using CLHEP::RandPoisson;
using CLHEP::RandExponential;


int G4TankEventAction::fgPrimaryId = 0;

G4TankEventAction::G4TankEventAction() :
  fTimer(new G4Timer)
{
  fRandomEngine =
    &RandomEngineRegistry::GetInstance().Get(RandomEngineRegistry::eDetector);
}

G4TankEventAction::~G4TankEventAction()
{
  delete fTimer;
}

void
G4TankEventAction::BeginOfEventAction(const G4Event* Evt)
{
  fgPrimaryId = Evt->GetPrimaryVertex()->GetPrimary()->GetPDGcode();
  fTimer->Start();

  G4TankTrackingAction::WCD_Track_p=0;

  G4TankTrackingAction::ASCII_Track=0;
  G4TankTrackingAction::ASCII_Track_p=0;
  
  G4TankTrackingAction::ASCII_Edep=0.;
  G4TankTrackingAction::ASCII_Edep_p=0.;
  G4TankTrackingAction::ASCII_Edep_Michel=0.;

  for (int i=0;i<27;i++) G4TankTrackingAction::ASCII_Edep_i[i]=0.; 

  G4TankTrackingAction::flagDecay=0;
  G4TankTrackingAction::MichelElectronKE=0.;
  G4TankTrackingAction::MichelElectronTheta=0.;
  G4TankTrackingAction::MichelElectronZ=0.;
  for (int idet=0;idet<2;idet++)
    {
      G4TankTrackingAction::MichelElectronNPE[idet]=0;
      G4TankTrackingAction::MichelElectronNPE_time[idet]=0.;
    }
}


void
G4TankEventAction::EndOfEventAction(const G4Event* event)
{
 
  // Reset WCD counters 
  G4TankTrackingAction::fNumC = 0;
  G4TankTrackingAction::fNumCDelta = 0;
  G4TankTrackingAction::fNumBounces = 0;

   for (int idet=0;idet<2;idet++)
    {
      if ( G4TankTrackingAction::MichelElectronNPE[idet]==0 ) continue;
      G4TankTrackingAction::MichelElectronNPE_time[idet]/=double(G4TankTrackingAction::MichelElectronNPE[idet]);
    }
 
  fTimer->Stop();
  
}