Changeset 18287
- Timestamp:
- 08/22/15 00:34:25 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Mars/fact/analysis/callisto.C
r18274 r18287 21 21 22 22 MStatusDisplay *d = new MStatusDisplay; 23 24 // ------------------------------------------------------ 25 26 // ------------------------------------------------------ 27 28 // Calib: 51 / 90 / 197 (20% TH) 29 // Data: 52 / 64 / 104 (20% TH) 30 31 // Extraction range in slices. It will always(!) contain the full range 32 // of integration 33 const int first_slice = 25; // 10ns 34 const int last_slice = 225; // 125ns 35 36 const double heighttm = 0.5; // IntegralAbs { 1.5pe * 9.6mV/pe } / IntegralRel { 0.5 } 37 38 Long_t max = 0; // All 39 Long_t max0 = max; // Time marker 40 Long_t max1 = max; // Light pulser 41 //Long_t max2 = 3000; // Calibration ratio 42 Long_t max3 = max; // Pedestal Rndm 43 Long_t max4 = max; // Pedestal Ext 44 Long_t max5 = max; // Data 45 46 // ========================= Result ================================== 47 48 //double scale = 0.1; 49 double scale = 0.1024; 50 51 // ====================================================== 52 53 if (map && gSystem->AccessPathName(map, kFileExists)) 54 { 55 gLog << err << "ERROR - Cannot access mapping file '" << map << "'" << endl; 56 return 1; 57 } 58 59 // The sequence file which defines the files for the analysis 60 MSequence seq(seqfile); 61 if (!seq.IsValid()) 62 { 63 gLog << err << "ERROR - Sequence '" << seqfile << "' invalid!" << endl; 64 return 2; 65 } 23 66 24 67 MBadPixelsCam badpixels; … … 38 81 // 1393 39 82 40 // ------------------------------------------------------ 41 42 // ------------------------------------------------------ 43 44 // Calib: 51 / 90 / 197 (20% TH) 45 // Data: 52 / 64 / 104 (20% TH) 46 47 // Extraction range in slices. It will always(!) contain the full range 48 // of integration 49 const int first_slice = 25; // 10ns 50 const int last_slice = 225; // 125ns 51 52 const double heighttm = 0.5; // IntegralAbs { 1.5pe * 9.6mV/pe } / IntegralRel { 0.5 } 53 54 Long_t max = 0; // All 55 Long_t max0 = max; // Time marker 56 Long_t max1 = max; // Light pulser 57 //Long_t max2 = 3000; // Calibration ratio 58 Long_t max3 = max; // Pedestal Rndm 59 Long_t max4 = max; // Pedestal Ext 60 Long_t max5 = max; // Data 61 62 // ========================= Result ================================== 63 64 //double scale = 0.1; 65 double scale = 0.1024; 66 67 // ====================================================== 68 69 if (map && gSystem->AccessPathName(map, kFileExists)) 70 { 71 gLog << err << "ERROR - Cannot access mapping file '" << map << "'" << endl; 72 return 1; 73 } 74 75 // The sequence file which defines the files for the analysis 76 MSequence seq(seqfile); 77 if (!seq.IsValid()) 78 { 79 gLog << err << "ERROR - Sequence '" << seqfile << "' invalid!" << endl; 80 return 2; 83 // broken drs-board: 84 // https://www.fact-project.org/logbook/showthread.php?tid=3521 85 // in the analysis effect visible 10.1. - 25.5.2015 86 // according to logbook: 87 // first report of problem: 88 // https://www.fact-project.org/logbook/showthread.php?tid=2772&pid=15756#pid15756 89 // repair action: 90 // https://www.fact-project.org/logbook/showthread.php?tid=3298 91 // affected pixels: 92 // SW: 1193 1194 1195 1391 1392 1393 1304 1305 1306 93 if (seq.GetNight().GetNightAsInt()>20141114 && seq.GetNight().GetNightAsInt()<20150526) 94 { 95 gLog << "Exclude broken drs-board..." << endl; 96 badpixels[1193].SetUnsuitable(MBadPixelsPix::kUnsuitable); 97 badpixels[1194].SetUnsuitable(MBadPixelsPix::kUnsuitable); 98 badpixels[1195].SetUnsuitable(MBadPixelsPix::kUnsuitable); 99 badpixels[1391].SetUnsuitable(MBadPixelsPix::kUnsuitable); 100 badpixels[1392].SetUnsuitable(MBadPixelsPix::kUnsuitable); 101 badpixels[1393].SetUnsuitable(MBadPixelsPix::kUnsuitable); 102 badpixels[1304].SetUnsuitable(MBadPixelsPix::kUnsuitable); 103 badpixels[1305].SetUnsuitable(MBadPixelsPix::kUnsuitable); 104 badpixels[1306].SetUnsuitable(MBadPixelsPix::kUnsuitable); 81 105 } 82 106
Note:
See TracChangeset
for help on using the changeset viewer.