Changeset 8897 for trunk/MagicSoft/Mars/mjobs

Timestamp:

05/24/08 15:45:57 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mjobs/MJCalibrateSignal.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mjobs/MJCalibrateSignal.cc

@@ -494 +494 @@
 
     //------------------------------
+    //
+    // this is the filter to find pedestal events. For the Siegen FADC
+    // these are all events which are not calibration events because
+    // the pedestal is extracted from the lo-gain signal. For MUX
+    // data this are artifiially triggered events with the pedestal
+    // trigger flag, and for safty without Lvl1 or Lvl2 flag)
+    //
+    // For the time "before" the trigger pattern all events (Siegen FADC)
+    // can be considered to be pedestal, because this was also the time
+    // without artifially calibration events
+    //
     MFTriggerPattern fped;
     fped.SetDefault(kTRUE);
@@ -505 +516 @@
 
     //------------------------------
-
-    MFCosmics fcosmics;
-    fcosmics.SetNamePedestalCam("MPedestalFundamental");
-    MContinue contcos(&fcosmics, "ContTrigEvts");
-    contcos.SetInverted();
-
-    //MMcPedestalCopy pcopy;
+    //
+    // Apply a filter against cosmics (this is to make sure that the
+    // trigger system was working properly and no empty events survive)
+    // For every event 5% of the pixel must not be empty. In PostProcess
+    // an error is raised if more than 50% of the events were skipped.
+    //
+    MFCosmics fcosmicscal;
+    fcosmicscal.SetNamePedestalCam("MPedestalFundamental"); //CORRECT?
+    fcosmicscal.SetMaxEmptyPixels(0.05);
+    fcosmicscal.SetMaxAcceptedFraction(0.5);
+
+    MContinue contcoscal(&fcosmicscal, "ContCosmicsCal");
+    contcoscal.SetInverted();
+
+    //------------------------------
+    /*
+     MFCosmics fcosmicsped; // Def=20% empty pixels
+     fcosmicsped.SetNamePedestalCam("MPedestalFundamental"); //CORRECT?
+     MContinue contcosped(&fcosmicsped, "ContCosmicsPed");
+     contcosped.SetInverted();
+    */
+
+    //------------------------------
+    //
+    // Thie signal extractors
+    //
     MTaskEnv taskenv1("ExtractSignal");
     MTaskEnv taskenv2("ExtractTime");
@@ -529 +559 @@
     calib.AddPedestal("FromExtractorRndm");
     calib.SetPedestalFlag(MCalibrateData::kEvent);
+
     //----------------------------------------------------------
 
@@ -740 +771 @@
         tlist2.AddToList(&fcalib);     // MFTriggerPattern
         tlist2.AddToList(&tlist4);
+
+        tlist4.AddToList(&taskenv3);
+        tlist4.AddToList(&contcoscal); // MContinue/ContCosmicsCal
         if (IsUsePINDiode())
             tlist4.AddToList(&pinext); // MExtractPINDiode
         if (IsUseBlindPixel())
             tlist4.AddToList(&bldext); // MExtractBlindPixel
-        tlist4.AddToList(&taskenv3);
         tlist4.AddToList(&pacalc);     // MCalibCalcFromPast
         /*
@@ -777 +810 @@
       tlist2.AddToList(&taskenv1);
 
-    // remove all events which definitly don#t have a signal
-    // using MFCosmics (ContTrigEvts)
-    tlist2.AddToList(&contcos);      
+    // remove all events which definitly don't have a signal
+    // using MFCosmics (ContCosmicsPed)
+    // tlist2.AddToList(&contcosped);
 
     // Extract arrival time (if a dedicated extrator given)