Ignore:
Timestamp:
05/27/05 16:33:19 (20 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mcalib
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/mcalib/MCalibrateData.cc

    r7053 r7099  
    135135using namespace std;
    136136
    137 const Float_t MCalibrateData::fgCalibConvMinLimit = 0.01;
    138 const Float_t MCalibrateData::fgCalibConvMaxLimit = 5.;
     137const Float_t MCalibrateData::gkCalibConvMinLimit = 0.01;
     138const Float_t MCalibrateData::gkCalibConvMaxLimit = 5.;
     139
     140const MCalibrateData::CalibrationMode_t MCalibrateData::gkDefault = kFfactor;
     141
    139142// --------------------------------------------------------------------------
    140143//
    … …  
    441444{
    442445
    443     *fLog << inf << GetDescriptor()
    444           << ": Updating Conversion Factors... " << endl;
     446    *fLog << inf << GetDescriptor() << ": Updating Conversion Factors... " << endl;
    445447   
    446448    fCalibConsts.Reset();
    … …  
    495497            if (updatecam)
    496498              {
    497                 MCalibrationChargePix &upix = (MCalibrationChargePix&)(*updatecam)[pixidx];
     499                const MCalibrationChargePix &upix = (MCalibrationChargePix&)(*updatecam)[pixidx];
     500
    498501                //
    499502                // Correct for the possible change in amplification of the individual pixels chain
    500503                //
    501504                const Float_t pixmean = upix.GetConvertedMean();
    502                 calibUpdate = (pixmean == 0.) ? 1. : pix.GetConvertedMean() / pixmean;
     505                calibUpdate = pixmean==0 ? 1 : pix.GetConvertedMean()/pixmean;
     506
    503507                //
    504508                // Correct for global shifts in light emission
    505509                //
    506                 MCalibrationChargePix &ugpix = (MCalibrationChargePix&)updatecam->GetAverageArea(0);
     510                const MCalibrationChargePix &ugpix = (MCalibrationChargePix&)updatecam->GetAverageArea(0);
    507511
    508512                const Float_t globmean = avpix.GetConvertedMean();
    509                 calibUpdate = (globmean == 0.) ? 1. : ugpix.GetConvertedMean() / globmean;
    510 
    511                 MBadPixelsPix         &ubad  = (MBadPixelsPix&)        updatecam->GetAverageBadArea(0);
     513                calibUpdate = globmean==0 ? 1 : ugpix.GetConvertedMean()/globmean;
     514
     515                MBadPixelsPix &ubad  = (MBadPixelsPix&)updatecam->GetAverageBadArea(0);
    512516                if (ubad.IsUncalibrated(MBadPixelsPix::kChargeIsPedestal))
    513517                  {
    … …  
    581585                calibUpdate = 1.;
    582586                break;
     587
     588              default:
     589                  break;
    583590              } /* switch calibration mode */
    584591          } /* if(fCalibrationMode!=kNone) */
    … …  
    755762        out << "   " << GetUniqueName() << ".EnablePedestalType(MCalibrateData::kRun)" << endl;
    756763
    757     if (fCalibrationMode != kDefault)
     764    if (fCalibrationMode != gkDefault)
    758765    {
    759766        out << "   " << GetUniqueName() << ".SetCalibrationMode(MCalibrateData::";
    … …  
    865872    {
    866873      *fLog << all
    867             << Form("%s%3i","Pixel: ",pixidx)
    868             << Form("%s%4.2f","  CalibConst: ",fCalibConsts[pixidx])
    869             << Form("%s%4.2f","  F-Factor: ",fCalibFFactors[pixidx])
    870             << Form("%s%4.2f","  HiLoConv: ",fHiLoConv[pixidx])
     874            << "Pixel: " << Form("%3i",pixidx)
     875            << "  CalibConst: " << Form("%4.2f",fCalibConsts[pixidx])
     876            << "  F-Factor: "   << Form("%4.2f",fCalibFFactors[pixidx])
     877            << "  HiLoConv: "   << Form("%4.2f",fHiLoConv[pixidx])
    871878            << endl;
    872879    }

  • trunk/MagicSoft/Mars/mcalib/MCalibrateData.h

    r6855 r7099  
    2323class MCalibrateData : public MTask
    2424{
     25public:
     26    enum CalibrationMode_t
     27    {
     28        kSkip       =  0,
     29        kNone       =  1,
     30        kFlatCharge =  2,
     31        kBlindPixel =  3,
     32        kFfactor    =  4,
     33        kPinDiode   =  5,
     34        kCombined   =  6,
     35        kDummy      =  7
     36    };
     37
     38    enum PedestalType_t
     39    {
     40        kNo    = BIT(0),
     41        kRun   = BIT(1),
     42        kEvent = BIT(2)
     43    };
     44
     45    enum SignalType_t
     46    {
     47        kPhe,
     48        kPhot
     49    };
     50
     51    static const CalibrationMode_t gkDefault;
     52
    2553private:
    26 
    27   static const Float_t  fgCalibConvMinLimit; //! Minimum limit for conv. factor
    28   static const Float_t  fgCalibConvMaxLimit; //! Maximum limit for conv. factor
     54  static const Float_t gkCalibConvMinLimit; //! Minimum limit for conv. factor
     55  static const Float_t gkCalibConvMaxLimit; //! Maximum limit for conv. factor
    2956
    3057  Float_t  fCalibConvMinLimit;           // Minimum limit for conv. factor
    … …  
    3966  MCalibConstCam        *fCalibConstCam; //! Temporary calib consts storage
    4067 
    41   UShort_t fCalibrationMode;             // Flag defining the calibration mode (CalibrationMode_t)
     68  CalibrationMode_t fCalibrationMode;    // Flag defining the calibration mode (CalibrationMode_t)
    4269  Byte_t   fPedestalFlag;                // Flags defining to calibrate the pedestal each event or each run
    4370  Byte_t   fSignalType;                  // Flag defining the signal type (kphot or kphe)
    … …  
    6390 
    6491public:
    65 
    66     enum CalibrationMode_t
    67     {
    68         kSkip       =  0,
    69         kNone       =  1,
    70         kFlatCharge =  2,
    71         kBlindPixel =  3,
    72         kFfactor    =  4,
    73         kPinDiode   =  5,
    74         kCombined   =  6,
    75         kDummy      =  7
    76     };
    77 
    78   static const CalibrationMode_t kDefault = kFfactor;
    79  
    80   enum PedestalType_t
    81   {
    82       kNo    = BIT(0),
    83       kRun   = BIT(1),
    84       kEvent = BIT(2)
    85   };
    86  
    87   enum SignalType_t
    88   {
    89       kPhe,
    90       kPhot
    91   };
    92  
    93   MCalibrateData(CalibrationMode_t calmode=kDefault,
    94                  const char *name=NULL, const char *title=NULL);
     92    MCalibrateData(CalibrationMode_t calmode=gkDefault,
     93                   const char *name=NULL, const char *title=NULL);
    9594 
    9695  void   AddPedestal(const char *name="Cam");
    … …  
    106105  Bool_t TestPedestalFlag(PedestalType_t i) const { return fPedestalFlag&i ? kTRUE : kFALSE; }
    107106 
    108   void   SetCalibrationMode ( CalibrationMode_t calmode=kDefault ) { fCalibrationMode=calmode; }
     107  void   SetCalibrationMode ( CalibrationMode_t calmode=gkDefault) { fCalibrationMode=calmode; }
    109108  void   SetSignalType      ( SignalType_t      sigtype=kPhe    ) { fSignalType     =sigtype; } 
    110109
    111   void   SetCalibConvMinLimit( const Float_t f=fgCalibConvMinLimit ) { fCalibConvMinLimit = f; }
    112   void   SetCalibConvMaxLimit( const Float_t f=fgCalibConvMaxLimit ) { fCalibConvMaxLimit = f; }
     110  void   SetCalibConvMinLimit( const Float_t f=gkCalibConvMinLimit ) { fCalibConvMinLimit = f; }
     111  void   SetCalibConvMaxLimit( const Float_t f=gkCalibConvMaxLimit ) { fCalibConvMaxLimit = f; }
    113112   
    114113  Bool_t UpdateConversionFactors( const MCalibrationChargeCam *updatecam=NULL);
    … …  
    118117
    119118#endif   /* MCalibrateData */
    120 
    121 
    122 
    123 
    124 
    125 

  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

    r7069 r7099  
    181181//  ClassVersion 2:
    182182//   - Float_t fPheErrLimit;
    183 //   + Float_t fPheErrLowerLimit;                   // Lower limit acceptance nr. phe's w.r.t. area idx mean (in sigmas)
    184 //   + Float_t fPheErrUpperLimit;                   // Upper limit acceptance nr. phe's w.r.t. area idx mean (in sigmas)
     183//   + Float_t fPheErrLowerLimit; // Lower limit acceptance nr. phe's w.r.t. area idx mean (in sigmas)
     184//   + Float_t fPheErrUpperLimit; // Upper limit acceptance nr. phe's w.r.t. area idx mean (in sigmas)
     185//
     186//  ClassVersion 3:
     187//   + Bool_t  fUseExtractorRes;  // Include extractor resolution in F-Factor method
    185188//
    186189//////////////////////////////////////////////////////////////////////////////
    … …  
    280283//
    281284MCalibrationChargeCalc::MCalibrationChargeCalc(const char *name, const char *title)
    282     : fGeom(NULL), fSignal(NULL), fCalibPattern(NULL)
     285    : fUseExtractorRes(kFALSE),
     286    fGeom(NULL), fSignal(NULL), fCalibPattern(NULL), fExtractor(NULL)
    283287{
    284288       
    … …  
    661665      return kTRUE;
    662666  }
    663  
     667
    664668  if (fNumProcessed == 0)
    665669    return kTRUE;
    … …  
    10171021            << endl;
    10181022      bad.SetUncalibrated(MBadPixelsPix::kDeviatingNumPhes);
     1023      bad.SetUnsuitable(MBadPixelsPix::kUnsuitableRun);
     1024      return kFALSE;
     1025    }
     1026
     1027  if (cal.GetPheFFactorMethod() < 0.)
     1028    {
     1029      bad.SetUncalibrated(MBadPixelsPix::kDeviatingNumPhes);
     1030      bad.SetUnsuitable(MBadPixelsPix::kUnsuitableRun);
     1031      cal.SetFFactorMethodValid(kFALSE);
     1032      return kFALSE;
     1033    }
     1034
     1035  if (!cal.CalcConvFFactor())
     1036    {
     1037      *fLog << warn << "Could not calculate the Conv. FADC counts to Phes in ";
     1038      *fLog << what << ":           " << Form("%4i", cal.GetPixId()) << endl;
     1039      bad.SetUncalibrated(MBadPixelsPix::kDeviatingNumPhes);
     1040      return kFALSE;
     1041    }
     1042
     1043  if (!fUseExtractorRes)
     1044    return kTRUE;
     1045
     1046  if (!fExtractor)
     1047    {
     1048      *fLog << err << "Extractor resolution has been chosen, but not extractor is set. Cannot calibrate" << endl;
     1049      return kFALSE;
     1050    }
     1051
     1052  const Float_t resinphes = cal.IsHiGainSaturation()
     1053    ? cal.GetPheFFactorMethod()*fExtractor->GetResolutionPerPheLoGain()
     1054    : cal.GetPheFFactorMethod()*fExtractor->GetResolutionPerPheHiGain();
     1055
     1056  Float_t resinfadc = cal.IsHiGainSaturation()
     1057    ? resinphes/cal.GetMeanConvFADC2Phe()/cal.GetConversionHiLo()
     1058    : resinphes/cal.GetMeanConvFADC2Phe();
     1059
     1060  if (resinfadc > 1.5*cal.GetPedRms() )
     1061    {
     1062      *fLog << warn << " Extractor Resolution: " << resinfadc << " bigger than Pedestal RMS " << cal.GetPedRms() << endl;
     1063      resinfadc = cal.GetPedRms();
     1064    }
     1065
     1066  if (!cal.CalcReducedSigma(resinfadc))
     1067    {
     1068        *fLog << warn << "Could not calculate the reduced sigma in " << what;
     1069        *fLog << ":                        " << Form("%4i",cal.GetPixId()) << endl;
     1070        bad.SetUncalibrated( MBadPixelsPix::kChargeSigmaNotValid );
     1071        return kFALSE;
     1072    }
     1073
     1074  if (!cal.CalcFFactor())
     1075    {
     1076        *fLog << warn << "Could not calculate the F-Factor in " << what;
     1077        *fLog << ":                             " << Form("%4i",cal.GetPixId()) << endl;
     1078      bad.SetUncalibrated(MBadPixelsPix::kDeviatingNumPhes);
     1079      bad.SetUnsuitable(MBadPixelsPix::kUnsuitableRun);
    10191080      return kFALSE;
    10201081    }
    … …  
    22992360      rc = kTRUE;
    23002361    }
     2362  if (IsEnvDefined(env, prefix, "UseExtractorRes", print))
     2363    {
     2364      SetUseExtractorRes(GetEnvValue(env, prefix, "UseExtractorRes", fUseExtractorRes));
     2365      rc = kTRUE;
     2366    }
    23012367
    23022368  return rc;

  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.h

    r7028 r7099  
    22#define MARS_MCalibrationChargeCalc
    33
    4 /////////////////////////////////////////////////////////////////////////////
    5 //                                                                         //
    6 // MCalibrationChargeCalc                                                   //
    7 //                                                                         //
    8 // Integrates the time slices of the all pixels of a calibration event     //
    9 // and substract the pedestal value                                        //
    10 //                                                                         //
    11 /////////////////////////////////////////////////////////////////////////////
    12 
    134#ifndef MARS_MTask
    145#include "MTask.h"
    156#endif
    167
     8#ifndef ROOT_TArrayC
     9#include <TArrayC.h>
     10#endif
     11
    1712#ifndef MARS_MBadPixelsPix
    1813#include "MBadPixelsPix.h"
    … …  
    2116#ifndef MARS_MCalibrationCam
    2217#include "MCalibrationCam.h"
    23 #endif
    24 
    25 #ifndef ROOT_TArrayC
    26 #include <TArrayC.h>
    2718#endif
    2819
    … …  
    4435class MBadPixelsIntensityCam;
    4536class MBadPixelsCam;
     37class MExtractor;
    4638
    4739class MCalibrationChargeCalc : public MTask
    … …  
    9183
    9284  Int_t   fNumProcessed;                       // Number of processed events (for Intensity calibration)
     85  Bool_t  fUseExtractorRes;                    // Include extractor resolution in F-Factor method
    9386 
    9487  // Pointers
    … …  
    108101  MCalibrationPattern            *fCalibPattern;   //! Calibration DM pattern
    109102  MPedestalCam                   *fPedestals;      //! Pedestals all pixels (calculated previously from ped.file)
     103  MExtractor                     *fExtractor;      //! Signal Extractor
    110104
    111105  // enums
    … …  
    211205      : CLRBIT(fFlags, kDebug); }
    212206
     207  void SetExtractor(MExtractor   *ext) { fExtractor=ext; }
    213208  void SetPedestals(MPedestalCam *cam) { fPedestals=cam; }
    214209                                                                         
    … …  
    223218  void SetPheErrUpperLimit  ( const Float_t f=fgPheErrUpperLimit       ) { fPheErrUpperLimit  = f;     }
    224219  void SetPulserColor       ( const MCalibrationCam::PulserColor_t col ) { fPulserColor       = col;   }
    225 
    226   ClassDef(MCalibrationChargeCalc, 2)   // Task calculating Calibration Containers and Quantum Efficiencies
     220  void SetUseExtractorRes(Bool_t b=kTRUE) { fUseExtractorRes = b; }
     221
     222  ClassDef(MCalibrationChargeCalc, 3)   // Task calculating Calibration Containers and Quantum Efficiencies
    227223};
    228224

  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

    r7013 r7099  
    674674// in GetRSigma() and GetRSigmaErr()
    675675//
    676 Bool_t MCalibrationChargePix::CalcReducedSigma()
     676Bool_t MCalibrationChargePix::CalcReducedSigma(const Float_t extractorres)
    677677{
    678678
    … …  
    687687  const Float_t pedRmsSquare    = IsHiGainSaturation() ? fLoGainPedRmsSquare    : fPedRms*fPedRms;
    688688  const Float_t pedRmsSquareVar = IsHiGainSaturation() ? fLoGainPedRmsSquareVar : 0.25*fPedVar*pedRmsSquare;
     689  const Float_t extractorresVar = extractorres * extractorres;
    689690
    690691  if (IsDebug())
    … …  
    696697            << " PedRmsSquare: " << pedRmsSquare
    697698            << " pedRmsSquareVar: " << pedRmsSquareVar
     699            << " extractorresVar: " << extractorresVar
    698700            << endl;
    699701    }
    … …  
    705707  // Calculate the reduced sigmas
    706708  //
    707   fRSigmaSquare = sigmaSquare - pedRmsSquare;
     709  fRSigmaSquare = sigmaSquare - pedRmsSquare - extractorresVar;
    708710
    709711  if (IsDebug())
    … …  
    764766Bool_t MCalibrationChargePix::CalcFFactor()
    765767{
    766 
    767   if (fRSigmaSquare < 0.)
    768     return kFALSE;
     768    SetFFactorMethodValid(kFALSE);
     769
     770    if (fRSigmaSquare < 0.)
     771        return kFALSE;
    769772 
    770773  //

  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h

    r7013 r7099  
    5353
    5454  MCalibrationChargePix(const char *name=NULL, const char *title=NULL);
    55   ~MCalibrationChargePix() {}
    5655 
    5756  void Clear(Option_t *o="");
    … …  
    111110
    112111  // Calculations
    113   void   CalcLoGainPedestal       ( const Float_t logainsamples);
    114   Bool_t CalcReducedSigma  (  );
    115   Bool_t CalcFFactor       (  );
    116   Bool_t CalcConvFFactor   (  ); 
    117   Bool_t CalcMeanFFactor   ( const Float_t nphotons, const Float_t nphotonsrelvar );
     112  void   CalcLoGainPedestal(const Float_t logainsamples);
     113  Bool_t CalcReducedSigma(const Float_t extractorres=0);
     114  Bool_t CalcFFactor();
     115  Bool_t CalcConvFFactor();
     116  Bool_t CalcMeanFFactor(const Float_t nphotons, const Float_t nphotonsrelvar);
    118117 
    119118  ClassDef(MCalibrationChargePix, 3)    // Container Charge Calibration Results Pixel
Note: See TracChangeset for help on using the changeset viewer.