Changeset 16377


Ignore:
Timestamp:
05/28/13 13:40:12 (12 years ago)
Author:
tbretz
Message:
Added the updated formula for the decrease of energy threshold with zenith angle directly using the result from Roger's thesis
Location:
trunk/FACT++/src
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/FACT++/src/makedata.cc

    r16282 r16377  
    188188
    189189        // Relative  energy threshold prediction
    190         const double cs = cos((90+hrz.alt)*M_PI/180);
    191         const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
     190        //const double cs = cos((90+hrz.alt)*M_PI/180);
     191        //const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
     192        const double ratio = pow(cos((90-hrz.alt)*M_PI/180, -2.664);
    192193
    193194        // Add points to curve

  • trunk/FACT++/src/makeplots.cc

    r16282 r16377  
    309309
    310310            // Relative  energy threshold prediction
    311             const double cs = cos((90+hrz.alt)*M_PI/180);
    312             const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
     311            //const double cs = cos((90+hrz.alt)*M_PI/180);
     312            //const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
     313            const double ratio = pow(cos((90-hrz.alt)*M_PI/180, -2.664);
    313314
    314315            // Add points to curve
Note: See TracChangeset for help on using the changeset viewer.