From 3c33ec3b9fc513794e7e683620fbce99aedb9bcb Mon Sep 17 00:00:00 2001
From: erkn <erik@vonbrun.org>
Date: Tue, 9 Aug 2022 08:12:19 +0200
Subject: [PATCH] refactor MCPL file read part of FileSource

A new constructor is added to FileSource that reads from
a given MCPL-file.
---
 src/source.cpp | 67 --------------------------------------------------
 1 file changed, 67 deletions(-)

diff --git a/src/source.cpp b/src/source.cpp
index 99ec3a74d..e43f03929 100644
--- a/src/source.cpp
+++ b/src/source.cpp
@@ -437,73 +437,6 @@ CustomSourceWrapper::~CustomSourceWrapper()
 }
 
 
-MCPLFileSource::~MCPLFileSource(){
-  mcpl_close_file(mcpl_file);
-}
-
-SourceSite MCPLFileSource::sample(uint64_t* seed) const
-{
-  size_t i_site = sites_.size() * prn(seed);
-  return sites_[i_site];
-}
-
-void MCPLFileSource::read_source_bank(vector<SourceSite> &sites_)
-{
-  sites_.resize(n_sites);
-  for (int i=0;i<n_sites;i++){
-    sites_[i]=read_single_particle();
-  }
-}
-
-SourceSite MCPLFileSource::read_single_particle() const
-{
-  SourceSite omc_particle_;
-  const mcpl_particle_t *mcpl_particle;
-  // extract particle from mcpl-file
-  mcpl_particle=mcpl_read(mcpl_file);
-  // check if it is a neutron, photon, electron, or positron. Otherwise skip.
-  int pdg=mcpl_particle->pdgcode;
-  while ( pdg!=2112 && pdg!=22 && pdg!=11 && pdg!=-11) {
-    mcpl_particle=mcpl_read(mcpl_file);
-    pdg=mcpl_particle->pdgcode;
-    //should check for file exhaustion. This could happen if particles are other than
-    //neutrons, photons, electrons, or positrons.
-  }
-
-  switch(pdg){
-    case 2112:
-      omc_particle_.particle=ParticleType::neutron;
-      break;
-    case 22:
-      omc_particle_.particle=ParticleType::photon;
-      break;
-    case 11:
-      omc_particle_.particle=ParticleType::electron;
-      break;
-    case -11:
-      omc_particle_.particle=ParticleType::positron;
-      break;
-  }
-
-  //particle is good, convert to openmc-formalism
-  omc_particle_.r.x=mcpl_particle->position[0];
-  omc_particle_.r.y=mcpl_particle->position[1];
-  omc_particle_.r.z=mcpl_particle->position[2];
-
-  omc_particle_.u.x=mcpl_particle->direction[0];
-  omc_particle_.u.y=mcpl_particle->direction[1];
-  omc_particle_.u.z=mcpl_particle->direction[2];
-
-  //mcpl stores kinetic energy in MeV
-  omc_particle_.E=mcpl_particle->ekin*1e6;
-  //mcpl stores time in ms
-  omc_particle_.time=mcpl_particle->time*1e-3;
-  omc_particle_.wgt=mcpl_particle->weight;
-  omc_particle_.delayed_group=0;
-  return omc_particle_;
-}
-#endif //OPENMC_MCPL
-
 //==============================================================================
 // Non-member functions
 //==============================================================================