KVINDRACsITotLightReconstructor.h

 
#ifndef __KVINDRACSITOTLIGHTRECONSTRUCTOR_H
#define __KVINDRACSITOTLIGHTRECONSTRUCTOR_H
 
#include "KVDetectorSignal.h"
 
#include <KVINDRADetector.h>
 
#define NOT_CALCULATED 0
#define CALCULATED_WITH_GAIN_CORRECTION 1
#define NO_GAIN_CORRECTION 4
#define NEGATIVE_PEDESTAL_CORRECTED_VALUE 15
#define CALCULATION_NOT_CONVERGED 12
 
class KVINDRACsITotLightReconstructor : public KVDetectorSignal {
   KVDetectorSignal* fFast = nullptr;
   KVDetectorSignal* fSlow = nullptr;
 
    mutable UInt_t fLumTotStatus;
    Double_t tau;
   Double_t c2;
   Double_t p0;
   Double_t p1;
   Double_t eps;
   /* Cette variable n'est pas utilisee, et ne figure pas dans e.g. le fortran de la 4eme campagne
    Float_t pre=0.4318;
   */
   Double_t c1;
   Double_t tau0;
   Double_t tau1;
   Double_t tau2;
 
   void Calculate() const;
 
   UInt_t GetStatusLumiere() const;
   Bool_t LightIsGood() const;
 
public:
   KVINDRACsITotLightReconstructor() = default;
   KVINDRACsITotLightReconstructor(const KVINDRADetector* det)
      : KVDetectorSignal("TotLight", det),
        fFast(det->GetDetectorSignal("R_PedCor")),
        fSlow(det->GetDetectorSignal("L_PedCor")),
        fLumTotStatus(NOT_CALCULATED), tau(20)
   {
      p0 = 400;
      p1 = 900;
      eps = 1.e-4;
      /* Cette variable n'est pas utilisee, et ne figure pas dans e.g. le fortran de la 4eme campagne
       Float_t pre=0.4318;
      */
      c1 = 1.5;
      tau0 = 390.;
      tau1 = 1590.;
      tau2 = 3090.;
      int ring = det->GetRingNumber();
      int module = det->GetModuleNumber();
      /****************************************************************************
      * Definition of PM rise-time constant tau
      * Comment based on a note from Bernard Borderie:
      *   rings 11-16 use PM bases with a different component, thus tau is 60 nsec
      *   ring 16, module 5 is an exception, it has a tau of 20 nsec !!
      *   ring 5, module 11 is also an exception, it has a tau of 60 nsec !!
      *   (all others have tau=20ns)
      */
      if (ring >= 11 && ring <= 16)
         tau = 60.;
      if (ring == 16 && module == 5)
         tau = 20.;
      if (ring == 5 && module == 11)
         tau = 60.;
      /****************************************************************************/
      if (ring > 3) {
         c2 = 2.2;
      }
      else {
         c2 = 3.3;
      }
      if (det) SetTitle(Form("Signal %s calculated from signals %s and %s of detector %s", GetName(), fFast->GetName(), fSlow->GetName(), GetDetector()->GetName()));
   }
 
   Double_t GetValue(const KVNameValueList& = "") const override
   {
      Calculate();
      return KVDetectorSignal::GetValue();
   }
   void Reset() override
   {
      set_value(0);
      fLumTotStatus = NOT_CALCULATED;
   }
   Int_t GetStatus(const TString&) const override;
 
   Bool_t IsRaw() const override
   {
      return kFALSE;
   }
 
   ClassDefOverride(KVINDRACsITotLightReconstructor, 1) //Calculation of light output from fast and slow components for INDRA CsI detectors
};
 
#endif