| Geant4 Cross Reference |
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 // G4PionRadiativeDecayChannel << 26 // ------------------------------------------------------------
>> 27 // GEANT 4 class header file
27 // 28 //
28 // Class decription: << 29 // History:
>> 30 // 30 July 2007 P.Gumplinger
>> 31 // Samples Radiative Pion Decay
>> 32 // Reference:
>> 33 // TRIUMF/PIENU Technote:
>> 34 // "Inclusion of pi->enug in the Monte Carlo"
>> 35 // by M. Blecher
29 // 36 //
30 // This class describes radiative pion decay k << 37 // ------------------------------------------------------------
31 // Rate is for gammas > 100keV. << 38 //
32 // NOTE: it gives incorrect energy spectrum fo << 39 //
33 << 40 //
34 // Author: P.Gumplinger, 30 July 2007 << 41 //
35 // Reference: M. Blecher, TRIUMF/PIENU Technot << 42
36 // "Inclusion of pi->enug in the Mo << 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
37 // ------------------------------------------- << 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
38 #ifndef G4PionRadiativeDecayChannel_hh << 45
39 #define G4PionRadiativeDecayChannel_hh 1 << 46 #ifndef G4PionRadiativeDecayChannel_h
>> 47 #define G4PionRadiativeDecayChannel_h 1
40 48
>> 49 #include "globals.hh"
>> 50 #include "Randomize.hh"
41 #include "G4ThreeVector.hh" 51 #include "G4ThreeVector.hh"
42 #include "G4VDecayChannel.hh" 52 #include "G4VDecayChannel.hh"
43 #include "Randomize.hh" <<
44 #include "globals.hh" <<
45 53
46 class G4PionRadiativeDecayChannel : public G4V 54 class G4PionRadiativeDecayChannel : public G4VDecayChannel
47 { 55 {
48 public: << 56 // Class Decription
49 G4PionRadiativeDecayChannel(const G4String << 57 // This class describes radiative pion decay kinemtics, but
50 ~G4PionRadiativeDecayChannel() override = << 58 // gives incorrect energy spectrum for neutrino
51 << 59
52 G4DecayProducts* DecayIt(G4double) overrid << 60 public: // With Description
53 << 61
54 protected: << 62 //Constructors
55 G4PionRadiativeDecayChannel() = default; << 63 G4PionRadiativeDecayChannel(const G4String& theParentName,
56 << 64 G4double theBR);
57 // Copy constructor and assignment operato << 65 // Destructor
58 G4PionRadiativeDecayChannel(const G4PionRa << 66 virtual ~G4PionRadiativeDecayChannel();
59 G4PionRadiativeDecayChannel& operator=(con << 67
>> 68 public: // With Description
>> 69
>> 70 virtual G4DecayProducts *DecayIt(G4double);
>> 71
>> 72 private:
>> 73
>> 74 G4double beta;
>> 75
>> 76 G4double cib;
>> 77 G4double csdp;
>> 78 G4double csdm;
>> 79 G4double cif;
>> 80 G4double cig;
>> 81
>> 82 G4double xl, yl, xu, yu, d2wmax;
>> 83
>> 84 G4double D2W(const G4double x, const G4double y);
>> 85
60 }; 86 };
>> 87
>> 88 inline G4double G4PionRadiativeDecayChannel::D2W(const G4double x,
>> 89 const G4double y)
>> 90 {
>> 91 G4double d2w = cib*(1.-y)*(1.+((1.-x)*(1.-x)))/((x*x)*(x+y-1.)) +
>> 92 csdp*(1.-x)*((x+y-1.)*(x+y-1.)) +
>> 93 csdm*(1.-x)*((1.-y)*(1.-y)) +
>> 94 cif*(x-1.)*(1.-y)/x +
>> 95 cig*(1.-y)*(1.-x+(x*x)/(x+y-1.))/x;
>> 96 return d2w;
>> 97 }
61 98
62 #endif 99 #endif
63 100