Changeset 3360 for trunk

Timestamp:

02/29/04 04:21:54 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• macros/calibration.C (modified) (12 diffs)
• mjobs/MJCalibration.cc (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -7 +7 @@
    
    * mjobs/MJPedestal.cc
+   * mjobs/MJCalibration.cc
+   * macros/calibration.C
     
      - put in call to MHCamera::Projection the name of teh calling 

trunk/MagicSoft/Mars/macros/calibration.C

@@ -28 +28 @@
 {
 
-    const char *inpath = "/remote/home/pc2/operator/Crab20040214/";
+//    const char *inpath = "/remote/home/pc2/operator/NewCalBoxTest/";
+//    const char *inpath = "/remote/home/pc2/operator/Crab20040214/";
+    const char *inpath = "./";
     //const TString inpath = "/mnt/Data/rootdata/CrabNebula/2004_02_10/";
     //const TString inpath = "/mnt/Data/rootdata/CrabNebula/2004_01_26/";
@@ -36 +38 @@
     MRunIter cruns;
 
-    pruns.AddRuns(15487,15494,inpath);
-    cruns.AddRuns(15495,15508,inpath);
+    pruns.AddRun(14400,inpath);
+    cruns.AddRuns(14401,14409,inpath);
+//    pruns.AddRuns(15487,15494,inpath);
+//    cruns.AddRuns(15495,15508,inpath);
+//    pruns.AddRun(19816,inpath);  // CalboxTest
+//    cruns.AddRun(19820,inpath); // 3blue
+//    cruns.AddRun(19820,inpath);   // 1blue
 
     gStyle->SetOptStat(1111);
@@ -43 +50 @@
 
     MStatusDisplay *display = new MStatusDisplay;
-    display->SetUpdateTime(500);
+    display->SetUpdateTime(3000);
     display->Resize(850,700);
 
@@ -95 +102 @@
     MHCalibrationChargePINDiode histpin;
     MHCalibrationChargeBlindPix histblind;
+//    histblind.SetSinglePheCut(80);
     //
     // As long, as we don't have digital modules,
@@ -110 +118 @@
     plist.AddToList(&histtime);
     plist.AddToList(&histcharge);
-    plist.AddToList(&histpin);
+//    plist.AddToList(&histpin);
     plist.AddToList(&histblind);
 
@@ -120 +128 @@
     MExtractPINDiode       pincalc;
     MExtractBlindPixel     blindcalc;
+//    blindcalc.SetRange(11,16);
+
     MArrivalTimeCalc2      timecalc;
     MCalibrationChargeCalc calcalc;
@@ -125 +135 @@
     
     MFillH filltime( "MHCalibrationRelTimeCam"    , "MArrivalTimeCam");
-    MFillH fillpin  ("MHCalibrationChargePINDiode", "MExtractedSignalPINDiode");
+//   MFillH fillpin  ("MHCalibrationChargePINDiode", "MExtractedSignalPINDiode");
     MFillH fillblind("MHCalibrationChargeBlindPix", "MExtractedSignalBlindPixel");
     MFillH fillcam  ("MHCalibrationChargeCam"     , "MExtractedSignalCam");
@@ -139 +149 @@
     //
     MFCosmics            cosmics;
+    cosmics.SetMaxEmptyPixels(0.4);
     MContinue            cont(&cosmics);
 
@@ -145 +156 @@
     tlist.AddToList(&sigcalc2);
     tlist.AddToList(&blindcalc);
-    tlist.AddToList(&pincalc);
+//    tlist.AddToList(&pincalc);
     //
     // In case, you want to skip the cosmics rejection, 
@@ -189 +200 @@
     //
     histblind.DrawClone("all");
-    histcharge[5].DrawClone("time");
+    histcharge[5].DrawClone("all");
+    histcharge(5).DrawClone("all");
     histtime[5].DrawClone("fourierevents");
 
@@ -469 +481 @@
       c.cd(i+2*j);
       gPad->SetBorderMode(0);
-      TH1D *obj2 = (TH1D*)obj1->Projection();
+      TH1D *obj2 = (TH1D*)obj1->Projection(obj1.GetName());
       
 //      obj2->Sumw2();
@@ -558 +570 @@
         }
       
+        // Just to get the right (maximum) binning
+        TH1D *half[2];
+        half[0] = obj1->Projection("Sector 6-1-2");
+        half[1] = obj1->Projection("Sector 3-4-5");
+
+        half[0]->Reset();
+        half[1]->Reset();
+
+        TH1D *dummy = obj1->Projection("Dummy");
+
+        for (int i=1; i<7/*obj1->GetGeomCam().GetNumSectors()*/; i++)
+            half[(i/3)%2]->Add(obj1->ProjectionS(i, "Dummy"));
+
+        for (int i=0; i<2; i++)
+        {
+            half[i]->SetLineColor(kRed+i);
+            half[i]->SetDirectory(0);
+            half[i]->SetBit(kCanDelete);
+            half[i]->Draw("same");
+        }
+
+        delete dummy;
+
       gPad->Modified();
       gPad->Update();

trunk/MagicSoft/Mars/mjobs/MJCalibration.cc

@@ -78 +78 @@
 {
 
-  TH1D *obj2 = (TH1D*)obj1->Projection();
+  TH1D *obj2 = (TH1D*)obj1->Projection(obj1->GetName());
   obj2->Draw();
   obj2->SetBit(kCanDelete);