Changeset 1365 for trunk/WuerzburgSoft/Thomas/mphys
- Timestamp:
- 06/14/02 09:08:54 (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc
r1364 r1365 12 12 ClassImp(MParticle); 13 13 14 /*********************** 14 /************************************************** 15 15 * 16 * calc r from z:16 * H0 = 50./3.0857e19; // [km / Mpc s] --> [s^-1] 17 17 * 18 * R = 2*c/H0 *(1+z - sqrt(1+z)) 19 * 20 * z = -0.5 * (3 + R' +/- sqrt(R'+1)) R' = R*H0/c 21 * 22 * H0 = 50./3.0857e19; // [km / Mpc s] --> [s^-1] 23 * 24 ************************ 25 */ 26 #include <iostream.h> 27 #include <iomanip.h> 18 **************************************************/ 28 19 29 20 Double_t MParticle::ZofR(Double_t *x, Double_t *k) … … 50 41 const Double_t R = 1./(1-r*H0/c/2); // [1] 51 42 52 53 43 return R*R - 1; 54 44 } … … 67 57 const Double_t R = c/H0 * 2 * (z1 - sqrt(z1)); // [m] 68 58 69 return R * pc/ly/1e3; // [kpc]59 return R * pc/ly/1e3; // [kpc] 70 60 */ 71 61 Double_t z1 = x[0] + 1; … … 79 69 const Double_t R = c/H0 * 2 * (1 - 1./sqrt(z1)); // [m] 80 70 81 return R * pc/ly/1e3;// [kpc]71 return R * pc/ly/1e3; // [kpc] 82 72 } 83 73 … … 107 97 108 98 return E*E / (exp(EkT)-1.); // [GeV^2] 109 110 99 } 111 100 … … 115 104 // 116 105 // default constructor 117 // creates an a list of histograms for all pixels and both gain channels118 106 // 119 120 //121 // set the name and title of this object122 //123 124 /*125 fName = name ? name : "MParticle";126 fTitle = title ? title : "Container for a particle";127 */128 107 } 129 108 … … 159 138 fPsi = atan2(r(1), r(0)); 160 139 161 if (fTheta*2 > TMath::Pi()) 140 /* 141 if (fTheta*2 > TMath::Pi()) 162 142 fTheta = fabs(fTheta-TMath::Pi()); 143 */ 163 144 } 164 145
Note:
See TracChangeset
for help on using the changeset viewer.