Changeset 5460


Ignore:
Timestamp:
11/23/04 10:00:51 (20 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/Changelog

    r5459 r5460  
    3838     - stores the histogram RMS in MCalibrationPix in function
    3939       FitHiGainHists and FitLoGainHists
     40
     41   * mjobs/MJCalibration.cc
     42     - display in the datacheck the RMS per Mean instead of the reduced
     43       sigma per mean which is anyhow reflected by the number of photo-
     44       electrons.
    4045
    4146
  • trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc

    r5239 r5460  
    11931193  pix.SetHiGainMean       ( hist.GetMean()      );
    11941194  pix.SetHiGainMeanVar    ( hist.GetMeanErr() * hist.GetMeanErr()   );
     1195  pix.SetHiGainRms        ( hist.GetHistRms()   );
    11951196  pix.SetHiGainSigma      ( hist.GetSigma()     );
    11961197  pix.SetHiGainSigmaVar   ( hist.GetSigmaErr()* hist.GetSigmaErr()  );
     
    12751276  pix.SetLoGainMean       ( hist.GetMean()      );
    12761277  pix.SetLoGainMeanVar    ( hist.GetMeanErr()  * hist.GetMeanErr()   );
     1278  pix.SetLoGainRms        ( hist.GetHistRms()   );
    12771279  pix.SetLoGainSigma      ( hist.GetSigma()     );
    12781280  pix.SetLoGainSigmaVar   ( hist.GetSigmaErr() * hist.GetSigmaErr()  );
  • trunk/MagicSoft/Mars/mhcalib/MHGausEvents.h

    r5135 r5460  
    122122        TH1I    *GetHPowerProbability()        { return fHPowerProbability;  }
    123123  const TH1I    *GetHPowerProbability()  const { return fHPowerProbability;  }
     124  const Double_t GetHistRms()            const { return fHGausHist.GetRMS(); }
    124125  const Double_t GetMean()               const { return fMean;               }
    125126  const Double_t GetMeanErr()            const { return fMeanErr;            }
  • trunk/MagicSoft/Mars/mjobs/MJCalibration.cc

    r5403 r5460  
    132132#include "MCalibrationChargeCalc.h"
    133133#include "MCalibColorSet.h"
     134#include "MCalibrationRelTimeCam.h"
     135#include "MCalibrationRelTimeCalc.h"
    134136
    135137#include "MHGausEvents.h"
     
    139141#include "MHCalibrationChargePINDiode.h"
    140142#include "MHCalibrationRelTimeCam.h"
    141 #include "MCalibrationRelTimeCam.h"
    142 #include "MCalibrationRelTimeCalc.h"
    143143
    144144#include "MReadMarsFile.h"
     145#include "MPedCalcPedRun.h"
    145146#include "MRawFileRead.h"
    146147#include "MGeomApply.h"
     
    341342    MHCamera disp34(geomcam, "TimeOscillating",   "Oscillating Pixels");
    342343    MHCamera disp35(geomcam, "TotalConv",         "Conversion Factor to photons");
     344    MHCamera disp36(geomcam, "RMSperMean",        "Charge histogram RMS per Mean");
    343345   
    344346    MCalibrationChargeCam  *cam    = NULL;
     
    425427    disp29.SetCamContent(*cam, 27);
    426428
     429    // RMS and Mean
     430    disp36.SetCamContent(*cam,32);
     431
    427432    disp1.SetYTitle("Q [FADC counts]");
    428433    disp2.SetYTitle("\\sigma_{Q} [FADC counts]");
     
    463468
    464469    disp35.SetYTitle("Conv.Factor [Ph/FADC counts]");
     470
     471    disp36.SetYTitle("Charge RMS/<Q> [1]");
    465472
    466473    for (UInt_t i=0;i<geomcam.GetNumPixels();i++)
     
    569576     
    570577      //
    571       // REDUCED SIGMAS
     578      // RMS per Charge
    572579      //
    573580     
     
    575582      gPad->SetBorderMode(0);
    576583      gPad->SetTicks();
    577       MHCamera *obj3=(MHCamera*)disp4.DrawCopy("hist");
     584      MHCamera *obj3=(MHCamera*)disp36.DrawCopy("hist");
    578585      //
    579586      // for the datacheck, fix the ranges!!
    580587      //
    581       obj3->SetMinimum(0.);
     588      //      obj3->SetMinimum(0.);
    582589      // obj3->SetMaximum(fgChargeMax);
    583590      //
     
    599606      gPad->SetBorderMode(0);
    600607      if (geomcam.InheritsFrom("MGeomCamMagic"))
    601         DisplayDoubleProject(&disp4, "dead", "noisy");
     608        DisplayDoubleProject(&disp36, "dead", "noisy");
    602609     
    603610      //
     
    15581565    plist.AddToList(&pedcam);
    15591566    plist.AddToList(&fBadPixels);
    1560 
     1567       
    15611568    //
    15621569    // Calibration Results containers
     
    16301637        calcalc.SetOutputPath(fPathOut);
    16311638        calcalc.SetOutputFile(Form("%s-ChargeCalibStat.txt",(const char*)fRuns->GetRunsAsFileName()));
     1639        timecalc.SetOutputPath(fPathOut);
     1640        timecalc.SetOutputFile(Form("%s-ChargeCalibStat.txt",(const char*)fRuns->GetRunsAsFileName()));
    16321641    }
    16331642
     
    16681677    filltme.SetDrawOption(drawoption.Data());
    16691678  */
    1670     //
    1671     // Apply a filter against cosmics
    1672     // (will have to be needed in the future
    1673     // when the calibration hardware-trigger is working)
    1674     //
    1675     MFCosmics cosmics;
    1676     MContinue cont(&cosmics);
    1677 
    16781679    MTaskEnv taskenv("ExtractSignal");
    16791680    taskenv.SetDefault(fExtractor ? fExtractor : &extract2);
     
    16811682    tlist.AddToList(&merge);
    16821683    tlist.AddToList(&apply);
     1684
     1685    MPedCalcPedRun           pedcalc;
     1686    pedcalc.SetWindowSize(fExtractor
     1687                          ? (Int_t)TMath::Nint(fExtractor->GetNumHiGainSamples())
     1688                          : (Int_t)TMath::Nint(extract2.GetNumHiGainSamples())
     1689                          );
     1690    pedcalc.SetRange(fExtractor
     1691                     ? fExtractor->GetHiGainFirst()
     1692                     : extract2.GetHiGainFirst(),
     1693                     fExtractor
     1694                     ? fExtractor->GetHiGainLast()
     1695                     : extract2.GetHiGainLast()
     1696                     );
     1697
     1698    if (IsIntensity())
     1699      tlist.AddToList(&pedcalc);
     1700
     1701    MCalibColorSet colorset;
     1702    tlist.AddToList(&colorset);
     1703
    16831704    tlist.AddToList(&taskenv);
    16841705    tlist.AddToList(&pinext);
     
    16951716   
    16961717    //
    1697     // FIXME: This is not yet implemented in the classes!!!
    1698     //
     1718    // Apply a filter against cosmics
     1719    // (will have to be needed in the future
     1720    // when the calibration hardware-trigger is working)
     1721    //
     1722    MFCosmics cosmics;
     1723    MContinue cont(&cosmics);
     1724
    16991725    if (fColor == MCalibrationCam::kCT1)
    17001726      tlist.AddToList(&cont);
    1701 
    1702     MCalibColorSet colorset;
    1703     tlist.AddToList(&colorset);
    17041727
    17051728    MCalibColorSteer steer;
Note: See TracChangeset for help on using the changeset viewer.