nam/muon-tomograpy
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

change muon energy

Browse Source
This commit is contained in:
Nam Tran
2021-01-15 11:44:40 -06:00
parent 8e081d10e1
commit 4c6a00e681
2 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
src/DetectorConstruction.cc
View File

@@ -51,12 +51,12 @@ G4VPhysicalVolume *DetectorConstruction::Construct() {
G4double layerSizeY = nScBar * scSizeX + nScBar * scGapX; G4double layerSizeY = nScBar * scSizeX + nScBar * scGapX;
G4double layerSizeZ = scSizeZ + scGapZ; G4double layerSizeZ = scSizeZ + scGapZ;
G4double concreteX = 20 * cm; G4double concreteX = 30 * cm;
G4double concreteY = 20 * cm; G4double concreteY = 30 * cm;
G4double concreteZ = 20 * cm; G4double concreteZ = 30 * cm;
G4double targetX = 5 * cm; G4double targetX = 10 * cm;
G4double targetY = 5 * cm; G4double targetY = 10 * cm;
G4double targetZ = 2 * cm; G4double targetZ = 5 * cm;
G4double worldSizeX = 1.2 * layerSizeX; G4double worldSizeX = 1.2 * layerSizeX;
G4double worldSizeY = 1.2 * layerSizeY; G4double worldSizeY = 1.2 * layerSizeY;
@@ -155,7 +155,7 @@ G4VPhysicalVolume *DetectorConstruction::Construct() {
logicScX->SetVisAttributes(new G4VisAttributes(G4Colour(.0, 1.0, 1., 1.0))); logicScX->SetVisAttributes(new G4VisAttributes(G4Colour(.0, 1.0, 1., 1.0)));
logicScY->SetVisAttributes(new G4VisAttributes(G4Colour(0.0, 0.0, 1., 0.4))); logicScY->SetVisAttributes(new G4VisAttributes(G4Colour(0.0, 0.0, 1., 0.4)));
logConcrete->SetVisAttributes(new G4VisAttributes(G4Colour(.8, .3, .3, 0.3))); logConcrete->SetVisAttributes(new G4VisAttributes(G4Colour(.8, .3, .3, 0.3)));
logTarget->SetVisAttributes(new G4VisAttributes(G4Colour(1., .3, .3, 1.0))); logTarget->SetVisAttributes(new G4VisAttributes(G4Colour(.3, 1., .3, 1.0)));
// always return the physical World // always return the physical World
// //

2
src/PrimaryGeneratorAction.cc
View File

@@ -25,7 +25,7 @@ PrimaryGeneratorAction::PrimaryGeneratorAction()
particleTable->FindParticle(particleName = "mu+"); particleTable->FindParticle(particleName = "mu+");
fParticleGun->SetParticleDefinition(particle); fParticleGun->SetParticleDefinition(particle);
fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0., 0., -1.)); fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0., 0., -1.));
fParticleGun->SetParticleEnergy(6. * GeV); fParticleGun->SetParticleEnergy(4. * GeV);
} }
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......