KVDataPatch_INDRA_camp4_ring10_He_isotopic_identification.cpp

#include "KVDataPatch_INDRA_camp4_ring10_He_isotopic_identification.h"
#include "KVINDRA.h"
 
ClassImp(KVDataPatch_INDRA_camp4_ring10_He_isotopic_identification)
 
 
 
KVDataPatch_INDRA_camp4_ring10_He_isotopic_identification::
KVDataPatch_INDRA_camp4_ring10_He_isotopic_identification()
{
   // Default constructor
   SetName(ClassName());
   SetTitle(Class()->GetTitle());
}
 
 
 
 
void KVDataPatch_INDRA_camp4_ring10_He_isotopic_identification::ApplyToParticle(KVReconstructedNucleus* rnuc)
{
   // For Z=2 particles identified in CSI_R_L on ring 10, we check the mass attribution
   // using new identification grids drawn by Seon Ho Nam
 
   if (rnuc->GetStoppingDetector()->GetIndex() / 100 == 10 && rnuc->GetZ() == 2 && rnuc->GetIDCode() == KVINDRA::ID_CSI_PSA) {
      auto old_A = rnuc->GetA();
      auto old_isAmeasured = rnuc->IsAMeasured();
      auto idt = rnuc->GetIdentifyingTelescope();
      auto idr = rnuc->GetIdentificationResult(idt);
      if (idt->Identify(idr)) {
         rnuc->SetIdentification(idr, idt);
         rnuc->SetParameter("HE_ISOTOPE_CORRECTION", true);
         rnuc->SetParameter("HE_ISOTOPE_CORRECTION.OLD_A", old_A);
         rnuc->SetParameter("HE_ISOTOPE_CORRECTION.OLD_ISAMEASURED", old_isAmeasured);
      }
   }
}
 
 
 
void KVDataPatch_INDRA_camp4_ring10_He_isotopic_identification::PrintPatchInfo() const
{
   std::cout << "Correct bad attribution of isotopic mass of He nuclei in ring 10" << std::endl;
   std::cout << std::endl;
   std::cout << "In original DST generation, many low energy 4He nuclei were" << std::endl;
   std::cout << "incorrectly identified as 3He. Using new identification grids" << std::endl;
   std::cout << "for Z=2 nuclei we re-assign the masses for thes particles." << std::endl;
   std::cout << "We give single KVNucleus::kVedaMass A to all" << std::endl;
}