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

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents):
Revision 1.67 by wulsin, Tue May 28 21:35:43 2013 UTC vs.
Revision 1.73 by lantonel, Thu May 30 18:09:58 2013 UTC

+  //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");
-+
+    vector<edm::ParameterSet> histogramList_  (histogramSets_.at(currentHistogramSet).getParameter<vector<edm::ParameterSet> >("histograms"));
+    for(uint currentHistogram = 0; currentHistogram != histogramList_.size(); currentHistogram++){
+    //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");
->
+    }
+    vector<string> subSubDirNames;//subdirectories in each channel.
+    //get list of cuts for this channel
+    vector<edm::ParameterSet> cuts_  (channels_.at(currentChannel).getParameter<vector<edm::ParameterSet> >("cuts"));
-<
->
+    //loop over and parse all cuts
+    for(uint currentCut = 0; currentCut != cuts_.size(); currentCut++){
+      int numberRequiredInt = atoi(numberRequiredVector.at(1).c_str());
+      tempCut.numberRequired = numberRequiredInt;// number of objects required to pass the cut
+      tempCut.eventComparativeOperator = numberRequiredVector.at(0);// comparison to make
-<
->
+      //Set up vectors to store the directories and subDIrectories for each channel.
-<
+      string subSubDirName;
->
+      string subSubDirName;
+      string tempCutName;
+      if(cuts_.at(currentCut).exists("alias")){
+        tempCutName = cuts_.at(currentCut).getParameter<string> ("alias");
+        tempCutName = cutName;
+        subSubDirName = "After " + numberRequiredString + " " + collectionString + " with " + cutString + " Cut Applied";
-<
->
+      subSubDirNames.push_back(subSubDirName);
+      tempCut.name = tempCutName;
+      tempChannel.cuts.push_back(tempCut);
-<
->
+    }//end loop over cuts
-<
->
+    vector<map<string, TH1D*>> oneDHistsTmp;
-<
+    vector<map<string, TH2D*>> twoDHistsTmp;
->
+    vector<map<string, TH2D*>> twoDHistsTmp;
+    //book a directory in the output file with the name of the channel
+    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++){
->
+       for( uint currentDir=0; currentDir != subSubDirNames.size(); currentDir++){
+            TFileDirectory subSubDir = subDir.mkdir( subSubDirNames[currentDir] );
+            directories.push_back(subSubDir);
+            map<string, TH1D*> oneDhistoMap;
+            map<string, TH2D*> twoDhistoMap;
+            twoDHistsTmp.push_back(twoDhistoMap);
-<
+      oneDHists_.push_back(oneDHistsTmp);
->
+      oneDHists_.push_back(oneDHistsTmp);
+      twoDHists_.push_back(twoDHistsTmp);
-<
+   //only set up directories with names of the channels if GetPlotsAfterEachCut_ is false.
->
+   //only set up directories with names of the channels if GetPlotsAfterEachCut_ is false.
+   else{
+      map<string, TH1D*> oneDhistoMap;
+      oneDHistsTmp.push_back(oneDhistoMap);
+  double masterScaleFactor = 1.0;
+  //get pile-up event weight
-<
+  if(doPileupReweighting_ && datasetType_ != "data") masterScaleFactor *= puWeight_->at (events->at (0).numTruePV);
->
+  if (doPileupReweighting_ && datasetType_ != "data") {
->
+    //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;
+  if (datasetType_ == "signalMC" && regex_match (dataset_, regex ("stop.*to.*_.*mm.*")))
+    cumulativeFlags.clear ();
+    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;
+    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)){
-–
+        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);
-–
+      }
-–
+    }
-<
+    //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);
->
+    //filling histograms
->
+    for(uint currentCut = 0; currentCut != oneDHists_.at(currentChannelIndex).size(); currentCut++){//loop over all the directories in each channel.
-<
+     } // end for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++)
->
+      uint currentDir;
->
+      if(!GetPlotsAfterEachCut_){ currentDir =  currentChannel.cuts.size() - oneDHists_.at(currentChannelIndex).size();}//if GetPlotsAfterEachCut_ is true, set currentDir point to the last cut.
->
+      else currentDir = currentCut;
->
->
+      if(!eventPassedPreviousCuts.at(currentDir)) continue;
->
->
->
+      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);
->
+        }
->
+      }
->
->
+      //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;
->
->
->
+        vector<bool> lastCutFlags = cumulativeFlags.at(currentObject).at(currentDir);
->
->
+        //count the number of objects passing all cuts
->
+        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++)
-<
+  }
->
+  } //end loop over channel
+  masterCutFlow_->fillCutFlow(masterScaleFactor);
+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){
+  else if(variable == "id2") value = object->id2;
+  else if(variable == "evt") value = object->evt;
+  else if(variable == "puScaleFactor"){
-<
+    if(datasetType_ != "data")
->
+    if(doPileupReweighting_ && datasetType_ != "data")
+      value = puWeight_->at (events->at (0).numTruePV);
+    else
+      value = 1.0;

Diff Legend

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

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents): Revision 1.67 by wulsin, Tue May 28 21:35:43 2013 UTC vs. Revision 1.73 by lantonel, Thu May 30 18:09:58 2013 UTC

Diff Legend

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents):
Revision 1.67 by wulsin, Tue May 28 21:35:43 2013 UTC vs.
Revision 1.73 by lantonel, Thu May 30 18:09:58 2013 UTC