5 |
|
#include <TLegend.h> |
6 |
|
#include <iostream> |
7 |
|
#include "../interface/CutsAndHistos.h" |
8 |
< |
#include "../interface/Histos.h" |
9 |
< |
#include "../plugins/Cuts/Cuts.h" |
8 |
> |
#include "../plugins/Histos.h" |
9 |
> |
#include "../plugins/Cuts/Cuts.hpp" |
10 |
|
#include "TF1.h" |
11 |
|
#include "TH1.h" |
12 |
|
#include "../interface/samples.hpp" |
13 |
|
#include "../interface/ntupleReader.hpp" |
14 |
< |
#include "sampleCollection.h" |
14 |
> |
#include "../test/SideBandAnalysis-Pt50To100/sampleSideBand.h" |
15 |
|
#include "Riostream.h" |
16 |
|
#include "TMatrixD.h" |
17 |
|
#include "TVectorD.h" |
18 |
|
#include "TDecompChol.h" |
19 |
|
#include "TDecompSVD.h" |
20 |
< |
#include "setTDRStyle.C" |
21 |
< |
#include "controlRegions.h" |
20 |
> |
#include "../plugins/setTDRStyle.C" |
21 |
> |
#include "../interface/controlRegions.h" |
22 |
|
|
23 |
|
#define ZeeL 4683.5 |
24 |
|
#define fA 0.46502 |
27 |
|
int main(int argc, char **argv) |
28 |
|
{ |
29 |
|
|
30 |
< |
typedef std::vector<Sample> samplesCollection; |
30 |
> |
typedef std::vector<Sample> sampleCollection; |
31 |
|
setTDRStyle(); |
32 |
|
bool verbose_ = true; |
33 |
|
TCanvas * c1 = new TCanvas("c1","c1", 600,600); |
39 |
|
signalRegionCutFlow.add( new VMassCutMin(75.) ); |
40 |
|
signalRegionCutFlow.add( new VMassCutMax(105.) ); |
41 |
|
signalRegionCutFlow.add( new JetPtCut(20.) ); |
42 |
– |
signalRegionCutFlow.add( new HPtCut(100.) ); |
42 |
|
signalRegionCutFlow.add( new VPtCut(100.) ); |
43 |
< |
signalRegionCutFlow.add( new JetBtagCut(0.5, 0.898) ); |
44 |
< |
signalRegionCutFlow.add( new HVdPhiCut(2.9) ); |
45 |
< |
signalRegionCutFlow.add( new JetVeto(2) ); |
47 |
< |
signalRegionCutFlow.add( new HMassCut(100.) ); |
48 |
< |
signalRegionCutFlow.add( new HMassCutMax(130.) ); |
43 |
> |
signalRegionCutFlow.add( new JetBtagCut(0.244, 0.244) ); |
44 |
> |
signalRegionCutFlow.add( new HMassCut(80.) ); |
45 |
> |
signalRegionCutFlow.add( new HMassCutMax(120.) ); |
46 |
|
|
47 |
< |
samplesCollection samples = Nov10thDiJetPtUpdatedFall11OKSlim(); |
47 |
> |
sampleCollection samples = trees(); |
48 |
|
std::vector< controlRegion > cutFlowCR; |
49 |
|
|
50 |
|
for(int i=0; i<signalRegionCutFlow.size(); ++i){ |