]>
Commit | Line | Data |
---|---|---|
b94a914b | 1 | // $Id$ |
2 | // Category: physics | |
3 | // | |
4 | // According to ExN06PhysicsList (geant4 1.1) | |
5 | ||
6 | #include "TG4PhysicsConstructorOptical.h" | |
7 | ||
8 | #include <G4ParticleDefinition.hh> | |
9 | #include <G4ProcessManager.hh> | |
10 | #include <G4Cerenkov.hh> | |
11 | #include <G4OpAbsorption.hh> | |
12 | #include <G4OpRayleigh.hh> | |
13 | #include <G4OpBoundaryProcess.hh> | |
14 | ||
15 | TG4PhysicsConstructorOptical::TG4PhysicsConstructorOptical(const G4String& name) | |
16 | : G4VPhysicsConstructor(name) | |
17 | { | |
18 | // | |
19 | SetVerboseLevel(1); | |
20 | } | |
21 | ||
22 | TG4PhysicsConstructorOptical::~TG4PhysicsConstructorOptical() { | |
23 | // | |
24 | } | |
25 | ||
26 | // protected methods | |
27 | ||
28 | void TG4PhysicsConstructorOptical::ConstructParticle() | |
29 | { | |
30 | // The particles are constructed in the | |
31 | // TG4ModularPhysicsList. | |
32 | // --- | |
33 | } | |
34 | ||
35 | void TG4PhysicsConstructorOptical::ConstructProcess() | |
36 | { | |
37 | // Constructs optical processes. | |
38 | // According to ExN06PhysicsList.cc. | |
39 | // (geant4 1.1) | |
40 | // --- | |
41 | ||
42 | G4Cerenkov* theCerenkovProcess = new G4Cerenkov("Cerenkov"); | |
43 | G4OpAbsorption* theAbsorptionProcess = new G4OpAbsorption(); | |
44 | G4OpRayleigh* theRayleighScatteringProcess = new G4OpRayleigh(); | |
45 | G4OpBoundaryProcess* theBoundaryProcess = new G4OpBoundaryProcess(); | |
46 | ||
47 | theCerenkovProcess->DumpPhysicsTable(); | |
48 | //theAbsorptionProcess->DumpPhysicsTable(); | |
49 | //theRayleighScatteringProcess->DumpPhysicsTable(); | |
50 | ||
51 | // add verbose | |
52 | //theCerenkovProcess->SetVerboseLevel(1); | |
53 | //theAbsorptionProcess->SetVerboseLevel(1); | |
54 | //theRayleighScatteringProcess->SetVerboseLevel(1); | |
55 | //theBoundaryProcess->SetVerboseLevel(1); | |
56 | ||
57 | G4int maxNumPhotons = 300; | |
58 | ||
59 | theCerenkovProcess->SetTrackSecondariesFirst(true); | |
60 | theCerenkovProcess->SetMaxNumPhotonsPerStep(maxNumPhotons); | |
61 | ||
62 | //G4OpticalSurfaceModel themodel = unified; | |
63 | // model from GEANT3 | |
64 | G4OpticalSurfaceModel themodel = glisur; | |
65 | theBoundaryProcess->SetModel(themodel); | |
66 | ||
67 | theParticleIterator->reset(); | |
68 | while( (*theParticleIterator)() ){ | |
69 | G4ParticleDefinition* particle = theParticleIterator->value(); | |
70 | G4ProcessManager* processManager = particle->GetProcessManager(); | |
71 | G4String particleName = particle->GetParticleName(); | |
72 | if (theCerenkovProcess->IsApplicable(*particle)) { | |
73 | processManager->AddContinuousProcess(theCerenkovProcess); | |
74 | } | |
75 | if (particleName == "opticalphoton") { | |
76 | G4cout << " AddDiscreteProcess to OpticalPhoton " << G4endl; | |
77 | processManager->AddDiscreteProcess(theAbsorptionProcess); | |
78 | processManager->AddDiscreteProcess(theRayleighScatteringProcess); | |
79 | processManager->AddDiscreteProcess(theBoundaryProcess); | |
80 | } | |
81 | } | |
82 | ||
83 | if (verboseLevel>0) | |
84 | G4cout << "### " << namePhysics << " physics constructed." << G4endl; | |
85 | } | |
86 |