ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents):
Revision 1.72 by lantonel, Wed May 29 20:15:01 2013 UTC vs.
Revision 1.73 by lantonel, Thu May 30 18:09:58 2013 UTC

#	Line 57 | Line 57 | OSUAnalysis::OSUAnalysis (const edm::Par
57
58		//always get vertex collection so we can assign the primary vertex in the event
59		objectsToGet.push_back("primaryvertexs");
60	–
61	–	//always make the plot of number of primary vertices (to check pile-up reweighting)
60		objectsToPlot.push_back("primaryvertexs");
61	+	objectsToCut.push_back("primaryvertexs");
62	+
63
64		//always get the MC particles to do GEN-matching
65		objectsToGet.push_back("mcparticles");
#	Line 108 | Line 108 | OSUAnalysis::OSUAnalysis (const edm::Par
108
109		}
110
111	+
112		vector<edm::ParameterSet> histogramList_  (histogramSets_.at(currentHistogramSet).getParameter<vector<edm::ParameterSet> >("histograms"));
113
114		for(uint currentHistogram = 0; currentHistogram != histogramList_.size(); currentHistogram++){
#	Line 811 | Line 812 | OSUAnalysis::analyze (const edm::Event &
812		scaleFactor *= electronScaleFactor_;
813
814		cutFlows_.at(currentChannelIndex)->fillCutFlow(scaleFactor);
815	<	if(!(GetPlotsAfterEachCut_ || eventPassedAllCuts)) continue;
815	>
816
817		if (printEventInfo_) {
818		// Write information about event to screen, for testing purposes.
#	Line 821 | Line 822 | OSUAnalysis::analyze (const edm::Event &
822		<< "  event=" << events->at(0).evt
823		<< endl;
824		}
824	–	//filling histograms
825	–	for(uint currentCut = 0; currentCut != oneDHists_.at(currentChannelIndex).size(); currentCut++){//loop over all the directories in each channel.
826	–	uint currentDir;
827	–	if(!GetPlotsAfterEachCut_){ currentDir =  currentChannel.cuts.size() - oneDHists_.at(currentChannelIndex).size();}//if GetPlotsAfterEachCut_ is true, set currentDir point to the lat cut.
828	–	else{
829	–	currentDir = currentCut;
830	–	}
831	–	if(eventPassedPreviousCuts.at(currentDir)){
832	–	for (uint histogramIndex = 0; histogramIndex != histograms.size(); histogramIndex++){
833	–	histogram currentHistogram = histograms.at(histogramIndex);
834	–
835	–	if(currentHistogram.inputVariables.size() == 1){
836	–	TH1D* histo;
837	–	histo = oneDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
838	–	if     (currentHistogram.inputCollection == "jets")            fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
839	–	else if(currentHistogram.inputCollection == "secondary jets")  fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
840	–	else if(currentHistogram.inputCollection == "muons")           fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
841	–	else if(currentHistogram.inputCollection == "secondary muons") fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
842	–	else if(currentHistogram.inputCollection == "secondary electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
843	–	else if(currentHistogram.inputCollection == "muon-muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
844	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
845	–	cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
846	–	else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
847	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
848	–	cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
849	–	else if(currentHistogram.inputCollection == "electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
850	–	else if(currentHistogram.inputCollection == "electron-electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(),\
851	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir),\
852	–	cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
853	–	else if(currentHistogram.inputCollection == "jet-jet pairs") fill1DHistogram(histo,currentHistogram,jets.product(),jets.product(),\
854	–	cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir),\
855	–	cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
856	–	else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
857	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
858	–	cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
859	–	else if(currentHistogram.inputCollection == "electron-muon pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),muons.product(), \
860	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir),
861	–	cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
862	–	else if(currentHistogram.inputCollection == "electron-jet pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),jets.product(), \
863	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
864	–	cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
865	–	else if(currentHistogram.inputCollection == "muon-jet pairs") fill1DHistogram(histo,currentHistogram, muons.product(),jets.product(), \
866	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
867	–	cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
868	–	else if(currentHistogram.inputCollection == "electron-track pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),tracks.product(),
869	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
870	–	cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
871	–	else if(currentHistogram.inputCollection == "muon-track pairs") fill1DHistogram(histo,currentHistogram, muons.product(),tracks.product(),
872	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
873	–	cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
874	–	else if(currentHistogram.inputCollection == "muon-tau pairs") fill1DHistogram(histo,currentHistogram, muons.product(),taus.product(),
875	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
876	–	cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
877	–	else if(currentHistogram.inputCollection == "tau-tau pairs") fill1DHistogram(histo,currentHistogram, taus.product(),taus.product(),
878	–	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
879	–	cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
880	–	else if(currentHistogram.inputCollection == "tau-track pairs") fill1DHistogram(histo,currentHistogram, taus.product(),tracks.product(),
881	–	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
882	–	cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
883	–	else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),trigobjs.product(),
884	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
885	–	cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
886	–	else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill1DHistogram(histo,currentHistogram, muons.product(),trigobjs.product(),
887	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
888	–	cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
889	–
890	–	else if(currentHistogram.inputCollection == "events") fill1DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
891	–	else if(currentHistogram.inputCollection == "taus") fill1DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
892	–	else if(currentHistogram.inputCollection == "mets") fill1DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
893	–	else if(currentHistogram.inputCollection == "tracks") fill1DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
894	–	else if(currentHistogram.inputCollection == "genjets") fill1DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
895	–	else if(currentHistogram.inputCollection == "mcparticles") fill1DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
896	–	else if(currentHistogram.inputCollection == "primaryvertexs") fill1DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
897	–	else if(currentHistogram.inputCollection == "bxlumis") fill1DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
898	–	else if(currentHistogram.inputCollection == "photons") fill1DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
899	–	else if(currentHistogram.inputCollection == "superclusters") fill1DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
900	–	else if(currentHistogram.inputCollection == "trigobjs") fill1DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
901	–	else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill1DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
902	–	}
903	–	else if(currentHistogram.inputVariables.size() == 2){
904	–	TH2D* histo;
905	–	histo = twoDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
906	–	if     (currentHistogram.inputCollection == "jets")            fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
907	–	else if(currentHistogram.inputCollection == "secondary jets")  fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
908	–	else if(currentHistogram.inputCollection == "muons")           fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
909	–	else if(currentHistogram.inputCollection == "secondary muons") fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
910	–	else if(currentHistogram.inputCollection == "muon-muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
911	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
912	–	cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
913	–	else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
914	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
915	–	cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
916	–	else if(currentHistogram.inputCollection == "electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
917	–	else if(currentHistogram.inputCollection == "secondary electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
918	–	else if(currentHistogram.inputCollection == "electron-electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
919	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir), \
920	–	cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
921	–	else if(currentHistogram.inputCollection == "jet-jet pairs") fill2DHistogram(histo,currentHistogram,jets.product(),jets.product(), \
922	–	cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
923	–	cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
924	–	else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
925	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
926	–	cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
927	–	else if(currentHistogram.inputCollection == "electron-muon pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),muons.product(), \
928	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
929	–	cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
930	–	else if(currentHistogram.inputCollection == "electron-jet pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),jets.product(), \
931	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
932	–	cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
933	–	else if(currentHistogram.inputCollection == "muon-jet pairs") fill2DHistogram(histo,currentHistogram,muons.product(),jets.product(), \
934	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
935	–	cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
936	–	else if(currentHistogram.inputCollection == "electron-track pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),tracks.product(),
937	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
938	–	cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
939	–	else if(currentHistogram.inputCollection == "muon-track pairs") fill2DHistogram(histo,currentHistogram,muons.product(),tracks.product(),
940	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
941	–	cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
942	–	else if(currentHistogram.inputCollection == "muon-tau pairs") fill2DHistogram(histo,currentHistogram,muons.product(),taus.product(),
943	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
944	–	cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
945	–	else if(currentHistogram.inputCollection == "tau-tau pairs") fill2DHistogram(histo,currentHistogram,taus.product(),taus.product(),
946	–	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
947	–	cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
948	–	else if(currentHistogram.inputCollection == "tau-track pairs") fill2DHistogram(histo,currentHistogram,taus.product(),tracks.product(),
949	–	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
950	–	cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
951	–	else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),trigobjs.product(),
952	–	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
953	–	cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
954	–	else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill2DHistogram(histo,currentHistogram,muons.product(),trigobjs.product(),
955	–	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
956	–	cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
957	–	else if(currentHistogram.inputCollection == "events") fill2DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
958	–	else if(currentHistogram.inputCollection == "taus") fill2DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
959	–	else if(currentHistogram.inputCollection == "mets") fill2DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
960	–	else if(currentHistogram.inputCollection == "tracks") fill2DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
961	–	else if(currentHistogram.inputCollection == "track-event pairs") fill2DHistogram(histo,currentHistogram,tracks.product(),events.product(),
962	–	cumulativeFlags.at("tracks").at(currentDir),cumulativeFlags.at("events").at(currentDir),
963	–	cumulativeFlags.at("track-event pairs").at(currentDir),scaleFactor);
964	–	else if(currentHistogram.inputCollection == "genjets") fill2DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
965	–	else if(currentHistogram.inputCollection == "mcparticles") fill2DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
966	–	else if(currentHistogram.inputCollection == "primaryvertexs") fill2DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
967	–	else if(currentHistogram.inputCollection == "bxlumis") fill2DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
968	–	else if(currentHistogram.inputCollection == "photons") fill2DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
969	–	else if(currentHistogram.inputCollection == "superclusters") fill2DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
970	–	else if(currentHistogram.inputCollection == "trigobjs") fill2DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
971	–	else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill2DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
972	–	}
973	–	}
974	–
825
826
977	–	//fills histograms with the sizes of collections
978	–	for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++){
827
828	<	string currentObject = objectsToPlot.at(currentObjectIndex);
829	<	string objectToPlot = "";
982	<
983	<	// Name of objectToPlot here must match the name specified in OSUAnalysis::OSUAnalysis().
984	<	if(currentObject == "muon-muon pairs")                         objectToPlot = "dimuonPairs";
985	<	else if(currentObject == "electron-electron pairs")            objectToPlot = "dielectronPairs";
986	<	else if(currentObject == "electron-muon pairs")                objectToPlot = "electronMuonPairs";
987	<	else if(currentObject == "electron-jet pairs")                 objectToPlot = "electronJetPairs";
988	<	else if(currentObject == "muon-jet pairs")                     objectToPlot = "muonJetPairs";
989	<	else if(currentObject == "jet-jet pairs")                      objectToPlot = "dijetPairs";
990	<	else if(currentObject == "secondary jets")                     objectToPlot = "secondaryJets";
991	<	else if(currentObject == "electron-track pairs")               objectToPlot = "electronTrackPairs";
992	<	else if(currentObject == "muon-track pairs")                   objectToPlot = "muonTrackPairs";
993	<	else if(currentObject == "muon-tau pairs")                     objectToPlot = "muonTauPairs";
994	<	else if(currentObject == "tau-tau pairs")                      objectToPlot = "ditauPairs";
995	<	else if(currentObject == "tau-track pairs")                    objectToPlot = "tauTrackPairs";
996	<	else if(currentObject == "track-event pairs")                  objectToPlot = "trackEventPairs";
997	<	else if(currentObject == "muon-secondary muon pairs")          objectToPlot = "muonSecondaryMuonPairs";
998	<	else if(currentObject == "secondary muons")                    objectToPlot = "secondaryMuons";
999	<	else if(currentObject == "electron-secondary electron pairs")  objectToPlot = "electronSecondaryElectronPairs";
1000	<	else if(currentObject == "secondary electrons")                objectToPlot = "secondaryElectrons";
1001	<	else if(currentObject == "electron-trigobj pairs")             objectToPlot = "electronTrigobjPairs";
1002	<	else if(currentObject == "muon-trigobj pairs")                 objectToPlot = "muonTrigobjPairs";
1003	<	else objectToPlot = currentObject;
1004	<
1005	<	string tempCurrentObject = objectToPlot;
1006	<	tempCurrentObject.at(0) = toupper(tempCurrentObject.at(0));
1007	<	string histoName = "num" + tempCurrentObject;
1008	<
1009	<	//set position of primary vertex in event, in order to calculate quantities relative to it
1010	<	if(find(objectsToCut.begin(), objectsToCut.end(), currentObject) != objectsToCut.end()) {
1011	<	vector<bool> lastCutFlags = cumulativeFlags.at(currentObject).at(currentDir);
1012	<	int numToPlot = 0;
1013	<	for (uint currentFlag = 0; currentFlag != lastCutFlags.size(); currentFlag++){
1014	<	if(lastCutFlags.at(currentFlag)) numToPlot++;
1015	<	}
1016	<	if(objectToPlot == "primaryvertexs"){
1017	<	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
1018	<	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
1019	<	}
1020	<	else {
1021	<	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(numToPlot,scaleFactor);
1022	<	}
1023	<	}
1024	<	else if(objectToPlot == "jets")           oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(jets->size(),scaleFactor);
1025	<	else if(objectToPlot == "secondaryJets")  oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(jets->size(),scaleFactor);
1026	<	else if(objectToPlot == "muons")          oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size(),scaleFactor);
1027	<	else if(objectToPlot == "secondaryMuons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size(),scaleFactor);
1028	<	else if(objectToPlot == "dimuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
1029	<	else if(objectToPlot == "muonSecondaryMuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*(muons->size()-1)/2,scaleFactor);
1030	<	else if(objectToPlot == "electrons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size(),scaleFactor);
1031	<	else if(objectToPlot == "secondaryElectrons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size(),scaleFactor);
1032	<	else if(objectToPlot == "dielectronPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
1033	<	else if(objectToPlot == "electronSecondaryElectronPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*(electrons->size()-1)/2,scaleFactor);
1034	<	else if(objectToPlot == "electronMuonPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*muons->size(),scaleFactor);
1035	<	else if(objectToPlot == "electronJetPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*jets->size(),scaleFactor);
1036	<	else if(objectToPlot == "muonJetPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*jets->size(),scaleFactor);
1037	<	else if(objectToPlot == "electronTrackPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*tracks->size(),scaleFactor);
1038	<	else if(objectToPlot == "electronTrigobjPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(electrons->size()*trigobjs->size(),scaleFactor);
1039	<	else if(objectToPlot == "muonTrigobjPairs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(muons->size()*trigobjs->size(),scaleFactor);
1040	<	else if(objectToPlot == "events") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(events->size(),scaleFactor);
1041	<	else if(objectToPlot == "taus") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(taus->size(),scaleFactor);
1042	<	else if(objectToPlot == "mets") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(mets->size(),scaleFactor);
1043	<	else if(objectToPlot == "tracks") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(tracks->size(),scaleFactor);
1044	<	else if(objectToPlot == "genjets") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(genjets->size(),scaleFactor);
1045	<	else if(objectToPlot == "mcparticles") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(mcparticles->size(),scaleFactor);
1046	<	else if(objectToPlot == "bxlumis") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(bxlumis->size(),scaleFactor);
1047	<	else if(objectToPlot == "photons") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(photons->size(),scaleFactor);
1048	<	else if(objectToPlot == "superclusters") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(superclusters->size(),scaleFactor);
1049	<	else if(objectToPlot == "trigobjs") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(trigobjs->size(),scaleFactor);
1050	<	else if(objectToPlot == "primaryvertexs"){
1051	<	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
1052	<	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
1053	<	}
1054	<	if(objectToPlot == "stops" && datasetType_ == "signalMC") oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(stops->size(),scaleFactor);
828	>	//filling histograms
829	>	for(uint currentCut = 0; currentCut != oneDHists_.at(currentChannelIndex).size(); currentCut++){//loop over all the directories in each channel.
830
831	<	} // end for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++)
831	>	uint currentDir;
832	>	if(!GetPlotsAfterEachCut_){ currentDir =  currentChannel.cuts.size() - oneDHists_.at(currentChannelIndex).size();}//if GetPlotsAfterEachCut_ is true, set currentDir point to the last cut.
833	>	else currentDir = currentCut;
834	>
835	>	if(!eventPassedPreviousCuts.at(currentDir)) continue;
836	>
837	>
838	>	for (uint histogramIndex = 0; histogramIndex != histograms.size(); histogramIndex++){
839	>	histogram currentHistogram = histograms.at(histogramIndex);
840	>
841	>	if(currentHistogram.inputVariables.size() == 1){
842	>	TH1D* histo;
843	>	histo = oneDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
844	>
845	>	if     (currentHistogram.inputCollection == "jets")            fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
846	>	else if(currentHistogram.inputCollection == "secondary jets")  fill1DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
847	>	else if(currentHistogram.inputCollection == "muons")           fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
848	>	else if(currentHistogram.inputCollection == "secondary muons") fill1DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
849	>	else if(currentHistogram.inputCollection == "secondary electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
850	>	else if(currentHistogram.inputCollection == "muon-muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
851	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
852	>	cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
853	>	else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill1DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
854	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
855	>	cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
856	>	else if(currentHistogram.inputCollection == "electrons") fill1DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
857	>	else if(currentHistogram.inputCollection == "electron-electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
858	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir), \
859	>	cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
860	>	else if(currentHistogram.inputCollection == "jet-jet pairs") fill1DHistogram(histo,currentHistogram,jets.product(),jets.product(), \
861	>	cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
862	>	cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
863	>	else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill1DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
864	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
865	>	cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
866	>	else if(currentHistogram.inputCollection == "electron-muon pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),muons.product(), \
867	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir),
868	>	cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
869	>	else if(currentHistogram.inputCollection == "electron-jet pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),jets.product(), \
870	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
871	>	cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
872	>	else if(currentHistogram.inputCollection == "muon-jet pairs") fill1DHistogram(histo,currentHistogram, muons.product(),jets.product(), \
873	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
874	>	cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
875	>	else if(currentHistogram.inputCollection == "electron-track pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),tracks.product(),
876	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
877	>	cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
878	>	else if(currentHistogram.inputCollection == "muon-track pairs") fill1DHistogram(histo,currentHistogram, muons.product(),tracks.product(),
879	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
880	>	cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
881	>	else if(currentHistogram.inputCollection == "muon-tau pairs") fill1DHistogram(histo,currentHistogram, muons.product(),taus.product(),
882	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
883	>	cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
884	>	else if(currentHistogram.inputCollection == "tau-tau pairs") fill1DHistogram(histo,currentHistogram, taus.product(),taus.product(),
885	>	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
886	>	cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
887	>	else if(currentHistogram.inputCollection == "tau-track pairs") fill1DHistogram(histo,currentHistogram, taus.product(),tracks.product(),
888	>	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
889	>	cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
890	>	else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill1DHistogram(histo,currentHistogram, electrons.product(),trigobjs.product(),
891	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
892	>	cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
893	>	else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill1DHistogram(histo,currentHistogram, muons.product(),trigobjs.product(),
894	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
895	>	cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
896	>
897	>	else if(currentHistogram.inputCollection == "events") fill1DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
898	>	else if(currentHistogram.inputCollection == "taus") fill1DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
899	>	else if(currentHistogram.inputCollection == "mets") fill1DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
900	>	else if(currentHistogram.inputCollection == "tracks") fill1DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
901	>	else if(currentHistogram.inputCollection == "genjets") fill1DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
902	>	else if(currentHistogram.inputCollection == "mcparticles") fill1DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
903	>	else if(currentHistogram.inputCollection == "primaryvertexs") fill1DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
904	>	else if(currentHistogram.inputCollection == "bxlumis") fill1DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
905	>	else if(currentHistogram.inputCollection == "photons") fill1DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
906	>	else if(currentHistogram.inputCollection == "superclusters") fill1DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
907	>	else if(currentHistogram.inputCollection == "trigobjs") fill1DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
908	>	else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill1DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
909	>	}
910	>	else if(currentHistogram.inputVariables.size() == 2){
911	>	TH2D* histo;
912	>	histo = twoDHists_.at(currentChannelIndex).at(currentCut).at(currentHistogram.name);
913	>
914	>	if     (currentHistogram.inputCollection == "jets")            fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("jets").at(currentDir),scaleFactor);
915	>	else if(currentHistogram.inputCollection == "secondary jets")  fill2DHistogram(histo,currentHistogram,jets.product(),cumulativeFlags.at("secondary jets").at(currentDir),scaleFactor);
916	>	else if(currentHistogram.inputCollection == "muons")           fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("muons").at(currentDir),scaleFactor);
917	>	else if(currentHistogram.inputCollection == "secondary muons") fill2DHistogram(histo,currentHistogram,muons.product(),cumulativeFlags.at("secondary muons").at(currentDir),scaleFactor);
918	>	else if(currentHistogram.inputCollection == "muon-muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
919	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
920	>	cumulativeFlags.at("muon-muon pairs").at(currentDir),scaleFactor);
921	>	else if(currentHistogram.inputCollection == "muon-secondary muon pairs") fill2DHistogram(histo,currentHistogram,muons.product(),muons.product(), \
922	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("secondary muons").at(currentDir), \
923	>	cumulativeFlags.at("muon-secondary muon pairs").at(currentDir),scaleFactor);
924	>	else if(currentHistogram.inputCollection == "electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("electrons").at(currentDir),scaleFactor);
925	>	else if(currentHistogram.inputCollection == "secondary electrons") fill2DHistogram(histo,currentHistogram,electrons.product(),cumulativeFlags.at("secondary electrons").at(currentDir),scaleFactor);
926	>	else if(currentHistogram.inputCollection == "electron-electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
927	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("electrons").at(currentDir), \
928	>	cumulativeFlags.at("electron-electron pairs").at(currentDir),scaleFactor);
929	>	else if(currentHistogram.inputCollection == "jet-jet pairs") fill2DHistogram(histo,currentHistogram,jets.product(),jets.product(), \
930	>	cumulativeFlags.at("jets").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
931	>	cumulativeFlags.at("jet-jet pairs").at(currentDir),scaleFactor);
932	>	else if(currentHistogram.inputCollection == "electron-secondary electron pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),electrons.product(), \
933	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("secondary electrons").at(currentDir), \
934	>	cumulativeFlags.at("electron-secondary electron pairs").at(currentDir),scaleFactor);
935	>	else if(currentHistogram.inputCollection == "electron-muon pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),muons.product(), \
936	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("muons").at(currentDir), \
937	>	cumulativeFlags.at("electron-muon pairs").at(currentDir),scaleFactor);
938	>	else if(currentHistogram.inputCollection == "electron-jet pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),jets.product(), \
939	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
940	>	cumulativeFlags.at("electron-jet pairs").at(currentDir),scaleFactor);
941	>	else if(currentHistogram.inputCollection == "muon-jet pairs") fill2DHistogram(histo,currentHistogram,muons.product(),jets.product(), \
942	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
943	>	cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
944	>	else if(currentHistogram.inputCollection == "electron-track pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),tracks.product(),
945	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
946	>	cumulativeFlags.at("electron-track pairs").at(currentDir),scaleFactor);
947	>	else if(currentHistogram.inputCollection == "muon-track pairs") fill2DHistogram(histo,currentHistogram,muons.product(),tracks.product(),
948	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
949	>	cumulativeFlags.at("muon-track pairs").at(currentDir),scaleFactor);
950	>	else if(currentHistogram.inputCollection == "muon-tau pairs") fill2DHistogram(histo,currentHistogram,muons.product(),taus.product(),
951	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
952	>	cumulativeFlags.at("muon-tau pairs").at(currentDir),scaleFactor);
953	>	else if(currentHistogram.inputCollection == "tau-tau pairs") fill2DHistogram(histo,currentHistogram,taus.product(),taus.product(),
954	>	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("taus").at(currentDir),
955	>	cumulativeFlags.at("tau-tau pairs").at(currentDir),scaleFactor);
956	>	else if(currentHistogram.inputCollection == "tau-track pairs") fill2DHistogram(histo,currentHistogram,taus.product(),tracks.product(),
957	>	cumulativeFlags.at("taus").at(currentDir),cumulativeFlags.at("tracks").at(currentDir),
958	>	cumulativeFlags.at("tau-track pairs").at(currentDir),scaleFactor);
959	>	else if(currentHistogram.inputCollection == "electron-trigobj pairs") fill2DHistogram(histo,currentHistogram,electrons.product(),trigobjs.product(),
960	>	cumulativeFlags.at("electrons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
961	>	cumulativeFlags.at("electron-trigobj pairs").at(currentDir),scaleFactor);
962	>	else if(currentHistogram.inputCollection == "muon-trigobj pairs") fill2DHistogram(histo,currentHistogram,muons.product(),trigobjs.product(),
963	>	cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("trigobjs").at(currentDir),
964	>	cumulativeFlags.at("muon-trigobj pairs").at(currentDir),scaleFactor);
965	>	else if(currentHistogram.inputCollection == "events") fill2DHistogram(histo,currentHistogram,events.product(),cumulativeFlags.at("events").at(currentDir),scaleFactor);
966	>	else if(currentHistogram.inputCollection == "taus") fill2DHistogram(histo,currentHistogram,taus.product(),cumulativeFlags.at("taus").at(currentDir),scaleFactor);
967	>	else if(currentHistogram.inputCollection == "mets") fill2DHistogram(histo,currentHistogram,mets.product(),cumulativeFlags.at("mets").at(currentDir),scaleFactor);
968	>	else if(currentHistogram.inputCollection == "tracks") fill2DHistogram(histo,currentHistogram,tracks.product(),cumulativeFlags.at("tracks").at(currentDir),scaleFactor);
969	>	else if(currentHistogram.inputCollection == "track-event pairs") fill2DHistogram(histo,currentHistogram,tracks.product(),events.product(),
970	>	cumulativeFlags.at("tracks").at(currentDir),cumulativeFlags.at("events").at(currentDir),
971	>	cumulativeFlags.at("track-event pairs").at(currentDir),scaleFactor);
972	>	else if(currentHistogram.inputCollection == "genjets") fill2DHistogram(histo,currentHistogram,genjets.product(),cumulativeFlags.at("genjets").at(currentDir),scaleFactor);
973	>	else if(currentHistogram.inputCollection == "mcparticles") fill2DHistogram(histo,currentHistogram,mcparticles.product(),cumulativeFlags.at("mcparticles").at(currentDir),scaleFactor);
974	>	else if(currentHistogram.inputCollection == "primaryvertexs") fill2DHistogram(histo,currentHistogram,primaryvertexs.product(),cumulativeFlags.at("primaryvertexs").at(currentDir),scaleFactor);
975	>	else if(currentHistogram.inputCollection == "bxlumis") fill2DHistogram(histo,currentHistogram,bxlumis.product(),cumulativeFlags.at("bxlumis").at(currentDir),scaleFactor);
976	>	else if(currentHistogram.inputCollection == "photons") fill2DHistogram(histo,currentHistogram,photons.product(),cumulativeFlags.at("photons").at(currentDir),scaleFactor);
977	>	else if(currentHistogram.inputCollection == "superclusters") fill2DHistogram(histo,currentHistogram,superclusters.product(),cumulativeFlags.at("superclusters").at(currentDir),scaleFactor);
978	>	else if(currentHistogram.inputCollection == "trigobjs") fill2DHistogram(histo,currentHistogram,trigobjs.product(),cumulativeFlags.at("trigobjs").at(currentDir),scaleFactor);
979	>	else if(currentHistogram.inputCollection == "stops" && datasetType_ == "signalMC") fill2DHistogram(histo,currentHistogram,stops.product(),cumulativeFlags.at("stops").at(currentDir),scaleFactor);
980	>	}
981	>	}
982	>
983	>	//fills histograms with the sizes of collections
984	>
985	>	for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++){
986	>
987	>	string currentObject = objectsToPlot.at(currentObjectIndex);
988	>	string objectToPlot = "";
989	>
990	>	// Name of objectToPlot here must match the name specified in OSUAnalysis::OSUAnalysis().
991	>	if(currentObject == "muon-muon pairs")                         objectToPlot = "dimuonPairs";
992	>	else if(currentObject == "electron-electron pairs")            objectToPlot = "dielectronPairs";
993	>	else if(currentObject == "electron-muon pairs")                objectToPlot = "electronMuonPairs";
994	>	else if(currentObject == "electron-jet pairs")                 objectToPlot = "electronJetPairs";
995	>	else if(currentObject == "muon-jet pairs")                     objectToPlot = "muonJetPairs";
996	>	else if(currentObject == "jet-jet pairs")                      objectToPlot = "dijetPairs";
997	>	else if(currentObject == "secondary jets")                     objectToPlot = "secondaryJets";
998	>	else if(currentObject == "electron-track pairs")               objectToPlot = "electronTrackPairs";
999	>	else if(currentObject == "muon-track pairs")                   objectToPlot = "muonTrackPairs";
1000	>	else if(currentObject == "muon-tau pairs")                     objectToPlot = "muonTauPairs";
1001	>	else if(currentObject == "tau-tau pairs")                      objectToPlot = "ditauPairs";
1002	>	else if(currentObject == "tau-track pairs")                    objectToPlot = "tauTrackPairs";
1003	>	else if(currentObject == "track-event pairs")                  objectToPlot = "trackEventPairs";
1004	>	else if(currentObject == "muon-secondary muon pairs")          objectToPlot = "muonSecondaryMuonPairs";
1005	>	else if(currentObject == "secondary muons")                    objectToPlot = "secondaryMuons";
1006	>	else if(currentObject == "electron-secondary electron pairs")  objectToPlot = "electronSecondaryElectronPairs";
1007	>	else if(currentObject == "secondary electrons")                objectToPlot = "secondaryElectrons";
1008	>	else if(currentObject == "electron-trigobj pairs")             objectToPlot = "electronTrigobjPairs";
1009	>	else if(currentObject == "muon-trigobj pairs")                 objectToPlot = "muonTrigobjPairs";
1010	>	else objectToPlot = currentObject;
1011	>
1012	>	string tempCurrentObject = objectToPlot;
1013	>	tempCurrentObject.at(0) = toupper(tempCurrentObject.at(0));
1014	>	string histoName = "num" + tempCurrentObject;
1015	>
1016	>
1017	>	vector<bool> lastCutFlags = cumulativeFlags.at(currentObject).at(currentDir);
1018	>
1019	>	//count the number of objects passing all cuts
1020	>	int numToPlot = 0;
1021	>	for (uint currentFlag = 0; currentFlag != lastCutFlags.size(); currentFlag++){
1022	>	if(lastCutFlags.at(currentFlag)) numToPlot++;
1023	>	}
1024	>
1025	>	if(objectToPlot == "primaryvertexs"){
1026	>	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"BeforePileupCorrection")->Fill(primaryvertexs->size());
1027	>	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName+"AfterPileupCorrection")->Fill(primaryvertexs->size(),scaleFactor);
1028	>	}
1029	>	else {
1030	>	oneDHists_.at(currentChannelIndex).at(currentCut).at(histoName)->Fill(numToPlot,scaleFactor);
1031	>	}
1032	>	} // end for (uint currentObjectIndex = 0; currentObjectIndex != objectsToPlot.size(); currentObjectIndex++)
1033		}
1034	<	}
1034	>
1035		} //end loop over channel
1036
1037		masterCutFlow_->fillCutFlow(masterScaleFactor);

Diff Legend

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