diff --git a/L1Trigger/TrackFindingTracklet/test/L1TrackNtuplePlot.C b/L1Trigger/TrackFindingTracklet/test/L1TrackNtuplePlot.C index d2a6bf1be3d62..5bbfa4dd4c15e 100644 --- a/L1Trigger/TrackFindingTracklet/test/L1TrackNtuplePlot.C +++ b/L1Trigger/TrackFindingTracklet/test/L1TrackNtuplePlot.C @@ -76,6 +76,9 @@ void L1TrackNtuplePlot(TString type, SetPlotStyle(); + cout.setf(ios::fixed); + cout.precision(2); + // ---------------------------------------------------------------------------------------------------------------- // define input options @@ -1021,7 +1024,6 @@ void L1TrackNtuplePlot(TString type, // ---------------------------------------------------------------------------------------------------------------- int nevt = tree->GetEntries(); - cout << "number of events = " << nevt << endl; // ---------------------------------------------------------------------------------------------------------------- // event loop @@ -3633,11 +3635,19 @@ void L1TrackNtuplePlot(TString type, c.SaveAs(DIR + type + "_trk_pt.pdf"); } + // Sample z0 resolution at a couple of rapidity points. + float etaSample1 = h2_resVsEta_z0_68->GetXaxis()->GetBinCenter(1); + float etaSample2 = h2_resVsEta_z0_68->GetXaxis()->GetBinCenter(0.8 * nETARANGE); + float z0ResSample1 = h2_resVsEta_z0_68->GetBinContent(1); + float z0ResSample2 = h2_resVsEta_z0_68->GetBinContent(0.8 * nETARANGE); + fout->Close(); // --------------------------------------------------------------------------------------------------------- //some printouts + cout << "number of events = " << nevt << endl; + float k = (float)n_match_eta1p0; float N = (float)n_all_eta1p0; if (std::abs(N) > 0) @@ -3691,7 +3701,11 @@ void L1TrackNtuplePlot(TString type, cout << "# tracks/event (no pt cut)= " << (float)ntrk / nevt << endl; cout << "# tracks/event (pt > " << std::max(TP_minPt, 2.0f) << ") = " << (float)ntrk_pt2 / nevt << endl; cout << "# tracks/event (pt > 3.0) = " << (float)ntrk_pt3 / nevt << endl; - cout << "# tracks/event (pt > 10.0) = " << (float)ntrk_pt10 / nevt << endl; + cout << "# tracks/event (pt > 10.0) = " << (float)ntrk_pt10 / nevt << endl << endl; + + // z0 resolution + cout << "z0 resolution = " << z0ResSample1 << "cm at |eta| = " << etaSample1 << endl; + cout << "z0 resolution = " << z0ResSample2 << "cm at |eta| = " << etaSample2 << endl; } void SetPlotStyle() { diff --git a/L1Trigger/TrackFindingTracklet/test/makeHists.csh b/L1Trigger/TrackFindingTracklet/test/makeHists.csh index d2be4335aab4e..69a5355341e0f 100755 --- a/L1Trigger/TrackFindingTracklet/test/makeHists.csh +++ b/L1Trigger/TrackFindingTracklet/test/makeHists.csh @@ -43,7 +43,7 @@ echo "MVA track quality Histograms written to MVA_plots/" # Run track performance plotting macro set plotMacro = $CMSSW_BASE/src/L1Trigger/TrackFindingTracklet/test/L1TrackNtuplePlot.C if (-e TrkPlots) rm -r TrkPlots -\root -b -q ${plotMacro}'("'${inputFileStem}'","'${dirName}'")' | tail -n 19 >! results.out +\root -b -q ${plotMacro}'("'${inputFileStem}'","'${dirName}'")' | tail -n 24 >! results.out cat results.out echo "Tracking performance summary written to results.out" echo "Track performance histograms written to TrkPlots/"