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