Changeset 1358
- Timestamp:
- 06/12/02 14:46:53 (23 years ago)
- Location:
- trunk/WuerzburgSoft/Thomas/mphys
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc
r1356 r1358 288 288 } 289 289 290 #include <iostream.h>291 292 290 MPhoton *MElectron::DoInvCompton(Double_t theta) 293 291 { -
trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc
r1357 r1358 120 120 const Double_t f = sqrt(sqrbetah*sqrbetae)*salpha; 121 121 122 cout << " {" << f << "," << E << "," << atan(tan1) << "," << atan(-tan2) << "} " << flush;122 // cout << " {" << f << "," << E << "," << atan(tan1) << "," << atan(-tan2) << "} " << flush; 123 123 124 124 MElectron &p0 = *new MElectron(E*(1.-f), gamma->GetZ()); -
trunk/WuerzburgSoft/Thomas/mphys/phys.C
r1357 r1358 212 212 if (phot.Integral(-log10(Eg)-8, -10.5)==0) 213 213 { 214 // ???????????????????215 hist.Fill(Eg, Eg*Eg*weight);216 angle.Fill(p->GetPsi()*kRad2Deg, p->GetTheta()*kRad2Deg, weight);217 hista.Fill(p->GetTheta()*kRad2Deg, weight);218 delete listg.Remove(p);219 214 cout << "z" << flush; 220 215 continue; 221 216 } 222 217 223 Double_t Ep = pow(10, phot.GetRandom()); 224 218 Double_t Ep = pow(10, phot.GetRandom()); 225 219 Double_t theta = RandomTheta(Eg, Ep); 226 220 if (theta==0)
Note:
See TracChangeset
for help on using the changeset viewer.