KVFAZIA.h

 
#ifndef __KVFAZIA_H
#define __KVFAZIA_H
 
#include "KVMultiDetArray.h"
 
#include <KVGeoImport.h>
#include <KVEnv.h>
#include <KVSignal.h>
#include "KVFAZIATrigger.h"
 
#if ROOT_VERSION_CODE <= ROOT_VERSION(5,32,0)
#include "TGeoMatrix.h"
#endif
 
class KVDetectorEvent;
#ifdef WITH_PROTOBUF
#ifndef __CINT__
namespace DAQ {
   class FzEvent;
}
#endif
#endif
 
class KVFAZIA : public KVMultiDetArray {
protected:
    TString fFGeoType;
    Double_t fFDist;
    Double_t fFThetaMin;
    Int_t fNblocks;
    Int_t fStartingBlockNumber;
    TString fCorrespondanceFile;
   KVString fDetectorLabels;
   KVString fSignalTypes;
    Double_t fImport_dTheta;
    Double_t fImport_dPhi;
    Double_t fImport_ThetaMin;
    Double_t fImport_ThetaMax;
    Double_t fImport_PhiMin;
    Double_t fImport_PhiMax;
    Double_t fImport_Xorg;
    Double_t fImport_Yorg;
    Double_t fImport_Zorg;
    int fQuartet[8][2];
    int fTelescope[8][2];
 
   Double_t fQH1risetime;
   Double_t fQ2risetime;
   Double_t fQ3slowrisetime;
   Double_t fQ3fastrisetime;
 
   KVFAZIATrigger fTrigger;
 
   double oldfaziats = -1.;
 
   TString fSiCsI_idtype = "SI2_CSI";
 
   void SetTriggerPatternsForDataSet(const TString& dataset);
   void SetTriggerPattern(uint16_t fp)
   {
      fTrigger.SetTriggerPattern(fp);
   }
 
   virtual void BuildFAZIA();
   virtual void GetGeometryParameters();
 
   virtual void BuildTarget();
   void GenerateCorrespondanceFile();
   virtual void SetNameOfDetectors(KVEnv& env);
 
   virtual void DefineStructureFormats(KVGeoImport&) {}
 
#ifdef WITH_MFM
   Bool_t handle_raw_data_event_mfmframe(const MFMCommonFrame&) override;
#endif
   void copy_fired_parameters_to_recon_param_list() override
   {
   }
 
   void PerformClosedROOTGeometryOperations() override;
 
   void CreateCorrespondence();
#ifdef WITH_PROTOBUF
   Bool_t handle_raw_data_event_protobuf(KVProtobufDataReader&) override;
#ifndef __CINT__
   Bool_t treat_event(const DAQ::FzEvent&);
#endif
   Double_t TreatEnergy(Int_t sigid, Int_t eid, UInt_t val);
#endif
   TString GetSignalName(Int_t bb, Int_t qq, Int_t tt, Int_t idsig);
 
   void ReadTriggerPatterns(KVExpDB* db);
 
   void prepare_to_handle_new_raw_data()
   {
      KVMultiDetArray::prepare_to_handle_new_raw_data();
      SetTriggerPattern(0);
   }
   void SetIDCodeForIDTelescope(KVIDTelescope*) const override;
 
   UChar_t GetNormalCalibrationCode() const override
   {
      return ECodes::NORMAL_CALIBRATION;
   }
public:
   enum IDCodes {
      NO_IDENTIFICATION = -1,    
      ID_GAMMA = 0,    
      ID_SI1_PSA = 11, 
      ID_SI1_SI2 = 12, 
      ID_SI1_CSI = 13, 
      ID_SI2_PSA = 22, 
      ID_SI2_CSI = 23, 
      ID_SI12_CSI = 123, 
      ID_CSI_PSA = 33, 
      ID_STOPPED_IN_FIRST_STAGE = 5, 
      ID_SI1_SI2_MAYBE_PUNCH_THROUGH = 120, 
      ID_SI1_SI2_PUNCH_THROUGH = 121, 
      ID_SI1_PUNCH_THROUGH = 110, 
      ID_INCOHERENT = 666 
   };
   Int_t GetIDCodeForParticlesStoppingInFirstStageOfTelescopes() const override
   {
      return IDCodes::ID_STOPPED_IN_FIRST_STAGE;
   }
   Short_t GetBadIDCode() const override
   {
      return IDCodes::NO_IDENTIFICATION;
   }
 
