KVDataPatch_E789_INDRANewCsICalib.h

#ifndef __KVDATAPATCH_E789_INDRANEWCSICALIB_H
#define __KVDATAPATCH_E789_INDRANEWCSICALIB_H
 
#include "KVDataPatch.h"
#include "KVEventReconstructor.h"
 
class KVDataPatch_E789_INDRANewCsICalib : public KVDataPatch {
   void correct_ecsi(KVReconstructedNucleus*);
   std::unique_ptr<KVEventReconstructor> fReconstructor;
   std::unique_ptr<KVReconstructedEvent> fReconEvent;
public:
   KVDataPatch_E789_INDRANewCsICalib()
      : KVDataPatch()
   {
      SetName(ClassName());
      SetTitle(Class()->GetTitle());
   }
 
   Bool_t IsRequired(TString dataset, TString datatype, Int_t,
                     TString dataseries, Int_t datareleasenumber, const TList*) override
   {
      auto required = (dataset == "INDRAFAZIA.E789") && (datatype == "recon") && (dataseries == "1.12" || dataseries == "1.13");
 
      if (required) {
         fReconEvent.reset(new KVReconstructedEvent);
         fReconstructor.reset(new KVEventReconstructor(gMultiDetArray, fReconEvent.get()));
      }
      return required;
   }
 
   Bool_t IsEventPatch() override
   {
      return false;
   }
   Bool_t IsParticlePatch() override
   {
      return true;
   }
   void ApplyToEvent(KVReconstructedEvent*) override {}
   void ApplyToParticle(KVReconstructedNucleus* rnuc) override
   {
      if (rnuc->InArray("INDRA") && rnuc->IsZMeasured() && (rnuc->GetZ() > 0) && rnuc->GetStoppingDetector()->IsType("CSI"))
         correct_ecsi(rnuc);
   }
 
   void PrintPatchInfo() const override;
 
   ClassDefOverride(KVDataPatch_E789_INDRANewCsICalib, 1) //Apply new CsI calibration for particles detected in INDRA
};
 
#endif