Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/decay/src/G4PionDecayMakeSpin.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/decay/src/G4PionDecayMakeSpin.cc (Version 11.3.0) and /processes/decay/src/G4PionDecayMakeSpin.cc (Version 9.1.p2)


  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 //                                                 27 //
 28                                                    28 
 29 #include "G4PionDecayMakeSpin.hh"                  29 #include "G4PionDecayMakeSpin.hh"
 30                                                    30 
 31 #include "G4Decay.hh"                              31 #include "G4Decay.hh"
 32 #include "G4DecayProducts.hh"                      32 #include "G4DecayProducts.hh"
 33                                                    33 
 34 #include "G4RandomDirection.hh"                    34 #include "G4RandomDirection.hh"
 35                                                    35 
 36 // constructor                                     36 // constructor
 37                                                    37 
 38 G4PionDecayMakeSpin::G4PionDecayMakeSpin(const     38 G4PionDecayMakeSpin::G4PionDecayMakeSpin(const G4String& processName)
 39                                : G4Decay(proce <<  39                                : G4Decay(processName) { }
 40 {                                              << 
 41   // set Process Sub Type                      << 
 42   SetProcessSubType(static_cast<int>(DECAY_Pio << 
 43                                                << 
 44 }                                              << 
 45                                                    40 
 46 G4PionDecayMakeSpin::~G4PionDecayMakeSpin() {      41 G4PionDecayMakeSpin::~G4PionDecayMakeSpin() { }
 47                                                    42 
 48 void G4PionDecayMakeSpin::DaughterPolarization     43 void G4PionDecayMakeSpin::DaughterPolarization(const G4Track& aTrack,
 49                                    G4DecayProd     44                                    G4DecayProducts* products)
 50 {                                                  45 {
 51   //  This routine deals only with particles t     46   //  This routine deals only with particles that can decay into a muon
 52   //                 pi+, pi-, K+, K- and K0_l     47   //                 pi+, pi-, K+, K- and K0_long
 53                                                    48 
 54   // get particle                                  49   // get particle
 55                                                    50 
 56   const G4DynamicParticle* aParticle = aTrack.     51   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
 57   const G4ParticleDefinition* aParticleDef = a     52   const G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();
 58                                                    53 
 59   G4ParticleDefinition* aMuonPlus =                54   G4ParticleDefinition* aMuonPlus =
 60             G4ParticleTable::GetParticleTable(     55             G4ParticleTable::GetParticleTable()->FindParticle("mu+");
 61   G4ParticleDefinition* aMuonMinus =               56   G4ParticleDefinition* aMuonMinus =
 62             G4ParticleTable::GetParticleTable(     57             G4ParticleTable::GetParticleTable()->FindParticle("mu-");
 63   G4ParticleDefinition* aPionPlus =                58   G4ParticleDefinition* aPionPlus = 
 64             G4ParticleTable::GetParticleTable(     59             G4ParticleTable::GetParticleTable()->FindParticle("pi+");
 65   G4ParticleDefinition* aPionMinus =               60   G4ParticleDefinition* aPionMinus =
 66             G4ParticleTable::GetParticleTable(     61             G4ParticleTable::GetParticleTable()->FindParticle("pi-");
 67   G4ParticleDefinition* aKaonPlus =                62   G4ParticleDefinition* aKaonPlus =
 68             G4ParticleTable::GetParticleTable(     63             G4ParticleTable::GetParticleTable()->FindParticle("kaon+");
 69   G4ParticleDefinition* aKaonMinus =               64   G4ParticleDefinition* aKaonMinus =
 70             G4ParticleTable::GetParticleTable(     65             G4ParticleTable::GetParticleTable()->FindParticle("kaon-");
 71   G4ParticleDefinition* aKaon0Long =               66   G4ParticleDefinition* aKaon0Long =
 72             G4ParticleTable::GetParticleTable(     67             G4ParticleTable::GetParticleTable()->FindParticle("kaon0L");
 73   G4ParticleDefinition* aNeutrinoMu =              68   G4ParticleDefinition* aNeutrinoMu =
 74             G4ParticleTable::GetParticleTable(     69             G4ParticleTable::GetParticleTable()->FindParticle("nu_mu");
 75   G4ParticleDefinition* aAntiNeutrinoMu =          70   G4ParticleDefinition* aAntiNeutrinoMu =
 76             G4ParticleTable::GetParticleTable(     71             G4ParticleTable::GetParticleTable()->FindParticle("anti_nu_mu");
 77                                                    72 
 78   if( aParticleDef == aPionPlus   ||               73   if( aParticleDef == aPionPlus   ||
 79       aParticleDef == aPionMinus  ||               74       aParticleDef == aPionMinus  || 
 80       aParticleDef == aKaonPlus   ||               75       aParticleDef == aKaonPlus   ||
 81       aParticleDef == aKaonMinus  ||               76       aParticleDef == aKaonMinus  ||
 82       aParticleDef == aKaon0Long ) {               77       aParticleDef == aKaon0Long ) {
 83   } else {                                         78   } else {
 84       return;                                      79       return;
 85   }                                                80   }
 86                                                    81 
 87   G4DynamicParticle* aMuon = nullptr;          <<  82   G4DynamicParticle* aMuon = NULL;
 88                                                    83 
 89   G4double emu(0), eneutrino(0);                   84   G4double emu(0), eneutrino(0);
 90   G4ThreeVector p_muon, p_neutrino;                85   G4ThreeVector p_muon, p_neutrino;
 91                                                    86 
 92   G4int numberOfSecondaries = products->entrie     87   G4int numberOfSecondaries = products->entries();
 93                                                    88 
 94   if (numberOfSecondaries > 0) {                   89   if (numberOfSecondaries > 0) {
 95     for (G4int index=0; index < numberOfSecond <<  90     for (G4int index=0; index < numberOfSecondaries; index++)
 96       G4DynamicParticle* aSecondary = (*produc <<  91     {
 97       const G4ParticleDefinition* aSecondaryDe <<  92         G4DynamicParticle* aSecondary = (*products)[index];
 98                                                <<  93         G4ParticleDefinition* aSecondaryDef = aSecondary->GetDefinition();
 99       if (aSecondaryDef == aMuonPlus ||        <<  94 
100     aSecondaryDef == aMuonMinus ) {            <<  95         if (aSecondaryDef == aMuonPlus ||
101   //          Muon+ or Muon-                   <<  96             aSecondaryDef == aMuonMinus ) {
102   aMuon = aSecondary;                          <<  97             //          Muon+ or Muon-
103   emu    = aSecondary->GetTotalEnergy();       <<  98             aMuon = aSecondary;
104   p_muon = aSecondary->GetMomentum();          <<  99             emu    = aSecondary->GetTotalEnergy();
105       } else if (aSecondaryDef == aNeutrinoMu  << 100             p_muon = aSecondary->GetMomentum();
106      aSecondaryDef == aAntiNeutrinoMu ) {      << 101         } else if (aSecondaryDef == aNeutrinoMu ||
107   //          Muon-Neutrino / Muon-Anti-Neutri << 102                    aSecondaryDef == aAntiNeutrinoMu ) {
108   eneutrino  = aSecondary->GetTotalEnergy();   << 103             //          Muon-Neutrino / Muon-Anti-Neutrino
109   p_neutrino = aSecondary->GetMomentum();      << 104             eneutrino  = aSecondary->GetTotalEnergy();
110       }                                        << 105             p_neutrino = aSecondary->GetMomentum();
                                                   >> 106         }
111     }                                             107     }
112   }                                               108   }
113                                                   109 
114   //  This routine deals only with decays with    110   //  This routine deals only with decays with a 
115   //  muon and mu-(anti)neutrinos in the final    111   //  muon and mu-(anti)neutrinos in the final state
116   if (aMuon == nullptr) return;                << 112 
                                                   >> 113   if (!aMuon) return;
117   if (eneutrino==0||emu==0) return;               114   if (eneutrino==0||emu==0) return;
118                                                   115 
119   G4ThreeVector spin(0,0,0);                      116   G4ThreeVector spin(0,0,0);
120                                                   117 
121   const G4DynamicParticle* theParentParticle =    118   const G4DynamicParticle* theParentParticle = products->GetParentParticle();
122                                                   119 
123   G4double amass = theParentParticle->GetMass(    120   G4double amass = theParentParticle->GetMass();
124   G4double emmu = aMuonPlus->GetPDGMass();        121   G4double emmu = aMuonPlus->GetPDGMass();
125                                                   122 
126   if (numberOfSecondaries == 2 ) {                123   if (numberOfSecondaries == 2 ) {
                                                   >> 124 
127      G4double scale = - (eneutrino - ( p_muon     125      G4double scale = - (eneutrino - ( p_muon * p_neutrino )/(emu+emmu));
128                                                   126 
129      p_muon = scale * p_muon;                     127      p_muon = scale * p_muon;
130      p_neutrino = emmu * p_neutrino;              128      p_neutrino = emmu * p_neutrino;
131      spin = p_muon + p_neutrino;                  129      spin = p_muon + p_neutrino;
132                                                   130 
133      scale = 2./(amass*amass-emmu*emmu);          131      scale = 2./(amass*amass-emmu*emmu);
134      spin = scale * spin;                         132      spin = scale * spin;
135                                                   133 
136      if (aParticle->GetCharge() < 0.0) spin =     134      if (aParticle->GetCharge() < 0.0) spin = -spin;
137                                                   135 
138   } else {                                        136   } else {
                                                   >> 137 
139        spin = G4RandomDirection();                138        spin = G4RandomDirection();
140                                                   139 
141   }                                               140   }
142                                                   141 
143   spin = spin.unit();                             142   spin = spin.unit();
144                                                   143 
145   aMuon->SetPolarization(spin.x(),spin.y(),spi    144   aMuon->SetPolarization(spin.x(),spin.y(),spin.z());
146                                                   145 
147   return;                                         146   return; 
148 }                                                 147 }
149                                                << 
150 void G4PionDecayMakeSpin::ProcessDescription(s << 
151 {                                              << 
152   outFile << GetProcessName()                  << 
153     << ": Decay of mesons that can decay into  << 
154     << " i.e. pi+, pi-, K+, K- and K0_long \n" << 
155     << " kinematics of daughters are dertermin << 
156     << " polarization of daughter particles ar << 
157 }                                              << 
158                                                << 
159                                                   148