| 1 | void plotlpsimple()
|
|---|
| 2 | {
|
|---|
| 3 | MParList plist;
|
|---|
| 4 |
|
|---|
| 5 | MDrsCalibration calib;
|
|---|
| 6 | if (!calib.ReadFits("/loc_data/raw/2011/10/30/20111030_023.drs.fits"))
|
|---|
| 7 | return;
|
|---|
| 8 | plist.AddToList(&calib);
|
|---|
| 9 |
|
|---|
| 10 | MTaskList tlist;
|
|---|
| 11 | plist.AddToList(&tlist);
|
|---|
| 12 |
|
|---|
| 13 | MRawFitsRead read1;
|
|---|
| 14 | read1.AddFile("/loc_data/raw/2011/10/30/20111030_028.fits");
|
|---|
| 15 | read1.LoadMap("../FACT++/FACTmap111030.txt");
|
|---|
| 16 | tlist.AddToList(&read1);
|
|---|
| 17 |
|
|---|
| 18 | MGeomApply apply;
|
|---|
| 19 | tlist.AddToList(&apply);
|
|---|
| 20 |
|
|---|
| 21 | MDrsCalibApply drscal;
|
|---|
| 22 | tlist.AddToList(&drscal);
|
|---|
| 23 |
|
|---|
| 24 | MExtractTimeAndChargeSpline extractor;
|
|---|
| 25 | extractor.SetRange(180, 240); // 0, 1023
|
|---|
| 26 | extractor.SetRiseTimeHiGain(13);
|
|---|
| 27 | extractor.SetFallTimeHiGain(23);
|
|---|
| 28 | extractor.SetChargeType(MExtralgoSpline::kIntegralRel); // 1.4V-1.8V T=290
|
|---|
| 29 | extractor.SetSaturationLimit(600000);
|
|---|
| 30 | extractor.SetNoiseCalculation(kFALSE);
|
|---|
| 31 | tlist.AddToList(&extractor);
|
|---|
| 32 |
|
|---|
| 33 | MHCamEvent evt2(0, "Extra'd", "Extracted Calibration Signal;;S [cnts/sl]");
|
|---|
| 34 | MHCamEvent evt9(4, "ArrTm", "Extracted ArrivalTime;;T");
|
|---|
| 35 |
|
|---|
| 36 | MHSectorVsTime histrms;
|
|---|
| 37 | MHSectorVsTime histtm;
|
|---|
| 38 | histrms.SetTitle("Pulser Amplitude (maximum extracted with spline);;Amplitude [mV]");
|
|---|
| 39 | histrms.SetType(0);
|
|---|
| 40 | histtm.SetTitle("Pulser arrival bin (leading edge extracted with spline);;Arrival time [slice]");
|
|---|
| 41 | histtm.SetType(0);
|
|---|
| 42 |
|
|---|
| 43 | MFillH fill2(&evt2, "MExtractedSignalCam", "FillExtractedSignal");
|
|---|
| 44 | MFillH fill9(&evt9, "MArrivalTimeCam", "FillArrivalTime");
|
|---|
| 45 | MFillH fillrms(&histrms, "MExtractedSignalCam", "FillAmplVsTime");
|
|---|
| 46 | MFillH filltm(&histtm, "MArrivalTimeCam", "FillATVsTime");
|
|---|
| 47 |
|
|---|
| 48 | tlist.AddToList(&fill2);
|
|---|
| 49 | tlist.AddToList(&fill9);
|
|---|
| 50 | tlist.AddToList(&fillrms);
|
|---|
| 51 | tlist.AddToList(&filltm);
|
|---|
| 52 |
|
|---|
| 53 | MStatusDisplay *d = new MStatusDisplay;
|
|---|
| 54 |
|
|---|
| 55 | MEvtLoop loop;
|
|---|
| 56 | loop.SetParList(&plist);
|
|---|
| 57 | loop.SetDisplay(d);
|
|---|
| 58 | loop.Eventloop();
|
|---|
| 59 | }
|
|---|
