KVDetectionSimulator.h

 
#ifndef __KVDETECTIONSIMULATOR_H
#define __KVDETECTIONSIMULATOR_H
 
#include "KVBase.h"
#include "KVMultiDetArray.h"
#include "KVNucleusEvent.h"
#include "KVDetectorEvent.h"
#include "KVTarget.h"
#include "KVRangeTableGeoNavigator.h"
 
class KVDetectionSimulator : public KVBase {
 
   KVMultiDetArray* fArray = nullptr;        
    KVDetectorEvent fHitGroups;
   Bool_t           fCalcTargELoss = kTRUE;  
   Bool_t           fGeoFilter = kFALSE;     
   TString          fDetectionFrame;         
 
   KVRangeTableGeoNavigator* get_array_navigator() const
   {
      return static_cast<KVRangeTableGeoNavigator*>(fArray->GetNavigator());
   }
   KVNameValueList PropagateParticle(KVNucleus*);
public:
   KVDetectionSimulator() {}
   KVDetectionSimulator(KVMultiDetArray* a, Double_t cut_off = 1.e-3);
 
   void SetGeometricFilterMode()
   {
      SetIncludeTargetEnergyLoss(kFALSE);
      fGeoFilter = kTRUE;
   }
 
   void SetArray(KVMultiDetArray* a, Double_t e_cut_off = 1.e-3)
   {
      fArray = a;
      a->Clear();
      a->SetSimMode(kTRUE);
      get_array_navigator()->SetCutOffKEForPropagation(e_cut_off);
      if (GetTarget()) GetTarget()->SetOutgoing();
   }
   KVMultiDetArray* GetArray() const
   {
      return fArray;
   }
   KVTarget* GetTarget() const
   {
      return fArray->GetTarget();
   }
   void SetIncludeTargetEnergyLoss(Bool_t y = kTRUE)
   {
      fCalcTargELoss = y;
   }
   Bool_t IncludeTargetEnergyLoss() const
   {
      return fCalcTargELoss;
   }
   Double_t GetMinKECutOff() const
   {
      return get_array_navigator()->GetCutOffKEForPropagation();
   }
   void SetMinKECutOff(Double_t cutoff)
   {
      get_array_navigator()->SetCutOffKEForPropagation(cutoff);
   }
 
   void DetectEvent(KVEvent* event, const Char_t* detection_frame = "");
 
   void ClearHitGroups()
   {
      fHitGroups.Clear();
   }
 
   TString GetDetectionFrame() const
   {
 
      return fDetectionFrame;
   }
 
   ClassDef(KVDetectionSimulator, 0) //Simulate detection of particles or events in a detector array
};
 
#endif