Changeset 8356


Ignore:
Timestamp:
03/03/07 22:06:24 (18 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mhcalib/MHCalibrationRelTimeCam.cc

    r8339 r8356  
    373373
    374374      histhi.AddSaturated(1);
    375 
     375/*
    376376      // If lo-gain arrival time has been extracted successfully,
    377377      // the hi-gain has saturateed and the lo-gain is switched on
     
    391391          fSumsectorlo[sector] += time;
    392392          fNumsectorlo[sector] ++;
    393       }
     393      }*/
    394394  }
    395395
     
    737737void MHCalibrationRelTimeCam::DrawDataCheckPixel(MHCalibrationPix &pix, const Float_t refline)
    738738{
    739 
    740739    if (pix.IsEmpty())
    741740        return;
     741
     742    TVirtualPad *pad = gPad;
     743    pad->Divide(1,2, 1e-10, 1e-10);
     744    pad->cd(1);
     745
     746    gPad->SetBorderMode(0);
     747
     748    gPad->SetTicks();
     749    if (!pix.IsEmpty() && !pix.IsOnlyOverflow() && !pix.IsOnlyUnderflow())
     750        gPad->SetLogy();
     751
     752    TH1F *hist = pix.GetHGausHist();
    742753 
    743   TVirtualPad *newpad = gPad;
    744   newpad->Divide(1,2);
    745   newpad->cd(1);
    746  
    747   gPad->SetTicks();
    748   if (!pix.IsEmpty() && !pix.IsOnlyOverflow() && !pix.IsOnlyUnderflow())
    749     gPad->SetLogy();
    750 
    751   TH1F *hist = pix.GetHGausHist();
    752  
    753   TH1F *null = new TH1F("Null",hist->GetTitle(),100,0.,pix.GetLast());
    754 
    755   null->SetMaximum(1.1*hist->GetMaximum());
    756   null->SetDirectory(NULL);
    757   null->SetBit(kCanDelete);
    758   null->SetStats(kFALSE);
    759   //
    760   // set the labels bigger
    761   //
    762   TAxis *xaxe = null->GetXaxis();
    763   TAxis *yaxe = null->GetYaxis();
    764   xaxe->CenterTitle();
    765   yaxe->CenterTitle();   
    766   xaxe->SetTitleSize(0.07);
    767   yaxe->SetTitleSize(0.07);   
    768   xaxe->SetTitleOffset(0.65);
    769   yaxe->SetTitleOffset(0.55);   
    770   xaxe->SetLabelSize(0.06);
    771   yaxe->SetLabelSize(0.06);   
    772   xaxe->SetTitle(hist->GetXaxis()->GetTitle());
    773   yaxe->SetTitle(hist->GetYaxis()->GetTitle()); 
    774 
    775   null->Draw();
    776   hist->Draw("same");
    777 
    778   gStyle->SetOptFit();
    779 
    780   TF1 *fit = pix.GetFGausFit();
    781 
    782   if (fit)
    783   {
    784     switch ( fColor )
    785       {
    786       case MCalibrationCam::kGREEN:
    787         fit->SetLineColor(kGreen);
    788         break;
    789       case MCalibrationCam::kBLUE:
    790         fit->SetLineColor(kBlue);
    791         break;
    792       case MCalibrationCam::kUV: 
    793         fit->SetLineColor(106);
    794         break;
    795       case MCalibrationCam::kCT1:
    796         fit->SetLineColor(006);
    797         break;
    798       default:
    799         fit->SetLineColor(kRed);
    800       }
    801     fit->Draw("same");
    802   }
    803 
    804   DisplayRefLines(null,refline);
    805 
    806   newpad->cd(2);
    807   gPad->SetTicks();
    808 
    809   TH1F *null2 = new TH1F("Null2",hist->GetTitle(),100,0.,pix.GetEvents()->GetSize()/pix.GetEventFrequency());
    810 
    811   null2->SetMinimum(pix.GetMean()-10.*pix.GetSigma());
    812   null2->SetMaximum(pix.GetMean()+10.*pix.GetSigma());
    813   null2->SetDirectory(NULL);
    814   null2->SetBit(kCanDelete);
    815   null2->SetStats(kFALSE);
    816   //
    817   // set the labels bigger
    818   //
    819   TAxis *xaxe2 = null2->GetXaxis();
    820   TAxis *yaxe2 = null2->GetYaxis();
    821   xaxe2->CenterTitle();
    822   yaxe2->CenterTitle();   
    823   xaxe2->SetTitleSize(0.07);
    824   yaxe2->SetTitleSize(0.07);   
    825   xaxe2->SetTitleOffset(0.65);
    826   yaxe2->SetTitleOffset(0.55);   
    827   xaxe2->SetLabelSize(0.06);
    828   yaxe2->SetLabelSize(0.06);   
    829  
    830   pix.CreateGraphEvents();
    831   TGraph *gr = pix.GetGraphEvents();
    832   if (gr)
    833   {
    834       xaxe2->SetTitle(gr->GetXaxis()->GetTitle());
    835       yaxe2->SetTitle(gr->GetYaxis()->GetTitle());
    836   }
    837 
    838   null2->Draw();
    839 
    840   pix.DrawEvents("same");
    841 
    842   return;
     754    //
     755    // set the labels bigger
     756    //
     757    TAxis *xaxe = hist->GetXaxis();
     758    TAxis *yaxe = hist->GetYaxis();
     759    xaxe->CenterTitle();
     760    yaxe->CenterTitle();
     761    xaxe->SetTitleSize(0.07);
     762    yaxe->SetTitleSize(0.07);
     763    xaxe->SetTitleOffset(0.65);
     764    yaxe->SetTitleOffset(0.55);
     765    xaxe->SetLabelSize(0.06);
     766    yaxe->SetLabelSize(0.06);
     767    xaxe->SetTitle(hist->GetXaxis()->GetTitle());
     768    yaxe->SetTitle(hist->GetYaxis()->GetTitle());
     769    xaxe->SetRange(hist->GetMaximumBin()-1500, hist->GetMaximumBin()+1500);
     770
     771    hist->Draw();
     772
     773    gStyle->SetOptFit();
     774
     775    TF1 *fit = pix.GetFGausFit();
     776
     777    if (fit)
     778    {
     779        switch (fColor)
     780        {
     781        case MCalibrationCam::kGREEN:
     782            fit->SetLineColor(kGreen);
     783            break;
     784        case MCalibrationCam::kBLUE:
     785            fit->SetLineColor(kBlue);
     786            break;
     787        case MCalibrationCam::kUV:
     788            fit->SetLineColor(106);
     789            break;
     790        case MCalibrationCam::kCT1:
     791            fit->SetLineColor(006);
     792            break;
     793        default:
     794            fit->SetLineColor(kRed);
     795        }
     796        fit->Draw("same");
     797    }
     798
     799    DisplayRefLines(hist, refline);
     800
     801    pad->cd(2);
     802    gPad->SetBorderMode(0);
     803    gPad->SetTicks();
     804
     805    pix.CreateGraphEvents();
     806    TGraph *gr = pix.GetGraphEvents();
     807    if (gr)
     808    {
     809        TH1F *null2 = gr->GetHistogram();
     810
     811        null2->SetMinimum(pix.GetMean()-10.*pix.GetSigma());
     812        null2->SetMaximum(pix.GetMean()+10.*pix.GetSigma());
     813
     814        //
     815        // set the labels bigger
     816        //
     817        TAxis *xaxe2 = null2->GetXaxis();
     818        TAxis *yaxe2 = null2->GetYaxis();
     819        xaxe2->CenterTitle();
     820        yaxe2->CenterTitle();
     821        xaxe2->SetTitleSize(0.07);
     822        yaxe2->SetTitleSize(0.07);
     823        xaxe2->SetTitleOffset(0.65);
     824        yaxe2->SetTitleOffset(0.55);
     825        xaxe2->SetLabelSize(0.06);
     826        yaxe2->SetLabelSize(0.06);
     827        xaxe2->SetRangeUser(0.,pix.GetEvents()->GetSize()/pix.GetEventFrequency());
     828    }
     829
     830    pix.DrawEvents();
    843831}
    844832
    845833void  MHCalibrationRelTimeCam::DisplayRefLines(const TH1F *hist, const Float_t refline) const
    846834{
    847 
    848   TGraph *gr = new TGraph(2);
    849   gr->SetPoint(0,refline,0.);
    850   gr->SetPoint(1,refline,hist->GetMaximum());
    851   gr->SetBit(kCanDelete);
    852   gr->SetLineColor(kGreen);
    853   gr->SetLineStyle(2);
    854   gr->SetLineWidth(3);
    855   gr->Draw("L");
    856 
    857   TLegend *leg = new TLegend(0.75,0.7,0.99,0.99);
    858   leg->SetBit(kCanDelete);
    859   leg->AddEntry(gr,"Trigger Calibration","l");
    860   leg->Draw();
     835    TLine *line = new TLine(refline, 0, refline, hist->GetMaximum());
     836    line->SetLineColor(kGreen);
     837    line->SetLineStyle(2);
     838    line->SetLineWidth(3);
     839    line->SetBit(kCanDelete);
     840    line->Draw();
     841
     842    TLegend *leg = new TLegend(0.75,0.01,0.99,0.3);
     843    leg->SetBit(kCanDelete);
     844    leg->AddEntry(line, "Trigger Calibration", "l");
     845    leg->Draw();
    861846}
    862847
Note: See TracChangeset for help on using the changeset viewer.