   TString GetIDCodeMeaning(Int_t idc) const override
   {
      switch (idc) {
         case NO_IDENTIFICATION:
            return "NO_IDENTIFICATION";
         case ID_GAMMA:
            return "ID_GAMMA";
         case ID_SI1_PSA:
            return "ID_SI1_PSA";
         case ID_SI2_PSA:
            return "ID_SI2_PSA";
         case ID_SI1_SI2:
            return "ID_SI1_SI2";
         case ID_SI2_CSI:
            return "ID_SI2_CSI";
         case ID_SI1_CSI:
            return "ID_SI1_CSI";
         case ID_SI12_CSI:
            return "ID_SI12_CSI";
         case ID_CSI_PSA:
            return "ID_CSI_PSA";
         case ID_STOPPED_IN_FIRST_STAGE:
            return "ID_STOPPED_IN_FIRST_STAGE";
         case ID_SI1_SI2_MAYBE_PUNCH_THROUGH:
            return "ID_SI1_SI2_MAYBE_PUNCH_THROUGH";
         case ID_SI1_SI2_PUNCH_THROUGH:
            return "ID_SI1_SI2_PUNCH_THROUGH";
         case ID_SI1_PUNCH_THROUGH:
            return "ID_SI1_PUNCH_THROUGH";
         default:
            return Form("(unknown:%d)", idc);
      }
   }
   enum ECodes {
      NO_CALIBRATION_ATTEMPTED = 0, 
      NORMAL_CALIBRATION = 1,       
      SOME_ENERGY_LOSSES_CALCULATED = 2, 
      WARNING_CSI_MAX_ENERGY = 3,    
      ENERGY_LOSSES_TENTATIVELY_CALCULATED = 4    
   };
   TString GetECodeMeaning(Int_t ec) const override
   {
      switch (ec) {
         case NO_CALIBRATION_ATTEMPTED:
            return "NO_CALIBRATION_ATTEMPTED";
         case NORMAL_CALIBRATION:
            return "NORMAL_CALIBRATION";
         case SOME_ENERGY_LOSSES_CALCULATED:
            return "SOME_ENERGY_LOSSES_CALCULATED";
         case WARNING_CSI_MAX_ENERGY:
            return "WARNING_CSI_MAX_ENERGY";
         case ENERGY_LOSSES_TENTATIVELY_CALCULATED:
            return "ENERGY_LOSSES_TENTATIVELY_CALCULATED (for Z=1,2 stopped in uncalibrated CsI)";
         default:
            return Form("(unknown:%d)", ec);
      }
   }
 
   KVFAZIA(const Char_t* title = "");
   virtual ~KVFAZIA();
   void AddDetectorLabel(const Char_t* label);
 
   void Build(Int_t = -1) override;
 
   void GetDetectorEvent(KVDetectorEvent* detev, const TSeqCollection* dets) override;
   Int_t GetNumberOfBlocks() const
   {
      return fNblocks;
   }
   void IncludeTargetInGeometry(Bool_t include = kTRUE)
   {
      fBuildTarget = include;
   }
 
   KVString GetDetectorLabels() const
   {
      return fDetectorLabels;
   }
   const Char_t* GetSignalTypes() const
   {
      return fSignalTypes.Data();
   }
 
   void SetGeometryImportParameters(Double_t dt = 0.25, Double_t dp = 1.0, Double_t tmin = 2., Double_t pmin = 0, Double_t tmax = 20., Double_t pmax = 360.,
                                    Double_t xorg = 0, Double_t yorg = 0, Double_t zorg = 0)
   {
      fImport_dPhi = dp;
      fImport_dTheta = dt;
      fImport_PhiMax = pmax;
      fImport_PhiMin = pmin;
      fImport_ThetaMax = tmax;
      fImport_ThetaMin = tmin;
      fImport_Xorg = xorg;
      fImport_Yorg = yorg;
      fImport_Zorg = zorg;
   }
   void FillDetectorList(KVReconstructedNucleus* rnuc, KVHashList* DetList, const KVString& DetNames) override;
 
   KVGroupReconstructor* GetReconstructorForGroup(const KVGroup*) const override;
   Double_t GetSetupParameter(const Char_t* parname);
 
   const KVFAZIATrigger& GetTrigger() const
   {
      return fTrigger;
   }
   void SetRawDataFromReconEvent(const KVReconstructedEvent& recev) override;
   void MakeCalibrationTables(KVExpDB*) override;
 
   std::string GetTriggerForCurrentRun() const;
   void DeduceIdentificationTelescopesFromGeometry() override;
 
   static TGeoHMatrix* GetFAZIAVolumePositioningMatrix(Double_t distance, Double_t theta, Double_t phi, TGeoTranslation* postTrans = nullptr)
   {
 
      TGeoRotation rot1, rot2;
      TGeoTranslation trans;
      phi += 90;
      rot1.SetAngles(-phi, 0, 0) ;
      rot2.SetAngles(phi, theta, 0) ;
      trans.SetDz(distance) ;
      TGeoHMatrix h;
      if (postTrans) h = (*postTrans) * rot2 * trans * rot1 ;
      else h = rot2 * trans * rot1;
      TGeoHMatrix* ph = new TGeoHMatrix(h);
      return ph;
   }
   void AssociateTrajectoriesAndNodes() override;
 
   TString GetSiCsIIDType() const
   {
      return fSiCsI_idtype;
   }
 
   ClassDefOverride(KVFAZIA, 1) //Base class for description of the FAZIA set up
};
 
R__EXTERN KVFAZIA* gFazia;
 
#endif