In This Package:

G4ToHistoryUtensils.cc

Go to the documentation of this file.

#include "G4ToHistoryUtensils.h"
#include "G4Track.hh"
#include "G4ParticleDefinition.hh"
#include "G4VProcess.hh"

using namespace DayaBay;


void SimProcessFromG4Process(const G4VProcess* g4proc, 
                                           SimProcess& outProc)
{
  if(g4proc){
    switch(g4proc->GetProcessType()) {
      case fTransportation      : outProc = SimProcess(SimProcess::kTransportation      ,g4proc->GetProcessName()); return;
      case fElectromagnetic     : outProc = SimProcess(SimProcess::kElectromagnetic     ,g4proc->GetProcessName()); return;
      case fOptical             : outProc = SimProcess(SimProcess::kOptical             ,g4proc->GetProcessName()); return;
      case fHadronic            : outProc = SimProcess(SimProcess::kHadronic            ,g4proc->GetProcessName()); return;
      case fPhotolepton_hadron  : outProc = SimProcess(SimProcess::kPhotolepton_hadron  ,g4proc->GetProcessName()); return;
      case fDecay               : outProc = SimProcess(SimProcess::kDecay               ,g4proc->GetProcessName()); return;
      case fGeneral             : outProc = SimProcess(SimProcess::kGeneral             ,g4proc->GetProcessName()); return;
      case fParameterisation    : outProc = SimProcess(SimProcess::kParameterisation    ,g4proc->GetProcessName()); return;
      case fUserDefined         : outProc = SimProcess(SimProcess::kUserDefined         ,g4proc->GetProcessName()); return;
      default:
      outProc = SimProcess(SimProcess::kUnknown,g4proc->GetProcessName());
    }
  }
}


int pdgFromG4Track(const G4Track* g4track)
{
  int code = g4track->GetDefinition()->GetPDGEncoding();
  // if(code==0) {
  //   if(g4track->GetDefinition()->GetParticleName()=="opticalphoton") code =20022;
  // }
  return code;
}

---