APEEstimation/macros/commandsDrawPlot_thesis.C

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