[ViewVC] Diff of: cvsroot/UserCode/MitHzz4l/NonMCBackground/src/plotSignalFakes.cc

Comparing UserCode/MitHzz4l/NonMCBackground/src/plotSignalFakes.cc (file contents):
Revision 1.4 by dkralph, Tue Jul 31 16:46:22 2012 UTC vs.
Revision 1.7 by dkralph, Mon Dec 17 17:25:55 2012 UTC

+#include "TCanvas.h"
+#include "TChain.h"
-+
+#include "TStopwatch.h"
+#include "TString.h"
+#include "TStyle.h"
+#include "TH1D.h"
+#include "TNtuple.h"
-+
+#include "TROOT.h"
+#include "CPlot.h"
+#include "FOArgs.h"
+#include "Various.h"
+#include "CommonDefs.h"
+#include "PlotHeaders.h"
-+
+#include "HistHeaders.h"
+#include "TriggerUtils.h"
+#include "JetDefs.h"
+#include "JetInfoStruct.h"
-+
+#include "FR_struct.h"
-+
+#include "SelectionFuncs.h"
-+
+#include "SampleWeight.h"
-+
+#include "MitStyleRemix.h"
+#include "TMVA/Reader.h"
+#include "TMVA/Tools.h"
+using namespace RooFit;
+using namespace mithep;
-–
+TH1D* hpu_2011;
-–
+TH1D* hpu_2012;
-–
+//----------------------------------------------------------------------------------------
-–
+class EvtsOfInterest
-–
+{
-–
+public:
-–
+  EvtsOfInterest(TString config="");
-–
+  bool has(unsigned run, unsigned evt);
-–
+  vector<unsigned> runs,evts;
-–
+};
-–
+EvtsOfInterest::EvtsOfInterest(TString config)
-–
+{
-–
+  if(config!="") {
-–
+    ifstream ifs(config);
-–
+    assert(ifs.is_open());
-–
+    string line;
-–
+    while(getline(ifs,line)) {
-–
+      stringstream ss(line);
-–
+      unsigned run,evt;
-–
+      ss >> run >> evt;
-–
+      runs.push_back(run);
-–
+      evts.push_back(evt);
-–
+    }
-–
+    ifs.close();
-–
+  }
-–
+}
-–
+bool EvtsOfInterest::has(unsigned run, unsigned evt)
-–
+{
-–
+  for(unsigned i=0; i<runs.size(); i++) {
-–
+    if(runs[i]==run && evts[i]==evt)
-–
+      return true;
-–
+  }
-–
+  return false;
-–
+}
+//----------------------------------------------------------------------------------------
+TCanvas *can;
-–
+//----------------------------------------------------------------------------------------
-–
+class lepFrs
-–
+{
-–
+public:
-–
+  double fwgt_3,fwgt_4;
-–
+  double fwgt_lo_3,fwgt_lo_4;
-–
+  double fwgt_hi_3,fwgt_hi_4;
-–
+};
-–
+//----------------------------------------------------------------------------------------
-–
+class dr_struct
-–
+{
-–
+public:
-–
+  dr_struct();
-–
+  void check();
-–
+  float l1j1,l2j1,l3j1,l4j1;
-–
+  float l1j2,l2j2,l3j2,l4j2;
-–
+  float l1j3,l2j3,l3j3,l4j3;
-–
+  float l1j4,l2j4,l3j4,l4j4;
-–
+};
-–
+dr_struct::dr_struct()
-–
+{
-–
+  l1j1 = l2j1 = l3j1 = l4j1 = l1j2 = l2j2 = l3j2 = l4j2 = l1j3 = l2j3 = l3j3 = l4j3 = l1j4 = l2j4 = l3j4 = l4j4 = -1;
-–
+}
-–
+void dr_struct::check()
-–
+{
-–
+  if(fabs(l1j1) < 0.2) { cout << "l1j1 too small: " << l1j1 << endl; }
-–
+  if(fabs(l2j1) < 0.2) { cout << "l2j1 too small: " << l2j1 << endl; }
-–
+  if(fabs(l3j1) < 0.2) { cout << "l3j1 too small: " << l3j1 << endl; }
-–
+  if(fabs(l4j1) < 0.2) { cout << "l4j1 too small: " << l4j1 << endl; }
-–
+  if(fabs(l1j2) < 0.2) { cout << "l1j2 too small: " << l1j2 << endl; }
-–
+  if(fabs(l2j2) < 0.2) { cout << "l2j2 too small: " << l2j2 << endl; }
-–
+  if(fabs(l3j2) < 0.2) { cout << "l3j2 too small: " << l3j2 << endl; }
-–
+  if(fabs(l4j2) < 0.2) { cout << "l4j2 too small: " << l4j2 << endl; }
-–
+  if(fabs(l1j3) < 0.2) { cout << "l1j3 too small: " << l1j3 << endl; }
-–
+  if(fabs(l2j3) < 0.2) { cout << "l2j3 too small: " << l2j3 << endl; }
-–
+  if(fabs(l3j3) < 0.2) { cout << "l3j3 too small: " << l3j3 << endl; }
-–
+  if(fabs(l4j3) < 0.2) { cout << "l4j3 too small: " << l4j3 << endl; }
-–
+  if(fabs(l1j4) < 0.2) { cout << "l1j4 too small: " << l1j4 << endl; }
-–
+  if(fabs(l2j4) < 0.2) { cout << "l2j4 too small: " << l2j4 << endl; }
-–
+  if(fabs(l3j4) < 0.2) { cout << "l3j4 too small: " << l3j4 << endl; }
-–
+  if(fabs(l4j4) < 0.2) { cout << "l4j4 too small: " << l4j4 << endl; }
-–
+}
+bool findGoodJets(vector<SimpleLepton> &goodJets, filestuff *fs,
+                  SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4);
-<
+void fillHist(CSample *cs, TString channel, TString type, TString var, double val, double wgt, double wgt_lo=0, double wgt_hi=0);
-<
+void fillAllHists(FOFlags ctrl, CSample *cs, TString channel, TString type, filestuff *fs, KinematicsStruct kine,
-<
+                  SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4,
-<
+                  double wgt, double wgt_lo=0, double wgt_hi=0);
-<
+void fillAllJetHists(FOFlags ctrl,  CSample *cs, TString channel, TString type, filestuff *fs, KinematicsStruct kine,
-<
+                      SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4,
-<
+                      FusionMva &fusion, vector<SimpleLepton> &goodJets, JetInfoStruct ji,
-<
+                      double wgt, double wgt_lo=0, double wgt_hi=0);
-<
+void makeHTML(FOFlags &ctrl, TString type, TString plotLabel);
->
+// void fillHist(CSample *cs, TString channel, TString type, TString var, double val, double wgt, double wgt_lo=0, double wgt_hi=0);
->
+// void fillAllHists(FOFlags ctrl, CSample *cs, TString channel, TString type, filestuff *fs, KinematicsStruct kine, Angles angles,
->
+//                SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4,
->
+//                double wgt, double wgt_lo=0, double wgt_hi=0);
->
+// void fillAllJetHists(FOFlags ctrl,  CSample *cs, TString channel, TString type, filestuff *fs, KinematicsStruct kine,
->
+//                    SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4,
->
+//                    FusionMva &fusion, vector<SimpleLepton> &goodJets, JetInfoStruct ji,
->
+//                    double wgt, double wgt_lo=0, double wgt_hi=0);
->
+// void makeHTML(FOFlags &ctrl, TString type, TString plotLabel, TString fullOutDir);
+map<TString,map<TString,TH1D*>* > init_hists(FOFlags &ctrl, TString str="");
+bool passHlt(FOFlags &ctrl, TrigInfo ti, InfoStruct *info, unsigned lep1matchBits,
+             unsigned lep2matchBits, unsigned lep3matchBits,
+             unsigned lep4matchBits);
-–
+dr_struct fill_dr_struct(vector<SimpleLepton> jets, SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4);
+void init_cuts(vector<TString> &cutstrs, map<TString,int> &cutvec);
-–
+void fillLeptonFakeWeights(lepFrs &fwgts, FR_struct *fr, SimpleLepton lep3, SimpleLepton lep4);
+void resetCutVect(map<TString,int> &cutvec) {
+  for(map<TString,int>::iterator it=cutvec.begin(); it!=cutvec.end(); it++)
+    cutvec[(*it).first] = 0;
+map<TString,TH1D*> setHistSet(CSample *cs, TString type, TString var);
-–
+void shiftOverflows(map<TString,TH1D*> &hset);
-–
+void scaleHists(map<TString,TH1D*> &hset);
+void fillFakeTuple(CSample *cs, filestuff *fs, EventData evtdata, unsigned ientry, double fillweight);
-+
+void incrementSsofCounters(FOFlags &ctrl, filestuff *fs, double minMz2,
-+
+                           pair<int,int> *best_z_indices, double wgt);
-+
+map<TString,double> init_counter(TString fname="");
-+
+map<TString,double> initSsofRatios(vector<CSample*> &samplev);
-+
+double getSsofWgt(TString channel, map<TString,double> &ssofRatios);
+//----------------------------------------------------------------------------------------
+int main(int argc, char** argv)
-+
+  // Three usage modes:
-+
+  //   ctrl.faketype=="SR":    plot events in signal region, if requested including ntuples of the extrapolation from 2P2F and 3P1F regions
-+
+  //   ctrl.faketype=="2P2F":  plot 2P2F region, and if requested write out ntuple which gives extrapolation to SR (ntuple would then be read by "SR", above)
-+
+  //   ctrl.faketype=="3P1F":  plot 3P1F region (including the extrapolation to there from the 2P2F region), and if requested write out ntuple with SR prediction
-+
+  //
-+
+  // ctrl.makeFakeTuples:  write out zznt-type ntuples to be read by plotH4l or by SR mode
-+
+  // ctrl.hiStatFakes:     same as makeFakeTuples, but reverse charge and flavor req's to get a larger sample for BDT training (or to plot the WrongFlavorCharge control region)
-+
+  // ctrl.ssof:            look only at ssof region
-+
+  // ctrl.writessofratio:  write the OS / SS ratio to text files that can be read later in ssof mode
-+
+  // ctrl.plotWholeSample: plot the whole fake sample without weighting with the fake rate
-+
-+
+  TStopwatch watch;
-+
+  SetStyle();
+  initPUWeights();
-+
+  ControlFlags dummyCtrl;
-+
+  vector<TString> cutstrs;
+  map<TString,int> cutvec;
+  init_cuts(cutstrs, cutvec);
-<
+  EvtsOfInterest eoiAdish("/tmp/adish-uniq");
-<
+  EvtsOfInterest eoiMe("/tmp/me-uniq");
->
+  bool dummyStr(true); // don't add the yield numbers to legends
+  // arguments...
+  FOFlags ctrl;
+  parse_foargs( argc, argv, ctrl );
+  ctrl.dump();
-–
+  bool makeFakeTuples = (ctrl.extraArgs.Contains("makeFakeTuples"));
+  FusionMva fusion(ctrl.uncert);//"./weights/againstZZ-fusion_BDTG.weights.xml");
+  bool makeJetTuple=false; // jet tuple for vbf mva training
-<
+  assert(ctrl.faketype=="SR"   || // plot events in signal region, if requested including ntuples of the extrapolation from 2P2F and 3P1F regions
-<
+         ctrl.faketype=="2P2F" || // plot 2P2F region, and if requested write out ntuple which gives extrapolation to SR (ntuple would then be read by "SR", above)
-<
+         ctrl.faketype=="3P1F");  // plot 3P1F region (including the extrapolation to there from the 2P2F region), and if requested write out ntuple with SR prediction
->
+  assert(ctrl.faketype=="SR" || ctrl.faketype=="2P2F" || ctrl.faketype=="3P1F");
+  if(ctrl.heavyFlavor) assert(ctrl.faketype=="2P2F");
+  if(ctrl.ssof) assert(ctrl.faketype=="2P2F");
-<
+  can = new TCanvas("can","can");
->
+  can = new TCanvas("can","can",700,500);
-<
+  FR_struct fr2011 = initFRs(ctrl.mufakefile2011,ctrl.elefakefile2011);
-<
+  FR_struct fr2012 = initFRs(ctrl.mufakefile2012,ctrl.elefakefile2012);
->
+  FR_struct fr2011(ctrl.mufakefile2011,ctrl.elefakefile2011);
->
+  FR_struct fr2012(ctrl.mufakefile2012,ctrl.elefakefile2012);
+  TString cmsswpath(CMSSW_BASE + TString("/src"));
+  string xspath = (string(cmsswpath)+"/MitPhysics/data/xs.dat");
+  TString ntupledir(""),label(""),plotLabel(""),jsonFile("");
+  int puTarget;
+  readConfigFile(ctrl.config, ntupledir, label, plotLabel, jsonFile, puTarget, samplev, &ctrl, init_hists);
-+
+  map<TString,double> ssofRatios;
-+
+  if(ctrl.ssof && !ctrl.writessofratio && !ctrl.plotWholeSample) ssofRatios = initSsofRatios(samplev);
-+
+  Angles angles; // NOTE: this is not the same as fs->angles -- the one in fs is filled with whichever choice of leptons was made in applyZPlusX
+  KinematicsStruct kine; // NOTE: this is not the same as fs->kine -- the one in fs is filled with whichever choice of leptons was made in applyZPlusX
+  JetInfoStruct ji;
+    AngleTuple *jettuple;
+    if(makeJetTuple) {
+      jettuple = new AngleTuple( (const char*)(TString(cs->name)+"_pass.root").Data(), (const char*)"zznt");
-+
+      jettuple->makeAngleBranch(angles);
+      jettuple->makeKinematicsBranch(kine);
+      jettuple->makeJetInfoBranch(ji);
+      const char *str1 = "evt";
+    resetCutVect(cutvec);
+    for(unsigned ifs=0; ifs<(cs->fsv).size(); ifs++) {
+      filestuff *fs = (cs->fsv)[ifs];
-<
+      if(makeFakeTuples && fs->dataset_!="fakes") { // write a copy of the input trees for events that will be used for the SR fake prediction
->
+      if(ctrl.makeFakeTuples && fs->dataset_!="fakes") { // write a copy of the input trees for events that will be used for the SR fake prediction
+        TString fakefilename(fs->fname_);
-<
+        fakefilename.ReplaceAll(".root","-fakes.root");
->
+        fakefilename.ReplaceAll(".root","-"+ctrl.faketype+(ctrl.ssof ? "-ssof" : "")+"-fakes.root");
->
+        if(ctrl.hiStatFakes) fakefilename.ReplaceAll("-fakes.root","-fakes-histat.root");
+        fs->makeOutputFile(fakefilename);
+      cout << "\t" << fs->fname_; cout.flush();
-<
+      FR_struct *fr = (fs->era_==2011) ? &fr2011 : &fr2012;
-<
+      unsigned nDuplSkipped=0;
->
+      unsigned nDuplSkipped=0,nBadZs=0;
->
+      watch.Start();
+      for(unsigned ientry=0; ientry<fs->getentries("FOtree"); ientry++) {
+        fs->getentry(ientry,"FOs","FOtree");
+        fs->getentry(ientry,"Zleptons","FOtree");
+        fs->getentry(ientry,"info","zznt");
-<
+        // if(eoiAdish.has(fs->info->run, fs->info->evt) &&
-<
+        //    !eoiMe.has(fs->info->run, fs->info->evt)) {
-<
+        //   cout << "found an event: " << fs->info->run << " " << fs->info->evt << endl;
-<
+        //   ctrl.debug = true;
-<
+        // } else
-<
+        //   ctrl.debug = false;
-<
->
+        if(!(ientry%15000)) {
->
+          watch.Stop();
->
+          cout << "\ntime: " << watch.RealTime() << "  (entry: " << ientry << " )" << endl;
->
+          watch.Start();
->
+        }
+        cutvec["start"] = fs->total_entries_;
+        if(fs->isdata_) {// && !(fs->dataset_=="fakes")) {
+          // WARNING: when I write the fakeTuple from here, I fill an event *twice* if it has two fake z2 candidates, so event numbers are duplicated.
+          // in plotH4l I *don't* want to remove these, but here I *do*, because in this file I again loop over all z2 candidates
+          setMinMaxRun(fs->info->run, minRun, maxRun);
+          bool dupl = takeCareOfDuplicateEvents(fs->info->run, fs->info->evt, runEvtv, nDuplSkipped);
-<
+          if(dupl) { if(ctrl.debug) cout << " failing dupl" << endl; continue; }
-<
+          if(!fs->rlrm_.HasRunLumi(pair<unsigned,unsigned> (fs->info->run, fs->info->lumi))) { if(ctrl.debug) cout << " failing rlrm" << endl; continue; }
->
+          if(dupl) continue;
->
+          if(!fs->rlrm_.HasRunLumi(pair<unsigned,unsigned> (fs->info->run, fs->info->lumi))) continue;
+        cutvec["rlrmAndDupl"] += 1;
-–
+        double wgt=1;
-–
+        if(!fs->isdata_) {
-–
+          double xsWgt = fs->lumi_*xstab.Get(fs->dataset_)/fs->total_entries_;
-–
+          double puWgt = getPUWeight(fs->era_,fs->info->npu);
-–
+          if(fs->era_==2012)
-–
+            puWgt = weightTrue2012(fs->info->npu);
-–
+          wgt = xsWgt*puWgt;
-–
+        }
-–
+        if(fs->dataset_=="fakes") { // if we're reading in the fake ntuple that we wrote out earlier, use the event weight that was stored (this should be the same as what we'd get if we calculated it below)
-–
+          fs->getentry(ientry,"weights","zznt");
-–
+          wgt = fs->weights->w;
-–
+        }
-–
+        unsigned npass = fs->passingL->size();
+        unsigned nfail = fs->failingL->size();
+        // look for a z1
+        pair<int,int> best_z_indices;
+        double best_mz1 = findZ1(fs->passingL,best_z_indices,40);
-<
+        if(best_mz1<=0) {
-<
+          cout << "WARNING: best z1 mass: " << best_mz1 << endl;
-<
+        }
-<
+        if(ctrl.heavyFlavor && best_mz1<60) { if(ctrl.debug) cout << " failing mz1" << endl; continue; }
->
+        if(best_mz1<=0) nBadZs++;
->
+        if(best_mz1 < 40 || best_mz1 > 120) continue;
->
+        if(ctrl.heavyFlavor && best_mz1<60) continue;
+        EventData evtdata;
+        evtdata.Z1leptons.push_back((*fs->passingL)[best_z_indices.first]);
+        evtdata.Z1leptons.push_back((*fs->passingL)[best_z_indices.second]);
+        vector<pair<SimpleLepton,SimpleLepton> > z2cands;
+        vector<TString> types;
+        TString signFlavor;
-<
+        if(ctrl.heavyFlavor) signFlavor = "";
-<
+        else if(ctrl.ssof)   signFlavor = "SS-SF";
-<
+        else                 signFlavor = "OS-SF";
->
+        if(ctrl.heavyFlavor)            signFlavor = "";
->
+        else if(ctrl.ssof)              signFlavor = "SS-SF";
->
+        else if(ctrl.hiStatFakes)       signFlavor = "signalExclude";
->
+        else                            signFlavor = "OS-SF";
+        double minMz2 = ctrl.heavyFlavor ? 0 : 12;
+        int hiPtPassingZ2 = findZ2CandidatesPassFail(fs->passingL, z2cands, types, signFlavor, minMz2, "pass",&best_z_indices); // signal region
+        int hiPtFailingZ2 = findZ2CandidatesPassFail(fs->failingL, z2cands, types, signFlavor, minMz2, "fail"); // 2P+2F
+        int hiPtPfZ2      = findZ2CandidatesPassFail(fs->passingL, z2cands, types, signFlavor, minMz2, "PF", &best_z_indices, fs->failingL); // 3P+1F
-<
+        if(ctrl.debug) {
-<
+          cout << "z2 cands: " << endl;
-<
+          for(unsigned iz2=0; iz2<z2cands.size(); iz2++)
-<
+            cout << " ---" << types[iz2] << "'ing z2" << endl;
->
+        double wgt=1;
->
+        if(z2cands.size() > 0) {
->
+          if(!fs->isdata_) {
->
+            fs->getentry(ientry,"weights","zznt");
->
+            double xsWgt = fs->lumi_*xstab.Get(fs->dataset_)/fs->total_entries_;
->
+            double puWgt = getPUWeight(fs->era_, fs->era_==2011 ? "/025/" : "/029/", fs->info->npu);
->
+            // if(fs->era_==2012)
->
+            //   puWgt = hpu_2012_me->GetBinContent(hpu_2012_me->FindBin(fs->info->npu));
->
+            // else if(fs->era_==2011)
->
+            //   puWgt = hpu_2011_me->GetBinContent(hpu_2011_me->FindBin(fs->info->npu));
->
+            // else assert(0);
->
+            wgt = xsWgt * puWgt;//* fs->weights->won * fs->weights->woff;
->
+            if(wgt!=wgt) cout << "xsWgt: " << xsWgt << " puWgt: " << puWgt << endl;
->
+          }
->
+          if(fs->dataset_=="fakes") { // if we're reading in the fake ntuple that we wrote out earlier, use the event weight that was stored (this should be the same as what we'd get if we calculated it below)
->
+            // fs->getentry(ientry,"weights","zznt");
->
+            wgt = fs->weights->w;
->
+          }
->
+          if(ctrl.writessofratio) incrementSsofCounters(ctrl, fs, minMz2, &best_z_indices, wgt);
+        for(unsigned iz2=0; iz2<z2cands.size(); iz2++) {
+          evtdata.Z2leptons.clear();
+          evtdata.Z2leptons.push_back(lep3);
+          evtdata.Z2leptons.push_back(lep4);
-+
+          fillAngles(evtdata,angles);
+          fillKinematics(evtdata,kine);
+          lepFrs fwgts;
-<
+          fillLeptonFakeWeights(fwgts, fr, lep3, lep4);
->
+          fillLeptonFakeWeights(fwgts, (fs->era_==2011) ? &fr2011 : &fr2012, lep3, lep4);
+          TString channel = getChannel(lep1,lep2,lep3,lep4);
+          // bool hltPass = fs->isdata_ ? passHlt(ctrl,ti,fs->info,lep1.scID,lep2.scID,lep3.scID,lep4.scID) : true;
+          if(ctrl.heavyFlavor) {
-<
+            if(fabs(lep1.dz) > 0.01 || fabs(lep2.dz) > 0.01)            { cout << " failing XXX" << endl; continue; }
->
+            if(fabs(lep1.dz) > 0.01 || fabs(lep2.dz) > 0.01)            continue;
+            if( !(fabs(lep3.ip3dSig) > 2 && fabs(lep4.ip3dSig) > 8) &&
-<
+                !(fabs(lep4.ip3dSig) > 2 && fabs(lep3.ip3dSig) > 8)  )  { cout << " failing XXX" << endl; continue; }
->
+                !(fabs(lep4.ip3dSig) > 2 && fabs(lep3.ip3dSig) > 8)  )  continue;
+            if( !(fabs(lep3.d0) > .002 && fabs(lep4.d0) > .01) &&
-<
+                !(fabs(lep4.d0) > .002 && fabs(lep3.d0) > .01)  )       { cout << " failing XXX" << endl; continue; }
-<
+            if(fabs(lep3.d0) > 0.1 || fabs(lep4.d0) > 0.1)              { cout << " failing XXX" << endl; continue; }
-<
+            if(fabs(lep3.dz) > 0.2 || fabs(lep4.dz) > 0.2)              { cout << " failing XXX" << endl; continue; }
->
+                !(fabs(lep4.d0) > .002 && fabs(lep3.d0) > .01)  )       continue;
->
+            if(fabs(lep3.d0) > 0.1 || fabs(lep4.d0) > 0.1)              continue;
->
+            if(fabs(lep3.dz) > 0.2 || fabs(lep4.dz) > 0.2)              continue;
+          } else {
-<
+            if(fabs(lep3.ip3dSig) >= 4 || fabs(lep4.ip3dSig) >= 4)      { if(ctrl.debug) cout << " failing sip" << endl; continue; } // |sip| < 4 is applied to lep1, lep2 in ZPlusX
-<
+            if(!leptonDrReqs(lep1,lep2,lep3,lep4))                      { if(ctrl.debug) cout << " failing dr" << endl; continue; }
-<
+            if(kine.m4l < 100)                                          { if(ctrl.debug) cout << " failing m4l 100" << endl; continue; }
-<
+            if(!resonanceKilling(lep1,lep2,lep3,lep4))                  { if(ctrl.debug) cout << " failing reskill" << endl; continue; }
->
+            if(fabs(lep3.ip3dSig) >= 4 || fabs(lep4.ip3dSig) >= 4)      continue; // |sip| < 4 is applied to lep1, lep2 in ZPlusX
->
+            if(!leptonDrReqs(lep1,lep2,lep3,lep4))                      continue;
->
+            if(kine.m4l < 100)                                          continue;
->
+            if(!resonanceKilling(lep1,lep2,lep3,lep4))                  continue;
-<
+          if(!finalLeptonPtReqs(lep1,lep2,lep3,lep4))                   { if(ctrl.debug) cout << " failing lepton pt " << endl; continue; }
->
+          if(!finalLeptonPtReqs(lep1,lep2,lep3,lep4))                   continue;
+          cutvec["4lsele"] += 1;
+            if(fs->jets->size() > 1) cutvec["twoJets"] += 1;
+            bool hazJetz = findGoodJets(goodJets, fs, lep1, lep2, lep3, lep4);
+            if(goodJets.size() > 1) cutvec["twoJetsAfter"] += 1;
-<
+            if(!hazJetz) { cout << " failing XXX" << endl; continue; }
->
+            if(!hazJetz) continue;
+            ControlFlags ctrlTmp;
+            fillJetInfo(goodJets,ji,ctrlTmp);
+            if(makeJetTuple) {
+            fusion.setValues(ji,kine);
-<
->
+          cutvec["filling"] += 1;
+          if(types[iz2]=="pass" || ((fs->dataset_=="fakes") && (types[iz2]=="fail"))) { // for the fakes ntuple, we plot the "failing" z2s as passing ones, i.e. put them in hObs, but with the proper weight
+            // llll events
-<
+            if((fs->dataset_!="fakes") && (iz2!=hiPtPassingZ2)) { cout << " failing XXX" << endl; continue; } // for non-fakes, we only want the  highest-pt "passing" z2 (for the fakes these are by design failing z2s, and we want all of 'em)
->
+            if((fs->dataset_!="fakes") && (iz2!=hiPtPassingZ2)) continue; // for non-fakes, we only want the  highest-pt "passing" z2 (for the fakes these are by design failing z2s, and we want all of 'em)
+            cutvec["fillPass"] += 1;
-<
+            fillAllHists( ctrl, cs, channel, "obs", fs, kine, lep1, lep2, lep3, lep4, wgt);
-<
+            if(doJets) fillAllJetHists( ctrl, cs, channel, "obs", fs, kine, lep1, lep2, lep3, lep4, fusion, goodJets, ji, wgt);
->
+            fillAllHists( ctrl, cs, channel, "obs", fs, kine, angles, lep1, lep2, lep3, lep4, wgt);
->
+            if(doJets) fillAllJetHists( ctrl, cs, channel, "obs", fs, kine, lep1, lep2, lep3, lep4, /*fusion,*/ /*goodJets, */ji, wgt);
+          } else if(types[iz2]=="fail") { // lljj events
-<
+            if(iz2!=hiPtFailingZ2 && ctrl.plotWholeSample) { cout << " failing XXX" << endl; continue; }
->
+            if(iz2!=hiPtFailingZ2 && ctrl.plotWholeSample) continue;
+            double centr=0,lo=0,hi=0;
+            if(ctrl.faketype=="2P2F") {
+              if(ctrl.plotWholeSample) { // make plots without weighting with the FR
+                centr = wgt*fwgts.fwgt_3*fwgts.fwgt_4;
+                lo    = wgt*fwgts.fwgt_lo_3*fwgts.fwgt_lo_4;
+                hi    = wgt*fwgts.fwgt_hi_3*fwgts.fwgt_hi_4;
-+
+                if(ctrl.ssof && !ctrl.writessofratio && !ctrl.plotWholeSample) {
-+
+                  double ssof_wgt = getSsofWgt(channel,ssofRatios);
-+
+                  centr *= ssof_wgt;
-+
+                  lo    *= ssof_wgt;
-+
+                  hi    *= ssof_wgt;
-+
+                }
+            } else if(ctrl.faketype=="3P1F") { // apply the weights for extrapolation from 2P2F region to 3P1F region
+              centr = wgt*(fwgts.fwgt_3    + fwgts.fwgt_4);
+              hi    = wgt*(fwgts.fwgt_hi_3 + fwgts.fwgt_hi_4);
-<
+            // runfile << fixed << setprecision(0) << fs->info->run << " " << fixed << setprecision(0) << fs->info->evt << endl;
-<
-<
+            if(makeFakeTuples && fs->dataset_!="fakes") {
->
+            if(ctrl.makeFakeTuples && fs->dataset_!="fakes") {
+              double fillwgt;
+              bool dofill=true;
+              if(ctrl.faketype=="2P2F") {
-<
+                assert(fs->isdata_);
->
+                // assert(fs->isdata_);
+                fillwgt = centr;
+              } else if(ctrl.faketype=="3P1F") {
+                fillwgt = - wgt*(fwgts.fwgt_3*fwgts.fwgt_4 + fwgts.fwgt_4*fwgts.fwgt_3);
+              fillFakeTuple(cs,fs,evtdata,ientry,fillwgt);
+            cutvec["fillFail"] += 1;
-<
+            fillAllHists( ctrl, cs, channel, "pred", fs, kine, lep1, lep2, lep3, lep4, centr, lo, hi);
-<
+            if(doJets) fillAllJetHists( ctrl, cs, channel, "pred", fs, kine, lep1, lep2, lep3, lep4, fusion, goodJets, ji, centr, lo, hi);
->
+            fillAllHists( ctrl, cs, channel, "pred", fs, kine, angles, lep1, lep2, lep3, lep4, centr, lo, hi);
->
+            if(doJets) fillAllJetHists( ctrl, cs, channel, "pred", fs, kine, lep1, lep2, lep3, lep4, /*fusion,*/ /*goodJets,*/ ji, centr, lo, hi);
+          } else if(types[iz2]=="PF") { // lllj events
+            if(fs->dataset_!="fakes" && iz2!=hiPtPfZ2) continue; // shouldn't have this when I extrapolate to signal region
-<
+            if(makeFakeTuples && fs->dataset_!="fakes" && ctrl.faketype=="3P1F") {
->
+            if(ctrl.makeFakeTuples && fs->dataset_!="fakes" && ctrl.faketype=="3P1F") {
+              assert(lep3.isLoose || lep4.isLoose);
+              assert(!(lep3.isLoose && lep4.isLoose));
+              double fwgt = lep3.isLoose ? fwgts.fwgt_4 : fwgts.fwgt_3;
+            cutvec["fillPF"] += 1;
-<
+            fillAllHists( ctrl, cs, channel, "PF", fs, kine, lep1, lep2, lep3, lep4, wgt);
-<
+            if(doJets) fillAllJetHists( ctrl, cs, channel, "PF", fs, kine, lep1, lep2, lep3, lep4, fusion, goodJets, ji, wgt);
->
+            fillAllHists( ctrl, cs, channel, "PF", fs, kine, angles, lep1, lep2, lep3, lep4, wgt);
->
+            if(doJets) fillAllJetHists( ctrl, cs, channel, "PF", fs, kine, lep1, lep2, lep3, lep4, /*fusion,*/ /*goodJets,*/ ji, wgt);
+          } else assert(0);
-<
+      cout << "\t\t" << setw(7) << nDuplSkipped << " duplicate events skipped" << endl;
-<
+      if(makeFakeTuples) {
->
+      if(ctrl.writessofratio) fs->writeSsofRatios();
->
+      cout << "\t\t" << setw(7) << nDuplSkipped << " duplicate events skipped (" << nBadZs << " bad Zs)" << endl;
->
+      if(ctrl.makeFakeTuples) {
+        fs->outFotuple->getFile()->cd();
+        fs->outFotuple->getTree()->Write();
+        fs->outtuple->WriteClose();
-+
+      // fs->del();
+    for(unsigned icut=0; icut<cutstrs.size(); icut++) {
+      cout << setw(22) << cutstrs[icut] << setw(22) << cutvec[cutstrs[icut]] << endl;
+  vector<TString> typev;
+  typev.push_back("all");
-<
+  bool plotAllCats=false;
-<
+  if(plotAllCats) {
->
+  if(ctrl.plotAllCats) {
+    typev.push_back("4e");
+    typev.push_back("4m");
+    typev.push_back("2e2m");
-+
+  ofstream txtOutFile("fakes.txt",ios_base::app);
-+
+  txtOutFile << "---> " << ctrl.faketype << " " << plotLabel << endl;
+  // loop over channels
+  for(unsigned itype=0; itype<typev.size(); itype++) {
+    TString type(typev[itype]);
-<
+    gSystem->mkdir(ctrl.outdir+"/"+plotLabel+"/"+type,true);
-<
+    TFile runHistFile(ctrl.outdir+"/"+plotLabel+"/"+type+"/runs.root","recreate");
->
+    TString fullOutDir(ctrl.outdir+"/"+plotLabel+"/"+type);
->
+    gSystem->mkdir(fullOutDir,true);
->
+    TFile runHistFile(fullOutDir+"/runs.root","recreate");
+    // loop over variables
+    map<TString,TH1D*>::iterator it_v;
+    for(it_v=(*(samplev[0]->hists)[type+"_obs"]).begin(); it_v!=(*(samplev[0]->hists)[type+"_obs"]).end(); it_v++) {
+      TString var((*it_v).first);
-<
+      CPlot cplot(var,"",(*(samplev[0]->hists)[type+"_obs"])[var]->GetXaxis()->GetTitle(),"events",ctrl.outdir+"/"+plotLabel+"/"+type+"/plots");
->
+      CPlot cplot(var,TString(""),(*(samplev[0]->hists)[type+"_obs"])[var]->GetXaxis()->GetTitle(),TString("events"),TString(fullOutDir+"/plots"));
+      double fakeCounter=0;
+      double ymax=0;
+      // loop over csamples
+          scaleHists(hs);
+        if(cs->isdata && !(cs->name=="fakes")) {
-<
+          if(ctrl.faketype=="SR" || (ctrl.heavyFlavor && !ctrl.plotWholeSample)) {
-<
+            cplot.AddHist1D(hs["obs"],cs->legend+": "+integral_str(hs["obs"],0),"E");
->
+          // add observation (llll events) in an 'E' hists
->
+          if(ctrl.faketype=="SR" || (ctrl.heavyFlavor && !ctrl.plotWholeSample) || (ctrl.ssof && !ctrl.plotWholeSample)) {
->
+            cplot.AddHist1D(hs["obs"],cs->legend+": "+integral_str(hs["obs"],0,dummyStr),"E");
+            ymax = max(ymax,hs["obs"]->GetMaximum());
+          if(ctrl.faketype=="2P2F") {
+            if(ctrl.plotWholeSample) {
-<
+              cplot.AddHist1D(hs["pred"],cs->legend+": "+integral_str(hs["pred"],0),"E");
->
+              // add 'pred' (lljj events with weight *1*) as an 'E' hist
->
+              cplot.AddHist1D(hs["pred"],cs->legend+": "+integral_str(hs["pred"],0,dummyStr),"E");
+              ymax = max(ymax,hs["pred"]->GetMaximum());
+            } else {
-<
+              cplot.AddHist1D(hs["pred"],   "FR predic: "+integral_str(hs["pred"],3),   "hist",kRed);
-<
+              cplot.AddHist1D(hs["pred_lo"],  "stat lo: "+integral_str(hs["pred_lo"],2),"hist",kRed,kDashed);
-<
+              cplot.AddHist1D(hs["pred_hi"],  "stat hi: "+integral_str(hs["pred_hi"],2),"hist",kRed,kDashed);
->
+              // add 'pred' (lljj events with weight w1*w2) as red-line hists
->
+              cplot.AddHist1D(hs["pred"],    TString("FR predic: "+integral_str(hs["pred"],3,dummyStr)),   "hist",kRed);
->
+              cplot.AddHist1D(hs["pred_lo"], TString(  "stat lo: "+integral_str(hs["pred_lo"],2,dummyStr)),"hist",kRed,kDashed);
->
+              cplot.AddHist1D(hs["pred_hi"], TString(  "stat hi: "+integral_str(hs["pred_hi"],2,dummyStr)),"hist",kRed,kDashed);
->
+              if(var=="m4l") {
->
+                txtOutFile << "  " << type << ": " << integral_str(hs["pred"],3) << " (" << integral_str(hs["pred_lo"],2) << "," << integral_str(hs["pred_hi"],2) << ")" << endl;
->
+              }
+              ymax = max(ymax,histMax(hs["pred"],hs["pred_lo"],hs["pred_hi"]));
+          if(ctrl.faketype=="3P1F") {
-<
+            double ratio = 0;//integrateHist(hs["PF"]) / integrateHist(hs["pred"]);
-<
+            cplot.AddHist1D(hs["PF"],cs->legend+": "+integral_str(hs["PF"],0)+", ratio: "+f_to_a(ratio),"E",cs->color);
->
+            // add 'PF' (lllj events with weight 1) as an 'E' hist
->
+            // double ratio = 0;//integrateHist(hs["PF"]) / integrateHist(hs["pred"]);
->
+            // cplot.AddHist1D(hs["PF"],TString(cs->legend+": "+integral_str(hs["PF"],0,dummyStr)+", ratio: "+f_to_a(ratio)),"E",cs->color);
->
+            cplot.AddHist1D(hs["PF"],TString(cs->legend+": "+integral_str(hs["PF"],0,dummyStr)),"E",cs->color);
+            if(!ctrl.plotWholeSample) {
-<
+              cplot.AddToStack(hs["pred"],cs->legend+"(2P2F extrap.): "+integral_str(hs["pred"],1), kRed);
->
+              // add 'pred' (lljj events with weight ~w1) as a red stacked hist
->
+              cplot.AddToStack(hs["pred"],TString(cs->legend+"(2P2F extrap.): "+integral_str(hs["pred"],1,dummyStr)), 823, -1);
+              ymax = max(ymax,histMax(hs["PF"],cplot.GetStack(),hs["pred"]));
+            } else {
+              ymax = max(ymax,hs["PF"]->GetMaximum());
+          if(ctrl.faketype=="3P1F")   mchist = hs["PF"];
+          assert(mchist);
+          bool stackMc = true;
-<
+          if(stackMc) cplot.AddToStack(mchist,cs->legend+": "+integral_str(mchist,2),cs->color);
-<
+          else        cplot.AddHist1D(mchist,cs->legend +": "+integral_str(mchist,3),"Ehist",cs->color);//cs->color);
->
+          if(stackMc) cplot.AddToStack(mchist,TString(cs->legend+": "+integral_str(mchist,2,dummyStr)),cs->color,-1);
->
+          else        cplot.AddHist1D(mchist,TString(cs->legend +": "+integral_str(mchist,3,dummyStr)),"Ehist",cs->color);//cs->color);
+          if(cs->name=="zz") fakeCounter -= integrateHist(hs["PF"]);
+          ymax = max(ymax,histMax(mchist,cplot.GetStack()));
+        if(!var.Contains("fusionMVA"))
+          cplot.SetYRange(0,1.2*ymax);
-<
+        if(ctrl.faketype=="3P1F") cplot.AddTextBox("diff: "+f_to_a(fakeCounter),.7,.5,.8,.6);
->
+        if(!dummyStr && ctrl.faketype=="3P1F") cplot.AddTextBox(TString("diff: "+f_to_a(fakeCounter)),.7,.5,.8,.6);
-+
+      TString txt("#bf{CMS Preliminary           L = 5 ( + 12 ) fb^{-1} at #sqrt{s} = 7 ( + 8 ) TeV}"); // #sqrt{s} = 7 TeV, L = 5 fb^{-1}; #sqrt{s} = 8 TeV, L = 12 fb^{-1}");
-+
+      cplot.AddTextBox(txt,0.16,0.95,0.85,0.99,0);
+      cplot.Draw(can,true,"png");
+      if(var.Contains("fusionMVA")) {
+        cplot.SetLogy();
+    runHistFile.Close();
-<
+    makeHTML(ctrl,type,plotLabel);
->
+    makeHTML(ctrl,type,plotLabel,fullOutDir);
-+
+  txtOutFile.close();
-+
-+
+  fr2011.FR_struct::~FR_struct();
-+
+  fr2012.FR_struct::~FR_struct();
+//----------------------------------------------------------------------------------------
+map<TString,map<TString,TH1D*>* > init_hists(FOFlags &ctrl, TString str)
+  map<TString,TH1D*> *h_all_PF_hi       = new map<TString,TH1D*>;      hists["all_PF_hi"]       = h_all_PF_hi;
+  map<TString,map<TString,TH1D*>* >::iterator it_h;
-<
+  for(it_h=hists.begin(); it_h!=hists.end(); it_h++) {
-<
+    (*((*it_h).second))["npv"]          = new TH1D(TString("npv") +"_"+(*it_h).first+str,";#bf{N PV};",         50,-0.5,49.5);      (*((*it_h).second))["npv"]->Sumw2();
-<
+    (*((*it_h).second))["run"]          = new TH1D(TString("run") +"_"+(*it_h).first+str,";#bf{run};",          50,160800,196535);   (*((*it_h).second))["run"]->Sumw2();
-<
+    (*((*it_h).second))["mZ1"]          = new TH1D(TString("mZ1") +"_"+(*it_h).first+str,";#bf{mZ1 [GeV]};",    20,40,120);          (*((*it_h).second))["mZ1"]->Sumw2();
-<
+    (*((*it_h).second))["mZ2"]          = new TH1D(TString("mZ2") +"_"+(*it_h).first+str,";#bf{mZ2 [GeV]};",    30,0,115);           (*((*it_h).second))["mZ2"]->Sumw2();
-<
+    (*((*it_h).second))["mZ2_lo"]       = new TH1D(TString("mZ2_lo") +"_"+(*it_h).first+str,";#bf{mZ2 [GeV]};", 30,0,10);           (*((*it_h).second))["mZ2_lo"]->Sumw2();
-<
+    (*((*it_h).second))["m4l_lo"]       = new TH1D(TString("m4l_lo") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};", 35,70,120);          (*((*it_h).second))["m4l_lo"]->Sumw2();
-<
+    (*((*it_h).second))["m4l_med"]      = new TH1D(TString("m4l_med") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",35,100,180);          (*((*it_h).second))["m4l_med"]->Sumw2();
-<
+    (*((*it_h).second))["m4l"]          = new TH1D(TString("m4l") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",    35,100,600);          (*((*it_h).second))["m4l"]->Sumw2();
-<
+    (*((*it_h).second))["Z1pt"]         = new TH1D(TString("Z1pt") +"_"+(*it_h).first+str,";#bf{Z1pt [GeV]};",  20,0,200);           (*((*it_h).second))["Z1pt"]->Sumw2();
-<
+    (*((*it_h).second))["ZZpt"]         = new TH1D(TString("ZZpt") +"_"+(*it_h).first+str,";#bf{ZZpt [GeV]};",  30,0,200);           (*((*it_h).second))["ZZpt"]->Sumw2();
-<
+    (*((*it_h).second))["met"]          = new TH1D(TString("met") +"_"+(*it_h).first+str,";#bf{met [GeV]};",    30,0,100);           (*((*it_h).second))["met"]->Sumw2();
-<
-<
+    (*((*it_h).second))["pt_l1"]        = new TH1D(TString("pt_l1") +"_"+(*it_h).first+str,";#bf{l1 p_{T} [GeV]};",     37,0,120);  (*((*it_h).second))["pt_l1"]->Sumw2();
-<
+    (*((*it_h).second))["pt_l2"]        = new TH1D(TString("pt_l2") +"_"+(*it_h).first+str,";#bf{l2 p_{T} [GeV]};",     37,0,65);  (*((*it_h).second))["pt_l2"]->Sumw2();
-<
+    (*((*it_h).second))["pt_l3"]        = new TH1D(TString("pt_l3") +"_"+(*it_h).first+str,";#bf{l3 p_{T} [GeV]};",     37,0,75);  (*((*it_h).second))["pt_l3"]->Sumw2();
-<
+    (*((*it_h).second))["pt_l4"]        = new TH1D(TString("pt_l4") +"_"+(*it_h).first+str,";#bf{l4 p_{T} [GeV]};",     37,0,50);  (*((*it_h).second))["pt_l4"]->Sumw2();
-<
-<
+    (*((*it_h).second))["eta_l1"]       = new TH1D(TString("eta_l1")+"_"+(*it_h).first+str,";#bf{l1 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l1"]->Sumw2();
-<
+    (*((*it_h).second))["eta_l2"]       = new TH1D(TString("eta_l2")+"_"+(*it_h).first+str,";#bf{l2 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l2"]->Sumw2();
-<
+    (*((*it_h).second))["eta_l3"]       = new TH1D(TString("eta_l3")+"_"+(*it_h).first+str,";#bf{l3 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l3"]->Sumw2();
-<
+    (*((*it_h).second))["eta_l4"]       = new TH1D(TString("eta_l4")+"_"+(*it_h).first+str,";#bf{l4 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l4"]->Sumw2();
-<
-<
+    (*((*it_h).second))["ip3ds_l1"]     = new TH1D(TString("ip3ds_l1") +"_"+(*it_h).first+str,";#bf{ip3ds l1};",37,-3,3);       (*((*it_h).second))["ip3ds_l1"]->Sumw2();
-<
+    (*((*it_h).second))["ip3ds_l2"]     = new TH1D(TString("ip3ds_l2") +"_"+(*it_h).first+str,";#bf{ip3ds l2};",37,-3,3);       (*((*it_h).second))["ip3ds_l2"]->Sumw2();
-<
+    (*((*it_h).second))["ip3ds_l3_lo"]  = new TH1D(TString("ip3ds_l3_lo") +"_"+(*it_h).first+str,";#bf{ip3ds l3_lo};",25,-6,6); (*((*it_h).second))["ip3ds_l3_lo"]->Sumw2();
-<
+    (*((*it_h).second))["ip3ds_l4_lo"]  = new TH1D(TString("ip3ds_l4_lo") +"_"+(*it_h).first+str,";#bf{ip3ds l4_lo};",25,-6,6); (*((*it_h).second))["ip3ds_l4_lo"]->Sumw2();
-<
+    (*((*it_h).second))["ip3ds_l3"]     = new TH1D(TString("ip3ds_l3") +"_"+(*it_h).first+str,";#bf{ip3ds l3};",55,-40,40);     (*((*it_h).second))["ip3ds_l3"]->Sumw2();
-<
+    (*((*it_h).second))["ip3ds_l4"]     = new TH1D(TString("ip3ds_l4") +"_"+(*it_h).first+str,";#bf{ip3ds l4};",55,-40,40);     (*((*it_h).second))["ip3ds_l4"]->Sumw2();
-<
-<
+    (*((*it_h).second))["d0_l1"]        = new TH1D(TString("d0_l1") +"_"+(*it_h).first+str,";#bf{d0 l1};",50,-.008,.008);      (*((*it_h).second))["d0_l1"]->Sumw2();
-<
+    (*((*it_h).second))["d0_l2"]        = new TH1D(TString("d0_l2") +"_"+(*it_h).first+str,";#bf{d0 l2};",50,-.008,.008);      (*((*it_h).second))["d0_l2"]->Sumw2();
-<
+    (*((*it_h).second))["d0_l3_lo"]     = new TH1D(TString("d0_l3_lo") +"_"+(*it_h).first+str,";#bf{d0 l3_lo};",25,-.03,.03);  (*((*it_h).second))["d0_l3_lo"]->Sumw2();
-<
+    (*((*it_h).second))["d0_l4_lo"]     = new TH1D(TString("d0_l4_lo") +"_"+(*it_h).first+str,";#bf{d0 l4_lo};",25,-.03,.03);  (*((*it_h).second))["d0_l4_lo"]->Sumw2();
-<
+    (*((*it_h).second))["d0_l3"]        = new TH1D(TString("d0_l3") +"_"+(*it_h).first+str,";#bf{d0 l3};",50,-.1,.1);          (*((*it_h).second))["d0_l3"]->Sumw2();
-<
+    (*((*it_h).second))["d0_l4"]        = new TH1D(TString("d0_l4") +"_"+(*it_h).first+str,";#bf{d0 l4};",50,-.1,.1);          (*((*it_h).second))["d0_l4"]->Sumw2();
-<
-<
+    (*((*it_h).second))["dz_l1"]        = new TH1D(TString("dz_l1") +"_"+(*it_h).first+str,";#bf{dz l1};",50,-.01,.01);   (*((*it_h).second))["dz_l1"]->Sumw2();
-<
+    (*((*it_h).second))["dz_l2"]        = new TH1D(TString("dz_l2") +"_"+(*it_h).first+str,";#bf{dz l2};",50,-.01,.01);   (*((*it_h).second))["dz_l2"]->Sumw2();
-<
+    (*((*it_h).second))["dz_l3"]        = new TH1D(TString("dz_l3") +"_"+(*it_h).first+str,";#bf{dz l3};",50,-.2,.2);   (*((*it_h).second))["dz_l3"]->Sumw2();
-<
+    (*((*it_h).second))["dz_l4"]        = new TH1D(TString("dz_l4") +"_"+(*it_h).first+str,";#bf{dz l4};",50,-.2,.2);   (*((*it_h).second))["dz_l4"]->Sumw2();
-<
-<
+    (*((*it_h).second))["dR_l3l4_lo"]   = new TH1D(TString("dR_l3l4_lo") +"_"+(*it_h).first+str,";#bf{#Delta R l3,l4};",75,0,.2);    (*((*it_h).second))["dR_l3l4_lo"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l3l4"]      = new TH1D(TString("dR_l3l4") +"_"+(*it_h).first+str,";#bf{#Delta R l3,l4};",75,0,6);        (*((*it_h).second))["dR_l3l4"]->Sumw2();
-<
-<
+    // jets
-<
+    (*((*it_h).second))["pt_allJet"]    = new TH1D(TString("pt_allJet") +"_"+(*it_h).first+str,";#bf{jet 1 p_{T} [GeV]};",37,0,400);  (*((*it_h).second))["pt_allJet"]->Sumw2();
-<
+    (*((*it_h).second))["eta_allJet"]   = new TH1D(TString("eta_allJet") +"_"+(*it_h).first+str,";#bf{jet 1 p_{T} [GeV]};",37,-4.8,4.8);  (*((*it_h).second))["eta_allJet"]->Sumw2();
-<
+    (*((*it_h).second))["pt_j1"]        = new TH1D(TString("pt_j1") +"_"+(*it_h).first+str,";#bf{jet 1 p_{T} [GeV]};",  37,0,400);  (*((*it_h).second))["pt_j1"]->Sumw2();
-<
+    (*((*it_h).second))["pt_j2"]        = new TH1D(TString("pt_j2") +"_"+(*it_h).first+str,";#bf{jet 2 p_{T} [GeV]};",  37,0,125);  (*((*it_h).second))["pt_j2"]->Sumw2();
-<
+    (*((*it_h).second))["pt_j3"]        = new TH1D(TString("pt_j3") +"_"+(*it_h).first+str,";#bf{jet 3 p_{T} [GeV]};",  37,0,125);  (*((*it_h).second))["pt_j3"]->Sumw2();
-<
+    (*((*it_h).second))["pt_j4"]        = new TH1D(TString("pt_j4") +"_"+(*it_h).first+str,";#bf{jet 4 p_{T} [GeV]};",  37,0,125);  (*((*it_h).second))["pt_j4"]->Sumw2();
-<
+    (*((*it_h).second))["eta_j1"]       = new TH1D(TString("eta_j1")+"_"+(*it_h).first+str,";#bf{jet 1 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j1"]->Sumw2();
-<
+    (*((*it_h).second))["eta_j2"]       = new TH1D(TString("eta_j2")+"_"+(*it_h).first+str,";#bf{jet 2 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j2"]->Sumw2();
-<
+    (*((*it_h).second))["eta_j3"]       = new TH1D(TString("eta_j3")+"_"+(*it_h).first+str,";#bf{jet 3 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j3"]->Sumw2();
-<
+    (*((*it_h).second))["eta_j4"]       = new TH1D(TString("eta_j4")+"_"+(*it_h).first+str,";#bf{jet 4 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j4"]->Sumw2();
-<
+    (*((*it_h).second))["phi_j1"]       = new TH1D(TString("phi_j1")+"_"+(*it_h).first+str,";#bf{jet 1 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j1"]->Sumw2();
-<
+    (*((*it_h).second))["phi_j2"]       = new TH1D(TString("phi_j2")+"_"+(*it_h).first+str,";#bf{jet 2 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j2"]->Sumw2();
-<
+    (*((*it_h).second))["phi_j3"]       = new TH1D(TString("phi_j3")+"_"+(*it_h).first+str,";#bf{jet 3 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j3"]->Sumw2();
-<
+    (*((*it_h).second))["phi_j4"]       = new TH1D(TString("phi_j4")+"_"+(*it_h).first+str,";#bf{jet 4 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j4"]->Sumw2();
-<
+    (*((*it_h).second))["idmva_j1"]     = new TH1D(TString("idmva_j1")+"_"+(*it_h).first+str,";#bf{jet 1 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j1"]->Sumw2();
-<
+    (*((*it_h).second))["idmva_j2"]     = new TH1D(TString("idmva_j2")+"_"+(*it_h).first+str,";#bf{jet 2 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j2"]->Sumw2();
-<
+    (*((*it_h).second))["idmva_j3"]     = new TH1D(TString("idmva_j3")+"_"+(*it_h).first+str,";#bf{jet 3 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j3"]->Sumw2();
-<
+    (*((*it_h).second))["idmva_j4"]     = new TH1D(TString("idmva_j4")+"_"+(*it_h).first+str,";#bf{jet 4 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j4"]->Sumw2();
-<
+    (*((*it_h).second))["mjj"]          = new TH1D(TString("mjj")+"_"+(*it_h).first+str,";#bf{m_{jj}};",                25,10,1550);     (*((*it_h).second))["mjj"]->Sumw2();
-<
+    (*((*it_h).second))["dEta"]         = new TH1D(TString("dEta")+"_"+(*it_h).first+str,";#bf{#Delta #eta};",          25,-6,6);       (*((*it_h).second))["dEta"]->Sumw2();
-<
+    (*((*it_h).second))["etaProd"]      = new TH1D(TString("etaProd")+"_"+(*it_h).first+str,";#bf{#eta_{1} * #eta_{2}};",25,-12,12);      (*((*it_h).second))["etaProd"]->Sumw2();
-<
+    (*((*it_h).second))["njets"]        = new TH1D(TString("njets")+"_"+(*it_h).first+str,";#bf{N Jets};",      7,-0.5,6.5);      (*((*it_h).second))["njets"]->Sumw2();
-<
-<
+    (*((*it_h).second))["dR_l1j1"]      = new TH1D(TString("dR_l1j1") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 1};",75,0,6);        (*((*it_h).second))["dR_l1j1"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l2j1"]      = new TH1D(TString("dR_l2j1") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 1};",75,0,6);        (*((*it_h).second))["dR_l2j1"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l3j1"]      = new TH1D(TString("dR_l3j1") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 1};",75,0,6);        (*((*it_h).second))["dR_l3j1"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l4j1"]      = new TH1D(TString("dR_l4j1") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 1};",75,0,6);        (*((*it_h).second))["dR_l4j1"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l1j2"]      = new TH1D(TString("dR_l1j2") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 2};",75,0,6);        (*((*it_h).second))["dR_l1j2"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l2j2"]      = new TH1D(TString("dR_l2j2") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 2};",75,0,6);        (*((*it_h).second))["dR_l2j2"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l3j2"]      = new TH1D(TString("dR_l3j2") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 2};",75,0,6);        (*((*it_h).second))["dR_l3j2"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l4j2"]      = new TH1D(TString("dR_l4j2") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 2};",75,0,6);        (*((*it_h).second))["dR_l4j2"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l1j3"]      = new TH1D(TString("dR_l1j3") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 3};",75,0,6);        (*((*it_h).second))["dR_l1j3"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l2j3"]      = new TH1D(TString("dR_l2j3") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 3};",75,0,6);        (*((*it_h).second))["dR_l2j3"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l3j3"]      = new TH1D(TString("dR_l3j3") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 3};",75,0,6);        (*((*it_h).second))["dR_l3j3"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l4j3"]      = new TH1D(TString("dR_l4j3") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 3};",75,0,6);        (*((*it_h).second))["dR_l4j3"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l1j4"]      = new TH1D(TString("dR_l1j4") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 4};",75,0,6);        (*((*it_h).second))["dR_l1j4"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l2j4"]      = new TH1D(TString("dR_l2j4") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 4};",75,0,6);        (*((*it_h).second))["dR_l2j4"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l3j4"]      = new TH1D(TString("dR_l3j4") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 4};",75,0,6);        (*((*it_h).second))["dR_l3j4"]->Sumw2();
-<
+    (*((*it_h).second))["dR_l4j4"]      = new TH1D(TString("dR_l4j4") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 4};",75,0,6);        (*((*it_h).second))["dR_l4j4"]->Sumw2();
-<
-<
+    (*((*it_h).second))["fusionMVA"]            = new TH1D(TString("fusionMVA")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,-1,1);      (*((*it_h).second))["fusionMVA"]->Sumw2();
-<
+    (*((*it_h).second))["fusionMVA_lo"]         = new TH1D(TString("fusionMVA_lo")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,-1,-.83);      (*((*it_h).second))["fusionMVA_lo"]->Sumw2();
-<
+    (*((*it_h).second))["fusionMVA_med"]                = new TH1D(TString("fusionMVA_med")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,-.5,1);      (*((*it_h).second))["fusionMVA_med"]->Sumw2();
-<
+    (*((*it_h).second))["fusionMVA_hi"]         = new TH1D(TString("fusionMVA_hi")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,.75,1);      (*((*it_h).second))["fusionMVA_hi"]->Sumw2();
-<
+    (*((*it_h).second))["m4l_jet_lo"]           = new TH1D(TString("m4l_jet_lo") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",     45,100,180);          (*((*it_h).second))["m4l_jet_lo"]->Sumw2();
-<
+    (*((*it_h).second))["m4l_jet"]              = new TH1D(TString("m4l_jet") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",        45,100,600);          (*((*it_h).second))["m4l_jet"]->Sumw2();
-<
+  }
->
+  allocate_hists(hists,str);
->
+  // for(it_h=hists.begin(); it_h!=hists.end(); it_h++) {
->
+  //   (*((*it_h).second))["npv"]               = new TH1D(TString("npv") +"_"+(*it_h).first+str,";#bf{N PV};",         50,-0.5,49.5);      (*((*it_h).second))["npv"]->Sumw2();
->
+  //   (*((*it_h).second))["run"]               = new TH1D(TString("run") +"_"+(*it_h).first+str,";#bf{run};",          50,160800,203744/*196535*/);   (*((*it_h).second))["run"]->Sumw2();
->
+  //   (*((*it_h).second))["mZ1"]               = new TH1D(TString("mZ1") +"_"+(*it_h).first+str,";#bf{mZ1 [GeV]};",    20,40,120);          (*((*it_h).second))["mZ1"]->Sumw2();
->
+  //   (*((*it_h).second))["mZ2"]               = new TH1D(TString("mZ2") +"_"+(*it_h).first+str,";#bf{mZ2 [GeV]};",    30,0,115);           (*((*it_h).second))["mZ2"]->Sumw2();
->
+  //   (*((*it_h).second))["mZ2_lo"]    = new TH1D(TString("mZ2_lo") +"_"+(*it_h).first+str,";#bf{mZ2 [GeV]};", 30,0,10);           (*((*it_h).second))["mZ2_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["m4l_lo"]    = new TH1D(TString("m4l_lo") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};", 35,70,120);          (*((*it_h).second))["m4l_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["m4l_med"]   = new TH1D(TString("m4l_med") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",35,100,180);          (*((*it_h).second))["m4l_med"]->Sumw2();
->
+  //   (*((*it_h).second))["m4l"]               = new TH1D(TString("m4l") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",    35,100,450);          (*((*it_h).second))["m4l"]->Sumw2();
->
+  //   // (*((*it_h).second))["Z1pt"]           = new TH1D(TString("Z1pt") +"_"+(*it_h).first+str,";#bf{Z1pt [GeV]};",  20,0,200);           (*((*it_h).second))["Z1pt"]->Sumw2();
->
+  //   // (*((*it_h).second))["ZZpt"]           = new TH1D(TString("ZZpt") +"_"+(*it_h).first+str,";#bf{ZZpt [GeV]};",  30,0,200);           (*((*it_h).second))["ZZpt"]->Sumw2();
->
+  //   (*((*it_h).second))["met"]               = new TH1D(TString("met") +"_"+(*it_h).first+str,";#bf{met [GeV]};",    30,0,100);           (*((*it_h).second))["met"]->Sumw2();
->
->
+  //   (*((*it_h).second))["pt_l1"]     = new TH1D(TString("pt_l1") +"_"+(*it_h).first+str,";#bf{l1 p_{T} [GeV]};",     37,0,120);  (*((*it_h).second))["pt_l1"]->Sumw2();
->
+  //   (*((*it_h).second))["pt_l2"]     = new TH1D(TString("pt_l2") +"_"+(*it_h).first+str,";#bf{l2 p_{T} [GeV]};",     37,0,65);  (*((*it_h).second))["pt_l2"]->Sumw2();
->
+  //   (*((*it_h).second))["pt_l3"]     = new TH1D(TString("pt_l3") +"_"+(*it_h).first+str,";#bf{l3 p_{T} [GeV]};",     37,0,75);  (*((*it_h).second))["pt_l3"]->Sumw2();
->
+  //   (*((*it_h).second))["pt_l4"]     = new TH1D(TString("pt_l4") +"_"+(*it_h).first+str,";#bf{l4 p_{T} [GeV]};",     37,0,50);  (*((*it_h).second))["pt_l4"]->Sumw2();
->
->
+  //   (*((*it_h).second))["eta_l1"]    = new TH1D(TString("eta_l1")+"_"+(*it_h).first+str,";#bf{l1 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l1"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_l2"]    = new TH1D(TString("eta_l2")+"_"+(*it_h).first+str,";#bf{l2 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l2"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_l3"]    = new TH1D(TString("eta_l3")+"_"+(*it_h).first+str,";#bf{l3 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l3"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_l4"]    = new TH1D(TString("eta_l4")+"_"+(*it_h).first+str,";#bf{l4 #eta};",            25,-2.5,2.5);     (*((*it_h).second))["eta_l4"]->Sumw2();
->
->
+  //   (*((*it_h).second))["ip3ds_l1"]  = new TH1D(TString("ip3ds_l1") +"_"+(*it_h).first+str,";#bf{ip3ds l1};",37,-3,3);       (*((*it_h).second))["ip3ds_l1"]->Sumw2();
->
+  //   (*((*it_h).second))["ip3ds_l2"]  = new TH1D(TString("ip3ds_l2") +"_"+(*it_h).first+str,";#bf{ip3ds l2};",37,-3,3);       (*((*it_h).second))["ip3ds_l2"]->Sumw2();
->
+  //   (*((*it_h).second))["ip3ds_l3_lo"]       = new TH1D(TString("ip3ds_l3_lo") +"_"+(*it_h).first+str,";#bf{ip3ds l3_lo};",25,-6,6); (*((*it_h).second))["ip3ds_l3_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["ip3ds_l4_lo"]       = new TH1D(TString("ip3ds_l4_lo") +"_"+(*it_h).first+str,";#bf{ip3ds l4_lo};",25,-6,6); (*((*it_h).second))["ip3ds_l4_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["ip3ds_l3"]  = new TH1D(TString("ip3ds_l3") +"_"+(*it_h).first+str,";#bf{ip3ds l3};",55,-40,40);     (*((*it_h).second))["ip3ds_l3"]->Sumw2();
->
+  //   (*((*it_h).second))["ip3ds_l4"]  = new TH1D(TString("ip3ds_l4") +"_"+(*it_h).first+str,";#bf{ip3ds l4};",55,-40,40);     (*((*it_h).second))["ip3ds_l4"]->Sumw2();
->
->
+  //   (*((*it_h).second))["d0_l1"]     = new TH1D(TString("d0_l1") +"_"+(*it_h).first+str,";#bf{d0 l1};",50,-.008,.008);      (*((*it_h).second))["d0_l1"]->Sumw2();
->
+  //   (*((*it_h).second))["d0_l2"]     = new TH1D(TString("d0_l2") +"_"+(*it_h).first+str,";#bf{d0 l2};",50,-.008,.008);      (*((*it_h).second))["d0_l2"]->Sumw2();
->
+  //   (*((*it_h).second))["d0_l3_lo"]  = new TH1D(TString("d0_l3_lo") +"_"+(*it_h).first+str,";#bf{d0 l3_lo};",25,-.03,.03);  (*((*it_h).second))["d0_l3_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["d0_l4_lo"]  = new TH1D(TString("d0_l4_lo") +"_"+(*it_h).first+str,";#bf{d0 l4_lo};",25,-.03,.03);  (*((*it_h).second))["d0_l4_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["d0_l3"]     = new TH1D(TString("d0_l3") +"_"+(*it_h).first+str,";#bf{d0 l3};",50,-.1,.1);          (*((*it_h).second))["d0_l3"]->Sumw2();
->
+  //   (*((*it_h).second))["d0_l4"]     = new TH1D(TString("d0_l4") +"_"+(*it_h).first+str,";#bf{d0 l4};",50,-.1,.1);          (*((*it_h).second))["d0_l4"]->Sumw2();
->
->
+  //   (*((*it_h).second))["dz_l1"]     = new TH1D(TString("dz_l1") +"_"+(*it_h).first+str,";#bf{dz l1};",50,-.01,.01);   (*((*it_h).second))["dz_l1"]->Sumw2();
->
+  //   (*((*it_h).second))["dz_l2"]     = new TH1D(TString("dz_l2") +"_"+(*it_h).first+str,";#bf{dz l2};",50,-.01,.01);   (*((*it_h).second))["dz_l2"]->Sumw2();
->
+  //   (*((*it_h).second))["dz_l3"]     = new TH1D(TString("dz_l3") +"_"+(*it_h).first+str,";#bf{dz l3};",50,-.2,.2);   (*((*it_h).second))["dz_l3"]->Sumw2();
->
+  //   (*((*it_h).second))["dz_l4"]     = new TH1D(TString("dz_l4") +"_"+(*it_h).first+str,";#bf{dz l4};",50,-.2,.2);   (*((*it_h).second))["dz_l4"]->Sumw2();
->
->
+  //   (*((*it_h).second))["dR_l3l4_lo"]        = new TH1D(TString("dR_l3l4_lo") +"_"+(*it_h).first+str,";#bf{#Delta R l3,l4};",75,0,.2);    (*((*it_h).second))["dR_l3l4_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["dR_l3l4"]   = new TH1D(TString("dR_l3l4") +"_"+(*it_h).first+str,";#bf{#Delta R l3,l4};",75,0,6);        (*((*it_h).second))["dR_l3l4"]->Sumw2();
->
->
+  //   // jets
->
+  //   (*((*it_h).second))["pt_allJet"] = new TH1D(TString("pt_allJet") +"_"+(*it_h).first+str,";#bf{jet 1 p_{T} [GeV]};",37,0,400);  (*((*it_h).second))["pt_allJet"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_allJet"]        = new TH1D(TString("eta_allJet") +"_"+(*it_h).first+str,";#bf{jet 1 p_{T} [GeV]};",37,-4.8,4.8);  (*((*it_h).second))["eta_allJet"]->Sumw2();
->
+  //   (*((*it_h).second))["pt_j1"]     = new TH1D(TString("pt_j1") +"_"+(*it_h).first+str,";#bf{jet 1 p_{T} [GeV]};",  37,0,400);  (*((*it_h).second))["pt_j1"]->Sumw2();
->
+  //   (*((*it_h).second))["pt_j2"]     = new TH1D(TString("pt_j2") +"_"+(*it_h).first+str,";#bf{jet 2 p_{T} [GeV]};",  37,0,125);  (*((*it_h).second))["pt_j2"]->Sumw2();
->
+  //   (*((*it_h).second))["pt_j3"]     = new TH1D(TString("pt_j3") +"_"+(*it_h).first+str,";#bf{jet 3 p_{T} [GeV]};",  37,0,125);  (*((*it_h).second))["pt_j3"]->Sumw2();
->
+  //   (*((*it_h).second))["pt_j4"]     = new TH1D(TString("pt_j4") +"_"+(*it_h).first+str,";#bf{jet 4 p_{T} [GeV]};",  37,0,125);  (*((*it_h).second))["pt_j4"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_j1"]    = new TH1D(TString("eta_j1")+"_"+(*it_h).first+str,";#bf{jet 1 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j1"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_j2"]    = new TH1D(TString("eta_j2")+"_"+(*it_h).first+str,";#bf{jet 2 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j2"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_j3"]    = new TH1D(TString("eta_j3")+"_"+(*it_h).first+str,";#bf{jet 3 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j3"]->Sumw2();
->
+  //   (*((*it_h).second))["eta_j4"]    = new TH1D(TString("eta_j4")+"_"+(*it_h).first+str,";#bf{jet 4 #eta};",         55,-4.8,4.8);     (*((*it_h).second))["eta_j4"]->Sumw2();
->
+  //   (*((*it_h).second))["phi_j1"]    = new TH1D(TString("phi_j1")+"_"+(*it_h).first+str,";#bf{jet 1 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j1"]->Sumw2();
->
+  //   (*((*it_h).second))["phi_j2"]    = new TH1D(TString("phi_j2")+"_"+(*it_h).first+str,";#bf{jet 2 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j2"]->Sumw2();
->
+  //   (*((*it_h).second))["phi_j3"]    = new TH1D(TString("phi_j3")+"_"+(*it_h).first+str,";#bf{jet 3 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j3"]->Sumw2();
->
+  //   (*((*it_h).second))["phi_j4"]    = new TH1D(TString("phi_j4")+"_"+(*it_h).first+str,";#bf{jet 4 #phi};",         25,-3.2,3.2);     (*((*it_h).second))["phi_j4"]->Sumw2();
->
+  //   (*((*it_h).second))["idmva_j1"]  = new TH1D(TString("idmva_j1")+"_"+(*it_h).first+str,";#bf{jet 1 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j1"]->Sumw2();
->
+  //   (*((*it_h).second))["idmva_j2"]  = new TH1D(TString("idmva_j2")+"_"+(*it_h).first+str,";#bf{jet 2 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j2"]->Sumw2();
->
+  //   (*((*it_h).second))["idmva_j3"]  = new TH1D(TString("idmva_j3")+"_"+(*it_h).first+str,";#bf{jet 3 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j3"]->Sumw2();
->
+  //   (*((*it_h).second))["idmva_j4"]  = new TH1D(TString("idmva_j4")+"_"+(*it_h).first+str,";#bf{jet 4 MVA};",        25,-1.1,1.1);     (*((*it_h).second))["idmva_j4"]->Sumw2();
->
+  //   (*((*it_h).second))["mjj"]               = new TH1D(TString("mjj")+"_"+(*it_h).first+str,";#bf{m_{jj}};",                25,10,1550);     (*((*it_h).second))["mjj"]->Sumw2();
->
+  //   (*((*it_h).second))["dEta"]              = new TH1D(TString("dEta")+"_"+(*it_h).first+str,";#bf{#Delta #eta};",          25,-6,6);       (*((*it_h).second))["dEta"]->Sumw2();
->
+  //   (*((*it_h).second))["etaProd"]   = new TH1D(TString("etaProd")+"_"+(*it_h).first+str,";#bf{#eta_{1} * #eta_{2}};",25,-12,12);      (*((*it_h).second))["etaProd"]->Sumw2();
->
+  //   (*((*it_h).second))["njets"]     = new TH1D(TString("njets")+"_"+(*it_h).first+str,";#bf{N Jets};",      7,-0.5,6.5);      (*((*it_h).second))["njets"]->Sumw2();
->
->
+  //   // (*((*it_h).second))["dR_l1j1"]        = new TH1D(TString("dR_l1j1") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 1};",75,0,6);        (*((*it_h).second))["dR_l1j1"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l2j1"]        = new TH1D(TString("dR_l2j1") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 1};",75,0,6);        (*((*it_h).second))["dR_l2j1"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l3j1"]        = new TH1D(TString("dR_l3j1") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 1};",75,0,6);        (*((*it_h).second))["dR_l3j1"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l4j1"]        = new TH1D(TString("dR_l4j1") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 1};",75,0,6);        (*((*it_h).second))["dR_l4j1"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l1j2"]        = new TH1D(TString("dR_l1j2") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 2};",75,0,6);        (*((*it_h).second))["dR_l1j2"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l2j2"]        = new TH1D(TString("dR_l2j2") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 2};",75,0,6);        (*((*it_h).second))["dR_l2j2"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l3j2"]        = new TH1D(TString("dR_l3j2") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 2};",75,0,6);        (*((*it_h).second))["dR_l3j2"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l4j2"]        = new TH1D(TString("dR_l4j2") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 2};",75,0,6);        (*((*it_h).second))["dR_l4j2"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l1j3"]        = new TH1D(TString("dR_l1j3") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 3};",75,0,6);        (*((*it_h).second))["dR_l1j3"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l2j3"]        = new TH1D(TString("dR_l2j3") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 3};",75,0,6);        (*((*it_h).second))["dR_l2j3"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l3j3"]        = new TH1D(TString("dR_l3j3") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 3};",75,0,6);        (*((*it_h).second))["dR_l3j3"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l4j3"]        = new TH1D(TString("dR_l4j3") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 3};",75,0,6);        (*((*it_h).second))["dR_l4j3"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l1j4"]        = new TH1D(TString("dR_l1j4") +"_"+(*it_h).first+str,";#bf{#Delta R l1,jet 4};",75,0,6);        (*((*it_h).second))["dR_l1j4"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l2j4"]        = new TH1D(TString("dR_l2j4") +"_"+(*it_h).first+str,";#bf{#Delta R l2,jet 4};",75,0,6);        (*((*it_h).second))["dR_l2j4"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l3j4"]        = new TH1D(TString("dR_l3j4") +"_"+(*it_h).first+str,";#bf{#Delta R l3,jet 4};",75,0,6);        (*((*it_h).second))["dR_l3j4"]->Sumw2();
->
+  //   // (*((*it_h).second))["dR_l4j4"]        = new TH1D(TString("dR_l4j4") +"_"+(*it_h).first+str,";#bf{#Delta R l4,jet 4};",75,0,6);        (*((*it_h).second))["dR_l4j4"]->Sumw2();
->
->
+  //   int nbins=35;
->
+  //   (*((*it_h).second))["costheta1"] = new TH1D(TString("costheta1") +"_"+(*it_h).first+str,";#bf{costheta1};",                nbins,-1,1);          (*((*it_h).second))["costheta1"]->Sumw2();
->
+  //   (*((*it_h).second))["costheta2"] = new TH1D(TString("costheta2") +"_"+(*it_h).first+str,";#bf{costheta2};",                nbins,-1,1);          (*((*it_h).second))["costheta2"]->Sumw2();
->
+  //   (*((*it_h).second))["costhetastar"]      = new TH1D(TString("costhetastar") +"_"+(*it_h).first+str,";#bf{costhetastar};",          nbins,-1,1);          (*((*it_h).second))["costhetastar"]->Sumw2();
->
+  //   (*((*it_h).second))["Phi"]               = new TH1D(TString("Phi") +"_"+(*it_h).first+str,";#bf{Phi};",                            nbins,-3.15,3.15);          (*((*it_h).second))["Phi"]->Sumw2();
->
+  //   (*((*it_h).second))["Phi1"]              = new TH1D(TString("Phi1") +"_"+(*it_h).first+str,";#bf{Phi1};",                          nbins,-3.15,3.15);          (*((*it_h).second))["Phi1"]->Sumw2();
->
->
+  //   (*((*it_h).second))["pt4l"]              = new TH1D(TString("pt4l") +"_"+(*it_h).first+str,";#bf{pt4l/m4l};",                      nbins-7,0,0.5);          (*((*it_h).second))["pt4l"]->Sumw2();
->
+  //   (*((*it_h).second))["pt4l_log"]  = new TH1D(TString("pt4l_log") +"_"+(*it_h).first+str,";#bf{pt4l/m4l};",                  nbins,0,1.5);          (*((*it_h).second))["pt4l_log"]->Sumw2();
->
+  //   (*((*it_h).second))["y4l"]               = new TH1D(TString("y4l") +"_"+(*it_h).first+str,";#bf{y4l};",                            nbins,-2.4,2.4);      (*((*it_h).second))["y4l"]->Sumw2();
->
+  //   (*((*it_h).second))["Z1pt"]              = new TH1D(TString("Z1pt") +"_"+(*it_h).first+str,";#bf{Z1pt/m4l};",                      nbins,0,0.5);          (*((*it_h).second))["Z1pt"]->Sumw2();
->
+  //   (*((*it_h).second))["Z1pt_log"]  = new TH1D(TString("Z1pt_log") +"_"+(*it_h).first+str,";#bf{Z1pt/m4l};",                  nbins,0,1.5);          (*((*it_h).second))["Z1pt_log"]->Sumw2();
->
+  //   (*((*it_h).second))["Z2pt"]              = new TH1D(TString("Z2pt") +"_"+(*it_h).first+str,";#bf{Z2pt/m4l};",                      nbins,0,0.5);          (*((*it_h).second))["Z2pt"]->Sumw2();
->
+  //   (*((*it_h).second))["Z2pt_log"]  = new TH1D(TString("Z2pt_log") +"_"+(*it_h).first+str,";#bf{Z2pt/m4l};",                  nbins,0,1.5);          (*((*it_h).second))["Z2pt_log"]->Sumw2();
->
->
+  //   (*((*it_h).second))["ZZdotZ1"]   = new TH1D(TString("ZZdotZ1") +"_"+(*it_h).first+str,";#bf{ZZ #bullet Z1/(m4l*mZ1)};",    nbins,-.2,2.5);          (*((*it_h).second))["ZZdotZ1"]->Sumw2();
->
+  //   (*((*it_h).second))["ZZdotZ2"]   = new TH1D(TString("ZZdotZ2") +"_"+(*it_h).first+str,";#bf{ZZ #bullet Z2/(m4l*mZ2)};",    nbins,-1,3);          (*((*it_h).second))["ZZdotZ2"]->Sumw2();
->
+  //   (*((*it_h).second))["ZZdotZ1_log"]       = new TH1D(TString("ZZdotZ1_log") +"_"+(*it_h).first+str,";#bf{ZZ #bullet Z1/(m4l*mZ1)};",nbins,-1.5,15);          (*((*it_h).second))["ZZdotZ1_log"]->Sumw2();
->
+  //   (*((*it_h).second))["ZZdotZ2_log"]       = new TH1D(TString("ZZdotZ2_log") +"_"+(*it_h).first+str,";#bf{ZZ #bullet Z2/(m4l*mZ2)};",nbins,-1.5,15);          (*((*it_h).second))["ZZdotZ2_log"]->Sumw2();
->
+  //   (*((*it_h).second))["dphi1"]     = new TH1D(TString("dphi1") +"_"+(*it_h).first+str,";#bf{dphi1};",                        nbins,-1,1);          (*((*it_h).second))["dphi1"]->Sumw2();
->
+  //   (*((*it_h).second))["dphi2"]     = new TH1D(TString("dphi2") +"_"+(*it_h).first+str,";#bf{dphi2};",                        nbins,-1,1);          (*((*it_h).second))["dphi2"]->Sumw2();
->
+  //   (*((*it_h).second))["dphi1_log"] = new TH1D(TString("dphi1_log") +"_"+(*it_h).first+str,";#bf{dphi1};",                    nbins,-1,1);          (*((*it_h).second))["dphi1_log"]->Sumw2();
->
+  //   (*((*it_h).second))["dphi2_log"] = new TH1D(TString("dphi2_log") +"_"+(*it_h).first+str,";#bf{dphi2};",                    nbins,-1,1);          (*((*it_h).second))["dphi2_log"]->Sumw2();
->
->
+  //   (*((*it_h).second))["fusionMVA"]         = new TH1D(TString("fusionMVA")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,-1,1);      (*((*it_h).second))["fusionMVA"]->Sumw2();
->
+  //   (*((*it_h).second))["fusionMVA_lo"]              = new TH1D(TString("fusionMVA_lo")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,-1,-.83);      (*((*it_h).second))["fusionMVA_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["fusionMVA_med"]             = new TH1D(TString("fusionMVA_med")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,-.5,1);      (*((*it_h).second))["fusionMVA_med"]->Sumw2();
->
+  //   (*((*it_h).second))["fusionMVA_hi"]              = new TH1D(TString("fusionMVA_hi")+"_"+(*it_h).first+str,";#bf{MVA output};",        100,.75,1);      (*((*it_h).second))["fusionMVA_hi"]->Sumw2();
->
+  //   (*((*it_h).second))["m4l_jet_lo"]                = new TH1D(TString("m4l_jet_lo") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",     45,100,180);          (*((*it_h).second))["m4l_jet_lo"]->Sumw2();
->
+  //   (*((*it_h).second))["m4l_jet"]           = new TH1D(TString("m4l_jet") +"_"+(*it_h).first+str,";#bf{m4l [GeV]};",        45,100,600);          (*((*it_h).second))["m4l_jet"]->Sumw2();
->
+  // }
+  return hists;
+//--------------------------------------------------------------------------------------------------
-<
+void makeHTML(FOFlags &ctrl, TString type, TString plotLabel)
-<
+{
-<
+  TString tmpString(plotLabel);
-<
+  tmpString.ReplaceAll("/"," ");
-<
+  TString title(ctrl.faketype+", "+tmpString+", "+type);
-<
+  TString htmlfname(ctrl.outdir+"/"+plotLabel+"/"+type+"/plots.html");
-<
+  ofstream htmlfile(htmlfname);
-<
-<
+  htmlfile << "<!DOCTYPE html" << endl;
-<
+  htmlfile << "    PUBLIC \"-//W3C//DTD HTML 3.2//EN\">" << endl;
-<
+  htmlfile << "<html>" << endl;
-<
-<
+  htmlfile << "<head><title>"+title+"</title></head>" << endl;
-<
+  htmlfile << "<body bgcolor=\"000000\">" << endl;
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">"+title+"</h3>" << endl;
-<
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">boson kinematics</h3>" << endl;
-<
+  htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_lo.png\"><img src=\"plots/m4l_lo.png\" alt=\"plots/m4l_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_med.png\"><img src=\"plots/m4l_med.png\" alt=\"plots/m4l_med.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l.png\"><img src=\"plots/m4l.png\" alt=\"plots/m4l.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mZ1.png\"><img src=\"plots/mZ1.png\" alt=\"plots/mZ1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mZ2_lo.png\"><img src=\"plots/mZ2_lo.png\" alt=\"plots/mZ2_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mZ2.png\"><img src=\"plots/mZ2.png\" alt=\"plots/mZ2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Z1pt.png\"><img src=\"plots/Z1pt.png\" alt=\"plots/Z1pt.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ZZpt.png\"><img src=\"plots/ZZpt.png\" alt=\"plots/ZZpt.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/met.png\"><img src=\"plots/met.png\" alt=\"plots/met.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "</table>" << endl;
-<
+  htmlfile << "<hr />" << endl;
-<
+  htmlfile << "control: " << endl;
-<
+  htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/run.png\"><img src=\"plots/run.png\" alt=\"plots/run.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/npv.png\"><img src=\"plots/npv.png\" alt=\"plots/npv.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "</table>" << endl;
-<
-<
+  htmlfile << "</body>" << endl;
-<
+  htmlfile << "</html>" << endl;
-<
+  htmlfile.close();
-<
-<
+  TString leptonfname(htmlfname);
-<
+  leptonfname.ReplaceAll("plots.html","leptonplots.html");
-<
+  htmlfile.open(leptonfname);
-<
+  cout << "opening: " << leptonfname << endl;
-<
-<
+  htmlfile << "<!DOCTYPE html" << endl;
-<
+  htmlfile << "    PUBLIC \"-//W3C//DTD HTML 3.2//EN\">" << endl;
-<
+  htmlfile << "<html>" << endl;
-<
-<
+  htmlfile << "<head><title>"+title+"</title></head>" << endl;
-<
+  htmlfile << "<body bgcolor=\"000000\">" << endl;
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">"+title+"</h3>" << endl;
-<
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">lepton plots</h3>" << endl;
-<
+  htmlfile << "<hr />" << endl;
-<
+  htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
-<
+  htmlfile << "</table>" << endl;
-<
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">Z1 leptons:</h3>" << endl;
-<
+  htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l1.png\"><img src=\"plots/pt_l1.png\" alt=\"plots/pt_l1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l2.png\"><img src=\"plots/pt_l2.png\" alt=\"plots/pt_l2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l1.png\"><img src=\"plots/eta_l1.png\" alt=\"plots/eta_l1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l2.png\"><img src=\"plots/eta_l2.png\" alt=\"plots/eta_l2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l1.png\"><img src=\"plots/dz_l1.png\" alt=\"plots/dz_l1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l2.png\"><img src=\"plots/dz_l2.png\" alt=\"plots/dz_l2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l1.png\"><img src=\"plots/ip3ds_l1.png\" alt=\"plots/ip3ds_l1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l2.png\"><img src=\"plots/ip3ds_l2.png\" alt=\"plots/ip3ds_l2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l1.png\"><img src=\"plots/d0_l1.png\" alt=\"plots/d0_l1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l2.png\"><img src=\"plots/d0_l2.png\" alt=\"plots/d0_l2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  // htmlfile << "<tr>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3l4_lo.png\"><img src=\"plots/dR_l3l4_lo.png\" alt=\"plots/dR_l3l4_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3l4.png\"><img src=\"plots/dR_l3l4.png\" alt=\"plots/dR_l3l4.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a><</a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a><</a></td>" << endl;
-<
+  // htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "</table>" << endl;
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">extra leptons (l3 and l4): </h3>" << endl;
-<
+  htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l3.png\"><img src=\"plots/pt_l3.png\" alt=\"plots/pt_l3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l4.png\"><img src=\"plots/pt_l4.png\" alt=\"plots/pt_l4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l3.png\"><img src=\"plots/eta_l3.png\" alt=\"plots/eta_l3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l4.png\"><img src=\"plots/eta_l4.png\" alt=\"plots/eta_l4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l3_lo.png\"><img src=\"plots/ip3ds_l3_lo.png\" alt=\"plots/ip3ds_l3_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l4_lo.png\"><img src=\"plots/ip3ds_l4_lo.png\" alt=\"plots/ip3ds_l4_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l3.png\"><img src=\"plots/ip3ds_l3.png\" alt=\"plots/ip3ds_l3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l4.png\"><img src=\"plots/ip3ds_l4.png\" alt=\"plots/ip3ds_l4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
->
+// void makeHTML(FOFlags &ctrl, TString type, TString plotLabel, TString fullOutDir)
->
+// {
->
+//   TString tmpString(plotLabel);
->
+//   tmpString.ReplaceAll("/"," ");
->
+//   TString title(ctrl.faketype+", "+tmpString+", "+type);
->
+//   TString htmlfname(fullOutDir+"/plots.html");
->
+//   ofstream htmlfile(htmlfname);
->
->
+//   htmlfile << "<!DOCTYPE html" << endl;
->
+//   htmlfile << "    PUBLIC \"-//W3C//DTD HTML 3.2//EN\">" << endl;
->
+//   htmlfile << "<html>" << endl;
->
->
+//   htmlfile << "<head><title>"+title+"</title></head>" << endl;
->
+//   htmlfile << "<body bgcolor=\"000000\">" << endl;
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">"+title+"</h3>" << endl;
->
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">boson kinematics</h3>" << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_lo.png\"><img src=\"plots/m4l_lo.png\" alt=\"plots/m4l_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_med.png\"><img src=\"plots/m4l_med.png\" alt=\"plots/m4l_med.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l.png\"><img src=\"plots/m4l.png\" alt=\"plots/m4l.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mZ1.png\"><img src=\"plots/mZ1.png\" alt=\"plots/mZ1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mZ2_lo.png\"><img src=\"plots/mZ2_lo.png\" alt=\"plots/mZ2_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mZ2.png\"><img src=\"plots/mZ2.png\" alt=\"plots/mZ2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Z1pt.png\"><img src=\"plots/Z1pt.png\" alt=\"plots/Z1pt.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ZZpt.png\"><img src=\"plots/ZZpt.png\" alt=\"plots/ZZpt.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/met.png\"><img src=\"plots/met.png\" alt=\"plots/met.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "</table>" << endl;
->
+//   htmlfile << "<hr />" << endl;
->
+//   htmlfile << "control: " << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/run.png\"><img src=\"plots/run.png\" alt=\"plots/run.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/npv.png\"><img src=\"plots/npv.png\" alt=\"plots/npv.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "</table>" << endl;
->
->
+//   htmlfile << "</body>" << endl;
->
+//   htmlfile << "</html>" << endl;
->
+//   htmlfile.close();
->
->
+//   TString anglefname(htmlfname);
->
+//   anglefname.ReplaceAll("plots.html","angleplots.html");
->
+//   htmlfile.open(anglefname);
->
->
+//   htmlfile << "<!DOCTYPE html" << endl;
->
+//   htmlfile << "    PUBLIC \"-//W3C//DTD HTML 3.2//EN\">" << endl;
->
+//   htmlfile << "<html>" << endl;
->
->
+//   htmlfile << "<head><title>"+title+"</title></head>" << endl;
->
+//   htmlfile << "<body bgcolor=\"000000\">" << endl;
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">"+title+"</h3>" << endl;
->
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">Angles & Co.</h3>" << endl;
->
+//   htmlfile << "<hr />" << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
+//   htmlfile << "</table>" << endl;
->
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">Angles:</h3>" << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/costheta1.png\"><img src=\"plots/costheta1.png\" alt=\"plots/costheta1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/costheta2.png\"><img src=\"plots/costheta2.png\" alt=\"plots/costheta2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/costhetastar.png\"><img src=\"plots/costhetastar.png\" alt=\"plots/costhetastar.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Phi.png\"><img src=\"plots/Phi.png\" alt=\"plots/Phi.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Phi1.png\"><img src=\"plots/Phi1.png\" alt=\"plots/Phi1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "</table>" << endl;
->
+//   htmlfile << "<hr />" << endl;
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">Boson pt variables: " << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/y4l.png\"><img src=\"plots/y4l.png\" alt=\"plots/y4l.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt4l.png\"><img src=\"plots/pt4l.png\" alt=\"plots/pt4l.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt4l_log.png\"><img src=\"plots/pt4l_log.png\" alt=\"plots/pt4l_log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<tr>" << endl;
->
->
+//   htmlfile << "</tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Z1pt.png\"><img src=\"plots/Z1pt.png\" alt=\"plots/Z1pt.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Z2pt.png\"><img src=\"plots/Z2pt.png\" alt=\"plots/Z2pt.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Z1pt_log.png\"><img src=\"plots/Z1pt_log.png\" alt=\"plots/Z1pt_log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/Z2pt_log.png\"><img src=\"plots/Z2pt_log.png\" alt=\"plots/Z2pt_log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ZZdotZ1.png\"><img src=\"plots/ZZdotZ1.png\" alt=\"plots/ZZdotZ1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ZZdotZ2.png\"><img src=\"plots/ZZdotZ2.png\" alt=\"plots/ZZdotZ2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ZZdotZ1_log.png\"><img src=\"plots/ZZdotZ1_log.png\" alt=\"plots/ZZdotZ1_log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ZZdotZ2_log.png\"><img src=\"plots/ZZdotZ2_log.png\" alt=\"plots/ZZdotZ2_log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dphi1.png\"><img src=\"plots/dphi1.png\" alt=\"plots/dphi1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dphi2.png\"><img src=\"plots/dphi2.png\" alt=\"plots/dphi2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dphi1_log.png\"><img src=\"plots/dphi1_log.png\" alt=\"plots/dphi1_log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dphi2_log.png\"><img src=\"plots/dphi2_log.png\" alt=\"plots/dphi2_log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
->
+//   htmlfile << "</table>" << endl;
->
->
+//   htmlfile << "</body>" << endl;
->
+//   htmlfile << "</html>" << endl;
->
+//   htmlfile.close();
->
->
+//   TString leptonfname(htmlfname);
->
+//   leptonfname.ReplaceAll("plots.html","leptonplots.html");
->
+//   htmlfile.open(leptonfname);
->
->
+//   htmlfile << "<!DOCTYPE html" << endl;
->
+//   htmlfile << "    PUBLIC \"-//W3C//DTD HTML 3.2//EN\">" << endl;
->
+//   htmlfile << "<html>" << endl;
->
->
+//   htmlfile << "<head><title>"+title+"</title></head>" << endl;
->
+//   htmlfile << "<body bgcolor=\"000000\">" << endl;
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">"+title+"</h3>" << endl;
->
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">lepton plots</h3>" << endl;
->
+//   htmlfile << "<hr />" << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
+//   htmlfile << "</table>" << endl;
->
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">Z1 leptons:</h3>" << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l1.png\"><img src=\"plots/pt_l1.png\" alt=\"plots/pt_l1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l2.png\"><img src=\"plots/pt_l2.png\" alt=\"plots/pt_l2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l1.png\"><img src=\"plots/eta_l1.png\" alt=\"plots/eta_l1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l2.png\"><img src=\"plots/eta_l2.png\" alt=\"plots/eta_l2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l1.png\"><img src=\"plots/ip3ds_l1.png\" alt=\"plots/ip3ds_l1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l2.png\"><img src=\"plots/ip3ds_l2.png\" alt=\"plots/ip3ds_l2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l1.png\"><img src=\"plots/d0_l1.png\" alt=\"plots/d0_l1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l2.png\"><img src=\"plots/d0_l2.png\" alt=\"plots/d0_l2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l1.png\"><img src=\"plots/dz_l1.png\" alt=\"plots/dz_l1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l2.png\"><img src=\"plots/dz_l2.png\" alt=\"plots/dz_l2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   // htmlfile << "<tr>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3l4_lo.png\"><img src=\"plots/dR_l3l4_lo.png\" alt=\"plots/dR_l3l4_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3l4.png\"><img src=\"plots/dR_l3l4.png\" alt=\"plots/dR_l3l4.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a><</a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a><</a></td>" << endl;
->
+//   // htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "</table>" << endl;
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">extra leptons (l3 and l4): </h3>" << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l3.png\"><img src=\"plots/pt_l3.png\" alt=\"plots/pt_l3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_l4.png\"><img src=\"plots/pt_l4.png\" alt=\"plots/pt_l4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l3.png\"><img src=\"plots/eta_l3.png\" alt=\"plots/eta_l3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_l4.png\"><img src=\"plots/eta_l4.png\" alt=\"plots/eta_l4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l3_lo.png\"><img src=\"plots/ip3ds_l3_lo.png\" alt=\"plots/ip3ds_l3_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l4_lo.png\"><img src=\"plots/ip3ds_l4_lo.png\" alt=\"plots/ip3ds_l4_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l3.png\"><img src=\"plots/ip3ds_l3.png\" alt=\"plots/ip3ds_l3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/ip3ds_l4.png\"><img src=\"plots/ip3ds_l4.png\" alt=\"plots/ip3ds_l4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l3_lo.png\"><img src=\"plots/d0_l3_lo.png\" alt=\"plots/d0_l3_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l4_lo.png\"><img src=\"plots/d0_l4_lo.png\" alt=\"plots/d0_l4_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l3.png\"><img src=\"plots/d0_l3.png\" alt=\"plots/d0_l3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l4.png\"><img src=\"plots/d0_l4.png\" alt=\"plots/d0_l4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l3.png\"><img src=\"plots/dz_l3.png\" alt=\"plots/dz_l3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l4.png\"><img src=\"plots/dz_l4.png\" alt=\"plots/dz_l4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l3_lo.png\"><img src=\"plots/d0_l3_lo.png\" alt=\"plots/d0_l3_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l4_lo.png\"><img src=\"plots/d0_l4_lo.png\" alt=\"plots/d0_l4_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l3.png\"><img src=\"plots/d0_l3.png\" alt=\"plots/d0_l3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/d0_l4.png\"><img src=\"plots/d0_l4.png\" alt=\"plots/d0_l4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l3.png\"><img src=\"plots/dz_l3.png\" alt=\"plots/dz_l3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dz_l4.png\"><img src=\"plots/dz_l4.png\" alt=\"plots/dz_l4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
-<
+  htmlfile << "</table>" << endl;
->
+//   htmlfile << "</table>" << endl;
-<
+  htmlfile << "</body>" << endl;
-<
+  htmlfile << "</html>" << endl;
-<
+  htmlfile.close();
-<
-<
+  TString jetfname(htmlfname);
-<
+  jetfname.ReplaceAll("plots.html","jetplots.html");
-<
+  htmlfile.open(jetfname);
-<
-<
+  htmlfile << "<!DOCTYPE html" << endl;
-<
+  htmlfile << "    PUBLIC \"-//W3C//DTD HTML 3.2//EN\">" << endl;
-<
+  htmlfile << "<html>" << endl;
-<
-<
+  htmlfile << "<head><title>"+title+"</title></head>" << endl;
-<
+  htmlfile << "<body bgcolor=\"000000\">" << endl;
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">"+title+"</h3>" << endl;
-<
-<
+  htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">jet plots</h3>" << endl;
-<
+  htmlfile << "<hr />" << endl;
-<
+  htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j1.png\"><img src=\"plots/pt_j1.png\" alt=\"plots/pt_j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j2.png\"><img src=\"plots/pt_j2.png\" alt=\"plots/pt_j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j2.png\"><img src=\"plots/pt_j2.png\" alt=\"plots/pt_j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j4.png\"><img src=\"plots/pt_j4.png\" alt=\"plots/pt_j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j1.png\"><img src=\"plots/eta_j1.png\" alt=\"plots/eta_j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j2.png\"><img src=\"plots/eta_j2.png\" alt=\"plots/eta_j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j3.png\"><img src=\"plots/eta_j3.png\" alt=\"plots/eta_j3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j4.png\"><img src=\"plots/eta_j4.png\" alt=\"plots/eta_j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j1.png\"><img src=\"plots/idmva_j1.png\" alt=\"plots/idmva_j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j2.png\"><img src=\"plots/idmva_j2.png\" alt=\"plots/idmva_j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j3.png\"><img src=\"plots/idmva_j3.png\" alt=\"plots/idmva_j3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j4.png\"><img src=\"plots/idmva_j4.png\" alt=\"plots/idmva_j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j1.png\"><img src=\"plots/phi_j1.png\" alt=\"plots/phi_j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j2.png\"><img src=\"plots/phi_j2.png\" alt=\"plots/phi_j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j3.png\"><img src=\"plots/phi_j3.png\" alt=\"plots/phi_j3.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j4.png\"><img src=\"plots/phi_j4.png\" alt=\"plots/phi_j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mjj.png\"><img src=\"plots/mjj.png\" alt=\"plots/mjj.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dEta.png\"><img src=\"plots/dEta.png\" alt=\"plots/dEta.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/etaProd.png\"><img src=\"plots/etaProd.png\" alt=\"plots/etaProd.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/njets.png\"><img src=\"plots/njets.png\" alt=\"plots/njets.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  // htmlfile << "<tr>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j1.png\"><img src=\"plots/dR_l1j1.png\" alt=\"plots/dR_l1j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j1.png\"><img src=\"plots/dR_l2j1.png\" alt=\"plots/dR_l2j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j1.png\"><img src=\"plots/dR_l3j1.png\" alt=\"plots/dR_l3j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j1.png\"><img src=\"plots/dR_l4j1.png\" alt=\"plots/dR_l4j1.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "</tr>" << endl;
-<
-<
+  // htmlfile << "<tr>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j2.png\"><img src=\"plots/dR_l1j2.png\" alt=\"plots/dR_l1j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j2.png\"><img src=\"plots/dR_l2j2.png\" alt=\"plots/dR_l2j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j2.png\"><img src=\"plots/dR_l3j2.png\" alt=\"plots/dR_l3j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j2.png\"><img src=\"plots/dR_l4j2.png\" alt=\"plots/dR_l4j2.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "</tr>" << endl;
-<
-<
+  // htmlfile << "<tr>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j3.png\"><img src=\"plots/dR_l1j3.png\" alt=\"plots/dR_l1j3.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j3.png\"><img src=\"plots/dR_l2j3.png\" alt=\"plots/dR_l2j3.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j3.png\"><img src=\"plots/dR_l3j3.png\" alt=\"plots/dR_l3j3.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j3.png\"><img src=\"plots/dR_l4j3.png\" alt=\"plots/dR_l4j3.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "</tr>" << endl;
-<
-<
+  // htmlfile << "<tr>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j4.png\"><img src=\"plots/dR_l1j4.png\" alt=\"plots/dR_l1j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j4.png\"><img src=\"plots/dR_l2j4.png\" alt=\"plots/dR_l2j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j4.png\"><img src=\"plots/dR_l3j4.png\" alt=\"plots/dR_l3j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j4.png\"><img src=\"plots/dR_l4j4.png\" alt=\"plots/dR_l4j4.png\" width=\"100%\"></a></td>" << endl;
-<
+  // htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA.png\"><img src=\"plots/fusionMVA.png\" alt=\"plots/fusionMVA.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_lo.png\"><img src=\"plots/fusionMVA_lo.png\" alt=\"plots/fusionMVA_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_hi.png\"><img src=\"plots/fusionMVA_hi.png\" alt=\"plots/fusionMVA_hi.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA-log.png\"><img src=\"plots/fusionMVA-log.png\" alt=\"plots/fusionMVA-log.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_lo-log.png\"><img src=\"plots/fusionMVA_lo-log.png\" alt=\"plots/fusionMVA_lo-log.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_hi-log.png\"><img src=\"plots/fusionMVA_hi-log.png\" alt=\"plots/fusionMVA_hi-log.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
-<
-<
+  htmlfile << "<tr>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_jet_lo.png\"><img src=\"plots/m4l_jet_lo.png\" alt=\"plots/m4l_jet_lo.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_jet.png\"><img src=\"plots/m4l_jet.png\" alt=\"plots/m4l_jet.png\" width=\"100%\"></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
-<
+  htmlfile << "</tr>" << endl;
->
+//   htmlfile << "</body>" << endl;
->
+//   htmlfile << "</html>" << endl;
->
+//   htmlfile.close();
->
->
+//   TString jetfname(htmlfname);
->
+//   jetfname.ReplaceAll("plots.html","jetplots.html");
->
+//   htmlfile.open(jetfname);
->
->
+//   htmlfile << "<!DOCTYPE html" << endl;
->
+//   htmlfile << "    PUBLIC \"-//W3C//DTD HTML 3.2//EN\">" << endl;
->
+//   htmlfile << "<html>" << endl;
->
->
+//   htmlfile << "<head><title>"+title+"</title></head>" << endl;
->
+//   htmlfile << "<body bgcolor=\"000000\">" << endl;
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">"+title+"</h3>" << endl;
->
->
+//   htmlfile << "<h3 style=\"text-align:left; color:DD6600;\">jet plots</h3>" << endl;
->
+//   htmlfile << "<hr />" << endl;
->
+//   htmlfile << "<table border=\"0\" cellspacing=\"5\" width=\"100%\">" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j1.png\"><img src=\"plots/pt_j1.png\" alt=\"plots/pt_j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j2.png\"><img src=\"plots/pt_j2.png\" alt=\"plots/pt_j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j2.png\"><img src=\"plots/pt_j2.png\" alt=\"plots/pt_j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/pt_j4.png\"><img src=\"plots/pt_j4.png\" alt=\"plots/pt_j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j1.png\"><img src=\"plots/eta_j1.png\" alt=\"plots/eta_j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j2.png\"><img src=\"plots/eta_j2.png\" alt=\"plots/eta_j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j3.png\"><img src=\"plots/eta_j3.png\" alt=\"plots/eta_j3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/eta_j4.png\"><img src=\"plots/eta_j4.png\" alt=\"plots/eta_j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j1.png\"><img src=\"plots/idmva_j1.png\" alt=\"plots/idmva_j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j2.png\"><img src=\"plots/idmva_j2.png\" alt=\"plots/idmva_j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j3.png\"><img src=\"plots/idmva_j3.png\" alt=\"plots/idmva_j3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/idmva_j4.png\"><img src=\"plots/idmva_j4.png\" alt=\"plots/idmva_j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j1.png\"><img src=\"plots/phi_j1.png\" alt=\"plots/phi_j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j2.png\"><img src=\"plots/phi_j2.png\" alt=\"plots/phi_j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j3.png\"><img src=\"plots/phi_j3.png\" alt=\"plots/phi_j3.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/phi_j4.png\"><img src=\"plots/phi_j4.png\" alt=\"plots/phi_j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/mjj.png\"><img src=\"plots/mjj.png\" alt=\"plots/mjj.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dEta.png\"><img src=\"plots/dEta.png\" alt=\"plots/dEta.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/etaProd.png\"><img src=\"plots/etaProd.png\" alt=\"plots/etaProd.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/njets.png\"><img src=\"plots/njets.png\" alt=\"plots/njets.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   // htmlfile << "<tr>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j1.png\"><img src=\"plots/dR_l1j1.png\" alt=\"plots/dR_l1j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j1.png\"><img src=\"plots/dR_l2j1.png\" alt=\"plots/dR_l2j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j1.png\"><img src=\"plots/dR_l3j1.png\" alt=\"plots/dR_l3j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j1.png\"><img src=\"plots/dR_l4j1.png\" alt=\"plots/dR_l4j1.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "</tr>" << endl;
->
->
+//   // htmlfile << "<tr>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j2.png\"><img src=\"plots/dR_l1j2.png\" alt=\"plots/dR_l1j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j2.png\"><img src=\"plots/dR_l2j2.png\" alt=\"plots/dR_l2j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j2.png\"><img src=\"plots/dR_l3j2.png\" alt=\"plots/dR_l3j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j2.png\"><img src=\"plots/dR_l4j2.png\" alt=\"plots/dR_l4j2.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "</tr>" << endl;
->
->
+//   // htmlfile << "<tr>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j3.png\"><img src=\"plots/dR_l1j3.png\" alt=\"plots/dR_l1j3.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j3.png\"><img src=\"plots/dR_l2j3.png\" alt=\"plots/dR_l2j3.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j3.png\"><img src=\"plots/dR_l3j3.png\" alt=\"plots/dR_l3j3.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j3.png\"><img src=\"plots/dR_l4j3.png\" alt=\"plots/dR_l4j3.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "</tr>" << endl;
->
->
+//   // htmlfile << "<tr>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l1j4.png\"><img src=\"plots/dR_l1j4.png\" alt=\"plots/dR_l1j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l2j4.png\"><img src=\"plots/dR_l2j4.png\" alt=\"plots/dR_l2j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l3j4.png\"><img src=\"plots/dR_l3j4.png\" alt=\"plots/dR_l3j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/dR_l4j4.png\"><img src=\"plots/dR_l4j4.png\" alt=\"plots/dR_l4j4.png\" width=\"100%\"></a></td>" << endl;
->
+//   // htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA.png\"><img src=\"plots/fusionMVA.png\" alt=\"plots/fusionMVA.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_lo.png\"><img src=\"plots/fusionMVA_lo.png\" alt=\"plots/fusionMVA_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_hi.png\"><img src=\"plots/fusionMVA_hi.png\" alt=\"plots/fusionMVA_hi.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA-log.png\"><img src=\"plots/fusionMVA-log.png\" alt=\"plots/fusionMVA-log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_lo-log.png\"><img src=\"plots/fusionMVA_lo-log.png\" alt=\"plots/fusionMVA_lo-log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/fusionMVA_hi-log.png\"><img src=\"plots/fusionMVA_hi-log.png\" alt=\"plots/fusionMVA_hi-log.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
->
->
+//   htmlfile << "<tr>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_jet_lo.png\"><img src=\"plots/m4l_jet_lo.png\" alt=\"plots/m4l_jet_lo.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a target=\"_blank\" href=\"plots/m4l_jet.png\"><img src=\"plots/m4l_jet.png\" alt=\"plots/m4l_jet.png\" width=\"100%\"></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "<td width=\"25%\"><a></a></td>" << endl;
->
+//   htmlfile << "</tr>" << endl;
-<
+  htmlfile << "</table>" << endl;
->
+//   htmlfile << "</table>" << endl;
-<
+  htmlfile << "<hr />" << endl;
->
+//   htmlfile << "<hr />" << endl;
-<
+  htmlfile << "<hr />" << endl;
->
+//   htmlfile << "<hr />" << endl;
-<
+  htmlfile << "</body>" << endl;
-<
+  htmlfile << "</html>" << endl;
-<
+  htmlfile.close();
-<
+}
-<
+//----------------------------------------------------------------------------------------
-<
+void fillHist(CSample *cs, TString channel, TString type, TString var, double val, double wgt, double wgt_lo, double wgt_hi)
-<
+{
-<
+  (*(cs->hists)["all_"+type])[var]->Fill(           val,   wgt);
-<
+  if(channel!="")
-<
+     (*(cs->hists)[channel+"_"+type])[var]->Fill(   val,   wgt);
-<
+  if(type=="pred" || type=="PF") {
-<
+    (*(cs->hists)["all_"+type+"_lo"])[var]->Fill(   val,   wgt_lo);
-<
+    (*(cs->hists)["all_"+type+"_hi"])[var]->Fill(   val,   wgt_hi);
-<
+    if(channel!="") {
-<
+      (*(cs->hists)[channel+"_"+type+"_lo"])[var]->Fill(   val,   wgt_lo);
-<
+      (*(cs->hists)[channel+"_"+type+"_hi"])[var]->Fill(   val,   wgt_hi);
-<
+    }
-<
+  }
-<
+}
-<
+//----------------------------------------------------------------------------------------
-<
+void fillAllHists(FOFlags ctrl, CSample *cs, TString channel, TString type, filestuff *fs, KinematicsStruct kine,
-<
+                  SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4,
-<
+                  double wgt, double wgt_lo, double wgt_hi)
-<
+{
-<
+  fillHist( cs, channel, type, "npv"         , fs->info->npv  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "run"         , fs->info->run  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "mZ1"         , kine.mZ1       , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "mZ2"         , kine.mZ2       , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "mZ2_lo"      , kine.mZ2       , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "m4l_lo"      , kine.m4l       , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "m4l_med"     , kine.m4l       , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "m4l"         , kine.m4l       , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "Z1pt"        , kine.Z1pt      , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "ZZpt"        , kine.ZZpt      , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "met"         , fs->info->met  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "ip3ds_l1"    , lep1.ip3dSig   , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "ip3ds_l2"    , lep2.ip3dSig   , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "ip3ds_l3"    , lep3.ip3dSig   , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "ip3ds_l4"    , lep4.ip3dSig   , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "ip3ds_l3_lo" , lep3.ip3dSig   , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "ip3ds_l4_lo" , lep4.ip3dSig   , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "d0_l1"       , lep1.d0        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "d0_l2"       , lep2.d0        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "d0_l3"       , lep3.d0        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "d0_l4"       , lep4.d0        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "d0_l3_lo"    , lep3.d0        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "d0_l4_lo"    , lep4.d0        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dz_l1"       , lep1.dz        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dz_l2"       , lep2.dz        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dz_l3"       , lep3.dz        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dz_l4"       , lep4.dz        , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "pt_l1"       , lep1.vec.Pt()  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "pt_l2"       , lep2.vec.Pt()  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "pt_l3"       , lep3.vec.Pt()  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "pt_l4"       , lep4.vec.Pt()  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_l1"      , lep1.vec.Eta() , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_l2"      , lep2.vec.Eta() , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_l3"      , lep3.vec.Eta() , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_l4"      , lep4.vec.Eta() , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l3l4_lo"  , dr(lep3,lep4)  , wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l3l4"     , dr(lep3,lep4)  , wgt, wgt_lo, wgt_hi);
-<
+}
-<
+//----------------------------------------------------------------------------------------
-<
+void fillAllJetHists( FOFlags ctrl, CSample *cs, TString channel, TString type, filestuff *fs, KinematicsStruct kine,
-<
+                      SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4,
-<
+                      FusionMva &fusion, vector<SimpleLepton> &goodJets, JetInfoStruct ji,
-<
+                      double wgt, double wgt_lo, double wgt_hi)
-<
+{
-<
+  double fMVAval = (ctrl.uncert=="") ? 0 : fusion.reader->EvaluateMVA("BDTG");
-<
+  dr_struct drs = fill_dr_struct(goodJets,lep1,lep2,lep3,lep4);
-<
+  drs.check();
-<
+  fillHist( cs, channel, type, "pt_j1"  ,     ji.ptJet1,        wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "pt_j2"  ,     ji.ptJet2,        wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "pt_j3"  ,     ji.ptJet3,        wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "pt_j4"  ,     ji.ptJet4,        wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_j1" ,     ji.etaJet1,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_j2" ,     ji.etaJet2,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_j3" ,     ji.etaJet3,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "eta_j4" ,     ji.etaJet4,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "phi_j1" ,     ji.phiJet1,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "phi_j2" ,     ji.phiJet2,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "phi_j3" ,     ji.phiJet3,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "phi_j4" ,     ji.phiJet4,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "idmva_j1",    ji.mvaJet1,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "idmva_j2",    ji.mvaJet2,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "idmva_j3",    ji.mvaJet3,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "idmva_j4",    ji.mvaJet4,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "mjj"    ,     ji.mjj,           wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dEta"   ,     ji.dEta,          wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "etaProd",     ji.etaProd,       wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "njets"    ,   goodJets.size(),  wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l1j1",     drs.l1j1,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l2j1",     drs.l2j1,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l3j1",     drs.l3j1,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l4j1",     drs.l4j1,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l1j2",     drs.l1j2,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l2j2",     drs.l2j2,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l3j2",     drs.l3j2,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l4j2",     drs.l4j2,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l1j3",     drs.l1j3,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l2j3",     drs.l2j3,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l3j3",     drs.l3j3,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l4j3",     drs.l4j3,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l1j4",     drs.l1j4,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l2j4",     drs.l2j4,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l3j4",     drs.l3j4,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "dR_l4j4",     drs.l4j4,         wgt, wgt_lo, wgt_hi);
-<
-<
+  fillHist( cs, channel, type, "fusionMVA_lo", fMVAval,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "fusionMVA_med",fMVAval,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "fusionMVA_hi", fMVAval,         wgt, wgt_lo, wgt_hi);
-<
+  fillHist( cs, channel, type, "fusionMVA",    fMVAval,         wgt, wgt_lo, wgt_hi);
-<
+  if(fMVAval > -.3) {
-<
+    fillHist( cs, channel, type, "m4l_jet_lo",  kine.m4l  ,     wgt, wgt_lo, wgt_hi);
-<
+    fillHist( cs, channel, type, "m4l_jet",     kine.m4l  ,     wgt, wgt_lo, wgt_hi);
-<
+  }
-<
+}
->
+//   htmlfile << "</body>" << endl;
->
+//   htmlfile << "</html>" << endl;
->
+//   htmlfile.close();
->
+// }
+//----------------------------------------------------------------------------------------
+bool passHlt(FOFlags &ctrl, TrigInfo ti, InfoStruct *info, unsigned lep1matchBits,
+             unsigned lep2matchBits, unsigned lep3matchBits,
+  return pass;
+//----------------------------------------------------------------------------------------
-–
+void fillLeptonFakeWeights(lepFrs &fwgts, FR_struct *fr, SimpleLepton lep3, SimpleLepton lep4)
-–
+{
-–
+  fwgts.fwgt_3    = get_fake_weight("",lep3,*fr);
-–
+  fwgts.fwgt_lo_3 = get_fake_weight("lo",lep3,*fr);
-–
+  fwgts.fwgt_hi_3 = get_fake_weight("hi",lep3,*fr);
-–
+  fwgts.fwgt_4    = get_fake_weight("",lep4,*fr);
-–
+  fwgts.fwgt_lo_4 = get_fake_weight("lo",lep4,*fr);
-–
+  fwgts.fwgt_hi_4 = get_fake_weight("hi",lep4,*fr);
-–
+}
-–
+//----------------------------------------------------------------------------------------
-–
+dr_struct fill_dr_struct(vector<SimpleLepton> jets, SimpleLepton lep1, SimpleLepton lep2, SimpleLepton lep3, SimpleLepton lep4)
-–
+{
-–
+  dr_struct drs;
-–
+  if(jets.size() > 0) {
-–
+    drs.l1j1 = dr(lep1,jets[0]);
-–
+    drs.l2j1 = dr(lep2,jets[0]);
-–
+    drs.l3j1 = dr(lep3,jets[0]);
-–
+    drs.l4j1 = dr(lep4,jets[0]);
-–
+  }
-–
+  if(jets.size() > 1) {
-–
+    drs.l1j2 = dr(lep1,jets[1]);
-–
+    drs.l2j2 = dr(lep2,jets[1]);
-–
+    drs.l3j2 = dr(lep3,jets[1]);
-–
+    drs.l4j2 = dr(lep4,jets[1]);
-–
+  }
-–
+  if(jets.size() > 2) {
-–
+    drs.l1j3 = dr(lep1,jets[2]);
-–
+    drs.l2j3 = dr(lep2,jets[2]);
-–
+    drs.l3j3 = dr(lep3,jets[2]);
-–
+    drs.l4j3 = dr(lep4,jets[2]);
-–
+  }
-–
+  if(jets.size() > 3) {
-–
+    drs.l1j4 = dr(lep1,jets[3]);
-–
+    drs.l2j4 = dr(lep2,jets[3]);
-–
+    drs.l3j4 = dr(lep3,jets[3]);
-–
+    drs.l4j4 = dr(lep4,jets[3]);
-–
+  }
-–
+  return drs;
-–
+}
-–
+//----------------------------------------------------------------------------------------
+void init_cuts( vector<TString> &cutstrs, map<TString,int> &cutvec)
+  cutstrs.push_back("start"  );    cutvec["start"]  = 0;
+  return hset;
+//----------------------------------------------------------------------------------------
-–
+void shiftOverflows(map<TString,TH1D*> &hset)
-–
+{
-–
+  for(map<TString,TH1D*>::iterator it=hset.begin(); it!=hset.end(); it++)
-–
+    shiftOverflows((*it).second);
-–
+}
-–
+//----------------------------------------------------------------------------------------
-–
+void scaleHists(map<TString,TH1D*> &hset)
-–
+{
-–
+  for(map<TString,TH1D*>::iterator it=hset.begin(); it!=hset.end(); it++)
-–
+    (*it).second->Scale(1./integrateHist((*it).second));
-–
+}
-–
+//----------------------------------------------------------------------------------------
+void fillFakeTuple(CSample *cs, filestuff *fs, EventData evtdata, unsigned ientry, double fillweight)
+  fs->getentry(ientry,"","zznt"); // make sure we've got this entry for all the branches
-+
+  fillAngles(evtdata,*fs->angles);
+  fillKinematics(evtdata,*fs->kine); // reminder: kine has the fake information, while fs->kine *had*, until this line, whatever ZPlusX put in it
+  fs->weights->w = fillweight;
+  fs->outtuple->Fill();
+  fs->outFotuple->Fill();
-+
+//----------------------------------------------------------------------------------------
-+
+void incrementSsofCounters(FOFlags &ctrl, filestuff *fs, double minMz2,
-+
+                           pair<int,int> *best_z_indices,
-+
+                           double wgt)
-+
+{
-+
+  assert(ctrl.ssof); // only makes sense to write this out when we're thinking about ssof method
-+
+  vector<pair<SimpleLepton,SimpleLepton> > z2cands;
-+
+  vector<TString> types;
-+
-+
+  int hiPtFailZ2 = findZ2CandidatesPassFail(fs->failingL, z2cands, types, "XX-SF", minMz2, "fail");
-+
+  if(hiPtFailZ2 == -1) return; // no z2 found
-+
+  SimpleLepton lep1 = (*fs->passingL)[best_z_indices->first];
-+
+  SimpleLepton lep2 = (*fs->passingL)[best_z_indices->second];
-+
+  SimpleLepton lep3 = z2cands[hiPtFailZ2].first;
-+
+  SimpleLepton lep4 = z2cands[hiPtFailZ2].second;
-+
+  TString channel = getChannel(lep1,lep2,lep3,lep4);
-+
+  TString sign = (lep3.charge == lep4.charge) ? "SS" : "OS";
-<
+          // double mjj=0,dEta=0,etaProd=0,nCentralJets=0;
-<
+          // if(jet1 && jet2) {
-<
+          //   cutvec["twoJets"] += 1;
-<
+          //   if(ctrl.debug) cout << "looping through " << goodJets.size() << " good jets" << endl;
-<
+          //   if(ctrl.debug) {cout << "jet 1 (" << jet1 << "): "; jet1->print();}
-<
+          //   if(ctrl.debug) {cout << "jet 2 (" << jet2 << "): "; jet2->print();}
-<
+          //   for(unsigned ijet=0; ijet<goodJets.size(); ijet++) {
-<
+          //     SimpleLepton *jet = goodJets[ijet];
-<
+          //     if(jet == jet1) {
-<
+          //    if(ctrl.debug) cout << "  " << jet << " this is jet 1" << endl;
-<
+          //    continue;
-<
+          //     }
-<
+          //     if(jet == jet2) {
-<
+          //    if(ctrl.debug) cout << "  " << jet << " this is jet 2" << endl;
-<
+          //    continue;
-<
+          //     }
-<
+          //     double eta = jet->vec.Eta();
-<
+          //     double eta1 = jet1->vec.Eta();
-<
+          //     double eta2 = jet2->vec.Eta();
-<
+          //     if(eta1 > eta2) {
-<
+          //    if(eta > eta2 && eta < eta1)
-<
+          //      nCentralJets++;
-<
+          //     } else if(eta2 > eta1) {
-<
+          //    if(eta > eta1 && eta < eta2)
-<
+          //      nCentralJets++;
-<
+          //     }
-<
+          //     if(ctrl.debug) cout << "  nCentralJets: " << nCentralJets << endl;
-<
+          //   }
-<
+          //   double mjjMin = 400;
-<
+          //   double dEtaMin = 4;
-<
+          //   TLorentzVector dijet(jet1->vec + jet2->vec);
-<
+          //   mjj = dijet.M();
-<
+          //   if(ctrl.debug) cout << "setting mjj: " << mjj << "(using " << jet1->vec.Pt() << " " << jet2->vec.Pt() << endl;
-<
+          //   dEta = jet1->vec.Eta() - jet2->vec.Eta();
-<
+          //   etaProd = jet1->vec.Eta()*jet2->vec.Eta();
-<
+          //   bool isVbf = (mjj > mjjMin) && (fabs(dEta) > dEtaMin) && (etaProd < 0) && (nCentralJets==0);
-<
+          // }
-<
->
+  fs->counters_[sign+"_"+channel] += wgt;
->
+  fs->counters_[sign] += wgt;
->
+}
->
+//----------------------------------------------------------------------------------------
->
+double getSsofWgt(TString channel, map<TString,double> &ssofRatios)
->
+{
->
+  double returnval = ssofRatios[channel];
->
+  if(returnval < 0.1 || returnval > 3) {
->
+    cout << "\nERROR! ssof ratio out of bounds for " << channel << ": " << returnval << endl;
->
+    assert(0);
->
+  }
->
+  return returnval;
->
+}
->
+//----------------------------------------------------------------------------------------
->
+map<TString,double> initSsofRatios(vector<CSample*> &samplev)
->
+{
->
+  double OS,SS;
->
+  double OS_4e,SS_4e;
->
+  double OS_4m,SS_4m;
->
+  double OS_2e2m,SS_2e2m;
->
+  cout << "Initializing ssof ratios: " << endl;
->
+  for(unsigned isam=0; isam<samplev.size(); isam++) {
->
+    CSample *cs = samplev[isam];
->
+    for(unsigned ifs=0; ifs<(cs->fsv).size(); ifs++) {
->
+      filestuff *fs = (cs->fsv)[ifs];
->
+      if(fs->useSsofRatio_) {
->
+        cout << "  adding from: " << fs->fname_ << endl;
->
+        cout << "        " << fs->counters_["ratio"] << setw(8) << fixed << setprecision(2) << fs->counters_["ratio_4e"] << setw(8) << fixed << setprecision(2) << fs->counters_["ratio_4m"] << setw(8) << fixed << setprecision(2) << fs->counters_["ratio_2e2m"] << endl;
->
+        OS      += fs->counters_["OS"];
->
+        SS      += fs->counters_["SS"];
->
+        OS_4e   += fs->counters_["OS_4e"];
->
+        SS_4e   += fs->counters_["SS_4e"];
->
+        OS_4m   += fs->counters_["OS_4m"];
->
+        SS_4m   += fs->counters_["SS_4m"];
->
+        OS_2e2m += fs->counters_["OS_2e2m"];
->
+        SS_2e2m += fs->counters_["SS_2e2m"];
->
+      }
->
+    }
->
+  }
->
+  assert(OS!=0 && SS!=0 && OS_4e!=0 && SS_4e!=0 && OS_4m!=0 && SS_4m!=0 && OS_2e2m && SS_2e2m!=0);
->
+  map<TString,double> ssofRatios;
->
+  ssofRatios[""]     = OS/SS;
->
+  ssofRatios["4e"]   = OS_4e/SS_4e;
->
+  ssofRatios["4m"]   = OS_4m/SS_4m;
->
+  ssofRatios["2e2m"] = OS_2e2m/SS_2e2m;
->
+  cout << "returning: " << ssofRatios[""] << setw(8) << fixed << setprecision(2) << ssofRatios["4e"] << setw(8) << fixed << setprecision(2) << ssofRatios["4m"] << setw(8) << fixed << setprecision(2) << ssofRatios["2e2m"] << endl;
->
+  return ssofRatios;
->
+}
->

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing UserCode/MitHzz4l/NonMCBackground/src/plotSignalFakes.cc (file contents): Revision 1.4 by dkralph, Tue Jul 31 16:46:22 2012 UTC vs. Revision 1.7 by dkralph, Mon Dec 17 17:25:55 2012 UTC

Diff Legend

Comparing UserCode/MitHzz4l/NonMCBackground/src/plotSignalFakes.cc (file contents):
Revision 1.4 by dkralph, Tue Jul 31 16:46:22 2012 UTC vs.
Revision 1.7 by dkralph, Mon Dec 17 17:25:55 2012 UTC