Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/g3tog4/src/G4gsmixt.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /g3tog4/src/G4gsmixt.cc (Version 11.3.0) and /g3tog4/src/G4gsmixt.cc (Version 6.0)


  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 //
                                                   >>  24 // $Id: G4gsmixt.cc,v 1.10 2003/06/16 16:50:56 gunter Exp $
                                                   >>  25 // GEANT4 tag $Name: geant4-06-00 $
 27 //                                                 26 //
 28 // by I.Hrivnacova, 27 Sep 99                      27 // by I.Hrivnacova, 27 Sep 99
 29                                                    28 
                                                   >>  29 #include "globals.hh"
 30 #include <iomanip>                                 30 #include <iomanip>
                                                   >>  31 #include <strstream>
 31 #include <iomanip>                                 32 #include <iomanip>
 32                                                    33 
 33 #include "globals.hh"                          << 
 34 #include "G4SystemOfUnits.hh"                  << 
 35                                                << 
 36 #include "G3toG4.hh"                               34 #include "G3toG4.hh"
 37 #include "G3EleTable.hh"                           35 #include "G3EleTable.hh"
 38 #include "G3MatTable.hh"                           36 #include "G3MatTable.hh"
 39 #include "G4Material.hh"                           37 #include "G4Material.hh"
 40 #include "G4Isotope.hh"                            38 #include "G4Isotope.hh"
 41                                                    39 
 42 void PG4gsmixt(G4String *tokens)               <<  40 void PG4gsmixt(G4String tokens[])
 43 {                                                  41 {
 44     // fill the parameter containers               42     // fill the parameter containers
 45     G3fillParams(tokens,PTgsmixt);                 43     G3fillParams(tokens,PTgsmixt);
 46                                                    44 
 47     // interpret the parameters                    45     // interpret the parameters
 48     G4String name = Spar[0].data();                46     G4String name = Spar[0].data();
 49     G4int imate = Ipar[0];                         47     G4int imate = Ipar[0];
 50     G4int nlmat = Ipar[1];                         48     G4int nlmat = Ipar[1];
 51     //G4double dens = Rpar[0]*g/cm3;               49     //G4double dens = Rpar[0]*g/cm3;
 52     G4double dens = Rpar[0];                       50     G4double dens = Rpar[0];
 53     G4double *a = Rpar + 1;                        51     G4double *a = Rpar + 1;
 54     G4double *z = Rpar + 1+std::abs(nlmat);    <<  52     G4double *z = Rpar + 1+abs(nlmat);
 55     G4double *wmat = Rpar + 1 + 2*std::abs(nlm <<  53     G4double *wmat = Rpar + 1 + 2*abs(nlmat);
 56 /*                                             <<  54 
 57     for (int i=0; i<std::abs(nlmat); i++){     <<  55     for (int i=0; i<abs(nlmat); i++){
 58       //Rpar[i]=Rpar[i]*g/mole;                    56       //Rpar[i]=Rpar[i]*g/mole;
 59       Rpar[i]=Rpar[i];                             57       Rpar[i]=Rpar[i];
 60     };                                             58     };
 61 */                                             << 
 62     G4gsmixt(imate,name,a,z,dens,nlmat,wmat);      59     G4gsmixt(imate,name,a,z,dens,nlmat,wmat);
 63 }                                                  60 }
 64                                                    61 
 65 // replaced with G3EleTable                        62 // replaced with G3EleTable 
 66 // only used G4Elements are created;               63 // only used G4Elements are created;
 67 // !! no checking of given A of the element;       64 // !! no checking of given A of the element;
 68 //                                                 65 //
 69 // extern G4Element* CreateElement(G4double ze     66 // extern G4Element* CreateElement(G4double zeff, G4double aeff, G4String matName);
 70                                                    67 
 71                                                    68 
 72 void G4gsmixt(G4int imate, G4String name, G4do     69 void G4gsmixt(G4int imate, G4String name, G4double* a, G4double* z,
 73               G4double dens, G4int nlmat, G4do     70               G4double dens, G4int nlmat, G4double* wmat)
 74 {                                                  71 {
 75   // in Geant3:                                    72   // in Geant3:
 76   // After a call with ratios by number (negat     73   // After a call with ratios by number (negative number of elements), 
 77   // the ratio array is changed to the ratio b     74   // the ratio array is changed to the ratio by weight, so all successive 
 78   // calls with the same array must specify th     75   // calls with the same array must specify the number of elements as 
 79   // positive                                      76   // positive 
 80   G4int i=0;                                       77   G4int i=0;
 81   if (nlmat<0) {                                   78   if (nlmat<0) {
 82     // in case of proportions given in atom co     79     // in case of proportions given in atom counts (nlmat<0),
 83     // the wmat[i] are converted to weight fra     80     // the wmat[i] are converted to weight fractions
 84     G4double aMol = 0.;                            81     G4double aMol = 0.;
 85     for (i=0; i<std::abs(nlmat); i++) {        <<  82     for (i=0; i<abs(nlmat); i++) { 
 86       // total molecular weight                    83       // total molecular weight 
 87       aMol += wmat[i]*a[i];                        84       aMol += wmat[i]*a[i]; 
 88     }                                              85     }  
 89     if (aMol == 0.) {                          <<  86     if (aMol == 0.)
 90       G4String text = "G4mixt: Total molecular <<  87       G4Exception("\nG4mixt: Total molecular weight in " +
 91       G4Exception("G4gsmixt()", "G3toG40016",  <<  88         name + " = 0.");       
 92       return;                                  <<  89     for (i=0; i<abs(nlmat); i++) {
 93     }                                          << 
 94     for (i=0; i<std::abs(nlmat); i++) {        << 
 95       // weight fractions                          90       // weight fractions
 96       wmat[i] = wmat[i]*a[i]/aMol;                 91       wmat[i] = wmat[i]*a[i]/aMol;
 97     }                                              92     }
 98   }                                                93   }
 99                                                    94 
100   // create material with given number of comp     95   // create material with given number of components
101   // (elements)                                    96   // (elements)
102                                                    97 
103   G4Material* material                             98   G4Material* material 
104     = new G4Material(name, dens*g/cm3, std::ab <<  99     = new G4Material(name, dens*g/cm3, abs(nlmat));
105   for (i=0; i<std::abs(nlmat); i++) {          << 100   for (i=0; i<abs(nlmat); i++) {
106     // add units                                  101     // add units
107     // G4Element* element = G4Element(z[i], a[    102     // G4Element* element = G4Element(z[i], a[i]*g/mole, name);
108     G4Element* element = G3Ele.GetEle(z[i]);      103     G4Element* element = G3Ele.GetEle(z[i]);
109     material->AddElement(element, wmat[i]);       104     material->AddElement(element, wmat[i]);    
110   }                                               105   }
111                                                   106 
112   // add the material to the List                 107   // add the material to the List
113   G3Mat.put(imate, material);                     108   G3Mat.put(imate, material);
114 }                                                 109 }
115                                                   110 
116 /*                                                111 /*
117 void G4gsmixt(G4int imate, G4String name, G4do    112 void G4gsmixt(G4int imate, G4String name, G4double a[], G4double z[],
118               G4double dens, G4int nlmat, G4do    113               G4double dens, G4int nlmat, G4double wmat[]){
119   G4int nmate = std::abs(nlmat);               << 114   G4int nmate = abs(nlmat);
120   G4String sname = name.strip(G4String::both);    115   G4String sname = name.strip(G4String::both);
121   G4double theDensity = dens*g/cm3;               116   G4double theDensity = dens*g/cm3;
122                                                   117 
123   G4Material* theMixture = new G4Material(name    118   G4Material* theMixture = new G4Material(name, dens, nmate); 
124   G4bool ok=true;                                 119   G4bool ok=true;
125   for (int i=0; i< nmate; i++){                   120   for (int i=0; i< nmate; i++){
126     G4Element* theElement = G3Ele.GetEle(z[i])    121     G4Element* theElement = G3Ele.GetEle(z[i]);
127     if (nlmat>0) {                                122     if (nlmat>0) {
128       G4double fractionmass = wmat[i];            123       G4double fractionmass = wmat[i];
129       ok = ok && std::abs(fractionmass)<=1.;   << 124       ok = ok && abs(fractionmass)<=1.;
130       theMixture->AddElement(theElement, fract    125       theMixture->AddElement(theElement, fractionmass);
131     } else if (nlmat<0) {                         126     } else if (nlmat<0) {
132       G4int natoms = wmat[i];                     127       G4int natoms = wmat[i];
133       ok = ok && wmat[i] == natoms;               128       ok = ok && wmat[i] == natoms;
134       theMixture->AddElement(theElement, natom    129       theMixture->AddElement(theElement, natoms);
135     } else {                                      130     } else {
136       ok=false;                                   131       ok=false;
137     }                                             132     }
138   }                                               133   }
139   if (ok) {                                       134   if (ok) {
140     G3Mat.put(imate, theMixture);                 135     G3Mat.put(imate, theMixture);
141   } else {                                        136   } else {
142     if (nlmat>0) {                                137     if (nlmat>0) {
143       G4cerr << "G4gsmixt: for mixture '" << n    138       G4cerr << "G4gsmixt: for mixture '" << name 
144        << "' some |weights|>1 : " << G4endl;      139        << "' some |weights|>1 : " << G4endl;
145       for (G4int i=0;i<nlmat; i++) {              140       for (G4int i=0;i<nlmat; i++) {
146   G4cerr << "Component " << std::setw(3) << i+    141   G4cerr << "Component " << std::setw(3) << i+1 << " fraction: "
147          << std::setw(10) << wmat[i] << G4endl    142          << std::setw(10) << wmat[i] << G4endl;
148       }                                           143       }
149     } else if (nlmat<0) {                         144     } else if (nlmat<0) {
150       G4cerr << "G4gsmixt: for mixture '" << n    145       G4cerr << "G4gsmixt: for mixture '" << name 
151        << "' some #natoms are non-integer: " <    146        << "' some #natoms are non-integer: " << G4endl;
152       for (G4int i=0;i<nlmat; i++) {              147       for (G4int i=0;i<nlmat; i++) {
153   G4cerr << "Component " << std::setw(3) << i+    148   G4cerr << "Component " << std::setw(3) << i+1 << " #atoms "
154          << std::setw(10) << wmat[i] << G4endl    149          << std::setw(10) << wmat[i] << G4endl;
155       }                                           150       }
156     } else {                                      151     } else {
157       G4cerr << "G4gsmixt: Number of component    152       G4cerr << "G4gsmixt: Number of components for mixture '" 
158        << name << "' (" << nlmat << ") not all    153        << name << "' (" << nlmat << ") not allowed." << G4endl;
159     }                                             154     }
160   }                                               155   }
161 }                                                 156 }
162 */                                                157 */
163                                                   158 
164                                                   159 
165                                                   160 
166                                                   161