DAS::Prefiring::Prefiring Namespace Reference

Functions
double	ansatz (double *x, double *p)
unique_ptr< TH1 >	Copy (const unique_ptr< TH1 > &h, double maxPt)

Function Documentation

ansatz()

double DAS::Prefiring::Prefiring::ansatz	(	double *	x,
		double *	p
	)

NOTE: this analytical is partly arbitrary.

{
    return 0.5*p[0]*(1+erf((log(*x)-log(p[1]))/(p[2])));
}

Copy()

unique_ptr<TH1> DAS::Prefiring::Prefiring::Copy	(	const unique_ptr< TH1 > &	h,
		double	maxPt
	)

{
    vector<double> edges;
    for (int i = 1; i <= h->GetNbinsX()+1; ++i) {
        auto edge = h->GetBinLowEdge(i);
        auto content = h->GetBinContent(i);
        cout << i << ' ' << edge << ' ' << content << '\n';
        edges.push_back(edge);
    }
    cout << flush;
    edges.back() = maxPt;
    cout << edges.back() << endl;
    auto g = make_unique<TH1D>("g", "", edges.size()-1, edges.data());
    for (int i = 1; i <= g->GetNbinsX(); ++i) {
        double content = h->GetBinContent(i),
               error   = h->GetBinError  (i);
        g->SetBinContent(i, content);
        //error = max(error, 0.2*content);
        g->SetBinError  (i, error  );
    }
    return g;
}