KVINDRAGroupReconstructor.h

 
#ifndef __KVINDRAGROUPRECONSTRUCTOR_H
#define __KVINDRAGROUPRECONSTRUCTOR_H
 
#include "KVGroupReconstructor.h"
#include "KVINDRA.h"
 
class KVINDRAGroupReconstructor : public KVGroupReconstructor {
 
   friend class KVINDRA;
 
protected:
    KVDetector* theChio;
   Double_t fECsI, fESi, fEChIo;
    bool print_part;
   static TString CSI_ID_TYPE;
 
   void SetBadCalibrationStatus(KVReconstructedNucleus* n)
   {
      n->SetIsUncalibrated();
      n->SetECode(KVINDRA::ECodes::BAD_CALIBRATION);
      n->SetEnergy(0);
   }
   void SetNoCalibrationStatus(KVReconstructedNucleus* n)
   {
      n->SetIsUncalibrated();
      n->SetECode(KVINDRA::ECodes::NO_CALIBRATION_ATTEMPTED);
      n->SetEnergy(0);
   }
 
   double DoBeryllium8Calibration(KVReconstructedNucleus* n);
   void CheckCsIEnergy(KVReconstructedNucleus* n);
   KVDetector* GetCsI(KVReconstructedNucleus* n)
   {
      return n->GetReconstructionTrajectory()->GetDetector("CSI");
   }
 
   Bool_t CalculateChIoDEFromResidualEnergy(KVReconstructedNucleus* n, Double_t ERES);
 
   void CalibrateCoherencyParticle(KVReconstructedNucleus*) override {}
 
public:
   KVINDRAGroupReconstructor(const KVGroup* g = nullptr) : KVGroupReconstructor(g)
   {
      if (g) theChio = g->GetDetectorByType("CI");
   }
   virtual ~KVINDRAGroupReconstructor() {}
 
   void SetGroup(const KVGroup* g) override
   {
      KVGroupReconstructor::SetGroup(g);
      theChio = g->GetDetectorByType("CI");
   }
 
   KVReconstructedNucleus* ReconstructTrajectory(const KVGeoDNTrajectory* traj, const KVGeoDetectorNode* node) override;
 
   void Identify() override;
   void IdentifyParticle(KVReconstructedNucleus& PART) override;
   void CalibrateParticle(KVReconstructedNucleus* PART) override;
 
   virtual bool DoCoherencyAnalysis(KVReconstructedNucleus&) = 0;
 
   virtual void DoCalibration(KVReconstructedNucleus*) {}
 
   ClassDefOverride(KVINDRAGroupReconstructor, 1) //Reconstruct particles in INDRA groups
};
 
#endif