Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/im_r_matrix/src/G4AngularDistributionPP.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/im_r_matrix/src/G4AngularDistributionPP.cc (Version 11.3.0) and /processes/hadronic/models/im_r_matrix/src/G4AngularDistributionPP.cc (Version 7.1)


  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 #include "G4AngularDistributionPP.hh"              23 #include "G4AngularDistributionPP.hh"
 27 #include "G4PhysicalConstants.hh"              << 
 28 #include "G4SystemOfUnits.hh"                  << 
 29 #include "Randomize.hh"                            24 #include "Randomize.hh"
 30 #include "G4ios.hh"                                25 #include "G4ios.hh"
 31                                                    26 
 32 // Initialization of static data arrays:           27 // Initialization of static data arrays:
 33 #include "G4AngularDistributionPPData.hh"          28 #include "G4AngularDistributionPPData.hh"
 34 #include "Randomize.hh"                            29 #include "Randomize.hh"
 35                                                    30 
 36                                                    31 
 37 G4double G4AngularDistributionPP::CosTheta(G4d <<  32 G4double G4AngularDistributionPP::CosTheta(G4double s, G4double m1, G4double m2) const
 38 {                                                  33 {
 39     G4int verboseLevel=1;                          34     G4int verboseLevel=1;
 40                                                    35 
 41     G4double ek= ((S - sqr(m_1) -sqr(m_2) )/(2 <<  36     G4double ek= ((s - sqr(m1) -sqr(m2) )/(2*m1) - m1 )/GeV   ;    // kinetic energy in GeV
 42                                                    37 
 43     // Find energy bin                             38     // Find energy bin
 44                                                    39 
 45     G4int je1 = 0;                                 40     G4int je1 = 0;
 46     G4int je2 = NENERGY - 1;                       41     G4int je2 = NENERGY - 1;
 47    G4int iterationsLeft=2*NENERGY +1;          << 
 48     do {                                           42     do {
 49       G4int midBin = (je1 + je2)/2;                43       G4int midBin = (je1 + je2)/2;
 50       if (ek < elab[midBin])                       44       if (ek < elab[midBin])
 51         je2 = midBin;                              45         je2 = midBin;
 52       else                                         46       else
 53         je1 = midBin;                              47         je1 = midBin;
 54     } while ( (je2 - je1) > 1 && --iterationsL <<  48     } while (je2 - je1 > 1);
 55    if ( iterationsLeft <= 0 ) {                << 
 56        G4Exception("G4AngularDistributionPP",  << 
 57              "Problem with energy bin (elab) d << 
 58                                                << 
 59    }                                           << 
 60     //    G4int j;                                 49     //    G4int j;
 61     //std::abs(ek-elab[je1]) < std::abs(ek-ela     50     //std::abs(ek-elab[je1]) < std::abs(ek-elab[je2]) ? j = je1 : j = je2;
 62     G4double delab = elab[je2] - elab[je1];        51     G4double delab = elab[je2] - elab[je1];
 63                                                    52 
 64     // Sample the angle                            53     // Sample the angle
 65                                                    54 
 66     G4double sample = G4UniformRand();         <<  55     G4float sample = G4UniformRand();
 67     G4int ke1 = 0;                                 56     G4int ke1 = 0;
 68     G4int ke2 = NANGLE - 1;                        57     G4int ke2 = NANGLE - 1;
 69     G4double dsig = sig[je2][0] - sig[je1][0];     58     G4double dsig = sig[je2][0] - sig[je1][0];
 70     G4double rc = dsig/delab;                      59     G4double rc = dsig/delab;
 71     G4double b = sig[je1][0] - rc*elab[je1];       60     G4double b = sig[je1][0] - rc*elab[je1];
 72     G4double sigint1 = rc*ek + b;                  61     G4double sigint1 = rc*ek + b;
 73     G4double sigint2 = 0.;                         62     G4double sigint2 = 0.;
 74                                                    63 
 75     if (verboseLevel > 1) G4cout << "sample="      64     if (verboseLevel > 1) G4cout << "sample=" << sample << G4endl
 76                                  << ek << " "      65                                  << ek << " " << ke1 << " " << ke2 << " "
 77                                  << sigint1 <<     66                                  << sigint1 << " " << sigint2 << G4endl;
 78                                                    67 
 79     iterationsLeft= 2*NANGLE +1;               <<  68     do {
 80      do {                                      << 
 81       G4int midBin = (ke1 + ke2)/2;                69       G4int midBin = (ke1 + ke2)/2;
 82       dsig = sig[je2][midBin] - sig[je1][midBi     70       dsig = sig[je2][midBin] - sig[je1][midBin];
 83       rc = dsig/delab;                             71       rc = dsig/delab;
 84       b = sig[je1][midBin] - rc*elab[je1];         72       b = sig[je1][midBin] - rc*elab[je1];
 85       G4double sigint = rc*ek + b;                 73       G4double sigint = rc*ek + b;
 86       if (sample < sigint) {                       74       if (sample < sigint) {
 87         ke2 = midBin;                              75         ke2 = midBin;
 88         sigint2 = sigint;                          76         sigint2 = sigint;
 89       }                                            77       }
 90       else {                                       78       else {
 91         ke1 = midBin;                              79         ke1 = midBin;
 92         sigint1 = sigint;                          80         sigint1 = sigint;
 93       }                                            81       }
 94       if (verboseLevel > 1)G4cout << ke1 << "      82       if (verboseLevel > 1)G4cout << ke1 << " " << ke2 << " "
 95                                   << sigint1 <     83                                   << sigint1 << " " << sigint2 << G4endl;
 96     } while ( (ke2 - ke1) > 1 && --iterationsL <<  84     } while (ke2 - ke1 > 1);
 97    if ( iterationsLeft <= 0 ) {                << 
 98        G4Exception("G4AngularDistributionPP",  << 
 99              "Problem with angular distributio << 
100    }                                           << 
101                                                    85 
102     // sigint1 and sigint2 should be recoverab     86     // sigint1 and sigint2 should be recoverable from above loop
103                                                    87 
104     //    G4double dsig = sig[je2][ke1] - sig[     88     //    G4double dsig = sig[je2][ke1] - sig[je1][ke1];
105     //    G4double rc = dsig/delab;                89     //    G4double rc = dsig/delab;
106     //    G4double b = sig[je1][ke1] - rc*elab     90     //    G4double b = sig[je1][ke1] - rc*elab[je1];
107     //    G4double sigint1 = rc*ek + b;            91     //    G4double sigint1 = rc*ek + b;
108                                                    92 
109     //    G4double dsig = sig[je2][ke2] - sig[     93     //    G4double dsig = sig[je2][ke2] - sig[je1][ke2];
110     //    G4double rc = dsig/delab;                94     //    G4double rc = dsig/delab;
111     //    G4double b = sig[je1][ke2] - rc*elab     95     //    G4double b = sig[je1][ke2] - rc*elab[je1];
112     //    G4double sigint2 = rc*ek + b;            96     //    G4double sigint2 = rc*ek + b;
113                                                    97 
114     dsig = sigint2 - sigint1;                      98     dsig = sigint2 - sigint1;
115     rc = 1./dsig;                                  99     rc = 1./dsig;
116     b = ke1 - rc*sigint1;                         100     b = ke1 - rc*sigint1;
117     G4double kint = rc*sample + b;                101     G4double kint = rc*sample + b;
118     G4double theta = (0.5 + kint)*pi/180.;        102     G4double theta = (0.5 + kint)*pi/180.;
119                                                   103 
120     //    G4int k;                                104     //    G4int k;
121     //std::abs(sample-sig[j][ke1]) < std::abs(    105     //std::abs(sample-sig[j][ke1]) < std::abs(sample-sig[j][ke2]) ? k = ke1 : k = ke2;
122     //    G4double theta = (0.5 + k)*pi/180.;     106     //    G4double theta = (0.5 + k)*pi/180.;
123                                                   107 
124     if (verboseLevel > 1) {                       108     if (verboseLevel > 1) {
125       G4cout << "   energy bin " << je1 << " e    109       G4cout << "   energy bin " << je1 << " energy=" << elab[je1] << G4endl;
126       G4cout << "   angle bin " << kint << " a    110       G4cout << "   angle bin " << kint << " angle=" << theta/degree << G4endl;
127     }                                             111     }
128     G4double costh= std::cos(theta);              112     G4double costh= std::cos(theta);
129 //    G4cout << "scattering angle hpw "<<costh    113 //    G4cout << "scattering angle hpw "<<costh<<" "<<s<<" "<<m1<<" "<<m2<<G4endl;
130     return costh;                                 114     return costh;
131 }                                                 115 }
132                                                   116 
133 G4double G4AngularDistributionPP::Phi() const     117 G4double G4AngularDistributionPP::Phi() const
134 {                                                 118 {
135     return twopi * G4UniformRand();               119     return twopi * G4UniformRand();
136 }                                                 120 }
137                                                   121