Geant4 LXR

Cross-Referencing   Geant4
Geant4/examples/advanced/purging_magnet/src/PurgMagPhysicsList.cc

   //
   // ********************************************************************
   // * License and Disclaimer                                           *
   // *                                                                  *
   // * The  Geant4 software  is  copyright of the Copyright Holders  of *
   // * the Geant4 Collaboration.  It is provided  under  the terms  and *
   // * conditions of the Geant4 Software License,  included in the file *
   // * LICENSE and available at  http://cern.ch/geant4/license .  These *
   // * include a list of copyright holders.                             *
  // *                                                                  *
  // * Neither the authors of this software system, nor their employing *
  // * institutes,nor the agencies providing financial support for this *
  // * work  make  any representation or  warranty, express or implied, *
  // * regarding  this  software system or assume any liability for its *
  // * use.  Please see the license in the file  LICENSE  and URL above *
  // * for the full disclaimer and the limitation of liability.         *
  // *                                                                  *
  // * This  code  implementation is the result of  the  scientific and *
  // * technical work of the GEANT4 collaboration.                      *
  // * By using,  copying,  modifying or  distributing the software (or *
  // * any work based  on the software)  you  agree  to acknowledge its *
  // * use  in  resulting  scientific  publications,  and indicate your *
  // * acceptance of all terms of the Geant4 Software license.          *
  // ********************************************************************
  //
  // Code developed by:
  //  S.Larsson
  //
  //    ********************************
  //    *                              *
  //    *    PurgMagPhysicsList.cc     *
  //    *                              *
  //    ********************************
  //
  
  #include "PurgMagPhysicsList.hh"
  #include "G4SystemOfUnits.hh"
  #include "G4ParticleDefinition.hh"
  #include "G4ParticleWithCuts.hh"
  #include "G4ProcessManager.hh"
  #include "G4ParticleTypes.hh"
  #include "G4ParticleTable.hh"
  #include "G4Material.hh"
  #include "G4UnitsTable.hh"
  #include "G4ios.hh"              
  #include "G4EmStandardPhysics_option4.hh"
  #include "G4VPhysicsConstructor.hh"
  #include "G4DecayPhysics.hh"
  
  PurgMagPhysicsList::PurgMagPhysicsList():  G4VUserPhysicsList()
  {
    defaultCutValue = 1*micrometer;
    cutForGamma     = defaultCutValue;
    cutForElectron  = defaultCutValue;
    cutForPositron  = defaultCutValue;
    cutForProton    = defaultCutValue;
  
    fEmPhysicsList = new G4EmStandardPhysics_option4();
    fDecPhysicsList = new G4DecayPhysics();
    SetVerboseLevel(1);
  }
  
  PurgMagPhysicsList::~PurgMagPhysicsList()
  {
   delete fDecPhysicsList;
   delete fEmPhysicsList;
  }
  
  void PurgMagPhysicsList::ConstructParticle()
  {
   fDecPhysicsList -> ConstructParticle();
  } 
  
  void PurgMagPhysicsList::ConstructProcess()
  {
    AddTransportation();
    fEmPhysicsList -> ConstructProcess();
    
  // Deexcitation
  // Both Fluorescence and Auger e- emission activated
  //G4VAtomDeexcitation* de = new G4UAtomicDeexcitation();
  //G4LossTableManager::Instance()->SetAtomDeexcitation(de);
  //de -> SetFluo(true);
  //de -> SetAuger(true);
  }
  
  void PurgMagPhysicsList::SetCuts()
  {
    if (verboseLevel >0){
      G4cout << "PurgMagPhysicsList::SetCuts:";
      G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
    }  
    
    // set cut values for gamma at first and for e- second and next for e+,
    // because some processes for e+/e- need cut values for gamma 
    SetCutValue(cutForGamma, "gamma");
    SetCutValue(cutForElectron, "e-");
    SetCutValue(cutForPositron, "e+");
    
   // set cut values for proton and anti_proton before all other hadrons
   // because some processes for hadrons need cut values for proton/anti_proton 
   SetCutValue(cutForProton, "proton");
   SetCutValue(cutForProton, "anti_proton");
   
   //  SetCutValueForOthers(defaultCutValue);
   
   if (verboseLevel>0) DumpCutValuesTable();
 }
 
 void PurgMagPhysicsList::SetGammaLowLimit(G4double lowcut)
 {
   if (verboseLevel >0){
     G4cout << "PurgMagPhysicsList::SetCuts:";
     G4cout << "Gamma cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
   }  
 
   // G4Gamma::SetEnergyRange(lowcut,1e5); 
   SetGELowLimit(lowcut);
 }
 
 void PurgMagPhysicsList::SetElectronLowLimit(G4double lowcut)
 {
   if (verboseLevel >0){
     
     G4cout << "PurgMagPhysicsList::SetCuts:";
     G4cout << "Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
   }  
 
   // G4Electron::SetEnergyRange(lowcut,1e5);
   SetGELowLimit(lowcut);
 }
 
 void PurgMagPhysicsList::SetPositronLowLimit(G4double lowcut)
 {
   if (verboseLevel >0){
     
     G4cout << "PurgMagPhysicsList::SetCuts:";
     G4cout << "Positron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
   }  
 
   G4cerr << "PurgMagPhysicsList::SetPositronLowLimit: Not currently able to set Positron LowLimit." << G4endl;
   G4Exception("PurgMagPhysicsList::SetPositronLowLimit()","PurMag001",
         FatalException,"Positron Low Limit: not implemented in PurgMagPhysicsList"); 
   //
   // G4Positron::SetEnergyRange(lowcut,1e5);
 }
 
 
 void PurgMagPhysicsList::SetProtonLowLimit(G4double lowcut)
 {
   if (verboseLevel >0){
     
     G4cout << "PurgMagPhysicsList::SetCuts:";
     G4cout << "Proton cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;  
   }  
 
   G4cerr << "PurgMagPhysicsList::SetProtonLowLimit: Not currently able to set Proton LowLimit." << G4endl;
   G4Exception("PurgMagPhysicsList::SetProtonLowLimit()","PurMag002",
         FatalException,"Proton Low Limit: not implemented in PurgMagPhysicsList"); 
   //
   // G4Proton::SetEnergyRange(lowcut,1e5);
   // G4AntiProton::SetEnergyRange(lowcut,1e5);
 }
 
 void PurgMagPhysicsList::SetGEPLowLimit(G4double lowcut)
 {
   if (verboseLevel >0){
     G4cout << "PurgMagPhysicsList::SetGEPLowLimit:";
     G4cout << "Gamma and Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
   }  
 
   // G4Gamma::SetEnergyRange(lowcut,1e5);
   // G4Electron::SetEnergyRange(lowcut,1e5);
   // G4Positron::SetEnergyRange(lowcut,1e5);
   this->SetGELowLimit(lowcut); 
 
   G4cerr << " SetGEPLowLimit : Uncertain whether setting Positron low limit " << G4endl;
 }
 
 void PurgMagPhysicsList::SetGELowLimit(G4double lowcut)
 {
   if (verboseLevel >0){
     G4cout << "PurgMagPhysicsList::SetGELowLimit:";
     G4cout << "Gamma and Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
   }  
  
   G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(lowcut,1e5);
 }
 void PurgMagPhysicsList::SetGammaCut(G4double val)
 {
   cutForGamma = val;
 }
 
 void PurgMagPhysicsList::SetElectronCut(G4double val)
 {
   cutForElectron = val;
 }
 
 void PurgMagPhysicsList::SetPositronCut(G4double val)
 {
   cutForPositron = val;
 }
 
 void PurgMagPhysicsList::SetProtonCut(G4double val)
 {
   cutForProton = val;
 }