MdcPerformance.h

Go to the documentation of this file.

#ifndef MDCPERFORMANCE_H
#define MDCPERFORMANCE_H
#include "GaudiKernel/Algorithm.h"
#include <string>
#include <vector>
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/HepPDT/ParticleDataTable.hh"
#include "GaudiKernel/NTuple.h"
#include "McTruth/McParticle.h"
#include "McTruth/MdcMcHit.h"
#include "RawDataProviderSvc/RawDataProviderSvc.h"
#include "MdcGeomSvc/IMdcGeomSvc.h"
#include "MdcGeomSvc/MdcGeoWire.h"
#include "MdcGeomSvc/MdcGeoLayer.h"
#include "MagneticField/IMagneticFieldSvc.h"
#include "BField/BField.h"
#include "MdcNavigation/MdcCheckUtil.h"
class MdcDetector;

class MdcPerformance:public Algorithm {
public:
  MdcPerformance (const std::string& name, ISvcLocator* pSvcLocator);
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();

private:
  StatusCode bookNTuple();
  StatusCode fillInit();
  StatusCode fillRawEvent();
  StatusCode getNoiseThroughRawDigi();
  StatusCode fillRecEvent();
  StatusCode fillRecTrack(const RecMdcTrack* tk);

  HepPDT::ParticleDataTable* m_particleTable;
  RawDataProviderSvc*        m_rawDataProviderSvc; 
  IMagneticFieldSvc*         m_pIMF;
  const MdcDetector *        m_gm;
  double Bz;

  int    havedigi[43][288]; 
  short  multiTdcCount[43][288];
  bool   m_doSag;
  bool   m_rawData;
  long   t_eventNo;
  long   t_runNo;
  double t_t0;
  double t_t0Stat;
  double t_t0Truth;
  int    t_timing;
  int    t_nTk;
  double t_recTkNum;
  int    t_mcTkNum;
  int    t_patRecTk ;
  int    t_mdcxRecTk;
  int    t_trkRecoTk;

  bool              m_keepFirstTdc; 
  bool              m_poca;
  int               m_maxMdcDigi;
  bool              m_keepBadTdc;
  bool              m_dropHot;
  bool              m_keepUnmatch;
  double            m_mass;
  MdcCheckUtil*     m_mdcCheckUtil;

  NTuple::Tuple*  m_tupleRawEvt;
  NTuple::Tuple*  m_tuple1;
  NTuple::Item<double>             m_t3_recTk;
  NTuple::Item<long>               m_t3_TrkReco;
  NTuple::Item<long>               m_t3_PatRec;
  NTuple::Item<double>             m_t3_evtNo;
  NTuple::Item<long>               m_t3_runNo;
  NTuple::Item<double>             m_t3_t0; 
  NTuple::Item<double>             m_t3_t0Truth; 
  NTuple::Item<long>               m_t3_t0Stat; 
  //trigger
  NTuple::Item<long>               m_trig_index;
  NTuple::Array<long>              m_trig_cond;
  NTuple::Array<long>              m_trig_chan;
  NTuple::Item<long>               m_trig_timewindow;
  NTuple::Item<long>               m_trig_timetype;

  NTuple::Item<long>               m_t3_nDigi;
  NTuple::Item<long>               m_t3_nDigiUnmatch;
  NTuple::Array<long>              m_t3_layer;
  NTuple::Array<long>              m_t3_wire;
  NTuple::Array<long>              m_t3_gwire;
  NTuple::Array<double>            m_t3_rt;
  NTuple::Array<double>            m_t3_rc;
  NTuple::Array<long>              m_t3_ovfl;

  NTuple::Item<double>             m_evtNo;
  NTuple::Item<double>             m_runNo;
  NTuple::Item<double>             m_d0;
  NTuple::Item<double>             m_phi0;
  NTuple::Item<double>             m_cpa;
  NTuple::Item<double>             m_z0;
  NTuple::Item<double>             m_tanl;
  NTuple::Item<double>             m_d0E;
  NTuple::Item<double>             m_phi0E;
  NTuple::Item<double>             m_cpaE;
  NTuple::Item<double>             m_z0E;
  NTuple::Item<double>             m_tanlE;
  NTuple::Item<double>             m_q;
  NTuple::Item<double>             m_recTkNum;
  NTuple::Item<double>             m_p;
  NTuple::Item<double>             m_pt;
  NTuple::Item<double>             m_pz;
  NTuple::Item<double>             m_pcms;

  NTuple::Item<long>               m_nHit;
  NTuple::Item<long>               m_nAct;
  NTuple::Item<long>               m_nSt;
  NTuple::Item<long>               m_nDof;
  NTuple::Item<double>             m_fiTerm;
  NTuple::Item<double>             m_chi2;
  NTuple::Item<double>             m_chi2Dof;
  NTuple::Item<double>             m_chi2Prob;
  NTuple::Item<long>               m_tkStat;


  NTuple::Item<long>               m_tkId;
  NTuple::Array<long>              m_layer;
  NTuple::Array<long>              m_wire;
  NTuple::Array<double>            m_gwire;
  NTuple::Array<double>            m_resid;
  NTuple::Array<double>            m_driftD;
  NTuple::Array<double>            m_driftT;
  NTuple::Array<double>            m_doca;
  NTuple::Array<double>            m_entra;
  NTuple::Array<double>            m_zhit;
  NTuple::Array<double>            m_chi2add;
  NTuple::Array<double>            m_flaglr;
  NTuple::Array<double>            m_Tdc;
  NTuple::Array<double>            m_Adc;
  NTuple::Array<long>            m_hitStat;


};

#endif

---