include/Geant4/G4MuonicAtomDecay.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 #ifndef G4MuonicAtomDecay_h
0027 #define G4MuonicAtomDecay_h 1
0028
0029 //
0030 //
0031 //
0032 //
0033 // ------------------------------------------------------------
0034 //  GEANT 4 class header file
0035 //
0036 //  History:
0037 //
0038 //  20170522 K L Genser first implementation somewhat based on
0039 //  G4RadioactiveDecay and code by V.Ivantchenko & K.Lynch
0040 //
0041 // ------------------------------------------------------------
0042
0043 #include "G4VRestDiscreteProcess.hh"
0044
0045 class G4MuonicAtom;
0046 class G4Ions;
0047 class G4HadronicInteraction;
0048 class G4HadFinalState;
0049
0050 class G4MuonicAtomDecay : public G4VRestDiscreteProcess
0051 {
0052
0053   // Class Description:
0054   //
0055   // MuonicAtom process in which Muon is captured by the nucleus or decays in orbit
0056
0057 public:  // with description
0058
0059   explicit G4MuonicAtomDecay(G4HadronicInteraction* hiptr=nullptr,
0060                              // G4HadronicInteraction* diptr=nullptr,
0061                              const G4String& processName="MuonicAtomDecay");
0062
0063   virtual ~G4MuonicAtomDecay();
0064
0065   virtual G4bool IsApplicable(const G4ParticleDefinition&);
0066
0067   //  virtual void PreparePhysicsTable(const G4ParticleDefinition&);
0068
0069   //  virtual void BuildPhysicsTable(const G4ParticleDefinition&);
0070
0071   inline void  SetVerboseLevel(G4int value) {verboseLevel = value;}
0072   // Sets the VerboseLevel
0073
0074   inline G4int GetVerboseLevel() const {return verboseLevel;}
0075   // Returns the VerboseLevel
0076
0077   virtual void ProcessDescription(std::ostream& outFile) const;
0078
0079   virtual G4double PostStepGetPhysicalInteractionLength(const G4Track& track,
0080                                                         G4double   previousStepSize,
0081                                                         G4ForceCondition* condition);
0082
0083   virtual G4double AtRestGetPhysicalInteractionLength(const G4Track&,
0084                                                       G4ForceCondition* );
0085
0086   inline virtual G4VParticleChange* AtRestDoIt(const G4Track& theTrack,
0087                                        const G4Step& theStep)
0088   {return DecayIt(theTrack, theStep);}
0089
0090   inline virtual G4VParticleChange* PostStepDoIt(const G4Track& theTrack,
0091                                          const G4Step& theStep)
0092   {return DecayIt(theTrack, theStep);}
0093
0094   virtual G4double GetMeanFreePath(const G4Track& aTrack,
0095                                    G4double   previousStepSize,
0096                                    G4ForceCondition* condition
0097                                    );
0098   //  Calculates from the macroscopic cross section a mean
0099   //  free path, the value is returned in units of distance.
0100
0101   virtual G4double GetMeanLifeTime(const G4Track& aTrack,G4ForceCondition* condition);
0102   //  Calculates the mean life-time (i.e. for decays) of the
0103   //  particle at rest due to the occurence of the given process,
0104   //  or converts the probability of interaction (i.e. for
0105   //  annihilation) into the life-time of the particle for the
0106   //  occurence of the given process.
0107
0108
0109 private:
0110   // hide copy constructor and assignment operator as private
0111   G4MuonicAtomDecay(const G4MuonicAtomDecay &right);
0112   G4MuonicAtomDecay& operator=(const G4MuonicAtomDecay &right);
0113
0114   G4VParticleChange* DecayIt(const G4Track& theTrack,
0115                              const G4Step&  theStep);
0116
0117   void FillResult(G4HadFinalState * aR, const G4Track & aT);
0118
0119   void DumpState(const G4Track&, const G4String&, G4ExceptionDescription&);
0120
0121   G4ParticleChange theTotalResult;
0122
0123   const G4double fMuMass;
0124
0125   G4HadronicInteraction* cmptr; // the capture model owned by its registry
0126   //  G4HadronicInteraction* dmptr; // the DIO model owned by its registry
0127
0128   G4int verboseLevel;
0129
0130
0131 };
0132 #endif