Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // 23 // 27 // ------------------------------------------- 24 // ------------------------------------------------------------------- 28 // GEANT4 Class file 25 // GEANT4 Class file 29 // 26 // 30 // For information related to this code c 27 // For information related to this code contact: 31 // 28 // 32 // File name: G4PartialWidthTable 29 // File name: G4PartialWidthTable 33 // 30 // 34 // Author: 31 // Author: 35 // 32 // 36 // Creation date: 15 April 1999 33 // Creation date: 15 April 1999 37 // 34 // 38 // Modifications: 35 // Modifications: 39 // 36 // 40 // ------------------------------------------- 37 // ------------------------------------------------------------------- 41 38 42 #include "globals.hh" 39 #include "globals.hh" 43 #include "G4SystemOfUnits.hh" << 44 #include "G4HadronicException.hh" 40 #include "G4HadronicException.hh" 45 #include "G4ios.hh" 41 #include "G4ios.hh" 46 #include "G4PartialWidthTable.hh" 42 #include "G4PartialWidthTable.hh" 47 #include "G4KineticTrack.hh" 43 #include "G4KineticTrack.hh" 48 44 >> 45 49 G4PartialWidthTable::G4PartialWidthTable(const 46 G4PartialWidthTable::G4PartialWidthTable(const G4double* energies, G4int entries) : nEnergies(entries) 50 { 47 { 51 // Fill the vector with tabulated energies 48 // Fill the vector with tabulated energies 52 G4int i; 49 G4int i; 53 for (i=0; i<entries; i++) 50 for (i=0; i<entries; i++) 54 { 51 { 55 G4double e = *(energies + i) * GeV; 52 G4double e = *(energies + i) * GeV; 56 energy.push_back(e); 53 energy.push_back(e); 57 } 54 } 58 } 55 } 59 56 60 57 61 G4PartialWidthTable::~G4PartialWidthTable() 58 G4PartialWidthTable::~G4PartialWidthTable() 62 { } 59 { } 63 60 64 61 65 G4bool G4PartialWidthTable::operator==(const G 62 G4bool G4PartialWidthTable::operator==(const G4PartialWidthTable &right) const 66 { 63 { 67 return (this == (G4PartialWidthTable *) &rig 64 return (this == (G4PartialWidthTable *) &right); 68 } 65 } 69 66 70 67 71 G4bool G4PartialWidthTable::operator!=(const G 68 G4bool G4PartialWidthTable::operator!=(const G4PartialWidthTable &right) const 72 { 69 { 73 return (this != (G4PartialWidthTable *) &rig 70 return (this != (G4PartialWidthTable *) &right); 74 } 71 } 75 72 76 73 77 G4int G4PartialWidthTable::NumberOfChannels() 74 G4int G4PartialWidthTable::NumberOfChannels() const 78 { 75 { 79 return (G4int)widths.size(); << 76 return widths.size(); 80 } 77 } 81 78 82 79 83 const G4PhysicsVector* G4PartialWidthTable::Wi 80 const G4PhysicsVector* G4PartialWidthTable::Width(const G4String& name1, const G4String& name2) const 84 { 81 { 85 // Returned pointer is not owned by the clie 82 // Returned pointer is not owned by the client >> 83 G4int i; 86 G4PhysicsVector* width = 0; 84 G4PhysicsVector* width = 0; 87 std::size_t n = 0; << 85 G4int n = 0; 88 std::size_t entries = widths.size(); << 86 G4int entries = widths.size(); 89 for (std::size_t i=0; i<entries; ++i) << 87 for (i=0; i<entries; i++) 90 { 88 { 91 if ( (daughter1[i] == name1 && daughter2 89 if ( (daughter1[i] == name1 && daughter2[i] == name2) || 92 (daughter2[i] == name1 && daughter1[i] == 90 (daughter2[i] == name1 && daughter1[i] == name2) ) 93 { 91 { 94 width = (G4PhysicsVector*) (widths[i]); 92 width = (G4PhysicsVector*) (widths[i]); 95 ++n; << 93 n++; 96 } 94 } 97 } 95 } 98 if (n > 1) throw G4HadronicException(__FILE_ 96 if (n > 1) throw G4HadronicException(__FILE__, __LINE__, "G4PartialWidthTable::Width - ambiguity"); 99 97 100 return width; 98 return width; 101 } 99 } 102 100 103 101 104 void G4PartialWidthTable::AddWidths(const G4do 102 void G4PartialWidthTable::AddWidths(const G4double* channelWidth, 105 const G4String& name1, const G4Str 103 const G4String& name1, const G4String& name2) 106 { 104 { 107 G4PhysicsFreeVector* width = new G4PhysicsFr 105 G4PhysicsFreeVector* width = new G4PhysicsFreeVector(nEnergies); 108 G4int i; 106 G4int i; 109 for (i=0; i<nEnergies; i++) 107 for (i=0; i<nEnergies; i++) 110 { 108 { 111 G4double value = *(channelWidth+ i) * Ge 109 G4double value = *(channelWidth+ i) * GeV; 112 G4double e = energy[i]; 110 G4double e = energy[i]; 113 width->PutValue(i,e,value); 111 width->PutValue(i,e,value); 114 } 112 } 115 113 116 widths.push_back(width); 114 widths.push_back(width); 117 daughter1.push_back(name1); 115 daughter1.push_back(name1); 118 daughter2.push_back(name2); 116 daughter2.push_back(name2); 119 117 120 return; 118 return; 121 } 119 } 122 120 123 121 124 void G4PartialWidthTable::Dump() const 122 void G4PartialWidthTable::Dump() const 125 { 123 { 126 std::size_t entries = widths.size(); << 124 G4int entries = widths.size(); 127 125 128 for (std::size_t i=0; i<entries; ++i) << 126 G4int i; >> 127 for (i=0; i<entries; i++) 129 { 128 { 130 G4cout << " Channel " << i << ": " 129 G4cout << " Channel " << i << ": " 131 << daughter1[i] << " " << daughte << 130 << daughter1[i] << " " << daughter2[i] << G4endl; 132 G4PhysicsFreeVector* width = widths[i]; 131 G4PhysicsFreeVector* width = widths[i]; 133 for (G4int j=0; j<nEnergies; ++j) << 132 G4int j; >> 133 for (j=0; j<nEnergies; j++) 134 { 134 { 135 G4bool dummy = false; 135 G4bool dummy = false; 136 G4double e = energy[i]; 136 G4double e = energy[i]; 137 G4double w = width->GetValue(e,dummy); 137 G4double w = width->GetValue(e,dummy); 138 G4cout << j << ") Energy = " << e << ", Wi 138 G4cout << j << ") Energy = " << e << ", Width = " << w << G4endl; 139 } 139 } 140 } 140 } 141 return; 141 return; 142 } 142 } 143 143