G4ExcitedDeltaConstructor.hh

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 //
00027 // $Id$
00028 //
00029 // 
00030 // --------------------------------------------------------------
00031 //      GEANT 4 class implementation file 
00032 //
00033 //      History: first implementation, based on object model of
00034 //      10 oct 1998  H.Kurashige
00035 // ---------------------------------------------------------------
00036 #ifndef G4ExcitedDeltaConstructor_h
00037 #define G4ExcitedDeltaConstructor_h 1
00038 
00039 #include "globals.hh"
00040 #include "G4ios.hh"
00041 #include "G4ExcitedBaryonConstructor.hh"
00042 
00043 class G4ExcitedDeltaConstructor: public G4ExcitedBaryonConstructor
00044 {
00045   //This class is a utility class for construction 
00046   //short lived particles
00047 
00048   public:
00049     G4ExcitedDeltaConstructor();
00050     virtual  ~G4ExcitedDeltaConstructor();
00051 
00052   protected:
00053     virtual  G4int    GetEncoding(G4int iIsoSpin3, G4int idxState);
00054 
00055   protected:  
00056     virtual  G4bool   Exist( G4int ){return true;}
00057 
00058     virtual  G4int    GetQuarkContents(G4int, G4int);
00059     virtual  G4String GetName(G4int iIso3, G4int iState);
00060     virtual  G4String GetMultipletName(G4int iState);
00061 
00062     virtual  G4double GetMass( G4int state, G4int iso);
00063     virtual  G4double GetWidth( G4int state, G4int iso);
00064     virtual  G4int    GetiSpin(G4int iState);
00065     virtual  G4int    GetiParity(G4int iState);
00066     virtual  G4int    GetEncodingOffset(G4int iState);
00067 
00068     virtual  G4DecayTable* CreateDecayTable(const G4String& name,
00069                                             G4int iIso3, G4int iState,
00070                                             G4bool fAnti = false);
00071   private:
00072     G4DecayTable* AddNGammaMode( G4DecayTable* table, const G4String& name,
00073                                     G4double br, G4int iIso3, G4bool fAnti);
00074     G4DecayTable* AddNPiMode( G4DecayTable* table, const G4String& name,
00075                                      G4double br, G4int iIso3, G4bool fAnti);
00076     G4DecayTable* AddNRhoMode( G4DecayTable* table, const G4String& name,
00077                                      G4double br, G4int iIso3, G4bool fAnti);
00078     G4DecayTable* AddDeltaPiMode( G4DecayTable* table, const G4String& name,
00079                                      G4double br, G4int iIso3, G4bool fAnti);
00080     G4DecayTable* AddNStarPiMode( G4DecayTable* table, const G4String& name,
00081                                      G4double br, G4int iIso3, G4bool fAnti);
00082 
00083   public:
00084     enum     { NStates = 9  };
00085   private:
00086     enum     { DeltaIsoSpin = 3 };
00087 
00088   private:
00089     static const char* name[ NStates ];
00090     static const G4double mass[ NStates ];
00091     static const G4double width[ NStates ];
00092     static const G4int    iSpin[ NStates ];
00093     static const G4int    iParity[ NStates ];
00094     static const G4int    encodingOffset[ NStates ];
00095 
00096   public:
00097     enum     { NumberOfDecayModes = 5};
00098   private:
00099     enum     { NGamma=0, NPi=1, NRho=2, DeltaPi=3, NStarPi=4 };
00100   private:
00101     static const G4double bRatio[ NStates ][ NumberOfDecayModes];
00102 };
00103 
00104 inline
00105  G4double G4ExcitedDeltaConstructor::GetMass(G4int iState, G4int)
00106 { 
00107   return mass[iState]; 
00108 }
00109 
00110 inline
00111  G4double G4ExcitedDeltaConstructor::GetWidth(G4int iState, G4int)
00112 {
00113   return width[iState];
00114 }
00115 
00116 inline
00117  G4int    G4ExcitedDeltaConstructor::GetiSpin(G4int iState)
00118 {
00119   return iSpin[iState];
00120 }
00121 
00122 inline
00123  G4int    G4ExcitedDeltaConstructor::GetiParity(G4int iState)
00124 {
00125   return iParity[iState];
00126 }
00127 
00128 inline
00129  G4int    G4ExcitedDeltaConstructor::GetEncodingOffset(G4int iState)
00130 {
00131   return encodingOffset[iState];
00132 }
00133 
00134 inline
00135  G4int  G4ExcitedDeltaConstructor::GetQuarkContents(G4int iQ, G4int iIso3)
00136 {
00137   // Quark contents
00138   //    iIso3 = +3 : uuu
00139   //    iIso3 = +1 : uud
00140   //    iIso3 = -1 : udd
00141   //    iIso3 = -3 : ddd
00142   G4int quark=0;
00143   if ( iQ == 0 ){
00144     if ( iIso3 == -3 ){
00145     // d-quark
00146       quark = 1;
00147     } else {
00148     // u-quark
00149       quark = 2;
00150     }
00151   } else if ( iQ == 2 ){
00152     if ( iIso3 == +3 ){
00153     // u-quark
00154       quark = 2;
00155     } else {
00156     // d-quark
00157       quark = 1;
00158     }
00159   } else {
00160     if (( iIso3 == -1 )||( iIso3 == -3 )) {
00161     // d-quark
00162       quark = 1;
00163     } else {
00164     // u-quark
00165       quark = 2;
00166     }
00167   }
00168   return quark;
00169 }
00170 
00171 inline 
00172  G4String G4ExcitedDeltaConstructor::GetMultipletName(G4int iState)
00173 {
00174   return name[iState];
00175 }
00176 
00177 inline 
00178  G4String G4ExcitedDeltaConstructor::GetName(G4int iIso3, G4int iState)
00179 {
00180   G4String particle = name[iState];
00181   if ( iIso3 == -3 ){
00182     particle += "-";
00183   } else if ( iIso3 == -1 ){
00184     particle += "0";
00185   } else if ( iIso3 == +1 ){
00186     particle += "+";
00187   } else {
00188     particle += "++";
00189   }
00190   return particle;
00191 }
00192 #endif
00193 
00194 
00195 
00196 
00197 
00198 
00199 
00200 
00201 
00202 

Generated on Mon May 27 17:48:13 2013 for Geant4 by  doxygen 1.4.7