In This Package:

EsPmtEffectPulseTool.h

Go to the documentation of this file.

#ifndef ESPMTEFFECTPULSETOOL_H
#define ESPMTEFFECTPULSETOOL_H 1

#include "ElecSim/IEsPulseTool.h"
#include "Event/ElecPmtPulse.h"
#include "GaudiAlg/GaudiTool.h"
#include "GaudiKernel/IRndmGenSvc.h"
#include "GaudiKernel/RndmGenerators.h"
#include "DataSvc/ICableSvc.h"
#include "DataSvc/ISimDataSvc.h"
#include "Context/TimeStamp.h"

namespace DayaBay{
  class SimHitCollection;
  class ElecPulseCollection;
}

class EsPmtEffectPulseTool : public GaudiTool , virtual public IEsPulseTool
{
public:  
  
    EsPmtEffectPulseTool(const std::string& type,
                     const std::string& name,
                     const IInterface* parent);
                 
    virtual ~EsPmtEffectPulseTool();
    
    virtual StatusCode generatePulses(DayaBay::SimHitCollection*,
                                      DayaBay::ElecPulseCollection*);
    virtual StatusCode initialize();
    virtual StatusCode finalize();

 private:

    // Property CableSvcName - Name of the cable service (default="StaticCableSvc")
    std::string m_cableSvcName;    
    /* Property SimDataSvcName - Name of the simulation data service
     * that provides access to pmt properties (default="StaticSimDataSvc")
     */
    std::string m_simDataSvcName;    

    // user-defined PDF and associated bin edges for pre- and after-pulses
    // PDFs defined wrt main pulse time
    // default is single count in [-50ns,0ns] and [0us,10us], respectively
    std::vector<double> m_prePulsePdf;
    std::vector<double> m_prePulseEdges;
    std::vector<double> m_afterPulsePdf;
    std::vector<double> m_afterPulseEdges;
    std::vector<double> m_afterPulseTime;
    
    bool m_enablePrePulse;
    bool m_enableAfterPulse;
    
    // Mode for afterpulse amplitude distribution:
    // SinglePE: Fixed single p.e. amplitude with Gaussian smearing
    // PDF: Use amplitude PDF with peak at single p.e and long tail
    std::string m_afterPulseAmpMode;
    std::string m_darkPulseAmpMode;
    
    // PDFs for afterpulse amplitude
    std::vector<double> m_afterPulseAmpPdf;
    std::vector<double> m_afterPulseAmpEdges;

    // Property: EnableNonlinearAfterpulse
    // Applies nonlinear afterpulsing rate model if true
    bool m_enableNonlinearAfterpulse;
    // Property: LinearAfterpulseThreshold
    // Upper limit of linear afterpulsing in number of PE.
    double m_linearAfterpulseThreshold;

    // Property: EnableVeryLongTimeSuppression
    // Enable suppression of hit times in the far future
    bool m_enableVeryLongTimeSuppression;
    // Property: VeryLongTime
    // Definition of very long time for hit time suppression
    double m_veryLongTime; 

    //Weight of the exponential contribution to the spe response function
    double m_expWeight;
    // Decay time of the exponential
    double m_speExpDecay;
    double m_speCutoff;

    // CableSvc
    ICableSvc *m_cableSvc;
    // SimDataSvc
    ISimDataSvc *m_simDataSvc;
    
    // Random numbers
    Rndm::Numbers m_uni,m_gauss,m_exp;
    Rndm::Numbers m_randPrePulseTime,m_randAfterPulseTime;

    // electronic simulation beginning and end times, defined by ElecHeader
    TimeStamp m_simTimeEarliest;
    TimeStamp m_simTimeLatest;

    // DWL: Size of PMT hits counting window
    double m_timeInterval;

    DayaBay::ElecPulse* makePrePulse  (DayaBay::ElecPulse* simpulse,
                                       const DayaBay::PmtSimData& pmtData);
    DayaBay::ElecPulse* makeAfterPulse(DayaBay::ElecPulse* simpulse,
                                       const DayaBay::PmtSimData& pmtData);
    DayaBay::ElecPulse* makeDarkPulse (const DayaBay::ElecChannelId& channelId,
                                       const DayaBay::PmtSimData& pmtData);
    float  PulseAmp(float weight, float gain, float gainFwhm);
    double ConvertPdfRand(double rand, std::vector<double> pdf, std::vector<double> edges);
    double ConvertPdfRand01(double rand, std::vector<double> pdf, std::vector<double> edges);
    int PoissonRand(double mean);
    // DWL: computing number of after-pulses based on the charge of the primary pulse
    double NumAfterPulse(const int numPmtHit);
    // Get afterpulse amplitude PDF
    void getAfterPulseAmpPdf(std::vector<double>& pdf);
    void getAfterPulseAmpEdges(std::vector<double>& edges);
};

#endif // ESPMTEFFECTPULSETOOL_H

---