Ignore:
Timestamp:
01/27/05 20:44:17 (20 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mjobs
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r6031 r6065  
    101101#include <TPaveText.h>
    102102#include <TPaveStats.h>
     103#include <TEnv.h>
    103104
    104105#include "MLog.h"
    … …  
    175176const Int_t MJCalibration::gkSpecialPixelsContInstallation = 34057;
    176177const Int_t MJCalibration::gkThirdBlindPixelInstallation   = 43308;
    177 
    178 const Double_t MJCalibration::fgConvFADC2PheMin           = 0.;
    179 const Double_t MJCalibration::fgConvFADC2PheMax           = 1.5;
    180 const Double_t MJCalibration::fgConvFADC2PhotMin          = 0.;
    181 const Double_t MJCalibration::fgConvFADC2PhotMax          = 10.;
    182 const Double_t MJCalibration::fgQEMin                     = 0.;
    183 const Double_t MJCalibration::fgQEMax                     = 0.3;
    184 const Double_t MJCalibration::fgArrivalTimeMin            = 1.;
    185 const Double_t MJCalibration::fgArrivalTimeMax            = 10.;
    186 
    187 const Float_t  MJCalibration::fgRefConvFADC2PheInner      = 0.14;       
    188 const Float_t  MJCalibration::fgRefConvFADC2PheOuter      = 0.4;       
    189 const Float_t  MJCalibration::fgRefConvFADC2PhotInner     = 0.8;
    190 const Float_t  MJCalibration::fgRefConvFADC2PhotOuter     = 3.8;
    191 const Float_t  MJCalibration::fgRefQEInner                = 0.18;           
    192 const Float_t  MJCalibration::fgRefQEOuter                = 0.12;
    193 const Float_t  MJCalibration::fgRefArrivalTimeInner       = 4.5;           
    194 const Float_t  MJCalibration::fgRefArrivalTimeOuter       = 5.0;
    195 const Float_t  MJCalibration::fgRefArrivalTimeRmsInner    = 0.5;           
    196 const Float_t  MJCalibration::fgRefArrivalTimeRmsOuter    = 0.5;
     178const TString MJCalibration::fgReferenceFile = "mjobs/calibrationref.rc";
    197179// --------------------------------------------------------------------------
    198180//
    … …  
    224206  SetNoStorage(kFALSE);
    225207  SetHistsStorage(kFALSE);
     208
     209  SetReferenceFile();
     210
     211  fConvFADC2PheMin           = 0.;
     212  fConvFADC2PheMax           = 1.5;
     213  fConvFADC2PhotMin          = 0.;
     214  fConvFADC2PhotMax          = 10.;
     215  fQEMin                     = 0.;
     216  fQEMax                     = 0.3;
     217  fArrivalTimeMin            = 1.;
     218  fArrivalTimeMax            = 10.;
     219 
     220  fRefConvFADC2PheInner      = 0.14;       
     221  fRefConvFADC2PheOuter      = 0.4;       
     222  fRefConvFADC2PhotInner     = 0.8;
     223  fRefConvFADC2PhotOuter     = 3.8;
     224  fRefQEInner                = 0.18;           
     225  fRefQEOuter                = 0.12;
     226  fRefArrivalTimeInner       = 4.5;           
     227  fRefArrivalTimeOuter       = 5.0;
     228  fRefArrivalTimeRmsInner    = 0.5;           
     229  fRefArrivalTimeRmsOuter    = 0.5;
     230
    226231}
    227232
    … …  
    526531      // for the datacheck, fix the ranges!!
    527532      //
    528       // obj1->SetMinimum(fgChargeMin);
    529       // obj1->SetMaximum(fgChargeMax);
     533      // obj1->SetMinimum(fChargeMin);
     534      // obj1->SetMaximum(fChargeMax);
    530535      //
    531536      // Set the datacheck sizes:
    … …  
    591596      //
    592597      //      obj3->SetMinimum(0.);
    593       // obj3->SetMaximum(fgChargeMax);
     598      // obj3->SetMaximum(fChargeMax);
    594599      //
    595600      // Set the datacheck sizes:
    … …  
    623628      // for the datacheck, fix the ranges!!
    624629      //
    625       // obj3->SetMinimum(fgChargeMin);
    626       // obj3->SetMaximum(fgChargeMax);
     630      // obj3->SetMinimum(fChargeMin);
     631      // obj3->SetMaximum(fChargeMax);
    627632      //
    628633      // Set the datacheck sizes:
    … …  
    658663      // for the datacheck, fix the ranges!!
    659664      //
    660       obj5->SetMinimum(fgConvFADC2PheMin);
    661       obj5->SetMaximum(fgConvFADC2PheMax);
     665      obj5->SetMinimum(fConvFADC2PheMin);
     666      obj5->SetMaximum(fConvFADC2PheMax);
    662667      //
    663668      // Set the datacheck sizes:
    … …  
    690695      // for the datacheck, fix the ranges!!
    691696      //
    692       obj6->SetMinimum(fgQEMin);
    693       obj6->SetMaximum(fgQEMax);
     697      obj6->SetMinimum(fQEMin);
     698      obj6->SetMaximum(fQEMax);
    694699      //
    695700      // Set the datacheck sizes:
    … …  
    723728      // for the datacheck, fix the ranges!!
    724729      //
    725       obj7->SetMinimum(fgConvFADC2PhotMin);
    726       obj7->SetMaximum(fgConvFADC2PhotMax);
     730      obj7->SetMinimum(fConvFADC2PhotMin);
     731      obj7->SetMaximum(fConvFADC2PhotMax);
    727732      //
    728733      // Set the datacheck sizes:
    … …  
    757762      // for the datacheck, fix the ranges!!
    758763      //
    759       obj10->SetMinimum(fgArrivalTimeMin);
    760       obj10->SetMaximum(fgArrivalTimeMax);
     764      obj10->SetMinimum(fArrivalTimeMin);
     765      obj10->SetMaximum(fArrivalTimeMax);
    761766      //
    762767      // Set the datacheck sizes:
    … …  
    789794      // for the datacheck, fix the ranges!!
    790795      //
    791       //      obj11->SetMinimum(fgArrivalTimeMin);
    792       //      obj11->SetMaximum(fgArrivalTimeMax);
     796      //      obj11->SetMinimum(fArrivalTimeMin);
     797      //      obj11->SetMaximum(fArrivalTimeMax);
    793798      //
    794799      // Set the datacheck sizes:
    … …  
    971976          // for the datacheck, fix the ranges!!
    972977          //
    973           //   obj10->SetMinimum(fgConvFADC2PheMin);
    974           //   obj10->SetMaximum(fgConvFADC2PheMax);
     978          //   obj10->SetMinimum(fConvFADC2PheMin);
     979          //   obj10->SetMaximum(fConvFADC2PheMax);
    975980          //
    976981          // Set the datacheck sizes:
    … …  
    10031008          // for the datacheck, fix the ranges!!
    10041009          //
    1005           //          obj11->SetMinimum(fgQEMin);
    1006           //          obj11->SetMaximum(fgQEMax);
     1010          //          obj11->SetMinimum(fQEMin);
     1011          //          obj11->SetMaximum(fQEMax);
    10071012          //
    10081013          // Set the datacheck sizes:
    … …  
    12311236
    12321237  if (what == 0)
    1233     l1 = line.DrawLine(0, fgRefQEInner, x, fgRefQEInner);     
     1238    l1 = line.DrawLine(0, fRefQEInner, x, fRefQEInner);     
    12341239  else if (what == 1)
    1235     l1 = line.DrawLine(0, fgRefConvFADC2PhotInner, x, fgRefConvFADC2PhotInner);     
     1240    l1 = line.DrawLine(0, fRefConvFADC2PhotInner, x, fRefConvFADC2PhotInner);     
    12361241  else if (what == 2)
    1237     l1 = line.DrawLine(0, fgRefConvFADC2PheInner, x, fgRefConvFADC2PheInner );
     1242    l1 = line.DrawLine(0, fRefConvFADC2PheInner, x, fRefConvFADC2PheInner );
    12381243  else if (what == 3)
    1239     l1 = line.DrawLine(0, fgRefArrivalTimeInner, x, fgRefArrivalTimeInner );
     1244    l1 = line.DrawLine(0, fRefArrivalTimeInner, x, fRefArrivalTimeInner );
    12401245  else if (what == 4)
    1241     l1 = line.DrawLine(0, fgRefArrivalTimeRmsInner, x, fgRefArrivalTimeRmsInner );
     1246    l1 = line.DrawLine(0, fRefArrivalTimeRmsInner, x, fRefArrivalTimeRmsInner );
    12421247 
    12431248  if (geom->InheritsFrom("MGeomCamMagic"))
    … …  
    12481253        {
    12491254        case 0:
    1250           line.DrawLine(x2, fgRefQEOuter, 398, fgRefQEOuter);     
     1255          line.DrawLine(x2, fRefQEOuter, 398, fRefQEOuter);     
    12511256          break;
    12521257        case 1:
    1253           line.DrawLine(x2, fgRefConvFADC2PhotOuter, 398, fgRefConvFADC2PhotOuter );     
     1258          line.DrawLine(x2, fRefConvFADC2PhotOuter, 398, fRefConvFADC2PhotOuter );     
    12541259          break;
    12551260        case 2:
    1256           line.DrawLine(x2, fgRefConvFADC2PheOuter, 398, fgRefConvFADC2PheOuter);
     1261          line.DrawLine(x2, fRefConvFADC2PheOuter, 398, fRefConvFADC2PheOuter);
    12571262          break;
    12581263        case 3:
    1259           line.DrawLine(x2, fgRefArrivalTimeOuter, 398, fgRefArrivalTimeOuter);
     1264          line.DrawLine(x2, fRefArrivalTimeOuter, 398, fRefArrivalTimeOuter);
    12601265          break;
    12611266        case 4:
    1262           line.DrawLine(x2, fgRefArrivalTimeRmsOuter, 398, fgRefArrivalTimeRmsOuter);
     1267          line.DrawLine(x2, fRefArrivalTimeRmsOuter, 398, fRefArrivalTimeRmsOuter);
    12631268          break;
    12641269        }
    … …  
    13821387    SetGeometry(GetEnv("Geometry", fGeometry));
    13831388
    1384     return kTRUE;
     1389    fReferenceFile = GetEnv("ReferenceFile",fReferenceFile.Data());
     1390
     1391    TEnv refenv(fReferenceFile);
     1392   
     1393    fConvFADC2PheMin           = refenv.GetValue("ConvFADC2PheMin",fConvFADC2PheMin);
     1394    fConvFADC2PheMax           = refenv.GetValue("ConvFADC2PheMax",fConvFADC2PheMax);
     1395    fConvFADC2PhotMin          = refenv.GetValue("ConvFADC2PhotMin",fConvFADC2PhotMin);
     1396    fConvFADC2PhotMax          = refenv.GetValue("ConvFADC2PhotMax",fConvFADC2PhotMax);
     1397    fQEMin                     = refenv.GetValue("QEMin",fQEMin);
     1398    fQEMax                     = refenv.GetValue("QEMax",fQEMax);
     1399    fArrivalTimeMin            = refenv.GetValue("ArrivalTimeMin",fArrivalTimeMin);
     1400    fArrivalTimeMax            = refenv.GetValue("ArrivalTimeMax",fArrivalTimeMax);
     1401                                                                       
     1402    fRefConvFADC2PheInner      = refenv.GetValue("RefConvFADC2PheInner",fRefConvFADC2PheInner);
     1403    fRefConvFADC2PheOuter      = refenv.GetValue("RefConvFADC2PheOuter",fRefConvFADC2PheOuter);
     1404    fRefConvFADC2PhotInner     = refenv.GetValue("RefConvFADC2PhotInner",fRefConvFADC2PhotInner);
     1405    fRefConvFADC2PhotOuter     = refenv.GetValue("RefConvFADC2PhotOuter",fRefConvFADC2PhotOuter);
     1406    fRefQEInner                = refenv.GetValue("RefQEInner",fRefQEInner);
     1407    fRefQEOuter                = refenv.GetValue("RefQEOuter",fRefQEOuter);
     1408    fRefArrivalTimeInner       = refenv.GetValue("RefArrivalTimeInner",fRefArrivalTimeInner);
     1409    fRefArrivalTimeOuter       = refenv.GetValue("RefArrivalTimeOuter",fRefArrivalTimeOuter);
     1410    fRefArrivalTimeRmsInner    = refenv.GetValue("RefArrivalTimeRmsInner",fRefArrivalTimeRmsInner);
     1411    fRefArrivalTimeRmsOuter    = refenv.GetValue("RefArrivalTimeRmsOuter",fRefArrivalTimeRmsOuter);
     1412                                 
     1413    return kTRUE;               
    13851414}
    13861415

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

    r5561 r6065  
    6464  static const Int_t gkThirdBlindPixelInstallation;     //! Run number upon which third blind pixel was installed 
    6565
    66   static const Double_t fgConvFADC2PheMin;              //! Histogram minimum for conversion factor to phes
    67   static const Double_t fgConvFADC2PheMax;              //! Histogram maximum for conversion factor to phes
    68   static const Double_t fgConvFADC2PhotMin;             //! Histogram minimum for conversion factor to phs
    69   static const Double_t fgConvFADC2PhotMax;             //! Histogram maixmum for conversion factor to phs
    70   static const Double_t fgQEMin;                        //! Histogram minimum for quantum efficiency
    71   static const Double_t fgQEMax;                        //! Histogram maximum for quantum efficiency
    72   static const Double_t fgArrivalTimeMin;               //! Histogram minimum for arrival times
    73   static const Double_t fgArrivalTimeMax;               //! Histogram maximum for arrival times
    74  
    75   static const Float_t  fgRefConvFADC2PheInner;         //! Reference value for the conversion factor to phes - inner
    76   static const Float_t  fgRefConvFADC2PheOuter;         //! Reference value for the conversion factor to phes - outer
    77   static const Float_t  fgRefConvFADC2PhotInner;        //! Reference value for the conversion factor to phs - inner
    78   static const Float_t  fgRefConvFADC2PhotOuter;        //! Reference value for the conversion factor to phs - outer
    79   static const Float_t  fgRefQEInner;                   //! Reference value for the quantum eff. cascades - inner
    80   static const Float_t  fgRefQEOuter;                   //! Reference value for the quantum eff. cascades - outer
    81   static const Float_t  fgRefArrivalTimeInner;          //! Reference value for the arrival time - inner
    82   static const Float_t  fgRefArrivalTimeOuter;          //! Reference value for the arrival time - outer
    83   static const Float_t  fgRefArrivalTimeRmsInner;       //! Reference value for the arrival time - inner
    84   static const Float_t  fgRefArrivalTimeRmsOuter;       //! Reference value for the arrival time - outer
     66  static const TString  fgReferenceFile;   //! default for fReferenceFile ("mjobs/calibrationref.rc")
     67
     68  Double_t fConvFADC2PheMin;               //! Histogram minimum for conversion factor to phes
     69  Double_t fConvFADC2PheMax;               //! Histogram maximum for conversion factor to phes
     70  Double_t fConvFADC2PhotMin;              //! Histogram minimum for conversion factor to phs
     71  Double_t fConvFADC2PhotMax;              //! Histogram maixmum for conversion factor to phs
     72  Double_t fQEMin;                         //! Histogram minimum for quantum efficiency
     73  Double_t fQEMax;                         //! Histogram maximum for quantum efficiency
     74  Double_t fArrivalTimeMin;                //! Histogram minimum for arrival times
     75  Double_t fArrivalTimeMax;                //! Histogram maximum for arrival times
     76                                           
     77  Float_t  fRefConvFADC2PheInner;          //! Reference value for the conversion factor to phes - inner
     78  Float_t  fRefConvFADC2PheOuter;          //! Reference value for the conversion factor to phes - outer
     79  Float_t  fRefConvFADC2PhotInner;         //! Reference value for the conversion factor to phs - inner
     80  Float_t  fRefConvFADC2PhotOuter;         //! Reference value for the conversion factor to phs - outer
     81  Float_t  fRefQEInner;                    //! Reference value for the quantum eff. cascades - inner
     82  Float_t  fRefQEOuter;                    //! Reference value for the quantum eff. cascades - outer
     83  Float_t  fRefArrivalTimeInner;           //! Reference value for the arrival time - inner
     84  Float_t  fRefArrivalTimeOuter;           //! Reference value for the arrival time - outer
     85  Float_t  fRefArrivalTimeRmsInner;        //! Reference value for the arrival time - inner
     86  Float_t  fRefArrivalTimeRmsOuter;        //! Reference value for the arrival time - outer
     87                                           
     88  TString  fReferenceFile;                 // File name containing the reference values
    8589
    8690  MRunIter       *fRuns;                                // Calibration files
    … …  
    213217  void SetHistsStorage ( const Bool_t b=kTRUE ) { b ? SETBIT(fStorage,kHistsStorage) : CLRBIT(fStorage,kHistsStorage); }
    214218 
     219  void SetReferenceFile( const TString ref=fgReferenceFile ) { fReferenceFile = ref; }
     220
    215221  // Precessing
    216222  Bool_t ReadCalibrationCam();
Note: See TracChangeset for help on using the changeset viewer.