triplePi0.cc File Reference

#include "triplePi0.h"
#include "expconfig/ExpConfig.h"
#include "utils/Combinatorics.h"
#include "utils/ParticleTools.h"
#include "base/vec/LorentzVec.h"
#include "base/Logger.h"
#include "utils/uncertainties/Interpolated.h"
#include "analysis/physics/Plotter.h"
#include "plot/CutTree.h"
#include "slowcontrol/SlowControlVariables.h"

Include dependency graph for triplePi0.cc:

Classes
class	triplePi0_Test

Typedefs
using	triplePi0_PlotBase = TreePlotterBase_t< triplePi0::PionProdTree >

Functions
triplePi0::fitRatings_t	applyTreeFit (utils::TreeFitter &fitter, const std::vector< utils::TreeFitter::tree_t > &intermediates, const utils::ProtonPhotonCombs::comb_t &protonSelection, const double Ebeam)

Variables
auto	getEgamma
auto	reducedChi2
auto	getTreeFitPhotonIndices
auto	getLorentzSumFitted

Typedef Documentation

triplePi0_PlotBase

using triplePi0_PlotBase = TreePlotterBase_t<triplePi0::PionProdTree>

Function Documentation

applyTreeFit()

triplePi0::fitRatings_t applyTreeFit	(	utils::TreeFitter &	fitter,
		const std::vector< utils::TreeFitter::tree_t > &	intermediates,
		const utils::ProtonPhotonCombs::comb_t &	protonSelection,
		const double	Ebeam
	)

Variable Documentation

getEgamma

auto getEgamma

Initial value:

= [] (const TParticleTree_t& tree)
{
    return tree->Get()->Ek();
}

getLorentzSumFitted

auto getLorentzSumFitted

Initial value:

= [](const vector<utils::TreeFitter::tree_t>& nodes)
{
    vector<TLorentzVector> acc;
    for ( const auto& node: nodes)
    {
        acc.push_back(node->Get().LVSum);
    }
    return acc;
}

getTreeFitPhotonIndices

auto getTreeFitPhotonIndices

Initial value:

= [] (const TParticleList& orig_Photons,
                                   const utils::TreeFitter& treeFitter)
{
    const auto allP = treeFitter.GetFitParticles();
    const vector<utils::Fitter::FitParticle> fitPhotons(allP.begin()+1,
                                                        allP.end());

    vector<unsigned> combination;
    for (unsigned iFit = 0; iFit < fitPhotons.size(); ++iFit)
    {
        for (unsigned iOrig = 0; iOrig < orig_Photons.size(); ++iOrig)
        {
            if ( fitPhotons.at(iFit).Particle == orig_Photons.at(iOrig))
            {
                combination.push_back(iOrig);
                continue;
            }
        }
    }
    return combination;
}

reducedChi2

auto reducedChi2

Initial value:

= [](const APLCON::Result_t& ares)
{
    return 1. * ares.ChiSquare / ares.NDoF;
}