KVINDRAReconNuc.h

/***************************************************************************
                          kvdetectedparticle.h  -  description
                             -------------------
    begin                : Thu Oct 10 2002
    copyright            : (C) 2002 by J.D. Frankland
    email                : frankland@ganil.fr
 
$Id: KVINDRAReconNuc.h,v 1.39 2009/04/03 14:28:37 franklan Exp $
 ***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
 
#ifndef KVINDRADETPART_H
#define KVINDRADETPART_H
 
 
#include "KVReconstructedNucleus.h"
#include "KVINDRADetector.h"
#include "KVTelescope.h"
#include "KVINDRACodes.h"
#include "KVINDRAIDTelescope.h"
 
class KVChIo;
class KVSilicon;
class KVSi75;
class KVSiLi;
class KVCsI;
 
 
 
class KVINDRAReconNuc: public KVReconstructedNucleus {
 
    mutable KVINDRACodes fCodes;
    Bool_t fCoherent;
    Bool_t fPileup;
    Bool_t fUseFullChIoEnergyForCalib;
    Float_t fECsI;
    Float_t fESi;
    Float_t fEChIo;
    Float_t fESi75;
    Float_t fESiLi;
    Bool_t fCorrectCalib;
 
    Bool_t fCoherentSi75SiLiCsI;
    Bool_t fPileupChIo;
    Bool_t fPileupSiLi;
    Bool_t fPileupSi75;
    Bool_t fIncludeEtalonsInCalibration;
 
   void CheckCsIEnergy();
 
   void SetBadCalibrationStatus()
   {
      SetECode(KVINDRACodeMask::kECode15);
      SetEnergy(-1.0);
   }
   void SetNoCalibrationStatus()
   {
      SetECode(KVINDRACodeMask::kECode0);
      SetEnergy(0.0);
   }
   void DoNeutronCalibration();
   void DoBeryllium8Calibration();
   void DoGammaCalibration();
   Bool_t CalculateSiliconDEFromResidualEnergy();
   void CalculateSiLiDEFromResidualEnergy(Double_t ERES);
   void CalculateSi75DEFromResidualEnergy(Double_t ERES);
   void CalculateChIoDEFromResidualEnergy(Double_t ERES);
 
public:
 
   Bool_t AreSiCsICoherent() const
   {
      return fCoherent;
   };
   Bool_t IsSiPileup() const
   {
      return fPileup;
   };
   Bool_t IsSi75Pileup() const
   {
      return fPileupSi75;
   };
   Bool_t IsSiLiPileup() const
   {
      return fPileupSiLi;
   };
   Bool_t IsChIoPileup() const
   {
      return fPileupChIo;
   };
   Bool_t UseFullChIoEnergyForCalib() const
   {
 
      return fUseFullChIoEnergyForCalib;
   };
   KVINDRAReconNuc();
   KVINDRAReconNuc(const KVINDRAReconNuc&);
   void init();
   virtual ~ KVINDRAReconNuc();
 
#if ROOT_VERSION_CODE >= ROOT_VERSION(3,4,0)
   virtual void Copy(TObject&) const;
#else
   virtual void Copy(TObject&);
#endif
   virtual void Clear(Option_t* t = "");
   void Print(Option_t* option = "") const;
 
   UInt_t GetRingNumber(void) const
   {
      if (GetTelescope()) {
         return GetTelescope()->GetRingNumber();
      }
      else {
         return 0;
      }
   };
   UInt_t GetModuleNumber(void) const
   {
      if (GetTelescope()) {
         return GetTelescope()->GetNumber();
      }
      else {
         return 0;
      }
   };
 
   virtual void Identify();
   virtual Bool_t CoherencySiCsI(KVIdentificationResult& theID);
   virtual Bool_t CoherencyChIoCsI(KVIdentificationResult& theID);
   virtual Bool_t CoherencyEtalons(KVIdentificationResult& theID);
   virtual Bool_t CoherencyChIoSiCsI(KVIdentificationResult);
   virtual void CalibrateRings1To9();
   virtual void CalibrateRings10To17();
   virtual void Calibrate();
 
   Float_t GetEnergyChIo()
   {
 
      return fEChIo;
   };
   Float_t GetEnergySi()
   {
 
      return fESi;
   };
   Float_t GetEnergySi75()
   {
 
      return fESi75;
   };
   Float_t GetEnergySiLi()
   {
 
      return fESiLi;
   };
   Float_t GetEnergyCsI()
   {
      return fECsI;
   };
 
   void SetEnergyChIo(Float_t val)
   {
 
      fEChIo = val;
   };
   void SetEnergySi(Float_t val)
   {
 
      fESi = val;
   };
   void SetEnergySi75(Float_t val)
   {
      fESi75 = val;
   };
   void SetEnergySiLi(Float_t val)
   {
      fESiLi = val;
   };
   void SetEnergyCsI(Float_t val)
   {
      fECsI = val;
   };
 
   KVChIo* GetChIo();
   KVSilicon* GetSi();
   KVSi75* GetSi75();
   KVSiLi* GetSiLi();
   KVCsI* GetCsI();
   Bool_t StoppedInChIo();
   Bool_t StoppedInSi();
   Bool_t StoppedInSi75();
   Bool_t StoppedInSiLi();
   Bool_t StoppedInCsI();
 
   KVINDRACodes& GetCodes() const
   {
      return fCodes;
   }
   void SetIDCode(UShort_t code_mask)
   {
      GetCodes().SetIDCode(code_mask);
      KVReconstructedNucleus::SetIDCode(GetCodes().GetVedaIDCode());
   }
 
   void SetECode(UChar_t code_mask)
   {
      GetCodes().SetECode(code_mask);
      KVReconstructedNucleus::SetECode(GetCodes().GetVedaECode());
   }
   Int_t GetIDCode() const
   {
      if (GetParameters()->HasIntParameter("IDCODE")) return GetParameters()->GetIntValue("IDCODE");
      return GetCodes().GetVedaIDCode();
   }
   Int_t GetECode() const
   {
      if (GetParameters()->HasIntParameter("ECODE")) return GetParameters()->GetIntValue("ECODE");
      return GetCodes().GetVedaECode();
   }
 
   Int_t GetIDSubCode(const Char_t* id_tel_type = "") const;
   const Char_t* GetIDSubCodeString(const Char_t* id_tel_type = "") const;
   KVINDRATelescope* GetTelescope() const
   {
      return (GetStoppingDetector() ?  dynamic_cast<KVINDRADetector*>(GetStoppingDetector())->GetTelescope() : 0);
   }
 
   ClassDef(KVINDRAReconNuc, 12) //Nucleus identified by INDRA array
};
 
#endif