Changeset 8355 for trunk/MagicSoft/Mars/mhcalib/MHCalibrationPulseTimeCam.cc

Timestamp:

03/03/07 22:01:13 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhcalib/MHCalibrationPulseTimeCam.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/mhcalib/MHCalibrationPulseTimeCam.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MHCalibrationPulseTimeCam.cc,v 1.32 2007-03-01 23:43:33 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MHCalibrationPulseTimeCam.cc,v 1.33 2007-03-03 22:01:13 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -105 +105 @@
 
 #include <TEnv.h>
+#include <TLine.h>
 #include <TGraph.h>
 #include <TLegend.h>
@@ -218 +219 @@
   cam->fLast                   = fLast;
   cam->fReferenceFile          = fReferenceFile;
+  cam->fInnerRefTime           = fInnerRefTime;
+  cam->fOuterRefTime           = fOuterRefTime;
 
   if (!IsAverageing())
@@ -301 +304 @@
         fHiGainArray->AddAt(new MHCalibrationPix(Form("%sHiGainPix%04d",fHistName.Data(),i),
                                                  Form("%s High Gain Pixel %4d",fHistTitle.Data(),i)),i);
-        
-        MHCalibrationPix &pix = (*this)[i];          
+
+        MHCalibrationPix &pix = (*this)[i];
         pix.SetBinning(fNbins, fFirst, fLast);
-        
+
         MBadPixelsPix &bad = fIntensBad ? (*fIntensBad)[i] : (*fBadPixels)[i];
         InitHists(pix,bad,i);
@@ -322 +325 @@
         fAverageHiGainAreas->AddAt(new MHCalibrationPix(Form("%sHiGainArea%d",fHistName.Data(),j),
                                                         Form("%s High Gain Area Idx %d",fHistTitle.Data(),j)),j);
-        
+
         MHCalibrationPix &pix = GetAverageHiGainArea(j);
         pix.SetBinning(fNbins, fFirst, fLast);
@@ -338 +341 @@
           MHCalibrationPix &pix = GetAverageHiGainSector(j);
           pix.SetBinning(fNbins, fFirst, fLast);
-          
+
           InitHists(pix,fIntensCam ? fIntensCam->GetAverageBadSector(j) : fCam->GetAverageBadSector(j),j);
       }
@@ -469 +472 @@
   for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
     {
-      
+
       MHCalibrationPix &hist = (*this)[i];
-      
+
       if (hist.IsExcluded())
         continue;
-      
+
       MCalibrationPix &pix    = calcam[i];
       CalcHists(hist,pix);
@@ -481 +484 @@
   if (!IsAverageing())
     return kTRUE;
-  
+
   for (Int_t j=0; j<fAverageHiGainAreas->GetSize(); j++)
     {
-      
       MHCalibrationPix &hist = GetAverageHiGainArea(j);
       MCalibrationPix  &pix  = calcam.GetAverageArea(j);
       CalcHists(hist,pix);
-  }
-  
+    }
+
   for (Int_t j=0; j<fAverageHiGainSectors->GetSize(); j++)
     {
@@ -496 +498 @@
       CalcHists(hist,pix);
     }
-  
+
   return kTRUE;
 }
@@ -570 +572 @@
   option.ToLower();
 
-  TVirtualPad *pad = gPad ? gPad : MH::MakeDefCanvas(this);  
+  TVirtualPad *pad = gPad ? gPad : MH::MakeDefCanvas(this);
   pad->SetBorderMode(0);
-  pad->Divide(1,nareas);
+  pad->Divide(1,nareas, 1e-10, 1e-10);
 
   //
   // Loop over inner and outer pixels
-  //  
-  for (Int_t i=0; i<nareas;i++) 
+  //
+  for (Int_t i=0; i<nareas;i++)
      {
-       
        pad->cd(i+1);
-       
+
        MHCalibrationPix &hipix = GetAverageHiGainArea(i);
        DrawDataCheckPixel(hipix,i ? fOuterRefTime : fInnerRefTime);
-    }      
+     }
 }
 
