| 4 |
|
#include <TGraphErrors.h> |
| 5 |
|
#include <TLegend.h> |
| 6 |
|
#include <iostream> |
| 7 |
< |
#include "CutsAndHistos.h" |
| 8 |
< |
#include "Histos.h" |
| 9 |
< |
#include "Cuts.h" |
| 7 |
> |
#include "../interface/CutsAndHistos.h" |
| 8 |
> |
#include "../plugins/Histos.h" |
| 9 |
> |
#include "../plugins/Cuts/Cuts.hpp" |
| 10 |
|
#include "TF1.h" |
| 11 |
|
#include "TH1.h" |
| 12 |
< |
#include "samples.h" |
| 13 |
< |
#include "ntupleReader.h" |
| 14 |
< |
#include "sampleCollection.h" |
| 12 |
> |
#include "../interface/samples.hpp" |
| 13 |
> |
#include "../interface/ntupleReader.hpp" |
| 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){ |
