Geant4 Cross Reference

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


  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 #include "globals.hh"                              25 #include "globals.hh"
 29 #include "G4ios.hh"                                26 #include "G4ios.hh"
 30 #include "G4KineticTrack.hh"                       27 #include "G4KineticTrack.hh"
 31 #include "G4XResonance.hh"                         28 #include "G4XResonance.hh"
 32 #include "Randomize.hh"                            29 #include "Randomize.hh"
 33 #include "G4Proton.hh"                             30 #include "G4Proton.hh"
 34                                                    31 
 35 G4XResonance::G4XResonance(const G4ParticleDef     32 G4XResonance::G4XResonance(const G4ParticleDefinition* in1, 
 36          const G4ParticleDefinition* in2,          33          const G4ParticleDefinition* in2,
 37          G4int iIsospinOut1, G4double spinOut1     34          G4int iIsospinOut1, G4double spinOut1, G4double massOut1,
 38          G4int iIsospinOut2, G4double spinOut2     35          G4int iIsospinOut2, G4double spinOut2, G4double massOut2,
 39          G4String subType1, G4String subType2,     36          G4String subType1, G4String subType2,
 40          const G4VXResonanceTable& sigmaTable)     37          const G4VXResonanceTable& sigmaTable) : 
 41   isoOut1(iIsospinOut1), iSpinOut1(spinOut1),      38   isoOut1(iIsospinOut1), iSpinOut1(spinOut1), mOut1(massOut1),
 42   isoOut2(iIsospinOut2), iSpinOut2(spinOut2),      39   isoOut2(iIsospinOut2), iSpinOut2(spinOut2), mOut2(massOut2)
 43                                                    40   
 44 {                                                  41 {
 45   table = sigmaTable.CrossSectionTable();          42   table = sigmaTable.CrossSectionTable();
 46   // Check if there is a valid cross section t     43   // Check if there is a valid cross section table for this channel
 47   if (table == 0)                                  44   if (table == 0)
 48     throw G4HadronicException(__FILE__, __LINE     45     throw G4HadronicException(__FILE__, __LINE__, "G4XResonance::G4XResonance - no cross section table available");
 49                                                    46   
 50   name = in1->GetParticleName() + in2->GetPart     47   name = in1->GetParticleName() + in2->GetParticleName() + " -> " + subType1 + subType2;
 51 }                                                  48 }
 52                                                    49 
 53                                                    50 
 54 G4XResonance::~G4XResonance()                      51 G4XResonance::~G4XResonance() 
 55 {                                              <<  52 { 
 56   delete table;                                    53   delete table;
 57   table = 0;                                       54   table = 0;
 58 }                                                  55 }
 59                                                    56 
 60                                                    57 
 61 G4bool G4XResonance::operator==(const G4XReson     58 G4bool G4XResonance::operator==(const G4XResonance &right) const
 62 {                                                  59 {
 63   return (this == (G4XResonance *) &right);        60   return (this == (G4XResonance *) &right);
 64 }                                                  61 }
 65                                                    62 
 66                                                    63 
 67 G4bool G4XResonance::operator!=(const G4XReson     64 G4bool G4XResonance::operator!=(const G4XResonance &right) const
 68 {                                                  65 {
 69   return (this != (G4XResonance *) &right);        66   return (this != (G4XResonance *) &right);
 70 }                                                  67 }
 71                                                    68 
 72                                                    69 
 73 G4String G4XResonance::Name() const                70 G4String G4XResonance::Name() const
 74 {                                                  71 {
 75   return name;                                     72   return name;
 76 }                                                  73 }
 77                                                    74 
 78 G4double G4XResonance::CrossSection(const G4Ki     75 G4double G4XResonance::CrossSection(const G4KineticTrack& trk1, const G4KineticTrack& trk2) const
 79 {                                                  76 {
 80   G4bool dummy = false;                            77   G4bool dummy = false;
 81   G4double sqrtS = (trk1.Get4Momentum() + trk2     78   G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
 82                                                    79 
 83   // pp -> trk1 + trk2 cross section               80   // pp -> trk1 + trk2 cross section
 84   G4double sigma = table->GetValue(sqrtS,dummy     81   G4double sigma = table->GetValue(sqrtS,dummy);
 85                                                    82   
 86   // Isospin correction                            83   // Isospin correction
 87   // from G4VXResonance;                           84   // from G4VXResonance; 
 88 //   G4cout << "SigmaCheck "<<sigma;               85 //   G4cout << "SigmaCheck "<<sigma;
 89   sigma *= IsospinCorrection(trk1,trk2,isoOut1     86   sigma *= IsospinCorrection(trk1,trk2,isoOut1,isoOut2,iSpinOut1,iSpinOut2); // from G4VXResonance
 90 //   G4cout << " "<<sigma<<G4endl;                 87 //   G4cout << " "<<sigma<<G4endl;
 91                                                    88 
 92   // Detailed balance                              89   // Detailed balance
 93   if (trk1.GetDefinition()->IsShortLived() ||      90   if (trk1.GetDefinition()->IsShortLived() || trk2.GetDefinition()->IsShortLived())
 94   {                                                91   {
 95    sigma *= DetailedBalance(trk1,trk2, isoOut1     92    sigma *= DetailedBalance(trk1,trk2, isoOut1,isoOut2, iSpinOut1,iSpinOut2, mOut1,mOut2); 
 96   }                                                93   }
 97                                                    94 
 98   return sigma;                                    95   return sigma;
 99                                                    96 
100 }                                                  97 }
101                                                    98 
102                                                    99