hps-java/src/main/java/org/lcsim/hps/users/phansson
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]