|
DAS
3.0
Das Analysis System
|
#include <ControlPlots.h>
|
| ControlPlots (TString Name) |
|
void | operator() (const std::vector< RecJet > &recjets, const double &evW, size_t iJEC=0, size_t iWgt=0) |
|
void | operator() (const std::vector< GenJet > &genjets, const double &evW, size_t iWgt=0) |
|
void | Write (TDirectory *D) |
|
◆ ControlPlots()
- Todo:
- TH3 pt eta rho plots TH2 pt_N TH2 ptmax_Dphi
15 cout << __func__ <<
'\t' << Name << endl;
◆ operator()() [1/2]
void operator() |
( |
const std::vector< GenJet > & |
genjets, |
|
|
const double & |
evW, |
|
|
size_t |
iWgt = 0 |
|
) |
| |
47 auto genjet0 = find_if(begin(
genjets), end(
genjets), inYacceptance);
48 genpt0_N->Fill(genjet0->p4.Pt(), N, evW * genjet0->weights.front());
51 for (
size_t i = 0; i <
genjets.size(); ++i) {
61 if ( j0.p4.Pt() < 100. || j1.p4.Pt() < 50. )
return;
63 auto ymax = max(j0.AbsRap(), j1.AbsRap());
64 auto djW = j0.weights.at(iWgt) * j1.weights.at(iWgt);
67 auto Y0 = j0.Rapidity(),
69 auto ystar = 0.5*abs(Y0-Y1),
70 yboost = 0.5*abs(Y0+Y1);
71 genMjj_ysyb->Fill(dijet.CorrP4().M(), yboost, ystar, evW * djW);
72 auto ptav = 0.5*(j0.p4.Pt() + j1.p4.Pt());
◆ operator()() [2/2]
void operator() |
( |
const std::vector< RecJet > & |
recjets, |
|
|
const double & |
evW, |
|
|
size_t |
iJEC = 0 , |
|
|
size_t |
iWgt = 0 |
|
) |
| |
89 auto recjet0 = find_if(begin(
recjets), end(
recjets), inYacceptance);
90 recpt0_N->Fill(recjet0->p4.Pt(), N, evW * recjet0->weights.front());
94 for (
size_t i = 0; i <
recjets.size(); ++i) {
100 if (
recjets.size() < 2)
return;
101 const auto& j0 =
recjets.at(0),
103 if ( j0.CorrPt(iJEC) < 100. || j1.CorrPt(iJEC) < 50. )
return;
104 auto ymax = max(j0.AbsRap(), j1.AbsRap());
105 auto djW = j0.weights.at(iWgt) * j1.weights.at(iWgt);
106 auto dijet = j0 + j1;
109 auto Y0 = j0.Rapidity(),
111 auto ystar = 0.5*abs(Y0-Y1),
112 yboost = 0.5*abs(Y0+Y1);
113 recMjj_ysyb->Fill(dijet.CorrP4().M(), yboost, ystar, evW * djW);
114 auto ptav = dijet.HT();
◆ Write()
void Write |
( |
TDirectory * |
D | ) |
|
121 auto d = D->mkdir(
name);
126 hs.insert(hs.end(), hs2.begin(), hs2.end());
129 TString n = h->GetName();
130 n.ReplaceAll(
name,
"");
◆ genMjj_y
std::optional<TH2F> genMjj_y |
◆ genMjj_ysyb
std::optional<TH3F> genMjj_ysyb |
◆ genpt
std::optional<TH1F> genpt |
◆ genpt0_N
std::optional<TH2F> genpt0_N |
◆ genpt_n
std::optional<TH2F> genpt_n |
◆ genpt_y
std::optional<TH2F> genpt_y |
◆ genptav_ysyb
std::optional<TH3F> genptav_ysyb |
◆ isMC
◆ name
◆ recMjj_y
std::optional<TH2F> recMjj_y |
◆ recMjj_ysyb
std::optional<TH3F> recMjj_ysyb |
◆ recpt
std::optional<TH1F> recpt |
◆ recpt0_N
std::optional<TH2F> recpt0_N |
◆ recpt_n
std::optional<TH2F> recpt_n |
◆ recpt_y
std::optional<TH2F> recpt_y |
◆ recptav_ysyb
std::optional<TH3F> recptav_ysyb |
◆ verbose
The documentation for this struct was generated from the following files:
std::optional< TH2F > genMjj_y
Definition: ControlPlots.h:26
static const std::vector< double > y_edges
Definition: binnings.h:36
std::optional< TH1F > recpt
Definition: ControlPlots.h:24
std::optional< TH3F > genptav_ysyb
Definition: ControlPlots.h:28
FourVector p4
raw four-momentum directly after reconstruction
Definition: PhysicsObject.h:47
static const std::vector< double > n_edges
Definition: binnings.h:37
std::optional< TH2F > genpt_y
Definition: ControlPlots.h:25
static const std::vector< double > pt_edges
Definition: binnings.h:33
static const std::vector< double > pt_av_edges
Definition: binnings.h:35
static const int nMjjBins
Definition: binnings.h:40
static const int maxMult
Definition: binnings.h:42
std::optional< TH3F > genMjj_ysyb
Definition: ControlPlots.h:27
static const double ymax
Definition: binnings.h:49
std::optional< TH2F > recpt_y
Definition: ControlPlots.h:25
static bool verbose
Definition: ControlPlots.h:20
float AbsRap(const Uncertainties::Variation &=Uncertainties::nominal) const
absolute rapidity
Definition: PhysicsObject.h:54
genjets
Definition: Ntupliser_cfg.py:268
std::optional< TH1F > genpt
Definition: ControlPlots.h:24
DAS::RecJet recjet
Definition: classes.h:15
std::optional< TH2F > recpt0_N
Definition: ControlPlots.h:29
static bool isMC
Definition: ControlPlots.h:19
std::optional< TH2F > genpt0_N
Definition: ControlPlots.h:29
recjets
Definition: Ntupliser_cfg.py:269
static const std::vector< double > Mjj_edges
Definition: binnings.h:34
float CorrPt(size_t i=0) const
corrected transverse momentum
Definition: PhysicsObject.h:52
static const int nPtBins
Definition: binnings.h:39
std::optional< TH3F > recMjj_ysyb
Definition: ControlPlots.h:27
std::optional< TH2F > genpt_n
Definition: ControlPlots.h:30
Weights weights
object weights
Definition: PhysicsObject.h:49
DAS::GenJet genjet
Definition: classes.h:14
static const int nPtDijetBins
Definition: binnings.h:41
std::optional< TH3F > recptav_ysyb
Definition: ControlPlots.h:28
const TString name
Definition: ControlPlots.h:22
std::optional< TH2F > recpt_n
Definition: ControlPlots.h:30
static const int nYbins
Definition: binnings.h:38
std::optional< TH2F > recMjj_y
Definition: ControlPlots.h:26