In This Package:

FlasherIdTool.h

Go to the documentation of this file.

/*
 * FlasherIdTool
 *
 * Identify PMT Flashing Events in the AD
 *
 * dandwyer@caltech.edu 2010/10/23
 */


#ifndef FLASHERIDTOOL_H
#define FLASHERIDTOOL_H


// Fast index for histograms
#define FLASHERQVST 0
#define FLASHERRINGVSCOLUMN 1
#define NHISTOGRAMS 2

#define NBOARDS 21
#define NCONNECTORS 17

#define MAXHISTS NBOARDS*NCONNECTORS*NHISTOGRAMS

#define MAXDT 4000

#include "GaudiAlg/GaudiTool.h"
#include "Context/TimeStamp.h"
#include "Context/Context.h"
#include "Conventions/Detectors.h"
#include "ProcessTools/IReadoutProcessor.h"
#include <vector>
#include <map>
#include <string>

class IStatisticsSvc;
class ICableSvc;
class IPmtGeomInfoSvc;
class TH1;

class FlasherIdTool : public GaudiTool,
                      virtual public IReadoutProcessor
{

 public:
    FlasherIdTool(const std::string& type,
                  const std::string& name, 
                  const IInterface* parent);
    virtual ~FlasherIdTool();

    // GaudiTool interface
    virtual StatusCode initialize();
    virtual StatusCode finalize();

    // IReadoutTool interface
    virtual StatusCode process(const DayaBay::ReadoutHeader* readoutHeader);

 private:
    TH1* getOrMakeHist(int run, const DayaBay::Detector& detector,
                       int board, int connector, int histogram);
    std::string getPath(int run, const char* detectorName,
                        int board, int connector, 
                        const char* histName);
    StatusCode initializeDistanceTable(const Context& context);
    unsigned int channelIndex(unsigned int board, unsigned int connector);
    double tdcOffset(int board, int connector);
    double transitTime(unsigned int board);

    // Properties
    double m_flasherThreshold;
    std::string m_detectorName;

    // Handle for statistics service
    IStatisticsSvc* m_statsSvc;
    ICableSvc* m_cableSvc;
    IPmtGeomInfoSvc* m_pmtInfoSvc;
    // Internal variables
    TimeStamp* m_firstTriggerTime;
    TimeStamp m_lastTriggerTime;
    std::map<int,TH1**> m_shortCuts;

    DayaBay::Detector m_detector;
    unsigned int m_nEvents;
    unsigned int m_maxTdcCalibEvents;
    double m_boardTdcOffset[NBOARDS];
    unsigned int m_boardNTdcs[NBOARDS];
    double m_topCalibTdcOffset;
    unsigned int m_topCalibNTdcs;
    double m_bottomCalibTdcOffset;
    unsigned int m_bottomCalibNTdcs;
    unsigned int m_flashBoard[NBOARDS*NCONNECTORS];
    unsigned int m_flashConnector[NBOARDS*NCONNECTORS];
    unsigned int m_flashRing[NBOARDS*NCONNECTORS];
    unsigned int m_flashColumn[NBOARDS*NCONNECTORS];
    float m_flashDeltaAdc[NBOARDS*NCONNECTORS];
    double m_flashTdc[NBOARDS*NCONNECTORS];
    double m_flashLikelihood[NBOARDS*NCONNECTORS];
    double m_adcSum;
    double m_dt;
    unsigned int m_nChannels;
    double m_minDt;
    double m_logLikelihood[MAXDT];
    double m_distance[NBOARDS*NCONNECTORS][NBOARDS*NCONNECTORS];
    unsigned int m_ring[NBOARDS*NCONNECTORS];
    unsigned int m_column[NBOARDS*NCONNECTORS];
};

#endif  // FLASHERIDTOOL_H

---