- Timestamp:
- 02/07/13 10:57:35 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/FACT++/src/makeplots.cc
r14723 r14882 309 309 310 310 // 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; 313 316 314 317 // Relative energy threshold prediction … … 334 337 335 338 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); 337 340 338 341 if (no_limits || (cur<max_current && 90-hrz.alt<max_zd))
Note:
See TracChangeset
for help on using the changeset viewer.