LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Added Si3N4 material, added G4hMultipleScattering process.

Browse Source
This commit is contained in:
suter_a
2020-04-20 14:41:16 +02:00
parent e9b675cc0a
commit 9cf33e2406
2 changed files with 7 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/musrDetectorConstruction.cc
View File

@ -1609,6 +1609,11 @@ void musrDetectorConstruction::DefineMaterials()
SiO2->AddElement(O, natoms=2); SiO2->AddElement(O, natoms=2);
SiO2->AddElement(Si, natoms=1); SiO2->AddElement(Si, natoms=1);
G4Material* Si3N4 = new G4Material("Si3N4", 3.2*CLHEP::g/CLHEP::cm3, ncomponents=2);
Si3N4->AddElement(Si, natoms=3);
Si3N4->AddElement(N, natoms=4);
G4Material* Al2O3 = new G4Material("Al2O3", 3.985*CLHEP::g/CLHEP::cm3, ncomponents=2); // saphire G4Material* Al2O3 = new G4Material("Al2O3", 3.985*CLHEP::g/CLHEP::cm3, ncomponents=2); // saphire
Al2O3->AddElement (Al, natoms=2); Al2O3->AddElement (Al, natoms=2);
Al2O3->AddElement (O, natoms=3); Al2O3->AddElement (O, natoms=3);

3
src/musrPhysicsList.cc
View File

@ -381,7 +381,8 @@ void musrPhysicsList::ConstructEM()
else if (stringProcessName=="G4StepLimiter") pManager->AddProcess(new G4StepLimiter,nr1,nr2,nr3); else if (stringProcessName=="G4StepLimiter") pManager->AddProcess(new G4StepLimiter,nr1,nr2,nr3);
else if (stringProcessName=="G4UserSpecialCuts") pManager->AddProcess(new G4UserSpecialCuts,nr1,nr2,nr3); else if (stringProcessName=="G4UserSpecialCuts") pManager->AddProcess(new G4UserSpecialCuts,nr1,nr2,nr3);
// else if (stringProcessName=="G4DecayWithSpin") pManager->AddProcess(new G4DecayWithSpin,nr1,nr2,nr3); // else if (stringProcessName=="G4DecayWithSpin") pManager->AddProcess(new G4DecayWithSpin,nr1,nr2,nr3);
// else if (stringProcessName=="G4hIonisation") pManager->AddProcess(new G4hIonisation,nr1,nr2,nr3); else if (stringProcessName=="G4hIonisation") pManager->AddProcess(new G4hIonisation,nr1,nr2,nr3);
else if (stringProcessName=="G4hMultipleScattering") pManager->AddProcess(new G4hMultipleScattering,nr1,nr2,nr3);
// else if (stringProcessName=="G4hLowEnergyIonisation") pManager->AddProcess(new G4hLowEnergyIonisation,nr1,nr2,nr3); // else if (stringProcessName=="G4hLowEnergyIonisation") pManager->AddProcess(new G4hLowEnergyIonisation,nr1,nr2,nr3);
else if (stringProcessName=="musrMuFormation") pManager->AddProcess(new musrMuFormation,nr1,nr2,nr3); else if (stringProcessName=="musrMuFormation") pManager->AddProcess(new musrMuFormation,nr1,nr2,nr3);
else if (stringProcessName=="musrMuEnergyLossLandau") pManager->AddProcess(new musrMuEnergyLossLandau,nr1,nr2,nr3); else if (stringProcessName=="musrMuEnergyLossLandau") pManager->AddProcess(new musrMuEnergyLossLandau,nr1,nr2,nr3);