// // ******************************************************************** // * 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. * // ******************************************************************** // // // $Id: XrayFluoPhysicsList.cc // GEANT4 tag $Name: xray_fluo-V03-02-00 // // Author: Elena Guardincerri (Elena.Guardincerri@ge.infn.it) // // History: // ----------- // 28 Nov 2001 Elena Guardincerri Created // // ------------------------------------------------------------------- #include "XrayFluoPhysicsList.hh" #include "XrayFluoPhysicsListMessenger.hh" #include "XrayFluoDetectorConstruction.hh" #include "XrayFluoPlaneDetectorConstruction.hh" #include "XrayFluoMercuryDetectorConstruction.hh" ///////////////////////////////////////// //#include "G4LeptonConstructor.hh" //#include "G4BosonConstructor.hh" //#include "G4MesonConstructor.hh" #include "G4BaryonConstructor.hh" ///////////////////////////////////////// #include "G4ParticleDefinition.hh" #include "G4ParticleWithCuts.hh" #include "G4ProcessManager.hh" #include "G4ParticleTypes.hh" #include "G4ParticleTable.hh" #include "G4ios.hh" #include "G4Region.hh" #include "G4RegionStore.hh" //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... XrayFluoPhysicsList::XrayFluoPhysicsList(XrayFluoDetectorConstruction* p) : G4VUserPhysicsList(),pDet(0),planeDet(0),mercuryDet(0) { pDet = p; // SetGELowLimit(250*eV); defaultCutValue = 10e-6*mm; cutForGamma = defaultCutValue; cutForElectron = defaultCutValue; cutForProton = 0.001*mm; SetVerboseLevel(1); physicsListMessenger = new XrayFluoPhysicsListMessenger(this); G4String regName = "SampleRegion"; G4double cutValue = 0.000001 * mm; G4Region* reg = G4RegionStore::GetInstance()->GetRegion(regName); G4ProductionCuts* cuts = new G4ProductionCuts; cuts->SetProductionCut(cutValue); reg->SetProductionCuts(cuts); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... XrayFluoPhysicsList::XrayFluoPhysicsList(XrayFluoPlaneDetectorConstruction* p) : G4VUserPhysicsList(),pDet(0),planeDet(0),mercuryDet(0) { planeDet = p; // SetGELowLimit(250*eV); defaultCutValue = 0.000001*mm; cutForGamma = defaultCutValue; cutForElectron = defaultCutValue; cutForProton = 0.001*mm; SetVerboseLevel(1); physicsListMessenger = new XrayFluoPhysicsListMessenger(this); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... XrayFluoPhysicsList::XrayFluoPhysicsList(XrayFluoMercuryDetectorConstruction* p) : G4VUserPhysicsList(),pDet(0),planeDet(0),mercuryDet(0) { mercuryDet = p; // SetGELowLimit(250*eV); defaultCutValue = 0.000001*mm; cutForGamma = defaultCutValue; cutForElectron = defaultCutValue; cutForProton = 0.001*mm; SetVerboseLevel(1); physicsListMessenger = new XrayFluoPhysicsListMessenger(this); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... XrayFluoPhysicsList::~XrayFluoPhysicsList() { delete physicsListMessenger; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::ConstructParticle() { // In this method, static member functions should be called // for all particles which you want to use. // This ensures that objects of these particle types will be // created in the program. ConstructBosons(); ConstructLeptons(); ConstructBarions(); ConstructIons(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::ConstructBosons() { // gamma G4Gamma::GammaDefinition(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::ConstructLeptons() { // leptons G4Electron::ElectronDefinition(); G4Positron::PositronDefinition(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::ConstructBarions() { G4BaryonConstructor baryon; baryon.ConstructParticle(); // G4Proton::ProtonDefinition(); } void XrayFluoPhysicsList::ConstructIons() { // Ions G4Alpha::AlphaDefinition(); } void XrayFluoPhysicsList::ConstructProcess() { AddTransportation(); ConstructEM(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... #include "G4LowEnergyCompton.hh" #include "G4LowEnergyGammaConversion.hh" #include "G4LowEnergyPhotoElectric.hh" #include "G4LowEnergyRayleigh.hh" // e+ #include "G4MultipleScattering.hh" #include "G4eIonisation.hh" #include "G4eBremsstrahlung.hh" #include "G4eplusAnnihilation.hh" #include "G4LowEnergyIonisation.hh" #include "G4LowEnergyBremsstrahlung.hh" #include "G4hLowEnergyIonisation.hh" //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::ConstructEM() { theParticleIterator->reset(); while( (*theParticleIterator)() ){ G4ParticleDefinition* particle = theParticleIterator->value(); G4ProcessManager* pmanager = particle->GetProcessManager(); G4String particleName = particle->GetParticleName(); if (particleName == "gamma") { // gamma pmanager->AddDiscreteProcess(new G4LowEnergyCompton); LePeprocess = new G4LowEnergyPhotoElectric(); LePeprocess->ActivateAuger(true); LePeprocess->SetCutForLowEnSecPhotons(0.250 * keV); LePeprocess->SetCutForLowEnSecElectrons(0.250 * keV); pmanager->AddDiscreteProcess(LePeprocess); pmanager->AddDiscreteProcess(new G4LowEnergyRayleigh("Rayleigh")); } else if (particleName == "e-") { //electron pmanager->AddProcess(new G4MultipleScattering,-1, 1,1); LeIoprocess = new G4LowEnergyIonisation("IONI"); LeIoprocess->ActivateAuger(true); //eIoProcess = new G4eIonisation("stdIONI"); LeIoprocess->SetCutForLowEnSecPhotons(0.1*keV); LeIoprocess->SetCutForLowEnSecElectrons(0.1*keV); pmanager->AddProcess(LeIoprocess, -1, 2, 2); LeBrprocess = new G4LowEnergyBremsstrahlung(); pmanager->AddProcess(LeBrprocess, -1, -1, 3); } else if (particleName == "e+") { //positron pmanager->AddProcess(new G4MultipleScattering,-1, 1,1); pmanager->AddProcess(new G4eIonisation, -1, 2,2); pmanager->AddProcess(new G4eBremsstrahlung, -1,-1,3); pmanager->AddProcess(new G4eplusAnnihilation, 0,-1,4); } else if (particleName == "proton") { //proton G4hLowEnergyIonisation* hIoni = new G4hLowEnergyIonisation; hIoni->SetFluorescence(true); pmanager->AddProcess(new G4MultipleScattering,-1,1,1); pmanager->AddProcess(hIoni,-1, 2,2); } else if ( particleName == "alpha" ) { pmanager->AddProcess(new G4MultipleScattering,-1,1,1); G4hLowEnergyIonisation* iIon = new G4hLowEnergyIonisation() ; pmanager->AddProcess(iIon,-1,2,2); } } } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::SetGELowLimit(G4double lowcut) { if (verboseLevel >0){ G4cout << "XrayFluoPhysicsList::SetCuts:"; G4cout << "Gamma and Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl; } G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(lowcut,1e5); } void XrayFluoPhysicsList::SetGammaLowLimit(G4double lowcut) { if (verboseLevel >0){ G4cout << "XrayFluoPhysicsList::SetCuts:"; G4cout << "Gamma cut in energy: " << lowcut*MeV << " (MeV)" << G4endl; } SetGELowLimit(lowcut); // G4Gamma::SetEnergyRange(lowcut,1e5); } void XrayFluoPhysicsList::SetElectronLowLimit(G4double lowcut) { if (verboseLevel >0){ G4cout << "XrayFluoPhysicsList::SetCuts:"; G4cout << "Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl; } SetGELowLimit(lowcut); // G4Electron::SetEnergyRange(lowcut,1e5); } void XrayFluoPhysicsList::SetGammaCut(G4double val) { ResetCuts(); cutForGamma = val; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::SetElectronCut(G4double val) { ResetCuts(); cutForElectron = val; } void XrayFluoPhysicsList::SetCuts(){ SetCutValue(cutForGamma,"gamma"); SetCutValue(cutForElectron,"e-"); SetCutValue(cutForElectron,"e+"); SetCutValue(cutForProton, "proton"); //SetCutValueForOthers(cutForProton); if (verboseLevel>0) DumpCutValuesTable(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::SetLowEnSecPhotCut(G4double cut){ G4cout<<"Low energy secondary photons cut is now set to: "<SetCutForLowEnSecPhotons(cut); LePeprocess->SetCutForLowEnSecPhotons(cut); LeIoprocess->SetCutForLowEnSecPhotons(cut); } void XrayFluoPhysicsList::SetLowEnSecElecCut(G4double cut){ G4cout<<"Low energy secondary electrons cut is now set to: "<SetCutForLowEnSecElectrons(cut); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::SetProtonCut(G4double val) { ResetCuts(); cutForProton = val; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... void XrayFluoPhysicsList::SetCutsByEnergy(G4double val) { G4ParticleTable* theXrayFluoParticleTable = G4ParticleTable::GetParticleTable(); G4Material* currMat = 0; if(pDet){ currMat = pDet->XrayFluoDetectorConstruction::GetSampleMaterial(); } else if(planeDet){ currMat = planeDet->XrayFluoPlaneDetectorConstruction::GetPlaneMaterial(); } else if(mercuryDet){ currMat = mercuryDet->GetMercuryMaterial(); } G4ParticleDefinition* part; G4double cut; part = theXrayFluoParticleTable->FindParticle("e-"); cut = G4EnergyLossTables::GetRange(part,val,currMat); SetCutValue(cut, "e-"); part = theXrayFluoParticleTable->FindParticle("proton"); cut = G4EnergyLossTables::GetRange(part,val,currMat); SetCutValue(cut, "proton"); // part = theXrayFluoParticleTable->FindParticle("gamma"); // cut = G4EnergyLossTables::GetRange(part,val,currMat); // SetCutValue(cut, "gamma"); }