Ignore:
Timestamp:
08/20/02 12:38:36 (22 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc

    r1507 r1510  
    5757    const Double_t o = x[0];
    5858
     59    const Double_t s1 = TMath::Pi()*2;
     60    const Double_t s2 = a*h*c/E0;             // [m]
     61
     62    if (o<1e-4)
     63        return s2*s2/s1 * 4/3;
     64
    5965    const Double_t o1  = o+1;
    6066    const Double_t o21 = o*2+1;
    6167
    62     const Double_t s1 = TMath::Pi()*2;
    63     const Double_t s2 = a*h*c/E0;             // [m]
    6468    const Double_t s3 = o1/(o*o*o);
    6569    const Double_t s4 = o*2*o1/o21;
     
    365369    return bcost;
    366370
    367     Double_t theta = acos(bcost/beta);
    368 
    369     return theta;
     371    /*
     372     Double_t theta = acos(bcost/beta);
     373     return theta;
     374     */
    370375}
    371376
Note: See TracChangeset for help on using the changeset viewer.