DAS::Photon Namespace Reference

Classes
class	ConversionVetoApplier
class	IDApplier

Functions
void	applyPhotonConversionVeto (const vector< fs::path > &inputs, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})
void	applyPhotonID (const vector< fs::path > &inputs, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})

Function Documentation

applyPhotonConversionVeto()

void DAS::Photon::applyPhotonConversionVeto	(	const vector< fs::path > &	inputs,
		const fs::path &	output,
		const pt::ptree &	config,
		const int	steering,
		const DT::Slice	slice = {1,0}
	)

Applies one of the photon conversion veto and (optionally) scale factors.

This function fetches the scale factor for the ID cut previously applied using applyPhotonID (automatically retrieved from the metainfo).

Todo:

ControlPlots

Parameters

inputs	input ROOT files (n-tuples)
output	output ROOT file (n-tuple)
config	config handled with `Darwin::Tools::options`
steering	parameters obtained from explicit options
slice	number and index of slice

                              {1,0} 
        )
{
    cout << __func__ << ' ' << slice << " start" << endl;
 
    unique_ptr<TChain> tIn = DT::GetChain(inputs);
    unique_ptr<TFile> fOut(DT_GetOutput(output));
    auto tOut = unique_ptr<TTree>(tIn->CloneTree(0));
 
    DT::MetaInfo metainfo(tOut);
    metainfo.Check(config);
 
    vector<RecPhoton> * recPhotons = nullptr;
    tIn->SetBranchAddress("recPhotons", &recPhotons);
 
    auto ID = metainfo.Get<string>("corrections", "photons", "ID", "WP");
    auto table = config.get<fs::path>("corrections.photons.conversion_veto.table");
    auto type = config.get<string>("corrections.photons.conversion_veto.type");
    metainfo.Set<string>("corrections", "photons", "conversion_veto", "type", type);
    metainfo.Set<fs::path>("corrections", "photons", "conversion_veto", "table", table);
    
    ConversionVetoApplier applier(table, ID, table != "/dev/null",
                          steering & DT::syst, type);
    for (const auto& name: applier.weightNames())
        metainfo.Set<string>("variations", RecPhoton::WeightVar, name);
 
 
    for (DT::Looper looper(tIn, slice); looper(); ++looper) {
        [[ maybe_unused]]
        static auto& cout = steering & DT::verbose ? ::cout : DT::dev_null;
 
        applier(*recPhotons);
 
        static const bool fillTree = steering & DT::fill;
        if (fillTree) tOut->Fill();
    }
 
    metainfo.Set<bool>("git", "complete", true);
    fOut->cd();
    tOut->Write();
 
    cout << __func__ << ' ' << slice << " stop" << endl;
}

applyPhotonID()

void DAS::Photon::applyPhotonID	(	const vector< fs::path > &	inputs,
		const fs::path &	output,
		const pt::ptree &	config,
		const int	steering,
		const DT::Slice	slice = {1,0}
	)

Todo:

ControlPlots

Parameters

inputs	input ROOT files (n-tuples)
output	output ROOT file (n-tuple)
config	config handled with `Darwin::Tools::options`
steering	parameters obtained from explicit options
slice	number and index of slice

                              {1,0} 
        )
{
    cout << __func__ << ' ' << slice << " start" << endl;
 
    unique_ptr<TChain> tIn = DT::GetChain(inputs);
    unique_ptr<TFile> fOut(DT_GetOutput(output));
    auto tOut = unique_ptr<TTree>(tIn->CloneTree(0));
 
    DT::MetaInfo metainfo(tOut);
    metainfo.Check(config);
 
    vector<RecPhoton> * recPhotons = nullptr;
    tIn->SetBranchAddress("recPhotons", &recPhotons);
 
    auto ID = config.get<string>("corrections.photons.ID.WP");
    auto path = config.get<fs::path>("corrections.photons.ID.table");
    metainfo.Set<string>("corrections", "photons", "ID", "WP", ID);
    metainfo.Set<fs::path>("corrections", "photons", "ID", "table", path);
    IDApplier applier(path, ID, path != "/dev/null", steering & DT::syst);
    for (const auto& name: applier.weightNames())
        metainfo.Set<string>("variations", RecPhoton::WeightVar, name);
 
 
    for (DT::Looper looper(tIn, slice); looper(); ++looper) {
        [[ maybe_unused]]
        static auto& cout = steering & DT::verbose ? ::cout : DT::dev_null;
 
        applier(*recPhotons);
 
        static const bool fillTree = steering & DT::fill;
        if (fillTree) tOut->Fill();
    }
 
    metainfo.Set<bool>("git", "complete", true);
    fOut->cd();
    tOut->Write();
 
    cout << __func__ << ' ' << slice << " stop" << endl;
}