applyEventMixing.cc File Reference

#include <cstdlib>
#include <cassert>
#include <filesystem>
#include <iostream>
#include <map>
#include <vector>
#include "Core/Objects/interface/Event.h"
#include "Core/Objects/interface/Di.h"
#include "Core/Objects/interface/Jet.h"
#include "Core/Objects/interface/Lepton.h"
#include "Core/CommonTools/interface/binnings.h"
#include "Core/CommonTools/interface/DASOptions.h"
#include "Core/CommonTools/interface/toolbox.h"
#include "Core/JetObservables/bin/common.h"
#include "Core/JetObservables/interface/fjcore.hh"
#include <TFile.h>
#include <TH1.h>
#include <TString.h>
#include "Math/VectorUtil.h"
#include <darwin.h>

Include dependency graph for applyEventMixing.cc:

Namespaces
	DAS
	DAS::TPS

Enumerations
enum	RemovalStrategy { DeltaR, AntikT, None }

Functions
RemovalStrategy	parseStrategy (const string &strategyStr, TTree *tree)
template<class Jet >
vector< fjcore::PseudoJet >	cluster (const vector< Jet > &jets, const double R)
template<class Jet >
bool	overlappingJets (const vector< Jet > &jets)
map< Long64_t, vector< Long64_t > >	buildRunToEntriesMap (TTree *tree, RecEvent *recEvent)
void	applyEventMixing (const vector< fs::path > &inputs1, const vector< fs::path > &inputs2, 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[]
	)

{
    try {
        DT::StandardInit();
 
        vector<fs::path> inputs1;
        vector<fs::path> inputs2;
        fs::path output;
 
        auto options = DAS::Options(
                            "Mix the jet collections of two events into one jet collection in one events."
                            "Event weights are multiplied according to the DPS pocket formula (modulo the "
                            "combinatorial prefactor and the sigma effective constant). Extra information "
                            "is untouched and should be handled by the user.\n"
                            "Overlapping jets can be treated by three different ways:\n"
                            " - DeltaR: drop events when facing azimuthally overlapping jets\n"
                            " - antikt: drop events when the jet multiplicity changes before and after clustering the jets with the anti-kt algorithm.\n"
                            " - none: do nothing\n",
                            DT::config | DT::split | DT::fill);
        options.inputs("inputs1", &inputs1, "input ROOT file(s) or directory")
               .inputs("inputs2", &inputs2, "input ROOT file(s) or directory")
               .output("output", &output, "output ROOT file")
               .arg<float>("minpt", "skims.jets.minpt", "Minimum pT of the last jet")
               .arg<string>("RemovalStrategy", "corrections.event_mixing.overlap", "'deltaR' or 'antikt' or 'none'");
        const auto& config = options(argc, argv);
        const auto& slice = options.slice();
        const int steering = options.steering();
 
        DAS::TPS::applyEventMixing(inputs1, inputs2, output, config, steering, slice);
    }
    catch (boost::exception& e) {
        DE::Diagnostic(e);
        return EXIT_FAILURE;
    }
 
    return EXIT_SUCCESS;
}