@@ -597 +598 @@
 void MHCalibrationPulseTimeCam::DrawDataCheckPixel(MHCalibrationPix &pix, const Float_t refline)
 {
-  
-  gPad->SetTicks();
-
-  TH1F *hist = pix.GetHGausHist();
-
-  TH1F *null = new TH1F("Null",hist->GetTitle(),100,
-                        pix.GetFirst() > 0. ? pix.GetFirst() : 0.,
-                        pix.GetLast() > pix.GetFirst()
-                        ? pix.GetLast() : pix.GetFirst()*2.);
-
-  null->SetMaximum(1.1*hist->GetMaximum());
-  null->SetDirectory(NULL);
-  null->SetBit(kCanDelete);
-  null->SetStats(kFALSE);
-  //
-  // set the labels bigger
-  //
-  TAxis *xaxe = null->GetXaxis();
-  TAxis *yaxe = null->GetYaxis();
-  xaxe->CenterTitle();
-  yaxe->CenterTitle();    
-  xaxe->SetTitleSize(0.06);
-  yaxe->SetTitleSize(0.076);    
-  xaxe->SetTitleOffset(0.6);
-  yaxe->SetTitleOffset(0.65);    
-  xaxe->SetLabelSize(0.06);
-  yaxe->SetLabelSize(0.06);    
-  xaxe->SetTitle(hist->GetXaxis()->GetTitle());
-  yaxe->SetTitle(hist->GetYaxis()->GetTitle());  
-
-  null->Draw();
-  hist->Draw("same");
-
-  DisplayRefLines(null,refline);
-
-  return;
-  
+    gPad->SetBorderMode(0);
+    gPad->SetTicks();
+
+    TH1F *hist = pix.GetHGausHist();
+
+    //
+    // set the labels bigger
+    //
+    TAxis *xaxe = hist->GetXaxis();
+    TAxis *yaxe = hist->GetYaxis();
+    xaxe->CenterTitle();
+    yaxe->CenterTitle();
+    xaxe->SetTitleSize(0.06);
+    yaxe->SetTitleSize(0.076);
+    xaxe->SetTitleOffset(0.6);
+    yaxe->SetTitleOffset(0.65);
+    xaxe->SetLabelSize(0.06);
+    yaxe->SetLabelSize(0.06);
+    xaxe->SetTitle(hist->GetXaxis()->GetTitle());
+    yaxe->SetTitle(hist->GetYaxis()->GetTitle());
+    xaxe->SetRange(hist->GetMaximumBin()-30, hist->GetMaximumBin()+30);
+
+    hist->Draw();
+
+    DisplayRefLines(hist, refline);
 }
 
 void  MHCalibrationPulseTimeCam::DisplayRefLines(const TH1F *hist, const Float_t refline) const
 {
-
-  TGraph *gr    = new TGraph(2);
-  gr->SetPoint(0,refline,0.);
-  gr->SetPoint(1,refline,hist->GetMaximum());
-  gr->SetBit(kCanDelete);
-  gr->SetLineColor(106);
-  gr->SetLineStyle(2);
-  gr->SetLineWidth(3);
-  gr->Draw("L");
-
-  TLegend *leg = new TLegend(0.8,0.35,0.99,0.65);
-  leg->SetBit(kCanDelete);
-  leg->AddEntry(gr,"Reference","l");
-
-  leg->Draw();
+    TLine *line = new TLine(refline, 0, refline, hist->GetMaximum());
+    line->SetLineColor(106);
+    line->SetLineStyle(2);
+    line->SetLineWidth(3);
+    line->SetBit(kCanDelete);
+    line->Draw();
+
+    TLegend *leg = new TLegend(0.8,0.35,0.99,0.65);
+    leg->AddEntry(line, "Reference", "l");
+    leg->SetBit(kCanDelete);
+    leg->Draw();
 }