applyPUcleaning.cc File Reference

#include <cstdlib>
#include <filesystem>
#include <iostream>
#include <unordered_map>
#include <utility>
#include <vector>
#include "Core/Objects/interface/Event.h"
#include "Core/Objects/interface/Jet.h"
#include "Core/CommonTools/interface/binnings.h"
#include "Core/CommonTools/interface/ControlPlots.h"
#include "Core/CommonTools/interface/DASOptions.h"
#include "Core/CommonTools/interface/toolbox.h"
#include "Core/PUstaubSauger/interface/Plots.h"
#include <boost/property_tree/ptree.hpp>
#include <boost/property_tree/info_parser.hpp>
#include <TH1.h>
#include <TString.h>
#include <darwin.h>

Include dependency graph for applyPUcleaning.cc:

Namespaces
	DAS
	DAS::PUstaub

Functions
unordered_map< unsigned long long, float >	GetMaxGenPts (const fs::path &input)
void	applyPUcleaning (const vector< fs::path > &inputs, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Note

156 GeV for Pythia UL18

{
    try {
        DT::StandardInit();
 
        vector<fs::path> inputs;
        fs::path output;
 
        auto options = DAS::Options(
                "Clean up events from anomalous pile-up in MC samples.",
                DT::config | DT::split | DT::Friend);
        options.inputs("inputs", &inputs, "input ROOT file(s) or directory")
               .output("output", &output, "output ROOT file")
               .arg<fs::path>("MBmaxgenpt", "corrections.PUcleaning.MBmaxgenpt",
                                "2-column file with MB event ID and max gen pt")
               .arg<float>("artifacts_minpt", "corrections.PUcleaning.artifacts_minpt",
                            "Partily arbitrary value to cut of overweighted reconstruction "
                            "artifacts (0 to deactivate)"); 
 
        const auto& config = options(argc, argv);
        const auto& slice = options.slice();
        const int steering = options.steering();
 
        DAS::PUstaub::applyPUcleaning(inputs, output, config, steering, slice);
    }
    catch (boost::exception& e) {
        DE::Diagnostic(e);
        return EXIT_FAILURE;
    }
 
    return EXIT_SUCCESS;
}