APEEstimation/macros/commandsDrawPlot.C

0001 {
0002
0003
0004
0005
0006 gROOT->SetStyle("Plain");
0007 gROOT->ForceStyle();
0008 gStyle->SetOptStat(0);
0009
0010 gStyle->SetPalette(1);
0011 //gStyle->SetNumberContours(20);  // Default: 20
0012
0013 double width = 600.;
0014
0015 gStyle->SetCanvasDefW(width);
0016 gStyle->SetCanvasDefH(width);
0017
0018 gStyle->SetPadLeftMargin(0.13);
0019 gStyle->SetPadRightMargin(0.11);
0020
0021 gStyle->SetPadTopMargin(0.10);
0022 gStyle->SetPadBottomMargin(0.13);
0023
0024 gStyle->SetTitleOffset(1.25,"Y");
0025 gStyle->SetTitleOffset(1.13,"X");
0026
0027
0028 TGaxis::SetMaxDigits(3);
0029
0030 gStyle->SetTitleX(0.25);        // move upper left corner of title box to specified value
0031 //gStyle->SetTitleY(0.99);        // move upper left corner of title box to specified value
0032
0033
0034 gStyle->SetTitleXSize(0.05);
0035 gStyle->SetTitleYSize(0.05);
0036 gStyle->SetTitleSize(0.05,"XY");
0037 gStyle->SetLabelSize(0.05,"XY");
0038
0039
0040
0041
0042 //++++++++++++++++++++++++++++++++++=====================================+++++++++++++++++++++++++++++++
0043
0044
0045 /*
0046
0047 gROOT->ProcessLine(".L tdrstyle.C");
0048 setTDRStyle();
0049
0050 */
0051
0052
0053 //++++++++++++++++++++++++++++++++++=====================================+++++++++++++++++++++++++++++++
0054
0055
0056
0057
0058
0059 gROOT->ProcessLine(".L DrawPlot.C");
0060
0061
0062
0063
0064 gROOT->ProcessLine("DrawPlot drawPlot1(0, true)");
0065
0066 drawPlot1.setLegendEntry("data (final #sigma_{align})","data (#sigma_{align}=0)","design");
0067 gStyle->SetOptStat(111110);
0068 //drawPlot1.setLegendCoordinate(0.69, 0.70, 0.99, 0.85);
0069 drawPlot1.setLegendCoordinate(0.69, 0.20, 0.99, 0.35);
0070 drawPlot1.drawPlot("ApeEstimatorSummary","h_weightX", false);
0071 drawPlot1.drawPlot("ApeEstimatorSummary","h_weightY", false);
0072 drawPlot1.setLegendCoordinate(0.69, 0.20, 0.99, 0.35);
0073 drawPlot1.drawPlot("ApeEstimatorSummary","h_residualWidthX1", false);
0074 drawPlot1.drawPlot("ApeEstimatorSummary","h_residualWidthY1", false);
0075 drawPlot1.drawPlot("ApeEstimatorSummary","h_rmsX", false);
0076 drawPlot1.drawPlot("ApeEstimatorSummary","h_rmsY", false);
0077
0078
0079 gROOT->ProcessLine("DrawPlot drawPlot2(0, false)");
0080
0081 drawPlot2.setLegendEntry("data (final #sigma_{align})","data (#sigma_{align}=0)","design");
0082 //drawPlot2.setLegendCoordinate(0.69, 0.70, 0.99, 0.85);
0083 drawPlot2.setLegendCoordinate(0.69, 0.20, 0.99, 0.35);
0084 gStyle->SetOptStat(111100);
0085 drawPlot2.drawPlot("ApeEstimator1","h_entriesX", false);
0086 drawPlot2.drawPlot("ApeEstimator1","h_entriesY", false);
0087 //drawPlot2.setLegendCoordinate(0.395, 0.845, 0.995, 0.995);
0088 gStyle->SetOptStat(111110);
0089 drawPlot2.drawPlot("ApeEstimator1","h_NorResX", false);
0090 drawPlot2.drawPlot("ApeEstimator1","h_NorResY", false);
0091 drawPlot2.drawPlot("ApeEstimator1","h_ResX", false);
0092 drawPlot2.drawPlot("ApeEstimator1","h_ResY", false);
0093
0094 /*
0095 // Additional plots for strip (mainly; some are common with pixel)
0096 drawPlot2.drawPlot("ApeEstimator2","h_ChargeStrip", false);
0097 drawPlot2.drawPlot("ApeEstimator2","h_SOverN", false);
0098 drawPlot2.drawPlot("ApeEstimator2","h_ChargeOnEdges", false);
0099 drawPlot2.drawPlot("ApeEstimator2","h_BaryStripX", false);
0100 drawPlot2.drawPlot("ApeEstimator2","h_WidthX", false);
0101 drawPlot2.drawPlot("ApeEstimator2","h_PhiSensX", false);
0102 drawPlot2.drawPlot("ApeEstimator2","h_sigmaXHit_1", false);
0103 drawPlot2.drawPlot("ApeEstimator2","h_sigmaXTrk_1", false);
0104
0105 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsChargeStrip", false);
0106 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsSOverN", false);
0107 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsChargeOnEdges", false);
0108 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsBaryStripX", false);
0109 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsPhiSensX", false);
0110
0111 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsD0Beamspot", false);
0112 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsDz", false);
0113 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsHitsPixel", false);
0114 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsHitsStrip", false);
0115 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsP", false);
0116 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsPhiSensX", false);
0117 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsPt", false);
0118 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXTrkVsTheta", false);
0119
0120 drawPlot2.drawPlot("ApeEstimator2","p_widthVsPhiSensX", false);
0121
0122 // Additional plots for pixel
0123 drawPlot2.drawPlot("ApeEstimator2","h_ChargePixel", false);
0124 drawPlot2.drawPlot("ApeEstimator2","h_ChargePixel_y", false);
0125 drawPlot2.drawPlot("ApeEstimator2","h_BaryStripY_y", false);
0126 drawPlot2.drawPlot("ApeEstimator2","h_ClusterProbXY", false);
0127 drawPlot2.drawPlot("ApeEstimator2","h_ClusterProbXY_y", false);
0128 drawPlot2.drawPlot("ApeEstimator2","h_LogClusterProb", false);
0129 drawPlot2.drawPlot("ApeEstimator2","h_LogClusterProb_y", false);
0130 drawPlot2.drawPlot("ApeEstimator2","h_PhiSensY_y", false);
0131 drawPlot2.drawPlot("ApeEstimator2","h_WidthY_y", false);
0132 drawPlot2.drawPlot("ApeEstimator2","h_sigmaYHit_1", false);
0133 drawPlot2.drawPlot("ApeEstimator2","h_sigmaYTrk_1", false);
0134
0135 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsChargePixel", false);
0136 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsClusterProbXY", false);
0137 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsLogClusterProb", false);
0138 drawPlot2.drawPlot("ApeEstimator2","p_sigmaXHitVsWidthX", false);
0139
0140 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYHitVsChargePixel", false);
0141 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYHitVsBaryStripY", false);
0142 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYHitVsClusterProbXY", false);
0143 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYHitVsLogClusterProb", false);
0144 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYHitVsWidthY", false);
0145 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYHitVsPhiSensY", false);
0146
0147 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsD0Beamspot", false);
0148 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsDz", false);
0149 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsHitsPixel", false);
0150 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsHitsStrip", false);
0151 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsP", false);
0152 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsPhiSensX", false);
0153 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsPt", false);
0154 drawPlot2.drawPlot("ApeEstimator2","p_sigmaYTrkVsTheta", false);
0155
0156 drawPlot2.drawPlot("ApeEstimator2","p_phiSensXVsBarycentreX", false);
0157 drawPlot2.drawPlot("ApeEstimator2","p_phiSensYVsBarycentreY", false);
0158 */
0159
0160
0161 //drawPlot2.setLegendCoordinate(0.395, 0.845, 0.995, 0.995);
0162 gStyle->SetOptStat(111110);
0163 drawPlot2.drawTrackPlot("ApeEstimator2","h_eta", false);
0164 drawPlot2.drawTrackPlot("ApeEstimator2","h_etaErr", false);
0165 drawPlot2.drawTrackPlot("ApeEstimator2","h_etaSig", false);
0166 drawPlot2.drawTrackPlot("ApeEstimator2","h_theta", false);
0167 drawPlot2.drawTrackPlot("ApeEstimator2","h_phi", false);
0168 drawPlot2.drawTrackPlot("ApeEstimator2","h_phiErr", false);
0169 drawPlot2.drawTrackPlot("ApeEstimator2","h_phiSig", false);
0170 drawPlot2.drawTrackPlot("ApeEstimator2","h_norChi2", false);
0171 drawPlot2.drawTrackPlot("ApeEstimator2","h_p", false);
0172 drawPlot2.drawTrackPlot("ApeEstimator2","h_pt", false);
0173 drawPlot2.drawTrackPlot("ApeEstimator2","h_ptErr", false);
0174 drawPlot2.drawTrackPlot("ApeEstimator2","h_ptSig", false);
0175 drawPlot2.drawTrackPlot("ApeEstimator2","h_charge", false);
0176 drawPlot2.drawTrackPlot("ApeEstimator2","h_d0Beamspot", false);
0177 drawPlot2.drawTrackPlot("ApeEstimator2","h_d0BeamspotErr", false);
0178 drawPlot2.drawTrackPlot("ApeEstimator2","h_d0BeamspotSig", false);
0179 drawPlot2.drawTrackPlot("ApeEstimator2","h_dz", false);
0180 drawPlot2.drawTrackPlot("ApeEstimator2","h_dzErr", false);
0181 drawPlot2.drawTrackPlot("ApeEstimator2","h_dzSig", false);
0182 drawPlot2.drawTrackPlot("ApeEstimator2","h_prob", false);
0183 drawPlot2.drawTrackPlot("ApeEstimator2","h_hitsPixel", false);
0184 drawPlot2.drawTrackPlot("ApeEstimator2","h_hitsStrip", false);
0185
0186 drawPlot2.drawEventPlot("ApeEstimator2","h_trackSizeGood", false);
0187
0188
0189
0190
0191
0192 gROOT->ProcessLine("DrawPlot drawPlot3(15, true)");
0193 //drawPlot3.thesisMode();
0194
0195 drawPlot3.setLegendEntry("data (final #sigma_{align})","data (#sigma_{align}=0)","design");
0196 drawPlot3.setLegendCoordinate(0.41, 0.27, 0.71, 0.42);
0197 gStyle->SetOptStat(111110);
0198 // drawPlot3.drawPlot("ApeEstimatorSummary","h_correctionX1", false);
0199 // drawPlot3.drawPlot("ApeEstimatorSummary","h_correctionX2", false);
0200 // drawPlot3.drawPlot("ApeEstimatorSummary","h_correctionY1", false);
0201 // drawPlot3.drawPlot("ApeEstimatorSummary","h_correctionY2", false);
0202 // drawPlot3.drawPlot("ApeEstimatorSummary","h_fitMeanX1", false);
0203 // drawPlot3.drawPlot("ApeEstimatorSummary","h_fitMeanX2", false);
0204 // drawPlot3.drawPlot("ApeEstimatorSummary","h_fitMeanY1", false);
0205 // drawPlot3.drawPlot("ApeEstimatorSummary","h_fitMeanY2", false);
0206 // drawPlot3.drawPlot("ApeEstimatorSummary","meanX", false);
0207 // drawPlot3.drawPlot("ApeEstimatorSummary","meanY", false);
0208 drawPlot3.setLegendCoordinate(0.69, 0.20, 0.99, 0.35);
0209 //drawPlot3.setLegendCoordinate(0.69, 0.70, 0.99, 0.85);
0210 drawPlot3.drawPlot("ApeEstimatorSummary","h_weightX", false);
0211 drawPlot3.drawPlot("ApeEstimatorSummary","h_weightY", false);
0212 drawPlot3.setLegendCoordinate(0.69, 0.20, 0.99, 0.35);
0213 drawPlot3.drawPlot("ApeEstimatorSummary","h_residualWidthX1", false);
0214 drawPlot3.drawPlot("ApeEstimatorSummary","h_residualWidthY1", false);
0215 //drawPlot3.drawPlot("ApeEstimatorSummary","h_residualWidthX2", false);
0216 //drawPlot3.drawPlot("ApeEstimatorSummary","h_residualWidthY2", false);
0217 drawPlot3.drawPlot("ApeEstimatorSummary","h_rmsX", false);
0218 drawPlot3.drawPlot("ApeEstimatorSummary","h_rmsY", false);
0219
0220
0221 gROOT->ProcessLine("DrawPlot drawPlot4(15, false)");
0222 //drawPlot4.thesisMode();
0223
0224 drawPlot4.setLegendEntry("data (final #sigma_{align})","data (#sigma_{align}=0)","design");
0225 drawPlot4.setLegendCoordinate(0.69, 0.20, 0.99, 0.35);
0226 //drawPlot4.setLegendCoordinate(0.69, 0.70, 0.99, 0.85);
0227 gStyle->SetOptStat(111100);
0228 drawPlot4.drawPlot("ApeEstimator1","h_entriesX", false);
0229 drawPlot4.drawPlot("ApeEstimator1","h_entriesY", false);
0230 //drawPlot4.setLegendCoordinate(0.395, 0.845, 0.995, 0.995);
0231 gStyle->SetOptStat(111110);
0232 drawPlot4.drawPlot("ApeEstimator1","h_NorResX");
0233 drawPlot4.drawPlot("ApeEstimator1","h_NorResY");
0234 drawPlot4.drawPlot("ApeEstimator1","h_ResX");
0235 drawPlot4.drawPlot("ApeEstimator1","h_ResY");
0236
0237 //drawPlot4.setLegendCoordinate(0.395, 0.845, 0.995, 0.995);
0238 gStyle->SetOptStat(111110);
0239 drawPlot4.drawTrackPlot("ApeEstimator2","h_eta");
0240 drawPlot4.drawTrackPlot("ApeEstimator2","h_etaErr");
0241 drawPlot4.drawTrackPlot("ApeEstimator2","h_etaSig");
0242 drawPlot4.drawTrackPlot("ApeEstimator2","h_theta");
0243 drawPlot4.drawTrackPlot("ApeEstimator2","h_phi");
0244 drawPlot4.drawTrackPlot("ApeEstimator2","h_phiErr");
0245 drawPlot4.drawTrackPlot("ApeEstimator2","h_phiSig");
0246 drawPlot4.drawTrackPlot("ApeEstimator2","h_norChi2");
0247 drawPlot4.drawTrackPlot("ApeEstimator2","h_p");
0248 drawPlot4.drawTrackPlot("ApeEstimator2","h_pt");
0249 drawPlot4.drawTrackPlot("ApeEstimator2","h_ptErr");
0250 drawPlot4.drawTrackPlot("ApeEstimator2","h_ptSig");
0251 drawPlot4.drawTrackPlot("ApeEstimator2","h_charge");
0252 drawPlot4.drawTrackPlot("ApeEstimator2","h_d0Beamspot");
0253 drawPlot4.drawTrackPlot("ApeEstimator2","h_d0BeamspotErr");
0254 drawPlot4.drawTrackPlot("ApeEstimator2","h_d0BeamspotSig");
0255 drawPlot4.drawTrackPlot("ApeEstimator2","h_dz");
0256 drawPlot4.drawTrackPlot("ApeEstimator2","h_dzErr");
0257 drawPlot4.drawTrackPlot("ApeEstimator2","h_dzSig");
0258 drawPlot4.drawTrackPlot("ApeEstimator2","h_prob");
0259 drawPlot4.drawTrackPlot("ApeEstimator2","h_hitsPixel");
0260 drawPlot4.drawTrackPlot("ApeEstimator2","h_hitsStrip");
0261
0262 drawPlot4.drawEventPlot("ApeEstimator2","h_trackSizeGood");
0263
0264
0265
0266 gROOT->ProcessLine(".q");
0267
0268
0269
0270 }
0271
0272
0273
0274
0275
0276