Geant4 Cross Reference

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


  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 "G4INCLNpiToNKKbChannel.hh"               38 #include "G4INCLNpiToNKKbChannel.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 #include <algorithm>                               44 #include <algorithm>
 45 #include "G4INCLPhaseSpaceGenerator.hh"            45 #include "G4INCLPhaseSpaceGenerator.hh"
 46                                                    46 
 47 namespace G4INCL {                                 47 namespace G4INCL {
 48                                                    48   
 49   const G4double NpiToNKKbChannel::angularSlop     49   const G4double NpiToNKKbChannel::angularSlope = 6.; // What is the exact effect? Sould be check
 50                                                    50   
 51   NpiToNKKbChannel::NpiToNKKbChannel(Particle      51   NpiToNKKbChannel::NpiToNKKbChannel(Particle *p1, Particle *p2)
 52     : particle1(p1), particle2(p2)                 52     : particle1(p1), particle2(p2)
 53     {}                                             53     {}
 54                                                    54   
 55   NpiToNKKbChannel::~NpiToNKKbChannel(){}          55   NpiToNKKbChannel::~NpiToNKKbChannel(){}
 56                                                    56   
 57   void NpiToNKKbChannel::fillFinalState(FinalS     57   void NpiToNKKbChannel::fillFinalState(FinalState *fs) {
 58                                                    58     
 59         // p pi+ -> p K+ K0b (1/2)                 59         // p pi+ -> p K+ K0b (1/2)
 60                                                    60         
 61         // p pi0 -> p K0 K0b (1/4)                 61         // p pi0 -> p K0 K0b (1/4)
 62         // p pi0 -> p K+ K- (1/4)                  62         // p pi0 -> p K+ K- (1/4)
 63         // p pi0 -> n K+ K0b (1)                   63         // p pi0 -> n K+ K0b (1)
 64                                                    64         
 65         // p pi- -> p K0 K- (1/2)                  65         // p pi- -> p K0 K- (1/2)
 66         // p pi- -> n K+ K- (1)                    66         // p pi- -> n K+ K- (1)
 67         // p pi- -> n K0 K0b (1)                   67         // p pi- -> n K0 K0b (1)
 68                                                    68     
 69     Particle *nucleon;                             69     Particle *nucleon;
 70     Particle *pion;                                70     Particle *pion;
 71                                                    71     
 72                                                    72     
 73     if(particle1->isNucleon()){                    73     if(particle1->isNucleon()){
 74       nucleon = particle1;                         74       nucleon = particle1;
 75       pion = particle2;                            75       pion = particle2;
 76     }                                              76     }
 77     else{                                          77     else{
 78       nucleon = particle2;                         78       nucleon = particle2;
 79       pion = particle1;                            79       pion = particle1;
 80     }                                              80     }
 81                                                    81     
 82     const G4double sqrtS = KinematicsUtils::to     82     const G4double sqrtS = KinematicsUtils::totalEnergyInCM(nucleon, pion);
 83                                                    83     
 84     const G4int iso = ParticleTable::getIsospi     84     const G4int iso = ParticleTable::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(pion->getType());
 85     G4double rdm = Random::shoot();                85     G4double rdm = Random::shoot();
 86                                                    86     
 87     ParticleType KaonType;                         87     ParticleType KaonType;
 88                                                    88     
 89     if(iso == 3 || iso == -3){                     89     if(iso == 3 || iso == -3){
 90       KaonType = ParticleTable::getKaonType(is     90       KaonType = ParticleTable::getKaonType(iso/3);
 91       pion->setType(ParticleTable::getAntiKaon     91       pion->setType(ParticleTable::getAntiKaonType(iso/3));
 92     }                                              92     }
 93     else if(pion->getType() == PiZero){            93     else if(pion->getType() == PiZero){
 94       if(rdm*6. < 1.){                             94       if(rdm*6. < 1.){
 95         KaonType = ParticleTable::getKaonType(     95         KaonType = ParticleTable::getKaonType(-iso);
 96         pion->setType(ParticleTable::getAntiKa     96         pion->setType(ParticleTable::getAntiKaonType(iso));
 97       }                                            97       }
 98       else if(rdm*6. < 2.){                        98       else if(rdm*6. < 2.){
 99         KaonType = ParticleTable::getKaonType(     99         KaonType = ParticleTable::getKaonType(iso);
100         pion->setType(ParticleTable::getAntiKa    100         pion->setType(ParticleTable::getAntiKaonType(-iso));
101       }                                           101       }
102       else{                                       102       else{
103         KaonType = ParticleTable::getKaonType(    103         KaonType = ParticleTable::getKaonType(iso);
104         pion->setType(ParticleTable::getAntiKa    104         pion->setType(ParticleTable::getAntiKaonType(iso));
105         nucleon->setType(ParticleTable::getNuc    105         nucleon->setType(ParticleTable::getNucleonType(-iso));
106       }                                           106       }
107     }                                             107     }
108     else{                                         108     else{
109       if(rdm*5. < 1.){                            109       if(rdm*5. < 1.){
110         KaonType = ParticleTable::getKaonType(    110         KaonType = ParticleTable::getKaonType(iso);
111         pion->setType(ParticleTable::getAntiKa    111         pion->setType(ParticleTable::getAntiKaonType(iso));
112       }                                           112       }
113       else if(rdm*5. < 3.){                       113       else if(rdm*5. < 3.){
114         KaonType = ParticleTable::getKaonType(    114         KaonType = ParticleTable::getKaonType(-iso);
115         pion->setType(ParticleTable::getAntiKa    115         pion->setType(ParticleTable::getAntiKaonType(iso));
116         nucleon->setType(ParticleTable::getNuc    116         nucleon->setType(ParticleTable::getNucleonType(iso));
117       }                                           117       }
118       else{                                       118       else{
119         KaonType = ParticleTable::getKaonType(    119         KaonType = ParticleTable::getKaonType(iso);
120         pion->setType(ParticleTable::getAntiKa    120         pion->setType(ParticleTable::getAntiKaonType(-iso));
121         nucleon->setType(ParticleTable::getNuc    121         nucleon->setType(ParticleTable::getNucleonType(iso));
122       }                                           122       }
123     }                                             123     }
124                                                   124     
125 #ifdef INCLXX_IN_GEANT4_MODE                   << 
126     // Erase the parent resonance information  << 
127     nucleon->setParentResonancePDGCode(0);     << 
128     nucleon->setParentResonanceID(0);          << 
129     pion->setParentResonancePDGCode(0);        << 
130     pion->setParentResonanceID(0);             << 
131 #endif                                         << 
132                                                << 
133     ParticleList list;                            125     ParticleList list;
134     list.push_back(nucleon);                      126     list.push_back(nucleon);
135     list.push_back(pion);                         127     list.push_back(pion);
136     const ThreeVector &rcol = pion->getPositio    128     const ThreeVector &rcol = pion->getPosition();
137     const ThreeVector zero;                       129     const ThreeVector zero;
138     Particle *kaon = new Particle(KaonType,zer    130     Particle *kaon = new Particle(KaonType,zero,rcol);
139     list.push_back(kaon);                         131     list.push_back(kaon);
140                                                   132     
141     PhaseSpaceGenerator::generateBiased(sqrtS,    133     PhaseSpaceGenerator::generateBiased(sqrtS, list, 0, angularSlope);
142                                                   134     
143     INCL_DEBUG("NpiToNKKb " << (kaon->getMomen    135     INCL_DEBUG("NpiToNKKb " << (kaon->getMomentum().theta()) * 180. / G4INCL::Math::pi << '\n');
144                                                   136     
145     fs->addModifiedParticle(nucleon);             137     fs->addModifiedParticle(nucleon);
146     fs->addModifiedParticle(pion);                138     fs->addModifiedParticle(pion);
147     fs->addCreatedParticle(kaon);                 139     fs->addCreatedParticle(kaon);
148                                                   140     
149   }                                               141   }
150 }                                                 142 }
151                                                   143