DPSCPair_factory.cc

Go to the documentation of this file.

// $Id$
//
//    File: DPSCPair_factory.cc
// Created: Tue Mar 24 21:35:49 EDT 2015
// Creator: nsparks (on Linux cua2.jlab.org 2.6.32-431.5.1.el6.x86_64 x86_64)
//


#include <iostream>
#include <iomanip>
#include <math.h>
using namespace std;

#include "DPSCPair_factory.h"
#include "DPSCHit.h"
using namespace jana;

inline bool DPSCPair_SortByTimeDifference(const DPSCPair* pair1, const DPSCPair* pair2)
{
  double tdiff1 = fabs(pair1->ee.first->t-pair1->ee.second->t);
  double tdiff2 = fabs(pair2->ee.first->t-pair2->ee.second->t);
  return (tdiff1<tdiff2);
}

//------------------
// init
//------------------
jerror_t DPSCPair_factory::init(void)
{
  DELTA_T_PAIR_MAX = 10.0; // ns
  gPARMS->SetDefaultParameter("PSCPair:DELTA_T_PAIR_MAX",DELTA_T_PAIR_MAX,
               "Maximum difference in ns between a pair of hits"
               " in left and right arm of coarse PS");
  return NOERROR;
}

//------------------
// brun
//------------------
jerror_t DPSCPair_factory::brun(jana::JEventLoop *eventLoop, int32_t runnumber)
{
  return NOERROR;
}

//------------------
// evnt
//------------------
jerror_t DPSCPair_factory::evnt(JEventLoop *loop, uint64_t eventnumber)
{
  // get coarse pair spectrometer hits
  vector<const DPSCHit*> hits;
  loop->Get(hits);
  // form PSC left-right hit pairs and sort by time difference
  pair<const DPSCHit*,const DPSCHit*> ee;
  if (hits.size()>1) {
    for (unsigned int i=0; i < hits.size()-1; i++) {
      for (unsigned int j=i+1; j < hits.size(); j++) {
   if (!hits[i]->has_TDC||!hits[j]->has_TDC) continue;
   if (!hits[i]->has_fADC||!hits[j]->has_fADC) continue;
   if (std::abs(hits[i]->arm-hits[j]->arm)==1&&fabs(hits[i]->t-hits[j]->t)<DELTA_T_PAIR_MAX) {
     if (hits[i]->arm==0) {
       ee.first = hits[i];
       ee.second = hits[j];
     }
     else if (hits[i]->arm==1) {
       ee.first = hits[j];
       ee.second = hits[i];
     }
     DPSCPair *pair = new DPSCPair;
     pair->ee = ee;
     _data.push_back(pair);
   }
      }
    }
  }
  sort(_data.begin(),_data.end(),DPSCPair_SortByTimeDifference);

  return NOERROR;
}

//------------------
// erun
//------------------
jerror_t DPSCPair_factory::erun(void)
{
  return NOERROR;
}

//------------------
// fini
//------------------
jerror_t DPSCPair_factory::fini(void)
{
  return NOERROR;
}