- Timestamp:
- 08/09/06 11:01:30 (18 years ago)
- Location:
- trunk/MagicSoft/Mars
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/Changelog
r7869 r7870 22 22 * mbase/MMath.cc: 23 23 - changed MedianDev to support 68% division instead of 0.5% 24 25 * datacenter/macros/fillsignal.C: 26 - adapted to use Median for Pulse Offset and Gain 24 27 25 28 -
trunk/MagicSoft/Mars/NEWS
r7867 r7870 10 10 11 11 - general: Added a function to MMath which returns the Median of 12 the distribution abs(y[i]-median) , which is somehow the median13 counterpart of the rms.12 the distribution abs(y[i]-median) but at 68:32 division instead 13 of 50:50, which is somehow the mediani counterpart of the rms. 14 14 15 15 - showplot: -
trunk/MagicSoft/Mars/datacenter/macros/fillsignal.C
r7858 r7870 213 213 Double_t meanoff = TMath::Nint(pullo->GetMean()*100.)/100.; 214 214 Double_t rmsoff = TMath::Nint(pullo->GetRMS() *100.)/100.; 215 */ 216 Double_t meanoff = TMath::Nint(hilooff->GetMean()*100.)/100.; 217 Double_t rmsoff = TMath::Nint(hilooff->GetRMS() *100.)/100.; 218 219 Double_t meancal = TMath::Nint(hilocal->GetMean()*100.)/100.; 220 Double_t rmscal = TMath::Nint(hilocal->GetRMS() *100.)/100.; 215 */ 216 217 // USE MEDIAN INSTEAD? GetQuantiles(Int_t nprobSum, Double_t *q, const Double_t *probSum=0); 218 219 Double_t medoff = TMath::Nint(hilooff->GetMedian() *10000)/10000.; 220 Double_t devoff = TMath::Nint(hilooff->GetMedianDev()*10000)/10000.; 221 222 Double_t medcal = TMath::Nint(hilocal->GetMedian() *100)/100.; 223 Double_t devcal = TMath::Nint(hilocal->GetMedianDev()*100)/100.; 221 224 222 225 //get sequence number from the filename … … 245 248 TString rmspullo =Form("%6.2f", rmslo); 246 249 */ 247 TString me anpuloff =Form("%6.2f", meanoff);248 TString rmspuloff =Form("%6.2f", rmsoff);249 TString me anhilocal =Form("%6.2f", meancal);250 TString rmshilocal =Form("%6.2f", rmscal);250 TString medpuloff =Form("%7.4f", medoff); 251 TString devpuloff =Form("%7.4f", devoff); 252 TString medhilocal =Form("%6.2f", medcal); 253 TString devhilocal =Form("%6.2f", devcal); 251 254 252 255 if (meanextpul<0 && rmsextpul<0) … … 270 273 // cout << " Mean ext.HiGain PulsePos " << meanpulhi << " +- " << rmspulhi << endl; 271 274 // cout << " Mean ext.LoGain PulsePos " << meanpullo << " +- " << rmspullo << endl; 272 cout << " Lo-Hi gain offset: " << meanpuloff << " +- " << rmspuloff << endl;273 cout << " Hi/Lo gain ratio: " << me anhilocal << " +- " << rmshilocal << endl;275 cout << " Lo-Hi gain offset: " << medpuloff << " +- " << devpuloff << endl; 276 cout << " Hi/Lo gain ratio: " << medhilocal << " +- " << devhilocal << endl; 274 277 cout << endl; 275 278 … … 285 288 //" fPulsePosHiMean=%s, fPulsePosHiRms=%s, " 286 289 //" fPulsePosLoMean=%s, fPulsePosLoRms=%s, " 287 " fPulsePosOffMe an=%s, fPulsePosOffRms=%s, "288 " fHiLoGainRatioMe an=%s, fHiLoGainRatioRms=%s"290 " fPulsePosOffMed=%s, fPulsePosOffDev=%s, " 291 " fHiLoGainRatioMed=%s, fHiLoGainRatioDev=%s " 289 292 " WHERE fSequenceFirst='%d' ", 290 293 meanrmsinner.Data(), meanrmsouter.Data(), … … 294 297 //meanpulhi.Data(), rmspulhi.Data(), 295 298 //meanpullo.Data(), rmspullo.Data(), 296 me anpuloff.Data(), rmspuloff.Data(),297 me anhilocal.Data(), rmshilocal.Data(),299 medpuloff.Data(), devpuloff.Data(), 300 medhilocal.Data(), devhilocal.Data(), 298 301 seq); 299 302 -
trunk/MagicSoft/Mars/mbase/MMath.cc
r7868 r7870 181 181 182 182 // FIXME: GausProb() is a workaround. It should be taken into account in Median! 183 return TMath::Median(n, arr);184 /*183 //return TMath::Median(n, arr); 184 185 185 // Sort distribution 186 186 Long64_t idx[n]; … … 199 199 200 200 return dev; 201 */202 201 } 203 202
Note:
See TracChangeset
for help on using the changeset viewer.