Changeset 3356 for trunk/MagicSoft


Ignore:
Timestamp:
02/28/04 20:43:30 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mcalib
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

    r3288 r3356  
    2525/////////////////////////////////////////////////////////////////////////////
    2626//                                                                         //
    27 // MCalibrationChargePix                                                         //
     27// MCalibrationChargePix                                                   //
    2828//                                                                         //
    29 // This is the storage container to hold informations about the pedestal   //
    30 // (offset) value of one Pixel (PMT).                                      //
     29// Storage container to hold informations about the calibration values     //
     30// values of one Pixel (PMT).                                              //
    3131//                                                                         //
    3232// The following values are initialized to meaningful values:
     
    3434// - The Electronic Rms to 1.5 per FADC slice
    3535// - The uncertainty about the Electronic RMS to 0.3 per slice
    36 // - The F-Factor is assumed to have been measured in Munich to 1.13 - 1.17. 
    37 //   with the Munich definition of the F-Factor, thus: 
     36// - The F-Factor is assumed to have been measured in Munich to 1.13 - 1.17.
     37//   with the Munich definition of the F-Factor, thus:
    3838//   F = Sigma(Out)/Mean(Out) * Mean(In)/Sigma(In)
    3939//   Mean F-Factor  = 1.15
     
    243243}
    244244 
     245
     246void  MCalibrationChargePix::SetMeanCharge( const Float_t f )
     247{
     248    if (IsHiGainSaturation())
     249        fLoGainMeanCharge = f;
     250    else
     251        fHiGainMeanCharge = f;
     252}
     253
     254void  MCalibrationChargePix::SetMeanChargeErr( const Float_t f )
     255{
     256    if (IsHiGainSaturation())
     257        fLoGainMeanChargeErr = f;
     258    else
     259        fHiGainMeanChargeErr = f;
     260
     261}
     262
     263void  MCalibrationChargePix::SetSigmaCharge( const Float_t f ) 
     264{
     265    if (IsHiGainSaturation())
     266        fLoGainSigmaCharge = f;
     267    else
     268        fHiGainSigmaCharge = f;
     269}
     270
     271
     272void  MCalibrationChargePix::SetSigmaChargeErr( const Float_t f )   
     273{
     274    if (IsHiGainSaturation())
     275        fLoGainSigmaChargeErr = f;
     276    else
     277        fHiGainSigmaChargeErr = f;
     278
     279}
     280
     281
    245282
    246283// --------------------------------------------------------------------------
     
    735772  fPheFFactorMethod = ffactorsquare * chargeSquare / rsigmaSquare;
    736773
     774  if (fPheFFactorMethod < 0.)
     775      return kFALSE;
     776
    737777  //
    738778  // Calculate the Error of Nphe
     
    751791  fMeanConversionFFactorMethod         =  fPheFFactorMethod / fAverageQE / GetMeanCharge();
    752792
    753   if (fMeanConversionFFactorMethod < 0.)
    754       return kFALSE;
    755793
    756794  //
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h

    r3282 r3356  
    167167  void SetLoGainSigmaChargeErr  ( const Float_t f ) { fLoGainSigmaChargeErr  = f; }
    168168
     169  void SetMeanCharge            ( const Float_t f );               
     170  void SetMeanChargeErr         ( const Float_t f );             
     171  void SetSigmaCharge           ( const Float_t f );               
     172  void SetSigmaChargeErr        ( const Float_t f );           
     173
    169174  void SetHiGainNumPickup       ( const Int_t   i ) { fHiGainNumPickup       = i; }
    170175  void SetLoGainNumPickup       ( const Int_t   i ) { fLoGainNumPickup       = i; }
  • trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc

    r3355 r3356  
    3131//
    3232// MHCalibrationChargeHiGainPix, MHCalibrationChargeLoGainPix for each pixel
    33 // and additionally two classes containing an average of the inner and the outer
    34 // pixels, respectively
     33// and additionally for an average of the inner and the outer pixels, respectively.
    3534//
    3635// By default, subsequently the hi-gain classes are treated unless
     
    4443// its errors and the fit probability are extracted. If none of these values are
    4544// NaN's and if the probability is bigger than fProbLimit (default: 0.5%), the fit
    46 // is declared valid. Otherwise, histogram means and RMS's are taken directly, but the
    47 // flag kFitValid is set to FALSE. Outliers of more than fPickUpLimit (default: 5) sigmas
     45// is declared valid. Otherwise, the fit is repeated within ranges of the previous mean
     46// +- 5 sigma. In case that this does not help, the histogram means and RMS's are taken directly,
     47// but the flag kFitValid is set to FALSE. Outliers of more than fPickUpLimit (default: 5) sigmas
    4848// from the mean are counted as PickUp events.
    4949//
    5050// Additionally, the slice number with the highest value is stored and a corresponding
    51 // histogram is filled. This histogram serves only for a rought cross-check if the
     51// histogram is filled. This histogram serves only for a rough cross-check if the
    5252// signal does not lie at the edges of chose extraction window.
    5353//
  • trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeHiGainPix.cc

    r3355 r3356  
    8080}
    8181
    82 void MHCalibrationChargeHiGainPix::Draw(const Option_t *opt)
    83 {
    84 
    85   TString option(opt);
    86   option.ToLower();
    87  
    88   Int_t win = 1;
    89 
    90   TVirtualPad *oldpad = gPad ? gPad : MH::MakeDefCanvas(this,600, 600);
    91   TVirtualPad *pad    = NULL;
    92 
    93   oldpad->SetBorderMode(0); 
    94 
    95   pad->SetTicks();
    96   pad->SetBorderMode(0);
    97    
    98 
    99   pad->Divide(1,win);
    100   pad->cd(1);
    101 
    102   if (!IsEmpty())
    103     pad->SetLogy();
    104 
    105   MHGausEvents::Draw(opt);
    106 
    107 }
    10882
    10983
    11084
     85
     86
     87
     88
     89
     90
Note: See TracChangeset for help on using the changeset viewer.