JEventProcessor_ST_Propagation_Time.h

Go to the documentation of this file.

// $Id$
//
//    File: JEventProcessor_ST_Propagation_Time.h
// Created: Thu Jan 14 12:16:20 EST 2016
// Creator: mkamel (on Linux ifarm1401 2.6.32-431.el6.x86_64 x86_64)
//

#ifndef _JEventProcessor_ST_Propagation_Time_
#define _JEventProcessor_ST_Propagation_Time_

#include <JANA/JEventProcessor.h>
// ROOT header files
#include <TMath.h>
#include <TDirectory.h>
#include <TH1.h>
#include <TH2.h>
#include <TTree.h>
// C++ header files
#include <stdint.h>
#include <vector>
#include <stdio.h>
// PID libraries
#include <PID/DEventRFBunch.h>
#include <PID/DParticleID.h>
#include <PID/DDetectorMatches.h>
#include <PID/DParticleID.h>
#include <PID/DChargedTrack.h>
// Tracking libraries
#include <TRACKING/DTrackTimeBased.h>
// RF libraries
#include <RF/DRFTDCDigiTime.h>
#include <RF/DRFTime_factory.h>
// ST libraries
#include <START_COUNTER/DSCHit.h>
#include <START_COUNTER/DSCDigiHit.h>
#include <START_COUNTER/DSCTDCDigiHit.h>
#include <START_COUNTER/DSCTruthHit.h>
// TOF libraries
#include <TOF/DTOFHit.h>
// DAQ/EPICS libraries
#include <DAQ/DEPICSvalue.h>
// Translation table libraries
#include <TTAB/DTTabUtilities.h>
#include <TTAB/DTranslationTable.h>
//
const Int_t NCHANNELS = 30;
class JEventProcessor_ST_Propagation_Time:public jana::JEventProcessor{
   public:
      JEventProcessor_ST_Propagation_Time();
      ~JEventProcessor_ST_Propagation_Time();
      const char* className(void){return "JEventProcessor_ST_Propagation_Time";}

   private:
      jerror_t init(void);                ///< Called once at program start.
      jerror_t brun(jana::JEventLoop *eventLoop, int32_t runnumber); ///< Called everytime a new run number is detected.
      jerror_t evnt(jana::JEventLoop *eventLoop, uint64_t eventnumber); ///< Called every event.
      jerror_t erun(void);                ///< Called everytime run number changes, provided brun has been called.
      jerror_t fini(void);                ///< Called after last event of last event source has been processed.
      // Define constants
      const DParticleID* dParticleID;
      DRFTime_factory *dRFTimeFactory;
      double z_target_center;  // Target center along z
      double dRFBunchPeriod;   // RF bunch period from CCDBr
      //////////////////////////////////
      //double locSCHitTime,       locSCTrackFlightTime,   locFlightTimeCorrectedSCTime;
      //double locTOFHitTime,      locTOFTrackFlightTime,  locFlightTimeCorrectedTOFTime;
      //double locCenteredRFTime,  locCenterToVertexRFTime,  locVertexRFTime, locPeriod;
      //double locTOFRFShiftedTime;
      ///double locSCzIntersection;
      //double locSCPropTime;
      //double sc_pos_soss, sc_pos_eoss, sc_pos_eobs, sc_pos_eons;
      //double  Corr_Time,Corr_Time_ss,Corr_Time_bs,Corr_Time_ns, Corr_Time_bn;
      //double time_lower_limit, time_upper_limit, z_lower_limit, z_upper_limit;
      //int NoBins_time, NoBins_z;

      vector<vector<DVector3> >sc_pos;   // SC geometry vector
      vector<vector<DVector3> >sc_norm;  
      // Grab match track detector parameters
      //DTOFHitMatchParams locTOFHitMatchParams;  // TOF
      //DSCHitMatchParams  locSCHitMatchParams;   // SC
      // Declare event vertex vector
      //DVector3 vertex;
      // Declare a vector which quantizes the point of the intersection of a charged particle 
      //   with a plane in the middle of the scintillator 
      //DVector3 IntersectionPoint;
      // Declare a vector which quantizes the unit vector of the charged particle track traversing
      //   through the scintillator with its origin at the intersection point
      //DVector3 IntersectionDir;
      // Grab the paramteres associated to a track matched to the ST
      //vector<DSCHitMatchParams> st_params;
      // Declare r and z coord. of track vertex, sc array index
      //double z_v, r_v;
      //int sc_index;
      double Photonspeed;
      // Cuts
      double trackingFOMCut;
      double pim_pmag_cut;
        double sc_angle_corr;
      //bool z_vertex_cut, r_vertex_cut;
      //bool foundTOF, foundSC, foundSCandTOF;
      //bool sc_match, sc_match_pid;
      // ******************* 2D histos **************
      TH2I **h2_PropTime_z_SS_chan;
      TH2I **h2_PropTime_z_BS_chan;
      TH2I **h2_PropTime_z_NS_chan;
      TH2I **h2_PpropTime_z;
      
      TH2I **h2_PropTimeCorr_z_SS_chan;
      TH2I **h2_PropTimeCorr_z_BS_chan;
      TH2I **h2_PropTimeCorr_z_NS_chan;
      
      TH2I **h2_CorrectedTime_z;
      //Define Calibration parameters variable called from CCDB
      vector<vector<double> >propagation_time_corr;
      vector<vector<double> >PTC_Boundary;
};
#endif // _JEventProcessor_ST_Propagation_Time_