Ignore:
Timestamp:
04/29/05 12:39:44 (20 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/manalysis/MEventRateCalc.cc

    r6455 r6985  
    155155Int_t MEventRateCalc::Process()
    156156{
    157     const ULong_t exec = GetNumExecutions()-1;
     157    const ULong_t exec = GetNumExecutions()-fNumFirstEvent-1;
    158158
    159159    //*fLog << all << fNumFirstEvent << " " << exec << endl;
     
    172172    const UInt_t cnt = n1<n2 ? n : n1-n2;
    173173
     174    if (cnt<n/10)
     175        return kTRUE;
     176
     177
    174178    // Calculate the rate
    175179    const Double_t rate = (Double_t)cnt/(fTimes[n1%n]-fTimes[n2%n]);
     
    190194    fRate->SetRate(exec>1?rate:0, cnt);
    191195    fRate->SetReadyToSave();
     196
     197    Double_t diff = fTimes[n1%n]- fTimes[n2%n];
     198    Double_t time = fTimes[n2%n]+ (cnt-n/10.)/(n-n/10.)*diff/2;
     199
     200    fTimeRate->SetMean(time, time);
     201
     202    //if (cnt>1000)
     203    //*fLog << all << *fTimeRate << " " << *fTime << " " << diff << " " << cnt << " " << n << " " << (float)(cnt-n/2)/n << endl;
    192204
    193205    // Store the corresponding time
Note: See TracChangeset for help on using the changeset viewer.