- Timestamp:
- 04/21/05 16:49:41 (20 years ago)
- Location:
- trunk/MagicSoft/Mars
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/Changelog
r6968 r6969 24 24 2005/04/21 Thomas Bretz 25 25 26 * mjobs/MJCalibrateSignal.cc: 26 * mhcalib/MHCalibrationRelTimeCam.cc: 27 - fixed usage of a uninitialized pointer 28 29 * mjobs/MJCalibrateSignal.cc, mjobs/MJCalibration.cc: 27 30 - fixed usage of MFTriggerPattern 28 31 -
trunk/MagicSoft/Mars/mhcalib/MHCalibrationRelTimeCam.cc
r6919 r6969 22 22 ! 23 23 \* ======================================================================== */ 24 24 25 ///////////////////////////////////////////////////////////////////////////// 25 // 26 // MHCalibrationRelTimeCam 27 // 28 // Fills the extracted relative arrival times of MArrivalTimeCam into 26 // 27 // MHCalibrationRelTimeCam 28 // 29 // Fills the extracted relative arrival times of MArrivalTimeCam into 29 30 // the MHCalibrationPix-classes MHCalibrationPix for every: 30 31 // 31 // - Pixel, stored in the TObjArray's MHCalibrationCam::fHiGainArray 32 // or MHCalibrationCam::fHiGainArray, respectively, depending if 32 // - Pixel, stored in the TObjArray's MHCalibrationCam::fHiGainArray 33 // or MHCalibrationCam::fHiGainArray, respectively, depending if 33 34 // MArrivalTimePix::IsLoGainUsed() is set. 34 35 // 35 // - Average pixel per AREA index (e.g. inner and outer for the MAGIC camera), 36 // stored in the TObjArray's MHCalibrationCam::fAverageHiGainAreas and 36 // - Average pixel per AREA index (e.g. inner and outer for the MAGIC camera), 37 // stored in the TObjArray's MHCalibrationCam::fAverageHiGainAreas and 37 38 // MHCalibrationCam::fAverageHiGainAreas 38 39 // 39 // - Average pixel per camera SECTOR (e.g. sectors 1-6 for the MAGIC camera), 40 // stored in the TObjArray's MHCalibrationCam::fAverageHiGainSectors 41 // and MHCalibrationCam::fAverageHiGainSectors 42 // 43 // Every relative time is calculated as the difference between the individual 44 // pixel arrival time and the one of pixel 1 (hardware number: 2). 45 // The relative times are filled into a histogram and an array, in order to perform 46 // a Fourier analysis (see MHGausEvents). The signals are moreover averaged on an 40 // - Average pixel per camera SECTOR (e.g. sectors 1-6 for the MAGIC camera), 41 // stored in the TObjArray's MHCalibrationCam::fAverageHiGainSectors 42 // and MHCalibrationCam::fAverageHiGainSectors 43 // 44 // Every relative time is calculated as the difference between the individual 45 // pixel arrival time and the one of pixel 1 (hardware number: 2). 46 // The relative times are filled into a histogram and an array, in order to perform 47 // a Fourier analysis (see MHGausEvents). The signals are moreover averaged on an 47 48 // event-by-event basis and written into the corresponding average pixels. 48 49 // 49 // The histograms are fitted to a Gaussian, mean and sigma with its errors 50 // and the fit probability are extracted. If none of these values are NaN's and 51 // if the probability is bigger than MHGausEvents::fProbLimit (default: 0.5%), 50 // The histograms are fitted to a Gaussian, mean and sigma with its errors 51 // and the fit probability are extracted. If none of these values are NaN's and 52 // if the probability is bigger than MHGausEvents::fProbLimit (default: 0.5%), 52 53 // the fit is declared valid. 53 // Otherwise, the fit is repeated within ranges of the previous mean 54 // +- MHCalibrationPix::fPickupLimit (default: 5) sigma (see MHCalibrationPix::RepeatFit())55 // In case this does not make the fit valid, the histogram means and RMS's are56 // taken directly (see MHCalibrationPix::BypassFit()) and the following flags are set:54 // Otherwise, the fit is repeated within ranges of the previous mean 55 // - MHCalibrationPix::fPickupLimit (default: 5) sigma (see MHCalibrationPix::RepeatFit()) 56 // In case this does not make the fit valid, the histogram means and RMS's are 57 // taken directly (see MHCalibrationPix::BypassFit()) and the following flags are set: 57 58 // - MBadPixelsPix::SetUncalibrated( MBadPixelsPix::kRelTimeNotFitted ) and 58 // - MBadPixelsPix::SetUnsuitable( MBadPixelsPix::kUnreliableRun ) 59 // 60 // Outliers of more than MHCalibrationPix::fPickupLimit (default: 5) sigmas 61 // from the mean are counted as Pickup events (stored in MHCalibrationPix::fPickup) 62 // 63 // The class also fills arrays with the signal vs. event number, creates a fourier 64 // spectrum (see MHGausEvents::CreateFourierSpectrum()) and investigates if the 65 // projected fourier components follow an exponential distribution. 66 // In case that the probability of the exponential fit is less than 59 // - MBadPixelsPix::SetUnsuitable( MBadPixelsPix::kUnreliableRun ) 60 // 61 // Outliers of more than MHCalibrationPix::fPickupLimit (default: 5) sigmas 62 // from the mean are counted as Pickup events (stored in MHCalibrationPix::fPickup) 63 // 64 // The class also fills arrays with the signal vs. event number, creates a fourier 65 // spectrum (see MHGausEvents::CreateFourierSpectrum()) and investigates if the 66 // projected fourier components follow an exponential distribution. 67 // In case that the probability of the exponential fit is less than 67 68 // MHGausEvents::fProbLimit (default: 0.5%), the following flags are set: 68 69 // - MBadPixelsPix::SetUncalibrated( MBadPixelsPix::kRelTimeOscillating ) and … … 83 84 // the number involved pixels in order to be able to compare it to the average of 84 85 // sigmas in the camera. 85 // 86 // 86 87 ///////////////////////////////////////////////////////////////////////////// 87 88 #include "MHCalibrationRelTimeCam.h" … … 136 137 const TString MHCalibrationRelTimeCam::gsHistYTitle = "Nr. events"; 137 138 const TString MHCalibrationRelTimeCam::fgReferenceFile = "mjobs/calibrationref.rc"; 139 138 140 // -------------------------------------------------------------------------- 139 141 // … … 808 810 pix.CreateGraphEvents(); 809 811 TGraph *gr = pix.GetGraphEvents(); 810 811 xaxe2->SetTitle(gr->GetXaxis()->GetTitle()); 812 yaxe2->SetTitle(gr->GetYaxis()->GetTitle()); 812 if (gr) 813 { 814 xaxe2->SetTitle(gr->GetXaxis()->GetTitle()); 815 yaxe2->SetTitle(gr->GetYaxis()->GetTitle()); 816 } 813 817 814 818 null2->Draw(); -
trunk/MagicSoft/Mars/mjobs/MJCalibration.cc
r6963 r6969 1849 1849 MTriggerPatternDecode trgpat; 1850 1850 MFTriggerPattern fcalib("CalibFilter"); 1851 fcalib.SetDefault(kFALSE); 1851 1852 fcalib.DenyCalibration(MFTriggerPattern::kPrescaled); 1852 1853
Note:
See TracChangeset
for help on using the changeset viewer.