Changeset 14882 for trunk


Ignore:
Timestamp:
02/07/13 10:57:35 (12 years ago)
Author:
tbretz
Message:
Added a new preliminary formula for current prediction based on a higher statistical bases.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/FACT++/src/makeplots.cc

    r14723 r14882  
    309309
    310310            // Current prediction
    311             const double lc = angle*hrzm.alt*pow(disk, 6)/360/360;
    312             const double cur = lc>0 ? 7.7+4942*lc : 7.7;
     311            const double cang = sin(angle   *M_PI/180);
     312            const double calt = sin(hrzm.alt*M_PI/180);
     313
     314            const double lc = calt>0 ? cang*sqrt(calt)*pow(disk, 3) : -1;
     315            const double cur = lc>0 ? 8+104.5*lc : 8;
    313316
    314317            // Relative  energy threshold prediction
     
    334337
    335338            if (no_limits || (cur<max_current && 90-hrz.alt<max_zd))
    336                 g3.SetPoint(g3.GetN(), axis, ratio*cur/7.7);
     339                g3.SetPoint(g3.GetN(), axis, ratio*cur/8);
    337340
    338341            if (no_limits || (cur<max_current && 90-hrz.alt<max_zd))
Note: See TracChangeset for help on using the changeset viewer.