- Timestamp:
- 02/23/05 17:10:12 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mtemp/mmpi/MSkyPlot.cc
r6190 r6667 634 634 Double_t dx, dy, beta, tanbeta, alphapar, distpar; 635 635 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 637 637 const Float_t fDistOffset = 0.9; 638 638 … … 725 725 (y_cam) (sin(gam) cos(gam)) (ytilde) (y_0) 726 726 */ 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; 729 729 730 730 731 731 /* end derotatiom */ 732 xsourcam = xsource;733 ysourcam = ysource;732 // xsourcam = xsource; 733 // ysourcam = ysource; 734 734 735 735 /* calculate ALPHA und DIST according to the source position */ … … 758 758 // apply cuts 759 759 logsize = log(fHillas->GetSize()); 760 lgsize = logsize-log3000 ;760 lgsize = logsize-log3000n; 761 761 lgsize2 = lgsize*lgsize; 762 762 dist2 = distpar*distpar - fDistOffset*fDistOffset;
Note:
See TracChangeset
for help on using the changeset viewer.