Ignore:
Timestamp:
02/23/05 17:10:12 (20 years ago)
Author:
mazin
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mtemp/mmpi/MSkyPlot.cc

    r6190 r6667  
    634634      Double_t dx, dy, beta, tanbeta, alphapar, distpar;
    635635      Double_t logsize, lgsize, lgsize2, dist2, hadr;
    636       const Double_t log3000 = log(3000.);
     636      const Double_t log3000n = log(3000.*0.18);  // now in phe, thus must shift the offset
    637637      const Float_t fDistOffset = 0.9;
    638638     
     
    725725                (y_cam)      (sin(gam)   cos(gam))   (ytilde)   (y_0)
    726726              */
    727               //xsourcam = - (cosgam * xsource - singam * ysource) + x_0;
    728               //ysourcam = - (singam * xsource + cosgam * ysource) + y_0;
     727              xsourcam = - (cosgam * xsource - singam * ysource) + x_0;
     728              ysourcam = - (singam * xsource + cosgam * ysource) + y_0;
    729729             
    730730             
    731731              /*    end derotatiom    */
    732                          xsourcam = xsource;
    733                          ysourcam = ysource;
     732              //           xsourcam = xsource;
     733              //           ysourcam = ysource;
    734734
    735735              /* calculate ALPHA und DIST according to the source position */
     
    758758                // apply cuts
    759759                logsize = log(fHillas->GetSize());
    760                 lgsize = logsize-log3000;
     760                lgsize = logsize-log3000n;
    761761                lgsize2 = lgsize*lgsize;
    762762                dist2 = distpar*distpar - fDistOffset*fDistOffset;
Note: See TracChangeset for help on using the changeset viewer.