phot::ScintTimeLAr Class Reference

#include <ScintTimeLAr.h>

Inheritance diagram for phot::ScintTimeLAr:

Public Member Functions
	ScintTimeLAr (fhicl::ParameterSet const &pset)
void	GenScintTime (bool is_fast, CLHEP::HepRandomEngine &engine)
Public Member Functions inherited from phot::ScintTime
	ScintTime ()
virtual	~ScintTime ()=default
double	GetScintTime () const

Private Member Functions
double	single_exp (double t, double tau2)
double	bi_exp (double t, double tau1, double tau2)

Private Attributes
int	LogLevel
double	SRTime
double	SDTime
double	FRTime
double	FDTime

Additional Inherited Members
Protected Attributes inherited from phot::ScintTime
double	timing

Detailed Description

Definition at line 27 of file ScintTimeLAr.h.

Constructor & Destructor Documentation

phot::ScintTimeLAr::ScintTimeLAr ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 14 of file ScintTimeLAr_tool.cc.

    : LogLevel{pset.get<int>("LogLevel")}
    , SRTime{pset.get<double>("SlowRisingTime", 0.0)}
    , SDTime{pset.get<double>("SlowDecayTime",  0.0)}
    , FRTime{pset.get<double>("FastRisingTime", 0.0)}
    , FDTime{pset.get<double>("FastDecayTime",  0.0)}
    {
        if ( LogLevel >= 1 ) 
        {
            std::cout << "ScintTimeLAr Tool configure:" << std::endl;
            std::cout << "Fast rising time: "   << FRTime
                      << ", Fast decay time: "  << FDTime
                      << ", Slow rising time: " << SRTime
                      << ", Slow decay time: "  << SDTime
                      << std::endl;
        }
    }

Member Function Documentation

double phot::ScintTimeLAr::bi_exp ( double  t, double  tau1, double  tau2  )

private

Definition at line 39 of file ScintTimeLAr_tool.cc.

    {
        return (((std::exp((-1.0 * t) / tau2) * (1.0 - std::exp((-1.0 * t) / tau1))) / tau2) / tau2) * (tau1 + tau2);
    }

void phot::ScintTimeLAr::GenScintTime ( bool  is_fast, CLHEP::HepRandomEngine &  engine  )

virtual

Implements phot::ScintTime.

Definition at line 49 of file ScintTimeLAr_tool.cc.

    {        
        double tau1;
        double tau2;
        
        if(is_fast == 1)  // fast scinitllation
        {
            tau1 = FRTime;
            tau2 = FDTime;
        }
        else
        {
            tau1 = SRTime;
            tau2 = SDTime;
        }
        
        CLHEP::RandFlat randflatscinttime{engine};
        
        if ((tau1 == 0.0) || (tau1 == -1.0))
        {
            timing = -tau2 * std::log(randflatscinttime());
            return;
        }
        
        //ran1, ran2 = random numbers for the algorithm        
        while (1)
        {
            auto ran1 = randflatscinttime();
            auto ran2 = randflatscinttime();
            auto d = (tau1 + tau2) / tau2;
            auto t = -tau2 * std::log(1 - ran1);
            auto g = d * single_exp(t, tau2);
            if (ran2 <= bi_exp(t, tau1, tau2) / g)
            {
                timing = t;
                return;
            }
        }
    }

double phot::ScintTimeLAr::single_exp ( double  t, double  tau2  )

private

Definition at line 33 of file ScintTimeLAr_tool.cc.

    {
        return std::exp((-1.0 * t) / tau2) / tau2;
    }

Member Data Documentation

double phot::ScintTimeLAr::FDTime

private

Definition at line 40 of file ScintTimeLAr.h.

double phot::ScintTimeLAr::FRTime

private

Definition at line 39 of file ScintTimeLAr.h.

int phot::ScintTimeLAr::LogLevel

private

Definition at line 34 of file ScintTimeLAr.h.

double phot::ScintTimeLAr::SDTime

private

Definition at line 38 of file ScintTimeLAr.h.

double phot::ScintTimeLAr::SRTime

private

Definition at line 37 of file ScintTimeLAr.h.

---

The documentation for this class was generated from the following files:

• larsim/larsim/PhotonPropagation/ScintTimeTools/ScintTimeLAr.h
• larsim/larsim/PhotonPropagation/ScintTimeTools/ScintTimeLAr_tool.cc