ZmmYDalitz.h

Go to the documentation of this file.

// SPDX-License-Identifier: GPLv3-or-later
//
// SPDX-FileCopyrightText: Louis Moureaux <louis.moureaux@cern.ch>
 
#pragma once
 
#include <vector>
#include <list>
 
#include <TUnfoldBinning.h>
#include <TH1.h>
#include <TH2.h>
 
#if !defined(__CLING__) || defined(__ROOTCLING__)
#include "Core/Objects/interface/Event.h"
#include "Core/Objects/interface/Lepton.h"
#include "Core/Objects/interface/Photon.h"
#include "Core/Unfolding/interface/Observable.h"
#include "Core/Unfolding/interface/ZmmYCommon.h"
#endif
 
namespace DAS::Unfolding::ZmmY {
 
// Phase space cuts
constexpr double minMu1Pt = 20;
constexpr double minMu2Pt = 15;
constexpr double maxMuAbsEta = 2.4;
constexpr double minGammaPt = 20;
constexpr double maxGammaAbsEta = 2.4;
constexpr double minMMuMu = 40;
constexpr double minMMuMuGamma = 76.1876, maxMMuMuGamma = 106.1876;
 
// Binning
constexpr int nM2RatioBins = 20;
constexpr int nMMuGamma2Bins = 10;
 
#if !defined(__CLING__) || defined(__ROOTCLING__)
struct Dalitz final : public Observable {
 
    Dalitz ();
 
    std::unique_ptr<Filler> getFiller (Darwin::Tools::Flow&) const override;
};
 
struct DalitzFiller final : public Filler {
    Dalitz obs; 
    EventInfo<RecMuon, RecPhoton> rInfo; 
 
    std::vector<GenMuon> * genMuons;
    std::vector<RecMuon> * recMuons;
    std::vector<GenPhoton> * genPhotons;
    std::vector<RecPhoton> * recPhotons;
    GenEvent * gEv;
    RecEvent * rEv;
 
    DalitzFiller (const Dalitz&, Darwin::Tools::Flow&);
 
    std::list<int> fillRec (DistVariation&) override;
 
    void match () override {}
 
    void fillMC (DistVariation&) override;
};
#endif
 
} // end of DAS::Unfolding::DrellYan namespace