Commit in hps-java/src/main/java/org/lcsim/hps/users/phansson on MAIN | |||
MultScatAnalysis.java | +1324 | -1095 | 1.3 -> 1.4 |
Cleaned up and split rate plots to separate driver
diff -u -r1.3 -r1.4 --- MultScatAnalysis.java 11 Jun 2012 23:00:48 -0000 1.3 +++ MultScatAnalysis.java 21 Jun 2012 15:35:44 -0000 1.4 @@ -56,12 +56,15 @@
private static int crystalCols; private static int crystalRows; int ecalClusterSel;
- boolean hide = false; -
+ boolean hide = true; +
private String outputPlotFileName = "test.aida"; private String trackCollectionName = "MatchedTracks"; private String ecalClusterCollectionName = "EcalClusters"; private String inputEcalHitCollection = "EcalCalHits";
+ private boolean doTriggerPart = false; + private boolean doTrackingPart = false; +
EcalTrackMatch trkMatchTool;
@@ -101,530 +104,688 @@
ecalClusterSel = 1;
- - plotterFrameTrig = new AIDAFrame(); - plotterFrameTrig.setTitle("Trigger"); - - IPlotter plotter_trig_other = af.createPlotterFactory().create(); - plotter_trig_other.createRegions(3,1,0); - plotter_trig_other.setTitle("Other"); - plotterFrameTrig.addPlotter(plotter_trig_other); - - - IHistogram ht1 = aida.histogram1D("trigger_count" , 7, 0, 7); - IHistogram ht2 = aida.histogram1D("trigger_bank" , 2, 0, 2); - IHistogram ht3 = aida.histogram1D("trigger_list" , 2, 0, 2); - plotter_trig_other.region(0).plot(ht1); - plotter_trig_other.region(1).plot(ht2); - plotter_trig_other.region(2).plot(ht3); - for(int i=0;i<3;++i) { - IPlotterStyle style = plotter_trig_other.region(i).style(); - //style.setParameter("hist2DStyle", "colorMap"); - //style.dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow"); - ((PlotterRegion) plotter_trig_other.region(i)).getPlot().setAllowUserInteraction(true); - ((PlotterRegion) plotter_trig_other.region(i)).getPlot().setAllowPopupMenus(true); - }
+ if(doTriggerPart) { + + plotterFrameTrig = new AIDAFrame(); + plotterFrameTrig.setTitle("Trigger"); + + IPlotter plotter_trig_other = af.createPlotterFactory().create(); + plotter_trig_other.createRegions(3,1,0); + plotter_trig_other.setTitle("Other"); + plotterFrameTrig.addPlotter(plotter_trig_other); + + + IHistogram ht1 = aida.histogram1D("trigger_count" , 7, 0, 7); + IHistogram ht2 = aida.histogram1D("trigger_bank" , 2, 0, 2); + IHistogram ht3 = aida.histogram1D("trigger_list" , 2, 0, 2); + plotter_trig_other.region(0).plot(ht1); + plotter_trig_other.region(1).plot(ht2); + plotter_trig_other.region(2).plot(ht3); + for(int i=0;i<3;++i) { + IPlotterStyle style = plotter_trig_other.region(i).style(); + //style.setParameter("hist2DStyle", "colorMap"); + //style.dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow"); + ((PlotterRegion) plotter_trig_other.region(i)).getPlot().setAllowUserInteraction(true); + ((PlotterRegion) plotter_trig_other.region(i)).getPlot().setAllowPopupMenus(true); + } + + + + + plotter_trig_tag = af.createPlotterFactory().create(); + plotter_trig_tag.createRegions(3,3,0); + plotter_trig_tag.setTitle("Tag&Probe Top"); + plotterFrameTrig.addPlotter(plotter_trig_tag); + IHistogram htag1 = aida.histogram1D("toptrig_cl_ecal_n_top" , 7, 0, 7); + IHistogram htag2 = aida.histogram1D("toptrig_cl_ecal_e_top" , 100, 0, 1500); + IHistogram htag3 = aida.histogram1D("toptrig_cl_ecal_emax_top" , 100, 0, 1500); + IHistogram htag4 = aida.histogram1D("toptrig_cl_ecal_e_bottom" , 100, 0, 1500); + IHistogram htag5 = aida.histogram1D("toptrig_cl_ecal_e_bottom_trig" , 100, 0, 1500); + IHistogram htag6 = aida.histogram1D("toptrigtag_cl_ecal_e_bottom" , 100, 0, 1500); + IHistogram htag7 = aida.histogram1D("toptrigtag_cl_ecal_e_bottom_trig" , 100, 0, 1500); + IHistogram htag8 = aida.histogram1D("toptrig_cl_ecal_e_bottom_trigeff" , 100, 0, 1500); + IHistogram htag9 = aida.histogram1D("toptrigtag_cl_ecal_e_bottom_trigeff" , 100, 0, 1500); + + plotter_trig_tag.region(0).plot(htag1); + plotter_trig_tag.region(3).plot(htag2); + plotter_trig_tag.region(6).plot(htag3); + plotter_trig_tag.region(1).plot(htag4); + plotter_trig_tag.region(2).plot(htag6); + plotter_trig_tag.region(4).plot(htag5); + plotter_trig_tag.region(5).plot(htag7); + plotter_trig_tag.region(7).plot(htag9); + plotter_trig_tag.region(8).plot(htag9); + + for(int i=0;i<9;++i) { + IPlotterStyle style = plotter_trig_tag.region(i).style(); + //style.setParameter("hist2DStyle", "colorMap"); + //style.dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow"); + ((PlotterRegion) plotter_trig_tag.region(i)).getPlot().setAllowUserInteraction(true); + ((PlotterRegion) plotter_trig_tag.region(i)).getPlot().setAllowPopupMenus(true); + } + + + IPlotter plotter_trig_tag2 = af.createPlotterFactory().create(); + plotter_trig_tag2.setTitle("Tag&Probe Top Counts"); + plotterFrameTrig.addPlotter(plotter_trig_tag2); + IHistogram htag10 = aida.histogram1D("toptrig_cl_ecal_n_bottom" , 7, 0, 7); + plotter_trig_tag2.createRegion().plot(htag10); + + if(!hide) { + plotterFrameTrig.pack(); + plotterFrameTrig.setVisible(true); + }
- - - - plotter_trig_tag = af.createPlotterFactory().create(); - plotter_trig_tag.createRegions(3,3,0); - plotter_trig_tag.setTitle("Tag&Probe Top"); - plotterFrameTrig.addPlotter(plotter_trig_tag); - IHistogram htag1 = aida.histogram1D("toptrig_cl_ecal_n_top" , 7, 0, 7); - IHistogram htag2 = aida.histogram1D("toptrig_cl_ecal_e_top" , 100, 0, 1500); - IHistogram htag3 = aida.histogram1D("toptrig_cl_ecal_emax_top" , 100, 0, 1500); - IHistogram htag4 = aida.histogram1D("toptrig_cl_ecal_e_bottom" , 100, 0, 1500); - IHistogram htag5 = aida.histogram1D("toptrig_cl_ecal_e_bottom_trig" , 100, 0, 1500); - IHistogram htag6 = aida.histogram1D("toptrigtag_cl_ecal_e_bottom" , 100, 0, 1500); - IHistogram htag7 = aida.histogram1D("toptrigtag_cl_ecal_e_bottom_trig" , 100, 0, 1500); - IHistogram htag8 = aida.histogram1D("toptrig_cl_ecal_e_bottom_trigeff" , 100, 0, 1500); - IHistogram htag9 = aida.histogram1D("toptrigtag_cl_ecal_e_bottom_trigeff" , 100, 0, 1500); - - plotter_trig_tag.region(0).plot(htag1); - plotter_trig_tag.region(3).plot(htag2); - plotter_trig_tag.region(6).plot(htag3); - plotter_trig_tag.region(1).plot(htag4); - plotter_trig_tag.region(2).plot(htag6); - plotter_trig_tag.region(4).plot(htag5); - plotter_trig_tag.region(5).plot(htag7); - plotter_trig_tag.region(7).plot(htag9); - plotter_trig_tag.region(8).plot(htag9); - - for(int i=0;i<9;++i) { - IPlotterStyle style = plotter_trig_tag.region(i).style(); - //style.setParameter("hist2DStyle", "colorMap"); - //style.dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow"); - ((PlotterRegion) plotter_trig_tag.region(i)).getPlot().setAllowUserInteraction(true); - ((PlotterRegion) plotter_trig_tag.region(i)).getPlot().setAllowPopupMenus(true); - } - - - IPlotter plotter_trig_tag2 = af.createPlotterFactory().create(); - plotter_trig_tag2.setTitle("Tag&Probe Top Counts"); - plotterFrameTrig.addPlotter(plotter_trig_tag2); - IHistogram htag10 = aida.histogram1D("toptrig_cl_ecal_n_bottom" , 7, 0, 7); - plotter_trig_tag2.createRegion().plot(htag10); - - if(!hide) { - plotterFrameTrig.pack(); - plotterFrameTrig.setVisible(true);
}
- -
plotterFrame = new AIDAFrame(); plotterFrame.setTitle("General");
- - IPlotter plotter_track_mom_pos = af.createPlotterFactory().create(); - plotter_track_mom_pos.createRegions(3,3,0); - plotter_track_mom_pos.setTitle("Track Momentum vs Position"); - plotter_track_mom_pos.style().statisticsBoxStyle().setVisible(false); - plotterFrame.addPlotter(plotter_track_mom_pos); - IHistogram hPzVsX_t = aida.histogram2D("Top track Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsX_b = aida.histogram2D("Bottom track Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsX_a = aida.histogram2D("Track Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsXqp_t = aida.histogram2D("Top track q>0 Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsXqp_b = aida.histogram2D("Bottom track q>0 Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsXqp_a = aida.histogram2D("Track q>0 Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsXqm_t = aida.histogram2D("Top track q<0 Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsXqm_b = aida.histogram2D("Bottom track q<0 Pz vs X" , 25, -500,500, 50, 0, 3500); - IHistogram hPzVsXqm_a = aida.histogram2D("Track q<0 Pz vs X" , 25, -500,500, 50, 0, 3500); - plotter_track_mom_pos.region(0).plot(hPzVsX_t); - plotter_track_mom_pos.region(1).plot(hPzVsX_b); - plotter_track_mom_pos.region(2).plot(hPzVsX_a); - plotter_track_mom_pos.region(3).plot(hPzVsXqp_t); - plotter_track_mom_pos.region(4).plot(hPzVsXqp_b); - plotter_track_mom_pos.region(5).plot(hPzVsXqp_a); - plotter_track_mom_pos.region(6).plot(hPzVsXqm_t); - plotter_track_mom_pos.region(7).plot(hPzVsXqm_b); - plotter_track_mom_pos.region(8).plot(hPzVsXqm_a); - for(int i=0;i<9;++i) { - plotter_track_mom_pos.region(i).style().setParameter("hist2DStyle", "colorMap"); - plotter_track_mom_pos.region(i).style().dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow"); - ((PlotterRegion) plotter_track_mom_pos.region(i)).getPlot().setAllowUserInteraction(true); - ((PlotterRegion) plotter_track_mom_pos.region(i)).getPlot().setAllowPopupMenus(true); - - } - - - IPlotter plotter_track_mom = af.createPlotterFactory().create(); - plotter_track_mom.createRegions(3,3,0); - plotter_track_mom.setTitle("Track Momentum"); - plotterFrame.addPlotter(plotter_track_mom); - IHistogram hPz_t = aida.histogram1D("Top track Pz" , 50, 0, 3500); - IHistogram hPz_b = aida.histogram1D("Bottom track Pz" , 50, 0, 3500); - IHistogram hPz_a = aida.histogram1D("Track Pz" , 50, 0, 3500); - IHistogram hmPz_t = aida.histogram1D("Matched top track Pz" , 50, 0, 3500); - IHistogram hmPz_b = aida.histogram1D("Matched bottom track Pz" , 50, 0, 3500); - IHistogram hmPz_a = aida.histogram1D("Matched track Pz" , 50, 0, 3500); - IHistogram hmsPz_t = aida.histogram1D("Matched sel top track Pz" , 50, 0, 3500); - IHistogram hmsPz_b = aida.histogram1D("Matched sel bottom track Pz" , 50, 0, 3500); - IHistogram hmsPz_a = aida.histogram1D("Matched sel track Pz" , 50, 0, 3500); - - - plotter_track_mom.region(0).plot(hPz_t); - plotter_track_mom.region(1).plot(hPz_b); - plotter_track_mom.region(2).plot(hPz_a); - plotter_track_mom.region(3).plot(hmPz_t); - plotter_track_mom.region(4).plot(hmPz_b); - plotter_track_mom.region(5).plot(hmPz_a); - plotter_track_mom.region(6).plot(hmsPz_t); - plotter_track_mom.region(7).plot(hmsPz_b); - plotter_track_mom.region(8).plot(hmsPz_a); - - IPlotter plotter_track_mom2 = af.createPlotterFactory().create(); - plotter_track_mom2.createRegions(3,3,0); - plotter_track_mom2.setTitle("Track Momentum Target Outliers X"); - plotterFrame.addPlotter(plotter_track_mom2); - IHistogram hPzTX_t = aida.histogram1D("Top track !target in X Pz" , 50, 0, 3500); - IHistogram hPzTX_b = aida.histogram1D("Bottom track !target in X Pz" , 50, 0, 3500); - IHistogram hPzTX_a = aida.histogram1D("Track !target in X Pz" , 50, 0, 3500); - IHistogram hmPzTX_t = aida.histogram1D("Matched top track !target in X Pz" , 50, 0, 3500); - IHistogram hmPzTX_b = aida.histogram1D("Matched bottom track !target in X Pz" , 50, 0, 3500); - IHistogram hmPzTX_a = aida.histogram1D("Matched track !target in X Pz" , 50, 0, 3500); - IHistogram hmsPzTX_t = aida.histogram1D("Matched sel top track !target in X Pz" , 50, 0, 3500); - IHistogram hmsPzTX_b = aida.histogram1D("Matched sel bottom track !target in X Pz" , 50, 0, 3500); - IHistogram hmsPzTX_a = aida.histogram1D("Matched sel track !target in X Pz" , 50, 0, 3500); - - - plotter_track_mom2.region(0).plot(hPzTX_t); - plotter_track_mom2.region(1).plot(hPzTX_b); - plotter_track_mom2.region(2).plot(hPzTX_a); - plotter_track_mom2.region(3).plot(hmPzTX_t); - plotter_track_mom2.region(4).plot(hmPzTX_b); - plotter_track_mom2.region(5).plot(hmPzTX_a); - plotter_track_mom2.region(6).plot(hmsPzTX_t); - plotter_track_mom2.region(7).plot(hmsPzTX_b); - plotter_track_mom2.region(8).plot(hmsPzTX_a); - - - IPlotter plotter_track_ext = af.createPlotterFactory().create(); - plotter_track_ext.createRegions(2,3,0); - plotter_track_ext.setTitle("Track @ Target"); - plotterFrame.addPlotter(plotter_track_ext); - IHistogram hTrkXAtConv_t = aida.histogram1D("Top track X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConv_t = aida.histogram1D("Top track Y @ -67cm" , 50, -20, 20); - IHistogram hTrkXAtConv_b = aida.histogram1D("Bottom track X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConv_b = aida.histogram1D("Bottom track Y @ -67cm" , 50, -20, 20); - IHistogram hTrkXAtConv_a = aida.histogram1D("Track X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConv_a = aida.histogram1D("Track Y @ -67cm" , 50, -20, 20); - plotter_track_ext.region(0).plot(hTrkXAtConv_t); - plotter_track_ext.region(1).plot(hTrkXAtConv_b); - plotter_track_ext.region(2).plot(hTrkXAtConv_a); - plotter_track_ext.region(3).plot(hTrkYAtConv_t); - plotter_track_ext.region(4).plot(hTrkYAtConv_b); - plotter_track_ext.region(5).plot(hTrkYAtConv_a); - - IHistogram hTrkHMXAtConv_t = aida.histogram1D("Top track Pz>1GeV X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConv_t = aida.histogram1D("Top track Pz>1GeV Y @ -67cm" , 50, -20, 20); - IHistogram hTrkHMXAtConv_b = aida.histogram1D("Bottom track Pz>1GeV X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConv_b = aida.histogram1D("Bottom track Pz>1GeV Y @ -67cm" , 50, -20, 20); - IHistogram hTrkHMXAtConv_a = aida.histogram1D("Track Pz>1GeV X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConv_a = aida.histogram1D("Track Pz>1GeV Y @ -67cm" , 50, -20, 20); - plotter_track_ext.region(0).plot(hTrkHMXAtConv_t,"mode=overlay"); - plotter_track_ext.region(1).plot(hTrkHMXAtConv_b,"mode=overlay"); - plotter_track_ext.region(2).plot(hTrkHMXAtConv_a,"mode=overlay"); - plotter_track_ext.region(3).plot(hTrkHMYAtConv_t,"mode=overlay"); - plotter_track_ext.region(4).plot(hTrkHMYAtConv_b,"mode=overlay"); - plotter_track_ext.region(5).plot(hTrkHMYAtConv_a,"mode=overlay"); - - IPlotter plotter_track_ext_entr = af.createPlotterFactory().create(); - plotter_track_ext_entr.createRegions(2,3,0); - plotter_track_ext_entr.setTitle("Track @ 0cm"); - plotterFrame.addPlotter(plotter_track_ext_entr); - IHistogram hTrkXAtEntr_t = aida.histogram1D("Top track X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntr_t = aida.histogram1D("Top track Y @ 0cm" , 50, -40, 40); - IHistogram hTrkXAtEntr_b = aida.histogram1D("Bottom track X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntr_b = aida.histogram1D("Bottom track Y @ 0cm" , 50, -40, 40); - IHistogram hTrkXAtEntr_a = aida.histogram1D("Track X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntr_a = aida.histogram1D("Track Y @ 0cm" , 50, -40, 40); - plotter_track_ext_entr.region(0).plot(hTrkXAtEntr_t); - plotter_track_ext_entr.region(1).plot(hTrkXAtEntr_b); - plotter_track_ext_entr.region(2).plot(hTrkXAtEntr_a); - plotter_track_ext_entr.region(3).plot(hTrkYAtEntr_t); - plotter_track_ext_entr.region(4).plot(hTrkYAtEntr_b); - plotter_track_ext_entr.region(5).plot(hTrkYAtEntr_a); - - - IHistogram hTrkHMXAtEntr_t = aida.histogram1D("Top track Pz>1GeV X @ 0cm" , 50, -100, 100); - IHistogram hTrkHMYAtEntr_t = aida.histogram1D("Top track Pz>1GeV Y @ 0cm" , 50, -40, 40); - IHistogram hTrkHMXAtEntr_b = aida.histogram1D("Bottom track Pz>1GeV X @ 0cm" , 50, -100, 100); - IHistogram hTrkHMYAtEntr_b = aida.histogram1D("Bottom track Pz>1GeV Y @ 0cm" , 50, -40, 40); - IHistogram hTrkHMXAtEntr_a = aida.histogram1D("Track Pz>1GeV X @ 0cm" , 50, -100, 100); - IHistogram hTrkHMYAtEntr_a = aida.histogram1D("Track Pz>1GeV Y @ 0cm" , 50, -40, 40); - plotter_track_ext_entr.region(0).plot(hTrkHMXAtEntr_t,"mode=overlay"); - plotter_track_ext_entr.region(1).plot(hTrkHMXAtEntr_b,"mode=overlay"); - plotter_track_ext_entr.region(2).plot(hTrkHMXAtEntr_a,"mode=overlay"); - plotter_track_ext_entr.region(3).plot(hTrkHMYAtEntr_t,"mode=overlay"); - plotter_track_ext_entr.region(4).plot(hTrkHMYAtEntr_b,"mode=overlay"); - plotter_track_ext_entr.region(5).plot(hTrkHMYAtEntr_a,"mode=overlay"); - - - IPlotter plotter_track_ext_coll = af.createPlotterFactory().create(); - plotter_track_ext_coll.createRegions(2,3,0); - plotter_track_ext_coll.setTitle("Track @ -150cm"); - plotterFrame.addPlotter(plotter_track_ext_coll); - IHistogram hTrkXAtColl_t = aida.histogram1D("Top track X @ -150cm" , 50, -100, 100); - IHistogram hTrkYAtColl_t = aida.histogram1D("Top track Y @ -150cm" , 50, -40, 40); - IHistogram hTrkXAtColl_b = aida.histogram1D("Bottom track X @ -150cm" , 50, -100, 100); - IHistogram hTrkYAtColl_b = aida.histogram1D("Bottom track Y @ -150cm" , 50, -40, 40); - IHistogram hTrkXAtColl_a = aida.histogram1D("Track X @ -150cm" , 50, -100, 100); - IHistogram hTrkYAtColl_a = aida.histogram1D("Track Y @ -150cm" , 50, -40, 40); - plotter_track_ext_coll.region(0).plot(hTrkXAtColl_t); - plotter_track_ext_coll.region(1).plot(hTrkXAtColl_b); - plotter_track_ext_coll.region(2).plot(hTrkXAtColl_a); - plotter_track_ext_coll.region(3).plot(hTrkYAtColl_t); - plotter_track_ext_coll.region(4).plot(hTrkYAtColl_b); - plotter_track_ext_coll.region(5).plot(hTrkYAtColl_a); - - - IHistogram hTrkHMXAtColl_t = aida.histogram1D("Top track Pz>1GeV X @ -150cm" , 50, -100, 100); - IHistogram hTrkHMYAtColl_t = aida.histogram1D("Top track Pz>1GeV Y @ -150cm" , 50, -40, 40); - IHistogram hTrkHMXAtColl_b = aida.histogram1D("Bottom track Pz>1GeV X @ -150cm" , 50, -100, 100); - IHistogram hTrkHMYAtColl_b = aida.histogram1D("Bottom track Pz>1GeV Y @ -150cm" , 50, -40, 40); - IHistogram hTrkHMXAtColl_a = aida.histogram1D("Track Pz>1GeV X @ -150cm" , 50, -100, 100); - IHistogram hTrkHMYAtColl_a = aida.histogram1D("Track Pz>1GeV Y @ -150cm" , 50, -40, 40); - plotter_track_ext_coll.region(0).plot(hTrkHMXAtColl_t,"mode=overlay"); - plotter_track_ext_coll.region(1).plot(hTrkHMXAtColl_b,"mode=overlay"); - plotter_track_ext_coll.region(2).plot(hTrkHMXAtColl_a,"mode=overlay"); - plotter_track_ext_coll.region(3).plot(hTrkHMYAtColl_t,"mode=overlay"); - plotter_track_ext_coll.region(4).plot(hTrkHMYAtColl_b,"mode=overlay"); - plotter_track_ext_coll.region(5).plot(hTrkHMYAtColl_a,"mode=overlay"); - - -
- IPlotter plotter_track_ext2 = af.createPlotterFactory().create(); - plotter_track_ext2.createRegions(2,3,0); - plotter_track_ext2.setTitle("Matched track @ Target "); - plotterFrame.addPlotter(plotter_track_ext2); - IHistogram hmTrkXAtConv_t = aida.histogram1D("Top matched track X @ -67cm" , 50, -100, 100); - IHistogram hmTrkYAtConv_t = aida.histogram1D("Top matched track Y @ -67cm" , 50, -20, 20); - IHistogram hmTrkXAtConv_b = aida.histogram1D("Bottom matched track X @ -67cm" , 50, -100, 100); - IHistogram hmTrkYAtConv_b = aida.histogram1D("Bottom matched track Y @ -67cm" , 50, -20, 20); - IHistogram hmTrkXAtConv_a = aida.histogram1D("Matched track X @ -67cm" , 50, -100, 100); - IHistogram hmTrkYAtConv_a = aida.histogram1D("Matched track Y @ -67cm" , 50, -20, 20); - plotter_track_ext2.region(0).plot(hmTrkXAtConv_t); - plotter_track_ext2.region(1).plot(hmTrkXAtConv_b); - plotter_track_ext2.region(2).plot(hmTrkXAtConv_a); - plotter_track_ext2.region(3).plot(hmTrkYAtConv_t); - plotter_track_ext2.region(4).plot(hmTrkYAtConv_b); - plotter_track_ext2.region(5).plot(hmTrkYAtConv_a); -
- - - - IPlotter plotter_track_ext3 = af.createPlotterFactory().create(); - plotter_track_ext3.createRegions(2,3,0); - plotter_track_ext3.setTitle("Matched sel track @ Target"); - plotterFrame.addPlotter(plotter_track_ext3); - IHistogram hmsTrkXAtConv_t = aida.histogram1D("Top sel matched track X @ -67cm" , 50, -100, 100); - IHistogram hmsTrkYAtConv_t = aida.histogram1D("Top sel matched track Y @ -67cm" , 50, -100, 100); - IHistogram hmsTrkXAtConv_b = aida.histogram1D("Bottom sel matched track X @ -67cm" , 50, -100, 100); - IHistogram hmsTrkYAtConv_b = aida.histogram1D("Bottom sel matched track Y @ -67cm" , 50, -100, 100); - IHistogram hmsTrkXAtConv_a = aida.histogram1D("Matched sel track X @ -67cm" , 50, -100, 100); - IHistogram hmsTrkYAtConv_a = aida.histogram1D("Matched sel track Y @ -67cm" , 50, -100, 100); - plotter_track_ext3.region(0).plot(hmsTrkXAtConv_t); - plotter_track_ext3.region(1).plot(hmsTrkXAtConv_b); - plotter_track_ext3.region(2).plot(hmsTrkXAtConv_a); - plotter_track_ext3.region(3).plot(hmsTrkYAtConv_t); - plotter_track_ext3.region(4).plot(hmsTrkYAtConv_b); - plotter_track_ext3.region(5).plot(hmsTrkYAtConv_a); - - - - - IPlotter plotter_track_ext4 = af.createPlotterFactory().create(); - plotter_track_ext4.createRegions(2,3,0); - plotter_track_ext4.setTitle("Track @ Target Pos. Charge"); - plotterFrame.addPlotter(plotter_track_ext4); - IHistogram hTrkXAtConvqp_t = aida.histogram1D("Top track q>0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConvqp_t = aida.histogram1D("Top track q>0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkXAtConvqp_b = aida.histogram1D("Bottom track q>0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConvqp_b = aida.histogram1D("Bottom track q>0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkXAtConvqp_a = aida.histogram1D("Track q>0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConvqp_a = aida.histogram1D("Track q>0 Y @ -67cm" , 50, -20, 20); - plotter_track_ext4.region(0).plot(hTrkXAtConvqp_t); - plotter_track_ext4.region(1).plot(hTrkXAtConvqp_b); - plotter_track_ext4.region(2).plot(hTrkXAtConvqp_a); - plotter_track_ext4.region(3).plot(hTrkYAtConvqp_t); - plotter_track_ext4.region(4).plot(hTrkYAtConvqp_b); - plotter_track_ext4.region(5).plot(hTrkYAtConvqp_a); - - - IHistogram hTrkHMXAtConvqp_t = aida.histogram1D("Top track Pz>1GeV q>0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConvqp_t = aida.histogram1D("Top track Pz>1GeV q>0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkHMXAtConvqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConvqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkHMXAtConvqp_a = aida.histogram1D("Track Pz>1GeV q>0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConvqp_a = aida.histogram1D("Track Pz>1GeV q>0 Y @ -67cm" , 50, -20, 20); - plotter_track_ext4.region(0).plot(hTrkHMXAtConvqp_t,"mode=overlay"); - plotter_track_ext4.region(1).plot(hTrkHMXAtConvqp_b,"mode=overlay"); - plotter_track_ext4.region(2).plot(hTrkHMXAtConvqp_a,"mode=overlay"); - plotter_track_ext4.region(3).plot(hTrkHMYAtConvqp_t,"mode=overlay"); - plotter_track_ext4.region(4).plot(hTrkHMYAtConvqp_b,"mode=overlay"); - plotter_track_ext4.region(5).plot(hTrkHMYAtConvqp_a,"mode=overlay"); - - - - IPlotter plotter_track_ext4_entr = af.createPlotterFactory().create(); - plotter_track_ext4_entr.createRegions(2,3,0); - plotter_track_ext4_entr.setTitle("Track @ 0cm Pos. Charge"); - plotterFrame.addPlotter(plotter_track_ext4_entr); - IHistogram hTrkXAtEntrqp_t = aida.histogram1D("Top track q>0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntrqp_t = aida.histogram1D("Top track q>0 Y @ 0cm" , 50, -40, 40); - IHistogram hTrkXAtEntrqp_b = aida.histogram1D("Bottom track q>0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntrqp_b = aida.histogram1D("Bottom track q>0 Y @ 0cm" , 50, -40, 40); - IHistogram hTrkXAtEntrqp_a = aida.histogram1D("Track q>0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntrqp_a = aida.histogram1D("Track q>0 Y @ 0cm" , 50, -40, 40); - plotter_track_ext4_entr.region(0).plot(hTrkXAtEntrqp_t); - plotter_track_ext4_entr.region(1).plot(hTrkXAtEntrqp_b); - plotter_track_ext4_entr.region(2).plot(hTrkXAtEntrqp_a); - plotter_track_ext4_entr.region(3).plot(hTrkYAtEntrqp_t); - plotter_track_ext4_entr.region(4).plot(hTrkYAtEntrqp_b); - plotter_track_ext4_entr.region(5).plot(hTrkYAtEntrqp_a); - - - IHistogram hTrkHMXAtEntrqp_t = aida.histogram1D("Top track Pz>1GeV q>0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkHMYAtEntrqp_t = aida.histogram1D("Top track Pz>1GeV q>0 Y @ 0cm" , 50, -40, 40); - IHistogram hTrkHMXAtEntrqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkHMYAtEntrqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 Y @ 0cm" , 50, -40, 40); - IHistogram hTrkHMXAtEntrqp_a = aida.histogram1D("Track Pz>1GeV q>0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkHMYAtEntrqp_a = aida.histogram1D("Track Pz>1GeV q>0 Y @ 0cm" , 50, -40, 40); - plotter_track_ext4_entr.region(0).plot(hTrkHMXAtEntrqp_t,"mode=overlay"); - plotter_track_ext4_entr.region(1).plot(hTrkHMXAtEntrqp_b,"mode=overlay"); - plotter_track_ext4_entr.region(2).plot(hTrkHMXAtEntrqp_a,"mode=overlay"); - plotter_track_ext4_entr.region(3).plot(hTrkHMYAtEntrqp_t,"mode=overlay"); - plotter_track_ext4_entr.region(4).plot(hTrkHMYAtEntrqp_b,"mode=overlay"); - plotter_track_ext4_entr.region(5).plot(hTrkHMYAtEntrqp_a,"mode=overlay"); - - - - IPlotter plotter_track_ext5 = af.createPlotterFactory().create(); - plotter_track_ext5.createRegions(2,3,0); - plotter_track_ext5.setTitle("Track @ Target Neg. Charge"); - plotterFrame.addPlotter(plotter_track_ext5); - IHistogram hTrkXAtConvqm_t = aida.histogram1D("Top track q<0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConvqm_t = aida.histogram1D("Top track q<0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkXAtConvqm_b = aida.histogram1D("Bottom track q<0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConvqm_b = aida.histogram1D("Bottom track q<0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkXAtConvqm_a = aida.histogram1D("Track q<0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkYAtConvqm_a = aida.histogram1D("Track q<0 Y @ -67cm" , 50, -20, 20); - plotter_track_ext5.region(0).plot(hTrkXAtConvqm_t); - plotter_track_ext5.region(1).plot(hTrkXAtConvqm_b); - plotter_track_ext5.region(2).plot(hTrkXAtConvqm_a); - plotter_track_ext5.region(3).plot(hTrkYAtConvqm_t); - plotter_track_ext5.region(4).plot(hTrkYAtConvqm_b); - plotter_track_ext5.region(5).plot(hTrkYAtConvqm_a); - - IHistogram hTrkHMXAtConvqm_t = aida.histogram1D("Top track Pz>1GeV q<0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConvqm_t = aida.histogram1D("Top track Pz>1GeV q<0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkHMXAtConvqm_b = aida.histogram1D("Bottom track Pz>1GeV q<0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConvqm_b = aida.histogram1D("Bottom track Pz>1GeV q<0 Y @ -67cm" , 50, -20, 20); - IHistogram hTrkHMXAtConvqm_a = aida.histogram1D("Track Pz>1GeV q<0 X @ -67cm" , 50, -100, 100); - IHistogram hTrkHMYAtConvqm_a = aida.histogram1D("Track Pz>1GeV q<0 Y @ -67cm" , 50, -20, 20); - plotter_track_ext5.region(0).plot(hTrkHMXAtConvqm_t,"mode=overlay"); - plotter_track_ext5.region(1).plot(hTrkHMXAtConvqm_b,"mode=overlay"); - plotter_track_ext5.region(2).plot(hTrkHMXAtConvqm_a,"mode=overlay"); - plotter_track_ext5.region(3).plot(hTrkHMYAtConvqm_t,"mode=overlay"); - plotter_track_ext5.region(4).plot(hTrkHMYAtConvqm_b,"mode=overlay"); - plotter_track_ext5.region(5).plot(hTrkHMYAtConvqm_a,"mode=overlay"); - - - IPlotter plotter_track_ext5_entr = af.createPlotterFactory().create(); - plotter_track_ext5_entr.createRegions(2,3,0); - plotter_track_ext5_entr.setTitle("Track @ 0cm Neg. Charge"); - plotterFrame.addPlotter(plotter_track_ext5_entr); - IHistogram hTrkXAtEntrqm_t = aida.histogram1D("Top track q<0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntrqm_t = aida.histogram1D("Top track q<0 Y @ 0cm" , 50, -40, 40); - IHistogram hTrkXAtEntrqm_b = aida.histogram1D("Bottom track q<0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntrqm_b = aida.histogram1D("Bottom track q<0 Y @ 0cm" , 50, -40, 40); - IHistogram hTrkXAtEntrqm_a = aida.histogram1D("Track q<0 X @ 0cm" , 50, -100, 100); - IHistogram hTrkYAtEntrqm_a = aida.histogram1D("Track q<0 Y @ 0cm" , 50, -40, 40); - plotter_track_ext5_entr.region(0).plot(hTrkXAtEntrqm_t); - plotter_track_ext5_entr.region(1).plot(hTrkXAtEntrqm_b); - plotter_track_ext5_entr.region(2).plot(hTrkXAtEntrqm_a); - plotter_track_ext5_entr.region(3).plot(hTrkYAtEntrqm_t); - plotter_track_ext5_entr.region(4).plot(hTrkYAtEntrqm_b); - plotter_track_ext5_entr.region(5).plot(hTrkYAtEntrqm_a); -
+ IPlotter plotter_track_mom_pos = af.createPlotterFactory().create(); + plotter_track_mom_pos.createRegions(3,3,0); + plotter_track_mom_pos.setTitle("Track Momentum vs Position"); + plotter_track_mom_pos.style().statisticsBoxStyle().setVisible(false); + plotterFrame.addPlotter(plotter_track_mom_pos); + IHistogram hPzVsX_t = aida.histogram2D("Top track Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsX_b = aida.histogram2D("Bottom track Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsX_a = aida.histogram2D("Track Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsXqp_t = aida.histogram2D("Top track q>0 Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsXqp_b = aida.histogram2D("Bottom track q>0 Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsXqp_a = aida.histogram2D("Track q>0 Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsXqm_t = aida.histogram2D("Top track q<0 Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsXqm_b = aida.histogram2D("Bottom track q<0 Pz vs X" , 25, -500,500, 50, 0, 3500); + IHistogram hPzVsXqm_a = aida.histogram2D("Track q<0 Pz vs X" , 25, -500,500, 50, 0, 3500); + plotter_track_mom_pos.region(0).plot(hPzVsX_t); + plotter_track_mom_pos.region(1).plot(hPzVsX_b); + plotter_track_mom_pos.region(2).plot(hPzVsX_a); + plotter_track_mom_pos.region(3).plot(hPzVsXqp_t); + plotter_track_mom_pos.region(4).plot(hPzVsXqp_b); + plotter_track_mom_pos.region(5).plot(hPzVsXqp_a); + plotter_track_mom_pos.region(6).plot(hPzVsXqm_t); + plotter_track_mom_pos.region(7).plot(hPzVsXqm_b); + plotter_track_mom_pos.region(8).plot(hPzVsXqm_a); + for(int i=0;i<9;++i) { + plotter_track_mom_pos.region(i).style().setParameter("hist2DStyle", "colorMap"); + plotter_track_mom_pos.region(i).style().dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow"); + ((PlotterRegion) plotter_track_mom_pos.region(i)).getPlot().setAllowUserInteraction(true); + ((PlotterRegion) plotter_track_mom_pos.region(i)).getPlot().setAllowPopupMenus(true); + + } + + + IPlotter plotter_track_mom = af.createPlotterFactory().create(); + plotter_track_mom.createRegions(3,3,0); + plotter_track_mom.setTitle("Track Momentum"); + plotterFrame.addPlotter(plotter_track_mom); + IHistogram hPz_t = aida.histogram1D("Top track Pz" , 50, 0, 3500); + IHistogram hPz_b = aida.histogram1D("Bottom track Pz" , 50, 0, 3500); + IHistogram hPz_a = aida.histogram1D("Track Pz" , 50, 0, 3500); + IHistogram hmPz_t = aida.histogram1D("Matched top track Pz" , 50, 0, 3500); + IHistogram hmPz_b = aida.histogram1D("Matched bottom track Pz" , 50, 0, 3500); + IHistogram hmPz_a = aida.histogram1D("Matched track Pz" , 50, 0, 3500); + IHistogram hmsPz_t = aida.histogram1D("Matched sel top track Pz" , 50, 0, 3500); + IHistogram hmsPz_b = aida.histogram1D("Matched sel bottom track Pz" , 50, 0, 3500); + IHistogram hmsPz_a = aida.histogram1D("Matched sel track Pz" , 50, 0, 3500); + + + plotter_track_mom.region(0).plot(hPz_t); + plotter_track_mom.region(1).plot(hPz_b); + plotter_track_mom.region(2).plot(hPz_a); + plotter_track_mom.region(3).plot(hmPz_t); + plotter_track_mom.region(4).plot(hmPz_b); + plotter_track_mom.region(5).plot(hmPz_a); + plotter_track_mom.region(6).plot(hmsPz_t); + plotter_track_mom.region(7).plot(hmsPz_b); + plotter_track_mom.region(8).plot(hmsPz_a); + + IPlotter plotter_track_mom2 = af.createPlotterFactory().create(); + plotter_track_mom2.createRegions(3,3,0); + plotter_track_mom2.setTitle("Track Momentum Target Outliers X"); + plotterFrame.addPlotter(plotter_track_mom2); + IHistogram hPzTX_t = aida.histogram1D("Top track !target in X Pz" , 50, 0, 3500); + IHistogram hPzTX_b = aida.histogram1D("Bottom track !target in X Pz" , 50, 0, 3500); + IHistogram hPzTX_a = aida.histogram1D("Track !target in X Pz" , 50, 0, 3500); + IHistogram hmPzTX_t = aida.histogram1D("Matched top track !target in X Pz" , 50, 0, 3500); + IHistogram hmPzTX_b = aida.histogram1D("Matched bottom track !target in X Pz" , 50, 0, 3500); + IHistogram hmPzTX_a = aida.histogram1D("Matched track !target in X Pz" , 50, 0, 3500); + IHistogram hmsPzTX_t = aida.histogram1D("Matched sel top track !target in X Pz" , 50, 0, 3500); + IHistogram hmsPzTX_b = aida.histogram1D("Matched sel bottom track !target in X Pz" , 50, 0, 3500); + IHistogram hmsPzTX_a = aida.histogram1D("Matched sel track !target in X Pz" , 50, 0, 3500); + + + plotter_track_mom2.region(0).plot(hPzTX_t); + plotter_track_mom2.region(1).plot(hPzTX_b); + plotter_track_mom2.region(2).plot(hPzTX_a); + plotter_track_mom2.region(3).plot(hmPzTX_t); + plotter_track_mom2.region(4).plot(hmPzTX_b); + plotter_track_mom2.region(5).plot(hmPzTX_a); + plotter_track_mom2.region(6).plot(hmsPzTX_t); + plotter_track_mom2.region(7).plot(hmsPzTX_b); + plotter_track_mom2.region(8).plot(hmsPzTX_a); + + + IPlotter plotter_track_ext = af.createPlotterFactory().create(); + plotter_track_ext.createRegions(2,3,0); + plotter_track_ext.setTitle("Track @ Target"); + plotterFrame.addPlotter(plotter_track_ext); + IHistogram hTrkXAtConv_t = aida.histogram1D("Top track X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConv_t = aida.histogram1D("Top track Y @ -67cm" , 50, -20, 20); + IHistogram hTrkXAtConv_b = aida.histogram1D("Bottom track X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConv_b = aida.histogram1D("Bottom track Y @ -67cm" , 50, -20, 20); + IHistogram hTrkXAtConv_a = aida.histogram1D("Track X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConv_a = aida.histogram1D("Track Y @ -67cm" , 50, -20, 20); + plotter_track_ext.region(0).plot(hTrkXAtConv_t); + plotter_track_ext.region(1).plot(hTrkXAtConv_b); + plotter_track_ext.region(2).plot(hTrkXAtConv_a); + plotter_track_ext.region(3).plot(hTrkYAtConv_t); + plotter_track_ext.region(4).plot(hTrkYAtConv_b); + plotter_track_ext.region(5).plot(hTrkYAtConv_a); + + IHistogram hTrkHMXAtConv_t = aida.histogram1D("Top track Pz>1GeV X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConv_t = aida.histogram1D("Top track Pz>1GeV Y @ -67cm" , 50, -20, 20); + IHistogram hTrkHMXAtConv_b = aida.histogram1D("Bottom track Pz>1GeV X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConv_b = aida.histogram1D("Bottom track Pz>1GeV Y @ -67cm" , 50, -20, 20); + IHistogram hTrkHMXAtConv_a = aida.histogram1D("Track Pz>1GeV X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConv_a = aida.histogram1D("Track Pz>1GeV Y @ -67cm" , 50, -20, 20); + plotter_track_ext.region(0).plot(hTrkHMXAtConv_t,"mode=overlay"); + plotter_track_ext.region(1).plot(hTrkHMXAtConv_b,"mode=overlay"); + plotter_track_ext.region(2).plot(hTrkHMXAtConv_a,"mode=overlay"); + plotter_track_ext.region(3).plot(hTrkHMYAtConv_t,"mode=overlay"); + plotter_track_ext.region(4).plot(hTrkHMYAtConv_b,"mode=overlay"); + plotter_track_ext.region(5).plot(hTrkHMYAtConv_a,"mode=overlay"); + + IPlotter plotter_track_ext_entr = af.createPlotterFactory().create(); + plotter_track_ext_entr.createRegions(2,3,0); + plotter_track_ext_entr.setTitle("Track @ 0cm"); + plotterFrame.addPlotter(plotter_track_ext_entr); + IHistogram hTrkXAtEntr_t = aida.histogram1D("Top track X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntr_t = aida.histogram1D("Top track Y @ 0cm" , 50, -40, 40); + IHistogram hTrkXAtEntr_b = aida.histogram1D("Bottom track X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntr_b = aida.histogram1D("Bottom track Y @ 0cm" , 50, -40, 40); + IHistogram hTrkXAtEntr_a = aida.histogram1D("Track X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntr_a = aida.histogram1D("Track Y @ 0cm" , 50, -40, 40); + plotter_track_ext_entr.region(0).plot(hTrkXAtEntr_t); + plotter_track_ext_entr.region(1).plot(hTrkXAtEntr_b); + plotter_track_ext_entr.region(2).plot(hTrkXAtEntr_a); + plotter_track_ext_entr.region(3).plot(hTrkYAtEntr_t); + plotter_track_ext_entr.region(4).plot(hTrkYAtEntr_b); + plotter_track_ext_entr.region(5).plot(hTrkYAtEntr_a); + + + IHistogram hTrkHMXAtEntr_t = aida.histogram1D("Top track Pz>1GeV X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntr_t = aida.histogram1D("Top track Pz>1GeV Y @ 0cm" , 50, -40, 40); + IHistogram hTrkHMXAtEntr_b = aida.histogram1D("Bottom track Pz>1GeV X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntr_b = aida.histogram1D("Bottom track Pz>1GeV Y @ 0cm" , 50, -40, 40); + IHistogram hTrkHMXAtEntr_a = aida.histogram1D("Track Pz>1GeV X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntr_a = aida.histogram1D("Track Pz>1GeV Y @ 0cm" , 50, -40, 40); + plotter_track_ext_entr.region(0).plot(hTrkHMXAtEntr_t,"mode=overlay"); + plotter_track_ext_entr.region(1).plot(hTrkHMXAtEntr_b,"mode=overlay"); + plotter_track_ext_entr.region(2).plot(hTrkHMXAtEntr_a,"mode=overlay"); + plotter_track_ext_entr.region(3).plot(hTrkHMYAtEntr_t,"mode=overlay"); + plotter_track_ext_entr.region(4).plot(hTrkHMYAtEntr_b,"mode=overlay"); + plotter_track_ext_entr.region(5).plot(hTrkHMYAtEntr_a,"mode=overlay"); + + + IPlotter plotter_track_ext_coll = af.createPlotterFactory().create(); + plotter_track_ext_coll.createRegions(2,3,0); + plotter_track_ext_coll.setTitle("Track @ -150cm"); + plotterFrame.addPlotter(plotter_track_ext_coll); + IHistogram hTrkXAtColl_t = aida.histogram1D("Top track X @ -150cm" , 50, -100, 100); + IHistogram hTrkYAtColl_t = aida.histogram1D("Top track Y @ -150cm" , 50, -40, 40); + IHistogram hTrkXAtColl_b = aida.histogram1D("Bottom track X @ -150cm" , 50, -100, 100); + IHistogram hTrkYAtColl_b = aida.histogram1D("Bottom track Y @ -150cm" , 50, -40, 40); + IHistogram hTrkXAtColl_a = aida.histogram1D("Track X @ -150cm" , 50, -100, 100); + IHistogram hTrkYAtColl_a = aida.histogram1D("Track Y @ -150cm" , 50, -40, 40); + plotter_track_ext_coll.region(0).plot(hTrkXAtColl_t); + plotter_track_ext_coll.region(1).plot(hTrkXAtColl_b); + plotter_track_ext_coll.region(2).plot(hTrkXAtColl_a); + plotter_track_ext_coll.region(3).plot(hTrkYAtColl_t); + plotter_track_ext_coll.region(4).plot(hTrkYAtColl_b); + plotter_track_ext_coll.region(5).plot(hTrkYAtColl_a); + + + IHistogram hTrkHMXAtColl_t = aida.histogram1D("Top track Pz>1GeV X @ -150cm" , 50, -100, 100); + IHistogram hTrkHMYAtColl_t = aida.histogram1D("Top track Pz>1GeV Y @ -150cm" , 50, -40, 40); + IHistogram hTrkHMXAtColl_b = aida.histogram1D("Bottom track Pz>1GeV X @ -150cm" , 50, -100, 100); + IHistogram hTrkHMYAtColl_b = aida.histogram1D("Bottom track Pz>1GeV Y @ -150cm" , 50, -40, 40); + IHistogram hTrkHMXAtColl_a = aida.histogram1D("Track Pz>1GeV X @ -150cm" , 50, -100, 100); + IHistogram hTrkHMYAtColl_a = aida.histogram1D("Track Pz>1GeV Y @ -150cm" , 50, -40, 40); + plotter_track_ext_coll.region(0).plot(hTrkHMXAtColl_t,"mode=overlay"); + plotter_track_ext_coll.region(1).plot(hTrkHMXAtColl_b,"mode=overlay"); + plotter_track_ext_coll.region(2).plot(hTrkHMXAtColl_a,"mode=overlay"); + plotter_track_ext_coll.region(3).plot(hTrkHMYAtColl_t,"mode=overlay"); + plotter_track_ext_coll.region(4).plot(hTrkHMYAtColl_b,"mode=overlay"); + plotter_track_ext_coll.region(5).plot(hTrkHMYAtColl_a,"mode=overlay"); + + + + + + IPlotter plotter_track_ext2 = af.createPlotterFactory().create(); + plotter_track_ext2.createRegions(2,3,0); + plotter_track_ext2.setTitle("Matched track @ Target "); + plotterFrame.addPlotter(plotter_track_ext2); + IHistogram hmTrkXAtConv_t = aida.histogram1D("Top matched track X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConv_t = aida.histogram1D("Top matched track Y @ -67cm" , 50, -20, 20); + IHistogram hmTrkXAtConv_b = aida.histogram1D("Bottom matched track X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConv_b = aida.histogram1D("Bottom matched track Y @ -67cm" , 50, -20, 20); + IHistogram hmTrkXAtConv_a = aida.histogram1D("Matched track X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConv_a = aida.histogram1D("Matched track Y @ -67cm" , 50, -20, 20); + plotter_track_ext2.region(0).plot(hmTrkXAtConv_t); + plotter_track_ext2.region(1).plot(hmTrkXAtConv_b); + plotter_track_ext2.region(2).plot(hmTrkXAtConv_a); + plotter_track_ext2.region(3).plot(hmTrkYAtConv_t); + plotter_track_ext2.region(4).plot(hmTrkYAtConv_b); + plotter_track_ext2.region(5).plot(hmTrkYAtConv_a); + + + + + + + IPlotter plotter_track_ext3 = af.createPlotterFactory().create(); + plotter_track_ext3.createRegions(2,3,0); + plotter_track_ext3.setTitle("Matched sel track @ Target"); + plotterFrame.addPlotter(plotter_track_ext3); + IHistogram hmsTrkXAtConv_t = aida.histogram1D("Top sel matched track X @ -67cm" , 50, -100, 100); + IHistogram hmsTrkYAtConv_t = aida.histogram1D("Top sel matched track Y @ -67cm" , 50, -100, 100); + IHistogram hmsTrkXAtConv_b = aida.histogram1D("Bottom sel matched track X @ -67cm" , 50, -100, 100); + IHistogram hmsTrkYAtConv_b = aida.histogram1D("Bottom sel matched track Y @ -67cm" , 50, -100, 100); + IHistogram hmsTrkXAtConv_a = aida.histogram1D("Matched sel track X @ -67cm" , 50, -100, 100); + IHistogram hmsTrkYAtConv_a = aida.histogram1D("Matched sel track Y @ -67cm" , 50, -100, 100); + plotter_track_ext3.region(0).plot(hmsTrkXAtConv_t); + plotter_track_ext3.region(1).plot(hmsTrkXAtConv_b); + plotter_track_ext3.region(2).plot(hmsTrkXAtConv_a); + plotter_track_ext3.region(3).plot(hmsTrkYAtConv_t); + plotter_track_ext3.region(4).plot(hmsTrkYAtConv_b); + plotter_track_ext3.region(5).plot(hmsTrkYAtConv_a); + + + + + IPlotter plotter_track_ext4 = af.createPlotterFactory().create(); + plotter_track_ext4.createRegions(2,3,0); + plotter_track_ext4.setTitle("Track @ Target Pos. Charge"); + plotterFrame.addPlotter(plotter_track_ext4); + IHistogram hTrkXAtConvqp_t = aida.histogram1D("Top track q>0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConvqp_t = aida.histogram1D("Top track q>0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkXAtConvqp_b = aida.histogram1D("Bottom track q>0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConvqp_b = aida.histogram1D("Bottom track q>0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkXAtConvqp_a = aida.histogram1D("Track q>0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConvqp_a = aida.histogram1D("Track q>0 Y @ -67cm" , 50, -20, 20); + plotter_track_ext4.region(0).plot(hTrkXAtConvqp_t); + plotter_track_ext4.region(1).plot(hTrkXAtConvqp_b); + plotter_track_ext4.region(2).plot(hTrkXAtConvqp_a); + plotter_track_ext4.region(3).plot(hTrkYAtConvqp_t); + plotter_track_ext4.region(4).plot(hTrkYAtConvqp_b); + plotter_track_ext4.region(5).plot(hTrkYAtConvqp_a); + + + IHistogram hTrkHMXAtConvqp_t = aida.histogram1D("Top track Pz>1GeV q>0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConvqp_t = aida.histogram1D("Top track Pz>1GeV q>0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkHMXAtConvqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConvqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkHMXAtConvqp_a = aida.histogram1D("Track Pz>1GeV q>0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConvqp_a = aida.histogram1D("Track Pz>1GeV q>0 Y @ -67cm" , 50, -20, 20); + plotter_track_ext4.region(0).plot(hTrkHMXAtConvqp_t,"mode=overlay"); + plotter_track_ext4.region(1).plot(hTrkHMXAtConvqp_b,"mode=overlay"); + plotter_track_ext4.region(2).plot(hTrkHMXAtConvqp_a,"mode=overlay"); + plotter_track_ext4.region(3).plot(hTrkHMYAtConvqp_t,"mode=overlay"); + plotter_track_ext4.region(4).plot(hTrkHMYAtConvqp_b,"mode=overlay"); + plotter_track_ext4.region(5).plot(hTrkHMYAtConvqp_a,"mode=overlay"); + + + + IPlotter plotter_track_ext4_entr = af.createPlotterFactory().create(); + plotter_track_ext4_entr.createRegions(2,3,0); + plotter_track_ext4_entr.setTitle("Track @ 0cm Pos. Charge"); + plotterFrame.addPlotter(plotter_track_ext4_entr); + IHistogram hTrkXAtEntrqp_t = aida.histogram1D("Top track q>0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntrqp_t = aida.histogram1D("Top track q>0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkXAtEntrqp_b = aida.histogram1D("Bottom track q>0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntrqp_b = aida.histogram1D("Bottom track q>0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkXAtEntrqp_a = aida.histogram1D("Track q>0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntrqp_a = aida.histogram1D("Track q>0 Y @ 0cm" , 50, -40, 40); + plotter_track_ext4_entr.region(0).plot(hTrkXAtEntrqp_t); + plotter_track_ext4_entr.region(1).plot(hTrkXAtEntrqp_b); + plotter_track_ext4_entr.region(2).plot(hTrkXAtEntrqp_a); + plotter_track_ext4_entr.region(3).plot(hTrkYAtEntrqp_t); + plotter_track_ext4_entr.region(4).plot(hTrkYAtEntrqp_b); + plotter_track_ext4_entr.region(5).plot(hTrkYAtEntrqp_a); + + + IHistogram hTrkHMXAtEntrqp_t = aida.histogram1D("Top track Pz>1GeV q>0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntrqp_t = aida.histogram1D("Top track Pz>1GeV q>0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkHMXAtEntrqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntrqp_b = aida.histogram1D("Bottom track Pz>1GeV q>0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkHMXAtEntrqp_a = aida.histogram1D("Track Pz>1GeV q>0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntrqp_a = aida.histogram1D("Track Pz>1GeV q>0 Y @ 0cm" , 50, -40, 40); + plotter_track_ext4_entr.region(0).plot(hTrkHMXAtEntrqp_t,"mode=overlay"); + plotter_track_ext4_entr.region(1).plot(hTrkHMXAtEntrqp_b,"mode=overlay"); + plotter_track_ext4_entr.region(2).plot(hTrkHMXAtEntrqp_a,"mode=overlay"); + plotter_track_ext4_entr.region(3).plot(hTrkHMYAtEntrqp_t,"mode=overlay"); + plotter_track_ext4_entr.region(4).plot(hTrkHMYAtEntrqp_b,"mode=overlay"); + plotter_track_ext4_entr.region(5).plot(hTrkHMYAtEntrqp_a,"mode=overlay"); + + + + IPlotter plotter_track_ext5 = af.createPlotterFactory().create(); + plotter_track_ext5.createRegions(2,3,0); + plotter_track_ext5.setTitle("Track @ Target Neg. Charge"); + plotterFrame.addPlotter(plotter_track_ext5); + IHistogram hTrkXAtConvqm_t = aida.histogram1D("Top track q<0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConvqm_t = aida.histogram1D("Top track q<0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkXAtConvqm_b = aida.histogram1D("Bottom track q<0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConvqm_b = aida.histogram1D("Bottom track q<0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkXAtConvqm_a = aida.histogram1D("Track q<0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkYAtConvqm_a = aida.histogram1D("Track q<0 Y @ -67cm" , 50, -20, 20); + plotter_track_ext5.region(0).plot(hTrkXAtConvqm_t); + plotter_track_ext5.region(1).plot(hTrkXAtConvqm_b); + plotter_track_ext5.region(2).plot(hTrkXAtConvqm_a); + plotter_track_ext5.region(3).plot(hTrkYAtConvqm_t); + plotter_track_ext5.region(4).plot(hTrkYAtConvqm_b); + plotter_track_ext5.region(5).plot(hTrkYAtConvqm_a); + + IHistogram hTrkHMXAtConvqm_t = aida.histogram1D("Top track Pz>1GeV q<0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConvqm_t = aida.histogram1D("Top track Pz>1GeV q<0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkHMXAtConvqm_b = aida.histogram1D("Bottom track Pz>1GeV q<0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConvqm_b = aida.histogram1D("Bottom track Pz>1GeV q<0 Y @ -67cm" , 50, -20, 20); + IHistogram hTrkHMXAtConvqm_a = aida.histogram1D("Track Pz>1GeV q<0 X @ -67cm" , 50, -100, 100); + IHistogram hTrkHMYAtConvqm_a = aida.histogram1D("Track Pz>1GeV q<0 Y @ -67cm" , 50, -20, 20); + plotter_track_ext5.region(0).plot(hTrkHMXAtConvqm_t,"mode=overlay"); + plotter_track_ext5.region(1).plot(hTrkHMXAtConvqm_b,"mode=overlay"); + plotter_track_ext5.region(2).plot(hTrkHMXAtConvqm_a,"mode=overlay"); + plotter_track_ext5.region(3).plot(hTrkHMYAtConvqm_t,"mode=overlay"); + plotter_track_ext5.region(4).plot(hTrkHMYAtConvqm_b,"mode=overlay"); + plotter_track_ext5.region(5).plot(hTrkHMYAtConvqm_a,"mode=overlay"); + + + IPlotter plotter_track_ext5_entr = af.createPlotterFactory().create(); + plotter_track_ext5_entr.createRegions(2,3,0); + plotter_track_ext5_entr.setTitle("Track @ 0cm Neg. Charge"); + plotterFrame.addPlotter(plotter_track_ext5_entr); + IHistogram hTrkXAtEntrqm_t = aida.histogram1D("Top track q<0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntrqm_t = aida.histogram1D("Top track q<0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkXAtEntrqm_b = aida.histogram1D("Bottom track q<0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntrqm_b = aida.histogram1D("Bottom track q<0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkXAtEntrqm_a = aida.histogram1D("Track q<0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkYAtEntrqm_a = aida.histogram1D("Track q<0 Y @ 0cm" , 50, -40, 40); + plotter_track_ext5_entr.region(0).plot(hTrkXAtEntrqm_t); + plotter_track_ext5_entr.region(1).plot(hTrkXAtEntrqm_b); + plotter_track_ext5_entr.region(2).plot(hTrkXAtEntrqm_a); + plotter_track_ext5_entr.region(3).plot(hTrkYAtEntrqm_t); + plotter_track_ext5_entr.region(4).plot(hTrkYAtEntrqm_b); + plotter_track_ext5_entr.region(5).plot(hTrkYAtEntrqm_a); + + + + + IHistogram hTrkHMXAtEntrqm_t = aida.histogram1D("Top track Pz>1GeV q<0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntrqm_t = aida.histogram1D("Top track Pz>1GeV q<0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkHMXAtEntrqm_b = aida.histogram1D("Bottom track Pz>1GeV q<0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntrqm_b = aida.histogram1D("Bottom track Pz>1GeV q<0 Y @ 0cm" , 50, -40, 40); + IHistogram hTrkHMXAtEntrqm_a = aida.histogram1D("Track Pz>1GeV q<0 X @ 0cm" , 50, -100, 100); + IHistogram hTrkHMYAtEntrqm_a = aida.histogram1D("Track Pz>1GeV q<0 Y @ 0cm" , 50, -40, 40); + plotter_track_ext5_entr.region(0).plot(hTrkHMXAtEntrqm_t,"mode=overlay"); + plotter_track_ext5_entr.region(1).plot(hTrkHMXAtEntrqm_b,"mode=overlay"); + plotter_track_ext5_entr.region(2).plot(hTrkHMXAtEntrqm_a,"mode=overlay"); + plotter_track_ext5_entr.region(3).plot(hTrkHMYAtEntrqm_t,"mode=overlay"); + plotter_track_ext5_entr.region(4).plot(hTrkHMYAtEntrqm_b,"mode=overlay"); + plotter_track_ext5_entr.region(5).plot(hTrkHMYAtEntrqm_a,"mode=overlay"); + + + + + + + IPlotter plotter_track_ext6 = af.createPlotterFactory().create(); + plotter_track_ext6.createRegions(2,3,0); + plotter_track_ext6.setTitle("Matched Track @ Target Pos. Charge"); + plotterFrame.addPlotter(plotter_track_ext6); + IHistogram hmTrkXAtConvqp_t = aida.histogram1D("Top matched track q>0 X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConvqp_t = aida.histogram1D("Top matched track q>0 Y @ -67cm" , 50, -20, 20); + IHistogram hmTrkXAtConvqp_b = aida.histogram1D("Bottom matched track q>0 X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConvqp_b = aida.histogram1D("Bottom matched track q>0 Y @ -67cm" , 50, -20, 20); + IHistogram hmTrkXAtConvqp_a = aida.histogram1D("Matched track q>0 X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConvqp_a = aida.histogram1D("Matched track q>0 Y @ -67cm" , 50, -20, 20); + plotter_track_ext6.region(0).plot(hmTrkXAtConvqp_t); + plotter_track_ext6.region(1).plot(hmTrkXAtConvqp_b); + plotter_track_ext6.region(2).plot(hmTrkXAtConvqp_a); + plotter_track_ext6.region(3).plot(hmTrkYAtConvqp_t); + plotter_track_ext6.region(4).plot(hmTrkYAtConvqp_b); + plotter_track_ext6.region(5).plot(hmTrkYAtConvqp_a); + + IPlotter plotter_track_ext7 = af.createPlotterFactory().create(); + plotter_track_ext7.createRegions(2,3,0); + plotter_track_ext7.setTitle("Matched Track @ Target Neg. Charge"); + plotterFrame.addPlotter(plotter_track_ext7); + IHistogram hmTrkXAtConvqm_t = aida.histogram1D("Top matched track q<0 X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConvqm_t = aida.histogram1D("Top matched track q<0 Y @ -67cm" , 50, -20, 20); + IHistogram hmTrkXAtConvqm_b = aida.histogram1D("Bottom matched track q<0 X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConvqm_b = aida.histogram1D("Bottom matched track q<0 Y @ -67cm" , 50, -20, 20); + IHistogram hmTrkXAtConvqm_a = aida.histogram1D("Matched track q<0 X @ -67cm" , 50, -100, 100); + IHistogram hmTrkYAtConvqm_a = aida.histogram1D("Matched track q<0 Y @ -67cm" , 50, -20, 20); + plotter_track_ext7.region(0).plot(hmTrkXAtConvqm_t); + plotter_track_ext7.region(1).plot(hmTrkXAtConvqm_b); + plotter_track_ext7.region(2).plot(hmTrkXAtConvqm_a); + plotter_track_ext7.region(3).plot(hmTrkYAtConvqm_t); + plotter_track_ext7.region(4).plot(hmTrkYAtConvqm_b); + plotter_track_ext7.region(5).plot(hmTrkYAtConvqm_a); + + + + + + IPlotter plotter_track_chi2 = af.createPlotterFactory().create(); + plotter_track_chi2.createRegions(3,3,0); + plotter_track_chi2.setTitle("Track Chi2"); + plotterFrame.addPlotter(plotter_track_chi2); + IHistogram hChi2_t = aida.histogram1D("Top track Chi2" , 50, 0, 15); + IHistogram hChi2_b = aida.histogram1D("Bottom track Chi2" , 50, 0, 15); + IHistogram hChi2_a = aida.histogram1D("Track Chi2" , 50, 0, 15); + IHistogram hmChi2_t = aida.histogram1D("Matched top track Chi2" , 50, 0, 15); + IHistogram hmChi2_b = aida.histogram1D("Matched bottom track Chi2" , 50, 0, 15); + IHistogram hmChi2_a = aida.histogram1D("Matched track Chi2" , 50, 0, 15); + IHistogram hmsChi2_t = aida.histogram1D("Matched sel top track Chi2" , 50, 0, 15); + IHistogram hmsChi2_b = aida.histogram1D("Matched sel bottom track Chi2" , 50, 0, 15); + IHistogram hmsChi2_a = aida.histogram1D("Matched sel track Chi2" , 50, 0, 15); + + plotter_track_chi2.region(0).plot(hChi2_t); + plotter_track_chi2.region(1).plot(hChi2_b); + plotter_track_chi2.region(2).plot(hChi2_a); + plotter_track_chi2.region(3).plot(hmChi2_t); + plotter_track_chi2.region(4).plot(hmChi2_b); + plotter_track_chi2.region(5).plot(hmChi2_a); + plotter_track_chi2.region(6).plot(hmsChi2_t); + plotter_track_chi2.region(7).plot(hmsChi2_b); + plotter_track_chi2.region(8).plot(hmsChi2_a); + + + + IPlotter plotter_track_chi2_2 = af.createPlotterFactory().create(); + plotter_track_chi2_2.createRegions(3,3,0); + plotter_track_chi2_2.setTitle("Track Chi2 Target Outliers X"); + plotterFrame.addPlotter(plotter_track_chi2_2); + IHistogram hChi2TX_t = aida.histogram1D("Top track !target in X Chi2" , 50, 0, 15); + IHistogram hChi2TX_b = aida.histogram1D("Bottom track !target in X Chi2" , 50, 0, 15); + IHistogram hChi2TX_a = aida.histogram1D("Track !target in X Chi2" , 50, 0, 15); + IHistogram hmChi2TX_t = aida.histogram1D("Matched top track !target in X Chi2" , 50, 0, 15); + IHistogram hmChi2TX_b = aida.histogram1D("Matched bottom track !target in X Chi2" , 50, 0, 15); + IHistogram hmChi2TX_a = aida.histogram1D("Matched track !target in X Chi2" , 50, 0, 15); + IHistogram hmsChi2TX_t = aida.histogram1D("Matched sel top track !target in X Chi2" , 50, 0, 15); + IHistogram hmsChi2TX_b = aida.histogram1D("Matched sel bottom track !target in X Chi2" , 50, 0, 15); + IHistogram hmsChi2TX_a = aida.histogram1D("Matched sel track !target in X Chi2" , 50, 0, 15); + + plotter_track_chi2_2.region(0).plot(hChi2TX_t); + plotter_track_chi2_2.region(1).plot(hChi2TX_b); + plotter_track_chi2_2.region(2).plot(hChi2TX_a); + plotter_track_chi2_2.region(3).plot(hmChi2TX_t); + plotter_track_chi2_2.region(4).plot(hmChi2TX_b); + plotter_track_chi2_2.region(5).plot(hmChi2TX_a); + plotter_track_chi2_2.region(6).plot(hmsChi2TX_t); + plotter_track_chi2_2.region(7).plot(hmsChi2TX_b); + plotter_track_chi2_2.region(8).plot(hmsChi2TX_a); + + + + + + IHistogram hmcls_t = aida.histogram1D("Cluster size matched top", 6, 0.0, 6.0); + IHistogram hmcls_b = aida.histogram1D("Cluster size matched bottom", 6, 0.0, 6.0); + IHistogram hmcls_a = aida.histogram1D("Cluster size matched", 6, 0.0, 6.0); + IHistogram hmscls_t = aida.histogram1D("Cluster size matched sel top", 6, 0.0, 6.0); + IHistogram hmscls_b = aida.histogram1D("Cluster size matched sel bottom", 6, 0.0, 6.0); + IHistogram hmscls_a = aida.histogram1D("Cluster size matched sel", 6, 0.0, 6.0); + + + + + + + IPlotter plotter_ecal_clsxm = af.createPlotterFactory().create(); + plotter_ecal_clsxm.createRegions(2,3,0); + plotter_ecal_clsxm.setTitle("Ecal Cluster size matched"); + plotterFrame.addPlotter(plotter_ecal_clsxm); + + plotter_ecal_clsxm.region(0).plot(hmcls_t); + plotter_ecal_clsxm.region(1).plot(hmcls_b); + plotter_ecal_clsxm.region(2).plot(hmcls_a); + plotter_ecal_clsxm.region(3).plot(hmscls_t); + plotter_ecal_clsxm.region(4).plot(hmscls_b); + plotter_ecal_clsxm.region(5).plot(hmscls_a); + + + + + + IHistogram hmHMcls_t = aida.histogram1D("Cluster size matched trk Pz>600MeV top", 6, 0.0, 6.0); + IHistogram hmHMcls_b = aida.histogram1D("Cluster size matched trk Pz>600MeV bottom", 6, 0.0, 6.0); + IHistogram hmHMcls_a = aida.histogram1D("Cluster size matched trk Pz>600MeV", 6, 0.0, 6.0); + IHistogram hmHM2cls_t = aida.histogram1D("Cluster size matched trk Pz>1000MeV top", 6, 0.0, 6.0); + IHistogram hmHM2cls_b = aida.histogram1D("Cluster size matched trk Pz>1000MeV bottom", 6, 0.0, 6.0); + IHistogram hmHM2cls_a = aida.histogram1D("Cluster size matched trk Pz>1000MeV", 6, 0.0, 6.0); + + + + + IPlotter plotter_ecal_clsxm2 = af.createPlotterFactory().create(); + plotter_ecal_clsxm2.createRegions(2,3,0); + plotter_ecal_clsxm2.setTitle("Ecal Cluster size matched Trk Pz");[truncated at 1000 lines; 1686 more skipped]
Use REPLY-ALL to reply to list
To unsubscribe from the LCD-CVS list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=LCD-CVS&A=1