[ViewVC] Diff of: cvsroot/UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents):
Revision 1.70 by qpython, Wed May 29 15:29:27 2013 UTC vs.
Revision 1.75 by wulsin, Fri May 31 09:15:04 2013 UTC

+  datasetType_ (cfg.getParameter<string> ("datasetType")),
+  channels_  (cfg.getParameter<vector<edm::ParameterSet> >("channels")),
+  histogramSets_ (cfg.getParameter<vector<edm::ParameterSet> >("histogramSets")),
-+
+  treeBranchSets_   (cfg.getParameter<vector<edm::ParameterSet> >("treeBranchSets")),
+  plotAllObjectsInPassingEvents_ (cfg.getParameter<bool> ("plotAllObjectsInPassingEvents")),
+  doPileupReweighting_ (cfg.getParameter<bool> ("doPileupReweighting")),
+  applyLeptonSF_ (cfg.getParameter<bool> ("applyLeptonSF")),
+  //always get vertex collection so we can assign the primary vertex in the event
+  objectsToGet.push_back("primaryvertexs");
-–
-–
+  //always make the plot of number of primary vertices (to check pile-up reweighting)
+  objectsToPlot.push_back("primaryvertexs");
-+
+  objectsToCut.push_back("primaryvertexs");
-+
+  //always get the MC particles to do GEN-matching
+  objectsToGet.push_back("mcparticles");
+  //always get the event collection to do pile-up reweighting
+  objectsToGet.push_back("events");
-+
-+
+  // Parse the tree variable definitions.
-+
+  for (uint iBranchSet = 0; iBranchSet<treeBranchSets_.size(); iBranchSet++) {
-+
+    string tempInputCollection = treeBranchSets_.at(iBranchSet).getParameter<string> ("inputCollection");
-+
+    if(tempInputCollection.find("pairs")!=string::npos) { cout << "Warning:  tree filling is not configured for pairs of objects, so will not work for collection: " << tempInputCollection << endl; }
-+
+    objectsToGet.push_back(tempInputCollection);
-+
+    objectsToCut.push_back(tempInputCollection);
-+
-+
+    //    vector<edm::ParameterSet> branchList(treeBranchSets_.at(iBranchSet).getParameter<vector<edm::ParameterSet> >("histograms"));
-+
+    vector<string> branchList(treeBranchSets_.at(iBranchSet).getParameter<vector<string> >("branches"));
-+
-+
+    for(uint iBranch = 0; iBranch<branchList.size(); iBranch++){
-+
-+
+      vector<double> defaultValue;
-+
+      defaultValue.push_back (-1.0);
-+
-+
+      BranchSpecs br;
-+
+      br.inputCollection = tempInputCollection;
-+
+      br.inputVariable = branchList.at(iBranch);
-+
+      TString newName = TString(br.inputCollection) + "_" + TString(br.inputVariable);
-+
+      br.name = string(newName.Data());
-+
-+
+      treeBranches_.push_back(br);
-+
-+
+    }
-+
-+
+  } // end   for (uint iBranchSet = 0; iBranchSet<treeBranchSets_.size(); iBranchSet++) {
-+
-+
+  //parse the histogram definitions
+  for(uint currentHistogramSet = 0; currentHistogramSet != histogramSets_.size(); currentHistogramSet++){
+      objectsToPlot.push_back(tempInputCollection);
+      objectsToCut.push_back(tempInputCollection);
-<
+    }
-<
+    else{//pair of objects, need to add the pair and the individual objects to the lists of things to Get/Plot/Cut
->
+    } else { //pair of objects, need to add the pair and the individual objects to the lists of things to Get/Plot/Cut
+      string obj1;
+      string obj2;
+      getTwoObjs(tempInputCollection, obj1, obj2);
+      objectsToGet.push_back(tempInputCollection);
+      objectsToGet.push_back(obj1);
+      objectsToGet.push_back(obj2ToGet);
-+
+    } // end else
-+
-–
+    }
+    vector<edm::ParameterSet> histogramList_  (histogramSets_.at(currentHistogramSet).getParameter<vector<edm::ParameterSet> >("histograms"));
+      histograms.push_back(tempHistogram);
-<
+  }
->
+  } //   for(uint currentHistogramSet = 0; currentHistogramSet != histogramSets_.size(); currentHistogramSet++){
->
+  //make unique vector of objects we need to plot (so we can book a histogram with the number of each object)
+  sort( objectsToPlot.begin(), objectsToPlot.end() );
+  objectsToPlot.erase( unique( objectsToPlot.begin(), objectsToPlot.end() ), objectsToPlot.end() );
+    //set triggers for this channel
+    vector<string> triggerNames;
+    triggerNames.clear();
-+
+    vector<string> triggerToVetoNames;
-+
+    triggerToVetoNames.clear();
-+
+    tempChannel.triggers.clear();
-+
+    tempChannel.triggersToVeto.clear();
+    if(channels_.at(currentChannel).exists("triggers")){
+      triggerNames   = channels_.at(currentChannel).getParameter<vector<string> >("triggers");
+      for(uint trigger = 0; trigger!= triggerNames.size(); trigger++)
+        tempChannel.triggers.push_back(triggerNames.at(trigger));
+      objectsToGet.push_back("triggers");
-<
->
+    if(channels_.at(currentChannel).exists("triggersToVeto")){
->
+      triggerToVetoNames = channels_.at(currentChannel).getParameter<vector<string> >("triggersToVeto");
->
+      for(uint trigger = 0; trigger!= triggerToVetoNames.size(); trigger++)
->
+        tempChannel.triggersToVeto.push_back(triggerToVetoNames.at(trigger));
->
+      objectsToGet.push_back("triggers");
->
+    }
+    TFileDirectory subDir = fs_->mkdir( channelName );
+    //loop over the cuts to set up subdirectory for each cut if GetPlotsAfterEachCut_ is true.
+    if(GetPlotsAfterEachCut_){
-<
+       for( uint currentDir=0; currentDir != subSubDirNames.size(); currentDir++){
-<
+            TFileDirectory subSubDir = subDir.mkdir( subSubDirNames[currentDir] );
-<
+            directories.push_back(subSubDir);
-<
+            map<string, TH1D*> oneDhistoMap;
-<
+            oneDHistsTmp.push_back(oneDhistoMap);
-<
+            map<string, TH2D*> twoDhistoMap;
-<
+            twoDHistsTmp.push_back(twoDhistoMap);
->
+      for( uint currentDir=0; currentDir != subSubDirNames.size(); currentDir++){
->
+        TFileDirectory subSubDir = subDir.mkdir( subSubDirNames[currentDir] );
->
+        directories.push_back(subSubDir);
->
+        map<string, TH1D*> oneDhistoMap;
->
+        oneDHistsTmp.push_back(oneDhistoMap);
->
+        map<string, TH2D*> twoDhistoMap;
->
+        twoDHistsTmp.push_back(twoDhistoMap);
+      oneDHists_.push_back(oneDHistsTmp);
+      twoDHists_.push_back(twoDHistsTmp);
-<
+   //only set up directories with names of the channels if GetPlotsAfterEachCut_ is false.
-<
+   else{
->
+    //only set up directories with names of the channels if GetPlotsAfterEachCut_ is false.
->
+    else{
+      map<string, TH1D*> oneDhistoMap;
+      oneDHistsTmp.push_back(oneDhistoMap);
+      map<string, TH2D*> twoDhistoMap;
+      oneDHists_.push_back(oneDHistsTmp);
+      twoDHists_.push_back(twoDHistsTmp);
+      directories.push_back(subDir);
-–
+   }
-–
+    //book all histograms included in the configuration
-–
+    for(uint currentDir = 0; currentDir != directories.size(); currentDir++){//loop over all the directories.
-–
+    for(uint currentHistogramIndex = 0; currentHistogramIndex != histograms.size(); currentHistogramIndex++){
-–
+      histogram currentHistogram = histograms.at(currentHistogramIndex);
-–
+      int numBinsElements = currentHistogram.bins.size();
-–
+      int numInputVariables = currentHistogram.inputVariables.size();
-–
+      int numBinEdgesX = currentHistogram.variableBinsX.size();
-–
+      int numBinEdgesY = currentHistogram.variableBinsY.size();
-–
-–
+      if(numBinsElements == 1){
-–
+        if(numBinEdgesX > 1){
-–
+          if(numBinEdgesY > 1)
-–
+            numBinsElements = 6;
-–
+          else
-–
+            numBinsElements = 3;
-–
+        }
-–
+      }
-–
+      if(numBinsElements != 3 && numBinsElements !=6){
-–
+        cout << "Error: Didn't find correct number of bin specifications for histogram named '" << currentHistogram.name << "'\n";
-–
+        exit(0);
-–
+      }
-–
+      else if((numBinsElements == 3 && numInputVariables !=1) || (numBinsElements == 6 && numInputVariables!=2)){
-–
+        cout << "Error: Didn't find correct number of input variables for histogram named '" << currentHistogram.name << "'\n";
-–
+        exit(0);
-–
+      }
-–
+      else if(numBinsElements == 3){
-–
+        if (currentHistogram.bins.size () == 3)
-–
+          oneDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH1D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, currentHistogram.bins.at(0), currentHistogram.bins.at(1), currentHistogram.bins.at(2));
-–
+        else
-–
+          {
-–
+            oneDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH1D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, numBinEdgesX - 1, currentHistogram.variableBinsX.data ());
-–
+          }
-–
+      }
-–
+      else if(numBinsElements == 6){
-–
+        if (currentHistogram.bins.size () == 6)
-–
+          twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH2D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, currentHistogram.bins.at(0), currentHistogram.bins.at(1), currentHistogram.bins.at(2),currentHistogram.bins.at(3),currentHistogram.bins.at(4),currentHistogram.bins.at(5));
-–
+        else
-–
+          twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH2D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, numBinEdgesX - 1, currentHistogram.variableBinsX.data (), numBinEdgesY - 1, currentHistogram.variableBinsY.data ());
-–
+      }
-–
-–
-–
+      if(currentHistogram.name.find("GenMatch")==string::npos) continue;
-–
-–
+      // bin      particle type
-–
+      // ---      -------------
-–
+      //  0        unmatched
-–
+      //  1        u
-–
+      //  2        d
-–
+      //  3        s
-–
+      //  4        c
-–
+      //  5        b
-–
+      //  6        t
-–
+      //  7        e
-–
+      //  8        mu
-–
+      //  9        tau
-–
+      // 10        nu
-–
+      // 11        g
-–
+      // 12        gamma
-–
+      // 13        Z
-–
+      // 14        W
-–
+      // 15        light meson
-–
+      // 16        K meson
-–
+      // 17        D meson
-–
+      // 18        B meson
-–
+      // 19        light baryon
-–
+      // 20        strange baryon
-–
+      // 21        charm baryon
-–
+      // 22        bottom baryon
-–
+      // 23        QCD string
-–
+      // 24        other
-–
-–
+      vector<TString> labelArray;
-–
+      labelArray.push_back("unmatched");
-–
+      labelArray.push_back("u");
-–
+      labelArray.push_back("d");
-–
+      labelArray.push_back("s");
-–
+      labelArray.push_back("c");
-–
+      labelArray.push_back("b");
-–
+      labelArray.push_back("t");
-–
+      labelArray.push_back("e");
-–
+      labelArray.push_back("#mu");
-–
+      labelArray.push_back("#tau");
-–
+      labelArray.push_back("#nu");
-–
+      labelArray.push_back("g");
-–
+      labelArray.push_back("#gamma");
-–
+      labelArray.push_back("Z");
-–
+      labelArray.push_back("W");
-–
+      labelArray.push_back("light meson");
-–
+      labelArray.push_back("K meson");
-–
+      labelArray.push_back("D meson");
-–
+      labelArray.push_back("B meson");
-–
+      labelArray.push_back("light baryon");
-–
+      labelArray.push_back("strange baryon");
-–
+      labelArray.push_back("charm baryon");
-–
+      labelArray.push_back("bottom baryon");
-–
+      labelArray.push_back("QCD string");
-–
+      labelArray.push_back("other");
-–
-–
+      for(int bin = 0; bin !=currentHistogram.bins.at(0); bin++){
-–
+        if(currentHistogram.name.find("GenMatchIdVsMotherId")==string::npos && currentHistogram.name.find("GenMatchIdVsGrandmotherId")==string::npos) {
-–
+          oneDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->SetBinLabel(bin+1,labelArray.at(bin));
-–
+        }
-–
+        else {
-–
+          twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetYaxis()->SetBinLabel(bin+1,labelArray.at(currentHistogram.bins.at(0)-bin-1));
-–
+          twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->SetBinLabel(bin+1,labelArray.at(bin));
-–
+        }
-–
+      }
-–
+      if(currentHistogram.name.find("GenMatchIdVsMotherId")!=string::npos || currentHistogram.name.find("GenMatchIdVsGrandmotherId")!=string::npos) {
-–
+        twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->CenterTitle();
-–
+        twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetYaxis()->CenterTitle();
-–
+      }
-<
+    // Book a histogram for the number of each object type to be plotted.
-<
+    // Name of objectToPlot here must match the name specified in OSUAnalysis::analyze().
-<
+    for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++){
-<
+      string currentObject = objectsToPlot.at(currentObjectIndex);
-<
+      int maxNum = 10;
-<
+      if(currentObject == "mcparticles") maxNum = 50;
-<
+      else if(currentObject == "primaryvertexs") maxNum = 50;
-<
-<
+      if(currentObject == "muon-muon pairs")                currentObject = "dimuonPairs";
-<
+      else if(currentObject == "electron-electron pairs")   currentObject = "dielectronPairs";
-<
+      else if(currentObject == "electron-muon pairs")       currentObject = "electronMuonPairs";
-<
+      else if(currentObject == "secondary jets")            currentObject = "secondaryJets";
-<
+      else if(currentObject == "jet-jet pairs")             currentObject = "dijetPairs";
-<
+      else if(currentObject == "electron-jet pairs")        currentObject = "electronJetPairs";
-<
+      else if(currentObject == "muon-jet pairs")            currentObject = "muonJetPairs";
-<
+      else if(currentObject == "track-event pairs")         currentObject = "trackEventPairs";
-<
+      else if(currentObject == "electron-track pairs")      currentObject = "electronTrackPairs";
-<
+      else if(currentObject == "muon-track pairs")          currentObject = "muonTrackPairs";
-<
+      else if(currentObject == "muon-tau pairs")            currentObject = "muonTauPairs";
-<
+      else if(currentObject == "tau-tau pairs")             currentObject = "ditauPairs";
-<
+      else if(currentObject == "tau-track pairs")           currentObject = "tauTrackPairs";
-<
+      else if(currentObject == "muon-secondary muon pairs") currentObject = "muonSecondaryMuonPairs";
-<
+      else if(currentObject == "secondary muons")           currentObject = "secondaryMuons";
-<
+      else if(currentObject == "electron-secondary electron pairs") currentObject = "electronSecondaryElectronPairs";
-<
+      else if(currentObject == "secondary electrons")       currentObject = "secondaryElectrons";
-<
+      else if(currentObject == "electron-trigobj pairs")    currentObject = "electronTrigobjPairs";
-<
+      else if(currentObject == "muon-trigobj pairs")        currentObject = "muonTrigobjPairs";
-<
-<
+      currentObject.at(0) = toupper(currentObject.at(0));
-<
+      string histoName = "num" + currentObject;
-<
-<
+      if(histoName == "numPrimaryvertexs"){
-<
+        string newHistoName = histoName + "BeforePileupCorrection";
-<
+        oneDHists_.at(currentChannel).at(currentDir)[newHistoName] = directories.at(currentDir).make<TH1D> (TString(newHistoName),channelLabel+" channel: Number of Selected "+currentObject+" Before Pileup Correction; # "+currentObject, maxNum, 0, maxNum);
-<
+        newHistoName = histoName + "AfterPileupCorrection";
-<
+        oneDHists_.at(currentChannel).at(currentDir)[newHistoName] = directories.at(currentDir).make<TH1D> (TString(newHistoName),channelLabel+" channel: Number of Selected "+currentObject+ " After Pileup Correction; # "+currentObject, maxNum, 0, maxNum);
->
->
+    //book all histograms included in the configuration
->
+    for(uint currentDir = 0; currentDir != directories.size(); currentDir++){//loop over all the directories.
->
->
+      TTree* newTree = directories.at(currentDir).make<TTree> (TString("BNTree_"+channelLabel), TString("BNTree_"+channelLabel));
->
+      BNTrees_.push_back(newTree);
->
->
+      for (uint iBranch = 0; iBranch < treeBranches_.size(); iBranch++){
->
+        BranchSpecs currentVar = treeBranches_.at(iBranch);
->
+        vector<float> newVec;
->
+        BNTreeBranchVals_[currentVar.name] = newVec;
->
+        BNTrees_.back()->Branch(TString(currentVar.name), &BNTreeBranchVals_.at(currentVar.name));
->
->
+      } // end for (uint iBranch = 0; iBranch < treeBranches_.size(); iBranch++){
->
->
->
+      for(uint currentHistogramIndex = 0; currentHistogramIndex != histograms.size(); currentHistogramIndex++){
->
->
+        histogram currentHistogram = histograms.at(currentHistogramIndex);
->
+        int numBinsElements = currentHistogram.bins.size();
->
+        int numInputVariables = currentHistogram.inputVariables.size();
->
+        int numBinEdgesX = currentHistogram.variableBinsX.size();
->
+        int numBinEdgesY = currentHistogram.variableBinsY.size();
->
->
+        if(numBinsElements == 1){
->
+          if(numBinEdgesX > 1){
->
+            if(numBinEdgesY > 1)
->
+              numBinsElements = 6;
->
+            else
->
+              numBinsElements = 3;
->
+          }
->
+        }
->
+        if(numBinsElements != 3 && numBinsElements !=6){
->
+          cout << "Error: Didn't find correct number of bin specifications for histogram named '" << currentHistogram.name << "'\n";
->
+          exit(0);
->
+        }
->
+        else if((numBinsElements == 3 && numInputVariables !=1) || (numBinsElements == 6 && numInputVariables!=2)){
->
+          cout << "Error: Didn't find correct number of input variables for histogram named '" << currentHistogram.name << "'\n";
->
+          exit(0);
->
+        }
->
+        else if(numBinsElements == 3){
->
+          if (currentHistogram.bins.size () == 3)
->
+            oneDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH1D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, currentHistogram.bins.at(0), currentHistogram.bins.at(1), currentHistogram.bins.at(2));
->
+          else
->
+            {
->
+              oneDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH1D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, numBinEdgesX - 1, currentHistogram.variableBinsX.data ());
->
+            }
->
+        }
->
+        else if(numBinsElements == 6){
->
+          if (currentHistogram.bins.size () == 6)
->
+            twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH2D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, currentHistogram.bins.at(0), currentHistogram.bins.at(1), currentHistogram.bins.at(2),currentHistogram.bins.at(3),currentHistogram.bins.at(4),currentHistogram.bins.at(5));
->
+          else
->
+            twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name] = directories.at(currentDir).make<TH2D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, numBinEdgesX - 1, currentHistogram.variableBinsX.data (), numBinEdgesY - 1, currentHistogram.variableBinsY.data ());
->
+        }
->
->
->
+        if(currentHistogram.name.find("GenMatch")==string::npos) continue;
->
->
+        // bin      particle type
->
+        // ---      -------------
->
+        //  0        unmatched
->
+        //  1        u
->
+        //  2        d
->
+        //  3        s
->
+        //  4        c
->
+        //  5        b
->
+        //  6        t
->
+        //  7        e
->
+        //  8        mu
->
+        //  9        tau
->
+        // 10        nu
->
+        // 11        g
->
+        // 12        gamma
->
+        // 13        Z
->
+        // 14        W
->
+        // 15        light meson
->
+        // 16        K meson
->
+        // 17        D meson
->
+        // 18        B meson
->
+        // 19        light baryon
->
+        // 20        strange baryon
->
+        // 21        charm baryon
->
+        // 22        bottom baryon
->
+        // 23        QCD string
->
+        // 24        other
->
->
+        vector<TString> labelArray;
->
+        labelArray.push_back("unmatched");
->
+        labelArray.push_back("u");
->
+        labelArray.push_back("d");
->
+        labelArray.push_back("s");
->
+        labelArray.push_back("c");
->
+        labelArray.push_back("b");
->
+        labelArray.push_back("t");
->
+        labelArray.push_back("e");
->
+        labelArray.push_back("#mu");
->
+        labelArray.push_back("#tau");
->
+        labelArray.push_back("#nu");
->
+        labelArray.push_back("g");
->
+        labelArray.push_back("#gamma");
->
+        labelArray.push_back("Z");
->
+        labelArray.push_back("W");
->
+        labelArray.push_back("light meson");
->
+        labelArray.push_back("K meson");
->
+        labelArray.push_back("D meson");
->
+        labelArray.push_back("B meson");
->
+        labelArray.push_back("light baryon");
->
+        labelArray.push_back("strange baryon");
->
+        labelArray.push_back("charm baryon");
->
+        labelArray.push_back("bottom baryon");
->
+        labelArray.push_back("QCD string");
->
+        labelArray.push_back("other");
->
->
+        for(int bin = 0; bin !=currentHistogram.bins.at(0); bin++){
->
+          if(currentHistogram.name.find("GenMatchIdVsMotherId")==string::npos && currentHistogram.name.find("GenMatchIdVsGrandmotherId")==string::npos) {
->
+            oneDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->SetBinLabel(bin+1,labelArray.at(bin));
->
+          }
->
+          else {
->
+            twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetYaxis()->SetBinLabel(bin+1,labelArray.at(currentHistogram.bins.at(0)-bin-1));
->
+            twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->SetBinLabel(bin+1,labelArray.at(bin));
->
+          }
->
+        }
->
+        if(currentHistogram.name.find("GenMatchIdVsMotherId")!=string::npos || currentHistogram.name.find("GenMatchIdVsGrandmotherId")!=string::npos) {
->
+          twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->CenterTitle();
->
+          twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetYaxis()->CenterTitle();
->
+        }
->
->
+      }  // end      for(uint currentHistogramIndex = 0; currentHistogramIndex != histograms.size(); currentHistogramIndex++){
->
->
->
+      // Book a histogram for the number of each object type to be plotted.
->
+      // Name of objectToPlot here must match the name specified in OSUAnalysis::analyze().
->
+      for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++){
->
+        string currentObject = objectsToPlot.at(currentObjectIndex);
->
+        int maxNum = 10;
->
+        if(currentObject == "mcparticles") maxNum = 50;
->
+        else if(currentObject == "primaryvertexs") maxNum = 50;
->
->
+        if(currentObject == "muon-muon pairs")                currentObject = "dimuonPairs";
->
+        else if(currentObject == "electron-electron pairs")   currentObject = "dielectronPairs";
->
+        else if(currentObject == "electron-muon pairs")       currentObject = "electronMuonPairs";
->
+        else if(currentObject == "secondary jets")            currentObject = "secondaryJets";
->
+        else if(currentObject == "jet-jet pairs")             currentObject = "dijetPairs";
->
+        else if(currentObject == "electron-jet pairs")        currentObject = "electronJetPairs";
->
+        else if(currentObject == "muon-jet pairs")            currentObject = "muonJetPairs";
->
+        else if(currentObject == "track-event pairs")         currentObject = "trackEventPairs";
->
+        else if(currentObject == "electron-track pairs")      currentObject = "electronTrackPairs";
->
+        else if(currentObject == "muon-track pairs")          currentObject = "muonTrackPairs";
->
+        else if(currentObject == "muon-tau pairs")            currentObject = "muonTauPairs";
->
+        else if(currentObject == "tau-tau pairs")             currentObject = "ditauPairs";
->
+        else if(currentObject == "tau-track pairs")           currentObject = "tauTrackPairs";
->
+        else if(currentObject == "muon-secondary muon pairs") currentObject = "muonSecondaryMuonPairs";
->
+        else if(currentObject == "secondary muons")           currentObject = "secondaryMuons";
->
+        else if(currentObject == "electron-secondary electron pairs") currentObject = "electronSecondaryElectronPairs";
->
+        else if(currentObject == "secondary electrons")       currentObject = "secondaryElectrons";
->
+        else if(currentObject == "electron-trigobj pairs")    currentObject = "electronTrigobjPairs";
->
+        else if(currentObject == "muon-trigobj pairs")        currentObject = "muonTrigobjPairs";
->
->
+        currentObject.at(0) = toupper(currentObject.at(0));
->
+        string histoName = "num" + currentObject;
->
->
+        if(histoName == "numPrimaryvertexs"){
->
+          string newHistoName = histoName + "BeforePileupCorrection";
->
+          oneDHists_.at(currentChannel).at(currentDir)[newHistoName] = directories.at(currentDir).make<TH1D> (TString(newHistoName),channelLabel+" channel: Number of Selected "+currentObject+" Before Pileup Correction; # "+currentObject, maxNum, 0, maxNum);
->
+          newHistoName = histoName + "AfterPileupCorrection";
->
+          oneDHists_.at(currentChannel).at(currentDir)[newHistoName] = directories.at(currentDir).make<TH1D> (TString(newHistoName),channelLabel+" channel: Number of Selected "+currentObject+ " After Pileup Correction; # "+currentObject, maxNum, 0, maxNum);
->
+        }
->
+        else
->
+          oneDHists_.at(currentChannel).at(currentDir)[histoName] = directories.at(currentDir).make<TH1D> (TString(histoName),channelLabel+" channel: Number of Selected "+currentObject+"; # "+currentObject, maxNum, 0, maxNum);
-<
+      else
-<
+        oneDHists_.at(currentChannel).at(currentDir)[histoName] = directories.at(currentDir).make<TH1D> (TString(histoName),channelLabel+" channel: Number of Selected "+currentObject+"; # "+currentObject, maxNum, 0, maxNum);
-<
+    }
-<
+   }//end of loop over directories
->
+    }//end of loop over directories
+    channels.push_back(tempChannel);
+    tempChannel.cuts.clear();
+  }//end loop over channels
+  //get pile-up event weight
+  if (doPileupReweighting_ && datasetType_ != "data") {
-<
+    if (events->at(0).numTruePV < 0) cout << "WARNING[OSUAnalysis::analyze]: Event has numTruePV<0.  This will give an incorrect pile-up weight." << endl;
-<
+else
-<
+    masterScaleFactor *= puWeight_->at (events->at (0).numTruePV);
->
+    //for "data" datasets, the numTruePV is always set to -1
->
+    if (events->at(0).numTruePV < 0) cout << "WARNING[OSUAnalysis::analyze]: Event has numTruePV<0.  Turning off pile-up reweighting." << endl;
->
+    else masterScaleFactor *= puWeight_->at (events->at (0).numTruePV);
+  stopCTauScaleFactor_ = 1.0;
+    counterMap passingCounter;
+    cumulativeFlags.clear ();
-+
+    for (map<string, vector<float>>::iterator iter = BNTreeBranchVals_.begin();
-+
+         iter != BNTreeBranchVals_.end(); iter++) {
-+
+      iter->second.clear();  // clear array
-+
+    }
-+
-+
+    bool triggerDecision = true;
-<
+    if(currentChannel.triggers.size() != 0){  //triggers specified
-<
+      triggerDecision = evaluateTriggers(currentChannel.triggers,triggers.product());
->
+    if(currentChannel.triggers.size() != 0 || currentChannel.triggersToVeto.size() != 0){  //triggers specified
->
+      triggerDecision = evaluateTriggers(currentChannel.triggers, currentChannel.triggersToVeto, triggers.product());
+      cutFlows_.at(currentChannelIndex)->at ("trigger") = triggerDecision;
+                                                                   "muon-muon pairs");
+        else if(currentObject == "muon-secondary muon pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,muons.product(),muons.product(), \
-<
+                                                                   cumulativeFlags.at("muons").at(flagsForPairCutsIndex), \
-<
+                                                                   cumulativeFlags.at("secondary muons").at(flagsForPairCutsIndex), \
-<
+                                                                   "muon-secondary muon pairs");
->
+                                                                             cumulativeFlags.at("muons").at(flagsForPairCutsIndex), \
->
+                                                                             cumulativeFlags.at("secondary muons").at(flagsForPairCutsIndex), \
->
+                                                                             "muon-secondary muon pairs");
+        else if(currentObject == "electron-secondary electron pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,electrons.product(),electrons.product(), \
-<
+                                                                   cumulativeFlags.at("electrons").at(flagsForPairCutsIndex), \
-<
+                                                                   cumulativeFlags.at("secondary electrons").at(flagsForPairCutsIndex), \
-<
+                                                                   "electron-secondary electron pairs");
->
+                                                                                     cumulativeFlags.at("electrons").at(flagsForPairCutsIndex), \
->
+                                                                                     cumulativeFlags.at("secondary electrons").at(flagsForPairCutsIndex), \
->
+                                                                                     "electron-secondary electron pairs");
+        else if(currentObject == "electron-electron pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,electrons.product(),electrons.product(), \
+                                                                           cumulativeFlags.at("electrons").at(flagsForPairCutsIndex), \
+                                                                       cumulativeFlags.at("muons").at(flagsForPairCutsIndex), \
+                                                                       "electron-muon pairs");
+        else if(currentObject == "jet-jet pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,jets.product(),jets.product(), \
-<
+                                                                           cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
-<
+                                                                           cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
-<
+                                                                           "jet-jet pairs");
->
+                                                                 cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
->
+                                                                 cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
->
+                                                                 "jet-jet pairs");
+        else if(currentObject == "electron-jet pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,electrons.product(),jets.product(), \
-<
+                                                                       cumulativeFlags.at("electrons").at(flagsForPairCutsIndex), \
-<
+                                                                       cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
-<
+                                                                       "electron-jet pairs");
->
+                                                                      cumulativeFlags.at("electrons").at(flagsForPairCutsIndex), \
->
+                                                                      cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
->
+                                                                      "electron-jet pairs");
+        else if(currentObject == "muon-jet pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,muons.product(),jets.product(), \
-<
+                                                                       cumulativeFlags.at("muons").at(flagsForPairCutsIndex), \
-<
+                                                                       cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
-<
+                                                                       "muon-jet pairs");
->
+                                                                  cumulativeFlags.at("muons").at(flagsForPairCutsIndex), \
->
+                                                                  cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
->
+                                                                  "muon-jet pairs");
+        else if(currentObject == "track-event pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,tracks.product(),events.product(),
+                                                                     cumulativeFlags.at("tracks").at(flagsForPairCutsIndex),
+                                                                     cumulativeFlags.at("events").at(flagsForPairCutsIndex),
+                                                                 cumulativeFlags.at("taus").at(flagsForPairCutsIndex),
+                                                                 "tau-tau pairs");
+        else if(currentObject == "tau-track pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,taus .product(),tracks.product(),
-<
+                                                                 cumulativeFlags.at("taus").at(flagsForPairCutsIndex),
-<
+                                                                 cumulativeFlags.at("tracks").at(flagsForPairCutsIndex),
-<
+                                                                 "tau-track pairs");
->
+                                                                   cumulativeFlags.at("taus").at(flagsForPairCutsIndex),
->
+                                                                   cumulativeFlags.at("tracks").at(flagsForPairCutsIndex),
->
+                                                                   "tau-track pairs");
+        else if(currentObject == "electron-trigobj pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,electrons.product(),trigobjs.product(),
+                                                                          cumulativeFlags.at("electrons").at(flagsForPairCutsIndex),
+                                                                          cumulativeFlags.at("trigobjs").at(flagsForPairCutsIndex),
+    //apply trigger (true if none were specified)
+    eventPassedAllCuts = eventPassedAllCuts && triggerDecision;
-–
+    for(uint currentCutIndex = 0; currentCutIndex != currentChannel.cuts.size(); currentCutIndex++){
+      //loop over all objects and count how many passed the cumulative selection up to this point
+    scaleFactor *= electronScaleFactor_;
+    cutFlows_.at(currentChannelIndex)->fillCutFlow(scaleFactor);
-–
+    if(!(GetPlotsAfterEachCut_ || eventPassedAllCuts)) continue;
+    if (printEventInfo_) {
+      // Write information about event to screen, for testing purposes.
+           << "  event=" << events->at(0).evt
+           << endl;
-+
-+
+    //filling histograms
+    for(uint currentCut = 0; currentCut != oneDHists_.at(currentChannelIndex).size(); currentCut++){//loop over all the directories in each channel.
-<
+        uint currentDir;
-<
+        if(!GetPlotsAfterEachCut_){ currentDir =  currentChannel.cuts.size() - oneDHists_.at(currentChannelIndex).size();}//if GetPlotsAfterEachCut_ is true, set currentDir point to the lat cut.
-<
+        else{
-<
+          currentDir = currentCut;
-<
+        }
-<
+        if(eventPassedPreviousCuts.at(currentDir)){
->
+      uint currentDir;
->
+      if(!GetPlotsAfterEachCut_){ currentDir =  currentChannel.cuts.size() - oneDHists_.at(currentChannelIndex).size();}//if GetPlotsAfterEachCut_ is true, set currentDir point to the lat cut.
->
+      else{
->
+        currentDir = currentCut;
->
+      }
->
+      if(eventPassedPreviousCuts.at(currentDir)){
->
->
+        // Assign BNTree variables
->
+        for (uint iBranch = 0; iBranch != treeBranches_.size(); iBranch++) {
->
+          BranchSpecs brSpecs = treeBranches_.at(iBranch);
->
+          string coll = brSpecs.inputCollection;
->
+          if (cumulativeFlags.count(coll) == 0) cout << "Error: no flags found for collection:  " << coll << ", will cause a seg fault" << endl;
->
->
+          if     (coll == "jets")                assignTreeBranch(brSpecs,jets.product(),          cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "secondary jets")      assignTreeBranch(brSpecs,jets.product(),          cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "muons")               assignTreeBranch(brSpecs,muons.product(),         cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "secondary muons")     assignTreeBranch(brSpecs,muons.product(),         cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "electrons")           assignTreeBranch(brSpecs,electrons.product(),     cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "secondary electrons") assignTreeBranch(brSpecs,electrons.product(),     cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "events")              assignTreeBranch(brSpecs,events.product(),        cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "taus")                assignTreeBranch(brSpecs,taus.product(),          cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "mets")                assignTreeBranch(brSpecs,mets.product(),          cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "tracks")              assignTreeBranch(brSpecs,tracks.product(),        cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "genjets")             assignTreeBranch(brSpecs,genjets.product(),       cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "mcparticles")         assignTreeBranch(brSpecs,mcparticles.product(),   cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "primaryvertexs")      assignTreeBranch(brSpecs,primaryvertexs.product(),cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "bxlumis")             assignTreeBranch(brSpecs,bxlumis.product(),       cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "photons")             assignTreeBranch(brSpecs,photons.product(),       cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "superclusters")       assignTreeBranch(brSpecs,superclusters.product(), cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "trigobjs")            assignTreeBranch(brSpecs,trigobjs.product(),      cumulativeFlags.at(coll).at(currentDir));
->
+          else if(coll == "stops"
->
+                  && datasetType_ == "signalMC") assignTreeBranch(brSpecs,stops.product(),         cumulativeFlags.at(coll).at(currentDir));
->
+        } // end loop over histograms
->
->
+        for (uint histogramIndex = 0; histogramIndex != histograms.size(); histogramIndex++){
+          histogram currentHistogram = histograms.at(histogramIndex);
-<
+      if(currentHistogram.inputVariables.size() == 1){
-<
+        TH1D* histo;
-<
+        histo = oneDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
-<
+        if     (currentHistogram.inputCollection == "jets")            fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "secondary jets")  fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muons")           fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "secondary muons") fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "secondary electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-<
+                                                                                       cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
-<
+                                                                                       cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-<
+                                                                                       cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
-<
+                                                                                       cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(),\
-<
+                                                                                               cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir),\
-<
+                                                                                               cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "jet-jet pairs") fill1DHistogram(histo,currentHistogram,jets.product(),jets.product(),\
-<
+                                                                                               cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir),\
-<
+                                                                                               cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-<
+                                                                                       cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
-<
+                                                                                       cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-muon pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),muons.product(), \
-<
+                                                                                           cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir),
-<
+                                                                                           cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-jet pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),jets.product(), \
-<
+                                                                                           cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
-<
+                                                                                           cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-jet pairs") fill1DHistogram(histo,currentHistogram, muons.product(),jets.product(), \
-<
+                                                                                           cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
-<
+                                                                                           cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-track pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),tracks.product(),
-<
+                                                                                            cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-<
+                                                                                            cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-track pairs") fill1DHistogram(histo,currentHistogram, muons.product(),tracks.product(),
-<
+                                                                                        cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-<
+                                                                                        cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-tau pairs") fill1DHistogram(histo,currentHistogram, muons.product(),taus.product(),
-<
+                                                                                      cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-<
+                                                                                      cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "tau-tau pairs") fill1DHistogram(histo,currentHistogram, taus.product(),taus.product(),
-<
+                                                                                     cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-<
+                                                                                     cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "tau-track pairs") fill1DHistogram(histo,currentHistogram, taus.product(),tracks.product(),
-<
+                                                                                     cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-<
+                                                                                     cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),trigobjs.product(),
-<
+                                                                                              cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-<
+                                                                                              cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill1DHistogram(histo,currentHistogram, muons.product(),trigobjs.product(),
-<
+                                                                                          cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-<
+                                                                                          cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
-<
-<
+        else if(currentHistogram.inputCollection == "events") fill1DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "taus") fill1DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "mets") fill1DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "tracks") fill1DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "genjets") fill1DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "mcparticles") fill1DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "primaryvertexs") fill1DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "bxlumis") fill1DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "photons") fill1DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "superclusters") fill1DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "trigobjs") fill1DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill1DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
-<
+      }
-<
+      else if(currentHistogram.inputVariables.size() == 2){
-<
+        TH2D* histo;
-<
+        histo = twoDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
-<
+        if     (currentHistogram.inputCollection == "jets")            fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "secondary jets")  fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muons")           fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "secondary muons") fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-<
+                                                                                       cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
-<
+                                                                                       cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-<
+                                                                                       cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
-<
+                                                                                       cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "secondary electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-<
+                                                                                               cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir), \
-<
+                                                                                               cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "jet-jet pairs") fill2DHistogram(histo,currentHistogram,jets.product(),jets.product(), \
-<
+                                                                                               cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
-<
+                                                                                               cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-<
+                                                                                       cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
-<
+                                                                                       cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-muon pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),muons.product(), \
-<
+                                                                                           cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
-<
+                                                                                           cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-jet pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),jets.product(), \
-<
+                                                                                           cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
-<
+                                                                                           cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-jet pairs") fill2DHistogram(histo,currentHistogram,muons.product(),jets.product(), \
-<
+                                                                                           cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
-<
+                                                                                           cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-track pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),tracks.product(),
-<
+                                                                                            cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-<
+                                                                                            cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-track pairs") fill2DHistogram(histo,currentHistogram,muons.product(),tracks.product(),
-<
+                                                                                        cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-<
+                                                                                        cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-tau pairs") fill2DHistogram(histo,currentHistogram,muons.product(),taus.product(),
-<
+                                                                                      cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-<
+                                                                                      cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "tau-tau pairs") fill2DHistogram(histo,currentHistogram,taus.product(),taus.product(),
-<
+                                                                                     cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-<
+                                                                                     cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "tau-track pairs") fill2DHistogram(histo,currentHistogram,taus.product(),tracks.product(),
-<
+                                                                                     cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-<
+                                                                                     cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),trigobjs.product(),
-<
+                                                                                              cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-<
+                                                                                              cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill2DHistogram(histo,currentHistogram,muons.product(),trigobjs.product(),
-<
+                                                                                          cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-<
+                                                                                          cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "events") fill2DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "taus") fill2DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "mets") fill2DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "tracks") fill2DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "track-event pairs") fill2DHistogram(histo,currentHistogram,tracks.product(),events.product(),
-<
+                                                                                         cumulativeFlags.at("tracks").at(currentDir),cumulativeFlags.at("events").at(currentDir),
-<
+                                                                                         cumulativeFlags.at("track-event pairs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "genjets") fill2DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "mcparticles") fill2DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "primaryvertexs") fill2DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "bxlumis") fill2DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "photons") fill2DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "superclusters") fill2DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "trigobjs") fill2DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
-<
+        else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill2DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
-<
+      }
-<
+    }
-<
->
+          if (cumulativeFlags.count(currentHistogram.inputCollection) == 0) cout << "Error: no flags found for collection:  " << currentHistogram.inputCollection << ", will cause a seg fault" << endl;
-+
+          if(currentHistogram.inputVariables.size() == 1){
-+
+            TH1D* histo;
-+
+            histo = oneDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
-+
+            if     (currentHistogram.inputCollection == "jets")            fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "secondary jets")  fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muons")           fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "secondary muons") fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "secondary electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-+
+                                                                                           cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
-+
+                                                                                           cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-+
+                                                                                                     cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
-+
+                                                                                                     cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(),\
-+
+                                                                                                   cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir),\
-+
+                                                                                                   cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "jet-jet pairs") fill1DHistogram(histo,currentHistogram,jets.product(),jets.product(),\
-+
+                                                                                         cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir),\
-+
+                                                                                         cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-+
+                                                                                                             cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
-+
+                                                                                                             cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-muon pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),muons.product(), \
-+
+                                                                                               cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir),
-+
+                                                                                               cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-jet pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),jets.product(), \
-+
+                                                                                              cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
-+
+                                                                                              cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-jet pairs") fill1DHistogram(histo,currentHistogram, muons.product(),jets.product(), \
-+
+                                                                                          cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
-+
+                                                                                          cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-track pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),tracks.product(),
-+
+                                                                                                cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-+
+                                                                                                cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-track pairs") fill1DHistogram(histo,currentHistogram, muons.product(),tracks.product(),
-+
+                                                                                            cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-+
+                                                                                            cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-tau pairs") fill1DHistogram(histo,currentHistogram, muons.product(),taus.product(),
-+
+                                                                                          cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-+
+                                                                                          cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "tau-tau pairs") fill1DHistogram(histo,currentHistogram, taus.product(),taus.product(),
-+
+                                                                                         cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-+
+                                                                                         cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "tau-track pairs") fill1DHistogram(histo,currentHistogram, taus.product(),tracks.product(),
-+
+                                                                                           cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-+
+                                                                                           cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),trigobjs.product(),
-+
+                                                                                                  cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-+
+                                                                                                  cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill1DHistogram(histo,currentHistogram, muons.product(),trigobjs.product(),
-+
+                                                                                              cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-+
+                                                                                              cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
-+
-+
+            else if(currentHistogram.inputCollection == "events") fill1DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "taus") fill1DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "mets") fill1DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "tracks") fill1DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "genjets") fill1DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "mcparticles") fill1DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "primaryvertexs") fill1DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "bxlumis") fill1DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "photons") fill1DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "superclusters") fill1DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "trigobjs") fill1DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill1DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
-+
+          }
-+
+          else if(currentHistogram.inputVariables.size() == 2){
-+
+            TH2D* histo;
-+
+            histo = twoDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
-+
+            if     (currentHistogram.inputCollection == "jets")            fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "secondary jets")  fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muons")           fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "secondary muons") fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-+
+                                                                                           cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
-+
+                                                                                           cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-+
+                                                                                                     cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
-+
+                                                                                                     cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "secondary electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-+
+                                                                                                   cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir), \
-+
+                                                                                                   cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "jet-jet pairs") fill2DHistogram(histo,currentHistogram,jets.product(),jets.product(), \
-+
+                                                                                         cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
-+
+                                                                                         cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-+
+                                                                                                             cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
-+
+                                                                                                             cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-muon pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),muons.product(), \
-+
+                                                                                               cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
-+
+                                                                                               cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-jet pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),jets.product(), \
-+
+                                                                                              cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
-+
+                                                                                              cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-jet pairs") fill2DHistogram(histo,currentHistogram,muons.product(),jets.product(), \
-+
+                                                                                          cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
-+
+                                                                                          cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-track pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),tracks.product(),
-+
+                                                                                                cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-+
+                                                                                                cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-track pairs") fill2DHistogram(histo,currentHistogram,muons.product(),tracks.product(),
-+
+                                                                                            cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-+
+                                                                                            cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-tau pairs") fill2DHistogram(histo,currentHistogram,muons.product(),taus.product(),
-+
+                                                                                          cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-+
+                                                                                          cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "tau-tau pairs") fill2DHistogram(histo,currentHistogram,taus.product(),taus.product(),
-+
+                                                                                         cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
-+
+                                                                                         cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "tau-track pairs") fill2DHistogram(histo,currentHistogram,taus.product(),tracks.product(),
-+
+                                                                                           cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
-+
+                                                                                           cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),trigobjs.product(),
-+
+                                                                                                  cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-+
+                                                                                                  cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill2DHistogram(histo,currentHistogram,muons.product(),trigobjs.product(),
-+
+                                                                                              cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
-+
+                                                                                              cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "events") fill2DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "taus") fill2DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "mets") fill2DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "tracks") fill2DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "track-event pairs") fill2DHistogram(histo,currentHistogram,tracks.product(),events.product(),
-+
+                                                                                             cumulativeFlags.at("tracks").at(currentDir),cumulativeFlags.at("events").at(currentDir),
-+
+                                                                                             cumulativeFlags.at("track-event pairs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "genjets") fill2DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "mcparticles") fill2DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "primaryvertexs") fill2DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "bxlumis") fill2DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "photons") fill2DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "superclusters") fill2DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "trigobjs") fill2DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
-+
+            else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill2DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
-+
+          }
-+
-+
+        }
-+
-+
-+
+        //fills histograms with the sizes of collections
-+
+        for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++){
-+
-+
+          string currentObject = objectsToPlot.at(currentObjectIndex);
-+
+          string objectToPlot = "";
-+
-+
+          // Name of objectToPlot here must match the name specified in OSUAnalysis::OSUAnalysis().
-+
+          if(currentObject == "muon-muon pairs")                         objectToPlot = "dimuonPairs";
-+
+          else if(currentObject == "electron-electron pairs")            objectToPlot = "dielectronPairs";
-+
+          else if(currentObject == "electron-muon pairs")                objectToPlot = "electronMuonPairs";
-+
+          else if(currentObject == "electron-jet pairs")                 objectToPlot = "electronJetPairs";
-+
+          else if(currentObject == "muon-jet pairs")                     objectToPlot = "muonJetPairs";
-+
+          else if(currentObject == "jet-jet pairs")                      objectToPlot = "dijetPairs";
-+
+          else if(currentObject == "secondary jets")                     objectToPlot = "secondaryJets";
-+
+          else if(currentObject == "electron-track pairs")               objectToPlot = "electronTrackPairs";
-+
+          else if(currentObject == "muon-track pairs")                   objectToPlot = "muonTrackPairs";
-+
+          else if(currentObject == "muon-tau pairs")                     objectToPlot = "muonTauPairs";
-+
+          else if(currentObject == "tau-tau pairs")                      objectToPlot = "ditauPairs";
-+
+          else if(currentObject == "tau-track pairs")                    objectToPlot = "tauTrackPairs";
-+
+          else if(currentObject == "track-event pairs")                  objectToPlot = "trackEventPairs";
-+
+          else if(currentObject == "muon-secondary muon pairs")          objectToPlot = "muonSecondaryMuonPairs";
-+
+          else if(currentObject == "secondary muons")                    objectToPlot = "secondaryMuons";
-+
+          else if(currentObject == "electron-secondary electron pairs")  objectToPlot = "electronSecondaryElectronPairs";
-+
+          else if(currentObject == "secondary electrons")                objectToPlot = "secondaryElectrons";
-+
+          else if(currentObject == "electron-trigobj pairs")             objectToPlot = "electronTrigobjPairs";
-+
+          else if(currentObject == "muon-trigobj pairs")                 objectToPlot = "muonTrigobjPairs";
-+
+          else objectToPlot = currentObject;
-+
-+
+          string tempCurrentObject = objectToPlot;
-+
+          tempCurrentObject.at(0) = toupper(tempCurrentObject.at(0));
-+
+          string histoName = "num" + tempCurrentObject;
-+
-+
+          //set position of primary vertex in event, in order to calculate quantities relative to it
-+
+          if(find(objectsToCut.begin(), objectsToCut.end(), currentObject) != objectsToCut.end()) {
-+
+            vector<bool> lastCutFlags = cumulativeFlags.at(currentObject).at(currentDir);
-+
+            int numToPlot = 0;
-+
+            for (uint currentFlag = 0; currentFlag != lastCutFlags.size(); currentFlag++){
-+
+              if(lastCutFlags.at(currentFlag)) numToPlot++;
-+
+            }
-+
+            if(objectToPlot == "primaryvertexs"){
-+
+              oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
-+
+              oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
-+
+            }
-+
+            else {
-+
+              oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(numToPlot,scaleFactor);
-+
+            }
-+
+          }
-+
+          else if(objectToPlot == "jets")           oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(jets->size(),scaleFactor);
-+
+          else if(objectToPlot == "secondaryJets")  oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(jets->size(),scaleFactor);
-+
+          else if(objectToPlot == "muons")          oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size(),scaleFactor);
-+
+          else if(objectToPlot == "secondaryMuons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size(),scaleFactor);
-+
+          else if(objectToPlot == "dimuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
-+
+          else if(objectToPlot == "muonSecondaryMuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
-+
+          else if(objectToPlot == "electrons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size(),scaleFactor);
-+
+          else if(objectToPlot == "secondaryElectrons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size(),scaleFactor);
-+
+          else if(objectToPlot == "dielectronPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
-+
+          else if(objectToPlot == "electronSecondaryElectronPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
-+
+          else if(objectToPlot == "electronMuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*muons->size(),scaleFactor);
-+
+          else if(objectToPlot == "electronJetPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*jets->size(),scaleFactor);
-+
+          else if(objectToPlot == "muonJetPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*jets->size(),scaleFactor);
-+
+          else if(objectToPlot == "electronTrackPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*tracks->size(),scaleFactor);
-+
+          else if(objectToPlot == "electronTrigobjPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*trigobjs->size(),scaleFactor);
-+
+          else if(objectToPlot == "muonTrigobjPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*trigobjs->size(),scaleFactor);
-+
+          else if(objectToPlot == "events") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(events->size(),scaleFactor);
-+
+          else if(objectToPlot == "taus") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(taus->size(),scaleFactor);
-+
+          else if(objectToPlot == "mets") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(mets->size(),scaleFactor);
-+
+          else if(objectToPlot == "tracks") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(tracks->size(),scaleFactor);
-+
+          else if(objectToPlot == "genjets") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(genjets->size(),scaleFactor);
-+
+          else if(objectToPlot == "mcparticles") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(mcparticles->size(),scaleFactor);
-+
+          else if(objectToPlot == "bxlumis") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(bxlumis->size(),scaleFactor);
-+
+          else if(objectToPlot == "photons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(photons->size(),scaleFactor);
-+
+          else if(objectToPlot == "superclusters") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(superclusters->size(),scaleFactor);
-+
+          else if(objectToPlot == "trigobjs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(trigobjs->size(),scaleFactor);
-+
+          else if(objectToPlot == "primaryvertexs"){
-+
+            oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
-+
+            oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
-+
+          }
-+
+          if(objectToPlot == "stops" && datasetType_ == "signalMC") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(stops->size(),scaleFactor);
-<
+    //fills histograms with the sizes of collections
-<
+    for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++){
-<
-<
+      string currentObject = objectsToPlot.at(currentObjectIndex);
-<
+      string objectToPlot = "";
-<
-<
+      // Name of objectToPlot here must match the name specified in OSUAnalysis::OSUAnalysis().
-<
+      if(currentObject == "muon-muon pairs")                         objectToPlot = "dimuonPairs";
-<
+      else if(currentObject == "electron-electron pairs")            objectToPlot = "dielectronPairs";
-<
+      else if(currentObject == "electron-muon pairs")                objectToPlot = "electronMuonPairs";
-<
+      else if(currentObject == "electron-jet pairs")                 objectToPlot = "electronJetPairs";
-<
+      else if(currentObject == "muon-jet pairs")                     objectToPlot = "muonJetPairs";
-<
+      else if(currentObject == "jet-jet pairs")                      objectToPlot = "dijetPairs";
-<
+      else if(currentObject == "secondary jets")                     objectToPlot = "secondaryJets";
-<
+      else if(currentObject == "electron-track pairs")               objectToPlot = "electronTrackPairs";
-<
+      else if(currentObject == "muon-track pairs")                   objectToPlot = "muonTrackPairs";
-<
+      else if(currentObject == "muon-tau pairs")                     objectToPlot = "muonTauPairs";
-<
+      else if(currentObject == "tau-tau pairs")                      objectToPlot = "ditauPairs";
-<
+      else if(currentObject == "tau-track pairs")                    objectToPlot = "tauTrackPairs";
-<
+      else if(currentObject == "track-event pairs")                  objectToPlot = "trackEventPairs";
-<
+      else if(currentObject == "muon-secondary muon pairs")          objectToPlot = "muonSecondaryMuonPairs";
-<
+      else if(currentObject == "secondary muons")                    objectToPlot = "secondaryMuons";
-<
+      else if(currentObject == "electron-secondary electron pairs")  objectToPlot = "electronSecondaryElectronPairs";
-<
+      else if(currentObject == "secondary electrons")                objectToPlot = "secondaryElectrons";
-<
+      else if(currentObject == "electron-trigobj pairs")             objectToPlot = "electronTrigobjPairs";
-<
+      else if(currentObject == "muon-trigobj pairs")                 objectToPlot = "muonTrigobjPairs";
-<
+      else objectToPlot = currentObject;
-<
-<
+      string tempCurrentObject = objectToPlot;
-<
+      tempCurrentObject.at(0) = toupper(tempCurrentObject.at(0));
-<
+      string histoName = "num" + tempCurrentObject;
-<
-<
+      //set position of primary vertex in event, in order to calculate quantities relative to it
-<
+      if(find(objectsToCut.begin(), objectsToCut.end(), currentObject) != objectsToCut.end()) {
-<
+        vector<bool> lastCutFlags = cumulativeFlags.at(currentObject).at(currentDir);
-<
+        int numToPlot = 0;
-<
+        for (uint currentFlag = 0; currentFlag != lastCutFlags.size(); currentFlag++){
-<
+          if(lastCutFlags.at(currentFlag)) numToPlot++;
-<
+        }
-<
+        if(objectToPlot == "primaryvertexs"){
-<
+          oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
-<
+          oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
-<
+        }
-<
+        else {
-<
+          oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(numToPlot,scaleFactor);
-<
+        }
->
+        } // end for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++)
-–
+      else if(objectToPlot == "jets")           oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(jets->size(),scaleFactor);
-–
+      else if(objectToPlot == "secondaryJets")  oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(jets->size(),scaleFactor);
-–
+      else if(objectToPlot == "muons")          oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size(),scaleFactor);
-–
+      else if(objectToPlot == "secondaryMuons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size(),scaleFactor);
-–
+      else if(objectToPlot == "dimuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
-–
+      else if(objectToPlot == "muonSecondaryMuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
-–
+      else if(objectToPlot == "electrons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size(),scaleFactor);
-–
+      else if(objectToPlot == "secondaryElectrons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size(),scaleFactor);
-–
+      else if(objectToPlot == "dielectronPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
-–
+      else if(objectToPlot == "electronSecondaryElectronPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
-–
+      else if(objectToPlot == "electronMuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*muons->size(),scaleFactor);
-–
+      else if(objectToPlot == "electronJetPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*jets->size(),scaleFactor);
-–
+      else if(objectToPlot == "muonJetPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*jets->size(),scaleFactor);
-–
+      else if(objectToPlot == "electronTrackPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*tracks->size(),scaleFactor);
-–
+      else if(objectToPlot == "electronTrigobjPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*trigobjs->size(),scaleFactor);
-–
+      else if(objectToPlot == "muonTrigobjPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*trigobjs->size(),scaleFactor);
-–
+      else if(objectToPlot == "events") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(events->size(),scaleFactor);
-–
+      else if(objectToPlot == "taus") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(taus->size(),scaleFactor);
-–
+      else if(objectToPlot == "mets") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(mets->size(),scaleFactor);
-–
+      else if(objectToPlot == "tracks") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(tracks->size(),scaleFactor);
-–
+      else if(objectToPlot == "genjets") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(genjets->size(),scaleFactor);
-–
+      else if(objectToPlot == "mcparticles") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(mcparticles->size(),scaleFactor);
-–
+      else if(objectToPlot == "bxlumis") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(bxlumis->size(),scaleFactor);
-–
+      else if(objectToPlot == "photons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(photons->size(),scaleFactor);
-–
+      else if(objectToPlot == "superclusters") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(superclusters->size(),scaleFactor);
-–
+      else if(objectToPlot == "trigobjs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(trigobjs->size(),scaleFactor);
-–
+      else if(objectToPlot == "primaryvertexs"){
-–
+        oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
-–
+        oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
-–
+      }
-–
+      if(objectToPlot == "stops" && datasetType_ == "signalMC") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(stops->size(),scaleFactor);
-–
-–
+     } // end for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++)
-<
+  }
->
->
+    BNTrees_.at(currentChannelIndex)->Fill();
->
+  } //end loop over channel
+  masterCutFlow_->fillCutFlow(masterScaleFactor);
-+
+} // end void OSUAnalysis::analyze (const edm::Event &event, const edm::EventSetup &setup)
+bool
-<
+OSUAnalysis::evaluateTriggers (vector<string> triggersToTest, const BNtriggerCollection* triggerCollection){
->
+OSUAnalysis::evaluateTriggers (vector<string> triggersToTest, vector<string> triggersToVeto, const BNtriggerCollection* triggerCollection){
-+
+  //initialize to false until a chosen trigger is passed
+  bool triggerDecision = false;
-<
+  for (BNtriggerCollection::const_iterator trigger = triggerCollection->begin (); trigger != triggerCollection->end (); trigger++)
-<
+    {
-<
+      if(trigger->pass != 1) continue;
-<
+      for(uint triggerName = 0; triggerName != triggersToTest.size(); triggerName++)
-<
+        {
-<
+          if(trigger->name.find(triggersToTest.at(triggerName))!=string::npos){
-<
+            triggerDecision = true;
-<
+          }
-<
+        }
-<
+    }
-<
+  return triggerDecision;
->
+  //loop over all triggers defined in the event
->
+  for (BNtriggerCollection::const_iterator trigger = triggerCollection->begin (); trigger != triggerCollection->end (); trigger++){
->
->
+    //we're only interested in triggers that actually passed
->
+    if(trigger->pass != 1) continue;
-+
+    //if the event passes one of the veto triggers, exit and return false
-+
+    for(uint triggerName = 0; triggerName != triggersToVeto.size(); triggerName++){
-+
+      if(trigger->name.find(triggersToVeto.at(triggerName))!=string::npos) return false;
-+
+    }
-+
+    //if the event passes one of the chosen triggers, set triggerDecision to true
-+
+    for(uint triggerName = 0; triggerName != triggersToTest.size(); triggerName++){
-+
+      if(trigger->name.find(triggersToTest.at(triggerName))!=string::npos) triggerDecision = true;
-+
+    }
-+
+  }
-+
+  //if none of the veto triggers fired:
-+
+  //return the OR of all the chosen triggers
-+
+  if (triggersToTest.size() != 0) return triggerDecision;
-+
+  //or if no triggers were defined return true
-+
+  else return true;
-+
+vector<string>
+OSUAnalysis::splitString (string inputString){
+    value = object->isGlobalMuon > 0                \
+      && object->isPFMuon > 0                        \
+      && object->normalizedChi2 < 10                \
-<
+                                  && object->numberOfValidMuonHits > 0        \
->
+      && object->numberOfValidMuonHits > 0        \
+      && object->numberOfMatchedStations > 1        \
+      && fabs(object->correctedD0Vertex) < 0.2        \
+      && fabs(object->correctedDZ) < 0.5        \
+    value = object->isGlobalMuon > 0                \
+      && object->isPFMuon > 0                        \
+      && object->normalizedChi2 < 10                \
-<
+                                  && object->numberOfValidMuonHits > 0        \
->
+      && object->numberOfValidMuonHits > 0        \
+      && object->numberOfMatchedStations > 1        \
+      && object->numberOfValidPixelHits > 0        \
+      && object->numberOfLayersWithMeasurement > 5;
+    d0Error = hypot (object->tkD0err, hypot (chosenVertex ()->xError, chosenVertex ()->yError));
+    dzError = hypot (object->tkDZerr, chosenVertex ()->zError);
+    value = fabs (object->correctedD0Vertex) > 0.2 || fabs (object->correctedDZ) > 0.5
-<
+         || fabs (object->correctedD0Vertex / d0Error) > 99.0
-<
+         || fabs (object->correctedDZ / dzError) > 99.0;
->
+      || fabs (object->correctedD0Vertex / d0Error) > 99.0
->
+      || fabs (object->correctedDZ / dzError) > 99.0;
+    value = (object->isStandAloneMuon && !object->isTrackerMuon && !object->isGlobalMuon) || !value;
+    d0Error = hypot (object->tkD0err, hypot (chosenVertex ()->xError, chosenVertex ()->yError));
+    dzError = hypot (object->tkDZerr, chosenVertex ()->zError);
+    value = fabs (object->correctedD0Vertex) > 0.2 || fabs (object->correctedDZ) > 0.5
-<
+         || fabs (object->correctedD0Vertex / d0Error) > 99.0
-<
+         || fabs (object->correctedDZ / dzError) > 99.0;
->
+      || fabs (object->correctedD0Vertex / d0Error) > 99.0
->
+      || fabs (object->correctedDZ / dzError) > 99.0;
+    value = !value;
+    for(BNprimaryvertexCollection::const_iterator vertex = primaryvertexs->begin (); vertex != primaryvertexs->end (); vertex++){
+      vertex_rank++;
+      int dist = sqrt((object->vx-vertex->x)*(object->vx-vertex->x) + \
-<
+                   (object->vy-vertex->y)*(object->vy-vertex->y) + \
-<
+                   (object->vz-vertex->z)*(object->vz-vertex->z));
->
+                      (object->vy-vertex->y)*(object->vy-vertex->y) + \
->
+                      (object->vz-vertex->z)*(object->vz-vertex->z));
+      if(abs(dist) < abs(minDistToVertex)){
+        value = vertex_rank;
+    string obj1, obj2;
+    getTwoObjs(currentCut.inputCollection, obj1, obj2);
+    if (inputType==obj1 ||
-<
+          inputType==obj2) {
->
+        inputType==obj2) {
+      // Do not add a cut to individualFlags or cumulativeFlags, if the cut is on a paired collection,
+      // and the inputType is a member of the pair.
+      // The cut will instead be applied when the setObjectFlags() is called for the paired collection.
-+
-+
+template <class InputCollection>
-+
+void OSUAnalysis::assignTreeBranch(BranchSpecs parameters, InputCollection inputCollection, vector<bool> flags){
-+
+  // This function is similar to fill1DHistogram(), but instead of filling a histogram it assigns a value to a variable for the BNTree
-+
-+
+  for (uint object = 0; object != inputCollection->size(); object++){
-+
-+
+    if (!plotAllObjectsInPassingEvents_ && !flags.at(object)) continue;
-+
-+
+    string inputVariable = parameters.inputVariable;
-+
+    string function = "";
-+
+    string stringValue = "";
-+
+    double value = valueLookup(&inputCollection->at(object), inputVariable, function, stringValue);
-+
+    BNTreeBranchVals_.at(parameters.name).push_back(value);
-+
-+
+  }
-+
+}
-+
-+
+template <class InputCollection>
+void OSUAnalysis::fill1DHistogram(TH1* histo, histogram parameters, InputCollection inputCollection,vector<bool> flags, double scaleFactor){
+      function = currentString.substr(0,currentString.find("("));//function comes before the "("
+      inputVariable = currentString.substr(currentString.find("(")+1);//get rest of string
+      inputVariable = inputVariable.substr(0,inputVariable.size()-1);//remove trailing ")"
-<
+    }
->
+    }
->
->
+    string stringValue = "";
->
+    double value = valueLookup(&inputCollection->at(object), inputVariable, function, stringValue);
->
+    histo->Fill(value,scaleFactor);
-–
+    string stringValue = "";
-–
+    double value = valueLookup(&inputCollection->at(object), inputVariable, function, stringValue);
-–
+    histo->Fill(value,scaleFactor);
+    if (printEventInfo_) {
+      // Write information about event to screen, for testing purposes.
+      cout << "  Info for event:  value for histogram " << histo->GetName() << ":  " << value << endl;

Diff Legend

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

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents): Revision 1.70 by qpython, Wed May 29 15:29:27 2013 UTC vs. Revision 1.75 by wulsin, Fri May 31 09:15:04 2013 UTC

Diff Legend

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents):
Revision 1.70 by qpython, Wed May 29 15:29:27 2013 UTC vs.
Revision 1.75 by wulsin, Fri May 31 09:15:04 2013 UTC