Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/particles/hadrons/mesons/src/G4JPsi.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 /particles/hadrons/mesons/src/G4JPsi.cc (Version 11.3.0) and /particles/hadrons/mesons/src/G4JPsi.cc (Version 3.2)


  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 //
                                                   >>  23 //
                                                   >>  24 // $Id: G4JPsi.cc,v 1.5.2.1 2001/06/28 19:10:56 gunter Exp $
                                                   >>  25 // GEANT4 tag $Name:  $
                                                   >>  26 //
                                                   >>  27 // 
 26 // -------------------------------------------     28 // ----------------------------------------------------------------------
 27 //      GEANT 4 class implementation file          29 //      GEANT 4 class implementation file
 28 //                                                 30 //
 29 //      Created                 Hisaya Kurashi     31 //      Created                 Hisaya Kurashige, 16 June 1997
 30 // *******************************************     32 // **********************************************************************
 31 //  New impelemenataion as an utility class  M <<  33 //  Change both methods to get the pointer into non-inlined H.Kurashige 4 Aug. 1998
 32 // ------------------------------------------- <<  34 // ------------------------------------------------------------
 33                                                    35 
 34 #include "G4JPsi.hh"                           <<  36 #include "g4std/fstream"
                                                   >>  37 #include "g4std/iomanip"
 35                                                    38 
 36 #include "G4ParticleTable.hh"                  <<  39 #include "G4JPsi.hh"
 37 #include "G4String.hh"                         << 
 38 #include "G4SystemOfUnits.hh"                  << 
 39                                                    40 
 40 G4JPsi* G4JPsi::theInstance = nullptr;         <<  41 #include "G4DecayTable.hh"
 41                                                    42 
 42 G4JPsi* G4JPsi::Definition()                   <<  43 // ######################################################################
                                                   >>  44 // ###                          JPsi                            ###
                                                   >>  45 // ######################################################################
                                                   >>  46 
                                                   >>  47 G4JPsi::G4JPsi(
                                                   >>  48        const G4String&     aName,        G4double            mass,
                                                   >>  49        G4double            width,        G4double            charge,   
                                                   >>  50        G4int               iSpin,        G4int               iParity,    
                                                   >>  51        G4int               iConjugation, G4int               iIsospin,   
                                                   >>  52        G4int               iIsospin3,    G4int               gParity,
                                                   >>  53        const G4String&     pType,        G4int               lepton,      
                                                   >>  54        G4int               baryon,       G4int               encoding,
                                                   >>  55        G4bool              stable,       G4double            lifetime,
                                                   >>  56        G4DecayTable        *decaytable )
                                                   >>  57  : G4VMeson( aName,mass,width,charge,iSpin,iParity,
                                                   >>  58              iConjugation,iIsospin,iIsospin3,gParity,pType,
                                                   >>  59              lepton,baryon,encoding,stable,lifetime,decaytable )
 43 {                                                  60 {
 44   if (theInstance != nullptr) return theInstan <<  61   SetParticleSubType("J/psi");
 45   const G4String name = "J/psi";               <<  62   // Anti-particle of J/Psi is J/Psi itself  
 46   // search in particle table]                 <<  63   SetAntiPDGEncoding(encoding);
 47   G4ParticleTable* pTable = G4ParticleTable::G <<  64 
 48   G4ParticleDefinition* anInstance = pTable->F << 
 49   if (anInstance == nullptr) {                 << 
 50     // create particle                         << 
 51     //                                         << 
 52     //    Arguments for constructor are as fol << 
 53     //               name             mass     << 
 54     //             2*spin           parity  C- << 
 55     //          2*Isospin       2*Isospin3     << 
 56     //               type    lepton number  ba << 
 57     //             stable         lifetime     << 
 58     //             shortlived      subType     << 
 59                                                << 
 60     // clang-format off                        << 
 61    anInstance = new G4ParticleDefinition(      << 
 62                  name,    3.096916*GeV,    0.0 << 
 63                     2,              -1,        << 
 64                     0,               0,        << 
 65               "meson",               0,        << 
 66                 false,          0.0*ns,        << 
 67                 false,         "J/psi",        << 
 68     // clang-format on                         << 
 69   }                                            << 
 70   theInstance = static_cast<G4JPsi*>(anInstanc << 
 71   return theInstance;                          << 
 72 }                                              << 
 73                                                    65 
 74 G4JPsi* G4JPsi::JPsiDefinition()               << 
 75 {                                              << 
 76   return Definition();                         << 
 77 }                                                  66 }
 78                                                    67 
 79 G4JPsi* G4JPsi::JPsi()                         <<  68 // ......................................................................
                                                   >>  69 // ...                 static member definitions                      ...
                                                   >>  70 // ......................................................................
                                                   >>  71 //     
                                                   >>  72 //    Arguments for constructor are as follows
                                                   >>  73 //               name             mass          width         charge
                                                   >>  74 //             2*spin           parity  C-conjugation
                                                   >>  75 //          2*Isospin       2*Isospin3       G-parity
                                                   >>  76 //               type    lepton number  baryon number   PDG encoding
                                                   >>  77 //             stable         lifetime    decay table 
                                                   >>  78 
                                                   >>  79 // In this version, charged pions are set to stable
                                                   >>  80 G4JPsi G4JPsi::theJPsi(
                                                   >>  81         "J/psi",     3.09688*GeV,     0.087*MeV,          0., 
                                                   >>  82         2,              -1,            -1,          
                                                   >>  83         0,               0,            -1,             
                                                   >>  84         "meson",               0,             0,         443,
                                                   >>  85     false,          0.0*ns,          NULL
                                                   >>  86 );
                                                   >>  87 
                                                   >>  88 G4JPsi*  G4JPsi::JPsiDefinition(){return &theJPsi;}
                                                   >>  89 G4JPsi*  G4JPsi::JPsi(){return &theJPsi;}
                                                   >>  90 // initialization for static cut values
                                                   >>  91 G4double   G4JPsi::theJPsiLengthCut = -1.0;
                                                   >>  92 G4double*  G4JPsi::theJPsiKineticEnergyCuts = NULL;
                                                   >>  93 
                                                   >>  94 // **********************************************************************
                                                   >>  95 // **************************** SetCuts *********************************
                                                   >>  96 // **********************************************************************
                                                   >>  97 //  In this version Input Cut Value is meaning less
                                                   >>  98 //  theKineticEnergyCuts for all materials are set to LowestEnergy
                                                   >>  99 
                                                   >> 100 void G4JPsi::SetCuts(G4double aCut)
 80 {                                                 101 {
 81   return Definition();                         << 102   theCutInMaxInteractionLength = aCut;
                                                   >> 103 
                                                   >> 104   const G4MaterialTable* materialTable = G4Material::GetMaterialTable();
                                                   >> 105   // Create the vector of cuts in energy
                                                   >> 106   // corresponding to the stopping range cut
                                                   >> 107   if(theKineticEnergyCuts) delete [] theKineticEnergyCuts;
                                                   >> 108   theKineticEnergyCuts = new G4double [materialTable->length()];
                                                   >> 109 
                                                   >> 110   // Build range vector for every material, convert cut into energy-cut,
                                                   >> 111   // fill theKineticEnergyCuts and delete the range vector
                                                   >> 112   for (size_t J=0; J<materialTable->length(); J++)
                                                   >> 113   {
                                                   >> 114     theKineticEnergyCuts[J] = LowestEnergy;
                                                   >> 115   }
                                                   >> 116   theJPsiLengthCut = theCutInMaxInteractionLength;  
                                                   >> 117   theJPsiKineticEnergyCuts = theKineticEnergyCuts;
 82 }                                                 118 }
                                                   >> 119 
                                                   >> 120 
                                                   >> 121 
                                                   >> 122 
                                                   >> 123 
 83                                                   124