Ignore:
Timestamp:
09/09/04 18:03:23 (20 years ago)
Author:
gaug
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/mjobs/MJCalibration.cc

    r4870 r4909  
    111111
    112112#include "MPedestalCam.h"
     113#include "MCalibColorSteer.h"
     114
     115#include "MCalibrationIntensityChargeCam.h"
     116#include "MCalibrationIntensityRelTimeCam.h"
     117#include "MCalibrationIntensityQECam.h"
     118
    113119#include "MCalibrationCam.h"
    114120#include "MCalibrationQECam.h"
    … …  
    122128#include "MCalibrationChargeBlindCamTwoNewStyle.h"
    123129#include "MCalibrationChargeCalc.h"
     130#include "MCalibColorSet.h"
    124131
    125132#include "MHGausEvents.h"
    … …  
    170177
    171178const Float_t  MJCalibration::fgRefConvFADC2PheInner      = 0.14;       
    172 const Float_t  MJCalibration::fgRefConvFADC2PheOuter      = 0.45;       
     179const Float_t  MJCalibration::fgRefConvFADC2PheOuter      = 0.4;       
    173180const Float_t  MJCalibration::fgRefConvFADC2PhotInner     = 0.8;
    174181const Float_t  MJCalibration::fgRefConvFADC2PhotOuter     = 3.8;
    … …  
    184191//
    185192// - Sets fRuns to 0, fExtractor to NULL, fTimeExtractor to NULL, fColor to kNONE,
    186 //   fDisplay to kNormalDisplay, fRelTime to kFALSE, fDataCheck to kFALSE,
     193//   fDisplay to kNormalDisplay, kRelTimes to kFALSE, kataCheck to kFALSE, kDebug to kFALSE,
     194//   kIntensity to kFALSE
    187195// - SetUseBlindPixel()
    188196// - SetUsePINDiode()
    … …  
    191199    : fRuns(0), fExtractor(NULL), fTimeExtractor(NULL),
    192200      fColor(MCalibrationCam::kNONE), fDisplayType(kNormalDisplay),
    193       fRelTimes(kFALSE), fDataCheck(kFALSE), fDebug(kFALSE), fGeometry("MGeomCamMagic")
     201      fGeometry("MGeomCamMagic")
    194202{
    195     fName  = name  ? name  : "MJCalibration";
    196     fTitle = title ? title : "Tool to create the calibration constants for one calibration run";
    197 
    198     SetUseBlindPixel();
    199     SetUsePINDiode();
     203
     204  fName  = name  ? name  : "MJCalibration";
     205  fTitle = title ? title : "Tool to create the calibration constants for one calibration run";
     206 
     207  SetUseBlindPixel();
     208  SetUsePINDiode();
     209 
     210  SetRelTimeCalibration( kFALSE );
     211  SetDataCheck         ( kFALSE );
     212  SetDebug             ( kFALSE );
     213  SetIntensity         ( kFALSE );
    200214}
    201215
    … …  
    297311    MHCamera disp34(geomcam, "TimeOscillating",   "Oscillating Pixels");
    298312    MHCamera disp35(geomcam, "TotalConv",         "Conversion Factor to photons");
    299 
     313   
     314    MCalibrationChargeCam  *cam    = NULL;
     315    MCalibrationQECam      *qecam  = NULL;
     316    MCalibrationRelTimeCam *relcam = NULL;
     317
     318    if (IsIntensity())
     319      {
     320        cam    = (MCalibrationChargeCam*) fIntensCalibCam.GetCam();
     321        qecam  = (MCalibrationQECam*)     fIntensQECam.GetCam();
     322        relcam = (MCalibrationRelTimeCam*)fIntensRelTimeCam.GetCam();
     323      }
     324    else
     325      {
     326        cam    = &fCalibrationCam;
     327        qecam  = &fQECam;
     328        relcam = &fRelTimeCam;
     329      }
     330   
    300331    // Fitted charge means and sigmas
    301     disp1.SetCamContent(fCalibrationCam,  0);
    302     disp1.SetCamError(  fCalibrationCam,  1);
    303     disp2.SetCamContent(fCalibrationCam,  2);
    304     disp2.SetCamError(  fCalibrationCam,  3);
     332    disp1.SetCamContent(*cam,  0);
     333    disp1.SetCamError(  *cam,  1);
     334    disp2.SetCamContent(*cam,  2);
     335    disp2.SetCamError(  *cam,  3);
    305336
    306337    // Reduced Sigmas and reduced sigmas per charge
    307     disp3.SetCamContent(fCalibrationCam,  5);
    308     disp3.SetCamError(  fCalibrationCam,  6);
    309     disp4.SetCamContent(fCalibrationCam,  7);
    310     disp4.SetCamError(  fCalibrationCam,  8);
     338    disp3.SetCamContent(*cam,  5);
     339    disp3.SetCamError(  *cam,  6);
     340    disp4.SetCamContent(*cam,  7);
     341    disp4.SetCamError(  *cam,  8);
    311342
    312343    // F-Factor Method
    313     disp5.SetCamContent(fCalibrationCam,  9);
    314     disp5.SetCamError(  fCalibrationCam, 10);
    315     disp6.SetCamContent(fCalibrationCam, 11);
    316     disp6.SetCamError(  fCalibrationCam, 12);
    317     disp7.SetCamContent(fCalibrationCam, 13);
    318     disp7.SetCamError(  fCalibrationCam, 14);
     344    disp5.SetCamContent(*cam,  9);
     345    disp5.SetCamError(  *cam, 10);
     346    disp6.SetCamContent(*cam, 11);
     347    disp6.SetCamError(  *cam, 12);
     348    disp7.SetCamContent(*cam, 13);
     349    disp7.SetCamError(  *cam, 14);
    319350
    320351    // Quantum Efficiencies
    321     disp8.SetCamContent (fQECam, 0 );
    322     disp8.SetCamError   (fQECam, 1 );
    323     disp9.SetCamContent (fQECam, 2 );
    324     disp9.SetCamError   (fQECam, 3 );
    325     disp10.SetCamContent(fQECam, 4 );
    326     disp10.SetCamError  (fQECam, 5 );
    327     disp11.SetCamContent(fQECam, 6 );
    328     disp11.SetCamError  (fQECam, 7 );
     352    disp8.SetCamContent (*qecam, 0 );
     353    disp8.SetCamError   (*qecam, 1 );
     354    disp9.SetCamContent (*qecam, 2 );
     355    disp9.SetCamError   (*qecam, 3 );
     356    disp10.SetCamContent(*qecam, 4 );
     357    disp10.SetCamError  (*qecam, 5 );
     358    disp11.SetCamContent(*qecam, 6 );
     359    disp11.SetCamError  (*qecam, 7 );
    329360
    330361    // Valid flags
    331     disp12.SetCamContent(fQECam, 8 );
    332     disp13.SetCamContent(fQECam, 9 );
    333     disp14.SetCamContent(fQECam, 10);
    334     disp15.SetCamContent(fQECam, 11);
     362    disp12.SetCamContent(*qecam, 8 );
     363    disp13.SetCamContent(*qecam, 9 );
     364    disp14.SetCamContent(*qecam, 10);
     365    disp15.SetCamContent(*qecam, 11);
    335366
    336367    // Conversion Hi-Lo
    337     disp16.SetCamContent(fCalibrationCam, 25);
    338     disp17.SetCamContent(fCalibrationCam, 16);
    339     disp17.SetCamError  (fCalibrationCam, 17);
    340     disp18.SetCamContent(fCalibrationCam, 18);
    341     disp18.SetCamError  (fCalibrationCam, 19);
     368    disp16.SetCamContent(*cam, 25);
     369    disp17.SetCamContent(*cam, 16);
     370    disp17.SetCamError  (*cam, 17);
     371    disp18.SetCamContent(*cam, 18);
     372    disp18.SetCamError  (*cam, 19);
    342373
    343374    // Pickup and Blackout
    344     disp19.SetCamContent(fCalibrationCam, 21);
    345     disp20.SetCamContent(fCalibrationCam, 22);
    346     disp21.SetCamContent(fCalibrationCam, 23);
    347     disp22.SetCamContent(fCalibrationCam, 24);
     375    disp19.SetCamContent(*cam, 21);
     376    disp20.SetCamContent(*cam, 22);
     377    disp21.SetCamContent(*cam, 23);
     378    disp22.SetCamContent(*cam, 24);
    348379
    349380    // Pixels with defects
    350     disp23.SetCamContent(fCalibrationCam, 20);
     381    disp23.SetCamContent(*cam, 20);
    351382    disp24.SetCamContent(fBadPixels, 6);
    352383    disp25.SetCamContent(fBadPixels, 7);
    … …  
    357388
    358389    // Arrival Times
    359     disp28.SetCamContent(fCalibrationCam, 26);
    360     disp28.SetCamError(  fCalibrationCam, 27);
    361     disp29.SetCamContent(fCalibrationCam, 27);
     390    disp28.SetCamContent(*cam, 26);
     391    disp28.SetCamError(  *cam, 27);
     392    disp29.SetCamContent(*cam, 27);
    362393
    363394    disp1.SetYTitle("Q [FADC counts]");
    … …  
    403434      {
    404435
    405         MCalibrationChargePix &pix = (MCalibrationChargePix&)fCalibrationCam[i];
    406         MCalibrationQEPix     &qe  = (MCalibrationQEPix&)    fQECam         [i];
     436        MCalibrationChargePix &pix = (MCalibrationChargePix&)(*cam)  [i];
     437        MCalibrationQEPix     &qe  = (MCalibrationQEPix&)    (*qecam)[i];
    407438
    408439        if (!pix.IsFFactorMethodValid())
    … …  
    417448   
    418449
    419     if (fRelTimes)
    420     {
    421         disp30.SetCamContent(fRelTimeCam,0);
    422         disp30.SetCamError(  fRelTimeCam,1);
    423         disp31.SetCamContent(fRelTimeCam,2);
    424         disp31.SetCamError(  fRelTimeCam,3);
    425         disp32.SetCamContent(fRelTimeCam,4);
     450    if (IsRelTimes())
     451    {
     452        disp30.SetCamContent(*relcam,0);
     453        disp30.SetCamError(  *relcam,1);
     454        disp31.SetCamContent(*relcam,2);
     455        disp31.SetCamError(  *relcam,3);
     456        disp32.SetCamContent(*relcam,4);
    426457        disp33.SetCamContent(fBadPixels,20);
    427458        disp34.SetCamContent(fBadPixels,21);
    … …  
    434465    }
    435466
    436 
    437    
    438467    if (fDisplayType == kDataCheckDisplay)
    439468    {
    … …  
    739768      if (geomcam.InheritsFrom("MGeomCamMagic"))
    740769        DisplayDoubleProject(&disp29, "", "jittering");
    741 
    742770
    743771      //
    … …  
    884912      obj9->Draw();
    885913
    886       if (fRelTimes)
     914      if (IsRelTimes())
    887915        {
    888916          // Rel. Times
    … …  
    955983        disp29.CamDraw(c16, 2, 2, 5);
    956984
    957         if (fRelTimes)
     985        if (IsRelTimes())
    958986        {
    959987            // Rel. Times
    … …  
    10521080        disp29.CamDraw(c30, 2, 2, 1);
    10531081
    1054         if (fRelTimes)
     1082        if (IsRelTimes())
    10551083        {
    10561084            // Rel. Times
    … …  
    11421170{
    11431171
    1144   const Float_t mean  = hist->GetFunction("gaus")->GetParameter(1);
    1145   const Float_t lolim = mean - 4.0*hist->GetFunction("gaus")->GetParameter(2);
    1146   const Float_t uplim = mean + 4.0*hist->GetFunction("gaus")->GetParameter(2);
     1172  const Int_t kNotDraw = 1<<9;
     1173  TF1 *f = hist->GetFunction("gaus");
     1174  f->ResetBit(kNotDraw);
     1175
     1176  const Float_t mean  = f->GetParameter(1);
     1177  const Float_t lolim = mean - 4.0*f->GetParameter(2);
     1178  const Float_t uplim = mean + 4.0*f->GetParameter(2);
    11471179  const Stat_t  dead  = hist->Integral(0,hist->FindBin(lolim)-1);
    11481180  const Stat_t  noisy = hist->Integral(hist->FindBin(uplim)+1,hist->GetNbinsX()+1);
    … …  
    13561388        SetNormalDisplay();
    13571389
    1358     SetRelTimeCalibration(GetEnv("RelTimeCalibration", fRelTimes));
    1359     SetDataCheck(GetEnv("Datacheck", fDataCheck));
    1360     SetDebug(GetEnv("Debug", fDebug));
     1390    SetRelTimeCalibration(GetEnv("RelTimeCalibration", IsRelTimes()));
     1391    SetDataCheck(GetEnv("Datacheck", IsDataCheck()));
     1392    SetDebug(GetEnv("Debug", IsDebug()));
     1393    SetIntensity(GetEnv("Intensity", IsIntensity()));
     1394
    13611395    SetUseBlindPixel(GetEnv("UseBlindPixel", IsUseBlindPixel()));
    13621396    SetUsePINDiode(GetEnv("UsePINDiode", IsUsePINDiode()));
    … …  
    14381472//   2) MCalibrationQECam     (fQECam)
    14391473//   3) MCalibrationChargeCam (fCalibrationCam)
    1440 //   4) MCalibrationRelTimeCam (fRelTimeCam)   (only if flag fRelTimes is chosen)
     1474//   4) MCalibrationRelTimeCam (fRelTimeCam)   (only if flag IsRelTimes() is chosen)
    14411475//   5) MBadPixelsCam         (fBadPixels)
    14421476//   6) MCalibrationChargePINDiode
    … …  
    14491483//   5)  MExtractPINDiode
    14501484//   6)  MExtractBlindPixel
    1451 //   7)  MExtractTime (only if flag fRelTimes is chosen)
     1485//   7)  MExtractTime (only if flag IsRelTimes() is chosen)
    14521486//   8)  MContinue(MFCosmics)
    1453 //   9)  MFillH("MHCalibrationChargePINDiode", "MExtractedSignalPINDiode")
    1454 //   10) MFillH("MHCalibrationChargeBlindCam", "MExtractedSignalBlindPixel")
    1455 //   11) MFillH("MHCalibrationChargeCam",      "MExtractedSignalCam")
    1456 //   12) MFillH("MHCalibrationChargeCam",      "MExtractedSignalCam")
     1487//   9)  MFillH("MHCalibrationChargePINDiode", "MExtractedSignalPINDiode",   "FillPINDiode")
     1488//   10) MFillH("MHCalibrationChargeBlindCam", "MExtractedSignalBlindPixel", "FillBlindCam")
     1489//   11) MFillH("MHCalibrationChargeCam",      "MExtractedSignalCam",        "FillChargeCam")
     1490//   12) MFillH("MHCalibrationChargeCam",      "MExtractedSignalCam",        "FillRelTime")
    14571491//   13) MCalibrationChargeCalc
    1458 //   14) MFillH("MHCalibrationRelTimeCam",     "MArrivalTimeCam") (only if flag fRelTimes is chosen)
     1492//   14) MFillH("MHCalibrationRelTimeCam",     "MArrivalTimeCam") (only if flag IsRelTimes() is chosen)
    14591493//   15) MCalibrationRelTimeCalc
    14601494// - Execute MEvtLoop
    … …  
    15071541    plist.AddToList(&tlist);
    15081542    plist.AddToList(this); // take care of fDisplay!
    1509 
    1510     MReadMarsFile read("Events");
    1511     MRawFileRead rawread(NULL);
    15121543
    15131544    MDirIter iter;
    … …  
    15281559    }
    15291560
    1530     if (fDataCheck)
     1561    //
     1562    // Input containers
     1563    //
     1564    plist.AddToList(&pedcam);
     1565    plist.AddToList(&fBadPixels);
     1566
     1567    //
     1568    // Calibration Results containers
     1569    //
     1570    if (IsIntensity())
     1571      {
     1572        plist.AddToList(&fIntensQECam);
     1573        plist.AddToList(&fIntensCalibCam);
     1574        //        plist.AddToList(&fIntensCalibrationBlindCam);
     1575        //        plist.AddToList(&fIntensCalibrationPINDiode);
     1576        plist.AddToList(&fIntensRelTimeCam);
     1577      }
     1578    else
     1579      {
     1580        plist.AddToList(&fQECam);
     1581        plist.AddToList(&fCalibrationCam);
     1582        plist.AddToList(&fCalibrationBlindCam);
     1583        plist.AddToList(&fCalibrationPINDiode);
     1584        plist.AddToList(&fRelTimeCam);
     1585      }
     1586
     1587    //
     1588    // Initialize two histogram containers which could be modified in this class
     1589    //
     1590    MHCalibrationChargeCam      chargecam;
     1591    MHCalibrationChargeBlindCam blindcam;
     1592    plist.AddToList(&chargecam);
     1593    plist.AddToList(&blindcam);
     1594
     1595    //
     1596    // Data Reading tasks
     1597    //
     1598    MReadMarsFile read("Events");
     1599    MRawFileRead rawread(NULL);
     1600
     1601    if (IsDataCheck())
    15311602    {
    15321603        rawread.AddFiles(fSequence.IsValid() ? iter : *fRuns);
    … …  
    15401611    }
    15411612
    1542     MHCalibrationChargeCam      chargecam;
    1543     MHCalibrationChargeBlindCam blindcam;
    1544 
    1545     plist.AddToList(&pedcam);
    1546     plist.AddToList(&chargecam);
    1547     plist.AddToList(&blindcam);
    1548     plist.AddToList(&fBadPixels);
    1549     plist.AddToList(&fQECam);
    1550     plist.AddToList(&fCalibrationCam);
    1551     plist.AddToList(&fCalibrationBlindCam);
    1552     plist.AddToList(&fCalibrationPINDiode);
    1553     plist.AddToList(&fRelTimeCam);
    1554 
     1613    //
     1614    // Other Tasks
     1615    //
    15551616    MGeomApply               apply;
    15561617    apply.SetGeometry(fGeometry);
    … …  
    15721633    }
    15731634
    1574     if (fDebug)
     1635    if (IsDebug())
    15751636    {
    15761637        chargecam.SetDebug();
    … …  
    15791640
    15801641    //
    1581     // As long as there are no DM's, have to colour by hand
    1582     //
    1583     chargecam.SetColor    (fColor);
    1584     calcalc.SetPulserColor(fColor);
    1585 
     1642    // Calibration histogramming
     1643    //
    15861644    MFillH fillpin("MHCalibrationChargePINDiode", "MExtractedSignalPINDiode",   "FillPINDiode");
    15871645    MFillH fillbnd("MHCalibrationChargeBlindCam", "MExtractedSignalBlindPixel", "FillBlindCam");
    … …  
    16131671    MContinue cont(&cosmics);
    16141672
     1673    MTaskEnv taskenv("ExtractSignal");
     1674    taskenv.SetDefault(fExtractor ? fExtractor : &extract2);
     1675
    16151676    tlist.AddToList(&merge);
    16161677    tlist.AddToList(&apply);
    1617 
    1618     MTaskEnv taskenv("ExtractSignal");
    1619     taskenv.SetDefault(fExtractor ? fExtractor : &extract2);
    1620 
    16211678    tlist.AddToList(&taskenv);
    16221679    tlist.AddToList(&pinext);
    … …  
    16261683    taskenv2.SetDefault(fTimeExtractor ? fTimeExtractor : &timespline);
    16271684
    1628     if (fRelTimes)
     1685    if (IsRelTimes())
    16291686        tlist.AddToList(&taskenv2);
    16301687
     1688    //
     1689    // FIXME: This is not yet implemented in the classes!!!
     1690    //
    16311691    if (fColor == MCalibrationCam::kCT1)
    16321692        tlist.AddToList(&cont);
    16331693
     1694    MCalibColorSet colorset;
     1695    tlist.AddToList(&colorset);
     1696
     1697    MCalibColorSteer steer;
     1698    if (IsIntensity())
     1699      tlist.AddToList(&steer);
     1700
    16341701    tlist.AddToList(&fillcam);
    16351702
    1636     if (fRelTimes)
     1703    if (IsRelTimes())
    16371704    {
    16381705        tlist.AddToList(&filltme);
    … …  
    16621729
    16631730    if (!WriteTasks(taskenv.GetTask(), taskenv2.GetTask()))
    1664         return kFALSE;
     1731      return kFALSE;
    16651732
    16661733    // Execute first analysis
    … …  
    16931760    if (fDisplay)
    16941761    {
    1695         /*
    16961762      if (IsUsePINDiode())
    16971763        {
    … …  
    17061772          cam->DrawClone(Form("nonew %s",drawoption.Data()));
    17071773        }
    1708       else if (fRelTimes)
     1774      else if (IsRelTimes())
    17091775        {
    17101776          MHCalibrationRelTimeCam *cam =
    … …  
    17131779        }
    17141780      else
    1715         {
    1716           MHCalibrationChargeCam *cam =
    1717             (MHCalibrationChargeCam*)plist.FindObject("MHCalibrationChargeCam");
    1718           cam->DrawClone(Form("nonew %s",drawoption.Data()));
    1719         } */
     1781        chargecam.DrawClone(Form("nonew %s",drawoption.Data()));
    17201782    }
    17211783
    … …  
    18961958    *fLog << inf << "ok." << endl;
    18971959
    1898     *fLog << inf << " - MCalibrationChargeCam..." << flush;
    1899     if (fCalibrationCam.Write()<=0)
    1900     {
    1901         *fLog << err << "Unable to write MCalibrationChargeCam to " << oname << endl;
    1902         return kFALSE;
    1903     }
    1904     *fLog << inf << "ok." << endl;
    1905 
     1960    if (IsIntensity())
     1961      {
     1962        *fLog << inf << " - MCalibrationIntensityChargeCam..." << flush;
     1963        if (fIntensCalibCam.Write()<=0)
     1964          {
     1965            *fLog << err << "Unable to write MCalibrationIntensityChargeCam to " << oname << endl;
     1966            return kFALSE;
     1967          }
     1968        *fLog << inf << "ok." << endl;
     1969
     1970        *fLog << inf << " - MCalibrationIntensityQECam..." << flush;
     1971        if (fIntensQECam.Write()<=0)
     1972          {
     1973            *fLog << err << "Unable to write MCalibrationIntensityQECam to " << oname << endl;
     1974            return kFALSE;
     1975          }
     1976        *fLog << inf << "ok." << endl;
     1977      }
     1978    else
     1979      {
     1980        *fLog << inf << " - MCalibrationChargeCam..." << flush;
     1981        if (fCalibrationCam.Write()<=0)
     1982          {
     1983            *fLog << err << "Unable to write MCalibrationChargeCam to " << oname << endl;
     1984            return kFALSE;
     1985          }
     1986        *fLog << inf << "ok." << endl;
     1987
     1988        *fLog << inf << " - MCalibrationQECam..." << flush;
     1989        if (fQECam.Write()<=0)
     1990          {
     1991            *fLog << err << "Unable to write MCalibrationQECam to " << oname << endl;
     1992            return kFALSE;
     1993          }
     1994        *fLog << inf << "ok." << endl;
     1995      }
     1996   
    19061997    *fLog << inf << " - MCalibrationChargeBlindCam..." << flush;
    19071998    if (fCalibrationBlindCam.Write()<=0)
    1908     {
     1999      {
    19092000        *fLog << err << "Unable to write MCalibrationChargeBlindCam to " << oname << endl;
    19102001        return kFALSE;
    1911     }
     2002      }
    19122003    *fLog << inf << "ok." << endl;
    1913 
     2004   
    19142005    *fLog << inf << " - MCalibrationChargePINDiode..." << flush;
    19152006    if (fCalibrationPINDiode.Write()<=0)
    1916     {
     2007      {
    19172008        *fLog << err << "Unable to write MCalibrationChargePINDiode to " << oname << endl;
    1918         return kFALSE;
    1919     }
    1920     *fLog << inf << "ok." << endl;
    1921 
    1922     *fLog << inf << " - MCalibrationQECam..." << flush;
    1923     if (fQECam.Write()<=0)
    1924     {
    1925         *fLog << err << "Unable to write MCalibrationQECam to " << oname << endl;
    19262009        return kFALSE;
    19272010    }
    … …  
    19482031    }
    19492032
    1950     if (fRelTimes)
    1951     {
    1952         *fLog << inf << " - MCalibrationRelTimeCam..." << flush;
    1953         if (fRelTimeCam.Write()<=0)
    1954         {
    1955             *fLog << err << "Unable to write MCalibrationRelTimeCam to " << oname << endl;
    1956             return kFALSE;
    1957         }
    1958         *fLog << inf << "ok." << endl;
     2033    if (IsRelTimes())
     2034    {
     2035      if (IsIntensity())
     2036        {
     2037          *fLog << inf << " - MCalibrationIntensityRelTimeCam..." << flush;
     2038          if (fIntensRelTimeCam.Write()<=0)
     2039            {
     2040              *fLog << err << "Unable to write MCalibrationIntensityRelTimeCam to " << oname << endl;
     2041              return kFALSE;
     2042            }
     2043          *fLog << inf << "ok." << endl;
     2044        }
     2045      else
     2046        {
     2047          *fLog << inf << " - MCalibrationRelTimeCam..." << flush;
     2048          if (fRelTimeCam.Write()<=0)
     2049            {
     2050              *fLog << err << "Unable to write MCalibrationRelTimeCam to " << oname << endl;
     2051              return kFALSE;
     2052            }
     2053          *fLog << inf << "ok." << endl;
     2054        }
     2055     
    19592056    }
    19602057
    … …  
    19722069  outer[0] = 1;
    19732070         
    1974   TArrayI s0(6);
    1975   s0[0] = 6;
    1976   s0[1] = 1;
    1977   s0[2] = 2;
    1978   s0[3] = 3;
    1979   s0[4] = 4;
    1980   s0[5] = 5;
    1981  
    19822071  TArrayI s1(3);
    19832072  s1[0] = 6;
    … …  
    19942083 
    19952084  TH1D *inout[2];
    1996   inout[0] = cam->ProjectionS(s0, inner, "Inner");
    1997   inout[1] = cam->ProjectionS(s0, outer, "Outer");
    1998   FixDataCheckHist(inout[0]);
    1999   FixDataCheckHist(inout[1]);
    2000  
    2001   inout[0]->SetTitle(Form("%s %s",cam->GetTitle(),"Inner"));
    2002   inout[1]->SetTitle(Form("%s %s",cam->GetTitle(),"Outer"));
    2003  
     2085
    20042086  for (int i=0; i<2; i++)
    20052087    {
    … …  
    20082090      gPad->SetTicks();
    20092091     
     2092      inout[i] = cam->ProjectionS(TArrayI(), TArrayI(1,&i), i==0 ? "Inner" : "Outer");
     2093      FixDataCheckHist(inout[i]);
     2094      inout[i]->SetTitle(Form("%s %s",cam->GetTitle(),i==0 ? "Inner" : "Outer"));
    20102095      inout[i]->SetDirectory(NULL);
    20112096      inout[i]->SetLineColor(kRed+i);
    20122097      inout[i]->SetBit(kCanDelete);
    20132098      inout[i]->Draw();
    2014       inout[i]->Fit("gaus", "Q");
     2099      //
     2100      // Display the outliers as dead and noisy pixels
     2101      //
     2102      if (!inout[i]->Fit("gaus","0Q"))
     2103        DisplayOutliers(inout[i],whatsmall,whatbig);
     2104
    20152105      gPad->Modified();
    20162106      gPad->Update();     
    … …  
    20242114
    20252115      TLegend *leg2 = new TLegend(0.55,0.4,0.99,0.6);
    2026       //
    2027       // Display the outliers as dead and noisy pixels
    2028       //
    2029       DisplayOutliers(inout[i],whatsmall,whatbig);
    20302116     
    20312117      //
    … …  
    20332119      //
    20342120      TH1D *half[2];
    2035       half[0] = cam->ProjectionS(s1, i==0 ? inner : outer , "Sector 6-1-2");
    2036       half[1] = cam->ProjectionS(s2, i==0 ? inner : outer , "Sector 3-4-5");
     2121      half[0] = cam->ProjectionS(s1, TArrayI(1,&i), "Sector 6-1-2");
     2122      half[1] = cam->ProjectionS(s2, TArrayI(1,&i), "Sector 3-4-5");
    20372123     
    20382124      for (int j=0; j<2; j++)
Note: See TracChangeset for help on using the changeset viewer.