| 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 //
27 /// \file ClassifiedDamage.cc
28 /// \brief Implementation of the ClassifiedDam
29
30 #include "ClassifiedDamage.hh"
31
32 //....oooOO0OOooo........oooOO0OOooo........oo
33
34 ClassifiedDamage::ClassifiedDamage()
35 {
36 fType = fNA;
37 fDamage.clear();
38 fBp_begin = 0;
39 fBp_end = 0;
40 fComplexity = -1;
41 fIncludeBase = false;
42 }
43
44 //....oooOO0OOooo........oooOO0OOooo........oo
45
46 void ClassifiedDamage::AddDamage(Damage pDmg)
47 {
48 fDamage.push_back(pDmg);
49 }
50
51 //....oooOO0OOooo........oooOO0OOooo........oo
52
53 void ClassifiedDamage::ComputeBp()
54 {
55 fBp_begin = fDamage[0].GetCopyNb();
56 fBp_end = fDamage[fDamage.size()-1].GetCopyN
57
58 fBp_barycenter=0;
59 for(auto it=fDamage.begin();it!=fDamage.end(
60 {
61 fBp_barycenter+=it->GetCopyNb();
62 }
63 fBp_barycenter = fBp_barycenter/fDamage.size
64 }
65
66 //....oooOO0OOooo........oooOO0OOooo........oo
67
68 void ClassifiedDamage::ComputePosition()
69 {
70 // TODO
71 }
72
73 //....oooOO0OOooo........oooOO0OOooo........oo
74
75 void ClassifiedDamage::ComputeComplexity()
76 {
77 if(fDamage.size()==0)
78 {
79 fComplexity = -1;
80 }
81 else
82 {
83 fComplexity = -1;
84 for(auto it=fDamage.begin();it!=fDamage.en
85 {
86 if((it->GetDamageType()==Damage::DamageT
87 {
88 fComplexity++;
89 }
90 else
91 {
92 if(fIncludeBase)
93 {
94 fComplexity++;
95 }
96 }
97 }
98 }
99 }
100
101 //....oooOO0OOooo........oooOO0OOooo........oo
102
103 void ClassifiedDamage::ComputeType()
104 {
105 bool firstStrandTouched = false;
106 bool secondStrandTouched = false;
107
108 if(fDamage.size()==0)
109 {
110 fType = fNA;
111 }
112 else
113 {
114 for(auto it=fDamage.begin();it!=fDamage.en
115 {
116 if((it->GetDamageType()==Damage::DamageT
117 {
118 int strand = it->GetStrand();
119 if(strand == 1)
120 {
121 firstStrandTouched = true;
122 }
123 if(strand == 2)
124 {
125 secondStrandTouched = true;
126 }
127 }
128
129 if (it->GetCause() == Damage::DamageCaus
130 fIsThereDirectContribution = true;
131 }
132
133 if (it->GetCause() == Damage::DamageCaus
134 fIsThereIndirectContribution = true;
135 }
136 }
137
138 if(firstStrandTouched && secondStrandTouch
139 {
140 fType = fDSB;
141 }
142 else
143 {
144 fType = fSSB;
145 }
146 }
147 }
148
149 //....oooOO0OOooo........oooOO0OOooo........oo
150
151 void ClassifiedDamage::Reset()
152 {
153 fType = fNA;
154 fDamage.clear();
155 fBp_begin = 0;
156 fBp_end = 0;
157 fComplexity = -1;
158 }
159