Ignore:
Timestamp:
08/28/05 01:40:18 (19 years ago)
Author:
Daniela Dorner
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/datacenter/macros/fillganymed.C

    r7280 r7328  
    105105    }
    106106
    107     Double_t exc = fit->GetEventsExcess();
    108     Double_t sig = fit->GetEventsSignal();
    109     Double_t bgd = fit->GetEventsBackground();
    110     Double_t S   = fit->GetSignificance();
    111     Double_t f   = fit->GetScaleFactor();
     107    Int_t exc = (Int_t)fit->GetEventsExcess();
     108    Int_t sig = (Int_t)fit->GetEventsSignal();
     109    Int_t bgd = (Int_t)fit->GetEventsBackground();
     110    Float_t S = fit->GetSignificance();
     111    TString signif = Form("%5.1f", S);
     112    Float_t f = fit->GetScaleFactor();
     113    TString scale = Form("%5.2f", f);
    112114
    113115    MStatusArray arr;
     
    126128
    127129
    128     Double_t tm = vstime->Integral();
     130    Int_t tm = (Int_t)vstime->Integral();
    129131
    130132    TString dataset = fname(TRegexp("ganymed[0-9]+[.]root$"));
     
    145147    cout << "  Total eff. on-time: " << tm << "s = " << tm/3600. << "h" << endl;
    146148
    147     cout << "  Excess       Rate:  " << exc*60/tm << " evts/min"  << endl;
    148     cout << "  Background   Rate:  " << bgd*60/tm << " evts/min"  << endl;
    149     cout << "  Significance Rate:  " << S/TMath::Sqrt(tm/3600.) << " sigma/sqrt(h)" << endl;
     149//    cout << "  Excess       Rate:  " << exc*60/tm << " evts/min"  << endl;
     150//    cout << "  Background   Rate:  " << bgd*60/tm << " evts/min"  << endl;
     151//    cout << "  Significance Rate:  " << S/TMath::Sqrt(tm/3600.) << " sigma/sqrt(h)" << endl;
    150152
    151153
     
    155157        query = Form("INSERT Ganymed SET"
    156158                     " fDataSetNumber=%d,"
    157                      " fExcessEvents=%f, "
    158                      " fBackgroundEvents=%f, "
    159                      " fSignalEvents=%f, "
     159                     " fExcessEvents=%d, "
     160                     " fBackgroundEvents=%d, "
     161                     " fSignalEvents=%d, "
    160162                     " fSignificance=%f, "
    161163                     " fScaleFactor=%f, "
    162                      " fEffOnTime=%f, "
    163                      " fExcessRate=%f, "
    164                      " fBackgroundRate=%f, "
    165                      " fSignificanceRate=%f ",
    166                      ds, exc, bgd, sig, S, f, tm,
    167                      exc*60/tm, bgd*60/tm,
    168                      S/TMath::Sqrt(tm/3600));
     164                     " fEffOnTime=%d ", //, "
     165//                     " fExcessRate=%f, "
     166//                     " fBackgroundRate=%f, "
     167//                     " fSignificanceRate=%f ",
     168                     ds, exc, bgd, sig, S, f, tm);//,
     169//                     exc*60/tm, bgd*60/tm,
     170//                     S/TMath::Sqrt(tm/3600));
    169171    }
    170172    else
    171173    {
    172174        query = Form("UPDATE Ganymed SET"
    173                      " fExcessEvents=%f, "
    174                      " fBackgroundEvents=%f, "
    175                      " fSignalEvents=%f, "
     175                     " fExcessEvents=%d, "
     176                     " fBackgroundEvents=%d, "
     177                     " fSignalEvents=%d, "
    176178                     " fSignificance=%f, "
    177179                     " fScaleFactor=%f, "
    178                      " fEffOnTime=%f, "
    179                      " fExcessRate=%f, "
    180                      " fBackgroundRate=%f, "
    181                      " fSignificanceRate=%f ",
     180                     " fEffOnTime=%d " //, "
     181//                     " fExcessRate=%f, "
     182//                     " fBackgroundRate=%f, "
     183//                     " fSignificanceRate=%f ",
    182184                     " WHERE fDataSetNumber=%d ",
    183                      exc, bgd, sig, S, f, tm,
    184                      exc*60/tm, bgd*60/tm,
    185                      S/TMath::Sqrt(tm/3600), ds);
    186     }
    187 
     185                     exc, bgd, sig, S, f, tm, ds);
     186//                     exc*60/tm, bgd*60/tm,
     187//                     S/TMath::Sqrt(tm/3600), ds);
     188    }
     189
     190    cout << "q: " << query << endl;
    188191    if (dummy)
    189192        return 0;
Note: See TracChangeset for help on using the changeset viewer.