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*) {}
 
public:
   KVINDRAGroupReconstructor() {}
   virtual ~KVINDRAGroupReconstructor() {}
 
   void SetGroup(KVGroup* g)
   {
      KVGroupReconstructor::SetGroup(g);
      theChio = g->GetDetectorByType("CI");
   }
 
   KVReconstructedNucleus* ReconstructTrajectory(const KVGeoDNTrajectory* traj, const KVGeoDetectorNode* node);
 
   void Identify();
   void IdentifyParticle(KVReconstructedNucleus& PART);
   void CalibrateParticle(KVReconstructedNucleus* PART);
 
   virtual bool DoCoherencyAnalysis(KVReconstructedNucleus&) = 0;
 
   virtual void DoCalibration(KVReconstructedNucleus*) {}
 
   ClassDef(KVINDRAGroupReconstructor, 1) //Reconstruct particles in INDRA groups
};
 
#endif