| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file StepMaxProcess.cc
27 /// \brief Implementation of the StepMaxProces
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo
31 //....oooOO0OOooo........oooOO0OOooo........oo
32
33 #include "StepMaxProcess.hh"
34
35 #include "HistoManager.hh"
36 #include "StepMaxMessenger.hh"
37
38 //....oooOO0OOooo........oooOO0OOooo........oo
39
40 StepMaxProcess::StepMaxProcess(const G4String&
41 : G4VDiscreteProcess(processName, type)
42 {
43 fMaxStep1 = fMaxStep2 = DBL_MAX;
44 fApplyMaxStep2 = true;
45 fMess = new StepMaxMessenger(this);
46 }
47
48 //....oooOO0OOooo........oooOO0OOooo........oo
49
50 StepMaxProcess::~StepMaxProcess()
51 {
52 delete fMess;
53 }
54
55 //....oooOO0OOooo........oooOO0OOooo........oo
56
57 G4bool StepMaxProcess::IsApplicable(const G4Pa
58 {
59 return (particle.GetPDGCharge() != 0. && !pa
60 }
61
62 //....oooOO0OOooo........oooOO0OOooo........oo
63
64 void StepMaxProcess::SetMaxStep1(G4double step
65 {
66 fMaxStep1 = step;
67 }
68
69 //....oooOO0OOooo........oooOO0OOooo........oo
70
71 void StepMaxProcess::ApplyMaxStep2(G4bool valu
72 {
73 fApplyMaxStep2 = value;
74 }
75
76 //....oooOO0OOooo........oooOO0OOooo........oo
77
78 G4double StepMaxProcess::PostStepGetPhysicalIn
79
80 {
81 // condition is set to "Not Forced"
82 *condition = NotForced;
83
84 if (fApplyMaxStep2) {
85 G4AnalysisManager* analysisManager = G4Ana
86 G4int ih = 10;
87 if (analysisManager->GetH1Activation(ih))
88 fMaxStep2 = analysisManager->GetH1Width(
89 return fMaxStep2;
90 }
91 else
92 return fMaxStep1;
93 }
94
95 //....oooOO0OOooo........oooOO0OOooo........oo
96
97 G4VParticleChange* StepMaxProcess::PostStepDoI
98 {
99 // do nothing
100 aParticleChange.Initialize(aTrack);
101 return &aParticleChange;
102 }
103
104 //....oooOO0OOooo........oooOO0OOooo........oo
105