// Co57SourceGen.cc
// Contact person: Valentina Lozza (valentina.lozza@tu-dresden.de)
// See Co57SourceGen.hh for more details
//———————————————————————//

#include <RAT/Co57SourceGen.hh>
#include <RAT/Detector.hh>
#include <RAT/Log.hh>
#include <RAT/string_utilities.hpp>
#include <RAT/DecayChain_Gen.hh>
#include <RAT/GLG4PosGen.hh>
#include <RAT/GLG4TimeGen.hh>
#include <RAT/GLG4StringUtil.hh>
#include <RAT/PosGen_Fill.hh>
#include <RAT/DB.hh>
#include <RAT/EventInfo.hh>
#include <RAT/GLG4StringUtil.hh>
#include <RAT/Factory.hh>

#include <G4UIcommand.hh>

#include <string>
#include <algorithm>
#include <TMath.h>

#include <globals.hh>
#include <G4Event.hh>
#include <G4VPhysicalVolume.hh>
#include <G4LogicalVolume.hh>
#include <G4TransportationManager.hh>
#include <G4LogicalVolumeStore.hh>
#include <G4PhysicalVolumeStore.hh>
#include <G4UnitsTable.hh>

#include <sstream>

namespace RAT {

  Co57SourceGen::Co57SourceGen() : DecayChain_Gen("Co57Source"),fTimeGenName("poisson"),fPosGenName("fill"),
				   fDecayStr("57CoChain"),fSourceName("Co57Source_sample_spot_logical")
  {

    DecayChain_Gen::SetState(fDecayStr + ":" + fPosGenName + ":" + fTimeGenName);

    // The sampler volume is hard-coded for this source.
    // In order to allow for the point generator we also want to store locally this volume
    // so one can compute local coordinates
    Log::Assert( Detector::FindPhysicalVolume(fSourceName) != NULL,
		 "Co57SourceGen: Did not find physical volume '" + fSourceName + "'. Please load the Co57Source geometry.");
    DecayChain_Gen::SetPosState(fSourceName);

  }

  void Co57SourceGen::SetState(G4String state)
  {

    if (state != "") {
      warn << "Warning: parameter '" << state << "' given for /generator/add co57source has no effect.\n";
    }

  }

  void Co57SourceGen::SetPosState(G4String /*state*/) {
    warn << "Warning: /generator/set/pos has no effect for the co57source generator.\n";
  }


} // namespace RAT