test_recombinationmodels.cxx File Reference

#include "WireCellGen/RecombinationModels.h"
#include "WireCellUtil/Units.h"
#include "WireCellUtil/Testing.h"
#include <vector>
#include <iostream>

Go to the source code of this file.

Functions
void	spin (IRecombinationModel &model, const std::string &label)
void	test_larsoft ()
void	test_mip ()
void	test_birks ()
void	test_box ()
int	main ()

Variables
const std::vector< double >	evec {0.0797872,0.0533235,0.000891013, 0.02}
const std::vector< double >	svec {0.0208263, 0.0311329, 0.000947381, 0.01}
const std::vector< int >	nvec {3380, 2259, 37, 847}

Function Documentation

int main

(

void

)

Definition at line 69 of file test_recombinationmodels.cxx.

{
    test_larsoft();
    test_mip();
    test_birks();
    test_box();
    return 0;
}

void spin	(	IRecombinationModel &	model,
		const std::string &	label
	)

Definition at line 29 of file test_recombinationmodels.cxx.

{
    cerr << "Model: " << label << endl;
    for (size_t ind=0; ind<nvec.size(); ++ind) {
        // values from larsoft dump
        double lsenergy = evec[ind]*units::MeV;
        double lsstep = svec[ind]*units::cm;
        double lscharge = nvec[ind]*(-1.0*units::eplus);
        double dEdX = lsenergy/lsstep;

        double gotcharge = model(lsenergy, lsstep);
        cerr << "\tdE=" << lsenergy/units::keV << " keV, "
             << "dX=" << lsstep/units::micrometer << " um, "
             << "dE/dX=" << dEdX/(units::MeV/units::cm) << " MeV/cm, "
             << "dQls=" << lscharge/(-1*units::eplus) << " ele, "
             << "dQ=" << gotcharge/(-1*units::eplus) << " ele\n";
    }

}

void test_birks

(

)

Definition at line 58 of file test_recombinationmodels.cxx.

{
    Gen::BirksRecombination model;
    spin(model, "birks");
}

void test_box

(

)

Definition at line 63 of file test_recombinationmodels.cxx.

{
    Gen::BoxRecombination model;
    spin(model, "box");
}

void test_larsoft

(

)

Definition at line 48 of file test_recombinationmodels.cxx.

{
    Gen::MipRecombination model(1.0);
    spin(model, "larsoft");
}

void test_mip

(

)

Definition at line 53 of file test_recombinationmodels.cxx.

{
    Gen::MipRecombination model;
    spin(model, "mip");
}

Variable Documentation

const std::vector<double> evec {0.0797872,0.0533235,0.000891013, 0.02}

Definition at line 25 of file test_recombinationmodels.cxx.

const std::vector<int> nvec {3380, 2259, 37, 847}

Definition at line 27 of file test_recombinationmodels.cxx.

const std::vector<double> svec {0.0208263, 0.0311329, 0.000947381, 0.01}

Definition at line 26 of file test_recombinationmodels.cxx.