KVINDRAFilterEtalonGroupReconstructor.h

#ifndef __KVINDRAFILTERETALONGROUPRECONSTRUCTOR_H
#define __KVINDRAFILTERETALONGROUPRECONSTRUCTOR_H
 
#include "KVINDRAFilterGroupReconstructor.h"
 
class KVINDRAFilterEtalonGroupReconstructor : public KVINDRAFilterGroupReconstructor {
   struct etalon_telescope {
      KVGeoDetectorNode* csi_etalon = nullptr;
      KVDetector* sili = nullptr;
      KVGeoDNTrajectory* etalon_trajectory = nullptr;
      KVGeoDNTrajectory* not_etalon_trajectory = nullptr;
   };
    std::map<int, etalon_telescope> etalons;
    std::map<const KVGeoDetectorNode*, int> etalon_csi_rings;
public:
   KVINDRAFilterEtalonGroupReconstructor(const KVGroup* g = nullptr);
 
   KVReconstructedNucleus* ReconstructTrajectory(const KVGeoDNTrajectory* traj, const KVGeoDetectorNode* node) override
   {
 
      KVReconstructedNucleus* rnuc = nullptr;
      if (auto ring = etalon_csi_rings[node]) { // is node the CSI behind the etalons?
         if (traj == etalons[ring].etalon_trajectory) { // are we on the etalon trajectory?
            if (node->GetDetector()->Fired(GetPartSeedCond()) && !node->GetDetector()->IsAnalysed() && etalons[ring].sili->Fired()) {
               ++nfireddets;
               rnuc = GetEventFragment()->AddParticle();
            }
            else
               rnuc = nullptr;
         }
         else { // on short CSI/CI/ trajectory
            if (node->GetDetector()->Fired(GetPartSeedCond()) && !node->GetDetector()->IsAnalysed()) {
               ++nfireddets;
               rnuc = GetEventFragment()->AddParticle();
            }
            else
               rnuc = nullptr;
         }
      }
      else // all other stopping detectors/trajectories use standard treatment
         rnuc = KVINDRAFilterGroupReconstructor::ReconstructTrajectory(traj, node);
      return rnuc;
   }
 
   ClassDefOverride(KVINDRAFilterEtalonGroupReconstructor, 0)
};
 
#endif