In This Package:

QMLFTool.h

Go to the documentation of this file.

#ifndef QMLFTOOL_H
#define QMLFTOOL_H 1

#include "ReconAlg/IReconTool.h"
#include "GaudiAlg/GaudiTool.h"
#include "Event/CalibReadoutPmtCrate.h"
#include "Conventions/Electronics.h"
#include <map>

namespace DayaBay{
  class CalibReadout;
  class RecTrigger;
  //class RecTriggerExt;
  class CalibReadoutPmtCrate;
}

class IPmtGeomInfoSvc;
class ICableSvc;
class ICalibDataSvc;

class OpPara;
class GeomPara;

class QMLFTool : public GaudiTool , virtual public IReconTool
{
public:  
  
    QMLFTool(const std::string& type,
                   const std::string& name,
                   const IInterface* parent);
                 
    virtual ~QMLFTool();

    virtual StatusCode reconstruct(const DayaBay::CalibReadout&,
                                   DayaBay::RecTrigger&);
    virtual StatusCode initialize();
    virtual StatusCode finalize();

private:

    static void NCLL_fcn(int& npar, double* grad, double& fval, double* xval, int iflag);

    StatusCode getOpPara(std::string &data_src);
    StatusCode getGeomPara(std::string &data_src);
    StatusCode printPara();

    StatusCode setUsedChannels(const DayaBay::CalibReadoutPmtCrate*,
                               DayaBay::RecTrigger*);
    
    typedef std::map<DayaBay::FeeChannelId, bool> UsedFeeChannelMap;
    typedef std::map<DayaBay::FeeChannelId, double> ExpChgMap;

public:

    // Property CableSvcName - Name of the cable service
    std::string m_cableSvcName;

    //CableSvc
    static ICableSvc*  m_cableSvc;

    // Property PmtGeomSvcName - Name of the PmtGeom service
    std::string m_pmtGeomSvcName;

    // Pmt Geometry Information Service
    static IPmtGeomInfoSvc* m_pmtGeomSvc;

    // Property CalibDataSvcName - Name of the CalibData service
    std::string m_pmtCalibDataSvcName;

    // Pmt CalibData Information Service
    static ICalibDataSvc* m_pmtCalibDataSvc;

    static DayaBay::CalibReadoutPmtCrate* m_calibCrate;

    //DayaBay::RecTriggerExt *m_recTrigExt;
    static  ExpChgMap m_expChgMap;

    // different recon mode
    // 1 : not consider charge smearing
    // others : take charge smearing (S.P.E spectrum) into account
    static int    m_LLFMode; 

    // A flag that whether to add time likelihood or not
    // if 0: not to add; if non-zero: add. Default is 0.
    static bool    m_addTimeLLF;

    // If add time likelihood, the charge cut for PMTs selection criteria
    // It is useful when m_addTimeLLF=1
    static double    m_pmtChgCut;

    // A flag that whether to record extended information from reconstruction
    // such as expected charge pattern, default is no (0)
    static bool m_recExtInfo; 

    // Optical parameters for reconstruction
    static OpPara *m_opPara;
    std::string   m_opLocation;

    // Geometry parameters for reconstruction
    static GeomPara *m_geomPara;
    std::string   m_geomLocation;

    double m_ATTL;
    double m_TopRef;
    double m_BotRef;
    double m_LightYield;
    double m_PhotoCathodeArea;
    
    // Geometry parameters
    double m_TopRefZpos;
    double m_BotRefZpos;

    // Used channels
    static UsedFeeChannelMap m_timeLLFUsedChannels;
    static UsedFeeChannelMap m_chgLLFUsedChannels;
};

#endif // QMLFTOOL_H

---