BesTrigL1.h

Go to the documentation of this file.

#ifndef BesTrigL1_h
#define BesTrigL1_h 1
 
// from Gaudi
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/INTuple.h"
#include "GaudiKernel/INTupleSvc.h"
#include "GaudiKernel/ITHistSvc.h"
#include "GaudiKernel/NTuple.h"
#include "MdcGeomSvc/IMdcGeomSvc.h"
#include "RawDataProviderSvc/IRawDataProviderSvc.h"
#include "Trigger/EmcTCFinder.h"
#include "Trigger/MdcTSF.h"
#include "Trigger/MucTrigHit.h"
#include "Trigger/TofHitCount.h"
#include <fstream>
 
#include "EmcCalibConstSvc/IEmcCalibConstSvc.h"
#include "RealizationSvc/IRealizationSvc.h"
#include "RealizationSvc/RealizationSvc.h"
 
#include "TH1F.h"
#include "TObjArray.h"
#include "TProfile.h"
#include "TROOT.h"
 
using namespace std;
 
class IBesGlobalTrigSvc;
class BesGlobalTrigSvc;
class TrigTOFT;
class TrigEACC;
class EmcWaveform;
 
class BesTrigL1 : public Algorithm {
public:
  /// Standard constructor
  BesTrigL1( const std::string& name, ISvcLocator* pSvcLocator );
  virtual ~BesTrigL1(){}; ///< Destructor
 
  virtual StatusCode initialize(); ///< Algorithm initialization
  virtual StatusCode execute();    ///< Algorithm execution
  virtual StatusCode finalize();   ///< Algorithm finalization
  void               findSETime( multimap<int, uint32_t, less<int>> mdc_hitmap,
                                 multimap<int, int, less<int>>      tof_hitmap,
                                 multimap<int, uint32_t, less<int>> emc_TC, double& stime, double& etime );
 
  void runAclock_mdc( int iclock, double stime,
                      multimap<int, uint32_t, less<int>> mdc_hitmap );
  void runAclock_tof( int iclock, double stime, int& idle_status,
                      std::multimap<int, int, less<int>> tof_hitmap );
  void runAclock_emc( int iclock, double stime, std::multimap<int, uint32_t, less<int>> emc_TC,
                      EmcWaveform* blockWave );
  void getEmcAnalogSig( EmcDigiCol* emcDigiCol, EmcWaveform ( &blockWave )[16],
                        multimap<int, uint32_t, less<int>>& emc_TC );
  void findEmcPeakTime( double& peak_time, EmcWaveform* blockWave );
 
  void stretchTrgCond( int nclock, int**& trgcond );
  void trgSAFDelay( int nclock, int**& trgcond );
  void trgGTLDelay( int nclock, int**& trgcond );
 
private:
  // define and initialize variables for output mdchit infor.
  double hit9[76], hit10[76], hit11[88], hit12[88];      // super layer 3
  double hit13[100], hit14[100], hit15[112], hit16[112]; // super layer 4
  double hit17[128], hit18[128], hit19[140], hit20[140]; // super layer 5
  double hit37[256], hit38[256], hit39[256], hit40[256]; // super layer 10
 
  BesGlobalTrigSvc*  m_pIBGT;
  IBesGlobalTrigSvc* m_tmpSvc;
 
  IRawDataProviderSvc* m_rawDataProviderSvc;
  IMdcGeomSvc*         m_MdcGeomSvc;
 
  // Reference to RealizationSvc
  RealizationSvc*    m_RealizationSvc;
  IEmcCalibConstSvc* emcCalibConstSvc;
 
  MdcTSF*      m_MdcTSF;
  TofHitCount* m_TofHitCount;
  EmcTCFinder* m_emcDigi;
  MucTrigHit*  m_mucDigi;
 
  ifstream    readin;
  ofstream    readout;
  ifstream    mdc_hit;
  std::string input;
  std::string output;
  std::string outEvtId;
  std::string indexfile;
  bool        mTrigRootFlag;
  int         m_runMode;
 
  bool ifpass;
  int  passNo;
  int  writeFile;
  int  ifoutEvtId;
 
  int totalEvent;
  int totalTracks;
 
  int clock_shift;
 
  NTuple::Tuple*      m_tuple;
  NTuple::Item<float> m_wire_x;
  NTuple::Item<float> m_wire_y;
  NTuple::Item<long>  m_wire_evtId;
  NTuple::Item<long>  m_delta_tdc;
 
  NTuple::Tuple*       m_tuple1;
  NTuple::Item<long>   m_RunId;
  NTuple::Item<long>   m_EventId;
  NTuple::Item<float>  m_mc_totE_all;
  NTuple::Item<float>  m_data_totE_all;
  NTuple::Item<float>  m_wetotE;
  NTuple::Item<float>  m_eetotE;
  NTuple::Item<long>   m_data_wetotE;
  NTuple::Item<long>   m_data_eetotE;
  NTuple::Item<long>   m_index_btc;
  NTuple::Array<float> m_btc_e;
  NTuple::Array<long>  m_data_btc;
  NTuple::Item<long>   m_cond_id;
  NTuple::Array<long>  m_mc_cond;
  NTuple::Array<long>  m_data_cond;
  NTuple::Item<long>   m_block_id;
  NTuple::Array<float> m_mc_blockE;
  NTuple::Array<float> m_data_blockE;
  NTuple::Array<float> m_R_blockE;
 
  NTuple::Tuple*      m_tuple2;
  NTuple::Item<long>  m_index2;
  NTuple::Array<long> m_fireLayer;
 
  NTuple::Item<long>  m_index3;
  NTuple::Array<long> m_hitLayer;
 
  NTuple::Item<long>  m_index4;
  NTuple::Array<long> m_hitSeg;
 
  NTuple::Item<long>   m_index5;
  NTuple::Array<float> m_costheta;
  NTuple::Array<float> m_phi;
  NTuple::Array<float> m_p;
  NTuple::Array<float> m_pdgcode;
 
  NTuple::Item<long>   m_index6;
  NTuple::Array<float> m_hitphi;
 
  NTuple::Item<long> m_nlayerEE;
  NTuple::Item<long> m_nlayerBR;
  NTuple::Item<long> m_nlayerWE;
  NTuple::Item<long> m_nlayerTotal;
 
  NTuple::Item<long> m_nhitEE;
  NTuple::Item<long> m_nhitWE;
  NTuple::Item<long> m_nhitBR;
  NTuple::Item<long> m_nhitTotal;
 
  NTuple::Item<float> m_mumcostheta;
  NTuple::Item<float> m_mumphi;
 
  NTuple::Item<long> m_trackfindall;
  NTuple::Item<long> m_trackfind3l;
  NTuple::Item<long> m_trackb;
  NTuple::Item<long> m_tracke;
  NTuple::Item<long> m_ntrack1;
  NTuple::Item<long> m_ntrack2;
  NTuple::Item<long> m_ntrack3;
  NTuple::Item<long> m_ngoodevent;
  NTuple::Item<long> m_ngoodtrack;
 
  NTuple::Tuple*     m_tuple3;
  NTuple::Item<long> m_evtId;
  NTuple::Item<long> m_condNOne[48];
  NTuple::Item<long> m_condNZero[48];
  int                nEvent;
 
  // get a handle on the Hist/TTree registration service
  ITHistSvc* m_thistsvc;
  TH1F*      m_trigCondi_MC;
  TH1F*      m_trigCondi_Data;
 
  TrigTOFT* toftrig;
  // TrigEACC* eacctrig;
};
#endif