Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroChemicalPotential.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 /processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroChemicalPotential.cc (Version 11.3.0) and /processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroChemicalPotential.cc (Version 10.3.p3)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
                                                   >>  27 // $Id: G4StatMFMacroChemicalPotential.cc 100379 2016-10-19 15:05:35Z gcosmo $
 27 //                                                 28 //
 28 // Hadronic Process: Nuclear De-excitations        29 // Hadronic Process: Nuclear De-excitations
 29 // by V. Lara                                      30 // by V. Lara
 30                                                    31 
 31 #include "G4StatMFMacroChemicalPotential.hh"       32 #include "G4StatMFMacroChemicalPotential.hh"
 32 #include "G4PhysicalConstants.hh"                  33 #include "G4PhysicalConstants.hh"
 33 #include "G4Pow.hh"                                34 #include "G4Pow.hh"
 34                                                    35 
 35 // operators definitions                           36 // operators definitions
 36 G4StatMFMacroChemicalPotential &                   37 G4StatMFMacroChemicalPotential & 
 37 G4StatMFMacroChemicalPotential::operator=(cons     38 G4StatMFMacroChemicalPotential::operator=(const G4StatMFMacroChemicalPotential & ) 
 38 {                                                  39 {
 39     throw G4HadronicException(__FILE__, __LINE <<  40     throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::operator= meant to not be accessable");
 40     return *this;                                  41     return *this;
 41 }                                                  42 }
 42                                                    43 
 43 G4bool G4StatMFMacroChemicalPotential::operato     44 G4bool G4StatMFMacroChemicalPotential::operator==(const G4StatMFMacroChemicalPotential & ) const 
 44 {                                                  45 {
 45     throw G4HadronicException(__FILE__, __LINE <<  46     throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::operator== meant to not be accessable");
 46     return false;                                  47     return false;
 47 }                                                  48 }
 48                                                    49 
 49                                                    50 
 50 G4bool G4StatMFMacroChemicalPotential::operato     51 G4bool G4StatMFMacroChemicalPotential::operator!=(const G4StatMFMacroChemicalPotential & ) const 
 51 {                                                  52 {
 52     throw G4HadronicException(__FILE__, __LINE <<  53     throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::operator!= meant to not be accessable");
 53     return true;                                   54     return true;
 54 }                                                  55 }
 55                                                    56 
 56 G4double G4StatMFMacroChemicalPotential::CalcC     57 G4double G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu(void) 
 57 //  Calculate Chemical potential \nu               58 //  Calculate Chemical potential \nu
 58 {                                                  59 {
 59   G4Pow* g4calc = G4Pow::GetInstance();            60   G4Pow* g4calc = G4Pow::GetInstance();
 60   G4double CP = G4StatMFParameters::GetCoulomb     61   G4double CP = G4StatMFParameters::GetCoulomb();
 61                                                    62 
 62   // Initial value for _ChemPotentialNu            63   // Initial value for _ChemPotentialNu 
 63   _ChemPotentialNu = (theZ/theA)*(8.0*G4StatMF     64   _ChemPotentialNu = (theZ/theA)*(8.0*G4StatMFParameters::GetGamma0()
 64           +2.0*CP*g4calc->Z23(theA))               65           +2.0*CP*g4calc->Z23(theA))
 65     - 4.0*G4StatMFParameters::GetGamma0();         66     - 4.0*G4StatMFParameters::GetGamma0();
 66                                                    67     
 67   G4double ChemPa = _ChemPotentialNu;              68   G4double ChemPa = _ChemPotentialNu;
 68   G4double ChemPb = 0.5*_ChemPotentialNu;          69   G4double ChemPb = 0.5*_ChemPotentialNu;
 69                                                    70     
 70   G4double fChemPa = this->operator()(ChemPa);     71   G4double fChemPa = this->operator()(ChemPa); 
 71   G4double fChemPb = this->operator()(ChemPb);     72   G4double fChemPb = this->operator()(ChemPb); 
 72                                                    73 
 73   if (fChemPa*fChemPb > 0.0) {                     74   if (fChemPa*fChemPb > 0.0) {    
 74     // bracketing the solution                     75     // bracketing the solution
 75     if (fChemPa < 0.0) {                           76     if (fChemPa < 0.0) {
 76       do {                                         77       do {
 77   ChemPb -= 1.5*std::abs(ChemPb-ChemPa);           78   ChemPb -= 1.5*std::abs(ChemPb-ChemPa);
 78   fChemPb = this->operator()(ChemPb);              79   fChemPb = this->operator()(ChemPb);   
 79   // Loop checking, 05-Aug-2015, Vladimir Ivan     80   // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
 80       } while (fChemPb < 0.0);                     81       } while (fChemPb < 0.0);
 81     } else {                                       82     } else {
 82       do {                                         83       do {
 83   ChemPb += 1.5*std::abs(ChemPb-ChemPa);           84   ChemPb += 1.5*std::abs(ChemPb-ChemPa);
 84   fChemPb = this->operator()(ChemPb);              85   fChemPb = this->operator()(ChemPb);
 85   // Loop checking, 05-Aug-2015, Vladimir Ivan     86   // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
 86       } while (fChemPb > 0.0);                     87       } while (fChemPb > 0.0);
 87     }                                              88     }
 88   }                                                89   }
 89                                                    90 
 90   G4Solver<G4StatMFMacroChemicalPotential> * t     91   G4Solver<G4StatMFMacroChemicalPotential> * theSolver =
 91     new G4Solver<G4StatMFMacroChemicalPotentia     92     new G4Solver<G4StatMFMacroChemicalPotential>(100,1.e-4);
 92   theSolver->SetIntervalLimits(ChemPa,ChemPb);     93   theSolver->SetIntervalLimits(ChemPa,ChemPb);
 93   //    if (!theSolver->Crenshaw(*this))           94   //    if (!theSolver->Crenshaw(*this)) 
 94   if (!theSolver->Brent(*this)){                   95   if (!theSolver->Brent(*this)){
 95     G4cout <<"G4StatMFMacroChemicalPotential:"     96     G4cout <<"G4StatMFMacroChemicalPotential:"<<" ChemPa="<<ChemPa
 96      <<" ChemPb="<<ChemPb<< G4endl;                97      <<" ChemPb="<<ChemPb<< G4endl;
 97     G4cout <<"G4StatMFMacroChemicalPotential:"     98     G4cout <<"G4StatMFMacroChemicalPotential:"<<" fChemPa="<<fChemPa
 98      <<" fChemPb="<<fChemPb<< G4endl;              99      <<" fChemPb="<<fChemPb<< G4endl;
 99     throw G4HadronicException(__FILE__, __LINE    100     throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu: I couldn't find the root.");
100   }                                               101   }
101   _ChemPotentialNu = theSolver->GetRoot();        102   _ChemPotentialNu = theSolver->GetRoot();
102   delete theSolver;                               103   delete theSolver;
103   return _ChemPotentialNu;                        104   return _ChemPotentialNu;
104 }                                                 105 }
105                                                   106 
106                                                   107 
107                                                   108 
108 G4double G4StatMFMacroChemicalPotential::CalcM    109 G4double G4StatMFMacroChemicalPotential::CalcMeanZ(const G4double nu)
109 {                                                 110 {
110   std::vector<G4VStatMFMacroCluster*>::iterato    111   std::vector<G4VStatMFMacroCluster*>::iterator i;
111   for (i= _theClusters->begin()+1; i != _theCl    112   for (i= _theClusters->begin()+1; i != _theClusters->end(); ++i) 
112     {                                             113     { 
113       (*i)->CalcZARatio(nu);                      114       (*i)->CalcZARatio(nu);
114     }                                             115     }
115   CalcChemicalPotentialMu(nu);                    116   CalcChemicalPotentialMu(nu);
116   // This is important, the Z over A ratio for    117   // This is important, the Z over A ratio for proton and neutron depends on the 
117   // chemical potential Mu, while for the firs    118   // chemical potential Mu, while for the first guess for Chemical potential mu 
118   // some values of Z over A ratio. This is th    119   // some values of Z over A ratio. This is the reason for that.
119   (*_theClusters->begin())->CalcZARatio(nu);      120   (*_theClusters->begin())->CalcZARatio(nu);
120                                                   121   
121   G4double MeanZ = 0.0;                           122   G4double MeanZ = 0.0;
122   G4int n = 1;                                    123   G4int n = 1;
123   for (i = _theClusters->begin(); i != _theClu    124   for (i = _theClusters->begin(); i != _theClusters->end(); ++i)
124     {                                             125     {
125       MeanZ += (n++) * (*i)->GetZARatio() * (*    126       MeanZ += (n++) * (*i)->GetZARatio() * (*i)->GetMeanMultiplicity(); 
126     }                                             127     }
127   return MeanZ;                                   128   return MeanZ;
128 }                                                 129 }
129                                                   130 
130 void G4StatMFMacroChemicalPotential::CalcChemi    131 void G4StatMFMacroChemicalPotential::CalcChemicalPotentialMu(const G4double nu)
131 //  Calculate Chemical potential \mu              132 //  Calculate Chemical potential \mu
132 // For that is necesary to calculate mean mult    133 // For that is necesary to calculate mean multiplicities
133 {                                                 134 {
134   G4StatMFMacroMultiplicity * theMultip = new     135   G4StatMFMacroMultiplicity * theMultip = new
135     G4StatMFMacroMultiplicity(theA,_Kappa,_Mea    136     G4StatMFMacroMultiplicity(theA,_Kappa,_MeanTemperature,nu,_theClusters);
136                                                   137   
137   _ChemPotentialMu = theMultip->CalcChemicalPo    138   _ChemPotentialMu = theMultip->CalcChemicalPotentialMu();
138   _MeanMultiplicity = theMultip->GetMeanMultip    139   _MeanMultiplicity = theMultip->GetMeanMultiplicity();
139                                                   140   
140   delete theMultip;                               141   delete theMultip;
141                                                   142   
142   return;                                         143   return; 
143 }                                                 144 }
144                                                   145