Changeset 1337 for trunk/examples/extended/radioactivedecay/exrdm/src
- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- Location:
- trunk/examples/extended/radioactivedecay/exrdm/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/radioactivedecay/exrdm/src/exrdmEventAction.cc
r1230 r1337 32 32 #include "G4EventManager.hh" 33 33 #include "G4HCofThisEvent.hh" 34 #include "G4TrajectoryContainer.hh"35 #include "G4Trajectory.hh"36 #include "G4VVisManager.hh"37 34 #include "G4UImanager.hh" 38 35 #include "G4UnitsTable.hh" … … 65 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 66 63 67 void exrdmEventAction::EndOfEventAction(const G4Event* evt)64 void exrdmEventAction::EndOfEventAction(const G4Event*) 68 65 { 69 66 //analysis 70 67 exrdmAnalysisManager::getInstance()->EndOfEvent(); 71 // visualisation72 #ifdef G4VIS_USE73 G4int event_id = evt->GetEventID();74 if (event_id < 100 && G4VVisManager::GetConcreteInstance()) {75 G4TrajectoryContainer * trajectoryContainer = evt->GetTrajectoryContainer();76 G4int n_trajectories = 0;77 if (trajectoryContainer) n_trajectories = trajectoryContainer->entries();78 for (G4int i=0; i<n_trajectories; i++) {79 G4Trajectory* trj = (G4Trajectory *)80 ((*(evt->GetTrajectoryContainer()))[i]);81 if (drawFlag == "all") {82 trj->DrawTrajectory(0);83 } else if (drawFlag == "charged" && trj->GetCharge() != 0.) {84 trj->DrawTrajectory(0);85 }86 }87 }88 #endif89 68 } 90 69 -
trunk/examples/extended/radioactivedecay/exrdm/src/exrdmPhysListEmLowEnergy.cc
r807 r1337 46 46 #include "G4eplusAnnihilation.hh" 47 47 48 #include "G4MultipleScattering.hh" 48 #include "G4eMultipleScattering.hh" 49 #include "G4MuMultipleScattering.hh" 50 #include "G4hMultipleScattering.hh" 49 51 50 52 #include "G4MuIonisation.hh" … … 88 90 } else if (particleName == "e-") { 89 91 90 pmanager->AddProcess(new G4 MultipleScattering, -1, 1,1);92 pmanager->AddProcess(new G4eMultipleScattering, -1, 1,1); 91 93 pmanager->AddProcess(new G4LowEnergyIonisation, -1, 2,2); 92 94 pmanager->AddProcess(new G4LowEnergyBremsstrahlung, -1,-3,3); … … 94 96 } else if (particleName == "e+") { 95 97 96 pmanager->AddProcess(new G4 MultipleScattering, -1, 1,1);98 pmanager->AddProcess(new G4eMultipleScattering, -1, 1,1); 97 99 pmanager->AddProcess(new G4eIonisation, -1, 2,2); 98 100 pmanager->AddProcess(new G4eBremsstrahlung, -1, 3,3); … … 102 104 particleName == "mu-" ) { 103 105 104 pmanager->AddProcess(new G4Mu ltipleScattering,-1, 1,1);106 pmanager->AddProcess(new G4MuMultipleScattering,-1, 1,1); 105 107 pmanager->AddProcess(new G4MuIonisation, -1, 2,2); 106 108 pmanager->AddProcess(new G4MuBremsstrahlung, -1, 3,3); … … 109 111 } else if (particleName == "GenericIon") { 110 112 111 pmanager->AddProcess(new G4 MultipleScattering, -1, 1,1);113 pmanager->AddProcess(new G4hMultipleScattering, -1, 1,1); 112 114 pmanager->AddProcess(new G4hLowEnergyIonisation, -1,2,2); 113 115 // pmanager->AddProcess(new G4ionIonisation, -1, 2,2); … … 117 119 (particle->GetParticleName() != "chargedgeantino")) { 118 120 119 pmanager->AddProcess(new G4 MultipleScattering,-1,1,1);121 pmanager->AddProcess(new G4hMultipleScattering,-1,1,1); 120 122 pmanager->AddProcess(new G4hLowEnergyIonisation, -1,2,2); 121 123 }
Note: See TracChangeset
for help on using the changeset viewer.