/** * Callable by scheduler. * * \param file Pointer to the std::string with the filename. */ void Console::restoreState(void * file) { std::string & f = *(std::string*)file; std::ifstream in(f.c_str()); std::string graphed_var; while(in >> graphed_var) { if (!s_variable_watcher.isGraphed(graphed_var)) { plotVariable(std::vector<std::string>(1,graphed_var)); } } }
void plotAllVariables(){ std::vector<std::string> categories; categories.push_back(std::string("All")); categories.push_back(std::string("NoEleMatch")); categories.push_back(std::string("woG")); categories.push_back(std::string("wGwoGSF")); categories.push_back(std::string("wGwGSFwoPFMVA")); categories.push_back(std::string("wGwGSFwPFMVA")); // categories.push_back(std::string("All")); // categories.push_back(std::string("TauNoGammas")); // categories.push_back(std::string("TauHasGammasNoGsfTrack")); // categories.push_back(std::string("TauHasGammasHasGsfTrackPFmvaBelow01")); // categories.push_back(std::string("TauHasGammasHasGsfTrackPFmvaOver01")); std::vector<std::string> variables; variables.push_back(std::string("Elec_EeOverPout")); variables.push_back(std::string("Elec_EgammaOverPdif")); variables.push_back(std::string("Elec_EtotOverPin")); variables.push_back(std::string("Elec_PFMvaOutput")); variables.push_back(std::string("Elec_Logsihih")); variables.push_back(std::string("Elec_DeltaEta")); variables.push_back(std::string("Elec_HoHplusE")); variables.push_back(std::string("Elec_GSFTrackResol")); variables.push_back(std::string("Elec_GSFTracklnPt")); variables.push_back(std::string("Elec_GSFTrackEta")); // variables.push_back(std::string("Tau_AbsEta")); variables.push_back(std::string("Tau_Eta")); variables.push_back(std::string("Tau_EtaAtEcalEntrance")); variables.push_back(std::string("Tau_EmFraction")); variables.push_back(std::string("Tau_HasGsf")); variables.push_back(std::string("Tau_HadrEoP")); variables.push_back(std::string("Tau_HadrHoP")); variables.push_back(std::string("Tau_NumChargedCands")); variables.push_back(std::string("Tau_NumGammaCands")); variables.push_back(std::string("Tau_HadrMva")); variables.push_back(std::string("Tau_GammaEnFrac")); variables.push_back(std::string("Tau_GammaEtaMom")); variables.push_back(std::string("Tau_GammaPhiMom")); variables.push_back(std::string("Tau_Pt")); variables.push_back(std::string("Tau_LeadHadronPt")); variables.push_back(std::string("Tau_VisMass")); variables.push_back(std::string("Elec_AbsEta")); variables.push_back(std::string("Elec_Pt")); variables.push_back(std::string("Elec_Fbrem")); variables.push_back(std::string("Elec_Chi2KF")); variables.push_back(std::string("Elec_Chi2GSF")); variables.push_back(std::string("Elec_EarlyBrem")); variables.push_back(std::string("Elec_LateBrem")); variables.push_back(std::string("Elec_NumHits")); // variables.push_back(std::string("Tau_GsfEleMatch")); std::map<std::string, std::string> xAxisTitles; xAxisTitles["Elec_AbsEta"] = "|#eta|(e)" ; xAxisTitles["Elec_Fbrem"] = "Fbrem"; xAxisTitles["Elec_Chi2KF"] = "Chi2KF"; xAxisTitles["Elec_EarlyBrem"] = "EarlyBrem"; xAxisTitles["Elec_LateBrem"] = "LateBrem"; xAxisTitles["Elec_EeOverPout"] = "E_{e}/P_{out}"; xAxisTitles["Elec_EgammaOverPdif"] = "E_{#gamma}/(P_{in}-P_{out})"; xAxisTitles["Elec_EtotOverPin"] = "E_{tot}/P_{in}"; xAxisTitles["Elec_NumHits"] = "NHits"; xAxisTitles["Tau_AbsEta"] = "|#eta|(#tau)"; xAxisTitles["Tau_Eta"] = "#eta(#tau)"; xAxisTitles["Tau_EtaAtEcalEntrance"]= "#eta(#tau at ECAL)"; xAxisTitles["Tau_GammaEnFrac"] = "GammaEnFrac(#tau)"; xAxisTitles["Tau_EmFraction"] = "EmFraction(#tau)"; xAxisTitles["Tau_HasGsf"] = "TauHasGsf"; xAxisTitles["Tau_HadrEoP"] = "TauHadrEoP"; xAxisTitles["Tau_HadrHoP"] = "TauHadrHoP"; xAxisTitles["Tau_NumChargedCands"] = "TauNumChargedHadrCands"; xAxisTitles["Tau_NumGammaCands"] = "TauNumGammaCands"; xAxisTitles["Tau_VisMass"] = "VisMass(#tau)"; xAxisTitles["Tau_GammaEtaMom"] = "Gamma #Delta#eta RMS(#tau)"; xAxisTitles["Tau_GammaPhiMom"] = "Gamma #Delta#phi RMS(#tau)"; xAxisTitles["Elec_Pt"] = "P_{T}(e)" ; xAxisTitles["Elec_PFMvaOutput"] = "PFmva(e)" ; xAxisTitles["Elec_Logsihih"] = "Logsihih(e)" ; xAxisTitles["Elec_DeltaEta"] = "#Delta #eta (SC-Track)" ; xAxisTitles["Elec_HoHplusE"] = "H/H+E" ; xAxisTitles["Elec_Chi2GSF"] = "Chi2GSF" ; xAxisTitles["Elec_GSFTrackResol"] = "GSFTrack Resol" ; xAxisTitles["Elec_GSFTracklnPt"] = "lnP_{T}(GSFTr)" ; xAxisTitles["Elec_GSFTrackEta"] = "#eta(GSFTr)" ; xAxisTitles["Tau_Pt"] = "P_{T}(#tau)"; xAxisTitles["Tau_LeadHadronPt"] = "P_{T}(Lead hadron #tau)"; xAxisTitles["Tau_HadrMva"] = "HadrMva(#tau)"; std::map<std::string, float> xMinValues; xMinValues["Elec_AbsEta"] = 0.; xMinValues["Elec_Fbrem"] = -0.2; xMinValues["Elec_Chi2KF"] = 0.; xMinValues["Elec_EarlyBrem"] = -2; xMinValues["Elec_LateBrem"] = -2; xMinValues["Elec_EeOverPout"] = 0.; xMinValues["Elec_EgammaOverPdif"] = 0.; xMinValues["Elec_EtotOverPin"] = 0.; xMinValues["Elec_NumHits"] = 0.; xMinValues["Tau_AbsEta"] = 0.; xMinValues["Tau_Eta"] = -3.; xMinValues["Tau_EtaAtEcalEntrance"] = -3.; xMinValues["Tau_GammaEnFrac"] = 0.; xMinValues["Tau_EmFraction"] = 0.; xMinValues["Tau_HasGsf"] = 0.; xMinValues["Tau_HadrEoP"] = 0.; xMinValues["Tau_HadrHoP"] = 0.; xMinValues["Tau_NumChargedCands"] = 0.; xMinValues["Tau_NumGammaCands"] = 0.; xMinValues["Tau_VisMass"] = 0.; xMinValues["Tau_GammaEtaMom"] = 0.; xMinValues["Tau_GammaPhiMom"] = 0.; xMinValues["Elec_Pt"] = 10.; xMinValues["Elec_PFMvaOutput"] = -1.; xMinValues["Elec_Logsihih"] = -13.; xMinValues["Elec_DeltaEta"] = 0.; xMinValues["Elec_HoHplusE"] = 0.; xMinValues["Elec_Chi2GSF"] = 0.; xMinValues["Elec_GSFTrackResol"] = 0.; xMinValues["Elec_GSFTracklnPt"] = 0.; xMinValues["Elec_GSFTrackEta"] = -3.; xMinValues["Tau_Pt"] = 10.; xMinValues["Tau_LeadHadronPt"] = 10.; xMinValues["Tau_HadrMva"] = -1.; std::map<std::string, float> xMaxValues; xMaxValues["Elec_AbsEta"] = 3; xMaxValues["Elec_Fbrem"] = 1; xMaxValues["Elec_Chi2KF"] = 5; xMaxValues["Elec_EarlyBrem"] = 2; xMaxValues["Elec_LateBrem"] = 2; xMaxValues["Elec_EeOverPout"] = 4.0; xMaxValues["Elec_EgammaOverPdif"] = 4; xMaxValues["Elec_EtotOverPin"] = 4; xMaxValues["Elec_NumHits"] = 30; xMaxValues["Tau_AbsEta"] = 3; xMaxValues["Tau_Eta"] = 3; xMaxValues["Tau_EtaAtEcalEntrance"] = 3; xMaxValues["Tau_GammaEnFrac"] = 1; xMaxValues["Tau_EmFraction"] = 1; xMaxValues["Tau_HasGsf"] = 2; xMaxValues["Tau_HadrEoP"] = 1; xMaxValues["Tau_HadrHoP"] = 1; xMaxValues["Tau_NumChargedCands"] = 5; xMaxValues["Tau_NumGammaCands"] = 5; xMaxValues["Tau_VisMass"] = 1.8; xMaxValues["Tau_GammaEtaMom"] = 3.; xMaxValues["Tau_GammaPhiMom"] = 5.; xMaxValues["Elec_Pt"] = 80.; xMaxValues["Elec_PFMvaOutput"] = 1.; xMaxValues["Elec_Logsihih"] = -2.; xMaxValues["Elec_DeltaEta"] = 0.05; xMaxValues["Elec_HoHplusE"] = 1.; xMaxValues["Elec_Chi2GSF"] = 5.; xMaxValues["Elec_GSFTrackResol"] = 1.; xMaxValues["Elec_GSFTracklnPt"] = 15.; xMaxValues["Elec_GSFTrackEta"] = 3.; xMaxValues["Tau_Pt"] = 150.; xMaxValues["Tau_LeadHadronPt"] = 150.; xMaxValues["Tau_HadrMva"] = 1.; std::map<std::string, int> nBins; nBins["Elec_AbsEta"] = 100; nBins["Elec_Fbrem"] = 100; nBins["Elec_Chi2KF"] = 50; nBins["Elec_EarlyBrem"] = 5; nBins["Elec_LateBrem"] = 5; nBins["Elec_EeOverPout"] = 100; nBins["Elec_EgammaOverPdif"] = 100; nBins["Elec_EtotOverPin"] = 100; nBins["Elec_NumHits"] = 30; nBins["Tau_AbsEta"] = 100; nBins["Tau_Eta"] = 100; nBins["Tau_EtaAtEcalEntrance"] = 100; nBins["Tau_GammaEnFrac"] = 20; nBins["Tau_EmFraction"] = 100; nBins["Tau_HasGsf"] = 2; nBins["Tau_HadrEoP"] = 20; nBins["Tau_HadrHoP"] = 20; nBins["Tau_NumChargedCands"] = 5; nBins["Tau_NumGammaCands"] = 5; nBins["Tau_VisMass"] = 34; nBins["Tau_GammaEtaMom"] = 30.; nBins["Tau_GammaPhiMom"] = 50.; nBins["Elec_Pt"] = 100.; nBins["Elec_PFMvaOutput"] = 100.; nBins["Elec_Logsihih"] = 50.; nBins["Elec_DeltaEta"] = 100.; nBins["Elec_HoHplusE"] = 100.; nBins["Elec_Chi2GSF"] = 50.; nBins["Elec_GSFTrackResol"] = 100.; nBins["Elec_GSFTracklnPt"] = 100.; nBins["Elec_GSFTrackEta"] = 100.; nBins["Tau_Pt"] = 100.; nBins["Tau_LeadHadronPt"] = 100.; nBins["Tau_HadrMva"] = 50.; for ( std::vector<std::string>::const_iterator category = categories.begin(); category != categories.end(); ++category ) { for ( std::vector<std::string>::const_iterator variable = variables.begin(); variable!= variables.end(); ++variable ) { plotVariable(variable->data(), category->data(), xAxisTitles[*variable].data(),"a.u.",xMinValues[*variable], xMaxValues[*variable], nBins[*variable], 0, 50, 15, 150,"Barrel" ); plotVariable(variable->data(), category->data(), xAxisTitles[*variable].data(),"a.u.",xMinValues[*variable], xMaxValues[*variable], nBins[*variable], 0, 50, 15, 150,"Endcap" ); } } }