KVDataPatch_INDRA_CorrectEtalonTrajectories.cpp

//Created by KVClassFactory on Fri Feb  7 13:42:32 2020
//Author: John Frankland,,,
 
#include "KVDataPatch_INDRA_CorrectEtalonTrajectories.h"
#include "KVDetectorStack.h"
#include "KVReconstructedEvent.h"
#include "KVINDRADetector.h"
 
ClassImp(KVDataPatch_INDRA_CorrectEtalonTrajectories)
 
 
 
KVDataPatch_INDRA_CorrectEtalonTrajectories::KVDataPatch_INDRA_CorrectEtalonTrajectories()
   : KVDataPatch(), fRandomAngles(false)
{
   // Default constructor
   SetName(ClassName());
   SetTitle(Class()->GetTitle());
}
 
 
 
 
void KVDataPatch_INDRA_CorrectEtalonTrajectories::ApplyToEvent(KVReconstructedEvent* e)
{
   // Will be called before ApplyToParticle
   //
   // Just need to know if event has mean or random angles
 
   fRandomAngles = e->HasRandomAngles();
}
 
 
 
 
void KVDataPatch_INDRA_CorrectEtalonTrajectories::ApplyToParticle(KVReconstructedNucleus* N)
{
   // Apply the following correction to INDRA data:
   //    all particles identified in etalon modules by either CSI_R_L or CI_CSI
   //    identification were treated as if they passed through the etalon detectors.
   //    we check the coherency of this with the energy of the particle and
   //    correct when necessary
 
   if (!N->IsOK()) return; //only apply to well-identified/calibrated particles
 
   int ring = N->GetStoppingDetector()->GetIndex() / 100;
 
   if (ring < 10) return;
 
   KVINDRADetector* csi;
   if (!(csi = (KVINDRADetector*)N->GetStoppingDetector())->IsType("CSI")) return;
 
   auto si75 = (KVINDRADetector*)N->GetReconstructionTrajectory()->GetDetector("SI75");
   auto sili = (KVINDRADetector*)N->GetReconstructionTrajectory()->GetDetector("SILI");
   if(si75 || sili)
   {
      // check particle's energy is larger than punch through for CI-SI75-SILI
      KVDetectorStack etalon_stack;
      etalon_stack.Add(csi->GetChIo());
      if(si75)
         etalon_stack.Add(si75);
      else
      {Info("ApplyToParticle", "si75=%p sili=%p", si75, sili); N->GetReconstructionTrajectory()->Print();}
      if(sili)
         etalon_stack.Add(sili);
      else
      {Info("ApplyToParticle", "si75=%p sili=%p", si75, sili); N->GetReconstructionTrajectory()->Print();}
      if (N->GetE() < etalon_stack.GetPunchThroughEnergy(N->GetZ(), N->GetA())) {
         // it is not possible: this particle went straight from CI to CSI
         // change trajectory
         auto traj = (KVGeoDNTrajectory*)csi->GetNode()->GetTrajectories()->FindObjectWithMethod("2", "GetN");
         N->ReplaceReconTraj(traj->GetName());
         // redefine particle angles from new trajectory
         if (fRandomAngles) N->GetAnglesFromReconstructionTrajectory("random");
         else N->GetAnglesFromReconstructionTrajectory("mean");
 
      }
   }
}
 
 
 
 
void KVDataPatch_INDRA_CorrectEtalonTrajectories::PrintPatchInfo() const
{
   std::cout <<  "Apply the following correction to INDRA data:\n";
   std::cout <<  "   all particles identified in etalon modules by either CSI_R_L or CI_CSI\n";
   std::cout <<  "   identification were treated as if they passed through the etalon detectors.\n";
   std::cout <<  "   we check the coherency of this with the energy of the particle and\n";
   std::cout <<  "   correct when necessary\n";
}