Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/eventgenerator/pythia/py8decayer/src/Py8DecayerPhysics.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 /examples/extended/eventgenerator/pythia/py8decayer/src/Py8DecayerPhysics.cc (Version 11.3.0) and /examples/extended/eventgenerator/pythia/py8decayer/src/Py8DecayerPhysics.cc (Version 11.0.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 //                                                 26 //
 27 /// \file eventgenerator/pythia/pythia8decayer     27 /// \file eventgenerator/pythia/pythia8decayer/src/Py8DecayerPhysics.cc
 28 /// \brief Implementation of the Py8DecayerPhy     28 /// \brief Implementation of the Py8DecayerPhysics class
 29 ///                                                29 ///
 30 /// \author J. Yarba; FNAL                         30 /// \author J. Yarba; FNAL
 31                                                    31 
 32 #include "Py8DecayerPhysics.hh"                    32 #include "Py8DecayerPhysics.hh"
 33                                                << 
 34 #include "Py8Decayer.hh"                           33 #include "Py8Decayer.hh"
 35                                                    34 
 36 #include "G4Decay.hh"                          << 
 37 #include "G4DecayTable.hh"                     << 
 38 #include "G4ParticleDefinition.hh"                 35 #include "G4ParticleDefinition.hh"
 39 #include "G4ProcessManager.hh"                     36 #include "G4ProcessManager.hh"
                                                   >>  37 #include "G4Decay.hh"
                                                   >>  38 #include "G4DecayTable.hh"
 40                                                    39 
 41 // factory                                         40 // factory
 42 //                                                 41 //
 43 #include "G4PhysicsConstructorFactory.hh"          42 #include "G4PhysicsConstructorFactory.hh"
 44 //                                                 43 //
 45 // register it with contructor factory             44 // register it with contructor factory
 46 //                                                 45 //
 47 G4_DECLARE_PHYSCONSTR_FACTORY(Py8DecayerPhysic     46 G4_DECLARE_PHYSCONSTR_FACTORY(Py8DecayerPhysics);
 48                                                    47 
 49 //....oooOO0OOooo........oooOO0OOooo........oo     48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 50                                                    49 
 51 Py8DecayerPhysics::Py8DecayerPhysics(G4int) :  <<  50 Py8DecayerPhysics::Py8DecayerPhysics(G4int)
                                                   >>  51   : G4VPhysicsConstructor("Py8DecayerPhysics")
                                                   >>  52 {
                                                   >>  53 }
 52                                                    54 
 53 //....oooOO0OOooo........oooOO0OOooo........oo     55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 54                                                    56 
 55 Py8DecayerPhysics::~Py8DecayerPhysics() {}     <<  57 Py8DecayerPhysics::~Py8DecayerPhysics() 
                                                   >>  58 {
                                                   >>  59 }
 56                                                    60 
 57 //....oooOO0OOooo........oooOO0OOooo........oo     61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 58                                                    62 
 59 void Py8DecayerPhysics::ConstructParticle()        63 void Py8DecayerPhysics::ConstructParticle()
 60 {                                                  64 {
 61   // Nothing needs to be done here             <<  65    // Nothing needs to be done here
 62 }                                                  66 }
 63                                                    67 
 64 //....oooOO0OOooo........oooOO0OOooo........oo     68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 65                                                    69 
 66 void Py8DecayerPhysics::ConstructProcess()         70 void Py8DecayerPhysics::ConstructProcess()
 67 {                                                  71 {
 68   // Adding external decayer to G4Decay proces <<  72 
 69   // G4Decay will use the external decayer if  <<  73    // Loop over all particles instantiated and add external decayer
 70   // assigned to an unstable particle and that <<  74    // to all decay processes where there's no decay table, plus tau's
 71   // have its decay table.                     <<  75 
 72                                                <<  76    // Create external decayer for G4Decay process
 73   // Loop over all particles instantiated and  <<  77    // NOTE: The extDecayer will be deleted in G4Decay destructor
 74   // decay table for tau's and B+/- so that th <<  78    Py8Decayer* extDecayer = new Py8Decayer();
 75   // the external decayer (Pythia8).           <<  79 
 76                                                <<  80    auto particleIterator=GetParticleIterator();
 77   // NOTE: The extDecayer will be deleted in G <<  81    particleIterator->reset();
 78                                                <<  82    while ((*particleIterator)())
 79   Py8Decayer* extDecayer = new Py8Decayer();   <<  83    {    
 80   G4bool setOnce = true;                       <<  84       G4ParticleDefinition* particle = particleIterator->value();
 81                                                <<  85       G4ProcessManager* pmanager = particle->GetProcessManager();    
 82   auto particleIterator = GetParticleIterator( <<  86 /*
 83   particleIterator->reset();                   <<  87       if ( verboseLevel > 1 ) {
 84   while ((*particleIterator)()) {              <<  88          G4cout << "Setting ext decayer for: " 
 85     G4ParticleDefinition* particle = particleI <<  89                 <<  particleIterator->value()->GetParticleName()
 86                                                <<  90                 << G4endl;
 87     // remove native/existing decay table for  <<  91       } 
 88     // a)tau's                                 <<  92 */    
 89     // b) B+/-                                 << 
 90     // so that G4Decay will use the external d << 
 91     if (std::abs(particle->GetPDGEncoding()) = << 
 92       if (particle->GetDecayTable()) {         << 
 93         delete particle->GetDecayTable();      << 
 94         particle->SetDecayTable(nullptr);      << 
 95         /*                                     << 
 96                   if ( verboseLevel > 1 ) {    << 
 97                      G4cout << "Use ext decaye << 
 98                         <<  particleIterator-> << 
 99                         << G4endl;             << 
100                   }                            << 
101         */                                     << 
102       }                                        << 
103     }                                          << 
104                                                << 
105     if (setOnce)                               << 
106     // One G4Decay object is shared by all uns << 
107     // Thus, we set the external decayer only  << 
108     {                                          << 
109       G4ProcessManager* pmanager = particle->G << 
110       G4ProcessVector* processVector = pmanage     93       G4ProcessVector* processVector = pmanager->GetProcessList();
111       for (size_t i = 0; i < processVector->le <<  94       for ( size_t i=0; i<processVector->length(); ++i ) 
112         G4Decay* decay = dynamic_cast<G4Decay* <<  95       {    
113         if (decay) {                           <<  96          G4Decay* decay = dynamic_cast<G4Decay*>((*processVector)[i]);
114           decay->SetExtDecayer(extDecayer);    <<  97          if ( decay ) 
115           setOnce = false;                     <<  98          {
116         }                                      <<  99             // remove native/existing decay table for
117       }                                        << 100             // a)tau's 
118     }                                          << 101             // b) B+/- 
119   }                                            << 102             // and replace with external decayer
                                                   >> 103             if ( std::abs(particle->GetPDGEncoding()) == 15 ||
                                                   >> 104                  std::abs(particle->GetPDGEncoding()) == 521 )
                                                   >> 105             {
                                                   >> 106                if ( particle->GetDecayTable() )
                                                   >> 107                {
                                                   >> 108                   delete particle->GetDecayTable();
                                                   >> 109                   particle->SetDecayTable(nullptr);
                                                   >> 110                }
                                                   >> 111                decay->SetExtDecayer(extDecayer);
                                                   >> 112             }
                                                   >> 113             // now set external decayer to all particles 
                                                   >> 114             // that don't yet have a decay table
                                                   >> 115             if ( !particle->GetDecayTable() )
                                                   >> 116             {
                                                   >> 117                decay->SetExtDecayer(extDecayer);
                                                   >> 118             }
                                                   >> 119          }
                                                   >> 120       }              
                                                   >> 121    }
120                                                   122 
121   return;                                      << 123    return;
                                                   >> 124    
122 }                                                 125 }
123                                                   126 
124 //....oooOO0OOooo........oooOO0OOooo........oo    127 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
125                                                   128