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

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents):
Revision 1.64 by ahart, Sun May 19 00:04:50 2013 UTC vs.
Revision 1.73 by lantonel, Thu May 30 18:09:58 2013 UTC

+  applyLeptonSF_ (cfg.getParameter<bool> ("applyLeptonSF")),
+  printEventInfo_ (cfg.getParameter<bool> ("printEventInfo")),
+  useTrackCaloRhoCorr_ (cfg.getParameter<bool> ("useTrackCaloRhoCorr")),
-<
+  stopCTau_ (cfg.getParameter<vector<double> > ("stopCTau"))
-<
->
+  stopCTau_ (cfg.getParameter<vector<double> > ("stopCTau")),
->
+  GetPlotsAfterEachCut_ (cfg.getParameter<bool> ("GetPlotsAfterEachCut"))
+  TH1::SetDefaultSumw2 ();
+  // Construct Cutflow Objects. These store the results of cut decisions and
+  // handle filling cut flow histograms.
+  masterCutFlow_ = new CutFlow (fs_);
-–
+  vector<TFileDirectory> directories;
+  //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");
->
+    }
+    //create cutFlow for this channel
+    cutFlows_.push_back (new CutFlow (fs_, channelName));
-+
+    vector<TFileDirectory> directories; //vector of directories in the output file.
-+
+    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"));
-–
+    //book a directory in the output file with the name of the channel
-–
+    TFileDirectory subDir = fs_->mkdir( channelName );
-–
+    directories.push_back(subDir);
-<
+    map<string, TH1D*> oneDhistoMap;
-<
+    oneDHists_.push_back(oneDhistoMap);
-<
+    map<string, TH2D*> twoDhistoMap;
-<
+    twoDHists_.push_back(twoDhistoMap);
->
+    //loop over and parse all cuts
->
+    for(uint currentCut = 0; currentCut != cuts_.size(); currentCut++){
->
+      cut tempCut;
->
+      //store input collection for cut
->
+      string tempInputCollection = cuts_.at(currentCut).getParameter<string> ("inputCollection");
->
+      if(tempInputCollection == "muon-electron pairs") tempInputCollection = "electron-muon pairs";
->
+      if(tempInputCollection == "jet-electron pairs") tempInputCollection = "electron-jet pairs";
->
+      if(tempInputCollection == "jet-muon pairs") tempInputCollection = "muon-jet pairs";
->
+      if(tempInputCollection == "event-track pairs")   tempInputCollection = "track-event pairs";
->
+      if(tempInputCollection == "secondary muon-muon pairs")   tempInputCollection = "muon-secondary muon pairs";
->
+      if(tempInputCollection == "secondary electron-electron pairs")   tempInputCollection = "electron-secondary electron pairs";
->
+      if(tempInputCollection == "trigobj-electron pairs")   tempInputCollection = "electron-trigobj pairs";
->
+      if(tempInputCollection == "trigobj-muon pairs")   tempInputCollection = "muon-trigobj pairs";
->
+      tempCut.inputCollection = tempInputCollection;
->
+      if(tempInputCollection.find("pairs")==string::npos){ //just a single object
->
+        if(tempInputCollection.find("secondary")!=string::npos){//secondary object
->
+          int spaceIndex = tempInputCollection.find(" ");
->
+          int secondWordLength = tempInputCollection.size() - spaceIndex;
->
+          objectsToGet.push_back(tempInputCollection.substr(spaceIndex+1,secondWordLength));
->
+        }
->
+        else{
->
+          objectsToGet.push_back(tempInputCollection);
->
+        }
->
+        objectsToCut.push_back(tempInputCollection);
->
+      }
->
+      else{//pair of objects, need to add them both to objectsToGet
->
+        string obj1;
->
+        string obj2;
->
+        getTwoObjs(tempInputCollection, obj1, obj2);
->
+        string obj2ToGet = getObjToGet(obj2);
->
+        objectsToCut.push_back(tempInputCollection);
->
+        objectsToCut.push_back(obj1);
->
+        objectsToCut.push_back(obj2);
->
+        objectsToGet.push_back(tempInputCollection);
->
+        objectsToGet.push_back(obj1);
->
+        objectsToGet.push_back(obj2ToGet);
->
->
+      }
->
->
->
->
+      //split cut string into parts and store them
->
+      string cutString = cuts_.at(currentCut).getParameter<string> ("cutString");
->
+      vector<string> cutStringVector = splitString(cutString);
->
+      if(cutStringVector.size()!=3 && cutStringVector.size() % 4 !=3){
->
+        cout << "Error: Didn't find the expected number elements in the following cut string: '" << cutString << "'\n";
->
+        exit(0);
->
+      }
->
+      tempCut.numSubcuts = (cutStringVector.size()+1)/4;
->
+      for(int subcutIndex = 0; subcutIndex != tempCut.numSubcuts; subcutIndex++){//loop over all the pieces of the cut combined using &,|
->
+        int indexOffset = 4 * subcutIndex;
->
+        string currentVariableString = cutStringVector.at(indexOffset);
->
+        if(currentVariableString.find("(")==string::npos){
->
+          tempCut.functions.push_back("");//no function was specified
->
+          tempCut.variables.push_back(currentVariableString);// variable to cut on
->
+        }
->
+        else{
->
+          tempCut.functions.push_back(currentVariableString.substr(0,currentVariableString.find("(")));//function comes before the "("
->
+          string tempVariable = currentVariableString.substr(currentVariableString.find("(")+1);//get rest of string
->
+          tempCut.variables.push_back(tempVariable.substr(0,tempVariable.size()-1));//remove trailing ")"
->
+        }
->
+        tempCut.comparativeOperators.push_back(cutStringVector.at(indexOffset+1));// comparison to make
->
+        tempCut.cutValues.push_back(atof(cutStringVector.at(indexOffset+2).c_str()));// threshold value to pass cut
->
+        tempCut.cutStringValues.push_back(cutStringVector.at(indexOffset+2));// string value to pass cut
->
+        if(subcutIndex != 0) tempCut.logicalOperators.push_back(cutStringVector.at(indexOffset-1)); // logical comparison (and, or)
->
+      }
->
->
+      //get number of objects required to pass cut for event to pass
->
+      string numberRequiredString = cuts_.at(currentCut).getParameter<string> ("numberRequired");
->
+      vector<string> numberRequiredVector = splitString(numberRequiredString);
->
+      if(numberRequiredVector.size()!=2){
->
+        cout << "Error: Didn't find two elements in the following number requirement string: '" << numberRequiredString << "'\n";
->
+        exit(0);
->
+      }
-+
+      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 tempCutName;
-+
+      if(cuts_.at(currentCut).exists("alias")){
-+
+        tempCutName = cuts_.at(currentCut).getParameter<string> ("alias");
-+
+        subSubDirName = "After " + cuts_.at(currentCut).getParameter<string> ("alias") + " Cut Applied";
-+
+      }
-+
+      else{
-+
+        //construct string for cutflow table
-+
+        bool plural = numberRequiredInt != 1;
-+
+        string collectionString = plural ? tempInputCollection : tempInputCollection.substr(0, tempInputCollection.size()-1);
-+
+        string cutName =  numberRequiredString + " " + collectionString + " with " + cutString;
-+
+        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;
-+
+    //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++){
-+
+            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{
-+
+      map<string, TH1D*> oneDhistoMap;
-+
+      oneDHistsTmp.push_back(oneDhistoMap);
-+
+      map<string, TH2D*> twoDhistoMap;
-+
+      twoDHistsTmp.push_back(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();
+      else if(numBinsElements == 3){
+        if (currentHistogram.bins.size () == 3)
-<
+          oneDHists_.at(currentChannel)[currentHistogram.name] = directories.at(currentChannel).make<TH1D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, currentHistogram.bins.at(0), currentHistogram.bins.at(1), currentHistogram.bins.at(2));
->
+          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)[currentHistogram.name] = directories.at(currentChannel).make<TH1D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, numBinEdgesX - 1, currentHistogram.variableBinsX.data ());
->
+            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)[currentHistogram.name] = directories.at(currentChannel).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));
->
+          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)[currentHistogram.name] = directories.at(currentChannel).make<TH2D> (TString(currentHistogram.name),channelLabel+" channel: "+currentHistogram.title, numBinEdgesX - 1, currentHistogram.variableBinsX.data (), numBinEdgesY - 1, currentHistogram.variableBinsY.data ());
->
+          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 ());
+      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)[currentHistogram.name]->GetXaxis()->SetBinLabel(bin+1,labelArray.at(bin));
->
+          oneDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->SetBinLabel(bin+1,labelArray.at(bin));
+        else {
-<
+          twoDHists_.at(currentChannel)[currentHistogram.name]->GetYaxis()->SetBinLabel(bin+1,labelArray.at(currentHistogram.bins.at(0)-bin-1));
-<
+          twoDHists_.at(currentChannel)[currentHistogram.name]->GetXaxis()->SetBinLabel(bin+1,labelArray.at(bin));
->
+          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)[currentHistogram.name]->GetXaxis()->CenterTitle();
-<
+        twoDHists_.at(currentChannel)[currentHistogram.name]->GetYaxis()->CenterTitle();
->
+        twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetXaxis()->CenterTitle();
->
+        twoDHists_.at(currentChannel).at(currentDir)[currentHistogram.name]->GetYaxis()->CenterTitle();
+      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 == "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 == "secondary electrons")       currentObject = "secondaryElectrons";
+      else if(currentObject == "electron-trigobj pairs")    currentObject = "electronTrigobjPairs";
+      else if(currentObject == "muon-trigobj pairs")        currentObject = "muonTrigobjPairs";
+      if(histoName == "numPrimaryvertexs"){
+        string newHistoName = histoName + "BeforePileupCorrection";
-<
+        oneDHists_.at(currentChannel)[newHistoName] = directories.at(currentChannel).make<TH1D> (TString(newHistoName),channelLabel+" channel: Number of Selected "+currentObject+" Before Pileup Correction; # "+currentObject, maxNum, 0, maxNum);
->
+        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)[newHistoName] = directories.at(currentChannel).make<TH1D> (TString(newHistoName),channelLabel+" channel: Number of Selected "+currentObject+ " After Pileup Correction; # "+currentObject, maxNum, 0, maxNum);
->
+        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)[histoName] = directories.at(currentChannel).make<TH1D> (TString(histoName),channelLabel+" channel: Number of Selected "+currentObject+"; # "+currentObject, maxNum, 0, maxNum);
->
+        oneDHists_.at(currentChannel).at(currentDir)[histoName] = directories.at(currentDir).make<TH1D> (TString(histoName),channelLabel+" channel: Number of Selected "+currentObject+"; # "+currentObject, maxNum, 0, maxNum);
-<
-<
-<
-<
-<
+    //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++){
-<
+      cut tempCut;
-<
+      //store input collection for cut
-<
+      string tempInputCollection = cuts_.at(currentCut).getParameter<string> ("inputCollection");
-<
+      if(tempInputCollection == "muon-electron pairs") tempInputCollection = "electron-muon pairs";
-<
+      if(tempInputCollection == "jet-electron pairs") tempInputCollection = "electron-jet pairs";
-<
+      if(tempInputCollection == "jet-muon pairs") tempInputCollection = "muon-jet pairs";
-<
+      if(tempInputCollection == "event-track pairs")   tempInputCollection = "track-event pairs";
-<
+      if(tempInputCollection == "secondary muon-muon pairs")   tempInputCollection = "muon-secondary muon pairs";
-<
+      if(tempInputCollection == "secondary electron-electron pairs")   tempInputCollection = "electron-secondary electron pairs";
-<
+      if(tempInputCollection == "trigobj-electron pairs")   tempInputCollection = "electron-trigobj pairs";
-<
+      if(tempInputCollection == "trigobj-muon pairs")   tempInputCollection = "muon-trigobj pairs";
-<
+      tempCut.inputCollection = tempInputCollection;
-<
+      if(tempInputCollection.find("pairs")==string::npos){ //just a single object
-<
+        if(tempInputCollection.find("secondary")!=string::npos){//secondary object
-<
+          int spaceIndex = tempInputCollection.find(" ");
-<
+          int secondWordLength = tempInputCollection.size() - spaceIndex;
-<
+          objectsToGet.push_back(tempInputCollection.substr(spaceIndex+1,secondWordLength));
-<
+        }
-<
+        else{
-<
+          objectsToGet.push_back(tempInputCollection);
-<
+        }
-<
+        objectsToCut.push_back(tempInputCollection);
-<
+      }
-<
+      else{//pair of objects, need to add them both to objectsToGet
-<
+        string obj1;
-<
+        string obj2;
-<
+        getTwoObjs(tempInputCollection, obj1, obj2);
-<
+        string obj2ToGet = getObjToGet(obj2);
-<
+        objectsToCut.push_back(tempInputCollection);
-<
+        objectsToCut.push_back(obj1);
-<
+        objectsToCut.push_back(obj2);
-<
+        objectsToGet.push_back(tempInputCollection);
-<
+        objectsToGet.push_back(obj1);
-<
+        objectsToGet.push_back(obj2ToGet);
-<
-<
+      }
-<
-<
-<
-<
+      //split cut string into parts and store them
-<
+      string cutString = cuts_.at(currentCut).getParameter<string> ("cutString");
-<
+      vector<string> cutStringVector = splitString(cutString);
-<
+      if(cutStringVector.size()!=3 && cutStringVector.size() % 4 !=3){
-<
+        cout << "Error: Didn't find the expected number elements in the following cut string: '" << cutString << "'\n";
-<
+        exit(0);
-<
+      }
-<
+      tempCut.numSubcuts = (cutStringVector.size()+1)/4;
-<
+      for(int subcutIndex = 0; subcutIndex != tempCut.numSubcuts; subcutIndex++){//loop over all the pieces of the cut combined using &,|
-<
+        int indexOffset = 4 * subcutIndex;
-<
+        string currentVariableString = cutStringVector.at(indexOffset);
-<
+        if(currentVariableString.find("(")==string::npos){
-<
+          tempCut.functions.push_back("");//no function was specified
-<
+          tempCut.variables.push_back(currentVariableString);// variable to cut on
-<
+        }
-<
+        else{
-<
+          tempCut.functions.push_back(currentVariableString.substr(0,currentVariableString.find("(")));//function comes before the "("
-<
+          string tempVariable = currentVariableString.substr(currentVariableString.find("(")+1);//get rest of string
-<
+          tempCut.variables.push_back(tempVariable.substr(0,tempVariable.size()-1));//remove trailing ")"
-<
+        }
-<
+        tempCut.comparativeOperators.push_back(cutStringVector.at(indexOffset+1));// comparison to make
-<
+        tempCut.cutValues.push_back(atof(cutStringVector.at(indexOffset+2).c_str()));// threshold value to pass cut
-<
+        tempCut.cutStringValues.push_back(cutStringVector.at(indexOffset+2));// string value to pass cut
-<
+        if(subcutIndex != 0) tempCut.logicalOperators.push_back(cutStringVector.at(indexOffset-1)); // logical comparison (and, or)
-<
+      }
-<
-<
+      //get number of objects required to pass cut for event to pass
-<
+      string numberRequiredString = cuts_.at(currentCut).getParameter<string> ("numberRequired");
-<
+      vector<string> numberRequiredVector = splitString(numberRequiredString);
-<
+      if(numberRequiredVector.size()!=2){
-<
+        cout << "Error: Didn't find two elements in the following number requirement string: '" << numberRequiredString << "'\n";
-<
+        exit(0);
-<
+      }
-<
-<
+      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
-<
-<
-<
+      string tempCutName;
-<
+      if(cuts_.at(currentCut).exists("alias")){
-<
+        tempCutName = cuts_.at(currentCut).getParameter<string> ("alias");
-<
+      }
-<
+      else{
-<
+        //construct string for cutflow table
-<
+        bool plural = numberRequiredInt != 1;
-<
+        string collectionString = plural ? tempInputCollection : tempInputCollection.substr(0, tempInputCollection.size()-1);
-<
+        string cutName =  numberRequiredString + " " + collectionString + " with " + cutString;
-<
+        tempCutName = cutName;
-<
+      }
-<
+      tempCut.name = tempCutName;
-<
-<
+      tempChannel.cuts.push_back(tempCut);
-<
-<
-<
+    }//end loop over cuts
-<
->
+   }//end of loop over directories
+    channels.push_back(tempChannel);
+    tempChannel.cuts.clear();
-–
+  }//end loop over channels
+  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;
+    //loop over all cuts
-–
+    for(uint currentCutIndex = 0; currentCutIndex != currentChannel.cuts.size(); currentCutIndex++){
+      cut currentCut = currentChannel.cuts.at(currentCutIndex);
-–
+      for(uint currentObjectIndex = 0; currentObjectIndex != objectsToCut.size(); currentObjectIndex++){
+        string currentObject = objectsToCut.at(currentObjectIndex);
+        int flagsForPairCutsIndex = max(int(currentCutIndex-1),0);
-<
+        if(currentObject == "jets") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,jets.product(),"jets");
-<
-<
+        else if(currentObject == "muons") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,muons.product(),"muons");
-<
->
+        if     (currentObject == "jets")            setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,jets.product(),"jets");
->
+        else if(currentObject == "secondary jets")  setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,jets.product(),"secondary jets");
->
+        else if(currentObject == "muons")           setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,muons.product(),"muons");
+        else if(currentObject == "secondary muons") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,muons.product(),"secondary muons");
+        else if(currentObject == "secondary electrons") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,electrons.product(),"secondary electrons");
+        else if(currentObject == "electrons") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,electrons.product(),"electrons");
+    }//end loop over all cuts
-–
-–
+    //use cumulative flags to apply cuts at event level
+    bool eventPassedAllCuts = true;
-<
->
+    //a vector to store cumulative cut descisions after each cut.
->
+    vector<bool> eventPassedPreviousCuts;
+    //apply trigger (true if none were specified)
+    eventPassedAllCuts = eventPassedAllCuts && triggerDecision;
+      bool cutDecision = evaluateComparison(numberPassing,currentCut.eventComparativeOperator,currentCut.numberRequired);
+      cutFlows_.at(currentChannelIndex)->at (currentCut.name) = cutDecision;
+      eventPassedAllCuts = eventPassedAllCuts && cutDecision;
-<
->
+      eventPassedPreviousCuts.push_back(eventPassedAllCuts);
-–
+    double scaleFactor = masterScaleFactor;
+    muonScaleFactor_ = electronScaleFactor_ = 1.0;
+    cutFlows_.at(currentChannelIndex)->fillCutFlow(scaleFactor);
-–
+    if(!eventPassedAllCuts)continue;
+    if (printEventInfo_) {
+      // Write information about event to screen, for testing purposes.
+           << endl;
-–
+    //filling 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(currentHistogram.name);
-–
-–
+        if(currentHistogram.inputCollection == "jets") fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muons") fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "secondary muons") fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "secondary electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-–
+                                                                                       cumulativeFlags.at("muons").back(),cumulativeFlags.at("muons").back(), \
-–
+                                                                                       cumulativeFlags.at("muon-muon pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-–
+                                                                                       cumulativeFlags.at("muons").back(),cumulativeFlags.at("secondary muons").back(), \
-–
+                                                                                       cumulativeFlags.at("muon-secondary muon pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(),\
-–
+                                                                                               cumulativeFlags.at("electrons").back(),cumulativeFlags.at("electrons").back(),\
-–
+                                                                                               cumulativeFlags.at("electron-electron pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "jet-jet pairs") fill1DHistogram(histo,currentHistogram,jets.product(),jets.product(),\
-–
+                                                                                               cumulativeFlags.at("jets").back(),cumulativeFlags.at("jets").back(),\
-–
+                                                                                               cumulativeFlags.at("jet-jet pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-–
+                                                                                       cumulativeFlags.at("electrons").back(),cumulativeFlags.at("secondary electrons").back(), \
-–
+                                                                                       cumulativeFlags.at("electron-secondary electron pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-muon pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),muons.product(), \
-–
+                                                                                           cumulativeFlags.at("electrons").back(),cumulativeFlags.at("muons").back(),
-–
+                                                                                           cumulativeFlags.at("electron-muon pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-jet pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),jets.product(), \
-–
+                                                                                           cumulativeFlags.at("electrons").back(),cumulativeFlags.at("jets").back(),
-–
+                                                                                           cumulativeFlags.at("electron-jet pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-jet pairs") fill1DHistogram(histo,currentHistogram, muons.product(),jets.product(), \
-–
+                                                                                           cumulativeFlags.at("muons").back(),cumulativeFlags.at("jets").back(),
-–
+                                                                                           cumulativeFlags.at("muon-jet pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-track pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),tracks.product(),
-–
+                                                                                            cumulativeFlags.at("electrons").back(),cumulativeFlags.at("tracks").back(),
-–
+                                                                                            cumulativeFlags.at("electron-track pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-track pairs") fill1DHistogram(histo,currentHistogram, muons.product(),tracks.product(),
-–
+                                                                                        cumulativeFlags.at("muons").back(),cumulativeFlags.at("tracks").back(),
-–
+                                                                                        cumulativeFlags.at("muon-track pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-tau pairs") fill1DHistogram(histo,currentHistogram, muons.product(),taus.product(),
-–
+                                                                                      cumulativeFlags.at("muons").back(),cumulativeFlags.at("taus").back(),
-–
+                                                                                      cumulativeFlags.at("muon-tau pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "tau-tau pairs") fill1DHistogram(histo,currentHistogram, taus.product(),taus.product(),
-–
+                                                                                     cumulativeFlags.at("taus").back(),cumulativeFlags.at("taus").back(),
-–
+                                                                                     cumulativeFlags.at("tau-tau pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "tau-track pairs") fill1DHistogram(histo,currentHistogram, taus.product(),tracks.product(),
-–
+                                                                                     cumulativeFlags.at("taus").back(),cumulativeFlags.at("tracks").back(),
-–
+                                                                                     cumulativeFlags.at("tau-track pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),trigobjs.product(),
-–
+                                                                                              cumulativeFlags.at("electrons").back(),cumulativeFlags.at("trigobjs").back(),
-–
+                                                                                              cumulativeFlags.at("electron-trigobj pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill1DHistogram(histo,currentHistogram, muons.product(),trigobjs.product(),
-–
+                                                                                          cumulativeFlags.at("muons").back(),cumulativeFlags.at("trigobjs").back(),
-–
+                                                                                          cumulativeFlags.at("muon-trigobj pairs").back(),scaleFactor);
-–
-–
+        else if(currentHistogram.inputCollection == "events") fill1DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "taus") fill1DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "mets") fill1DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "tracks") fill1DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "genjets") fill1DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "mcparticles") fill1DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "primaryvertexs") fill1DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "bxlumis") fill1DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "photons") fill1DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "superclusters") fill1DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "trigobjs") fill1DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill1DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").back(),scaleFactor);
-–
+      }
-–
+      else if(currentHistogram.inputVariables.size() == 2){
-–
+        TH2D* histo;
-–
+        histo = twoDHists_.at(currentChannelIndex).at(currentHistogram.name);
-–
-–
+        if(currentHistogram.inputCollection == "jets") fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muons") fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "secondary muons") fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-–
+                                                                                       cumulativeFlags.at("muons").back(),cumulativeFlags.at("muons").back(), \
-–
+                                                                                       cumulativeFlags.at("muon-muon pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
-–
+                                                                                       cumulativeFlags.at("muons").back(),cumulativeFlags.at("secondary muons").back(), \
-–
+                                                                                       cumulativeFlags.at("muon-secondary muon pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "secondary electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-–
+                                                                                               cumulativeFlags.at("electrons").back(),cumulativeFlags.at("electrons").back(), \
-–
+                                                                                               cumulativeFlags.at("electron-electron pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "jet-jet pairs") fill2DHistogram(histo,currentHistogram,jets.product(),jets.product(), \
-–
+                                                                                               cumulativeFlags.at("jets").back(),cumulativeFlags.at("jets").back(), \
-–
+                                                                                               cumulativeFlags.at("jet-jet pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
-–
+                                                                                       cumulativeFlags.at("electrons").back(),cumulativeFlags.at("secondary electrons").back(), \
-–
+                                                                                       cumulativeFlags.at("electron-secondary electron pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-muon pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),muons.product(), \
-–
+                                                                                           cumulativeFlags.at("electrons").back(),cumulativeFlags.at("muons").back(), \
-–
+                                                                                           cumulativeFlags.at("electron-muon pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-jet pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),jets.product(), \
-–
+                                                                                           cumulativeFlags.at("electrons").back(),cumulativeFlags.at("jets").back(), \
-–
+                                                                                           cumulativeFlags.at("electron-jet pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-jet pairs") fill2DHistogram(histo,currentHistogram,muons.product(),jets.product(), \
-–
+                                                                                           cumulativeFlags.at("muons").back(),cumulativeFlags.at("jets").back(), \
-–
+                                                                                           cumulativeFlags.at("muon-jet pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-track pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),tracks.product(),
-–
+                                                                                            cumulativeFlags.at("electrons").back(),cumulativeFlags.at("tracks").back(),
-–
+                                                                                            cumulativeFlags.at("electron-track pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-track pairs") fill2DHistogram(histo,currentHistogram,muons.product(),tracks.product(),
-–
+                                                                                        cumulativeFlags.at("muons").back(),cumulativeFlags.at("tracks").back(),
-–
+                                                                                        cumulativeFlags.at("muon-track pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-tau pairs") fill2DHistogram(histo,currentHistogram,muons.product(),taus.product(),
-–
+                                                                                      cumulativeFlags.at("muons").back(),cumulativeFlags.at("taus").back(),
-–
+                                                                                      cumulativeFlags.at("muon-tau pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "tau-tau pairs") fill2DHistogram(histo,currentHistogram,taus.product(),taus.product(),
-–
+                                                                                     cumulativeFlags.at("taus").back(),cumulativeFlags.at("taus").back(),
-–
+                                                                                     cumulativeFlags.at("tau-tau pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "tau-track pairs") fill2DHistogram(histo,currentHistogram,taus.product(),tracks.product(),
-–
+                                                                                     cumulativeFlags.at("taus").back(),cumulativeFlags.at("tracks").back(),
-–
+                                                                                     cumulativeFlags.at("tau-track pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),trigobjs.product(),
-–
+                                                                                              cumulativeFlags.at("electrons").back(),cumulativeFlags.at("trigobjs").back(),
-–
+                                                                                              cumulativeFlags.at("electron-trigobj pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill2DHistogram(histo,currentHistogram,muons.product(),trigobjs.product(),
-–
+                                                                                          cumulativeFlags.at("muons").back(),cumulativeFlags.at("trigobjs").back(),
-–
+                                                                                          cumulativeFlags.at("muon-trigobj pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "events") fill2DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "taus") fill2DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "mets") fill2DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "tracks") fill2DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "track-event pairs") fill2DHistogram(histo,currentHistogram,tracks.product(),events.product(),
-–
+                                                                                         cumulativeFlags.at("tracks").back(),cumulativeFlags.at("events").back(),
-–
+                                                                                         cumulativeFlags.at("track-event pairs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "genjets") fill2DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "mcparticles") fill2DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "primaryvertexs") fill2DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "bxlumis") fill2DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "photons") fill2DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "superclusters") fill2DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "trigobjs") fill2DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").back(),scaleFactor);
-–
+        else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill2DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").back(),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 == "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).back();
-<
+        int numToPlot = 0;
-<
+        for (uint currentFlag = 0; currentFlag != lastCutFlags.size(); currentFlag++){
-<
+          if(lastCutFlags.at(currentFlag)) numToPlot++;
-<
+        }
-<
+        if(objectToPlot == "primaryvertexs"){
-<
+          oneDHists_.at(currentChannelIndex).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
-<
+          oneDHists_.at(currentChannelIndex).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
-<
+        }
-<
+        else {
-<
+          oneDHists_.at(currentChannelIndex).at(histoName)->Fill(numToPlot,scaleFactor);
-<
+        }
-<
+      }
-<
+      else if(objectToPlot == "jets") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(jets->size(),scaleFactor);
-<
+      else if(objectToPlot == "muons") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(muons->size(),scaleFactor);
-<
+      else if(objectToPlot == "secondaryMuons") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(muons->size(),scaleFactor);
-<
+      else if(objectToPlot == "dimuonPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
-<
+      else if(objectToPlot == "muonSecondaryMuonPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
-<
+      else if(objectToPlot == "electrons") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size(),scaleFactor);
-<
+      else if(objectToPlot == "secondaryElectrons") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size(),scaleFactor);
-<
+      else if(objectToPlot == "dielectronPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
-<
+      else if(objectToPlot == "electronSecondaryElectronPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
-<
+      else if(objectToPlot == "electronMuonPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size()*muons->size(),scaleFactor);
-<
+      else if(objectToPlot == "electronJetPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size()*jets->size(),scaleFactor);
-<
+      else if(objectToPlot == "muonJetPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(muons->size()*jets->size(),scaleFactor);
-<
+      else if(objectToPlot == "electronTrackPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size()*tracks->size(),scaleFactor);
-<
+      else if(objectToPlot == "electronTrigobjPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(electrons->size()*trigobjs->size(),scaleFactor);
-<
+      else if(objectToPlot == "muonTrigobjPairs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(muons->size()*trigobjs->size(),scaleFactor);
-<
+      else if(objectToPlot == "events") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(events->size(),scaleFactor);
-<
+      else if(objectToPlot == "taus") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(taus->size(),scaleFactor);
-<
+      else if(objectToPlot == "mets") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(mets->size(),scaleFactor);
-<
+      else if(objectToPlot == "tracks") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(tracks->size(),scaleFactor);
-<
+      else if(objectToPlot == "genjets") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(genjets->size(),scaleFactor);
-<
+      else if(objectToPlot == "mcparticles") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(mcparticles->size(),scaleFactor);
-<
+      else if(objectToPlot == "bxlumis") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(bxlumis->size(),scaleFactor);
-<
+      else if(objectToPlot == "photons") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(photons->size(),scaleFactor);
-<
+      else if(objectToPlot == "superclusters") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(superclusters->size(),scaleFactor);
-<
+      else if(objectToPlot == "trigobjs") oneDHists_.at(currentChannelIndex).at(histoName)->Fill(trigobjs->size(),scaleFactor);
-<
+      else if(objectToPlot == "primaryvertexs"){
-<
+        oneDHists_.at(currentChannelIndex).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
-<
+        oneDHists_.at(currentChannelIndex).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
-<
+      }
-<
+      if(objectToPlot == "stops" && datasetType_ == "signalMC") oneDHists_.at(currentChannelIndex).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
+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 == "tightIDdisplaced"){
+    value = object->isGlobalMuon > 0                \
-+
+      && object->isPFMuon > 0                        \
+      && object->normalizedChi2 < 10                \
+                                  && object->numberOfValidMuonHits > 0        \
+      && object->numberOfMatchedStations > 1        \
+  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;
+  if(variable == "deltaPhi") value = fabs(deltaPhi(object1->phi,object2->phi));
+  else if(variable == "deltaEta") value = fabs(object1->eta - object2->eta);
-+
+  else if(variable == "absDeltaPt") value = fabs(object1->pt - object2->pt);
+  else if(variable == "deltaR") value = deltaR(object1->eta,object1->phi,object2->eta,object2->phi);
+  else if(variable == "invMass"){
+    TLorentzVector fourVector1(object1->px, object1->py, object1->pz, object1->energy);
+    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.64 by ahart, Sun May 19 00:04:50 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.64 by ahart, Sun May 19 00:04:50 2013 UTC vs.
Revision 1.73 by lantonel, Thu May 30 18:09:58 2013 UTC