AtPulseGADGET.h

Go to the documentation of this file.

/*****************************************************************/
/*    AtPulseTask: Simulates the ionized electrons which are */
/*    drifted and amplified by avalanche when reaching the mesh. */
/*    Log: Created in 24-10-2016                 */
/*    Author: Nathan Watwood (NSCL)              */
/*    ayyadlim@nscl.msu.edu                                      */
/*****************************************************************/
#ifndef AtPULSETASKGADGET_H
#define AtPULSETASKGADGET_H
 
#include "AtPulse.h"
#include "AtRawEvent.h" // for AtRawEvent
 
#include <Rtypes.h>
 
#include <memory> // for make_shared, shared_ptr
#include <vector>
 
class AtSimulatedPoint;
 
class AtPulseGADGET : public AtPulse {
 
private:
   Int_t AdjecentPads = 1; // Number of adjecent pads to be considered
   Int_t Items = 2 * AdjecentPads + 1;
   Double_t R = 10.0e6;        // Resistance (ohms)
   Double_t C = 2.8e-7;        // Capacitance[C] (Farad)
   Double_t t_amp = 1.32e-3;   // time amplitude (microseconds)
   Double_t Dc = 0.8673;       // Diffusion coefficent  (millimeter^2/microseconds)
   Double_t W = 2.2;           // widith (millimeter)
   Float_t SigmaPercent = 1.0; // Sigma percent for charge dispersion
 
   int fSkippy = 0; // keeps track Number of skipped charge disprsion points that have hit our veto pads
 
protected:
   Double_t ChargeDispersion(Double_t G, Double_t time, Double_t x0, Double_t y0, Double_t xi, Double_t yi);
   virtual bool AssignElectronsToPad(AtSimulatedPoint *point) override;
 
public:
   AtPulseGADGET(AtMapPtr map);
 
   ~AtPulseGADGET() = default;
   void SetSigmaPercent(Float_t sigma) { SigmaPercent = sigma; }
   void SetAdjecentPads(Int_t pads) { AdjecentPads = pads; };
 
   virtual AtRawEvent GenerateEvent(std::vector<AtSimulatedPoint *> &vec) override; 
   virtual std::shared_ptr<AtPulse> Clone() const override { return std::make_shared<AtPulseGADGET>(*this); }
};
 
#endif //#ifndef AtPULSETASKGADGET_H