KVDataPatch_E789_FAZIABadCsICalibFragments.h

#ifndef __KVDATAPATCH_E789_FAZIABADCSICALIBFRAGMENTS_H
#define __KVDATAPATCH_E789_FAZIABADCSICALIBFRAGMENTS_H
 
#include "KVDataPatch.h"
#include "KVReconstructedNucleus.h"
#include "KVFAZIA.h"
 
class KVDataPatch_E789_FAZIABadCsICalibFragments : public KVDataPatch {
public:
   KVDataPatch_E789_FAZIABadCsICalibFragments(): KVDataPatch()
   {
      SetName(ClassName());
      SetTitle(Class()->GetTitle());
   }
 
   Bool_t IsRequired(TString dataset, TString datatype, Int_t,
                     TString dataseries, Int_t datareleasenumber, const TList*) override
   {
      return (dataset == "INDRAFAZIA.E789") && (datatype == "recon") && (dataseries == "1.12") && (datareleasenumber == 5);
   }
 
   Bool_t IsEventPatch() override
   {
      return false;
   }
   Bool_t IsParticlePatch() override
   {
      return true;
   }
   void ApplyToEvent(KVReconstructedEvent*) override {}
   void ApplyToParticle(KVReconstructedNucleus* rnuc) override
   {
      if (rnuc->InArray("FAZIA") && rnuc->IsZMeasured() && rnuc->IsCalibrated()
            && rnuc->GetStoppingDetector()->IsType("CsI") && rnuc->GetZ() > 2)
         correct_ecsi(rnuc);
   }
 
   void PrintPatchInfo() const override;
 
   bool correct_ecsi(KVReconstructedNucleus* rnuc)
   {
 
      double ecsi, esi1, esi2;
      if ((ecsi = rnuc->GetParameters()->GetDoubleValue("FAZIA.ECSI")) > 0
            && (esi1 = rnuc->GetParameters()->GetDoubleValue("FAZIA.ESI1")) > 0
            && (esi2 = rnuc->GetParameters()->GetDoubleValue("FAZIA.ESI2")) > 0) {
         auto si1 = rnuc->GetReconstructionTrajectory()->GetDetector("SI1");
         auto si2 = rnuc->GetReconstructionTrajectory()->GetDetector("SI2");
         KVDetector si1si2("Si", si1->GetThickness() + si2->GetThickness());
         auto new_ecsi = si1si2.GetEResFromDeltaE(rnuc->GetZ(), rnuc->GetA(), esi1 + esi2);
         if (new_ecsi > 0) {
            rnuc->SetParameter("DATAPATCH.E789_FAZIABadCsICalibFragments.APPLIED", true);
            rnuc->SetParameter("FAZIA.ECSI_BAD", ecsi);
            rnuc->SetParameter("FAZIA.ECSI", -new_ecsi); // calculated CsI energy
            rnuc->SetECode(KVFAZIA::ECodes::SOME_ENERGY_LOSSES_CALCULATED);
            rnuc->SetEnergy(esi1 + esi2 + new_ecsi);
            rnuc->UpdateAllFrames();
            return true;
         }
      }
      return false;
   }
 
   ClassDefOverride(KVDataPatch_E789_FAZIABadCsICalibFragments, 1) //Correct mistaken use of Z=2 CsI calib for all fragments in FAZIA
};
 
#endif