Pion 0.0.2

Simple particle physics library
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106

use crate::AntiQuark;
use crate::Quark;

///Currently only supports proton and neutron
#[derive(Debug,Copy,Clone)]
pub enum Baryon{
     Proton(f64), 
     Neutron(f64),
}
/// Currently only supports anti-proton and anti-neutron
#[derive(Debug,Copy,Clone)]
pub enum AntiBaryon{
     Proton(f64), 
     Neutron(f64),
}


 impl Baryon{
     ///Compose proton or neutron using quarks
  pub const fn new(x: Quark, y: Quark, z: Quark, ke: f64)->Baryon{
      
          match (x,y,z) {
          
          (Quark::Up, Quark::Up, Quark::Down)   => return Baryon::Proton(ke),
           _                                    => return Baryon::Neutron(ke),
 
          }
       }
       
  pub const fn charge(&self)->f64{
         match self{
          Baryon::Proton(_)  => 1f64,
          Baryon::Neutron(_) => 0f64,
          } 
  }      
       // MeV
  pub  const fn rest_mass(&self)->f64{
       
        match self{
          Baryon::Proton(_)  => 938.27231,
          Baryon::Neutron(_) => 939.5656,
          } 
       }
       
  pub  fn rest_mass_kg(&self)->f64{
      self.rest_mass()*1.79E-30
    }    
       
   pub const fn kinetic(&self)-> f64{
          match self {
          Baryon::Proton(x)  => *x,
          Baryon::Neutron(x) => *x,
          }
   }
           //relativistic momentum
   pub  fn momentum(&self)->f64{
           self.rest_mass()*self.kinetic() / (1.0 - (self.kinetic()*self.kinetic())/89875517873681764f64).sqrt()
        }
        
  }    
  
  
  impl AntiBaryon{
     
  pub const fn new(x: AntiQuark, y: AntiQuark, z: AntiQuark, ke: f64)->AntiBaryon{
      
          match (x,y,z) {
          
          (AntiQuark::Up, AntiQuark::Up, AntiQuark::Down)   => return AntiBaryon::Proton(ke),
           _                                                => return AntiBaryon::Neutron(ke),
 
          }
       }
       
  pub const fn charge(&self)->f64{
         match self{
          AntiBaryon::Proton(_)  => 1f64,
          AntiBaryon::Neutron(_) => 0f64,
          } 
  }      
       /// MeV
  pub const fn rest_mass(&self)->f64{
       
        match self{
          AntiBaryon::Proton(_)  => 938.27231,
          AntiBaryon::Neutron(_) => 939.5656,
          } 
       }
       
  pub  fn rest_mass_kg(&self)->f64{
      self.rest_mass()*1.79E-30
    }    
       
   pub const fn kinetic(&self)-> f64{
          match self {
          AntiBaryon::Proton(x)  => *x,
          AntiBaryon::Neutron(x) => *x,
          }
   }
           //relativistic momentum
   pub  fn momentum(&self)->f64{
           self.rest_mass()*self.kinetic() / (1.0 - (self.kinetic()*self.kinetic())/89875517873681764f64).sqrt()
        }
        
  }