Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/inclxx/incl_physics/src/G4INCLOmegaNToPiNChannel.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/inclxx/incl_physics/src/G4INCLOmegaNToPiNChannel.cc (Version 11.3.0) and /processes/hadronic/models/inclxx/incl_physics/src/G4INCLOmegaNToPiNChannel.cc (Version 10.3.p1)


  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 // INCL++ intra-nuclear cascade model              26 // INCL++ intra-nuclear cascade model
 27 // Alain Boudard, CEA-Saclay, France               27 // Alain Boudard, CEA-Saclay, France
 28 // Joseph Cugnon, University of Liege, Belgium     28 // Joseph Cugnon, University of Liege, Belgium
 29 // Jean-Christophe David, CEA-Saclay, France       29 // Jean-Christophe David, CEA-Saclay, France
 30 // Pekka Kaitaniemi, CEA-Saclay, France, and H     30 // Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
 31 // Sylvie Leray, CEA-Saclay, France                31 // Sylvie Leray, CEA-Saclay, France
 32 // Davide Mancusi, CEA-Saclay, France              32 // Davide Mancusi, CEA-Saclay, France
 33 //                                                 33 //
 34 #define INCLXX_IN_GEANT4_MODE 1                    34 #define INCLXX_IN_GEANT4_MODE 1
 35                                                    35 
 36 #include "globals.hh"                              36 #include "globals.hh"
 37                                                    37 
 38 #include "G4INCLOmegaNToPiNChannel.hh"             38 #include "G4INCLOmegaNToPiNChannel.hh"
 39 #include "G4INCLKinematicsUtils.hh"                39 #include "G4INCLKinematicsUtils.hh"
 40 #include "G4INCLBinaryCollisionAvatar.hh"          40 #include "G4INCLBinaryCollisionAvatar.hh"
 41 #include "G4INCLRandom.hh"                         41 #include "G4INCLRandom.hh"
 42 #include "G4INCLGlobals.hh"                        42 #include "G4INCLGlobals.hh"
 43 #include "G4INCLLogger.hh"                         43 #include "G4INCLLogger.hh"
 44                                                    44 
 45 namespace G4INCL {                                 45 namespace G4INCL {
 46                                                    46 
 47     OmegaNToPiNChannel::OmegaNToPiNChannel(Par     47     OmegaNToPiNChannel::OmegaNToPiNChannel(Particle *p1, Particle *p2)
 48     : particle1(p1), particle2(p2)                 48     : particle1(p1), particle2(p2)
 49     {                                              49     {
 50                                                    50 
 51     }                                              51     }
 52                                                    52 
 53     OmegaNToPiNChannel::~OmegaNToPiNChannel(){     53     OmegaNToPiNChannel::~OmegaNToPiNChannel(){
 54                                                    54 
 55     }                                              55     }
 56                                                    56 
 57     void OmegaNToPiNChannel::fillFinalState(Fi     57     void OmegaNToPiNChannel::fillFinalState(FinalState *fs) {
 58         Particle * nucleon;                        58         Particle * nucleon;
 59         Particle * omega;                          59         Particle * omega;
 60         if(particle1->isNucleon()) {               60         if(particle1->isNucleon()) {
 61             nucleon = particle1;                   61             nucleon = particle1;
 62             omega = particle2;                     62             omega = particle2;
 63         } else {                                   63         } else {
 64             nucleon = particle2;                   64             nucleon = particle2;
 65             omega = particle1;                     65             omega = particle1;
 66         }                                          66         }
 67                                                    67 
 68     const G4double r2 = Random::shoot();           68     const G4double r2 = Random::shoot();
 69     if (nucleon->getType() == Neutron) {           69     if (nucleon->getType() == Neutron) {
 70       if (r2*3. < 2.) {                            70       if (r2*3. < 2.) {
 71         nucleon->setType(Proton);                  71         nucleon->setType(Proton);
 72         omega->setType(PiMinus);                   72         omega->setType(PiMinus);
 73       }                                            73       }
 74       else {                                       74       else {
 75         nucleon->setType(Neutron);                 75         nucleon->setType(Neutron);
 76         omega->setType(PiZero);                    76         omega->setType(PiZero);
 77       }                                            77       }
 78     }                                              78     }
 79     else {                                         79     else {
 80       if (r2*3. < 2.) {                            80       if (r2*3. < 2.) {
 81         nucleon->setType(Neutron);                 81         nucleon->setType(Neutron);
 82         omega->setType(PiPlus);                    82         omega->setType(PiPlus);
 83       }                                            83       }
 84       else {                                       84       else {
 85         nucleon->setType(Proton);                  85         nucleon->setType(Proton);
 86         omega->setType(PiZero);                    86         omega->setType(PiZero);
 87       }                                            87       }
 88     }                                              88     }
 89                                                    89     
 90         G4double sh=nucleon->getEnergy()+omega     90         G4double sh=nucleon->getEnergy()+omega->getEnergy();
 91     G4double mn=nucleon->getMass();                91     G4double mn=nucleon->getMass();
 92     G4double me=omega->getMass();                  92     G4double me=omega->getMass();
 93     G4double en=(sh*sh+mn*mn-me*me)/(2*sh);        93     G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
 94     nucleon->setEnergy(en);                        94     nucleon->setEnergy(en);
 95     G4double ee=std::sqrt(en*en-mn*mn+me*me);      95     G4double ee=std::sqrt(en*en-mn*mn+me*me);
 96     omega->setEnergy(ee);                          96     omega->setEnergy(ee);
 97     G4double pn=std::sqrt(en*en-mn*mn);            97     G4double pn=std::sqrt(en*en-mn*mn);
 98 /*// test isotropy                                 98 /*// test isotropy
 99     const G4double pi=std::acos(-1.0);         <<  99     const G4double pi=std::acos(-1.0);    
100     G4double x1;                                  100     G4double x1;
101     G4double u1;                                  101     G4double u1;
102     G4double fteta;                               102     G4double fteta;
103     G4double teta;                                103     G4double teta;
104     G4double fi;                                  104     G4double fi;
105                                                   105 
106     G4int passe1=0;                               106     G4int passe1=0;
107     while (passe1==0) {                           107     while (passe1==0) {
108       // Sample x from 0 to pi/2                  108       // Sample x from 0 to pi/2
109       x1=(pi/2.)*Random::shoot();                 109       x1=(pi/2.)*Random::shoot();
110       // Sample u from 0 to 1                     110       // Sample u from 0 to 1
111       u1=Random::shoot();                         111       u1=Random::shoot();
112       fteta=std::sin(x1);                         112       fteta=std::sin(x1);
113       // The condition                            113       // The condition
114       if (u1 < fteta) {                           114       if (u1 < fteta) {
115         teta=x1;                                  115         teta=x1;
116         passe1=1;                                 116         passe1=1;
117           if (Random::shoot() < 0.5) {            117           if (Random::shoot() < 0.5) {
118           teta=pi-teta;                           118           teta=pi-teta;
119         }                                         119         }
120       }                                           120       }
121       }                                           121       }
122         fi=(2.0*pi)*Random::shoot();              122         fi=(2.0*pi)*Random::shoot();
123                                                   123     
124 // end test                                       124 // end test   
125         ThreeVector mom_nucleon(                  125         ThreeVector mom_nucleon(
126                                 pn*std::sin(te    126                                 pn*std::sin(teta)*std::cos(fi),
127                                 pn*std::sin(te    127                                 pn*std::sin(teta)*std::sin(fi),
128                                 pn*std::cos(te    128                                 pn*std::cos(teta)
129                                 );                129                                 );
130     if (Random::shoot() < 0.5) {                  130     if (Random::shoot() < 0.5) {
131       nucleon->setMomentum(mom_nucleon);          131       nucleon->setMomentum(mom_nucleon);
132       omega->setMomentum(-mom_nucleon);           132       omega->setMomentum(-mom_nucleon);
133     }                                             133     }
134     else {                                        134     else {
135       omega->setMomentum(mom_nucleon);            135       omega->setMomentum(mom_nucleon);
136         nucleon->setMomentum(-mom_nucleon);       136         nucleon->setMomentum(-mom_nucleon);
137     }                                             137     }
138 */                                                138 */
139     ThreeVector mom_nucleon = Random::normVect    139     ThreeVector mom_nucleon = Random::normVector(pn);
140                                                   140     
141     nucleon->setMomentum(mom_nucleon);            141     nucleon->setMomentum(mom_nucleon);
142     omega->setMomentum(-mom_nucleon);             142     omega->setMomentum(-mom_nucleon);
143                                                   143         
144         fs->addModifiedParticle(nucleon);         144         fs->addModifiedParticle(nucleon);
145         fs->addModifiedParticle(omega);           145         fs->addModifiedParticle(omega);
146     }                                             146     }
147                                                   147 
148 }                                                 148 }
149                                                   149