Ignore:
Timestamp:
03/20/04 00:30:10 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/macros/calibration.C

    r3551 r3566  
    2727{
    2828
    29   //const TString inpath = "/mnt/Data/rootdata/CrabNebula/2004_02_10/";
    30   //const TString inpath = "/mnt/Data/rootdata/CrabNebula/2004_01_26/";
    3129  const TString inpath = "/mnt/Data/rootdata/CrabNebula/2004_01_27/";
    32   //const TString inpath = "/mnt/Data/rootdata/Miscellaneous/2003_12_19/";
    33   //  const TString inpath = "/mnt/Data/rootdata/Miscellaneous/2004_03_03/";
    34   //    const TString inpath = "/home/rootdata/BlindPixel/";
    35 
    36   //  const TString inpath = ".";
     30
    3731  MRunIter pruns;
    3832  MRunIter cruns;
    3933
    40   //  pruns.AddRun(20491,inpath);
    41   //  pruns.AddRun(20123,inpath);
    42 
    4334  pruns.AddRun(12386,inpath);
    44   //    pruns.AddRun(14400,inpath);
    45   //    cruns.AddRuns(14401,14409,inpath);
    46   //    pruns.AddRuns(15487,15494,inpath);
    47  //  cruns.AddRuns(12525,12529,inpath);
    4835  cruns.AddRun(12525,inpath);
    49   //    cruns.AddRuns(15495,15508,inpath);
    50   //  cruns.AddRun(20491,inpath);
    51   //  cruns.AddRun(20129,inpath);
    5236
    5337  gStyle->SetOptStat(1111);
     
    7761  //
    7862  /*
    79     MJCalibration calloop;
    80     calloop.SetInput(&cruns);
    81     calloop.SetDisplay(display);
    82     if (!calloop.Process(pedloop.GetPedestalCam()))
    83         return;
    84         //        #if 0
     63  MJCalibration calloop;
     64  calloop.SetInput(&cruns);
     65  calloop.SetDisplay(display);
     66  calloop.SetBadPixels(badcam);
     67  if (!calloop.Process(pedloop.GetPedestalCam()))
     68    return;
     69    #if 0
    8570  */
    8671  //
     
    206191  // print the most important results of all pixels to a file
    207192  //
     193  /*
    208194  MLog gauglog;
    209195  gauglog.SetOutputFile(Form("%s%s",calcam.GetName(),".txt"),1);
    210196  calcam.SetLogStream(&gauglog);
    211197  badcam.Print();
    212   //    calcam.Print();
    213198  calcam.SetLogStream(&gLog);
     199  */
    214200  //
    215201  // just one example how to get the plots of individual pixels
    216202  //
    217203  //    histblind.DrawClone("all");
    218   //histcharge[5].DrawClone("all");
     204  //    histcharge[5].DrawClone("all");
    219205  //    histcharge(5).DrawClone("all");
    220206  //    histtime[5].DrawClone("fourierevents");
     
    226212  MHCamera disp4  (geomcam, "Cal;RSigma",         "Reduced Sigmas");
    227213  MHCamera disp5  (geomcam, "Cal;RSigma/Charge",  "Reduced Sigma per Charge");
    228   MHCamera disp6  (geomcam, "Cal;FFactorPh",      "Nr. of Photo-electrons (F-Factor Method)");
     214  MHCamera disp6  (geomcam, "Cal;FFactorPhe",     "Nr. of Photo-electrons (F-Factor Method)");
    229215  MHCamera disp7  (geomcam, "Cal;FFactorConv",    "Conversion Factor to photons (F-Factor Method)");
    230216  MHCamera disp8  (geomcam, "Cal;FFactorFFactor", "Total F-Factor (F-Factor Method)");
     
    254240  MHCamera disp30 (geomcam, "time;Oscillating", "Oscillating Pixels");
    255241 
    256   MHCamera disp31 (geomcam, "Cal;AbsTimeMean", "Abs. Arrival Times");
    257   MHCamera disp32 (geomcam, "Cal;AbsTimeRms",  "RMS of Arrival Times");
     242  MHCamera disp31 (geomcam, "Cal;AbsTimeMean",  "Abs. Arrival Times");
     243  MHCamera disp32 (geomcam, "Cal;AbsTimeRms",   "RMS of Arrival Times");
    258244 
    259245  // Fitted charge means and sigmas
     
    294280  // Pixels with defects
    295281  disp13.SetCamContent(calcam,26);
    296   disp14.SetCamContent(calcam,27);
     282  disp14.SetCamContent(badcam,7);
    297283  disp15.SetCamContent(badcam,9);
    298284  disp16.SetCamContent(badcam,16);
    299285  disp17.SetCamContent(badcam,15);
    300   disp18.SetCamContent(calcam,28);
    301   disp19.SetCamContent(calcam,29);
     286  disp18.SetCamContent(calcam,27);
     287  disp19.SetCamContent(calcam,28);
    302288 
    303289  // Lo Gain calibration
    304   disp20.SetCamContent(calcam,30);
     290  disp20.SetCamContent(calcam,29);
    305291 
    306292  // Valid flags
     
    310296 
    311297  // Pedestals
    312   disp24.SetCamContent(calcam,31);
    313   disp24.SetCamError(  calcam,32);
    314   disp25.SetCamContent(calcam,33);
    315   disp25.SetCamError(  calcam,34);
     298  disp24.SetCamContent(calcam,30);
     299  disp24.SetCamError(  calcam,31);
     300  disp25.SetCamContent(calcam,32);
     301  disp25.SetCamError(  calcam,33);
    316302 
    317303  // Relative Times
     
    325311 
    326312  // Absolute Times
    327   disp31.SetCamContent(calcam,35);
     313  disp31.SetCamContent(calcam,34);
    328314  disp31.SetCamError(  calcam,35);
    329   disp32.SetCamContent(calcam,36);
     315  disp32.SetCamContent(calcam,35);
    330316 
    331317  disp1.SetYTitle("Charge [FADC units]");
     
    472458  CamDraw(c13,disp31,calcam,1,2,2);
    473459  CamDraw(c13,disp32,calcam,2,2,2);
     460#endif
    474461
    475462}
     
    516503  if (fit != 0)
    517504    {
    518 
    519505      c.cd(i+2*j);
    520506      gPad->SetBorderMode(0);
     
    541527          hprof[k]->Fit("pol1","Q","",min,max);
    542528          hprof[k]->GetFunction("pol1")->SetLineColor(kRed+k);
     529          hprof[k]->GetFunction("pol1")->SetLineWidth(1);
    543530        }
    544531     
Note: See TracChangeset for help on using the changeset viewer.