In This Package:

PreElecSimSvc.h

Go to the documentation of this file.

#ifndef PREELECSIMALG_H
#define PREELECSIMALG_H

#include "MixSimHit.h"

#include "Event/GenHeader.h"
#include "Event/SimHeader.h"
#include "Event/SimHit.h"

#include "Conventions/Electronics.h"

#include "MixInputSvc/IMixInputSvc.h"
#include "PreElecSimSvc/IPreElecSimSvc.h"

#include "GaudiKernel/Service.h"
#include "GaudiKernel/MsgStream.h"

#include <map>
#include <deque>
#include <vector>

using namespace std;

class PreElecSimSvc : public Service, 
                      virtual public IPreElecSimSvc
{
 public:
  PreElecSimSvc(const string& name, ISvcLocator* pSvcLocator);
  virtual ~PreElecSimSvc();

  virtual StatusCode initialize();
  virtual StatusCode queryInterface(const InterfaceID& riid, void** ppint);

  SimHitHeader* getSimHitHeader(vector<const DayaBay::SimHeader*>& shs, TimeStamp& hite, TimeStamp& hitl);

 private:
  //MsgStream
  mutable MsgStream         m_log;

  //Service for multi-events inputting
  IMixInputSvc*             m_mixInputSvc;

  //Parameters to dealing with SimHits
  //vector<double>    m_hitsgap;
  deque<TimeStamp>  m_hitsGap;     //for SimHeader : absolute time
  deque<MixSimHit>  m_tranSimHits; //for merged SimHits buffer
  deque<MixSimHit>  m_SimHits;     //real SimHits buffer

  //Fill Hits buffer and HitHeader with given SimHeader's
  StatusCode fillHitsBuffer(DayaBay::SimHeader* sh);

  //Search Hits buffer for a gap, update the gap vector
  StatusCode findGap();

  //Filling one or more SimHitHeader with proper SimHits
  //then push_back the SimHitHeader into m_propSimHitHeaders
  //Bookkeep references between SimHitHeader and SimHeaders
  StatusCode createHitHeaders();

  //SimHitHeader buffer to hold current SimHitHeader
  //This two containers' size must be the same!
  map<DayaBay::SimHitHeader*, vector<const DayaBay::SimHeader*> > m_HitHeaderToHeaders;
  deque<DayaBay::SimHitHeader*> m_propSimHitHeaders;

  //SimHit chunks' earliest and latest
  deque<TimeStamp> m_hitEarliests;
  deque<TimeStamp> m_hitLatests;
    
  //Minimum width of time gap that determines a chunk.
  //double    m_mintimeGap;
  TimeStamp m_minTimeGap;
  TimeStamp m_absTime;    //current SimHeader's timeStamp

  //Decide the ouput SimHit chunk's site
  bool                            m_ifDeal;
  deque<TimeStamp>::size_type     m_dealInd;
};

#endif //PREELECSIMALG_H

---