Changeset 7328 for trunk/MagicSoft/Mars/datacenter/macros
- Timestamp:
- 08/28/05 01:40:18 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/datacenter/macros/fillganymed.C
r7280 r7328 105 105 } 106 106 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); 112 114 113 115 MStatusArray arr; … … 126 128 127 129 128 Double_t tm =vstime->Integral();130 Int_t tm = (Int_t)vstime->Integral(); 129 131 130 132 TString dataset = fname(TRegexp("ganymed[0-9]+[.]root$")); … … 145 147 cout << " Total eff. on-time: " << tm << "s = " << tm/3600. << "h" << endl; 146 148 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; 150 152 151 153 … … 155 157 query = Form("INSERT Ganymed SET" 156 158 " fDataSetNumber=%d," 157 " fExcessEvents=% f, "158 " fBackgroundEvents=% f, "159 " fSignalEvents=% f, "159 " fExcessEvents=%d, " 160 " fBackgroundEvents=%d, " 161 " fSignalEvents=%d, " 160 162 " fSignificance=%f, " 161 163 " 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)); 169 171 } 170 172 else 171 173 { 172 174 query = Form("UPDATE Ganymed SET" 173 " fExcessEvents=% f, "174 " fBackgroundEvents=% f, "175 " fSignalEvents=% f, "175 " fExcessEvents=%d, " 176 " fBackgroundEvents=%d, " 177 " fSignalEvents=%d, " 176 178 " fSignificance=%f, " 177 179 " 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 ", 182 184 " 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; 188 191 if (dummy) 189 192 return 0;
Note:
See TracChangeset
for help on using the changeset viewer.