Changeset 6730

Timestamp:

03/03/05 22:18:05 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mjobs/MJPedestal.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -34 +34 @@
      - re-added MFEventSelector which had been removed by TB on the
        15.2.2005.
+     - added pulse position selector, as decided at the Berlin soft-
+       ware meeting.
    
    * mjobs/MJCalib.[h,cc]

trunk/MagicSoft/Mars/mjobs/MJPedestal.cc

@@ -62 +62 @@
 
 // Other basic classes
-#include "MExtractTimeAndChargeSlidingWindow.h"
+#include "MExtractTimeAndCharge.h"
 
 // parameter containers
@@ -78 +78 @@
 #include "MReadMarsFile.h"
 #include "MRawFileRead.h"
+#include "MRawEvtData.h"
 #include "MGeomApply.h"
 #include "MBadPixelsMerge.h"
@@ -83 +84 @@
 #include "MPedCalcPedRun.h"
 #include "MPedCalcFromLoGain.h"
+#include "MFTriggerPattern.h"
 
 #include <fstream>
@@ -978 +980 @@
 
     MFillH fillped(&hpedcam, "MPedestalCamIntermediate", "FillPedCam");
-    MFillH fillpul("MHCalibrationPulseTimeCam",   "MExtractedSignalCam", "FillPulseTime");
+    MFillH fillpul("MHCalibrationPulseTimeCam", "MRawEvtData", "FillPulseTime");
     //    MFillH fillped(&fPedestalHist, "MPedestalCamIntermediate", "FillPedCam");
     fillped.SetBit(MFillH::kDoNotDisplay);
@@ -996 +998 @@
       }
 
-    MExtractTimeAndChargeSlidingWindow exttest;
-
-    if (fPulsePosCheck && fExtractionType != kFundamental)
+    MFTriggerPattern fcalib("CalibFilter");
+    fcalib.RequireCalibration();
+    fcalib.SetInverted();
+
+    if (fPulsePosCheck)
       {
-        exttest.SetRange(0,25,0,0);
-        exttest.SetPedestals(&fPedestalCamIn);
-        tlist.AddToList(&exttest);
+        fillpul.SetFilter(&fcalib);
         tlist.AddToList(&fillpul);
       }
@@ -1186 +1188 @@
       }
     
-    if (fPulsePosCheck && fExtractionType != kFundamental)
+    if (fDisplayType!=kDisplayNone)
+        DisplayResult(plist);
+
+    if (!WriteResult())
+        return kFALSE;
+
+    if (fPulsePosCheck)
       {
+
+        MHCalibrationPulseTimeCam *hcam = (MHCalibrationPulseTimeCam*)plist.FindObject("MHCalibrationPulseTimeCam");
+        if (fPixelCheck)
+          {
+            hcam->DrawClone();
+            MHCalibrationPix &pix = (*hcam)[fCheckedPixId];
+            pix.DrawClone();
+          }
+        
         MCalibrationPulseTimeCam *cam = (MCalibrationPulseTimeCam*)plist.FindObject("MCalibrationPulseTimeCam");
+        MRawEvtData *data = (MRawEvtData*)plist.FindObject("MRawEvtData");
         if (cam)
           {
             const Float_t meanpulsetime = cam->GetAverageArea(0).GetHiGainMean();
-            const Float_t sigmapulsetime = cam->GetAverageArea(0).GetHiGainSigma();
-            *fLog << inf << "Mean pulse time (cosmics): " << meanpulsetime << "+-" << sigmapulsetime << endl;
+            const Float_t rmspulsetime = cam->GetAverageArea(0).GetHiGainRms();
+            *fLog << inf << "Mean pulse time (cosmics): " << meanpulsetime << "+-" << rmspulsetime << endl;
+            const Int_t newfirst = 1.5*rmspulsetime > 3. 
+              ? (Int_t)(meanpulsetime-1.5*rmspulsetime) 
+              : (Int_t)(meanpulsetime-3.);
+            const Int_t wshigain = fExtractor->InheritsFrom("MExtractTimeAndCharge") 
+              ? ((MExtractTimeAndCharge*)fExtractor)->GetWindowSizeHiGain()
+              : 6;
+            const Int_t wslogain = fExtractor->InheritsFrom("MExtractTimeAndCharge") 
+              ? ((MExtractTimeAndCharge*)fExtractor)->GetWindowSizeLoGain()
+              : 6;
+            const Int_t newlast  = 1.5*rmspulsetime > 3. 
+              ? (Int_t)(meanpulsetime+1.5*rmspulsetime)
+              : (Int_t)(meanpulsetime+3.);
+            *fLog << inf << "Try to set new range limits: ("<<newfirst<<","<<newlast+wshigain<<")"<<endl;
+            if (newfirst < 0)
+              {
+                *fLog << err << "Pulse is too much to the left, cannot go below 0! " << endl;
+                return kFALSE;
+              }
+            if (newlast+wshigain > data->GetNumHiGainSamples())
+              {
+                *fLog << err << "Pulse is too much to the right, cannot go beyond higain limits! " 
+                      << " Try to use a different extractor ... " << endl;
+                return kFALSE;
+              }
+            if (newlast+wslogain > data->GetNumLoGainSamples())
+              {
+                *fLog << err << "Pulse is too much to the right, cannot go beyond logain limits! " 
+                      << " Try to use a different extractor ... " << endl;
+                return kFALSE;
+              }
+            fExtractor->SetRange(newfirst,newlast+wshigain,newfirst,data->GetNumLoGainSamples());
+          }
+        else
+          {
+            *fLog << err << "Could not determine mean pulse position, abort... " << endl;
+            return kFALSE;
           }
       }
     
 
-    if (fDisplayType!=kDisplayNone)
-        DisplayResult(plist);
-
-    if (!WriteResult())
-        return kFALSE;
-
     *fLog << all << GetDescriptor() << ": Done." << endl;
     *fLog << endl << endl;