/home/alex/Work/Research/Km3NeT/PPM-DU/ppm-du-v2.3/triggeringEvent.cc

Go to the documentation of this file.

#include <cmath>
#include <algorithm>
#include <iostream>

#include "triggeringEvent.h"

void
AddDom2Event(std::vector<triggeredSignals>& vTriggeredDom,
             std::vector<triggeredSignals>& vEventList)
{
  vEventList.insert(vEventList.end(), vTriggeredDom.begin(), vTriggeredDom.end());
}

void
Do2DomTrigger(std::vector<triggeredSignals>& vTrigSignals,
               std::vector<triggeredSignals>& vFinalEvent,
               const int domTimeWindow)
{
  //std::cout << "test0" << std::endl;
  sort(vTrigSignals.begin(), vTrigSignals.end());
  //std::cout << "test1" << std::endl;
  std::vector<signals> vEventSignals;
  vEventSignals.clear();
  for (unsigned int dom1 = 0; dom1 < vTrigSignals.size() - 1; ++dom1) {
    if (vTrigSignals.size() == 0) break;
    const int domId1 = vTrigSignals[dom1].GetSignals()[0].GetDomId();
    int skipped = 0;
    for (unsigned int dom2 = dom1 + 1; dom2 < vTrigSignals.size(); ++dom2) {
      const int domId2 = vTrigSignals[dom2].GetSignals()[0].GetDomId();
      //std::cout << domId1 << ' ' << domId2 << ' ' << dom1 << ' ' << dom2 << ' ' << vTrigSignals.size() << std::endl;
      if (domId1 == domId2) continue;
      const long long int time1 = vTrigSignals[dom1].GetSignals()[0].GetStart();
      const long long int time2 = vTrigSignals[dom2].GetSignals()[0].GetStart();
      //std::cout << domId1 << ' ' << domId2 << ' ' << time2 - time1 << std::endl;
      if (fabs(time1 - time2) > domTimeWindow) break;
      for (unsigned int j = 0; j < vTrigSignals[dom1].GetSignals().size(); ++j)
        vEventSignals.push_back(vTrigSignals[dom1].GetSignals()[j]);
      for (unsigned int j = 0; j < vTrigSignals[dom2].GetSignals().size(); ++j)
        vEventSignals.push_back(vTrigSignals[dom2].GetSignals()[j]);
      triggeredSignals trigEv("");
      trigEv.SetSignals(vEventSignals);
      vFinalEvent.push_back(trigEv);
      vEventSignals.clear();
      skipped = dom2 - dom1 - 1;
      break;
    }
    if (dom1 + skipped >= vTrigSignals.size()) break;
    else dom1 += skipped;
  }  
}


void
DoEventTrigger(std::vector<triggeredSignals>& vEventList,
               std::vector<triggeredSignals>& vFinalEvent,
               const int domTimeWindow)
{
  sort(vEventList.begin(), vEventList.end());
  std::vector<signals> vEventSignals;
  vEventSignals.clear();
  for (unsigned int i = 0; i < vEventList.size() - 1; ++i) {
    if (vEventList.size() == 0) break;
    const int dom1 = vEventList[i].GetSignals()[0].GetDomId();
    const int dom2 = vEventList[i + 1].GetSignals()[0].GetDomId();
    if (dom1 != dom2) {
      const long long int time1 = vEventList[i].GetSignals()[0].GetStart();
      const long long int time2 = vEventList[i + 1].GetSignals()[0].GetStart();
      if (fabs(time1 - time2) < domTimeWindow) {
        for (unsigned int j = 0; j < vEventList[i].GetSignals().size(); ++j)
          vEventSignals.push_back(vEventList[i].GetSignals()[j]);
        for (unsigned int j = 0; j < vEventList[i + 1].GetSignals().size(); ++j)
          vEventSignals.push_back(vEventList[i + 1].GetSignals()[j]);
      } else {
        triggeredSignals trigEv("");
        trigEv.SetSignals(vEventSignals);
        vFinalEvent.push_back(trigEv);
        vEventSignals.clear();      
      }
    } else {
      triggeredSignals trigEv("");
      trigEv.SetSignals(vEventSignals);
      vFinalEvent.push_back(trigEv);
      vEventSignals.clear();      
    }
  }  
  //std::cout << vEventSignals.size() << std::endl;
  //triggeredSignals trigEv("");
  //trigEv.SetSignals(vEventSignals);
  //vFinalEvent.push_back(trigEv);
  //vEventSignals.clear();
  //std::cout << vFinalEvent.size() << std::endl;
}