Changeset 7069 for trunk/MagicSoft/Mars/mjobs

Timestamp:

05/20/05 11:59:31 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mjobs

Files:

• MJCalibration.cc (modified) (2 diffs)
• MJPedestal.cc (modified) (6 diffs)
• MJPedestal.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mjobs/MJCalibration.cc

@@ -1485 +1485 @@
 
   MCalibrationHiLoCam hilocam;
-  
-  *fLog << all << "Initializing High gain vs. Low gain intercalibration from " << fHiLoCalibFile << endl;
-  *fLog << all << endl;
-
   if (hilocam.Read()<=0)
   {
@@ -1500 +1496 @@
       return kFALSE;
   }
+
+  *fLog << all << "Hi-/Lo-Gain intercalibration constants read from " << fHiLoCalibFile << endl << endl;
 
   if (fCalibrationCam.GetSize() < 1)

trunk/MagicSoft/Mars/mjobs/MJPedestal.cc

@@ -822 +822 @@
 }
 
+Bool_t MJPedestal::PulsePosCheck(const MParList &plist) const
+{
+    if (fIsPixelCheck)
+    {
+        MHPedestalCam *hcam = (MHPedestalCam*)plist.FindObject("MHPedestalCam");
+        if (hcam)
+        {
+            MHPedestalPix &pix1 = (MHPedestalPix&)(*hcam)[fCheckedPixId];
+            pix1.DrawClone("");
+        }
+    }
+
+    if (!fIsPulsePosCheck)
+        return kTRUE;
+
+    Int_t numhigainsamples = 0;
+    Int_t numlogainsamples = 0;
+
+    Float_t meanpulsetime  = 0.;
+    Float_t rmspulsetime   = 0.;
+
+    if (IsUseMC())
+    {
+        //
+        // FIXME:
+        // The MC cannot run over the first 2000 pedestal events since almost all
+        // events are empty, therefore a pulse pos. check is not possible, either.
+        // For the moment, have to fix the problem hardcoded...
+        //
+        //            MMcEvt *evt = (MMcEvt*)plist.FindObject("MMcEvt");
+        //            const Float_t meanpulsetime = evt->GetFadcTimeJitter();
+        meanpulsetime = 4.5;
+        rmspulsetime  = 1.0;
+
+        numhigainsamples = 15;
+        numlogainsamples = 15;
+
+    }
+    else
+    {
+        if (fIsPixelCheck)
+        {
+            MHCalibrationPulseTimeCam *hcam = (MHCalibrationPulseTimeCam*)plist.FindObject("MHCalibrationPulseTimeCam");
+            if (!hcam)
+            {
+                *fLog << err << "MHCalibrationPulseTimeCam not found... abort." << endl;
+                return kFALSE;
+            }
+            hcam->DrawClone();
+            gPad->SaveAs(Form("%s/PulsePosTest_all.root",fPathOut.Data()));
+
+            MHCalibrationPix &pix = (*hcam)[fCheckedPixId];
+            pix.DrawClone();
+            gPad->SaveAs(Form("%s/PulsePosTest_Pixel%04d.root",fPathOut.Data(),fCheckedPixId));
+        }
+
+        MCalibrationPulseTimeCam *cam = (MCalibrationPulseTimeCam*)plist.FindObject("MCalibrationPulseTimeCam");
+        if (!cam)
+        {
+            *fLog << err << "MCalibrationPulseTimeCam not found... abort." << endl;
+            return kFALSE;
+        }
+
+        meanpulsetime = cam->GetAverageArea(0).GetHiGainMean();
+        rmspulsetime  = cam->GetAverageArea(0).GetHiGainRms();
+
+        MRawEvtData *data = (MRawEvtData*)plist.FindObject("MRawEvtData");
+        if (!data)
+        {
+            *fLog << err << "MRawEvtData not found... abort." << endl;
+            return kFALSE;
+        }
+
+        numhigainsamples = data->GetNumHiGainSamples();
+        numlogainsamples = data->GetNumLoGainSamples();
+    }
+
+    *fLog << all << "Mean pulse time (" << (IsUseMC()?"MC":"cosmics") << "): ";
+    *fLog << meanpulsetime << "+-" << rmspulsetime << endl;
+
+    const MExtractTimeAndCharge *ext = dynamic_cast<MExtractTimeAndCharge*>(fExtractor);
+
+    //
+    // Get the ranges for the new extractor setting
+    //
+    const Int_t newfirst = TMath::Nint(meanpulsetime-fExtractWinLeft);
+
+    Int_t wshigain = ext ? ext->GetWindowSizeHiGain() : 6;
+    if (wshigain > 6)
+        wshigain = 6;
+
+    Int_t wslogain = ext ? ext->GetWindowSizeLoGain() : 4;
+    if (wslogain > 4)
+        wslogain = 4;
+
+    const Int_t newlast  = TMath::Nint(meanpulsetime+fExtractWinRight);
+
+    *fLog << underline;
+    *fLog << "Try to set new range limits: (" << newfirst << "," << newlast;
+    *fLog << "+" << wshigain << "," << newfirst-1 << "," << newlast << "+";
+    *fLog << wslogain << ")" << endl;
+
+    //
+    // Check the ranges for the new extractor setting
+    //
+    if (newfirst < 0)
+    {
+        *fLog << err << "Pulse is too much to the left, cannot go below 0!" << endl;
+        return kFALSE;
+
+    }
+    if (newlast+wshigain > numhigainsamples+numlogainsamples-1)
+    {
+        *fLog << err << "Pulse is too much to the right, cannot go beyond limits: ";
+        *fLog << numhigainsamples << "+" << numlogainsamples << "-1" << endl;
+        *fLog << " Cannot extract at all!" << endl;
+        return kFALSE;
+    }
+    if (newlast+wslogain > numlogainsamples)
+    {
+        *fLog << err << "Pulse is too much to the right, cannot go beyond logain limits!" << endl;
+        *fLog << endl;
+        *fLog << "Try to use a different extractor (e.g. with a window size of only 4 sl.) or:" << endl;
+        *fLog << "Set the limit to a lower value (callisto.rc):" << endl;
+        *fLog << "    MJPedestalY2:ExtractWinRight: 4.5" << endl;
+        *fLog << "(ATTENTION, you will lose late cosmics pulses!)" << endl;
+        *fLog << endl;
+        return kFALSE;
+    }
+
+    //
+    // Set and store the new ranges
+    //
+    const Int_t hi0 = newfirst;
+    const Int_t hi1 = newlast+wshigain;
+
+    const Int_t lo0 = newfirst>0 ? newfirst-1 : newfirst;
+    const Int_t lo1 = numlogainsamples-1;
+
+    fExtractor->SetRange(hi0, hi1, lo0, lo1);
+    return kTRUE;
+}
+
 Bool_t MJPedestal::Process()
 {
@@ -989 +1132 @@
     }
     
