QwRaster.h

Go to the documentation of this file.

/**********************************************************\
* File: QwRaster.h                                         *
*                                                          *
* Author: J. Pan                                           *
* jpan@jlab.org                                            *
*                                                          *
* Tue Jun 29 00:16:14 CDT 2010                             *
\**********************************************************/

#ifndef __QWRASTER__
#define __QWRASTER__

// System headers
#include <vector>

// ROOT headers
#include "TTree.h"
#include "TFile.h"
#include "TProfile2D.h"

// Qweak headers
#include "VQwSubsystemTracking.h"
#include "MQwV775TDC.h"
#include "QwScaler_Channel.h"
#include "QwVQWK_Channel.h"
#include "QwPMT_Channel.h"

// Forward declarations
class QwDatabase;

class QwRaster: public VQwSubsystemTracking, public MQwSubsystemCloneable<QwRaster>
{
 private:
  /// Private default constructor (not implemented, will throw linker error on use)
  QwRaster();

 public:
   /// Constructor with name
   QwRaster(const TString& name);
   /// Virtual destructor
   virtual ~QwRaster();

   // VQwSubsystem methods
   void ProcessOptions(QwOptions &options); //Handle command line options

   Int_t LoadEventCuts(TString filename) { return 0; };
   Bool_t ApplySingleEventCuts() { return kTRUE; };
   void PrintErrorCounters() const {};
   Bool_t CheckRunningAverages(Bool_t ) { return kTRUE; };

  /*  Member functions derived from VQwSubsystem. */
  Int_t LoadChannelMap(TString mapfile);
  Int_t GetDetectorInfo(std::vector< std::vector< QwDetectorInfo > > & detector_info){return 0;};
  Int_t LoadInputParameters(TString parameterfile);
  Int_t ProcessConfigurationBuffer(const UInt_t roc_id, const UInt_t bank_id, UInt_t* buffer, UInt_t num_words);
  void  InitializeChannel(TString name, TString datatosave);
  void  ClearEventData();

  Int_t ProcessEvBuffer(const UInt_t roc_id, const UInt_t bank_id, UInt_t* buffer, UInt_t num_words);
  void  ProcessEvent();

  using VQwSubsystem::ConstructHistograms;
  void  ConstructHistograms(TDirectory *folder, TString &prefix);
  void  FillHistograms();

  using VQwSubsystem::ConstructBranchAndVector;
  void  ConstructBranchAndVector(TTree *tree, TString &prefix, std::vector<Double_t> &values);
  void  FillTreeVector(std::vector<Double_t> &values) const;

  void  FillDB(QwDatabase *db, TString type) {return;};

  void  GetHitList(QwHitContainer & grandHitContainer){};
  void  ReportConfiguration();

  void PrintValue() { };
  void PrintInfo();

 protected:

  EQwModuleType fCurrentType;
  Bool_t fDEBUG;

  MQwV775TDC fQDCTDC;

  TString fRegion;  ///  Name of this subsystem (the region).

 protected:
  size_t fCurrentBankIndex;
  Int_t fCurrentSlot;
  Int_t fCurrentIndex;

  UInt_t fBankID[2];  //bank ID's of 2 different modules for QwRaster
                      //fBankID[0] for V792/V775 QDC_Bank
                      //fBankID[1] for SIS3801   SCA_Bank

 protected:
  static const UInt_t kMaxNumberOfModulesPerROC;
  static const UInt_t kMaxNumberOfChannelsPerModule;
  Int_t fNumberOfModules;

  std::vector< std::vector<Int_t> > fModuleIndex;  /// Module index, indexed by bank_index and slot_number

  std::vector< EQwModuleType > fModuleTypes;
  std::vector< std::vector< std::pair< EQwModuleType, Int_t> > > fModulePtrs; // Indexed by Module_index and Channel

  //    We need a mapping of module,channel into PMT index, ADC/TDC
  std::vector< std::vector<QwPMT_Channel> > fPMTs;  // for QDC/TDC and F1TDC

  std::vector<QwSIS3801D24_Channel> fSCAs;
  std::map<TString,size_t> fSCAs_map;
  std::vector<Int_t> fSCAs_offset;

  void FillRawWord(Int_t bank_index, Int_t slot_num, Int_t chan, UInt_t data);
  void  ClearAllBankRegistrations();
  Int_t RegisterROCNumber(const UInt_t roc_id);
  Int_t RegisterSubbank(const UInt_t bank_id);

  // Tells this object that it will decode data from the current bank
  Int_t RegisterSlotNumber(const UInt_t slot_id);
  EQwModuleType RegisterModuleType(TString moduletype);
  Int_t GetModuleIndex(size_t bank_index, size_t slot_num) const;
  Bool_t IsSlotRegistered(Int_t bank_index, Int_t slot_num) const {
    return (GetModuleIndex(bank_index,slot_num) != -1);
  };

  Int_t LinkChannelToSignal(const UInt_t chan, const TString &name);
  Int_t FindSignalIndex(const EQwModuleType modtype, const TString &name) const;

  //UInt_t GetEventcutErrorFlag() { return 0; };//return the error flag to the main routine

  // raster specified histograms
  TH2D* fRateMap;

 private:

  //static const Bool_t bStoreRawData;

  Double_t fPositionOffsetX;
  Double_t fPositionOffsetY;
  Double_t fChannel_Offset_X;
  Double_t fChannel_Offset_Y;
  Double_t fCal_Factor_QDC_X;
  Double_t fCal_Factor_QDC_Y;

  Double_t fPositionX_ADC;
  Double_t fPositionY_ADC;

  Double_t fbpm_3h07a_slope;
  Double_t fbpm_3h07a_intercept;
  Double_t fbpm_3h09b_slope;
  Double_t fbpm_3h09b_intercept;

  Double_t fbpm_3h07a_xp;
  Double_t fbpm_3h07a_xm;
  Double_t fbpm_3h07a_yp;
  Double_t fbpm_3h07a_ym;
  Double_t fbpm_3h09b_xp;
  Double_t fbpm_3h09b_xm;
  Double_t fbpm_3h09b_yp;
  Double_t fbpm_3h09b_ym;

  Double_t fbpm_3h07a_pos_x;
  Double_t fbpm_3h07a_pos_y;
  Double_t fbpm_3h09b_pos_x;
  Double_t fbpm_3h09b_pos_y;

};


#endif