Ignore:
Timestamp:
07/24/02 15:36:10 (22 years ago)
Author:
bigongia
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Simulation/Detector/ReflectorII/ph2cph.c

    r869 r1431  
    237237    + 2.0*SQR(rCT[0])*xCT[2] + 2.0*SQR(rCT[1])*xCT[2]);
    238238
    239   c = 2*rCT[0]*rCT[2]*x[0]*x[2] + 2*rCT[1]*rCT[2]*x[1]*x[2]
    240     - SQR(rCT[2])*SQR(x[0]) - SQR(rCT[2])*SQR(x[1])
    241     - SQR(rCT[0])*SQR(x[2]) - SQR(rCT[1])*SQR(x[2]);
     239  /*  FIXED Lines below, May 2002, AM : formerly (up to V0.4)
     240   *  there was a confusion between telescope coordinates xCT and
     241   *  the original coordinates x. Thanks to T. Hengstebeck for
     242   *  reporting the bug.
     243   */
     244
     245  c = 2*rCT[0]*rCT[2]*xCT[0]*xCT[2] + 2*rCT[1]*rCT[2]*xCT[1]*xCT[2]
     246    - SQR(rCT[2])*SQR(xCT[0]) - SQR(rCT[2])*SQR(xCT[1])
     247    - SQR(rCT[0])*SQR(xCT[2]) - SQR(rCT[1])*SQR(xCT[2]);
     248
    242249
    243250  if ( fabs(a) < 1.e-6 ) {
     
    583590  Debug("@19 xcam-AD %f %f \n", xcam[0], xcam[1]);
    584591  /* CBC */
    585 
     592 
    586593  /*
    587594  ++
Note: See TracChangeset for help on using the changeset viewer.