Ignore:
Timestamp:
07/31/01 10:46:18 (23 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mhist/MHMcRate.cc

    r896 r900  
    5151    fBackSim=-1.0;          // Number of simulated showers for the background
    5252
    53     fTrigRate= -1.0;        // Trigger rate in Hz
    54     fTrigRateError= -1.0;   // Estimated error for the trigger rate in Hz
     53    fTriggerRate= -1.0;        // Trigger rate in Hz
     54    fTriggerRateError= -1.0;   // Estimated error for the trigger rate in Hz
    5555}
    5656
     
    6767    fFlux0=-1.0;            // dn/dE = fFlux0 * E^{-a}
    6868
    69     fShowRate= -1.0;        // Showers rate in Hz
    70     fShowRateError=0.0;     // Estimated error of shower rate in Hz
     69    fShowerRate= -1.0;        // Showers rate in Hz
     70    fShowerRateError=0.0;     // Estimated error of shower rate in Hz
    7171}
    7272
     
    8585    fFlux0=-1.0;            // dn/dE = fFlux0 * E^{-a}
    8686
    87     fShowRate= showrate;               // Showers rate in Hz
    88     fShowRateError=sqrt(showrate);     // Estimated error of shower rate in Hz
     87    fShowerRate= showrate;               // Showers rate in Hz
     88    fShowerRateError=sqrt(showrate);     // Estimated error of shower rate in Hz
    8989}
    9090
     
    103103    fFlux0=flux0;           // dn/dE = fFlux0 * E^{-a}
    104104
    105     fShowRate= -1.0;
    106     fShowRateError=0.0;
     105    fShowerRate= -1.0;
     106    fShowerRateError=0.0;
    107107}
    108108
     
    130130//  update the limits for energy, theta, phi and impact parameter
    131131//
    132 void MHMcRate::Boundaries(Float_t energy, Float_t theta,
    133                           Float_t phi, Float_t impact)
     132void MHMcRate::UpdateBoundaries(Float_t energy, Float_t theta,
     133                                Float_t phi, Float_t impact)
    134134{
    135135  // It updates the limit values
     
    163163    const Double_t epowmin = pow(fEnergyMin, specidx);
    164164
    165     fShowRate = fFlux0/specidx*(epowmax-epowmin);
     165    fShowerRate = fFlux0/specidx*(epowmax-epowmin);
    166166
    167167    if (fPartId!=1)
    168         fShowRate *= (fPhiMax-fPhiMin)*(cos(fThetaMax)-cos(fThetaMin));
     168        fShowerRate *= (fPhiMax-fPhiMin)*(cos(fThetaMax)-cos(fThetaMin));
    169169
    170170    const Double_t impactdiff = fImpactMax-fImpactMin;
    171171
    172     fShowRate *= TMath::Pi()*(impactdiff/100.0*impactdiff/100.0);
    173 
    174     fShowRateError = sqrt(fShowRate);
     172    fShowerRate *= TMath::Pi()*(impactdiff/100.0*impactdiff/100.0);
     173
     174    fShowerRateError = sqrt(fShowerRate);
    175175
    176176    const Double_t anal2 = 1.0-anal*160.0e-9;
     
    179179    // Then the trigger rate and its error is evaluated
    180180    if(fBackTrig<0){
    181         fTrigRateError = sqrt((trig*fShowRate*fShowRate/(simu*simu)) +
     181        fTriggerRateError = sqrt((trig*fShowerRate*fShowerRate/(simu*simu)) +
    182182                              (anal2*anal2*1/(fBackSim*back2*back2)));
    183183        fBackTrig=0;
    184184    }
    185185    else
    186         fTrigRateError = sqrt((trig*fShowRate*fShowRate/(simu*simu)) +
     186        fTriggerRateError = sqrt((trig*fShowerRate*fShowerRate/(simu*simu)) +
    187187                              (anal2*anal2*fBackTrig/(back2*back2)));
    188188
    189     fTrigRate = trig*fShowRate/simu + anal2*fBackTrig/back2;
     189    fTriggerRate = trig*fShowerRate/simu + anal2*fBackTrig/back2;
    190190}
    191191
     
    194194//  print the trigger rate
    195195//
    196 void MHMcRate::Print()
    197 {
    198     *fLog << "Incident rate " << fShowRate << " Hz " << endl;
    199     *fLog << "Trigger Rate " << fTrigRate << " +- " << fTrigRateError << " Hz" << endl;
     196void MHMcRate::Print(Option_t *)
     197{
     198    *fLog << "Incident rate " << fShowerRate << " Hz " << endl;
     199    *fLog << "Trigger Rate " << fTriggerRate << " +- " << fTriggerRateError << " Hz" << endl;
    200200}
    201201
     
    204204//  draw the trigger rate
    205205//
    206 void MHMcRate::Draw(Option_t* option)
     206void MHMcRate::Draw(Option_t *)
    207207{
    208208  *fLog << "To be iplemented" << endl;
Note: See TracChangeset for help on using the changeset viewer.