DistVariation.h

Go to the documentation of this file.

#pragma once
 
#include <iostream>
#include <list>
#include <memory>
 
#include "Core/Objects/interface/Variation.h"
 
#include <TString.h>
 
#include <MetaInfo.h>
 
class TH1;
class TH2;
class TDirectory;
class TUnfoldBinning;
 
namespace DAS::Unfolding {
 
struct DistVariation : public DAS::Uncertainties::Variation {
 
    static bool isMC; 
 
    static TUnfoldBinning * genBinning, 
                          * recBinning; 
 
    std::unique_ptr<TH1> rec, 
                         tmp, 
                         gen, 
                         missNoMatch, 
                         missOut, 
                         fakeNoMatch, 
                         fakeOut; 
    // TODO: include miss entries in RM
    // TODO: use TUnfoldBinning to decorrelate categories of miss/fake entries
    std::unique_ptr<TH2> cov, 
                         RM; 
 
    ~DistVariation ();
 
    DistVariation (const TString& group, const TString& name, size_t index = 0, int bit = 0);
 
    DistVariation (const DistVariation& v) :
        DAS::Uncertainties::Variation(v),
        rec        (std::unique_ptr<TH1>(dynamic_cast<TH1*>(v.rec        ->Clone()))),
        tmp        (std::unique_ptr<TH1>(dynamic_cast<TH1*>(v.tmp        ->Clone()))),
        gen        (std::unique_ptr<TH1>(dynamic_cast<TH1*>(v.gen        ->Clone()))),
        missNoMatch(std::unique_ptr<TH1>(dynamic_cast<TH1*>(v.missNoMatch->Clone()))),
        missOut    (std::unique_ptr<TH1>(dynamic_cast<TH1*>(v.missOut    ->Clone()))),
        fakeNoMatch(std::unique_ptr<TH1>(dynamic_cast<TH1*>(v.fakeNoMatch->Clone()))),
        fakeOut    (std::unique_ptr<TH1>(dynamic_cast<TH1*>(v.fakeOut    ->Clone()))),
        cov        (std::unique_ptr<TH2>(dynamic_cast<TH2*>(v.cov        ->Clone()))),
        RM         (std::unique_ptr<TH2>(dynamic_cast<TH2*>(v.RM         ->Clone())))
    { }
 
    DistVariation (DistVariation&& v) :
        DAS::Uncertainties::Variation(v),
        rec(std::move(v.rec)), tmp(std::move(v.tmp)), gen(std::move(v.gen)),
        missNoMatch(std::move(v.missNoMatch)), missOut(std::move(v.missOut)),
        fakeNoMatch(std::move(v.fakeNoMatch)), fakeOut(std::move(v.fakeOut)),
        cov(std::move(v.cov)), RM(std::move(v.RM))
    { }
 
    void Write (TDirectory * d);
};
 
void fillCov (const DistVariation&, const std::list<int>&); // TODO: make this function a static member of Observable
 
std::vector<DistVariation> GetVariations
            (Darwin::Tools::MetaInfo&, 
             bool = false, 
             std::ostream& = std::cout);
 
} // end of namespace DAS::Unfolding