//
/// \file PhysicsList.cc
/// \brief Implementation of the PhysicsList class
//

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

#include "PhysicsList.hh"

#include "G4SystemOfUnits.hh"
#include "G4UnitsTable.hh"
#include "G4NuclideTable.hh"
#include "G4LossTableManager.hh"
#include "G4UAtomicDeexcitation.hh"
#include "G4NuclearLevelData.hh"
#include "G4DeexPrecoParameters.hh"
#include "G4PhysListUtil.hh"
#include "G4EmParameters.hh"
#include "G4GenericIon.hh"
#include "G4Radioactivation.hh"
#include "G4PhysicsListHelper.hh"
#include "G4EmBuilder.hh"         // For G4EmBuilder::ConstructMinimalEmSet()
#include "G4IonConstructor.hh"    // For G4IonConstructor

// Hadronic + EM Physics
#include "HadronElasticPhysicsHP.hh"
#include "G4HadronPhysicsFTFP_BERT_HP.hh"
#include "G4IonPhysicsXS.hh"
#include "G4StoppingPhysics.hh"
#include "GammaNuclearPhysics.hh"
#include "G4EmStandardPhysics_option4.hh"
#include "G4IonElasticPhysics.hh"
#include "RadioactiveDecayPhysics.hh"
#include "G4RadioactiveDecayPhysics.hh"
// Decay physics
#include "G4DecayPhysics.hh"

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

PhysicsList::PhysicsList()
:G4VModularPhysicsList()
{
  G4int verb = 1;
  SetVerboseLevel(verb);
  
  //add new units
  //
  new G4UnitDefinition( "mm2/g",  "mm2/g", "Surface/Mass", mm2/g);
  new G4UnitDefinition( "um2/mg", "um2/mg","Surface/Mass", um*um/mg);
  
  // mandatory for G4NuclideTable
  //
  const G4double meanLife = 1*picosecond, halfLife = meanLife*std::log(2);
  G4NuclideTable::GetInstance()->SetThresholdOfHalfLife(halfLife);
  G4NuclideTable::GetInstance()->SetMeanLifeThreshold(meanLife);
  G4NuclideTable::GetInstance()->SetLevelTolerance(1.0 * eV);

  
  // Hadron Elastic scattering
  RegisterPhysics( new HadronElasticPhysicsHP(verb) );
  ////RegisterPhysics( new G4HadronElasticPhysicsXS(verb) );
    
  // Hadron Inelastic Physics
  RegisterPhysics( new G4HadronPhysicsFTFP_BERT_HP(verb));
  //RegisterPhysics( new G4HadronPhysicsQGSP_BIC(verb));
  ////RegisterPhysics( new G4HadronPhysicsQGSP_BIC_AllHP(verb));
  ////RegisterPhysics( new G4HadronInelasticQBBC(verb));
  ////RegisterPhysics( new G4HadronPhysicsINCLXX(verb));
  
  // Ion Elastic scattering
  //
  RegisterPhysics( new G4IonElasticPhysics(verb));
  
  // Ion Inelastic physics
  RegisterPhysics( new G4IonPhysicsXS(verb));
  ////RegisterPhysics( new G4IonINCLXXPhysics(verb));
  
  // stopping Particles
  RegisterPhysics( new G4StoppingPhysics(verb));
      
  // Gamma-Nuclear Physics
  RegisterPhysics( new GammaNuclearPhysics("gamma"));
  ////RegisterPhysics( new GammaNuclearPhysicsLEND("gamma"));
    
  // EM physics
 // RegisterPhysics(new ElectromagneticPhysics());
  RegisterPhysics(new G4EmStandardPhysics_option4());
  
  // Decay
  RegisterPhysics(new G4DecayPhysics());

  // Radioactive decay
  RegisterPhysics(new RadioactiveDecayPhysics());
   // RegisterPhysics(new G4RadioactiveDecayPhysics());  

  G4EmParameters* emParams = G4EmParameters::Instance();
  emParams->SetDefaults();
  emParams->SetFluo(true);
  emParams->SetAuger(true);
  emParams->SetAugerCascade(true);
  emParams->SetPixe(false); // Optional
  emParams->SetDeexcitationIgnoreCut(true);
  emParams->SetApplyCuts(true);

  // Setup nuclear deexcitation parameters
  G4DeexPrecoParameters* deexParams = G4NuclearLevelData::GetInstance()->GetParameters();
  deexParams->SetCorrelatedGamma(true);
  deexParams->SetStoreAllLevels(true);
  deexParams->SetInternalConversionFlag(true);
  deexParams->SetIsomerProduction(true);
  deexParams->SetMaxLifeTime(DBL_MAX); // Adjustable
}



PhysicsList::~PhysicsList()
{ }



void PhysicsList::SetCuts()
{
  SetCutValue(1*um, "proton");
  SetCutValue(1*um, "e-");
  SetCutValue(1*um, "e+");
  SetCutValue(1*um, "gamma");      
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......