|
DAS
3.0
Das Analysis System
|
#include <HTn.h>
◆ HTnFiller()
Constructor.
30 ,
rEv(reader,
"recEvent")
◆ fillMC()
See Filler::fillMC
Reimplemented from Filler.
130 BOOST_THROW_EXCEPTION( runtime_error(__func__ +
" should only be called for MC"s) );
133 auto [recn,recW] = getMultiplicity<RecJet>(
recJets, v);
135 recW *=
rEv->Weight(v);
139 auto [genn,genW] = getMultiplicity<GenJet>(*
genJets, v);
141 genW = (*gEv)->Weight(v);
144 if (igen > 0) v.
gen->Fill(igen, genW);
147 if (irec > 0 && igen > 0) { v.
RM->Fill(igen, irec, genW * recW );
149 else if (irec == 0 && igen > 0) v.
missOut->Fill(igen, genW );
150 else if (irec > 0 && igen == 0) v.
fakeOut->Fill( irec, genW * recW );
154 if (irec > 0) v.
fakeNoMatch->Fill(irec, genW * recW);
◆ fillRec()
See Filler::fillRec
Reimplemented from Filler.
110 auto dijet = selection<RecJet>(
recJets, v);
111 if (!dijet)
return {};
113 auto [n,
w] = getMultiplicity(
recJets, v);
115 if (i == 0)
return {};
118 if (
obs.
isMC)
w *= (*gEv)->Weight(v);
119 w *= dijet.Weight(v);
◆ match()
Match the two pairs of leading jets (if any) and set matched
member.
Reimplemented from Filler.
38 auto match = [
this](
size_t i,
size_t j) {
41 using ROOT::Math::VectorUtil::DeltaR;
42 auto DR = DeltaR(g, r);
◆ fakes
◆ genJets
std::optional<TTreeReaderArray<GenJet> > genJets |
◆ gEv
std::optional<TTreeReaderValue<GenEvent> > gEv |
◆ matched
std::optional<bool> matched |
◆ matches
◆ misses
◆ obs
Backreference to the observable.
◆ recJets
TTreeReaderArray<RecJet> recJets |
◆ rEv
The documentation for this struct was generated from the following files:
- /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Unfolding/interface/HTn.h
- /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Unfolding/src/HTn.cc
double getBinNumber(const Di< Jet, Jet > &MNjets, const Uncertainties::Variation &v, TUnfoldBinning *bng)
Definition: MNjets.cc:47
double Weight(const Uncertainties::Variation &v=Uncertainties::nominal) const override
Definition: Di.h:73
std::unique_ptr< TH1 > rec
reconstructed-level distribution
Definition: DistVariation.h:32
TTreeReaderValue< RecEvent > rEv
Definition: HTn.h:76
HTn obs
Backreference to the observable.
Definition: HTn.h:71
static const float w
Definition: common.h:47
TUnfoldBinning * genBinning
particle-level binning
Definition: Observable.h:150
TTreeReaderArray< RecJet > recJets
Definition: HTn.h:74
std::optional< TTreeReaderValue< GenEvent > > gEv
Definition: HTn.h:75
std::optional< TTreeReaderArray< GenJet > > genJets
Definition: HTn.h:73
std::optional< bool > matched
Definition: HTn.h:90
static bool isMC
Definition: Observable.h:144
auto initOptionalBranch(TTreeReader &reader, const char *name)
Definition: Observable.h:39
DAS::RecDijet recdijet
Definition: classes.h:34
std::unique_ptr< TH1 > gen
generated-level distribution
Definition: DistVariation.h:34
std::unique_ptr< TH1 > tmp
temporary histogram help fill the covariance matrix
Definition: DistVariation.h:33
std::unique_ptr< TH1 > missNoMatch
losses (unmatched entries)
Definition: DistVariation.h:35
std::unique_ptr< TH2 > RM
response matrix
Definition: DistVariation.h:42
std::unique_ptr< TH1 > missOut
losses (migration out of phase space)
Definition: DistVariation.h:36
std::unique_ptr< TH1 > fakeNoMatch
background (unmatched entries)
Definition: DistVariation.h:37
ROOT::Math::LorentzVector< ROOT::Math::PtEtaPhiM4D< float > > FourVector
Definition: PhysicsObject.h:15
TUnfoldBinning * recBinning
detector-level binning
Definition: Observable.h:149
std::unique_ptr< TH1 > fakeOut
background (migration out of phase space)
Definition: DistVariation.h:38
DAS::GenDijet gendijet
Definition: classes.h:31
void match() override
Match the two pairs of leading jets (if any) and set matched member.
Definition: HTn.cc:34
static double maxDR
max Delta R
Definition: Observable.h:146