KVINDRAFilterEtalonGroupReconstructor.cpp

#include "KVINDRAFilterEtalonGroupReconstructor.h"
 
 
 
KVINDRAFilterEtalonGroupReconstructor::KVINDRAFilterEtalonGroupReconstructor(const KVGroup* g)
   : KVINDRAFilterGroupReconstructor(g)
{
   // Analysis of etalon group is carried out in order to identify:
   //   * the CSI behind each etalon telescope
   //   * the long ('etalon') trajectory for each CSI (CSI/SILI/SI75/CI/)
   //   * the short ('not_etalon') trajectory (CSI/CI)
   //   * the SILI detector on each etalon trajectory
   //
   TIter it_traj(GetGroup()->GetTrajectories());
   KVGeoDNTrajectory* traj;
   // first find trajectory with etalons + csi behind etalons
   while ((traj = (KVGeoDNTrajectory*)it_traj())) {
      traj->IterateFrom();
      auto det = traj->GetNextNode()->GetDetector();
      auto ring_number = dynamic_cast<KVINDRADetector*>(det)->GetRingNumber();
      if (traj->Contains(Form("SILI_%d", ring_number))) {
         etalons[ring_number].etalon_trajectory = traj;
         etalons[ring_number].csi_etalon = det->GetNode();
         etalon_csi_rings[det->GetNode()] = ring_number;
         etalons[ring_number].sili = traj->GetNode(Form("SILI_%d", ring_number))->GetDetector();
      }
   }
   it_traj.Reset();
   while ((traj = (KVGeoDNTrajectory*)it_traj())) {
      for (auto& etal : etalons) {
         auto& etalel = etal.second;
         if (traj != etalel.etalon_trajectory) {
            if (traj->BeginsAt(etalel.csi_etalon)) {
               etalel.not_etalon_trajectory = traj;
            }
         }
      }
   }
   // now reiterate over trajectories, put all except the etalon trajectories in temporary list
   TList trajectories;
   it_traj.Reset();
   while ((traj = (KVGeoDNTrajectory*)it_traj())) {
      bool not_an_etalon = true;
      for (auto& etal : etalons) {
         if (traj == etal.second.etalon_trajectory || traj == etal.second.not_etalon_trajectory) {
            not_an_etalon = false;
         }
      }
      if (not_an_etalon) trajectories.Add(traj);
   }
   // finally add etalon trajectories to list & to map
   for (auto& etal : etalons) {
      trajectories.Add(etal.second.etalon_trajectory);
      trajectories.Add(etal.second.not_etalon_trajectory);
   }
   // now replace group's list of trajectories with this one
   GetGroup()->ReplaceTrajectories(&trajectories);
}
 
 
ClassImp(KVINDRAFilterEtalonGroupReconstructor)