-    if (fIsUseHists)
+    if (fIsUseHists && fExtractor)
       {
         if (fExtractor->InheritsFrom("MExtractTimeAndCharge"))
@@ -1021 +1164 @@
     if (fExtractor)
     {
+      fExtractor->SetPedestals(&fPedestalCamIn);
 
       if (fExtractionType!=kFundamental)
@@ -1065 +1209 @@
         *fLog << "windows." << endl;
     }
-    
-    /*
-    switch (fExtractType)
-    {
-    case kUseData:
-        *fLog << all << "TYPE: USEDATA " << fExtractor << endl;
-        taskenv.SetDefault(&pedlogain);
-        tlist.AddToList(&taskenv);
-
-        if (!SetupExtractor(plist, pedlogain))
-        {
-            *fLog << all <<  "SETTING TO: " << fExtractor << " " << fExtractor->GetNumHiGainSamples() << endl;
-            fExtractor->Print();
-            pedlogain.SetExtractWindow(15, (Int_t)TMath::Nint(fExtractor->GetNumHiGainSamples()));
-        }
-        break;
-
-    case kUsePedRun:
-        *fLog << all <<  "TYPE: USEPEDRUN " << fExtractor << endl;
-        taskenv.SetDefault(&pedcalc);
-        tlist.AddToList(&taskenv);
-
-        if (!SetupExtractor(plist, pedcalc))
-            pedcalc.SetExtractWindow(fExtractor->GetHiGainFirst(), TMath::Nint(fExtractor->GetNumHiGainSamples()));
-        break;
-
-    case kUseHists:
-        if (!fExtractor)
-        {
-            *fLog << err << GetDescriptor() << " - ERROR: ";
-            *fLog << "Extraction Type is kUseHists, but no extractor was set" << endl;
-            return kFALSE;
-        }
-
-        tlist.AddToList(fExtractor);
-        tlist.AddToList(&fillped);
-        break;
-    } */
-
-    /*
-    if (!fPathIn.IsNull())
-    {
-        delete fExtractor;
-        fExtractor = 0;
-    }
-    */
-
 
     //
@@ -1135 +1232 @@
     tlist.PrintStatistics();
 
-    if (fIsPixelCheck)
-      {
-        MHPedestalCam *hcam = (MHPedestalCam*)plist.FindObject("MHPedestalCam");
-        if (hcam)
-          {
-            MHPedestalPix &pix1 = (MHPedestalPix&)(*hcam)[fCheckedPixId];
-            pix1.DrawClone("");
-          }
-      }
-    
     if (fDisplayType!=kDisplayNone)
         DisplayResult(plist);
@@ -1151 +1238 @@
         return kFALSE;
 
-    if (fIsPulsePosCheck)
-      {
-        
-        Int_t numhigainsamples = 0;
-        Int_t numlogainsamples = 0;
-        Float_t meanpulsetime  = 0.;
-        Float_t rmspulsetime   = 0.;
-
-        if (IsUseMC())
-          {
-            //
-            // FIXME:
-            // The MC cannot run over the first 2000 pedestal events since almost all
-            // events are empty, therefore a pulse pos. check is not possible, either. 
-            // For the moment, have to fix the problem hardcoded...
-            //
-            //            MMcEvt *evt = (MMcEvt*)plist.FindObject("MMcEvt");
-            //            const Float_t meanpulsetime = evt->GetFadcTimeJitter();
-            meanpulsetime = 4.5;
-            rmspulsetime  = 1.0;
-
-            *fLog << all << "Mean pulse time (MC): " << meanpulsetime << "+-" << rmspulsetime << endl;
-
-            numhigainsamples = 15;
-            numlogainsamples = 15;
-
-          }
-        else
-          {
-            MHCalibrationPulseTimeCam *hcam = (MHCalibrationPulseTimeCam*)plist.FindObject("MHCalibrationPulseTimeCam");
-            if (fIsPixelCheck)
-              {
-                hcam->DrawClone();
-                gPad->SaveAs(Form("%s/PulsePosTest_all.root",fPathOut.Data()));
-                MHCalibrationPix &pix = (*hcam)[fCheckedPixId];
-                pix.DrawClone();
-                gPad->SaveAs(Form("%s/PulsePosTest_Pixel%04d.root",fPathOut.Data(),fCheckedPixId));
-              }
-            
-            MCalibrationPulseTimeCam *cam = (MCalibrationPulseTimeCam*)plist.FindObject("MCalibrationPulseTimeCam");
-            if (!cam)
-              {
-                *fLog << err << "Could not determine mean pulse position, abort... " << endl;
-                return kFALSE;
-              }
-            
-            meanpulsetime = cam->GetAverageArea(0).GetHiGainMean();
-            rmspulsetime  = cam->GetAverageArea(0).GetHiGainRms();
-
-            *fLog << all << "Mean pulse time (cosmics): " << meanpulsetime << "+-" << rmspulsetime << endl;
-
-            MRawEvtData *data = (MRawEvtData*)plist.FindObject("MRawEvtData");
-
-            numhigainsamples = data->GetNumHiGainSamples();
-            numlogainsamples = data->GetNumLoGainSamples();
-          }
-        
-        //
-        // Get the ranges for the new extractor setting
-        //
-        const Int_t newfirst = TMath::Nint(meanpulsetime-fExtractWinLeft);
-        Int_t wshigain = fExtractor->InheritsFrom("MExtractTimeAndCharge") 
-          ?  ((MExtractTimeAndCharge*)fExtractor)->GetWindowSizeHiGain()
-          : 6;
-
-        if (wshigain > 6)
-          wshigain = 6;
-
-        Int_t wslogain = fExtractor->InheritsFrom("MExtractTimeAndCharge") 
-          ? ((MExtractTimeAndCharge*)fExtractor)->GetWindowSizeLoGain()
-          : 6;
-        if (wslogain > 4)
-          wslogain = 4;
-
-        const Int_t newlast  = TMath::Nint(meanpulsetime+fExtractWinRight);
-
-        *fLog << all << underline 
-              << "Try to set new range limits: ("<<newfirst<<","<<newlast<<"+"<<wshigain
-              <<","<<newfirst-1<<","<<newlast<<"+"<<wslogain<<")"<<endl;
-        //
-        // Check the ranges for the new extractor setting
-        //
-        if (newfirst < 0)
-          {
-            *fLog << err << "Pulse is too much to the left, cannot go below 0! " << endl;
-            return kFALSE;
-          }
-        if (newlast+wshigain > numhigainsamples+numlogainsamples-1)
-          {
-            *fLog << err << "Pulse is too much to the right, cannot go beyond limits: " 
-                  << numhigainsamples << "+" << numlogainsamples << "-1" << endl;
-            *fLog << " Cannot extract at all! ... " << endl;
-            return kFALSE;
-          }
-        if (newlast+wslogain > numlogainsamples)
-          {
-            *fLog << err << "Pulse is too much to the right, cannot go beyond logain limits! " << endl;
-            *fLog << endl;
-            *fLog << "Try to use a different extractor (e.g. with a window size of only 4 sl.) or:" << endl;
-            *fLog << "Set the limit to a lower value (callisto.rc: line 329): " << endl;
-            *fLog << "    MJPedestalY2:ExtractWinRight: 4.5           " << endl;
-            *fLog << "(ATTENTION, you will lose late cosmics pulses!)" << endl;
-            *fLog << endl;
-            return kFALSE;
-          }
-        //
-        // Set and store the new ranges 
-        //
-        fExtractor->SetRange(newfirst,newlast+wshigain,
-                             newfirst>0?newfirst-1:newfirst,numlogainsamples-1);
-      }
-    
+    if (!PulsePosCheck(plist))
+        return kFALSE;
 
     *fLog << all << GetDescriptor() << ": Done." << endl;

trunk/MagicSoft/Mars/mjobs/MJPedestal.h

@@ -26 +26 @@
 {
 private:
-
     static const TString  fgReferenceFile;    //! default for fReferenceFile ("pedestalref.rc")
     static const TString  fgBadPixelsFile;    //! default for fBadPixelsFile ("badpixels_0_559.rc")
@@ -83 +82 @@
     void   DisplayOutliers(TH1D *hist) const;
     void   FixDataCheckHist(TH1D *hist) const;
+    Bool_t PulsePosCheck(const MParList &plist) const;
 
     Bool_t CheckEnvLocal();