In This Package:

EsIdealFeeTool.h

Go to the documentation of this file.

#ifndef ESIDEALFEETOOL_H
#define ESIDEALFEETOOL_H 1

#include "ElecSim/IEsFrontEndTool.h"
#include "GaudiAlg/GaudiTool.h"
#include "Conventions/Electronics.h"
#include "Event/ElecPulseCollection.h"
#include "GaudiKernel/IRndmGenSvc.h"
#include "GaudiKernel/RndmGenerators.h"
#include "DataSvc/ICableSvc.h"
#include "DataSvc/ISimDataSvc.h"

namespace DayaBay{
  class ElecCrate;
  class ElecFeeCrate;
  class ElecFeeChannel;
  class ElecPmtPulse;
}

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

    virtual StatusCode generateSignals(DayaBay::ElecPulseCollection*,
                                       DayaBay::ElecCrate*);
    virtual StatusCode initialize();
    virtual StatusCode finalize();

 private:
    StatusCode loadResponse();
    StatusCode mapPulsesByChannel(
                    const DayaBay::ElecPulseCollection::PulseContainer& pulses,
                    std::map<DayaBay::ElecChannelId,
                    std::vector<DayaBay::ElecPmtPulse*> >& pulseMap,
                    int* hitCountSlices);
    StatusCode updateBoardSignals(DayaBay::ElecFeeCrate* crate);
    StatusCode updateEsumSignals(DayaBay::ElecFeeCrate* crate);

 private:
    virtual StatusCode generateOneChannel(const DayaBay::FeeChannelId&,
                                          std::vector<DayaBay::ElecPmtPulse*>&,
                                          DayaBay::ElecFeeChannel&,
                                          double simTime,
                                          const ServiceMode&);
    int convertClock(int clock, int inputFrequency, int outputFrequency);
    DayaBay::ESumComp::ESumComp_t boardEsumType(int boardId);
    void convolve(const DayaBay::AnalogSignal& signal,
                  const DayaBay::AnalogSignal& response,
                  DayaBay::AnalogSignal& output); 
    void shape(const DayaBay::AnalogSignal& rawSignal,
                  DayaBay::AnalogSignal& shapedSignal, 
                  const DayaBay::AnalogSignal& response, 
                  std::map<int,int> convolutionWindows); 
    void digitize(const DayaBay::AnalogSignal& signal,
                  DayaBay::DigitalSignal& output,
                  double range, int bits, double offset);
    void discriminate(const DayaBay::AnalogSignal& signal,
                      double threshold,
                      std::vector<int>& output,
                      DayaBay::Threshold::Threshold_t mode);
    double pmtPulse(double deltaT, int nPulse, bool addOvershoot);
    double overshoot(double deltaT, int nPulse);
    double ringing(double deltaT, double pulseMax);
    double satFactor(double amp);
    double pulseShaping(double deltaT);
    double esumShaping(double deltaT);
    double getPulseWidth(int nPulse);
    double getPulseForm(int nPulse);
    double getPulseShift(int nPulse);
    double getOvershootAmp(int nPulse);
    double noise();
    void sample(const DayaBay::AnalogSignal& signal,
                DayaBay::AnalogSignal& output,
                int inputFrequency, int outputFrequency);
 private:
    // Property CableSvcName - Name of the cable service
    std::string m_cableSvcName;
    // Property CableSvcName - Service name for electronics sim properties
    std::string m_simDataSvcName;
    // CableSvc
    ICableSvc *m_cableSvc;
    // CableSvc
    ISimDataSvc *m_simDataSvc;
    // Trigger Window Cycles
    int m_triggerWindowCycles;
    // Mean PMT gain (in units of 10^7)
    double m_gainFactor;
    // Simulation Frequency (in Hz)
    int m_simFrequency;
    // Esum Frequency (in Hz)
    int m_eSumFrequency;
    // Bool to turn on (true) or off (false) all nonlinear effects
    bool m_enableNonlinearity;
    // Bool to turn on (true) or off (false) high pass filter before discriminator
    bool m_enableAcCoupling;
    // Maximum number of pulses before nonlinear effects begin
    unsigned int m_linearityThreshold;
    
    // Bool to turn on (true) or off (false) pretrigger
    bool m_enablePretrigger;
    // Thresholds as set in TrigSim
    double m_eSumTriggerThreshold;
    int m_nHitTriggerThreshold;
    // Pretrigger threshold
    int m_pretriggerThreshold;
    // Bool to individually turn on (true) or off (false) noise
    bool m_enableNoise;
    // Bool to individually turn on (true) or off (false) dynamic single p.e. waveform model
    bool m_enableDynamicWaveform;
    // Bool to individually turn on (true) or off (false) saturation
    bool m_enableSaturation;
    // Bool to individually turn on (true) or off (false) ringing
    bool m_enableRinging;
    // Bool to individually turn on (true) or off (false) overshoot
    bool m_enableOvershoot;
    // Bool to turn on( true) or off (false) simulation of ESum waveforms
    bool m_enableESumTotal;
    bool m_enableESumH;
    bool m_enableESumL;
    // Bool to turn on (true) or off (false) fast simulation mode
    bool m_enableFastSimMode;
    // Pulse counting time slot width.
    double m_pulseCountSlotWidth;
    // Number of counted slots before and after the pulse.
    int m_pulseCountWindow;
    // Voltage amplitude of noise.
    double m_noiseAmp;
    // Spe pulse height in V at 1e7 gain
    double m_speAmp;
    // Factor to scale all channels' discriminator threshold
    double m_discThreshScale;
    
    double m_adcRange;
    double m_adcBits;
    double m_adcOffset;
    // Ideal PMT pulse shape.
    DayaBay::AnalogSignal m_pmtPulse;
    // Ideal PMT pulse shape without overshoot
    DayaBay::AnalogSignal m_pmtPulse_noOvershoot;
    // Ideal CR-(RC)^4 shaped PMT pulse including overshoot.
    DayaBay::AnalogSignal m_shapedPmtPulse;
    // Ideal CR-(RC)^4 shaped PMT pulse with no overshoot.
    DayaBay::AnalogSignal m_shapedPmtPulse_noOvershoot;
    // Single P.E. Ringing shape.
    DayaBay::AnalogSignal m_ringing;
    // Single P.E. Ringing shape, after CR-(RC)^4 shaping
    DayaBay::AnalogSignal m_shapedRinging;
    //ESum response after RC shaping
    DayaBay::AnalogSignal m_esumResponse;
    // CR-(RC)^4 shaping response.
    DayaBay::AnalogSignal m_shapingResponse;
    // Area of ideal CR-(RC)^4 shaped PMT pulse.
    double m_shapingSum;
    //Area of shaped ESum pulse
    double m_esumShapingSum;
    // PMT pulse peak amplitude.
    double m_pulseMax;
    
    // Random numbers.
    Rndm::Numbers m_gauss;
    // Transient buffers for signal processing
    DayaBay::AnalogSignal m_rawSignal;
    DayaBay::AnalogSignal m_esumL;
    DayaBay::AnalogSignal m_esumH;
    DayaBay::AnalogSignal m_esumTotal;
    DayaBay::DigitalSignal m_shapedEsumADC;
    DayaBay::AnalogSignal m_shapedEsumH;
    DayaBay::AnalogSignal m_shapedEsumL;
    DayaBay::AnalogSignal m_shapedEsumTotal;
    DayaBay::AnalogSignal m_discriminatorSignal;
    DayaBay::AnalogSignal m_shapedSignal;
    DayaBay::AnalogSignal m_tdcSignal;
    DayaBay::AnalogSignal m_adcAnalogSignal;
    DayaBay::DigitalSignal m_hitSignal;
    DayaBay::DigitalSignal m_adcHighGain;
    DayaBay::DigitalSignal m_adcLowGain;
    DayaBay::AnalogSignal m_energySignal;
    DayaBay::DigitalSignal m_hitSync;
    DayaBay::DigitalSignal m_hitHold;
    std::map<int,int> m_adcConvolutionWindows;
    std::map<int,int> m_esumConvolutionWindows;
};

#endif // ESIDEALFEETOOL_H

---