File indexing completed on 2024-04-06 12:05:40
0001 #include <iostream>
0002 #include <fstream>
0003 #include "TH1.h"
0004 #include "TH2.h"
0005 #include "TH1F.h"
0006 #include "TH2F.h"
0007 #include "TCanvas.h"
0008 #include "TROOT.h"
0009 #include <TMath.h>
0010 #include "TStyle.h"
0011 #include "TSystem.h"
0012 #include "TLegend.h"
0013 #include "TText.h"
0014 #include "TAxis.h"
0015 #include "TFile.h"
0016 #include "TLine.h"
0017 #include "TGraph.h"
0018 #include <TPaveText.h>
0019 #include <TPostScript.h>
0020 #include <TClass.h>
0021 #include <TLatex.h>
0022 #include <TFrame.h>
0023 #include <TPad.h>
0024
0025 #include "TGraphErrors.h"
0026 #include "TF1.h"
0027 #include "TArrow.h"
0028
0029
0030 using namespace std;
0031
0032
0033 int main(int argc, char *argv[]) {
0034
0035
0036 TGraph *gr1;
0037 TGraph *gr2;
0038 TGraph *gr3;
0039 TGraph *gr4;
0040 TGraph *gr5;
0041 TGraph *gr6;
0042 TGraph *gr7;
0043 TGraph *gr8;
0044 TGraph *gr9;
0045 TGraphErrors *ge1;
0046 TGraphErrors *ge2;
0047 TGraphErrors *ge3;
0048 TGraphErrors *ge4;
0049 TGraphErrors *ge5;
0050 TGraphErrors *ge6;
0051 TGraphErrors *ge7;
0052 TGraphErrors *ge8;
0053 TGraphErrors *ge9;
0054
0055
0056 const Int_t NMethods = 38;
0057 const Int_t NMethodsStart = 0;
0058 const Int_t NMethodsToUse = 37;
0059
0060 TString MethodName[NMethods] = {"PLOT0", "PLOT1", "PLOT2", "PLOT3", "PLOT4", "PLOT5", "PLOT6", "PLOT7",
0061 "PLOT8", "PLOT9", "PLOT10", "PLOT11", "PLOT12", "PLOT13", "PLOT14", "PLOT15",
0062 "PLOT16", "PLOT17", "PLOT18", "PLOT19", "PLOT20", "PLOT21", "PLOT22", "PLOT23",
0063 "PLOT24", "PLOT25", "PLOT26", "PLOT27", "PLOT28", "PLOT29", "PLOT30", "PLOT31",
0064 "PLOT32", "PLOT33", "PLOT34", "PLOT35", "PLOT36", "PLOT37"};
0065
0066
0067 const Int_t NMAX0 = 220;
0068
0069
0070 Float_t LU[NMAX0] = {
0071 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.2,
0072 1.3, 1.6, 1.9, 2.1, 2.4, 2.7, 3.0, 3.0, 3.0, 3.0, 3.0, 3.0, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0,
0073 6.5, 6.9, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.4, 7.9, 8.3, 8.7, 9.1, 9.5, 10.0, 10.5, 10.9, 11.3, 11.7,
0074 12.1, 12.5, 13.0, 13.4, 13.8, 14.2, 14.6, 15.0, 15.4, 15.8, 16.2, 16.6, 17.0, 17.4, 17.8, 18.2, 18.7, 19.1, 19.5,
0075 19.5, 19.5, 19.5, 19.5, 19.5, 19.5, 19.5, 19.8, 20.1, 20.4, 20.8, 21.1, 21.5, 21.9, 22.2, 22.5, 22.5, 22.5, 22.5,
0076 23.0, 23.5, 24.0, 24.5, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.2, 25.4, 26.0, 26.5, 27.0, 27.5,
0077 28.0, 28.5, 29.0, 29.5, 30.0, 30.5, 31.0, 31.5, 31.5, 31.5, 31.5, 31.5, 31.5, 31.5, 31.5, 31.5, 31.5, 31.5, 31.5,
0078 31.5, 31.5, 31.5, 31.5, 31.5, 32.6, 32.8, 33.3, 33.7, 34.0, 34.2, 34.2, 34.2, 34.2, 34.2, 34.2, 34.2, 34.2, 34.2,
0079 34.2, 34.9, 35.4, 35.9, 36.4, 36.9, 37.4, 37.8, 38.3, 38.6, 39.0, 39.4, 40.0, 40.5, 41.0, 41.3, 41.3, 41.3, 41.3,
0080 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3,
0081 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 41.3, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0082 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.};
0083
0084
0085 Float_t LU0[NMAX0] = {
0086 24.6, 50.4, 47.7, 28.0, 122.5, 0., 98.5, 0.02, 37., 256.6, 0., 0., 0., 0., 0., 0.,
0087 9.4, 5.1, 101.8, 130.6, 190., 198.8, 228.5, 320.3, 301., 338.6, 219.4, 0., 0., 0., 0., 0.,
0088 443., 464.6, 0., 365., 191.4, 158.8, 387.8, 182., 740.1, 0., 0., 0., 73.9, 472., 550., 0.,
0089 237.4, 345.1, 104.9, 0., 608.4, 340.8, 617.7, 0., 500., 532.6, 485., 360.6, 491.6, 289.7, 525., 0.,
0090 541.5, 481.4, 529.1, 172., 470., 611., 0., 673.9, 518.7, 756.5, 250.4, 61.3, 0., 0., 0., 0.,
0091 0., 0., 385.4, 0., 324.7, 367.4, 39.2, 412.7, 468., 298., 274., 0., 0., 0., 580.7, 482.4,
0092 200.2, 249.9, 407.3, 126.6, 40.3, 180., 0., 0., 0., 0., 238.7, 391.6, 284.9, 805.2, 513.6, 224.7,
0093 464.5, 265.2, 471.3, 857.4, 0., 286.5, 428.7, 0., 505.5, 374.3, 0., 0., 0., 0., 0., 0.,
0094 0., 0., 0., 0., 0., 0., 0., 0., 0., 17.4, 519.7, 469.4, 397.9, 0., 752., 0.,
0095 392.6, 0., 0., 24.9, 0., 10.4, 494.4, 441.7, 211.8, 152.6, 217.4, 600.7, 350.8, 434.6, 500., 127.6,
0096 400., 126.8, 399., 271.3, 819., 322.8, 364.6, 106.4, 350., 111.1, 385.3, 0.0, 0.9, 102.5, 130.3, 0.0,
0097 0.112, 0.071, 0.036, 1.771, 1.004, 12.40, 10.81, 3.80, 5.396, 8.971, 0.006, 11.17, 4.364, 0.671, 1.310, 7.819,
0098 8.927, 7.179, 8.289, 8.268, 4.788, 5.304, 5.443, 0.668, 1.128, 2.569, 7.931, 10.70, 7.068, 10.96, 10.29, 5.189,
0099 0.112, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.};
0100
0101 Float_t LU1[NMAX0];
0102 Float_t klumi = 1.076;
0103
0104 int ncounts = 1;
0105 LU1[ncounts - 1] = klumi * LU0[ncounts - 1] / 1000.;
0106
0107 cout << "Terminated: ind = " << ncounts - 1 << " LU1 = " << LU1[ncounts - 1] << " LU = " << LU[ncounts - 1]
0108 << endl;
0109
0110 for (int jj = 1; jj < NMAX0; jj++) {
0111 LU1[ncounts] = LU1[ncounts - 1] + klumi * LU0[jj] / 1000.;
0112
0113 ++ncounts;
0114 }
0115 cout << "LU1[NMAX0-1] = " << LU1[NMAX0 - 1] << endl;
0116
0117
0118
0119
0120
0121
0122
0123
0124
0125
0126
0127
0128
0129
0130 const Int_t ALF = 670;
0131
0132 Int_t DeltaDaysLumiTaken[ALF] = {
0133 1, 2, 3, 5, 5, 5, 7, 7, 8, 9, 9, 11, 18, 18, 18, 19, 20, 20, 21, 22, 23, 23, 24,
0134 25, 25, 26, 27, 33, 33, 35, 36, 37, 39, 39, 40, 41, 41, 42, 46, 47, 48, 50, 51, 52, 52, 53,
0135 54, 54, 55, 56, 56, 58, 59, 60, 61, 62, 64, 64, 66, 67, 68, 69, 70, 71, 73, 74, 74, 74, 75,
0136 75, 76, 77, 77, 84, 84, 84, 86, 87, 87, 88, 88, 89, 90, 91, 91, 92, 95, 96, 96, 96, 97, 98,
0137 99, 99, 100, 100, 101, 101, 102, 102, 103, 108, 108, 108, 109, 110, 110, 111, 111, 112, 113, 114, 114, 115, 115,
0138 116, 117, 117, 118, 119, 120, 122, 123, 139, 139, 139, 140, 140, 141, 142, 144, 144, 146, 150, 152, 152, 153, 153,
0139 154, 155, 156, 156, 157, 158, 158, 159, 160, 161, 161, 162, 163, 163, 163, 165, 165, 166, 167, 167, 168, 169, 170,
0140 170, 171, 185, 185, 187, 189, 190, 191, 191, 193, 193, 193, 194, 195, 195, 197, 197, 198, 198, 199, 200, 201, 201,
0141 201, 202, 202, 203, 203, 204, 205, 205, 205, 206, 206, 206, 207, 207, 208, 208, 209, 209, 339, 371, 379, 380, 381,
0142 382, 383, 383, 384, 384, 386, 389, 390, 391, 392, 398, 398, 398, 400, 400, 401, 401, 402, 402, 403, 403, 404, 405,
0143 406, 408, 409, 409, 410, 410, 411, 411, 412, 412, 413, 414, 414, 415, 415, 415, 425, 425, 427, 428, 430, 430, 431,
0144 431, 432, 432, 434, 434, 435, 435, 435, 435, 436, 436, 437, 438, 440, 443, 444, 444, 444, 445, 445, 445, 445, 446,
0145 447, 447, 449, 450, 451, 451, 452, 452, 453, 454, 454, 455, 455, 456, 456, 459, 461, 461, 461, 462, 462, 463, 463,
0146 464, 464, 464, 464, 465, 466, 466, 467, 468, 469, 469, 470, 471, 473, 474, 474, 475, 475, 476, 477, 477, 478, 479,
0147 479, 479, 480, 480, 481, 481, 482, 483, 483, 483, 484, 484, 485, 486, 486, 486, 487, 487, 488, 488, 489, 490, 490,
0148 490, 491, 500, 501, 501, 502, 502, 503, 503, 504, 505, 505, 506, 507, 508, 508, 509, 509, 509, 510, 511, 512, 512,
0149 513, 513, 514, 514, 515, 516, 516, 517, 517, 518, 519, 520, 520, 521, 522, 523, 523, 524, 525, 525, 526, 526, 527,
0150 527, 528, 529, 529, 530, 530, 531, 531, 532, 533, 534, 535, 535, 536, 536, 537, 538, 538, 539, 539, 540, 540, 541,
0151 542, 543, 545, 548, 548, 550, 551, 551, 551, 552, 553, 554, 554, 554, 707, 707, 708, 709, 710, 711, 711, 712, 712,
0152 714, 715, 715, 716, 716, 717, 717, 718, 718, 719, 719, 720, 721, 722, 722, 722, 723, 724, 724, 725, 725, 725, 726,
0153 726, 727, 727, 728, 729, 730, 730, 731, 732, 732, 733, 734, 735, 735, 736, 736, 736, 737, 738, 738, 739, 739, 740,
0154 740, 741, 742, 742, 743, 744, 745, 746, 747, 748, 748, 748, 750, 750, 751, 751, 752, 752, 753, 754, 755, 755, 756,
0155 757, 757, 759, 759, 760, 761, 761, 762, 762, 776, 777, 777, 778, 779, 779, 780, 781, 782, 783, 783, 784, 784, 785,
0156 785, 786, 787, 788, 789, 789, 790, 792, 792, 793, 794, 795, 795, 796, 796, 796, 797, 797, 798, 799, 800, 800, 801,
0157 801, 802, 802, 803, 803, 804, 811, 812, 813, 813, 814, 814, 815, 815, 816, 817, 817, 818, 819, 819, 819, 820, 821,
0158 821, 822, 823, 823, 824, 825, 825, 826, 826, 826, 827, 829, 829, 830, 831, 832, 835, 835, 836, 837, 838, 838, 839,
0159 839, 840, 840, 841, 841, 841, 842, 843, 843, 843, 844, 844, 845, 845, 845, 846, 846, 847, 847, 848, 849, 849, 850,
0160 850, 850, 851, 851, 852, 852, 853, 859, 866, 866, 867, 867, 868, 868, 868, 872, 872, 873, 873, 874, 874, 876, 876,
0161 877, 877, 879, 880, 881, 881, 882, 883, 887, 887, 888, 889, 889, 890, 890, 891, 892, 892, 893, 894, 894, 895, 895,
0162 896, 899, 909
0163
0164 };
0165
0166 const Int_t NMAXDAYS = DeltaDaysLumiTaken[ALF - 1];
0167 cout << "*********** NMAXDAYS = ************* " << NMAXDAYS << endl;
0168
0169 Float_t DeliveredLumi[ALF] = {
0170 28.0, 59.1, 54.9, 33.5, 31.0, 103.2, 31.4, 72.3, 0.0, 39.3, 0.1, 269.2, 9.7,
0171 0.0, 7.1, 103.1, 6.7, 137.5, 204.2, 158.6, 57.9, 249.3, 352.8, 8.0, 308.5, 353.9,
0172 235.9, 48.6, 435.7, 488.1, 406.5, 201.7, 198.5, 403.0, 203.4, 107.3, 144.1, 523.0, 85.6,
0173 547.4, 710.5, 254.4, 363.4, 19.3, 100.0, 121.2, 122.1, 557.0, 356.7, 171.3, 496.7, 538.8,
0174 564.1, 505.6, 381.3, 535.1, 311.8, 552.3, 570.9, 517.4, 549.5, 181.8, 540.6, 652.2, 514.2,
0175 197.6, 150.0, 409.9, 34.2, 485.3, 115.8, 163.7, 65.7, 0.0, 54.9, 366.4, 340.2, 36.9,
0176 353.0, 29.9, 21.3, 436.0, 150.3, 369.8, 313.9, 295.6, 0.0, 0.0, 182.2, 448.1, 519.3,
0177 223.5, 103.9, 38.7, 124.1, 430.9, 34.3, 108.3, 27.6, 14.8, 188.0, 0.0, 7.1, 241.2,
0178 447.0, 216.8, 96.4, 431.4, 412.0, 546.5, 232.0, 84.8, 394.7, 334.9, 77.1, 486.4, 96.4,
0179 251.9, 598.7, 339.4, 465.3, 542.1, 401.2, 0.1, 11.3, 8.5, 88.2, 495.1, 493.6, 412.8,
0180 415.9, 383.4, 405.9, 25.0, 0.2, 11.2, 38.4, 498.8, 467.6, 45.8, 213.1, 166.6, 263.8,
0181 6.5, 615.5, 374.4, 458.1, 0.0, 534.7, 132.6, 1.1, 234.2, 205.0, 229.8, 417.9, 297.5,
0182 443.3, 416.9, 339.1, 380.0, 111.1, 385.3, 0.0, 0.9, 0.0, 102.5, 130.3, 0.0, 111.8,
0183 71.2, 36.5, 1770.9, 10047.3, 12403.8, 10891.5, 3806.9, 5396.6, 8971.2, 5.8, 11172.5, 4364.3, 671.1,
0184 6.7, 1310.7, 7819.1, 8927.2, 7179.0, 8289.8, 8268.8, 4788.3, 5304.8, 5443.9, 668.1, 11286.1, 2569.6,
0185 7931.7, 10707.8, 7068.5, 10965.6, 10299.9, 5189.4, 112.7, 0.0, 0.0, 0.2, 0.1, 0.9, 0.6,
0186 8.3, 7.1, 3.9, 53.1, 63.3, 5.7, 85.2, 2.5, 50.1, 0.1, 43.3, 13.9, 132.2,
0187 70.6, 114.2, 209.8, 70.2, 140.9, 153.2, 254.1, 93.1, 380.7, 519.7, 427.4, 44.2, 214.1,
0188 347.0, 89.5, 8.6, 518.9, 390.4, 399.7, 435.9, 87.9, 390.3, 253.3, 12.3, 189.5, 0.0,
0189 0.0, 2.7, 26.6, 54.6, 4.9, 334.2, 56.8, 125.2, 226.6, 186.9, 330.1, 85.6, 99.2,
0190 17.4, 44.5, 15.0, 518.7, 460.9, 153.7, 111.6, 0.1, 0.0, 0.1, 2.1, 14.9, 54.8,
0191 10.8, 535.2, 523.1, 135.6, 522.3, 423.1, 110.5, 440.3, 487.8, 67.4, 523.4, 486.9, 353.8,
0192 561.3, 48.1, 105.1, 496.7, 427.1, 5.8, 22.9, 105.0, 123.1, 71.8, 149.4, 165.8, 178.1,
0193 70.2, 62.7, 102.8, 296.1, 74.6, 214.9, 107.2, 295.0, 254.9, 67.9, 136.2, 142.4, 20.3,
0194 92.0, 78.5, 304.2, 33.3, 309.3, 297.8, 10.4, 314.2, 122.3, 17.6, 166.5, 109.4, 112.6,
0195 146.8, 253.1, 189.9, 102.8, 101.3, 159.6, 136.2, 248.6, 324.2, 322.4, 74.5, 340.9, 14.8,
0196 13.8, 269.6, 83.3, 316.0, 339.0, 330.3, 111.0, 164.2, 16.3, 0.0, 3.5, 16.0, 149.6,
0197 344.8, 416.5, 432.5, 420.6, 285.3, 399.2, 457.1, 563.5, 390.3, 435.8, 150.7, 164.1, 134.0,
0198 389.8, 386.6, 73.4, 441.1, 437.7, 468.7, 237.8, 272.2, 161.1, 450.8, 300.1, 342.3, 502.0,
0199 463.7, 587.6, 0.0, 6.3, 501.9, 595.2, 335.4, 526.7, 494.5, 525.0, 420.7, 392.9, 521.9,
0200 48.8, 525.6, 472.8, 70.2, 546.6, 471.7, 519.0, 133.1, 501.9, 591.2, 552.9, 14.5, 279.2,
0201 49.6, 471.8, 255.3, 607.6, 213.0, 241.3, 162.3, 52.5, 50.7, 486.7, 555.8, 703.0, 712.6,
0202 743.8, 374.4, 680.8, 0.1, 4.2, 2.6, 22.7, 57.2, 36.6, 0.1, 1.2, 12.6, 0.3,
0203 0.1, 1.5, 13.6, 5.8, 9.4, 0.4, 63.2, 39.0, 63.2, 115.6, 7.2, 82.0, 78.1,
0204 181.2, 181.6, 95.0, 155.5, 264.8, 260.3, 354.6, 392.0, 37.6, 124.3, 509.7, 459.8, 332.3,
0205 365.0, 214.7, 25.7, 117.0, 123.7, 520.5, 317.6, 354.2, 305.0, 523.0, 141.9, 501.3, 259.6,
0206 554.5, 597.7, 575.8, 552.5, 32.3, 93.8, 204.3, 4.6, 342.6, 396.4, 402.1, 0.6, 493.0,
0207 189.5, 83.8, 436.6, 383.8, 490.6, 456.3, 558.2, 390.6, 45.0, 495.0, 172.5, 272.8, 8.9,
0208 496.8, 0.0, 407.9, 181.9, 67.5, 541.7, 518.8, 234.4, 54.7, 500.9, 500.6, 439.0, 358.9,
0209 449.0, 112.2, 527.2, 529.3, 487.2, 543.4, 118.7, 231.9, 0.1, 16.5, 46.4, 76.5, 13.9,
0210 400.0, 0.2, 0.4, 17.0, 0.0, 0.1, 0.4, 0.8, 0.1, 0.2, 2.4, 0.9, 1.3,
0211 229.7, 276.6, 602.5, 159.5, 613.2, 0.3, 793.8, 70.6, 405.4, 234.0, 95.4, 17.1, 100.1,
0212 486.3, 72.2, 342.3, 454.2, 115.0, 65.6, 529.3, 454.3, 308.6, 408.0, 424.5, 0.2, 0.0,
0213 0.2, 98.9, 156.3, 69.5, 116.9, 432.6, 452.6, 290.4, 441.8, 148.1, 506.5, 357.6, 39.0,
0214 78.0, 453.7, 473.8, 395.4, 514.5, 314.6, 458.4, 371.8, 214.6, 500.3, 322.9, 10.0, 411.6,
0215 754.0, 523.4, 572.2, 472.6, 535.4, 500.7, 117.1, 523.6, 446.2, 518.1, 80.0, 520.7, 527.0,
0216 162.5, 292.7, 521.7, 282.1, 107.8, 156.8, 535.8, 520.1, 13.9, 484.7, 16.0, 10.4, 248.2,
0217 361.3, 297.8, 332.0, 496.5, 528.9, 445.7, 510.0, 331.6, 436.2, 282.6, 188.7, 183.5, 495.2,
0218 481.5, 3.3, 429.6, 409.9, 0.0, 4.0, 70.2, 134.8, 398.5, 451.0, 0.2, 521.8, 497.0,
0219 413.0, 313.9, 337.4, 313.4, 268.8, 471.9, 328.5, 408.5, 462.0, 501.4, 79.7, 437.9, 537.3,
0220 196.0, 251.9, 0.0, 0.0, 0.2, 31.6, 533.4, 452.3, 546.8, 424.0, 54.2, 545.8, 550.8,
0221 412.6, 452.3, 127.9, 559.9, 449.9, 3.6, 0.0};
0222
0223
0224
0225 cout << "********************* treatment to obtaine massives of Days and DeliveredLumiForDays "
0226 "************************** "
0227 << endl;
0228 Float_t DeliveredLumiForDays[ALF];
0229 Int_t Days[ALF];
0230 for (int w9 = 0; w9 < ALF; w9++) {
0231 Days[w9] = 0.;
0232 DeliveredLumiForDays[w9] = 0.;
0233 }
0234 int w8 = 0;
0235 Days[w8] = DeltaDaysLumiTaken[0];
0236 DeliveredLumiForDays[w8] += DeliveredLumi[0];
0237
0238 for (int w9 = 1; w9 < ALF; w9++) {
0239 if (DeltaDaysLumiTaken[w9] != DeltaDaysLumiTaken[w9 - 1])
0240 w8++;
0241 Days[w8] = DeltaDaysLumiTaken[w9];
0242 if (DeltaDaysLumiTaken[w9] < 172)
0243 DeliveredLumiForDays[w8] += DeliveredLumi[w9];
0244 if (DeltaDaysLumiTaken[w9] > 171 && DeltaDaysLumiTaken[w9] < 210)
0245 DeliveredLumiForDays[w8] += DeliveredLumi[w9] / 1000.;
0246 if (DeltaDaysLumiTaken[w9] > 209)
0247 DeliveredLumiForDays[w8] += DeliveredLumi[w9];
0248
0249 }
0250
0251
0252
0253
0254
0255
0256
0257 Float_t LU2[NMAXDAYS];
0258 for (int w9 = 0; w9 < NMAXDAYS; w9++) {
0259 LU2[w9] = 0.;
0260 }
0261
0262 LU2[0] = DeliveredLumiForDays[0] / 1000.;
0263
0264
0265 int w8counts = 1;
0266 for (int jj = 1; jj < NMAXDAYS; jj++) {
0267 if (jj + 1 == Days[w8counts]) {
0268 LU2[jj] = LU2[jj - 1] + DeliveredLumiForDays[w8counts] / 1000.;
0269
0270 w8counts++;
0271 if (w8counts > w8)
0272 break;
0273 } else {
0274 LU2[jj] = LU2[jj - 1];
0275
0276 }
0277 }
0278
0279 cout << "before any correction: LU2[NMAXDAYS-1] = " << LU2[NMAXDAYS - 1] << endl;
0280
0281
0282
0283
0284 Float_t klumi2 = 1.0;
0285
0286 TH1F *lumiplot = new TH1F("lumiplot", "", NMAXDAYS + 1, 1., float(NMAXDAYS) + 1.);
0287 TH1F *lumiplot1 = new TH1F("lumiplot1", "", 172 + 1, 1., float(172) + 1.);
0288 TH1F *lumiplot2 = new TH1F("lumiplot2", "", 190 + 1, 1., float(190) + 1.);
0289 TH1F *lumiplot3 = new TH1F("lumiplot3", "", 195 + 1, 1., float(195) + 1.);
0290 for (int jj = 0; jj < NMAXDAYS; jj++) {
0291 if (jj < 172) {
0292 LU2[jj] *= klumi2;
0293 lumiplot1->Fill(float(jj + 1), LU2[jj]);
0294 }
0295 if (jj > 370) {
0296 lumiplot2->Fill(float(jj + 1 - 370), LU2[jj] - LU2[371]);
0297 }
0298 if (jj > 707) {
0299 lumiplot3->Fill(float(jj + 1 - 707), LU2[jj] - LU2[708]);
0300 }
0301 lumiplot->Fill(float(jj + 1), LU2[jj]);
0302
0303 }
0304 cout << "day 172(27.10.2016) LU2 = " << LU2[171] << endl;
0305
0306 cout << "since 15.05.2017 LU2 = " << LU2[555] - LU2[370] << endl;
0307 cout << "since 17.04.2018 LU2 = " << LU2[NMAXDAYS - 1] - LU2[706] << endl;
0308 cout << "after any correction: LU2[NMAXDAYS-1] = " << LU2[NMAXDAYS - 1] << endl;
0309
0310
0311
0312
0313
0314
0315
0316
0317 const Int_t NMAXP =
0318 134;
0319
0320 TString LASERrunsP[NMAXP] = {
0321 "272957", "273354", "273899", "274185", "274260", "274357", "274409", "274462", "274767", "275273", "275392",
0322 "275707", "275802", "275877", "276256", "276424", "276465", "276479", "276523", "276556", "276823", "277134",
0323 "277340", "278074", "278223", "278326", "278592", "278901", "279098", "279610", "279734", "279804", "280141",
0324 "280204", "280266", "280378", "282043", "282777", "282850", "282854", "282931", "283117", "283328", "283365",
0325 "283401", "283445", "283800", "283840", "283982", "285355",
0326
0327 "285701", "285878", "287385", "287739", "287765", "287766", "287866", "287867", "287868", "287869", "287870",
0328 "288174", "288180", "288181", "288182", "288988", "295678", "296810", "297084", "297526", "297573", "297632",
0329 "297692", "297850", "298104", "298295", "298794", "299461", "299971", "300181", "302109", "302256", "302412",
0330 "302431", "302434", "302587", "302614", "303249", "303853", "303928", "304232", "304425", "304601", "304633",
0331 "304682", "304748", "304860", "305127", "305217", "305472", "305702", "305874",
0332
0333 "315059", "315435", "315591", "315809", "316190", "316191", "316192", "316396", "316672", "316748", "316774",
0334 "316814", "316885", "316967", "317048", "317113", "317199", "317359", "317423", "317598", "317679", "317724",
0335 "317783", "317975", "320576", "321326", "321353", "321444", "321858", "321863", "323505", "324704"
0336
0337
0338
0339
0340
0341
0342
0343 };
0344
0345
0346
0347
0348
0349
0350 Float_t X0P[NMAXP] = {
0351 1., 4., 15., 20., 22., 25., 26., 28., 30., 38., 43., 46., 48., 50., 55., 57., 59., 59.,
0352 60., 62., 67., 74., 77., 85., 87., 90., 93., 99., 103., 109., 112., 114., 119., 120., 121., 123.,
0353 145., 154., 155., 155., 156., 157., 160., 161., 161., 162., 165., 166., 169., 190.,
0354
0355 197., 200., 277., 284., 284., 284., 287., 287., 287., 287., 287., 291., 291., 291., 291., 304., 387., 401.,
0356 404., 413., 414., 415., 416., 420., 422., 423., 428., 437., 444., 448., 479., 481., 485., 485., 485., 489.,
0357 490., 498., 504., 505., 511., 513., 515., 516., 518., 519., 521., 525., 527., 532., 535., 539.,
0358
0359 715., 721., 723., 728., 733., 733., 733., 737., 742., 743., 744., 745., 746., 748., 749., 750., 751., 755.,
0360 756., 760., 763., 764., 766., 770., 813., 829., 829., 832., 840., 840., 869., 889.
0361
0362
0363
0364 };
0365
0366
0367
0368
0369
0370
0371
0372
0373
0374
0375 const Int_t NMAXM = 121;
0376 TString LASERrunsM[NMAXM] = {
0377 "272967", "273211", "273460", "273562", "273897", "274355", "274408", "274765", "275017", "275272",
0378 "275357", "275706", "275874", "276252", "276420", "276470", "276478", "276522", "276532", "276860",
0379 "277133", "277174", "277869", "278071", "278221", "278325", "278590", "278898", "279097", "279701",
0380 "279802", "279909", "280138", "280202", "280207", "280264", "280376", "282042", "282243", "283325",
0381 "283364", "283400", "283443", "283799", "283837", "284735", "285353",
0382
0383 "285699", "286399", "287041", "288987", "290828", "293565", "295406", "295677", "295930", "295986",
0384 "296189", "296809", "297083", "297524", "297572", "297631", "297690", "297848", "298101", "298294",
0385 "298717", "298802", "299459", "299965", "300176", "302105", "302250", "302432", "302433", "302461",
0386 "302586", "302607", "303482", "303492", "303849", "303923", "304235", "304424", "304553", "304631",
0387 "304678", "304746", "304858", "305214", "305468", "305701", "305871",
0388
0389 "315056", "315586", "315808", "316118", "316193", "316194", "316253", "316527", "316528", "316772",
0390 "316812", "316884", "316966", "317044", "317109", "317198", "317310", "317356", "317424", "317597",
0391 "317723", "317782", "317967", "319607", "321030", "325265", "325634"
0392
0393
0394
0395
0396
0397
0398
0399
0400
0401
0402
0403
0404
0405
0406
0407
0408
0409
0410
0411
0412
0413
0414
0415 };
0416
0417
0418
0419
0420
0421
0422
0423 Float_t X0M[NMAXM] = {1., 3., 6., 9., 15., 25., 26., 30., 36., 38., 41., 46., 50., 55., 57., 59.,
0424 59., 60., 61., 68., 74., 75., 83., 85., 87., 90., 93., 99., 103., 111., 114., 116.,
0425 119., 120., 120., 121., 123., 145., 148., 159., 161., 161., 162., 165., 166., 181., 190.,
0426
0427 197., 207., 273., 304., 329., 364., 384., 387., 389., 390., 393., 401., 404., 413., 414., 415.,
0428 416., 420., 422., 423., 427., 428., 437., 444., 448., 479., 481., 485., 485., 486., 489., 490.,
0429 500., 500., 504., 505., 511., 513., 515., 516., 518., 519., 521., 527., 532., 535., 539.,
0430
0431 715., 723., 728., 732., 733., 733., 735., 739., 739., 744., 745., 746., 747., 748., 749., 750.,
0432 752., 754., 755., 759., 763., 765., 769., 795., 820., 899., 899.
0433
0434 };
0435
0436
0437
0438
0439
0440
0441
0442
0443
0444
0445
0446
0447
0448
0449
0450
0451
0452
0453
0454
0455
0456 const Int_t yearcase = 4;
0457
0458
0459
0460
0461 Int_t flaglumiuse = 2;
0462
0463
0464 cout << "!!********************* flaglumiuse = " << flaglumiuse << endl;
0465
0466
0467
0468
0469
0470 TCanvas *c1 = new TCanvas("c1", "Hcal4test", 10, 10, 900, 700);
0471
0472
0473 c1->Range(0, 0, 25, 18);
0474
0475 c1->SetFillColor(0);
0476 c1->Clear();
0477
0478 TLatex *t = new TLatex();
0479 t->SetTextFont(32);
0480
0481 t->SetTextSize(0.03);
0482 t->SetTextAlign(12);
0483
0484 t->DrawLatex(0.9,
0485 15.5,
0486 "Time dependence of full Amplitude averaged over events of each LASER run for some random HE channels");
0487
0488
0489
0490
0491
0492
0493
0494
0495
0496
0497 Float_t x1 = 0.00;
0498 Float_t x2 = 200.00;
0499 int Dref = 3;
0500 int DepthMax = 4;
0501 if (yearcase == 2)
0502 x2 = 200.00;
0503 if (yearcase == 3)
0504 x2 = 550.00;
0505 if (yearcase == 4)
0506 x2 = 200.00;
0507 if (yearcase == 5)
0508 x2 = 900.00;
0509
0510
0511
0512 Float_t y1 = 0.60;
0513 Float_t y2 = 1.01;
0514 Float_t y3 = 0.85;
0515 Float_t y4 = 1.01;
0516 Float_t y5 = 0.70;
0517 Float_t y6 = 1.01;
0518 if (yearcase == 2) {
0519 y1 = 0.60;
0520 y2 = 1.01;
0521 y3 = 0.85;
0522 y4 = 1.01;
0523 y5 = 0.70;
0524 y6 = 1.01;
0525 }
0526 if (yearcase == 3) {
0527 y1 = 0.40;
0528 y2 = 1.01;
0529 y3 = 0.75;
0530 y4 = 1.01;
0531 y5 = 0.60;
0532 y6 = 1.01;
0533 }
0534 if (yearcase == 4) {
0535 Dref = 4;
0536 DepthMax = 8;
0537 y1 = 0.60;
0538 y2 = 1.01;
0539 y3 = 0.85;
0540 y4 = 1.01;
0541 y5 = 0.70;
0542 y6 = 1.01;
0543 }
0544 if (yearcase > 4) {
0545 Dref = 4;
0546 DepthMax = 8;
0547 y1 = 0.30;
0548 y2 = 1.01;
0549 y3 = 0.65;
0550 y4 = 1.01;
0551 y5 = 0.50;
0552 y6 = 1.01;
0553 }
0554
0555
0556 Float_t xx1 = 0.12;
0557 Float_t xx2 = 0.31;
0558 Float_t yy1 = 0.17;
0559 Float_t yy2 = 0.38;
0560
0561
0562
0563
0564
0565 Int_t NMAXSIZE = NMAXM;
0566 if (NMAXP > NMAXM)
0567 NMAXSIZE = NMAXP;
0568
0569 Float_t X0[NMAXSIZE];
0570 TString LASERruns[NMAXSIZE];
0571 Double_t XX[NMAXSIZE];
0572 Double_t XXE[NMAXSIZE];
0573 const Int_t NDEP = 8;
0574 const Int_t NETA = 14;
0575 const Int_t NETA0 = 83;
0576 const Int_t NPHI0 = 73;
0577
0578 Double_t ZZM[NMAXSIZE][NDEP][NETA];
0579 Double_t ZZP[NMAXSIZE][NDEP][NETA];
0580 Double_t YY[NMAXSIZE];
0581 Double_t YYE[NMAXSIZE];
0582
0583 Double_t mapRADDAM_HE[NDEP][NETA0][NPHI0];
0584 Double_t mapRADDAM_HED2[NDEP][NETA0];
0585 Int_t mapRADDAM_HED20[NDEP][NETA0];
0586 Float_t gainRADDAM_HE[NDEP][NETA0][NPHI0];
0587 Double_t normM[NDEP][NETA0][NPHI0];
0588 Double_t normP[NDEP][NETA0][NPHI0];
0589
0590 Double_t THETA[13] = {17., 18., 19., 20., 21., 22., 23., 24., 25., 26., 27., 28., 29.};
0591 Double_t RRM[13];
0592 Double_t RRP[13];
0593 Double_t THETA7[11] = {18., 19., 20., 21., 22., 23., 24., 25., 26., 27., 28.};
0594 Double_t RRM7[11];
0595 Double_t RRP7[11];
0596 for (int ij = 0; ij < 13; ij++) {
0597 RRM[ij] = 0.;
0598 RRP[ij] = 0.;
0599 }
0600 for (int ij = 0; ij < 11; ij++) {
0601 RRM7[ij] = 0.;
0602 RRP7[ij] = 0.;
0603 }
0604
0605 TH1F *phidependenceM = new TH1F("phidependenceM", "", 72, 0., 72.);
0606 TH1F *phidependenceM0 = new TH1F("phidependenceM0", "", 72, 0., 72.);
0607 TH1F *phidependenceP = new TH1F("phidependenceP", "", 72, 0., 72.);
0608 TH1F *phidependenceP0 = new TH1F("phidependenceP0", "", 72, 0., 72.);
0609 TH1F *phidepdrawing30 = (TH1F *)phidependenceM->Clone("phidepdrawing30");
0610 TH1F *phidepdrawing31 = (TH1F *)phidependenceP->Clone("phidepdrawing31");
0611 TH1F *phiRRM = new TH1F("phiRRM", "", 100, 0., 5.);
0612 TH1F *phiRRP = new TH1F("phiRRP", "", 100, 0., 5.);
0613
0614 Int_t variant = 1;
0615
0616
0617
0618
0619 Int_t NMINmM = -1;
0620 Int_t NMAXmM = -1;
0621 Int_t NMAXaM = -1;
0622 Int_t NMINmP = -1;
0623 Int_t NMAXmP = -1;
0624 Int_t NMAXaP = -1;
0625 if (yearcase == 1) {
0626
0627 NMINmM = 0;
0628 NMAXmM = 46;
0629 NMAXaM = 47;
0630 NMINmP = 0;
0631 NMAXmP = 49;
0632 NMAXaP = 50;
0633 } else if (yearcase == 2) {
0634
0635
0636
0637 NMINmM = 53;
0638 NMAXmM = 93;
0639 NMAXaM = 41;
0640 NMINmP = 66;
0641 NMAXmP = 101;
0642 NMAXaP = 36;
0643 } else if (yearcase == 3) {
0644
0645 NMINmM = 0;
0646 NMAXmM = 93;
0647 NMAXaM = 94;
0648 NMINmP = 0;
0649 NMAXmP = 101;
0650 NMAXaP = 102;
0651 } else if (yearcase == 4) {
0652
0653
0654
0655 NMINmM = 94;
0656 NMAXmM = 120;
0657 NMAXaM = 27;
0658
0659 NMINmP = 102;
0660 NMAXmP = 133;
0661 NMAXaP = 32;
0662 } else if (yearcase == 5) {
0663
0664 NMINmM = 0;
0665 NMAXmM = -1;
0666 NMAXaM = -1;
0667 NMINmP = 0;
0668 NMAXmP = -1;
0669 NMAXaP = -1;
0670 }
0671
0672 if (NMINmM == -1 || NMAXmM == -1 || NMAXaM == -1 || NMINmP == -1 || NMAXmP == -1 || NMAXaP == -1)
0673 cout << " STOP !!!!!! NMINmM == -1 || NMAXmM == -1 || NMAXaM == -1 || NMINmP == -1 || NMAXmP == -1 || NMAXaP "
0674 "== -1"
0675 << endl;
0676
0677
0678
0679
0680
0681
0682
0683
0684 Int_t NMAX = 0;
0685
0686 for (int kdirection = 0; kdirection < 2; kdirection++) {
0687
0688
0689 if (kdirection == 0) {
0690 NMAX = NMAXaM;
0691
0692 for (int jjj = NMINmM; jjj <= NMAXmM; jjj++) {
0693 LASERruns[jjj - NMINmM] = LASERrunsM[jjj];
0694 X0[jjj - NMINmM] = X0M[jjj];
0695
0696 }
0697 }
0698
0699 else if (kdirection == 1) {
0700 NMAX = NMAXaP;
0701
0702 for (int jjj = NMINmP; jjj <= NMAXmP; jjj++) {
0703 LASERruns[jjj - NMINmP] = LASERrunsP[jjj];
0704 X0[jjj - NMINmP] = X0P[jjj];
0705
0706 }
0707 }
0708
0709
0710
0711
0712
0713 Int_t NP = NMAX;
0714 cout << "************************-=============================*************************=======__________++++++++++"
0715 "================================================= NP = "
0716 << NP << endl;
0717 cout << "***********************++++++++++================================================= kdirection = "
0718 << kdirection << endl;
0719
0720
0721
0722
0723
0724 Int_t flagcutamplitude = 0;
0725
0726 Double_t cutMinA = 10.;
0727 Double_t cutMaxA = 50000.;
0728 if (kdirection == 0) {
0729 cutMinA = 100.;
0730 cutMaxA = 3600.;
0731 }
0732
0733 Double_t kfactor = 2.6;
0734 Double_t k2factor = 1.5;
0735
0736 cout << "!!********************* flagcutamplitude= " << flagcutamplitude << endl;
0737 if (flagcutamplitude == 1)
0738 cout << " kfactor= " << kfactor << endl;
0739 if (flagcutamplitude == 2)
0740 cout << " k2factor= " << k2factor << endl;
0741
0742
0743
0744 Double_t cutMinM = 600.;
0745 if (kdirection == 0)
0746 cutMinM = 800.;
0747
0748
0749
0750
0751 Double_t CutrwithphiMin = 0.0;
0752 Double_t CutrwithphiMax = 6.0;
0753 if (kdirection == 0) {
0754 CutrwithphiMin = 0.0;
0755 CutrwithphiMax = 6.0;
0756 }
0757
0758 if (flagcutamplitude < 2) {
0759 CutrwithphiMin = 0.0;
0760 CutrwithphiMax = 6.0;
0761 if (kdirection == 0) {
0762 CutrwithphiMin = 0.0;
0763 CutrwithphiMax = 6.0;
0764 }
0765 }
0766
0767
0768
0769
0770
0771 Double_t CutrrwithphiMin = 0.0;
0772 Double_t CutrrwithphiMax = 4.0;
0773 if (kdirection == 0) {
0774 CutrrwithphiMin = 0.0;
0775 CutrrwithphiMax = 3.0;
0776 }
0777
0778
0779
0780
0781
0782
0783
0784
0785 Double_t CutrwithphiMinVar2 = 0.1;
0786 Double_t CutrwithphiMaxVar2 = 1.5;
0787 if (kdirection == 0) {
0788 CutrwithphiMinVar2 = 0.1;
0789 CutrwithphiMaxVar2 = 1.5;
0790 }
0791
0792
0793
0794
0795 for (int kk1 = 0; kk1 < NP; kk1++) {
0796 XX[kk1] = X0[kk1];
0797 int k10 = 0;
0798 int z1 = X0[k10];
0799 XX[kk1] -= z1;
0800 XXE[kk1] = 0.001 * XX[kk1];
0801 }
0802
0803
0804 Int_t flaggaincorrections = 0;
0805
0806 cout << " !!********************* flaggaincorrections = " << flaggaincorrections << endl;
0807
0808
0809
0810 if (flaglumiuse == 1) {
0811 for (int k1 = 0; k1 < NP; k1++) {
0812 int z1 = X0[k1];
0813 XX[k1] = LU1[z1];
0814 }
0815 }
0816
0817 if (flaglumiuse == 2) {
0818 for (int k1 = 0; k1 < NP; k1++) {
0819 int z1 = X0[k1];
0820 XX[k1] = LU2[z1];
0821 int k10 = 0;
0822 z1 = X0[k10];
0823 XX[k1] -= LU2[z1];
0824
0825 }
0826 }
0827
0828
0829 if (flaglumiuse > 0) {
0830 if (yearcase == 1)
0831 x2 = 45.0;
0832 if (yearcase == 2)
0833 x2 = 50.0;
0834 if (yearcase == 3)
0835 x2 = 95.0;
0836 if (yearcase == 4)
0837 x2 = 70.0;
0838 if (yearcase == 5)
0839 x2 = 165.0;
0840 }
0841
0842
0843 Int_t verbosity =
0844 0;
0845
0846
0847
0848
0849
0850
0851
0852
0853
0854
0855
0856
0857
0858
0859
0860
0861
0862
0863
0864
0865
0866
0867
0868
0869
0870
0871
0872
0873
0874
0875
0876
0877
0878
0879
0880
0881
0882
0883
0884
0885
0886
0887
0888
0889
0890
0891
0892
0893
0894
0895
0896
0897
0898
0899
0900
0901
0902
0903
0904
0905
0906
0907
0908
0909
0910
0911
0912
0913
0914
0915
0916
0917
0918
0919
0920
0921
0922
0923
0924
0925
0926
0927
0928
0929
0930
0931
0932
0933
0934
0935
0936
0937
0938
0939
0940
0941
0942
0943
0944
0945
0946
0947
0948
0949
0950
0951 for (int iii = 0; iii < DepthMax; iii++) {
0952 for (int jjj = 0; jjj <= 82; jjj++) {
0953 mapRADDAM_HED2[iii][jjj] = 0.;
0954 mapRADDAM_HED20[iii][jjj] = 0.;
0955 for (int kkk = 0; kkk <= 72; kkk++) {
0956 mapRADDAM_HE[iii][jjj][kkk] = 0.;
0957 normM[iii][jjj][kkk] = 0.;
0958 normP[iii][jjj][kkk] = 0.;
0959 gainRADDAM_HE[iii][jjj][kkk] = 0.;
0960 }
0961 }
0962 }
0963
0964
0965 TH1F *gainkoeff = new TH1F("gainkoeff", "", 100, 0.5, 1.5);
0966
0967
0968
0969
0970
0971
0972 TFile *hfile1 =
0973 new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_272847.root",
0974 "READ");
0975 TFile *hfile2 =
0976 new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_306249.root",
0977 "READ");
0978 TFile *hfile3 =
0979 new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_272847.root",
0980 "READ");
0981 TFile *hfile4 =
0982 new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_285349.root",
0983 "READ");
0984 TFile *hfile5 =
0985 new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_294054.root",
0986 "READ");
0987 TFile *hfile6 =
0988 new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_306249.root",
0989 "READ");
0990
0991
0992 TH2F *xwod1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl12_HE");
0993 TH2F *xwod2 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl12_HE");
0994 TH2F *xwod3 = (TH2F *)hfile1->Get("h_mapDepth3ADCAmpl12_HE");
0995 TH2F *xwo10 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
0996 TH2F *xwo20 = (TH2F *)hfile1->Get("h_mapDepth2_HE");
0997 TH2F *xwo30 = (TH2F *)hfile1->Get("h_mapDepth3_HE");
0998 TH2F *Xwof1 = (TH2F *)xwod1->Clone("Xwof1");
0999 TH2F *Xwof2 = (TH2F *)xwod2->Clone("Xwof2");
1000 TH2F *Xwof3 = (TH2F *)xwod3->Clone("Xwof3");
1001 Xwof1->Divide(xwod1, xwo10, 1, 1, "B");
1002 Xwof2->Divide(xwod2, xwo20, 1, 1, "B");
1003 Xwof3->Divide(xwod3, xwo30, 1, 1, "B");
1004
1005 TH2F *ywod1 = (TH2F *)hfile2->Get("h_mapDepth1ADCAmpl12_HE");
1006 TH2F *ywod2 = (TH2F *)hfile2->Get("h_mapDepth2ADCAmpl12_HE");
1007 TH2F *ywod3 = (TH2F *)hfile2->Get("h_mapDepth3ADCAmpl12_HE");
1008 TH2F *ywo10 = (TH2F *)hfile2->Get("h_mapDepth1_HE");
1009 TH2F *ywo20 = (TH2F *)hfile2->Get("h_mapDepth2_HE");
1010 TH2F *ywo30 = (TH2F *)hfile2->Get("h_mapDepth3_HE");
1011 TH2F *Ywof1 = (TH2F *)ywod1->Clone("Ywof1");
1012 TH2F *Ywof2 = (TH2F *)ywod2->Clone("Ywof2");
1013 TH2F *Ywof3 = (TH2F *)ywod3->Clone("Ywof3");
1014 Ywof1->Divide(ywod1, ywo10, 1, 1, "B");
1015 Ywof2->Divide(ywod2, ywo20, 1, 1, "B");
1016 Ywof3->Divide(ywod3, ywo30, 1, 1, "B");
1017
1018
1019 TH2F *xyod1 = (TH2F *)hfile3->Get("h_mapDepth1ADCAmpl12_HE");
1020 TH2F *xyod2 = (TH2F *)hfile3->Get("h_mapDepth2ADCAmpl12_HE");
1021 TH2F *xyod3 = (TH2F *)hfile3->Get("h_mapDepth3ADCAmpl12_HE");
1022 TH2F *xyo10 = (TH2F *)hfile3->Get("h_mapDepth1_HE");
1023 TH2F *xyo20 = (TH2F *)hfile3->Get("h_mapDepth2_HE");
1024 TH2F *xyo30 = (TH2F *)hfile3->Get("h_mapDepth3_HE");
1025 TH2F *Xyof1 = (TH2F *)xyod1->Clone("Xyof1");
1026 TH2F *Xyof2 = (TH2F *)xyod2->Clone("Xyof2");
1027 TH2F *Xyof3 = (TH2F *)xyod3->Clone("Xyof3");
1028 Xyof1->Divide(xyod1, xyo10, 1, 1, "B");
1029 Xyof2->Divide(xyod2, xyo20, 1, 1, "B");
1030 Xyof3->Divide(xyod3, xyo30, 1, 1, "B");
1031
1032 TH2F *yyod1 = (TH2F *)hfile4->Get("h_mapDepth1ADCAmpl12_HE");
1033 TH2F *yyod2 = (TH2F *)hfile4->Get("h_mapDepth2ADCAmpl12_HE");
1034 TH2F *yyod3 = (TH2F *)hfile4->Get("h_mapDepth3ADCAmpl12_HE");
1035 TH2F *yyo10 = (TH2F *)hfile4->Get("h_mapDepth1_HE");
1036 TH2F *yyo20 = (TH2F *)hfile4->Get("h_mapDepth2_HE");
1037 TH2F *yyo30 = (TH2F *)hfile4->Get("h_mapDepth3_HE");
1038 TH2F *Yyof1 = (TH2F *)yyod1->Clone("Yyof1");
1039 TH2F *Yyof2 = (TH2F *)yyod2->Clone("Yyof2");
1040 TH2F *Yyof3 = (TH2F *)yyod3->Clone("Yyof3");
1041 Yyof1->Divide(yyod1, yyo10, 1, 1, "B");
1042 Yyof2->Divide(yyod2, yyo20, 1, 1, "B");
1043 Yyof3->Divide(yyod3, yyo30, 1, 1, "B");
1044
1045
1046 TH2F *xzod1 = (TH2F *)hfile5->Get("h_mapDepth1ADCAmpl12_HE");
1047 TH2F *xzod2 = (TH2F *)hfile5->Get("h_mapDepth2ADCAmpl12_HE");
1048 TH2F *xzod3 = (TH2F *)hfile5->Get("h_mapDepth3ADCAmpl12_HE");
1049 TH2F *xzo10 = (TH2F *)hfile5->Get("h_mapDepth1_HE");
1050 TH2F *xzo20 = (TH2F *)hfile5->Get("h_mapDepth2_HE");
1051 TH2F *xzo30 = (TH2F *)hfile5->Get("h_mapDepth3_HE");
1052 TH2F *Xzof1 = (TH2F *)xzod1->Clone("Xzof1");
1053 TH2F *Xzof2 = (TH2F *)xzod2->Clone("Xzof2");
1054 TH2F *Xzof3 = (TH2F *)xzod3->Clone("Xzof3");
1055 Xzof1->Divide(xzod1, xzo10, 1, 1, "B");
1056 Xzof2->Divide(xzod2, xzo20, 1, 1, "B");
1057 Xzof3->Divide(xzod3, xzo30, 1, 1, "B");
1058
1059 TH2F *yzod1 = (TH2F *)hfile6->Get("h_mapDepth1ADCAmpl12_HE");
1060 TH2F *yzod2 = (TH2F *)hfile6->Get("h_mapDepth2ADCAmpl12_HE");
1061 TH2F *yzod3 = (TH2F *)hfile6->Get("h_mapDepth3ADCAmpl12_HE");
1062 TH2F *yzo10 = (TH2F *)hfile6->Get("h_mapDepth1_HE");
1063 TH2F *yzo20 = (TH2F *)hfile6->Get("h_mapDepth2_HE");
1064 TH2F *yzo30 = (TH2F *)hfile6->Get("h_mapDepth3_HE");
1065 TH2F *Yzof1 = (TH2F *)yzod1->Clone("Yzof1");
1066 TH2F *Yzof2 = (TH2F *)yzod2->Clone("Yzof2");
1067 TH2F *Yzof3 = (TH2F *)yzod3->Clone("Yzof3");
1068 Yzof1->Divide(yzod1, yzo10, 1, 1, "B");
1069 Yzof2->Divide(yzod2, yzo20, 1, 1, "B");
1070 Yzof3->Divide(yzod3, yzo30, 1, 1, "B");
1071
1072
1073
1074
1075 int gx = Xwof1->GetXaxis()->GetNbins();
1076 int gy = Xwof1->GetYaxis()->GetNbins();
1077 for (int jjj = 0; jjj < gx; jjj++) {
1078 for (int kkk = 0; kkk < gy; kkk++) {
1079 if (Xwof1->GetBinContent(jjj, kkk) != 0.) {
1080 gainRADDAM_HE[1][jjj][kkk + 1] = (Ywof1->GetBinContent(jjj, kkk) / Xwof1->GetBinContent(jjj, kkk));
1081 }
1082 if (Xwof2->GetBinContent(jjj, kkk) != 0.) {
1083 gainRADDAM_HE[2][jjj][kkk + 1] = (Ywof2->GetBinContent(jjj, kkk) / Xwof2->GetBinContent(jjj, kkk));
1084 }
1085 if (Xwof3->GetBinContent(jjj, kkk) != 0.) {
1086 gainRADDAM_HE[3][jjj][kkk + 1] = (Ywof3->GetBinContent(jjj, kkk) / Xwof3->GetBinContent(jjj, kkk));
1087 }
1088 }
1089 }
1090
1091 for (int jjj = 0; jjj < gx; jjj++) {
1092 for (int kkk = 0; kkk < gy; kkk++) {
1093 if (Xyof1->GetBinContent(jjj, kkk) != 0. &&
1094 fabs((Yyof1->GetBinContent(jjj, kkk) / Xyof1->GetBinContent(jjj, kkk)) - 1.) <
1095 fabs(gainRADDAM_HE[1][jjj][kkk + 1] - 1.)) {
1096 gainRADDAM_HE[1][jjj][kkk + 1] = (Yyof1->GetBinContent(jjj, kkk) / Xyof1->GetBinContent(jjj, kkk));
1097 }
1098 if (Xyof2->GetBinContent(jjj, kkk) != 0. &&
1099 fabs((Yyof2->GetBinContent(jjj, kkk) / Xyof2->GetBinContent(jjj, kkk)) - 1.) <
1100 fabs(gainRADDAM_HE[1][jjj][kkk + 1] - 1.)) {
1101 gainRADDAM_HE[2][jjj][kkk + 1] = (Yyof2->GetBinContent(jjj, kkk) / Xyof2->GetBinContent(jjj, kkk));
1102 }
1103 if (Xyof3->GetBinContent(jjj, kkk) != 0. &&
1104 fabs((Yyof3->GetBinContent(jjj, kkk) / Xyof3->GetBinContent(jjj, kkk)) - 1.) <
1105 fabs(gainRADDAM_HE[1][jjj][kkk + 1] - 1.)) {
1106 gainRADDAM_HE[3][jjj][kkk + 1] = (Yyof3->GetBinContent(jjj, kkk) / Xyof3->GetBinContent(jjj, kkk));
1107 }
1108 }
1109 }
1110
1111 for (int jjj = 0; jjj < gx; jjj++) {
1112 for (int kkk = 0; kkk < gy; kkk++) {
1113 if (Xzof1->GetBinContent(jjj, kkk) != 0. &&
1114 fabs((Yzof1->GetBinContent(jjj, kkk) / Xzof1->GetBinContent(jjj, kkk)) - 1.) <
1115 fabs(gainRADDAM_HE[1][jjj][kkk + 1] - 1.)) {
1116 gainRADDAM_HE[1][jjj][kkk + 1] = (Yzof1->GetBinContent(jjj, kkk) / Xzof1->GetBinContent(jjj, kkk));
1117 }
1118 if (Xzof2->GetBinContent(jjj, kkk) != 0. &&
1119 fabs((Yzof2->GetBinContent(jjj, kkk) / Xzof2->GetBinContent(jjj, kkk)) - 1.) <
1120 fabs(gainRADDAM_HE[1][jjj][kkk + 1] - 1.)) {
1121 gainRADDAM_HE[2][jjj][kkk + 1] = (Yzof2->GetBinContent(jjj, kkk) / Xzof2->GetBinContent(jjj, kkk));
1122 }
1123 if (Xzof3->GetBinContent(jjj, kkk) != 0. &&
1124 fabs((Yzof3->GetBinContent(jjj, kkk) / Xzof3->GetBinContent(jjj, kkk)) - 1.) <
1125 fabs(gainRADDAM_HE[1][jjj][kkk + 1] - 1.)) {
1126 gainRADDAM_HE[3][jjj][kkk + 1] = (Yzof3->GetBinContent(jjj, kkk) / Xzof3->GetBinContent(jjj, kkk));
1127 }
1128 }
1129 }
1130
1131
1132
1133 for (int iii = 1; iii < DepthMax; iii++) {
1134 for (int jjj = 1; jjj <= 82; jjj++) {
1135 for (int kkk = 1; kkk <= 72; kkk++) {
1136 if (fabs(gainRADDAM_HE[iii][jjj][kkk] - 1.) > 0.5)
1137 gainRADDAM_HE[iii][jjj][kkk] = 1.;
1138 if (flaggaincorrections == 0)
1139 gainRADDAM_HE[iii][jjj][kkk] = 1.;
1140 gainkoeff->Fill(gainRADDAM_HE[iii][jjj][kkk]);
1141 }
1142 }
1143 }
1144
1145
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159 TString fname = "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LASER_";
1160
1161
1162
1163 TH1F *Ceff1D = new TH1F("Ceff1D", "", 200, 0., 5000.);
1164 TH1F *Seff1D = new TH1F("Seff1D", "", 200, 0., 1000.);
1165 TH1F *Ceff1M = new TH1F("Ceff1M", "", 200, 0., 1300.);
1166 TH1F *Ceff1R = new TH1F("Ceff1R", "", 100, 0., 1000.);
1167
1168 TH1F *Meff1D = new TH1F("Meff1D", "", 100, 0., 4000.);
1169 TH1F *Meff2D = new TH1F("Meff2D", "", 100, 0., 4000.);
1170 TH1F *Meff3D = new TH1F("Meff3D", "", 100, 0., 4000.);
1171 TH1F *Meff4D = new TH1F("Meff4D", "", 100, 0., 4000.);
1172 TH1F *Meff5D = new TH1F("Meff5D", "", 100, 0., 4000.);
1173 TH1F *Meff6D = new TH1F("Meff6D", "", 100, 0., 4000.);
1174 TH1F *Meff7D = new TH1F("Meff7D", "", 100, 0., 4000.);
1175
1176
1177
1178
1179
1180 TH1F *rwithphi = new TH1F("rwithphi", "", 100, -1., 9.);
1181
1182
1183 TH1F *rrwithphiM = new TH1F("rrwithphiM", "", 100, -1., 4.);
1184 TH1F *rrwithphiP = new TH1F("rrwithphiP", "", 100, -1., 4.);
1185
1186 Int_t kmm = 0;
1187 Int_t kpp = 0;
1188
1189 Int_t mk1f = 0;
1190 Int_t pk1f = 0;
1191
1192
1193
1194
1195
1196
1197
1198 Int_t RUN3 = 0;
1199 Dref = 3;
1200 DepthMax = 4;
1201 for (int k1 = 0; k1 < NP; k1++) {
1202
1203
1204
1205 TString epsName = fname + LASERruns[k1] + ".root";
1206 TFile *hfile = new TFile(epsName, "READ");
1207
1208
1209 if (LASERruns[k1] > "305874")
1210 RUN3 = 1;
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220 TH2F *twod1 = (TH2F *)hfile->Get("h_mapDepth1ADCAmpl12_HE");
1221 TH2F *twod2 = (TH2F *)hfile->Get("h_mapDepth2ADCAmpl12_HE");
1222 TH2F *twod3 = (TH2F *)hfile->Get("h_mapDepth3ADCAmpl12_HE");
1223 TH2F *twod4 = (TH2F *)hfile->Get("h_mapDepth4ADCAmpl12_HE");
1224 TH2F *twod5 = (TH2F *)hfile->Get("h_mapDepth5ADCAmpl12_HE");
1225 TH2F *twod6 = (TH2F *)hfile->Get("h_mapDepth6ADCAmpl12_HE");
1226 TH2F *twod7 = (TH2F *)hfile->Get("h_mapDepth7ADCAmpl12_HE");
1227
1228 TH2F *t0010 = (TH2F *)hfile->Get("h_mapDepth1_HE");
1229 TH2F *t0020 = (TH2F *)hfile->Get("h_mapDepth2_HE");
1230 TH2F *t0030 = (TH2F *)hfile->Get("h_mapDepth3_HE");
1231 TH2F *t0040 = (TH2F *)hfile->Get("h_mapDepth4_HE");
1232 TH2F *t0050 = (TH2F *)hfile->Get("h_mapDepth5_HE");
1233 TH2F *t0060 = (TH2F *)hfile->Get("h_mapDepth6_HE");
1234 TH2F *t0070 = (TH2F *)hfile->Get("h_mapDepth7_HE");
1235 if (RUN3 = 1) {
1236 Dref = 4;
1237 DepthMax = 8;
1238 }
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271 TH2F *Ceff1 = (TH2F *)twod1->Clone("Ceff1");
1272 TH2F *Ceff2 = (TH2F *)twod2->Clone("Ceff2");
1273 TH2F *Ceff3 = (TH2F *)twod3->Clone("Ceff3");
1274 TH2F *Ceff4 = (TH2F *)twod4->Clone("Ceff4");
1275 TH2F *Ceff5 = (TH2F *)twod5->Clone("Ceff5");
1276 TH2F *Ceff6 = (TH2F *)twod6->Clone("Ceff6");
1277 TH2F *Ceff7 = (TH2F *)twod7->Clone("Ceff7");
1278
1279
1280
1281
1282
1283
1284
1285
1286 int kscale = 1;
1287 if (RUN3 == 1)
1288 kscale = 200;
1289 Ceff1->Divide(twod1, t0010, 1, kscale, "B");
1290 Ceff2->Divide(twod2, t0020, 1, kscale, "B");
1291 Ceff3->Divide(twod3, t0030, 1, kscale, "B");
1292 Ceff4->Divide(twod4, t0040, 1, kscale, "B");
1293 Ceff5->Divide(twod5, t0050, 1, kscale, "B");
1294 Ceff6->Divide(twod6, t0060, 1, kscale, "B");
1295 Ceff7->Divide(twod7, t0070, 1, kscale, "B");
1296
1297
1298
1299 int nx = Ceff1->GetXaxis()->GetNbins();
1300 int ny = Ceff1->GetYaxis()->GetNbins();
1301
1302 if (verbosity == 2) {
1303 cout << "************* nx = " << nx << " ny = " << ny << endl;
1304 for (int jjj = 0; jjj <= nx; jjj++) {
1305 for (int kkk = 0; kkk <= ny; kkk++) {
1306 if (Ceff1->GetBinContent(jjj, kkk) != 0.)
1307 cout << "*111111************ jjj = " << jjj << "************* kkk = " << kkk << endl;
1308 if (Ceff2->GetBinContent(jjj, kkk) != 0.)
1309 cout << "*******222****** jjj = " << jjj << "************* kkk = " << kkk << endl;
1310 if (Ceff3->GetBinContent(jjj, kkk) != 0.)
1311 cout << "**********3*** jjj = " << jjj << "************* kkk = " << kkk << endl;
1312 if (Ceff4->GetBinContent(jjj, kkk) != 0.)
1313 cout << "*144444************ jjj = " << jjj << "************* kkk = " << kkk << endl;
1314 if (Ceff5->GetBinContent(jjj, kkk) != 0.)
1315 cout << "*******555****** jjj = " << jjj << "************* kkk = " << kkk << endl;
1316 if (Ceff6->GetBinContent(jjj, kkk) != 0.)
1317 cout << "**********6*** jjj = " << jjj << "************* kkk = " << kkk << endl;
1318 if (Ceff7->GetBinContent(jjj, kkk) != 0.)
1319 cout << "**********7*** jjj = " << jjj << "************* kkk = " << kkk << endl;
1320 }
1321 }
1322 }
1323
1324
1325
1326 for (int jjj = 0; jjj < nx; jjj++) {
1327 for (int kkk = 0; kkk < ny; kkk++) {
1328
1329
1330 if (Ceff1->GetBinContent(jjj, kkk) != 0.) {
1331 Ceff1D->Fill(Ceff1->GetBinContent(jjj, kkk));
1332 Seff1D->Fill(Ceff1->GetBinContent(jjj, kkk));
1333 }
1334 if (Ceff2->GetBinContent(jjj, kkk) != 0.) {
1335 Ceff1D->Fill(Ceff2->GetBinContent(jjj, kkk));
1336 Seff1D->Fill(Ceff2->GetBinContent(jjj, kkk));
1337 }
1338 if (Ceff3->GetBinContent(jjj, kkk) != 0.) {
1339 Ceff1D->Fill(Ceff3->GetBinContent(jjj, kkk));
1340 Seff1D->Fill(Ceff3->GetBinContent(jjj, kkk));
1341 }
1342 if (Ceff4->GetBinContent(jjj, kkk) != 0.) {
1343 Ceff1D->Fill(Ceff4->GetBinContent(jjj, kkk));
1344 Seff1D->Fill(Ceff4->GetBinContent(jjj, kkk));
1345 }
1346 if (Ceff5->GetBinContent(jjj, kkk) != 0.) {
1347 Ceff1D->Fill(Ceff5->GetBinContent(jjj, kkk));
1348 Seff1D->Fill(Ceff5->GetBinContent(jjj, kkk));
1349 }
1350 if (Ceff6->GetBinContent(jjj, kkk) != 0.) {
1351 Ceff1D->Fill(Ceff6->GetBinContent(jjj, kkk));
1352 Seff1D->Fill(Ceff6->GetBinContent(jjj, kkk));
1353 }
1354 if (Ceff7->GetBinContent(jjj, kkk) != 0.) {
1355 Ceff1D->Fill(Ceff7->GetBinContent(jjj, kkk));
1356 Seff1D->Fill(Ceff7->GetBinContent(jjj, kkk));
1357 }
1358
1359 if (Ceff1->GetBinContent(jjj, kkk) > cutMinA) {
1360 Meff1D->Fill(Ceff1->GetBinContent(jjj, kkk));
1361 }
1362 if (Ceff2->GetBinContent(jjj, kkk) > cutMinA) {
1363 Meff2D->Fill(Ceff2->GetBinContent(jjj, kkk));
1364 }
1365 if (Ceff3->GetBinContent(jjj, kkk) > cutMinA) {
1366 Meff3D->Fill(Ceff3->GetBinContent(jjj, kkk));
1367 }
1368 if (Ceff4->GetBinContent(jjj, kkk) > cutMinA) {
1369 Meff4D->Fill(Ceff4->GetBinContent(jjj, kkk));
1370 }
1371 if (Ceff5->GetBinContent(jjj, kkk) > cutMinA) {
1372 Meff5D->Fill(Ceff5->GetBinContent(jjj, kkk));
1373 }
1374 if (Ceff6->GetBinContent(jjj, kkk) > cutMinA) {
1375 Meff6D->Fill(Ceff6->GetBinContent(jjj, kkk));
1376 }
1377 if (Ceff7->GetBinContent(jjj, kkk) > cutMinA) {
1378 Meff7D->Fill(Ceff7->GetBinContent(jjj, kkk));
1379 }
1380 }
1381 }
1382
1383 Ceff1M->Fill(Meff1D->GetMean());
1384 Ceff1M->Fill(Meff2D->GetMean());
1385 Ceff1M->Fill(Meff3D->GetMean());
1386 Ceff1M->Fill(Meff4D->GetMean());
1387 Ceff1M->Fill(Meff5D->GetMean());
1388 Ceff1M->Fill(Meff6D->GetMean());
1389 Ceff1M->Fill(Meff7D->GetMean());
1390
1391 if (verbosity == 9 && (Meff1D->GetMean() < cutMinM || Meff2D->GetMean() < cutMinM || Meff3D->GetMean() < cutMinM))
1392 cout << " kdirection = " << kdirection << "****** run out Mean = " << k1
1393 << " Meff1D->GetMean() = " << Meff1D->GetMean() << " cutMinM = " << cutMinM << endl;
1394 Ceff1R->Fill(Meff1D->GetRMS());
1395 Ceff1R->Fill(Meff2D->GetRMS());
1396 Ceff1R->Fill(Meff3D->GetRMS());
1397 Ceff1R->Fill(Meff4D->GetRMS());
1398 Ceff1R->Fill(Meff5D->GetRMS());
1399 Ceff1R->Fill(Meff6D->GetRMS());
1400 Ceff1R->Fill(Meff7D->GetRMS());
1401
1402
1403
1404
1405
1406
1407
1408
1409
1410
1411
1412
1413
1414
1415
1416
1417
1418 if (flagcutamplitude == 0) {
1419 for (int jjj = 0; jjj < nx; jjj++) {
1420 for (int kkk = 0; kkk < ny; kkk++) {
1421 if (Ceff1->GetBinContent(jjj, kkk) > cutMinA && Ceff1->GetBinContent(jjj, kkk) < cutMaxA)
1422 mapRADDAM_HE[1][jjj][kkk + 1] = Ceff1->GetBinContent(jjj, kkk);
1423 if (Ceff2->GetBinContent(jjj, kkk) > cutMinA && Ceff2->GetBinContent(jjj, kkk) < cutMaxA)
1424 mapRADDAM_HE[2][jjj][kkk + 1] = Ceff2->GetBinContent(jjj, kkk);
1425 if (Ceff3->GetBinContent(jjj, kkk) > cutMinA && Ceff3->GetBinContent(jjj, kkk) < cutMaxA)
1426 mapRADDAM_HE[3][jjj][kkk + 1] = Ceff3->GetBinContent(jjj, kkk);
1427 if (Ceff4->GetBinContent(jjj, kkk) > cutMinA && Ceff4->GetBinContent(jjj, kkk) < cutMaxA)
1428 mapRADDAM_HE[4][jjj][kkk + 1] = Ceff4->GetBinContent(jjj, kkk);
1429 if (Ceff5->GetBinContent(jjj, kkk) > cutMinA && Ceff5->GetBinContent(jjj, kkk) < cutMaxA)
1430 mapRADDAM_HE[5][jjj][kkk + 1] = Ceff5->GetBinContent(jjj, kkk);
1431 if (Ceff6->GetBinContent(jjj, kkk) > cutMinA && Ceff6->GetBinContent(jjj, kkk) < cutMaxA)
1432 mapRADDAM_HE[6][jjj][kkk + 1] = Ceff6->GetBinContent(jjj, kkk);
1433 if (Ceff7->GetBinContent(jjj, kkk) > cutMinA && Ceff7->GetBinContent(jjj, kkk) < cutMaxA)
1434 mapRADDAM_HE[7][jjj][kkk + 1] = Ceff7->GetBinContent(jjj, kkk);
1435 }
1436 }
1437 }
1438
1439
1440
1441
1442
1443
1444
1445 else if (flagcutamplitude == 1) {
1446 for (int jjj = 0; jjj < nx; jjj++) {
1447 for (int kkk = 0; kkk < ny; kkk++) {
1448 if (Ceff1->GetBinContent(jjj, kkk) > cutMinA &&
1449 Ceff1->GetBinContent(jjj, kkk) < kfactor * Meff1D->GetMean())
1450 mapRADDAM_HE[1][jjj][kkk + 1] = Ceff1->GetBinContent(jjj, kkk);
1451 if (Ceff2->GetBinContent(jjj, kkk) > cutMinA &&
1452 Ceff2->GetBinContent(jjj, kkk) < kfactor * Meff2D->GetMean())
1453 mapRADDAM_HE[2][jjj][kkk + 1] = Ceff2->GetBinContent(jjj, kkk);
1454 if (Ceff3->GetBinContent(jjj, kkk) > cutMinA &&
1455 Ceff3->GetBinContent(jjj, kkk) < kfactor * Meff3D->GetMean())
1456 mapRADDAM_HE[3][jjj][kkk + 1] = Ceff3->GetBinContent(jjj, kkk);
1457 if (Ceff4->GetBinContent(jjj, kkk) > cutMinA &&
1458 Ceff4->GetBinContent(jjj, kkk) < kfactor * Meff1D->GetMean())
1459 mapRADDAM_HE[4][jjj][kkk + 1] = Ceff4->GetBinContent(jjj, kkk);
1460 if (Ceff5->GetBinContent(jjj, kkk) > cutMinA &&
1461 Ceff5->GetBinContent(jjj, kkk) < kfactor * Meff2D->GetMean())
1462 mapRADDAM_HE[5][jjj][kkk + 1] = Ceff5->GetBinContent(jjj, kkk);
1463 if (Ceff6->GetBinContent(jjj, kkk) > cutMinA &&
1464 Ceff6->GetBinContent(jjj, kkk) < kfactor * Meff3D->GetMean())
1465 mapRADDAM_HE[6][jjj][kkk + 1] = Ceff6->GetBinContent(jjj, kkk);
1466 if (Ceff7->GetBinContent(jjj, kkk) > cutMinA &&
1467 Ceff7->GetBinContent(jjj, kkk) < kfactor * Meff3D->GetMean())
1468 mapRADDAM_HE[7][jjj][kkk + 1] = Ceff7->GetBinContent(jjj, kkk);
1469 }
1470 }
1471 }
1472 else if (flagcutamplitude == 2) {
1473 for (int jjj = 0; jjj < nx; jjj++) {
1474 for (int kkk = 0; kkk < ny; kkk++) {
1475 if (Ceff1->GetBinContent(jjj, kkk) > cutMinA &&
1476 Ceff1->GetBinContent(jjj, kkk) < (Meff1D->GetMean() + k2factor * Meff1D->GetRMS()))
1477 mapRADDAM_HE[1][jjj][kkk + 1] = Ceff1->GetBinContent(jjj, kkk);
1478 if (Ceff2->GetBinContent(jjj, kkk) > cutMinA &&
1479 Ceff2->GetBinContent(jjj, kkk) < (Meff2D->GetMean() + k2factor * Meff2D->GetRMS()))
1480 mapRADDAM_HE[2][jjj][kkk + 1] = Ceff2->GetBinContent(jjj, kkk);
1481 if (Ceff3->GetBinContent(jjj, kkk) > cutMinA &&
1482 Ceff3->GetBinContent(jjj, kkk) < (Meff3D->GetMean() + k2factor * Meff3D->GetRMS()))
1483 mapRADDAM_HE[3][jjj][kkk + 1] = Ceff3->GetBinContent(jjj, kkk);
1484 if (Ceff4->GetBinContent(jjj, kkk) > cutMinA &&
1485 Ceff4->GetBinContent(jjj, kkk) < (Meff4D->GetMean() + k2factor * Meff4D->GetRMS()))
1486 mapRADDAM_HE[4][jjj][kkk + 1] = Ceff4->GetBinContent(jjj, kkk);
1487 if (Ceff5->GetBinContent(jjj, kkk) > cutMinA &&
1488 Ceff5->GetBinContent(jjj, kkk) < (Meff5D->GetMean() + k2factor * Meff5D->GetRMS()))
1489 mapRADDAM_HE[5][jjj][kkk + 1] = Ceff5->GetBinContent(jjj, kkk);
1490 if (Ceff6->GetBinContent(jjj, kkk) > cutMinA &&
1491 Ceff6->GetBinContent(jjj, kkk) < (Meff6D->GetMean() + k2factor * Meff6D->GetRMS()))
1492 mapRADDAM_HE[6][jjj][kkk + 1] = Ceff6->GetBinContent(jjj, kkk);
1493 if (Ceff7->GetBinContent(jjj, kkk) > cutMinA &&
1494 Ceff7->GetBinContent(jjj, kkk) < (Meff7D->GetMean() + k2factor * Meff7D->GetRMS()))
1495 mapRADDAM_HE[7][jjj][kkk + 1] = Ceff7->GetBinContent(jjj, kkk);
1496 }
1497 }
1498 }
1499
1500
1501 for (int iii = 1; iii < DepthMax; iii++) {
1502 for (int jjj = 1; jjj <= 82; jjj++) {
1503 if ((kdirection == 0 && jjj > 41) || (kdirection == 1 && jjj < 42)) {
1504 for (int kkk = 1; kkk <= 72; kkk++) {
1505 mapRADDAM_HE[iii][jjj][kkk] = 0.;
1506 }
1507 }
1508 }
1509 }
1510
1511
1512 if (k1 > 36 && k1 < 45) {
1513 for (int iii = 1; iii < DepthMax; iii++) {
1514 for (int jjj = 1; jjj <= 82; jjj++) {
1515 for (int kkk = 1; kkk <= 72; kkk++) {
1516 if ((iii == 1) &&
1517 (jjj == 59 || jjj == 60 || jjj == 63 || jjj == 65 || jjj == 66 || jjj == 67 || jjj == 68) &&
1518 (kkk == 23 || kkk == 24)) {
1519 mapRADDAM_HE[iii][jjj][kkk] = 0.;
1520 }
1521 }
1522 }
1523 }
1524 }
1525
1526
1527
1528
1529
1530
1531 for (int iii = 1; iii < DepthMax; iii++) {
1532 for (int jjj = 1; jjj <= 82; jjj++) {
1533 for (int kkk = 1; kkk <= 72; kkk++) {
1534 if (mapRADDAM_HE[iii][jjj][kkk] != 0 && gainRADDAM_HE[iii][jjj][kkk] != 0)
1535 mapRADDAM_HE[iii][jjj][kkk] /= gainRADDAM_HE[iii][jjj][kkk];
1536 }
1537 }
1538 }
1539
1540
1541
1542
1543
1544 if (k1 == mk1f && kmm == 0 && kdirection == 0) {
1545 kmm = 1;
1546 cout << "* to be printed once only CHECK IT please ************ kmm = " << kmm << endl;
1547 cout << "NEGATIVE direction refRun !!!! read LASER run with k1= " << k1 << " run: " << LASERruns[k1] << endl;
1548 for (int iii = 1; iii < DepthMax; iii++) {
1549 for (int jjj = 1; jjj <= 82; jjj++) {
1550 for (int kkk = 1; kkk <= 72; kkk++) {
1551
1552 normM[iii][jjj][kkk] = mapRADDAM_HE[iii][jjj][kkk];
1553
1554
1555 }
1556 }
1557 }
1558
1559 }
1560
1561 if (k1 == pk1f && kpp == 0 && kdirection == 1) {
1562 kpp = 1;
1563 cout << "* to be printed once only CHECK IT please ************ kpp = " << kpp << endl;
1564 cout << "POSITIVE direction refRun !!!! read LASER run with k1= " << k1 << " run: " << LASERruns[k1] << endl;
1565 for (int iii = 1; iii < DepthMax; iii++) {
1566 for (int jjj = 1; jjj <= 82; jjj++) {
1567 for (int kkk = 1; kkk <= 72; kkk++) {
1568
1569 normP[iii][jjj][kkk] = mapRADDAM_HE[iii][jjj][kkk];
1570
1571
1572 }
1573 }
1574 }
1575
1576 }
1577
1578
1579
1580
1581
1582 for (int iii = 1; iii < DepthMax; iii++) {
1583 for (int jjj = 1; jjj <= 82; jjj++) {
1584 for (int kkk = 1; kkk <= 72; kkk++) {
1585 if (jjj < 42 && kdirection == 0) {
1586 if (normM[iii][jjj][kkk] != 0) {
1587 mapRADDAM_HE[iii][jjj][kkk] /= normM[iii][jjj][kkk];
1588 } else {
1589 mapRADDAM_HE[iii][jjj][kkk] = 0.;
1590 }
1591 }
1592 if (jjj > 41 && kdirection == 1) {
1593 if (normP[iii][jjj][kkk] != 0) {
1594 mapRADDAM_HE[iii][jjj][kkk] /= normP[iii][jjj][kkk];
1595 } else {
1596 mapRADDAM_HE[iii][jjj][kkk] = 0.;
1597 }
1598 }
1599 }
1600 }
1601 }
1602
1603
1604
1605
1606
1607
1608 if (variant == 2) {
1609 Double_t normchannelM = 0;
1610 Double_t normchannelP = 0;
1611
1612 if (kdirection == 0) {
1613 Int_t mcounter = 0;
1614 for (int kkk = 1; kkk <= 72; kkk++) {
1615 if (mapRADDAM_HE[Dref][26][kkk] > 0.) {
1616 normchannelM += mapRADDAM_HE[Dref][26][kkk];
1617 mcounter++;
1618 }
1619 }
1620 if (mcounter != 0)
1621 normchannelM /= mcounter;
1622 } else if (kdirection == 1) {
1623 Int_t pcounter = 0;
1624 for (int kkk = 1; kkk <= 72; kkk++) {
1625 if (mapRADDAM_HE[Dref][57][kkk] > 0.) {
1626 normchannelP += mapRADDAM_HE[Dref][57][kkk];
1627 pcounter++;
1628 }
1629 }
1630 if (pcounter != 0)
1631 normchannelP /= pcounter;
1632 }
1633
1634
1635
1636
1637 for (int iii = 1; iii < DepthMax; iii++) {
1638 for (int jjj = 1; jjj <= 82; jjj++) {
1639 for (int kkk = 1; kkk <= 72; kkk++) {
1640 if (mapRADDAM_HE[iii][jjj][kkk] > 0. && normchannelM != 0 && jjj < 42 && kdirection == 0) {
1641 mapRADDAM_HE[iii][jjj][kkk] /= normchannelM;
1642 }
1643 if (mapRADDAM_HE[iii][jjj][kkk] > 0. && normchannelP != 0 && jjj > 41 && kdirection == 1) {
1644 mapRADDAM_HE[iii][jjj][kkk] /= normchannelP;
1645 }
1646 }
1647 }
1648 }
1649
1650 }
1651
1652
1653
1654
1655
1656
1657 for (int iii = 1; iii < DepthMax; iii++) {
1658 for (int jjj = 1; jjj <= 82; jjj++) {
1659 for (int kkk = 1; kkk <= 72; kkk++) {
1660 if (mapRADDAM_HE[iii][jjj][kkk] > 0) {
1661
1662
1663 rwithphi->Fill(mapRADDAM_HE[iii][jjj][kkk]);
1664
1665
1666 if ((mapRADDAM_HE[iii][jjj][kkk] > CutrwithphiMin && mapRADDAM_HE[iii][jjj][kkk] < CutrwithphiMax &&
1667 variant != 2) ||
1668 (mapRADDAM_HE[iii][jjj][kkk] > CutrwithphiMinVar2 &&
1669 mapRADDAM_HE[iii][jjj][kkk] < CutrwithphiMaxVar2 && variant == 2)) {
1670
1671 Double_t RRrrRr = 0.;
1672 if (variant != 2) {
1673 if (jjj < 42 && kdirection == 0) {
1674 if (mapRADDAM_HE[Dref][26][kkk] != 0) {
1675 RRrrRr = mapRADDAM_HE[iii][jjj][kkk] / mapRADDAM_HE[Dref][26][kkk];
1676 rrwithphiM->Fill(RRrrRr);
1677 }
1678 } else if (jjj > 41 && kdirection == 1) {
1679 if (mapRADDAM_HE[Dref][57][kkk] != 0) {
1680 RRrrRr = mapRADDAM_HE[iii][jjj][kkk] / mapRADDAM_HE[Dref][57][kkk];
1681 rrwithphiP->Fill(RRrrRr);
1682 }
1683 }
1684 } else if (variant == 2) {
1685 if (jjj < 42 && kdirection == 0) {
1686 RRrrRr = mapRADDAM_HE[iii][jjj][kkk];
1687 rrwithphiM->Fill(RRrrRr);
1688 } else if (jjj > 41 && kdirection == 1) {
1689 RRrrRr = mapRADDAM_HE[iii][jjj][kkk];
1690 rrwithphiP->Fill(RRrrRr);
1691 }
1692
1693 }
1694
1695
1696 if ((RRrrRr > CutrrwithphiMin && RRrrRr < CutrrwithphiMax && variant != 2) || variant == 2) {
1697
1698 mapRADDAM_HED2[iii][jjj] += mapRADDAM_HE[iii][jjj][kkk];
1699
1700 mapRADDAM_HED20[iii][jjj]++;
1701
1702 if (verbosity == 3 && iii == Dref)
1703 cout << " depth = " << iii << " k1 = " << k1 << " eta = " << jjj
1704 << " mapRADDAM_HED2[iii][jjj] = " << mapRADDAM_HED2[iii][jjj]
1705 << " Nphi = " << mapRADDAM_HED20[iii][jjj]
1706 << " mapRADDAM_HE[iii][jjj][kkk = " << mapRADDAM_HE[iii][jjj][kkk] << endl;
1707
1708
1709 if (jjj < 42 && kdirection == 0) {
1710 phidependenceM->Fill(float(kkk - 1), mapRADDAM_HE[iii][jjj][kkk]);
1711 phidependenceM0->Fill(float(kkk - 1), 1.);
1712 phiRRM->Fill(mapRADDAM_HE[iii][jjj][kkk]);
1713 } else if (jjj > 41 && kdirection == 1) {
1714 phidependenceP->Fill(float(kkk - 1), mapRADDAM_HE[iii][jjj][kkk]);
1715 phidependenceP0->Fill(float(kkk - 1), 1.);
1716 phiRRP->Fill(mapRADDAM_HE[iii][jjj][kkk]);
1717 }
1718
1719 }
1720 }
1721
1722 }
1723 }
1724 }
1725 }
1726
1727
1728
1729
1730
1731
1732
1733
1734
1735
1736
1737
1738
1739
1740
1741
1742
1743
1744
1745
1746
1747
1748
1749
1750
1751 for (int iii = 1; iii < DepthMax; iii++) {
1752 for (int jjj = 1; jjj <= 82; jjj++) {
1753 if (mapRADDAM_HED20[iii][jjj] != 0) {
1754 if (verbosity == 1)
1755 cout << "****aver per N-phi_sectors********* depth = " << iii << " RunIndex = " << k1 << " eta = " << jjj
1756 << " Nphi = " << mapRADDAM_HED20[iii][jjj] << endl;
1757 mapRADDAM_HED2[iii][jjj] /= mapRADDAM_HED20[iii][jjj];
1758 if (verbosity == 3 && iii == Dref)
1759 cout << "****aver per N-phi_sectors********* depth = " << iii << " RunIndex = " << k1 << " eta = " << jjj
1760 << " Nphi = " << mapRADDAM_HED20[iii][jjj] << endl;
1761
1762
1763
1764
1765 if (jjj < 42) {
1766 if (kdirection == 0) {
1767 int jj = jjj - 13;
1768 if (jj < 0 || jj > 13)
1769 cout << "***ZZM ERROR *** depth = " << iii << " eta = " << jjj << " new-eta = " << jj << endl;
1770 ZZM[k1][iii][jj] = mapRADDAM_HED2[iii][jjj];
1771 if (verbosity == 1)
1772 cout << "****aver per N-phi_sectors********* ZZM[k1][iii][jj] = " << ZZM[k1][iii][jj] << endl;
1773 }
1774 } else {
1775 if (kdirection == 1) {
1776 int jj = jjj - 13;
1777 jj -= 44;
1778 if (jj < 0 || jj > 13)
1779 cout << "***ZZP ERROR *** depth = " << iii << " eta = " << jjj << " new-eta = " << jj << endl;
1780 ZZP[k1][iii][jj] = mapRADDAM_HED2[iii][jjj];
1781 if (verbosity == 1)
1782 cout << "****aver per N-phi_sectors****** ZZP[k1][iii][jj] = " << ZZP[k1][iii][jj] << endl;
1783 }
1784 }
1785
1786 }
1787
1788 }
1789 }
1790
1791
1792
1793
1794 for (int iii = 0; iii < DepthMax; iii++) {
1795 for (int jjj = 0; jjj <= 82; jjj++) {
1796 mapRADDAM_HED2[iii][jjj] = 0.;
1797 mapRADDAM_HED20[iii][jjj] = 0.;
1798 for (int kkk = 0; kkk <= 72; kkk++) {
1799 mapRADDAM_HE[iii][jjj][kkk] = 0.;
1800 }
1801 }
1802 }
1803
1804 hfile->Close();
1805
1806
1807
1808
1809
1810
1811
1812
1813
1814
1815
1816
1817
1818
1819
1820
1821
1822 }
1823 cout << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! DONE: last LASER run READ " << endl;
1824
1825
1826
1827
1828
1829
1830
1831
1832
1833 if (variant != 2) {
1834
1835
1836 for (int k1 = 0; k1 < NP; k1++) {
1837 if (verbosity == 1) {
1838 cout << " ********************************************************************time *************point = "
1839 << k1 << endl;
1840 }
1841 for (int iii = 1; iii < DepthMax; iii++) {
1842 for (int jj = 0; jj <= 13; jj++) {
1843 if (verbosity == 1) {
1844 cout << " **** DepthMax = " << DepthMax << " **** iii = " << iii << " **** jj = " << jj << endl;
1845 }
1846
1847
1848 if (kdirection == 1) {
1849 if (iii == Dref && jj == 0) {
1850 } else {
1851 if (ZZP[k1][Dref][0] != 0) {
1852 ZZP[k1][iii][jj] /= ZZP[k1][Dref][0];
1853 }
1854 if (verbosity == 1 && ZZP[k1][iii][jj] != 0) {
1855 cout << "Normalized P Positive Direction ***** depth = " << iii << " eta = " << jj
1856 << " ZZP[k1][iii][jj] = " << ZZP[k1][iii][jj] << endl;
1857 }
1858 if (verbosity == 3 && ZZP[k1][iii][jj] != 0 && iii == Dref) {
1859 cout << "Normalized P Positive Direction ***** depth = " << iii << " eta = " << jj
1860 << " ZZP[k1][iii][jj] = " << ZZP[k1][iii][jj] << " ZZP[k1][Dref][0] = " << ZZP[k1][Dref][0]
1861 << endl;
1862 }
1863
1864 }
1865 }
1866
1867 else if (kdirection == 0) {
1868
1869 if (iii == Dref && jj == 13) {
1870 } else {
1871 if (ZZM[k1][Dref][13] != 0) {
1872 ZZM[k1][iii][jj] /= ZZM[k1][Dref][13];
1873 }
1874
1875
1876 }
1877
1878 }
1879
1880 }
1881 }
1882 }
1883
1884 }
1885 cout << "DONE variant = " << variant << endl;
1886
1887
1888
1889
1890
1891
1892
1893
1894
1895
1896
1897
1898
1899
1900
1901
1902
1903
1904
1905
1906
1907
1908
1909
1910
1911
1912
1913
1914
1915
1916
1917
1918
1919 int k1last = NP - 1;
1920
1921
1922 if (kdirection == 1 && ZZP[k1last][2][13] < 0.1)
1923 k1last = NP - 2;
1924
1925
1926 if (kdirection == 0) {
1927
1928 cout << "******************** M L1 k1last = " << k1last
1929 << " ZZM[k1last][2][10] = " << ZZM[k1last][2][10] << endl;
1930 cout << "********************************************************************************************************"
1931 "********************************************************* "
1932 << endl;
1933
1934 cout << "********************* M L7 k1last = " << k1last
1935 << " ZZM[k1last][5][1] = " << ZZM[k1last][5][1] << endl;
1936 cout << "********************* M L7 k1last = " << k1last
1937 << " ZZM[k1last][5][2] = " << ZZM[k1last][5][2] << endl;
1938 cout << "********************* M L7 k1last = " << k1last
1939 << " ZZM[k1last][5][3] = " << ZZM[k1last][5][3] << endl;
1940 cout << "********************* M L7 k1last = " << k1last
1941 << " ZZM[k1last][4][4] = " << ZZM[k1last][4][4] << endl;
1942 cout << "********************* M L7 k1last = " << k1last
1943 << " ZZM[k1last][4][5] = " << ZZM[k1last][4][5] << endl;
1944 cout << "********************* M L7 k1last = " << k1last
1945 << " ZZM[k1last][4][6] = " << ZZM[k1last][4][6] << endl;
1946 cout << "********************* M L7 k1last = " << k1last
1947 << " ZZM[k1last][4][7] = " << ZZM[k1last][4][7] << endl;
1948 cout << "********************* M L7 k1last = " << k1last
1949 << " ZZM[k1last][4][8] = " << ZZM[k1last][4][8] << endl;
1950 cout << "********************* M L7 k1last = " << k1last
1951 << " ZZM[k1last][4][9] = " << ZZM[k1last][4][9] << endl;
1952 cout << "********************* M L7 k1last = " << k1last
1953 << " ZZM[k1last][4][10] = " << ZZM[k1last][4][10] << endl;
1954 cout << "********************* M L7 k1last = " << k1last
1955 << " ZZM[k1last][4][11] = " << ZZM[k1last][4][11] << endl;
1956 cout << "********************* M L7 k1last = " << k1last
1957 << " ZZM[k1last][3][12] = " << ZZM[k1last][3][12] << endl;
1958 cout << "****************************************************************************** " << endl;
1959 }
1960
1961
1962
1963
1964
1965
1966
1967
1968
1969
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981 if (kdirection == 0)
1982 for (int jj = 0; jj < 13; jj++) {
1983 if (LASERruns[k1last] > "305874") {
1984 RRM[12 - jj] = ZZM[k1last][2][jj];
1985 } else {
1986 RRM[12 - jj] = ZZM[k1last][1][jj];
1987 }
1988 }
1989 if (kdirection == 1)
1990 for (int jj = 1; jj < 14; jj++) {
1991 if (LASERruns[k1last] > "305874") {
1992 RRP[jj - 1] = ZZP[k1last][2][jj];
1993 } else {
1994 RRP[jj - 1] = ZZP[k1last][1][jj];
1995 }
1996 }
1997
1998
1999 if (kdirection == 0)
2000 for (int iii = 2; iii < DepthMax; iii++) {
2001 for (int jj = 1; jj < 12; jj++) {
2002 if (LASERruns[k1last] > "305874") {
2003 if ((iii == 3 && 13 - jj == 1) || (iii == 4 && 13 - jj > 1 && 13 - jj < 10) ||
2004 (iii == 5 && 13 - jj > 9 && 13 - jj < 13)) {
2005 RRM7[11 - jj] = ZZM[k1last][iii][jj];
2006
2007 }
2008 } else {
2009 if ((iii == 2 && jj > 2 && jj < 12) || (iii == 3 && jj > 0 && jj < 3))
2010 RRM7[11 - jj] = ZZM[k1last][iii][jj];
2011 }
2012 }
2013 }
2014
2015 if (kdirection == 1)
2016 for (int iii = 2; iii < DepthMax; iii++) {
2017 for (int jj = 1; jj < 13; jj++) {
2018 if (LASERruns[k1last] > "305874") {
2019 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13)) {
2020 RRP7[jj - 2] = ZZP[k1last][iii][jj];
2021 }
2022 } else {
2023 if ((iii == 2 && jj > 1 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
2024 RRP7[jj - 2] = ZZP[k1last][iii][jj];
2025 }
2026 }
2027 }
2028
2029
2030
2031
2032
2033
2034
2035
2036 cout << "**** DRAWING **** **** **** **** Draw Dependencies: *** " << endl;
2037
2038
2039 Int_t Method = NMethodsStart;
2040
2041
2042 while (Method < NMethodsToUse + 1) {
2043 if ((kdirection == 1 && Method != 9 && Method != 19 && Method != 24 && Method != 26 && Method != 28 &&
2044 Method != 32 && Method != 1 && Method != 2 && Method != 3 && Method != 4 && Method != 5 && Method != 6 &&
2045 Method != 7 && Method != 8 && Method != 30) ||
2046 (kdirection == 0 && Method != 0 && Method != 10 && Method != 22 && Method != 27 && Method != 29 &&
2047 Method != 33 && Method != 11 && Method != 12 && Method != 13 && Method != 14 && Method != 15 &&
2048 Method != 16 && Method != 17 && Method != 18 && Method != 31 && Method != 20 && Method != 21 &&
2049 Method != 23 && Method != 25 && Method != 34 && Method != 35 && Method != 36)) {
2050 cout << " **** START Method " << Method << endl;
2051
2052 if (Method == 0) {
2053 cout << "****Draw 0 Dependencies: P: draw <Aijk>_ev.run for all runs *** " << endl;
2054 gPad->SetGridy();
2055 gPad->SetGridx();
2056 gPad->SetLogy();
2057 gStyle->SetOptStat(101110);
2058
2059 Ceff1D->SetMarkerStyle(20);
2060 Ceff1D->SetMarkerSize(0.4);
2061 Ceff1D->GetYaxis()->SetLabelSize(0.04);
2062 Ceff1D->SetTitle("P: all runs & channels");
2063 Ceff1D->SetXTitle("<Aijk>");
2064 Ceff1D->SetYTitle("N");
2065 Ceff1D->SetMarkerColor(4);
2066 Ceff1D->SetLineColor(4);
2067
2068 Ceff1D->Draw("PL");
2069 }
2070
2071
2072 if (Method == 10) {
2073 cout << "****Draw 10 Dependencies: P: draw <Aijk>_ev.run for all runs *** " << endl;
2074 gPad->SetGridy();
2075 gPad->SetGridx();
2076 gPad->SetLogy();
2077 gStyle->SetOptStat(101110);
2078
2079 Seff1D->SetMarkerStyle(20);
2080 Seff1D->SetMarkerSize(0.4);
2081 Seff1D->GetYaxis()->SetLabelSize(0.04);
2082 Seff1D->SetTitle("P: all runs & channels");
2083 Seff1D->SetXTitle("<Aijk>");
2084 Seff1D->SetYTitle("N");
2085 Seff1D->SetMarkerColor(4);
2086 Seff1D->SetLineColor(4);
2087
2088 Seff1D->Draw("PL");
2089 }
2090
2091
2092 if (Method == 9) {
2093 cout << "****Draw 9 Dependencies: M: draw <Aijk>_ev.run for all runs *** " << endl;
2094 gPad->SetGridy();
2095 gPad->SetGridx();
2096 gPad->SetLogy();
2097 gStyle->SetOptStat(101110);
2098 c1->SetFillColor(0);
2099
2100 Ceff1D->SetMarkerStyle(20);
2101 Ceff1D->SetMarkerSize(0.4);
2102 Ceff1D->GetYaxis()->SetLabelSize(0.04);
2103 Ceff1D->SetTitle("M: all runs & channels");
2104 Ceff1D->SetXTitle("<Aijk>");
2105 Ceff1D->SetYTitle("N");
2106 Ceff1D->SetMarkerColor(4);
2107 Ceff1D->SetLineColor(4);
2108
2109 Ceff1D->Draw("PL");
2110 }
2111
2112
2113 if (Method == 19) {
2114 cout << "****Draw 19 Dependencies: M: draw <Aijk>_ev.run for all runs *** " << endl;
2115 gPad->SetGridy();
2116 gPad->SetGridx();
2117 gPad->SetLogy();
2118 gStyle->SetOptStat(101110);
2119 c1->SetFillColor(0);
2120
2121 Seff1D->SetMarkerStyle(20);
2122 Seff1D->SetMarkerSize(0.4);
2123 Seff1D->GetYaxis()->SetLabelSize(0.04);
2124 Seff1D->SetTitle("M: all runs & channels");
2125 Seff1D->SetXTitle("<Aijk>");
2126 Seff1D->SetYTitle("N");
2127 Seff1D->SetMarkerColor(4);
2128 Seff1D->SetLineColor(4);
2129
2130 Seff1D->Draw("PL");
2131 }
2132
2133
2134 if (Method == 20) {
2135 cout << "****Draw 20 Dependencies: P: draw mean <Aijk>_ev.run for all runs Depth3 *** " << endl;
2136 gPad->SetGridy();
2137 gPad->SetGridx();
2138 gPad->SetLogy();
2139 gStyle->SetOptStat(101110);
2140 Meff1D->SetMarkerStyle(20);
2141 Meff1D->SetMarkerSize(0.4);
2142 Meff1D->GetYaxis()->SetLabelSize(0.04);
2143 Meff1D->SetXTitle("for mean <Aijk>");
2144 Meff1D->SetYTitle("N");
2145 Meff1D->SetMarkerColor(4);
2146 Meff1D->SetLineColor(4);
2147 Meff1D->Draw("Error");
2148 Meff2D->SetMarkerColor(2);
2149 Meff2D->SetLineColor(2);
2150 Meff2D->Draw("ErrorSame");
2151 Meff3D->SetMarkerColor(3);
2152 Meff3D->SetLineColor(3);
2153 Meff3D->Draw("ErrorSame");
2154
2155 }
2156
2157
2158
2159
2160
2161
2162
2163
2164
2165
2166 if (Method == 1) {
2167 for (int k1 = 0; k1 < NP; k1++) {
2168 YY[k1] = 0.;
2169 }
2170 c1->Range(0, 0, 25, 18);
2171 c1->SetFillColor(0);
2172 c1->Clear();
2173 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
2174
2175 pad1->Draw();
2176 pad1->cd();
2177 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
2178
2179 pad1->DrawFrame(x1, y1, x2, y2);
2180
2181 pad1->GetFrame()->SetFillColor(0);
2182 pad1->Draw();
2183 gPad->SetGridy();
2184 gPad->SetGridx();
2185 cout << "****Draw 1 Dependencies: M r vs t for Depth1 and some eta 21 22 23 25 27 28 *** " << endl;
2186
2187
2188
2189 for (int k1 = 0; k1 < NP; k1++) {
2190 YY[k1] = ZZM[k1][1][1];
2191 }
2192 gr1 = new TGraph(NP, XX, YY);
2193 gr1->SetLineColor(1);
2194 gr1->SetMarkerColor(1);
2195 gr1->SetMarkerStyle(20);
2196 gr1->SetMarkerSize(1.2);
2197 gr1->SetTitle("ieta 28");
2198 gr1->SetFillStyle(0);
2199 gr1->Draw("P");
2200 for (int k1 = 0; k1 < NP; k1++) {
2201 YY[k1] = ZZM[k1][1][2];
2202 }
2203 gr2 = new TGraph(NP, XX, YY);
2204 gr2->SetLineColor(2);
2205 gr2->SetMarkerColor(2);
2206 gr2->SetMarkerStyle(20);
2207 gr2->SetMarkerSize(1.2);
2208 gr2->SetTitle("ieta 27");
2209 gr2->SetFillStyle(0);
2210 gr2->Draw("P");
2211 for (int k1 = 0; k1 < NP; k1++) {
2212 YY[k1] = ZZM[k1][1][4];
2213 }
2214 gr3 = new TGraph(NP, XX, YY);
2215 gr3->SetLineColor(3);
2216 gr3->SetMarkerColor(3);
2217 gr3->SetMarkerStyle(20);
2218 gr3->SetMarkerSize(1.2);
2219 gr3->SetTitle("ieta 25");
2220 gr3->SetFillStyle(0);
2221 gr3->Draw("P");
2222 for (int k1 = 0; k1 < NP; k1++) {
2223 YY[k1] = ZZM[k1][1][6];
2224 }
2225 gr4 = new TGraph(NP, XX, YY);
2226 gr4->SetLineColor(4);
2227 gr4->SetMarkerColor(4);
2228 gr4->SetMarkerStyle(20);
2229 gr4->SetMarkerSize(1.2);
2230 gr4->SetTitle("ieta 23");
2231 gr4->SetFillStyle(0);
2232 gr4->Draw("P");
2233 for (int k1 = 0; k1 < NP; k1++) {
2234 YY[k1] = ZZM[k1][1][7];
2235 }
2236 gr5 = new TGraph(NP, XX, YY);
2237 gr5->SetLineColor(5);
2238 gr5->SetMarkerColor(5);
2239 gr5->SetMarkerStyle(20);
2240 gr5->SetMarkerSize(1.2);
2241 gr5->SetTitle("ieta 22");
2242 gr5->SetFillStyle(0);
2243 gr5->Draw("P");
2244 for (int k1 = 0; k1 < NP; k1++) {
2245 YY[k1] = ZZM[k1][1][8];
2246 }
2247 gr6 = new TGraph(NP, XX, YY);
2248 gr6->SetLineColor(6);
2249 gr6->SetMarkerColor(6);
2250 gr6->SetMarkerStyle(20);
2251 gr6->SetMarkerSize(1.2);
2252 gr6->SetTitle("ieta 21");
2253 gr6->SetFillStyle(0);
2254 gr6->Draw("P");
2255
2256 TLegend leg(xx1, yy1, xx2, yy2, " HEM Depth1");
2257 leg.AddEntry(gr1);
2258 leg.AddEntry(gr2);
2259 leg.AddEntry(gr3);
2260 leg.AddEntry(gr4);
2261 leg.AddEntry(gr5);
2262 leg.AddEntry(gr6);
2263 leg.SetMargin(0.3);
2264 leg.SetBorderSize(1);
2265 leg.DrawClone("Same");
2266
2267 }
2268
2269 if (Method == 2) {
2270 for (int k1 = 0; k1 < NP; k1++) {
2271 YY[k1] = 0.;
2272 }
2273 c1->Range(0, 0, 25, 18);
2274 c1->SetFillColor(0);
2275 c1->Clear();
2276 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
2277 pad1->Draw();
2278 pad1->cd();
2279 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
2280
2281 pad1->DrawFrame(x1, y1, x2, y2);
2282
2283 pad1->GetFrame()->SetFillColor(0);
2284 pad1->Draw();
2285 gPad->SetGridy();
2286 gPad->SetGridx();
2287 cout << "****Draw 2 Dependencies: M r vs t for Depth1 and some eta 17 18 19 20 24 26 29 *** " << endl;
2288
2289
2290
2291 for (int k1 = 0; k1 < NP; k1++) {
2292 YY[k1] = ZZM[k1][1][0];
2293
2294
2295
2296
2297 }
2298 gr1 = new TGraph(NP, XX, YY);
2299 gr1->SetLineColor(1);
2300 gr1->SetMarkerColor(1);
2301 gr1->SetMarkerStyle(20);
2302 gr1->SetMarkerSize(1.2);
2303 gr1->SetTitle("ieta 29");
2304 gr1->SetFillStyle(0);
2305 gr1->Draw("P");
2306 for (int k1 = 0; k1 < NP; k1++) {
2307 YY[k1] = ZZM[k1][1][3];
2308 }
2309 gr2 = new TGraph(NP, XX, YY);
2310 gr2->SetLineColor(2);
2311 gr2->SetMarkerColor(2);
2312 gr2->SetMarkerStyle(20);
2313 gr2->SetMarkerSize(1.2);
2314 gr2->SetTitle("ieta 26");
2315 gr2->SetFillStyle(0);
2316 gr2->Draw("P");
2317 for (int k1 = 0; k1 < NP; k1++) {
2318 YY[k1] = ZZM[k1][1][5];
2319 }
2320 gr3 = new TGraph(NP, XX, YY);
2321 gr3->SetLineColor(3);
2322 gr3->SetMarkerColor(3);
2323 gr3->SetMarkerStyle(20);
2324 gr3->SetMarkerSize(1.2);
2325 gr3->SetTitle("ieta 24");
2326 gr3->SetFillStyle(0);
2327 gr3->Draw("P");
2328
2329
2330 for (int k1 = 0; k1 < NP; k1++) {
2331 YY[k1] = ZZM[k1][1][10];
2332 }
2333 gr4 = new TGraph(NP, XX, YY);
2334 gr4->SetLineColor(5);
2335 gr4->SetMarkerColor(5);
2336 gr4->SetMarkerStyle(20);
2337 gr4->SetMarkerSize(1.2);
2338 gr4->SetTitle("ieta 19");
2339 gr4->SetFillStyle(0);
2340 gr4->Draw("P");
2341
2342
2343 for (int k1 = 0; k1 < NP; k1++) {
2344 YY[k1] = ZZM[k1][1][12];
2345 }
2346 gr5 = new TGraph(NP, XX, YY);
2347 gr5->SetLineColor(7);
2348 gr5->SetMarkerColor(7);
2349 gr5->SetMarkerStyle(20);
2350 gr5->SetMarkerSize(1.2);
2351 gr5->SetTitle("ieta 17");
2352 gr5->SetFillStyle(0);
2353 gr5->Draw("P");
2354
2355 TLegend leg(xx1, yy1, xx2, yy2, " HEM Depth1");
2356 leg.AddEntry(gr1);
2357 leg.AddEntry(gr2);
2358 leg.AddEntry(gr3);
2359 leg.AddEntry(gr4);
2360 leg.AddEntry(gr5);
2361 leg.SetMargin(0.3);
2362 leg.SetBorderSize(1);
2363 leg.DrawClone("Same");
2364
2365 }
2366
2367
2368 if (Method == 3) {
2369 for (int k1 = 0; k1 < NP; k1++) {
2370 YY[k1] = 0.;
2371 }
2372 c1->Range(0, 0, 25, 18);
2373 c1->SetFillColor(0);
2374 c1->Clear();
2375 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
2376 pad1->Draw();
2377 pad1->cd();
2378 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
2379
2380 pad1->DrawFrame(x1, y3, x2, y4);
2381
2382 pad1->GetFrame()->SetFillColor(0);
2383 pad1->Draw();
2384 gPad->SetGridy();
2385 gPad->SetGridx();
2386 cout << "****Draw 3 Dependencies: M r vs t for Depth2 and some eta 21 22 23 25 27 28 *** " << endl;
2387
2388
2389
2390 for (int k1 = 0; k1 < NP; k1++) {
2391 YY[k1] = ZZM[k1][2][4];
2392 }
2393 gr1 = new TGraph(NP, XX, YY);
2394 gr1->SetLineColor(1);
2395 gr1->SetMarkerColor(1);
2396 gr1->SetMarkerStyle(20);
2397 gr1->SetMarkerSize(1.2);
2398 gr1->SetTitle("ieta 25");
2399 gr1->SetFillStyle(0);
2400 gr1->Draw("P");
2401 for (int k1 = 0; k1 < NP; k1++) {
2402 YY[k1] = ZZM[k1][2][6];
2403 }
2404 gr2 = new TGraph(NP, XX, YY);
2405 gr2->SetLineColor(2);
2406 gr2->SetMarkerColor(2);
2407 gr2->SetMarkerStyle(20);
2408 gr2->SetMarkerSize(1.2);
2409 gr2->SetTitle("ieta 23");
2410 gr2->SetFillStyle(0);
2411 gr2->Draw("P");
2412 for (int k1 = 0; k1 < NP; k1++) {
2413 YY[k1] = ZZM[k1][2][7];
2414 }
2415 gr3 = new TGraph(NP, XX, YY);
2416 gr3->SetLineColor(3);
2417 gr3->SetMarkerColor(3);
2418 gr3->SetMarkerStyle(20);
2419 gr3->SetMarkerSize(1.2);
2420 gr3->SetTitle("ieta 22");
2421 gr3->SetFillStyle(0);
2422 gr3->Draw("P");
2423 for (int k1 = 0; k1 < NP; k1++) {
2424 YY[k1] = ZZM[k1][2][8];
2425 }
2426 gr4 = new TGraph(NP, XX, YY);
2427 gr4->SetLineColor(4);
2428 gr4->SetMarkerColor(4);
2429 gr4->SetMarkerStyle(20);
2430 gr4->SetMarkerSize(1.2);
2431 gr4->SetTitle("ieta 21");
2432 gr4->SetFillStyle(0);
2433 gr4->Draw("P");
2434
2435 TLegend leg(xx1, yy1, xx2, yy2, " HEM Depth2");
2436 leg.AddEntry(gr1);
2437 leg.AddEntry(gr2);
2438 leg.AddEntry(gr3);
2439 leg.AddEntry(gr4);
2440 leg.SetMargin(0.3);
2441 leg.SetBorderSize(1);
2442 leg.DrawClone("Same");
2443
2444 }
2445
2446
2447 if (Method == 4) {
2448 for (int k1 = 0; k1 < NP; k1++) {
2449 YY[k1] = 0.;
2450 }
2451 c1->Range(0, 0, 25, 18);
2452 c1->SetFillColor(0);
2453 c1->Clear();
2454 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
2455 pad1->Draw();
2456 pad1->cd();
2457 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
2458
2459 pad1->DrawFrame(x1, y3, x2, y4);
2460
2461 pad1->GetFrame()->SetFillColor(0);
2462 pad1->Draw();
2463 gPad->SetGridy();
2464 gPad->SetGridx();
2465 cout << "****Draw 4 Dependencies: M r vs t for Depth2 and some eta 17 18 19 20 24 26 29 *** " << endl;
2466
2467
2468
2469 for (int k1 = 0; k1 < NP; k1++) {
2470 YY[k1] = ZZM[k1][2][3];
2471 }
2472 gr1 = new TGraph(NP, XX, YY);
2473 gr1->SetLineColor(1);
2474 gr1->SetMarkerColor(1);
2475 gr1->SetMarkerStyle(20);
2476 gr1->SetMarkerSize(1.2);
2477 gr1->SetTitle("ieta 26");
2478 gr1->SetFillStyle(0);
2479 gr1->Draw("P");
2480 for (int k1 = 0; k1 < NP; k1++) {
2481 YY[k1] = ZZM[k1][2][5];
2482 }
2483 gr2 = new TGraph(NP, XX, YY);
2484 gr2->SetLineColor(2);
2485 gr2->SetMarkerColor(2);
2486 gr2->SetMarkerStyle(20);
2487 gr2->SetMarkerSize(1.2);
2488 gr2->SetTitle("ieta 24");
2489 gr2->SetFillStyle(0);
2490 gr2->Draw("P");
2491 for (int k1 = 0; k1 < NP; k1++) {
2492 YY[k1] = ZZM[k1][2][9];
2493 }
2494 gr3 = new TGraph(NP, XX, YY);
2495 gr3->SetLineColor(3);
2496 gr3->SetMarkerColor(3);
2497 gr3->SetMarkerStyle(20);
2498 gr3->SetMarkerSize(1.2);
2499 gr3->SetTitle("ieta 20");
2500 gr3->SetFillStyle(0);
2501 gr3->Draw("P");
2502 for (int k1 = 0; k1 < NP; k1++) {
2503 YY[k1] = ZZM[k1][2][10];
2504 }
2505 gr4 = new TGraph(NP, XX, YY);
2506 gr4->SetLineColor(4);
2507 gr4->SetMarkerColor(4);
2508 gr4->SetMarkerStyle(20);
2509 gr4->SetMarkerSize(1.2);
2510 gr4->SetTitle("ieta 19");
2511 gr4->SetFillStyle(0);
2512 gr4->Draw("P");
2513 for (int k1 = 0; k1 < NP; k1++) {
2514 YY[k1] = ZZM[k1][2][11];
2515 }
2516 gr5 = new TGraph(NP, XX, YY);
2517 gr5->SetLineColor(5);
2518 gr5->SetMarkerColor(5);
2519 gr5->SetMarkerStyle(20);
2520 gr5->SetMarkerSize(1.2);
2521 gr5->SetTitle("ieta 18");
2522 gr5->SetFillStyle(0);
2523 gr5->Draw("P");
2524
2525
2526
2527 TLegend leg(xx1, yy1, xx2, yy2, " HEM Depth2");
2528 leg.AddEntry(gr1);
2529 leg.AddEntry(gr2);
2530 leg.AddEntry(gr3);
2531 leg.AddEntry(gr4);
2532 leg.AddEntry(gr5);
2533 leg.SetMargin(0.3);
2534 leg.SetBorderSize(1);
2535 leg.DrawClone("Same");
2536
2537 }
2538
2539
2540
2541
2542
2543 if (Method == 5) {
2544 for (int k1 = 0; k1 < NP; k1++) {
2545 YY[k1] = 0.;
2546 }
2547 c1->Range(0, 0, 25, 18);
2548 c1->SetFillColor(0);
2549 c1->Clear();
2550 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
2551 pad1->Draw();
2552 pad1->cd();
2553 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
2554
2555 pad1->DrawFrame(x1, y1, x2, y2);
2556
2557 pad1->GetFrame()->SetFillColor(0);
2558 pad1->Draw();
2559 gPad->SetGridy();
2560 gPad->SetGridx();
2561 cout << "****Draw 5 Dependencies: M r vs t for L1 and eta 21 22 23 25 27 28 **** " << endl;
2562
2563
2564
2565
2566
2567
2568
2569
2570
2571
2572
2573
2574
2575
2576
2577
2578
2579
2580 for (int k1 = 0; k1 < NP; k1++) {
2581 YY[k1] = 0.;
2582 if (LASERruns[k1] > "305874") {
2583 for (int iii = 2; iii < 3; iii++) {
2584 for (int jj = 0; jj < 13; jj++) {
2585 int LM1 = 0;
2586 if ((iii == 2 && jj < 13))
2587 LM1 = 1;
2588 if (LM1 == 1 && jj == 1) {
2589 YY[k1] = ZZM[k1][iii][jj];
2590 }
2591 }
2592 }
2593 } else {
2594 for (int iii = 1; iii < 2; iii++) {
2595 for (int jj = 0; jj < 13; jj++) {
2596 int LM1 = 0;
2597 if ((iii == 1 && jj < 13))
2598 LM1 = 1;
2599 if (LM1 == 1 && jj == 1) {
2600 YY[k1] = ZZM[k1][iii][jj];
2601 }
2602 }
2603 }
2604 }
2605 }
2606 gr1 = new TGraph(NP, XX, YY);
2607 gr1->SetLineColor(1);
2608 gr1->SetMarkerColor(1);
2609 gr1->SetMarkerStyle(20);
2610 gr1->SetMarkerSize(1.2);
2611 gr1->SetTitle("ieta 28");
2612 gr1->SetFillStyle(0);
2613 gr1->Draw("P");
2614
2615 for (int k1 = 0; k1 < NP; k1++) {
2616 YY[k1] = 0.;
2617 if (LASERruns[k1] > "305874") {
2618 for (int iii = 2; iii < 3; iii++) {
2619 for (int jj = 0; jj < 13; jj++) {
2620 int LM1 = 0;
2621 if ((iii == 2 && jj < 13))
2622 LM1 = 1;
2623 if (LM1 == 1 && jj == 2) {
2624 YY[k1] = ZZM[k1][iii][jj];
2625 }
2626 }
2627 }
2628 } else {
2629 for (int iii = 1; iii < 2; iii++) {
2630 for (int jj = 0; jj < 13; jj++) {
2631 int LM1 = 0;
2632 if ((iii == 1 && jj < 13))
2633 LM1 = 1;
2634 if (LM1 == 1 && jj == 2) {
2635 YY[k1] = ZZM[k1][iii][jj];
2636 }
2637 }
2638 }
2639 }
2640 }
2641 gr2 = new TGraph(NP, XX, YY);
2642 gr2->SetLineColor(2);
2643 gr2->SetMarkerColor(2);
2644 gr2->SetMarkerStyle(20);
2645 gr2->SetMarkerSize(1.2);
2646 gr2->SetTitle("ieta 27");
2647 gr2->SetFillStyle(0);
2648 gr2->Draw("P");
2649
2650 for (int k1 = 0; k1 < NP; k1++) {
2651 YY[k1] = 0.;
2652 if (LASERruns[k1] > "305874") {
2653 for (int iii = 2; iii < 3; iii++) {
2654 for (int jj = 0; jj < 13; jj++) {
2655 int LM1 = 0;
2656 if ((iii == 2 && jj < 13))
2657 LM1 = 1;
2658 if (LM1 == 1 && jj == 4) {
2659 YY[k1] = ZZM[k1][iii][jj];
2660 }
2661 }
2662 }
2663 } else {
2664 for (int iii = 1; iii < 2; iii++) {
2665 for (int jj = 0; jj < 13; jj++) {
2666 int LM1 = 0;
2667 if ((iii == 1 && jj < 13))
2668 LM1 = 1;
2669 if (LM1 == 1 && jj == 4) {
2670 YY[k1] = ZZM[k1][iii][jj];
2671 }
2672 }
2673 }
2674 }
2675 }
2676 gr3 = new TGraph(NP, XX, YY);
2677 gr3->SetLineColor(3);
2678 gr3->SetMarkerColor(3);
2679 gr3->SetMarkerStyle(20);
2680 gr3->SetMarkerSize(1.2);
2681 gr3->SetTitle("ieta 25");
2682 gr3->SetFillStyle(0);
2683 gr3->Draw("P");
2684
2685 for (int k1 = 0; k1 < NP; k1++) {
2686 YY[k1] = 0.;
2687 if (LASERruns[k1] > "305874") {
2688 for (int iii = 2; iii < 3; iii++) {
2689 for (int jj = 0; jj < 13; jj++) {
2690 int LM1 = 0;
2691 if ((iii == 2 && jj < 13))
2692 LM1 = 1;
2693 if (LM1 == 1 && jj == 6) {
2694 YY[k1] = ZZM[k1][iii][jj];
2695 }
2696 }
2697 }
2698 } else {
2699 for (int iii = 1; iii < 2; iii++) {
2700 for (int jj = 0; jj < 13; jj++) {
2701 int LM1 = 0;
2702 if ((iii == 1 && jj < 13))
2703 LM1 = 1;
2704 if (LM1 == 1 && jj == 6) {
2705 YY[k1] = ZZM[k1][iii][jj];
2706 }
2707 }
2708 }
2709 }
2710 }
2711 gr4 = new TGraph(NP, XX, YY);
2712 gr4->SetLineColor(4);
2713 gr4->SetMarkerColor(4);
2714 gr4->SetMarkerStyle(20);
2715 gr4->SetMarkerSize(1.2);
2716 gr4->SetTitle("ieta 23");
2717 gr4->SetFillStyle(0);
2718 gr4->Draw("P");
2719
2720 for (int k1 = 0; k1 < NP; k1++) {
2721 YY[k1] = 0.;
2722 if (LASERruns[k1] > "305874") {
2723 for (int iii = 2; iii < 3; iii++) {
2724 for (int jj = 0; jj < 13; jj++) {
2725 int LM1 = 0;
2726 if ((iii == 2 && jj < 13))
2727 LM1 = 1;
2728 if (LM1 == 1 && jj == 7) {
2729 YY[k1] = ZZM[k1][iii][jj];
2730 }
2731 }
2732 }
2733 } else {
2734 for (int iii = 1; iii < 2; iii++) {
2735 for (int jj = 0; jj < 13; jj++) {
2736 int LM1 = 0;
2737 if ((iii == 1 && jj < 13))
2738 LM1 = 1;
2739 if (LM1 == 1 && jj == 7) {
2740 YY[k1] = ZZM[k1][iii][jj];
2741 }
2742 }
2743 }
2744 }
2745 }
2746 gr5 = new TGraph(NP, XX, YY);
2747 gr5->SetLineColor(5);
2748 gr5->SetMarkerColor(5);
2749 gr5->SetMarkerStyle(20);
2750 gr5->SetMarkerSize(1.2);
2751 gr5->SetTitle("ieta 22");
2752 gr5->SetFillStyle(0);
2753 gr5->Draw("P");
2754
2755 for (int k1 = 0; k1 < NP; k1++) {
2756 YY[k1] = 0.;
2757 if (LASERruns[k1] > "305874") {
2758 for (int iii = 2; iii < 3; iii++) {
2759 for (int jj = 0; jj < 13; jj++) {
2760 int LM1 = 0;
2761 if ((iii == 2 && jj < 13))
2762 LM1 = 1;
2763 if (LM1 == 1 && jj == 8) {
2764 YY[k1] = ZZM[k1][iii][jj];
2765 }
2766 }
2767 }
2768 } else {
2769 for (int iii = 1; iii < 2; iii++) {
2770 for (int jj = 0; jj < 13; jj++) {
2771 int LM1 = 0;
2772 if ((iii == 1 && jj < 13))
2773 LM1 = 1;
2774 if (LM1 == 1 && jj == 8) {
2775 YY[k1] = ZZM[k1][iii][jj];
2776 }
2777 }
2778 }
2779 }
2780 }
2781 gr6 = new TGraph(NP, XX, YY);
2782 gr6->SetLineColor(6);
2783 gr6->SetMarkerColor(6);
2784 gr6->SetMarkerStyle(20);
2785 gr6->SetMarkerSize(1.2);
2786 gr6->SetTitle("ieta 21");
2787 gr6->SetFillStyle(0);
2788 gr6->Draw("P");
2789
2790 TLegend leg(xx1, yy1, xx2, yy2, " HEM L1");
2791 leg.AddEntry(gr1);
2792 leg.AddEntry(gr2);
2793 leg.AddEntry(gr3);
2794 leg.AddEntry(gr4);
2795 leg.AddEntry(gr5);
2796 leg.AddEntry(gr6);
2797 leg.SetMargin(0.3);
2798 leg.SetBorderSize(1);
2799 leg.DrawClone("Same");
2800
2801 }
2802
2803
2804
2805 if (Method == 6) {
2806 for (int k1 = 0; k1 < NP; k1++) {
2807 YY[k1] = 0.;
2808 }
2809 c1->Range(0, 0, 25, 18);
2810 c1->SetFillColor(0);
2811 c1->Clear();
2812 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
2813 pad1->Draw();
2814 pad1->cd();
2815 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
2816
2817 pad1->DrawFrame(x1, y1, x2, y2);
2818
2819 pad1->GetFrame()->SetFillColor(0);
2820 pad1->Draw();
2821 gPad->SetGridy();
2822 gPad->SetGridx();
2823 cout << "****Draw 6 Dependencies: M r vs t for L1 and eta 17 18 19 20 24 26 29 **** " << endl;
2824
2825
2826
2827 for (int k1 = 0; k1 < NP; k1++) {
2828 YY[k1] = 0.;
2829 if (LASERruns[k1] > "305874") {
2830 for (int iii = 2; iii < 3; iii++) {
2831 for (int jj = 0; jj < 13; jj++) {
2832 int LM1 = 0;
2833 if ((iii == 2 && jj < 13))
2834 LM1 = 1;
2835 if (LM1 == 1 && jj == 0) {
2836 YY[k1] = ZZM[k1][iii][jj];
2837 }
2838 }
2839 }
2840 } else {
2841 for (int iii = 1; iii < 2; iii++) {
2842 for (int jj = 0; jj < 13; jj++) {
2843 int LM1 = 0;
2844 if ((iii == 1 && jj < 13))
2845 LM1 = 1;
2846 if (LM1 == 1 && jj == 0) {
2847 YY[k1] = ZZM[k1][iii][jj];
2848 }
2849 }
2850 }
2851 }
2852 }
2853 gr1 = new TGraph(NP, XX, YY);
2854 gr1->SetLineColor(1);
2855 gr1->SetMarkerColor(1);
2856 gr1->SetMarkerStyle(20);
2857 gr1->SetMarkerSize(1.2);
2858 gr1->SetTitle("ieta 29");
2859 gr1->SetFillStyle(0);
2860 gr1->Draw("P");
2861
2862 for (int k1 = 0; k1 < NP; k1++) {
2863 YY[k1] = 0.;
2864 if (LASERruns[k1] > "305874") {
2865 for (int iii = 2; iii < 3; iii++) {
2866 for (int jj = 0; jj < 13; jj++) {
2867 int LM1 = 0;
2868 if ((iii == 2 && jj < 13))
2869 LM1 = 1;
2870 if (LM1 == 1 && jj == 3) {
2871 YY[k1] = ZZM[k1][iii][jj];
2872 }
2873 }
2874 }
2875 } else {
2876 for (int iii = 1; iii < 2; iii++) {
2877 for (int jj = 0; jj < 13; jj++) {
2878 int LM1 = 0;
2879 if ((iii == 1 && jj < 13))
2880 LM1 = 1;
2881 if (LM1 == 1 && jj == 3) {
2882 YY[k1] = ZZM[k1][iii][jj];
2883 }
2884 }
2885 }
2886 }
2887 }
2888 gr2 = new TGraph(NP, XX, YY);
2889 gr2->SetLineColor(2);
2890 gr2->SetMarkerColor(2);
2891 gr2->SetMarkerStyle(20);
2892 gr2->SetMarkerSize(1.2);
2893 gr2->SetTitle("ieta 26");
2894 gr2->SetFillStyle(0);
2895 gr2->Draw("P");
2896
2897 for (int k1 = 0; k1 < NP; k1++) {
2898 YY[k1] = 0.;
2899 if (LASERruns[k1] > "305874") {
2900 for (int iii = 2; iii < 3; iii++) {
2901 for (int jj = 0; jj < 13; jj++) {
2902 int LM1 = 0;
2903 if ((iii == 2 && jj < 13))
2904 LM1 = 1;
2905 if (LM1 == 1 && jj == 5) {
2906 YY[k1] = ZZM[k1][iii][jj];
2907 }
2908 }
2909 }
2910 } else {
2911 for (int iii = 1; iii < 2; iii++) {
2912 for (int jj = 0; jj < 13; jj++) {
2913 int LM1 = 0;
2914 if ((iii == 1 && jj < 13))
2915 LM1 = 1;
2916 if (LM1 == 1 && jj == 5) {
2917 YY[k1] = ZZM[k1][iii][jj];
2918 }
2919 }
2920 }
2921 }
2922 }
2923 gr3 = new TGraph(NP, XX, YY);
2924 gr3->SetLineColor(3);
2925 gr3->SetMarkerColor(3);
2926 gr3->SetMarkerStyle(20);
2927 gr3->SetMarkerSize(1.2);
2928 gr3->SetTitle("ieta 24");
2929 gr3->SetFillStyle(0);
2930 gr3->Draw("P");
2931
2932
2933
2934
2935 for (int k1 = 0; k1 < NP; k1++) {
2936 YY[k1] = 0.;
2937 if (LASERruns[k1] > "305874") {
2938 for (int iii = 2; iii < 3; iii++) {
2939 for (int jj = 0; jj < 13; jj++) {
2940 int LM1 = 0;
2941 if ((iii == 2 && jj < 13))
2942 LM1 = 1;
2943 if (LM1 == 1 && jj == 10) {
2944 YY[k1] = ZZM[k1][iii][jj];
2945 }
2946 }
2947 }
2948 } else {
2949 for (int iii = 1; iii < 2; iii++) {
2950 for (int jj = 0; jj < 13; jj++) {
2951 int LM1 = 0;
2952 if ((iii == 1 && jj < 13))
2953 LM1 = 1;
2954 if (LM1 == 1 && jj == 10) {
2955 YY[k1] = ZZM[k1][iii][jj];
2956 }
2957 }
2958 }
2959 }
2960 }
2961 gr4 = new TGraph(NP, XX, YY);
2962 gr4->SetLineColor(5);
2963 gr4->SetMarkerColor(5);
2964 gr4->SetMarkerStyle(20);
2965 gr4->SetMarkerSize(1.2);
2966 gr4->SetTitle("ieta 19");
2967 gr4->SetFillStyle(0);
2968 gr4->Draw("P");
2969
2970
2971
2972
2973 for (int k1 = 0; k1 < NP; k1++) {
2974 YY[k1] = 0.;
2975 if (LASERruns[k1] > "305874") {
2976 for (int iii = 2; iii < 3; iii++) {
2977 for (int jj = 0; jj < 13; jj++) {
2978 int LM1 = 0;
2979 if ((iii == 2 && jj < 13))
2980 LM1 = 1;
2981 if (LM1 == 1 && jj == 12) {
2982 YY[k1] = ZZM[k1][iii][jj];
2983 }
2984 }
2985 }
2986 } else {
2987 for (int iii = 1; iii < 2; iii++) {
2988 for (int jj = 0; jj < 13; jj++) {
2989 int LM1 = 0;
2990 if ((iii == 1 && jj < 13))
2991 LM1 = 1;
2992 if (LM1 == 1 && jj == 12) {
2993 YY[k1] = ZZM[k1][iii][jj];
2994 }
2995 }
2996 }
2997 }
2998 }
2999 gr5 = new TGraph(NP, XX, YY);
3000 gr5->SetLineColor(7);
3001 gr5->SetMarkerColor(7);
3002 gr5->SetMarkerStyle(20);
3003 gr5->SetMarkerSize(1.2);
3004 gr5->SetTitle("ieta 17");
3005 gr5->SetFillStyle(0);
3006 gr5->Draw("P");
3007
3008 TLegend leg(xx1, yy1, xx2, yy2, " HEM L1");
3009 leg.AddEntry(gr1);
3010 leg.AddEntry(gr2);
3011 leg.AddEntry(gr3);
3012 leg.AddEntry(gr4);
3013 leg.AddEntry(gr5);
3014 leg.SetMargin(0.3);
3015 leg.SetBorderSize(1);
3016 leg.DrawClone("Same");
3017 }
3018
3019
3020
3021 if (Method == 7) {
3022 for (int k1 = 0; k1 < NP; k1++) {
3023 YY[k1] = 0.;
3024 }
3025
3026
3027 cout << "**** Negative direction: take good runs eta 21 **** " << endl;
3028 c1->Range(0, 0, 25, 18);
3029 c1->SetFillColor(0);
3030 c1->Clear();
3031 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
3032 pad1->Draw();
3033 pad1->cd();
3034 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
3035
3036 pad1->DrawFrame(x1, y5, x2, y6);
3037
3038 pad1->GetFrame()->SetFillColor(0);
3039 pad1->Draw();
3040 gPad->SetGridy();
3041 gPad->SetGridx();
3042 cout << "****Draw 7 Dependencies: M r vs t for L7 and eta 21 22 23 25 27 28 **** " << endl;
3043
3044
3045
3046
3047
3048
3049
3050
3051 for (int k1 = 0; k1 < NP; k1++) {
3052 YY[k1] = 0.;
3053 if (LASERruns[k1] > "305874") {
3054 for (int iii = 3; iii < 6; iii++) {
3055 for (int jj = 1; jj < 13; jj++) {
3056 int LM7 = 0;
3057 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3058 LM7 = 1;
3059 if (LM7 == 1 && jj == 1) {
3060 YY[k1] = ZZM[k1][iii][jj];
3061 }
3062 }
3063 }
3064 } else {
3065 for (int iii = 2; iii < 4; iii++) {
3066 for (int jj = 1; jj < 13; jj++) {
3067 int LM7 = 0;
3068 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3069 LM7 = 1;
3070 if (LM7 == 1 && jj == 1) {
3071 YY[k1] = ZZM[k1][iii][jj];
3072 }
3073 }
3074 }
3075 }
3076 }
3077 gr1 = new TGraph(NP, XX, YY);
3078 gr1->SetLineColor(1);
3079 gr1->SetMarkerColor(1);
3080 gr1->SetMarkerStyle(20);
3081 gr1->SetTitle("ieta 28");
3082 gr1->SetFillStyle(0);
3083 gr1->SetMarkerSize(1.2);
3084 gr1->Draw("P");
3085
3086 for (int k1 = 0; k1 < NP; k1++) {
3087 YY[k1] = 0.;
3088 if (LASERruns[k1] > "305874") {
3089 for (int iii = 3; iii < 6; iii++) {
3090 for (int jj = 1; jj < 13; jj++) {
3091 int LM7 = 0;
3092 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3093 LM7 = 1;
3094 if (LM7 == 1 && jj == 2) {
3095 YY[k1] = ZZM[k1][iii][jj];
3096 }
3097 }
3098 }
3099 } else {
3100 for (int iii = 2; iii < 4; iii++) {
3101 for (int jj = 1; jj < 13; jj++) {
3102 int LM7 = 0;
3103 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3104 LM7 = 1;
3105 if (LM7 == 1 && jj == 2) {
3106 YY[k1] = ZZM[k1][iii][jj];
3107 }
3108 }
3109 }
3110 }
3111 }
3112 gr2 = new TGraph(NP, XX, YY);
3113 gr2->SetLineColor(2);
3114 gr2->SetMarkerColor(2);
3115 gr2->SetMarkerStyle(20);
3116 gr2->SetTitle("ieta 27");
3117 gr2->SetFillStyle(0);
3118 gr2->SetMarkerSize(1.2);
3119 gr2->Draw("P");
3120
3121 for (int k1 = 0; k1 < NP; k1++) {
3122 YY[k1] = 0.;
3123 if (LASERruns[k1] > "305874") {
3124 for (int iii = 3; iii < 6; iii++) {
3125 for (int jj = 1; jj < 13; jj++) {
3126 int LM7 = 0;
3127 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3128 LM7 = 1;
3129 if (LM7 == 1 && jj == 4) {
3130 YY[k1] = ZZM[k1][iii][jj];
3131 }
3132 }
3133 }
3134 } else {
3135 for (int iii = 2; iii < 4; iii++) {
3136 for (int jj = 1; jj < 13; jj++) {
3137 int LM7 = 0;
3138 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3139 LM7 = 1;
3140 if (LM7 == 1 && jj == 4) {
3141 YY[k1] = ZZM[k1][iii][jj];
3142 }
3143 }
3144 }
3145 }
3146 }
3147 gr3 = new TGraph(NP, XX, YY);
3148 gr3->SetLineColor(3);
3149 gr3->SetMarkerColor(3);
3150 gr3->SetMarkerStyle(20);
3151 gr3->SetTitle("ieta 25");
3152 gr3->SetFillStyle(0);
3153 gr3->SetMarkerSize(1.2);
3154 gr3->Draw("P");
3155
3156 for (int k1 = 0; k1 < NP; k1++) {
3157 YY[k1] = 0.;
3158 if (LASERruns[k1] > "305874") {
3159 for (int iii = 3; iii < 6; iii++) {
3160 for (int jj = 1; jj < 13; jj++) {
3161 int LM7 = 0;
3162 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3163 LM7 = 1;
3164 if (LM7 == 1 && jj == 6) {
3165 YY[k1] = ZZM[k1][iii][jj];
3166 }
3167 }
3168 }
3169 } else {
3170 for (int iii = 2; iii < 4; iii++) {
3171 for (int jj = 1; jj < 13; jj++) {
3172 int LM7 = 0;
3173 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3174 LM7 = 1;
3175 if (LM7 == 1 && jj == 6) {
3176 YY[k1] = ZZM[k1][iii][jj];
3177 }
3178 }
3179 }
3180 }
3181 }
3182 gr4 = new TGraph(NP, XX, YY);
3183 gr4->SetLineColor(4);
3184 gr4->SetMarkerColor(4);
3185 gr4->SetMarkerStyle(20);
3186 gr4->SetTitle("ieta 23");
3187 gr4->SetFillStyle(0);
3188 gr4->SetMarkerSize(1.2);
3189 gr4->Draw("P");
3190
3191 for (int k1 = 0; k1 < NP; k1++) {
3192 YY[k1] = 0.;
3193 if (LASERruns[k1] > "305874") {
3194 for (int iii = 3; iii < 6; iii++) {
3195 for (int jj = 1; jj < 13; jj++) {
3196 int LM7 = 0;
3197 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3198 LM7 = 1;
3199 if (LM7 == 1 && jj == 7) {
3200 YY[k1] = ZZM[k1][iii][jj];
3201 }
3202 }
3203 }
3204 } else {
3205 for (int iii = 2; iii < 4; iii++) {
3206 for (int jj = 1; jj < 13; jj++) {
3207 int LM7 = 0;
3208 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3209 LM7 = 1;
3210 if (LM7 == 1 && jj == 7) {
3211 YY[k1] = ZZM[k1][iii][jj];
3212 }
3213 }
3214 }
3215 }
3216 }
3217 gr5 = new TGraph(NP, XX, YY);
3218 gr5->SetLineColor(5);
3219 gr5->SetMarkerColor(5);
3220 gr5->SetMarkerStyle(20);
3221 gr5->SetTitle("ieta 22");
3222 gr5->SetFillStyle(0);
3223 gr5->SetMarkerSize(1.2);
3224 gr5->Draw("P");
3225
3226 for (int k1 = 0; k1 < NP; k1++) {
3227 YY[k1] = 0.;
3228 if (LASERruns[k1] > "305874") {
3229 for (int iii = 3; iii < 6; iii++) {
3230 for (int jj = 1; jj < 13; jj++) {
3231 int LM7 = 0;
3232 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3233 LM7 = 1;
3234 if (LM7 == 1 && jj == 8) {
3235 YY[k1] = ZZM[k1][iii][jj];
3236 }
3237 }
3238 }
3239 } else {
3240 for (int iii = 2; iii < 4; iii++) {
3241 for (int jj = 1; jj < 13; jj++) {
3242 int LM7 = 0;
3243 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3244 LM7 = 1;
3245 if (LM7 == 1 && jj == 8) {
3246 YY[k1] = ZZM[k1][iii][jj];
3247 }
3248 }
3249 }
3250 }
3251
3252
3253
3254 }
3255 gr6 = new TGraph(NP, XX, YY);
3256 gr6->SetLineColor(6);
3257 gr6->SetMarkerColor(6);
3258 gr6->SetMarkerStyle(20);
3259 gr6->SetTitle("ieta 21");
3260 gr6->SetFillStyle(0);
3261 gr6->SetMarkerSize(1.2);
3262 gr6->Draw("P");
3263
3264
3265 TLegend leg(xx1, yy1, xx2, yy2, " HEM L7");
3266 leg.AddEntry(gr1);
3267 leg.AddEntry(gr2);
3268 leg.AddEntry(gr3);
3269 leg.AddEntry(gr4);
3270 leg.AddEntry(gr5);
3271 leg.AddEntry(gr6);
3272 leg.SetMargin(0.3);
3273 leg.SetBorderSize(1);
3274 leg.DrawClone("Same");
3275 }
3276
3277
3278 if (Method == 8) {
3279 for (int k1 = 0; k1 < NP; k1++) {
3280 YY[k1] = 0.;
3281 }
3282 c1->Range(0, 0, 25, 18);
3283 c1->SetFillColor(0);
3284 c1->Clear();
3285 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
3286 pad1->Draw();
3287 pad1->cd();
3288 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
3289
3290 pad1->DrawFrame(x1, y5, x2, y6);
3291
3292 pad1->GetFrame()->SetFillColor(0);
3293 pad1->Draw();
3294 gPad->SetGridy();
3295 gPad->SetGridx();
3296 cout << "****Draw 8 Dependencies: M r vs t for L7 and eta 17 18 19 20 24 26 29 **** " << endl;
3297
3298
3299
3300
3301
3302
3303
3304
3305 for (int k1 = 0; k1 < NP; k1++) {
3306 YY[k1] = 0.;
3307 if (LASERruns[k1] > "305874") {
3308 for (int iii = 3; iii < 6; iii++) {
3309 for (int jj = 1; jj < 13; jj++) {
3310 int LM7 = 0;
3311 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3312 LM7 = 1;
3313 if (LM7 == 1 && jj == 3) {
3314 YY[k1] = ZZM[k1][iii][jj];
3315 }
3316 }
3317 }
3318 } else {
3319 for (int iii = 2; iii < 4; iii++) {
3320 for (int jj = 1; jj < 13; jj++) {
3321 int LM7 = 0;
3322 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3323 LM7 = 1;
3324 if (LM7 == 1 && jj == 3) {
3325 YY[k1] = ZZM[k1][iii][jj];
3326 }
3327 }
3328 }
3329 }
3330 }
3331 gr1 = new TGraph(NP, XX, YY);
3332 gr1->SetLineColor(1);
3333 gr1->SetMarkerColor(1);
3334 gr1->SetMarkerStyle(20);
3335 gr1->SetTitle("ieta 26");
3336 gr1->SetFillStyle(0);
3337 gr1->SetMarkerSize(1.2);
3338 gr1->Draw("P");
3339
3340 for (int k1 = 0; k1 < NP; k1++) {
3341 YY[k1] = 0.;
3342 if (LASERruns[k1] > "305874") {
3343 for (int iii = 3; iii < 6; iii++) {
3344 for (int jj = 1; jj < 13; jj++) {
3345 int LM7 = 0;
3346 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3347 LM7 = 1;
3348 if (LM7 == 1 && jj == 5) {
3349 YY[k1] = ZZM[k1][iii][jj];
3350 }
3351 }
3352 }
3353 } else {
3354 for (int iii = 2; iii < 4; iii++) {
3355 for (int jj = 1; jj < 13; jj++) {
3356 int LM7 = 0;
3357 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3358 LM7 = 1;
3359 if (LM7 == 1 && jj == 5) {
3360 YY[k1] = ZZM[k1][iii][jj];
3361 }
3362 }
3363 }
3364 }
3365 }
3366 gr2 = new TGraph(NP, XX, YY);
3367 gr2->SetLineColor(2);
3368 gr2->SetMarkerColor(2);
3369 gr2->SetMarkerStyle(20);
3370 gr2->SetTitle("ieta 24");
3371 gr2->SetFillStyle(0);
3372 gr2->SetMarkerSize(1.2);
3373 gr2->Draw("P");
3374
3375 for (int k1 = 0; k1 < NP; k1++) {
3376 YY[k1] = 0.;
3377 if (LASERruns[k1] > "305874") {
3378 for (int iii = 3; iii < 6; iii++) {
3379 for (int jj = 1; jj < 13; jj++) {
3380 int LM7 = 0;
3381 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3382 LM7 = 1;
3383 if (LM7 == 1 && jj == 9) {
3384 YY[k1] = ZZM[k1][iii][jj];
3385 }
3386 }
3387 }
3388 } else {
3389 for (int iii = 2; iii < 4; iii++) {
3390 for (int jj = 1; jj < 13; jj++) {
3391 int LM7 = 0;
3392 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3393 LM7 = 1;
3394 if (LM7 == 1 && jj == 9) {
3395 YY[k1] = ZZM[k1][iii][jj];
3396 }
3397 }
3398 }
3399 }
3400 }
3401 gr3 = new TGraph(NP, XX, YY);
3402 gr3->SetLineColor(3);
3403 gr3->SetMarkerColor(3);
3404 gr3->SetMarkerStyle(20);
3405 gr3->SetTitle("ieta 20");
3406 gr3->SetFillStyle(0);
3407 gr3->SetMarkerSize(1.2);
3408 gr3->Draw("P");
3409
3410 for (int k1 = 0; k1 < NP; k1++) {
3411 YY[k1] = 0.;
3412 if (LASERruns[k1] > "305874") {
3413 for (int iii = 3; iii < 6; iii++) {
3414 for (int jj = 1; jj < 13; jj++) {
3415 int LM7 = 0;
3416 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3417 LM7 = 1;
3418 if (LM7 == 1 && jj == 10) {
3419 YY[k1] = ZZM[k1][iii][jj];
3420 }
3421 }
3422 }
3423 } else {
3424 for (int iii = 2; iii < 4; iii++) {
3425 for (int jj = 1; jj < 13; jj++) {
3426 int LM7 = 0;
3427 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3428 LM7 = 1;
3429 if (LM7 == 1 && jj == 10) {
3430 YY[k1] = ZZM[k1][iii][jj];
3431 }
3432 }
3433 }
3434 }
3435 }
3436 gr4 = new TGraph(NP, XX, YY);
3437 gr4->SetLineColor(4);
3438 gr4->SetMarkerColor(4);
3439 gr4->SetMarkerStyle(20);
3440 gr4->SetTitle("ieta 19");
3441 gr4->SetFillStyle(0);
3442 gr4->SetMarkerSize(1.2);
3443 gr4->Draw("P");
3444
3445 for (int k1 = 0; k1 < NP; k1++) {
3446 YY[k1] = 0.;
3447 if (LASERruns[k1] > "305874") {
3448 for (int iii = 3; iii < 6; iii++) {
3449 for (int jj = 1; jj < 13; jj++) {
3450 int LM7 = 0;
3451 if ((iii == 3 && jj == 12) || (iii == 4 && jj > 3 && jj < 12) || (iii == 5 && jj > 0 && jj < 4))
3452 LM7 = 1;
3453 if (LM7 == 1 && jj == 11) {
3454 YY[k1] = ZZM[k1][iii][jj];
3455 }
3456 }
3457 }
3458 } else {
3459 for (int iii = 2; iii < 4; iii++) {
3460 for (int jj = 1; jj < 13; jj++) {
3461 int LM7 = 0;
3462 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
3463 LM7 = 1;
3464 if (LM7 == 1 && jj == 11) {
3465 YY[k1] = ZZM[k1][iii][jj];
3466 }
3467 }
3468 }
3469 }
3470 }
3471 gr5 = new TGraph(NP, XX, YY);
3472 gr5->SetLineColor(5);
3473 gr5->SetMarkerColor(5);
3474 gr5->SetMarkerStyle(20);
3475 gr5->SetTitle("ieta 18");
3476 gr5->SetFillStyle(0);
3477 gr5->SetMarkerSize(1.2);
3478 gr5->Draw("P");
3479
3480
3481
3482
3483 TLegend leg(xx1, yy1, xx2, yy2, " HEM L7");
3484 leg.AddEntry(gr1);
3485 leg.AddEntry(gr2);
3486 leg.AddEntry(gr3);
3487 leg.AddEntry(gr4);
3488 leg.AddEntry(gr5);
3489 leg.SetMargin(0.3);
3490 leg.SetBorderSize(1);
3491 leg.DrawClone("Same");
3492 }
3493
3494
3495
3496
3497
3498
3499
3500
3501
3502 if (Method == 11) {
3503 for (int k1 = 0; k1 < NP; k1++) {
3504 YY[k1] = 0.;
3505 }
3506 c1->Range(0, 0, 25, 18);
3507 c1->SetFillColor(0);
3508 c1->Clear();
3509 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
3510 pad1->Draw();
3511 pad1->cd();
3512 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
3513
3514 pad1->DrawFrame(x1, y1, x2, y2);
3515
3516 pad1->GetFrame()->SetFillColor(0);
3517 pad1->Draw();
3518 gPad->SetGridy();
3519 gPad->SetGridx();
3520 cout << "****Draw 11 Dependencies: P r vs t for Depth1 and some eta 21 22 23 25 27 28 *** " << endl;
3521
3522
3523
3524 for (int k1 = 0; k1 < NP; k1++) {
3525 YY[k1] = ZZP[k1][1][12];
3526 }
3527 gr1 = new TGraph(NP, XX, YY);
3528 gr1->SetLineColor(1);
3529 gr1->SetMarkerColor(1);
3530 gr1->SetMarkerStyle(20);
3531 gr1->SetMarkerSize(1.2);
3532 gr1->SetTitle("ieta 28");
3533 gr1->SetFillStyle(0);
3534 gr1->Draw("P");
3535 for (int k1 = 0; k1 < NP; k1++) {
3536 YY[k1] = ZZP[k1][1][11];
3537 }
3538 gr2 = new TGraph(NP, XX, YY);
3539 gr2->SetLineColor(2);
3540 gr2->SetMarkerColor(2);
3541 gr2->SetMarkerStyle(20);
3542 gr2->SetMarkerSize(1.2);
3543 gr2->SetTitle("ieta 27");
3544 gr2->SetFillStyle(0);
3545 gr2->Draw("P");
3546 for (int k1 = 0; k1 < NP; k1++) {
3547 YY[k1] = ZZP[k1][1][9];
3548 }
3549 gr3 = new TGraph(NP, XX, YY);
3550 gr3->SetLineColor(3);
3551 gr3->SetMarkerColor(3);
3552 gr3->SetMarkerStyle(20);
3553 gr3->SetMarkerSize(1.2);
3554 gr3->SetTitle("ieta 25");
3555 gr3->SetFillStyle(0);
3556 gr3->Draw("P");
3557 for (int k1 = 0; k1 < NP; k1++) {
3558 YY[k1] = ZZP[k1][1][7];
3559 }
3560 gr4 = new TGraph(NP, XX, YY);
3561 gr4->SetLineColor(4);
3562 gr4->SetMarkerColor(4);
3563 gr4->SetMarkerStyle(20);
3564 gr4->SetMarkerSize(1.2);
3565 gr4->SetTitle("ieta 23");
3566 gr4->SetFillStyle(0);
3567 gr4->Draw("P");
3568 for (int k1 = 0; k1 < NP; k1++) {
3569 YY[k1] = ZZP[k1][1][6];
3570 }
3571 gr5 = new TGraph(NP, XX, YY);
3572 gr5->SetLineColor(5);
3573 gr5->SetMarkerColor(5);
3574 gr5->SetMarkerStyle(20);
3575 gr5->SetMarkerSize(1.2);
3576 gr5->SetTitle("ieta 22");
3577 gr5->SetFillStyle(0);
3578 gr5->Draw("P");
3579 for (int k1 = 0; k1 < NP; k1++) {
3580 YY[k1] = ZZP[k1][1][5];
3581 }
3582 gr6 = new TGraph(NP, XX, YY);
3583 gr6->SetLineColor(6);
3584 gr6->SetMarkerColor(6);
3585 gr6->SetMarkerStyle(20);
3586 gr6->SetMarkerSize(1.2);
3587 gr6->SetTitle("ieta 21");
3588 gr6->SetFillStyle(0);
3589 gr6->Draw("P");
3590
3591 TLegend leg(xx1, yy1, xx2, yy2, " HEP Depth1");
3592 leg.AddEntry(gr1);
3593 leg.AddEntry(gr2);
3594 leg.AddEntry(gr3);
3595 leg.AddEntry(gr4);
3596 leg.AddEntry(gr5);
3597 leg.AddEntry(gr6);
3598 leg.SetMargin(0.3);
3599 leg.SetBorderSize(1);
3600 leg.DrawClone("Same");
3601
3602 }
3603
3604 if (Method == 12) {
3605 for (int k1 = 0; k1 < NP; k1++) {
3606 YY[k1] = 0.;
3607 }
3608 c1->Range(0, 0, 25, 18);
3609 c1->SetFillColor(0);
3610 c1->Clear();
3611 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
3612 pad1->Draw();
3613 pad1->cd();
3614 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
3615
3616 pad1->DrawFrame(x1, y1, x2, y2);
3617
3618 pad1->GetFrame()->SetFillColor(0);
3619 pad1->Draw();
3620 gPad->SetGridy();
3621 gPad->SetGridx();
3622 cout << "****Draw 12 Dependencies: P r vs t for Depth1 and some eta 17 18 19 20 24 26 29 *** " << endl;
3623
3624
3625
3626 for (int k1 = 0; k1 < NP; k1++) {
3627 YY[k1] = ZZP[k1][1][13];
3628
3629
3630
3631
3632 }
3633 gr1 = new TGraph(NP, XX, YY);
3634 gr1->SetLineColor(1);
3635 gr1->SetMarkerColor(1);
3636 gr1->SetMarkerStyle(20);
3637 gr1->SetMarkerSize(1.2);
3638 gr1->SetTitle("ieta 29");
3639 gr1->SetFillStyle(0);
3640 gr1->Draw("P");
3641 for (int k1 = 0; k1 < NP; k1++) {
3642 YY[k1] = ZZP[k1][1][10];
3643 }
3644 gr2 = new TGraph(NP, XX, YY);
3645 gr2->SetLineColor(2);
3646 gr2->SetMarkerColor(2);
3647 gr2->SetMarkerStyle(20);
3648 gr2->SetMarkerSize(1.2);
3649 gr2->SetTitle("ieta 26");
3650 gr2->SetFillStyle(0);
3651 gr2->Draw("P");
3652 for (int k1 = 0; k1 < NP; k1++) {
3653 YY[k1] = ZZP[k1][1][8];
3654 }
3655 gr3 = new TGraph(NP, XX, YY);
3656 gr3->SetLineColor(3);
3657 gr3->SetMarkerColor(3);
3658 gr3->SetMarkerStyle(20);
3659 gr3->SetMarkerSize(1.2);
3660 gr3->SetTitle("ieta 24");
3661 gr3->SetFillStyle(0);
3662 gr3->Draw("P");
3663
3664
3665 for (int k1 = 0; k1 < NP; k1++) {
3666 YY[k1] = ZZP[k1][1][3];
3667 }
3668 gr4 = new TGraph(NP, XX, YY);
3669 gr4->SetLineColor(5);
3670 gr4->SetMarkerColor(5);
3671 gr4->SetMarkerStyle(20);
3672 gr4->SetMarkerSize(1.2);
3673 gr4->SetTitle("ieta 19");
3674 gr4->SetFillStyle(0);
3675 gr4->Draw("P");
3676
3677
3678 for (int k1 = 0; k1 < NP; k1++) {
3679 YY[k1] = ZZP[k1][1][1];
3680 }
3681 gr5 = new TGraph(NP, XX, YY);
3682 gr5->SetLineColor(7);
3683 gr5->SetMarkerColor(7);
3684 gr5->SetMarkerStyle(20);
3685 gr5->SetMarkerSize(1.2);
3686 gr5->SetTitle("ieta 17");
3687 gr5->SetFillStyle(0);
3688 gr5->Draw("P");
3689
3690 TLegend leg(xx1, yy1, xx2, yy2, " HEP Depth1");
3691 leg.AddEntry(gr1);
3692 leg.AddEntry(gr2);
3693 leg.AddEntry(gr3);
3694 leg.AddEntry(gr4);
3695 leg.AddEntry(gr5);
3696 leg.SetMargin(0.3);
3697 leg.SetBorderSize(1);
3698 leg.DrawClone("Same");
3699
3700 }
3701
3702
3703 if (Method == 13) {
3704 for (int k1 = 0; k1 < NP; k1++) {
3705 YY[k1] = 0.;
3706 }
3707 c1->Range(0, 0, 25, 18);
3708 c1->SetFillColor(0);
3709 c1->Clear();
3710 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
3711 pad1->Draw();
3712 pad1->cd();
3713 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
3714
3715 pad1->DrawFrame(x1, y3, x2, y4);
3716
3717 pad1->GetFrame()->SetFillColor(0);
3718 pad1->Draw();
3719 gPad->SetGridy();
3720 gPad->SetGridx();
3721 cout << "****Draw 13 Dependencies: P r vs t for Depth2 and some eta 21 22 23 25 27 28 *** " << endl;
3722
3723
3724
3725 for (int k1 = 0; k1 < NP; k1++) {
3726 YY[k1] = ZZP[k1][2][9];
3727 }
3728 gr1 = new TGraph(NP, XX, YY);
3729 gr1->SetLineColor(1);
3730 gr1->SetMarkerColor(1);
3731 gr1->SetMarkerStyle(20);
3732 gr1->SetMarkerSize(1.2);
3733 gr1->SetTitle("ieta 25");
3734 gr1->SetFillStyle(0);
3735 gr1->Draw("P");
3736 for (int k1 = 0; k1 < NP; k1++) {
3737 YY[k1] = ZZP[k1][2][7];
3738 }
3739 gr2 = new TGraph(NP, XX, YY);
3740 gr2->SetLineColor(2);
3741 gr2->SetMarkerColor(2);
3742 gr2->SetMarkerStyle(20);
3743 gr2->SetMarkerSize(1.2);
3744 gr2->SetTitle("ieta 23");
3745 gr2->SetFillStyle(0);
3746 gr2->Draw("P");
3747 for (int k1 = 0; k1 < NP; k1++) {
3748 YY[k1] = ZZP[k1][2][6];
3749 }
3750 gr3 = new TGraph(NP, XX, YY);
3751 gr3->SetLineColor(3);
3752 gr3->SetMarkerColor(3);
3753 gr3->SetMarkerStyle(20);
3754 gr3->SetMarkerSize(1.2);
3755 gr3->SetTitle("ieta 22");
3756 gr3->SetFillStyle(0);
3757 gr3->Draw("P");
3758 for (int k1 = 0; k1 < NP; k1++) {
3759 YY[k1] = ZZP[k1][2][5];
3760 }
3761 gr4 = new TGraph(NP, XX, YY);
3762 gr4->SetLineColor(4);
3763 gr4->SetMarkerColor(4);
3764 gr4->SetMarkerStyle(20);
3765 gr4->SetMarkerSize(1.2);
3766 gr4->SetTitle("ieta 21");
3767 gr4->SetFillStyle(0);
3768 gr4->Draw("P");
3769
3770 TLegend leg(xx1, yy1, xx2, yy2, " HEP Depth2");
3771 leg.AddEntry(gr1);
3772 leg.AddEntry(gr2);
3773 leg.AddEntry(gr3);
3774 leg.AddEntry(gr4);
3775 leg.SetMargin(0.3);
3776 leg.SetBorderSize(1);
3777 leg.DrawClone("Same");
3778
3779 }
3780
3781
3782 if (Method == 14) {
3783 for (int k1 = 0; k1 < NP; k1++) {
3784 YY[k1] = 0.;
3785 }
3786 c1->Range(0, 0, 25, 18);
3787 c1->SetFillColor(0);
3788 c1->Clear();
3789 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
3790 pad1->Draw();
3791 pad1->cd();
3792 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
3793
3794 pad1->DrawFrame(x1, y3, x2, y4);
3795
3796 pad1->GetFrame()->SetFillColor(0);
3797 pad1->Draw();
3798 gPad->SetGridy();
3799 gPad->SetGridx();
3800 cout << "****Draw 14 Dependencies: P r vs t for Depth2 and some eta 17 18 19 20 24 26 29 *** " << endl;
3801
3802
3803
3804 for (int k1 = 0; k1 < NP; k1++) {
3805 YY[k1] = ZZP[k1][2][10];
3806 }
3807 gr1 = new TGraph(NP, XX, YY);
3808 gr1->SetLineColor(1);
3809 gr1->SetMarkerColor(1);
3810 gr1->SetMarkerStyle(20);
3811 gr1->SetMarkerSize(1.2);
3812 gr1->SetTitle("ieta 26");
3813 gr1->SetFillStyle(0);
3814 gr1->Draw("P");
3815 for (int k1 = 0; k1 < NP; k1++) {
3816 YY[k1] = ZZP[k1][2][8];
3817 }
3818 gr2 = new TGraph(NP, XX, YY);
3819 gr2->SetLineColor(2);
3820 gr2->SetMarkerColor(2);
3821 gr2->SetMarkerStyle(20);
3822 gr2->SetMarkerSize(1.2);
3823 gr2->SetTitle("ieta 24");
3824 gr2->SetFillStyle(0);
3825 gr2->Draw("P");
3826 for (int k1 = 0; k1 < NP; k1++) {
3827 YY[k1] = ZZP[k1][2][4];
3828 }
3829 gr3 = new TGraph(NP, XX, YY);
3830 gr3->SetLineColor(3);
3831 gr3->SetMarkerColor(3);
3832 gr3->SetMarkerStyle(20);
3833 gr3->SetMarkerSize(1.2);
3834 gr3->SetTitle("ieta 20");
3835 gr3->SetFillStyle(0);
3836 gr3->Draw("P");
3837 for (int k1 = 0; k1 < NP; k1++) {
3838 YY[k1] = ZZP[k1][2][3];
3839 }
3840 gr4 = new TGraph(NP, XX, YY);
3841 gr4->SetLineColor(4);
3842 gr4->SetMarkerColor(4);
3843 gr4->SetMarkerStyle(20);
3844 gr4->SetMarkerSize(1.2);
3845 gr4->SetTitle("ieta 19");
3846 gr4->SetFillStyle(0);
3847 gr4->Draw("P");
3848 for (int k1 = 0; k1 < NP; k1++) {
3849 YY[k1] = ZZP[k1][2][2];
3850 }
3851 gr5 = new TGraph(NP, XX, YY);
3852 gr5->SetLineColor(5);
3853 gr5->SetMarkerColor(5);
3854 gr5->SetMarkerStyle(20);
3855 gr5->SetMarkerSize(1.2);
3856 gr5->SetTitle("ieta 18");
3857 gr5->SetFillStyle(0);
3858 gr5->Draw("P");
3859
3860
3861
3862 TLegend leg(xx1, yy1, xx2, yy2, " HEP Depth2");
3863 leg.AddEntry(gr1);
3864 leg.AddEntry(gr2);
3865 leg.AddEntry(gr3);
3866 leg.AddEntry(gr4);
3867 leg.AddEntry(gr5);
3868 leg.SetMargin(0.3);
3869 leg.SetBorderSize(1);
3870 leg.DrawClone("Same");
3871
3872 }
3873
3874
3875
3876
3877 if (Method == 15) {
3878 for (int k1 = 0; k1 < NP; k1++) {
3879 YY[k1] = 0.;
3880 }
3881 c1->Range(0, 0, 25, 18);
3882 c1->SetFillColor(0);
3883 c1->Clear();
3884 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
3885 pad1->Draw();
3886 pad1->cd();
3887 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
3888
3889 pad1->DrawFrame(x1, y1, x2, y2);
3890
3891 pad1->GetFrame()->SetFillColor(0);
3892 pad1->Draw();
3893 gPad->SetGridy();
3894 gPad->SetGridx();
3895 cout << "****Draw 15 Dependencies: P r vs t for L1 and eta 21 22 23 25 27 28 **** " << endl;
3896
3897
3898
3899
3900
3901 for (int k1 = 0; k1 < NP; k1++) {
3902 YY[k1] = 0.;
3903 if (LASERruns[k1] > "305874") {
3904 for (int iii = 2; iii < 3; iii++) {
3905 for (int jj = 1; jj <= 13; jj++) {
3906 int LP1 = 0;
3907 if ((iii == 2 && jj > 0))
3908 LP1 = 1;
3909 if (LP1 == 1 && jj == 12) {
3910 YY[k1] = ZZP[k1][iii][jj];
3911 }
3912 }
3913 }
3914 } else {
3915 for (int iii = 1; iii < 2; iii++) {
3916 for (int jj = 1; jj <= 13; jj++) {
3917 int LP1 = 0;
3918 if ((iii == 1 && jj > 0))
3919 LP1 = 1;
3920 if (LP1 == 1 && jj == 12) {
3921 YY[k1] = ZZP[k1][iii][jj];
3922 }
3923 }
3924 }
3925 }
3926 }
3927 gr1 = new TGraph(NP, XX, YY);
3928 gr1->SetLineColor(1);
3929 gr1->SetMarkerColor(1);
3930 gr1->SetMarkerStyle(20);
3931 gr1->SetMarkerSize(1.2);
3932 gr1->SetTitle("ieta 28");
3933 gr1->SetFillStyle(0);
3934 gr1->Draw("P");
3935
3936
3937 for (int k1 = 0; k1 < NP; k1++) {
3938 YY[k1] = 0.;
3939 if (LASERruns[k1] > "305874") {
3940 for (int iii = 2; iii < 3; iii++) {
3941 for (int jj = 1; jj <= 13; jj++) {
3942 int LP1 = 0;
3943 if ((iii == 2 && jj > 0))
3944 LP1 = 1;
3945 if (LP1 == 1 && jj == 11) {
3946 YY[k1] = ZZP[k1][iii][jj];
3947 }
3948 }
3949 }
3950 } else {
3951 for (int iii = 1; iii < 2; iii++) {
3952 for (int jj = 1; jj <= 13; jj++) {
3953 int LP1 = 0;
3954 if ((iii == 1 && jj > 0))
3955 LP1 = 1;
3956 if (LP1 == 1 && jj == 11) {
3957 YY[k1] = ZZP[k1][iii][jj];
3958 }
3959 }
3960 }
3961 }
3962 }
3963 gr2 = new TGraph(NP, XX, YY);
3964 gr2->SetLineColor(2);
3965 gr2->SetMarkerColor(2);
3966 gr2->SetMarkerStyle(20);
3967 gr2->SetMarkerSize(1.2);
3968 gr2->SetTitle("ieta 27");
3969 gr2->SetFillStyle(0);
3970 gr2->Draw("P");
3971
3972 for (int k1 = 0; k1 < NP; k1++) {
3973 YY[k1] = 0.;
3974 if (LASERruns[k1] > "305874") {
3975 for (int iii = 2; iii < 3; iii++) {
3976 for (int jj = 1; jj <= 13; jj++) {
3977 int LP1 = 0;
3978 if ((iii == 2 && jj > 0))
3979 LP1 = 1;
3980 if (LP1 == 1 && jj == 9) {
3981 YY[k1] = ZZP[k1][iii][jj];
3982 }
3983 }
3984 }
3985 } else {
3986 for (int iii = 1; iii < 2; iii++) {
3987 for (int jj = 1; jj <= 13; jj++) {
3988 int LP1 = 0;
3989 if ((iii == 1 && jj > 0))
3990 LP1 = 1;
3991 if (LP1 == 1 && jj == 9) {
3992 YY[k1] = ZZP[k1][iii][jj];
3993 }
3994 }
3995 }
3996 }
3997 }
3998 gr3 = new TGraph(NP, XX, YY);
3999 gr3->SetLineColor(3);
4000 gr3->SetMarkerColor(3);
4001 gr3->SetMarkerStyle(20);
4002 gr3->SetMarkerSize(1.2);
4003 gr3->SetTitle("ieta 25");
4004 gr3->SetFillStyle(0);
4005 gr3->Draw("P");
4006
4007 for (int k1 = 0; k1 < NP; k1++) {
4008 YY[k1] = 0.;
4009 if (LASERruns[k1] > "305874") {
4010 for (int iii = 2; iii < 3; iii++) {
4011 for (int jj = 1; jj <= 13; jj++) {
4012 int LP1 = 0;
4013 if ((iii == 2 && jj > 0))
4014 LP1 = 1;
4015 if (LP1 == 1 && jj == 7) {
4016 YY[k1] = ZZP[k1][iii][jj];
4017 }
4018 }
4019 }
4020 } else {
4021 for (int iii = 1; iii < 2; iii++) {
4022 for (int jj = 1; jj <= 13; jj++) {
4023 int LP1 = 0;
4024 if ((iii == 1 && jj > 0))
4025 LP1 = 1;
4026 if (LP1 == 1 && jj == 7) {
4027 YY[k1] = ZZP[k1][iii][jj];
4028 }
4029 }
4030 }
4031 }
4032 }
4033 gr4 = new TGraph(NP, XX, YY);
4034 gr4->SetLineColor(4);
4035 gr4->SetMarkerColor(4);
4036 gr4->SetMarkerStyle(20);
4037 gr4->SetMarkerSize(1.2);
4038 gr4->SetTitle("ieta 23");
4039 gr4->SetFillStyle(0);
4040 gr4->Draw("P");
4041
4042 for (int k1 = 0; k1 < NP; k1++) {
4043 YY[k1] = 0.;
4044 if (LASERruns[k1] > "305874") {
4045 for (int iii = 2; iii < 3; iii++) {
4046 for (int jj = 1; jj <= 13; jj++) {
4047 int LP1 = 0;
4048 if ((iii == 2 && jj > 0))
4049 LP1 = 1;
4050 if (LP1 == 1 && jj == 6) {
4051 YY[k1] = ZZP[k1][iii][jj];
4052 }
4053 }
4054 }
4055 } else {
4056 for (int iii = 1; iii < 2; iii++) {
4057 for (int jj = 1; jj <= 13; jj++) {
4058 int LP1 = 0;
4059 if ((iii == 1 && jj > 0))
4060 LP1 = 1;
4061 if (LP1 == 1 && jj == 6) {
4062 YY[k1] = ZZP[k1][iii][jj];
4063 }
4064 }
4065 }
4066 }
4067 }
4068 gr5 = new TGraph(NP, XX, YY);
4069 gr5->SetLineColor(5);
4070 gr5->SetMarkerColor(5);
4071 gr5->SetMarkerStyle(20);
4072 gr5->SetMarkerSize(1.2);
4073 gr5->SetTitle("ieta 22");
4074 gr5->SetFillStyle(0);
4075 gr5->Draw("P");
4076
4077 for (int k1 = 0; k1 < NP; k1++) {
4078 YY[k1] = 0.;
4079 if (LASERruns[k1] > "305874") {
4080 for (int iii = 2; iii < 3; iii++) {
4081 for (int jj = 1; jj <= 13; jj++) {
4082 int LP1 = 0;
4083 if ((iii == 2 && jj > 0))
4084 LP1 = 1;
4085 if (LP1 == 1 && jj == 5) {
4086 YY[k1] = ZZP[k1][iii][jj];
4087 }
4088 }
4089 }
4090 } else {
4091 for (int iii = 1; iii < 2; iii++) {
4092 for (int jj = 1; jj <= 13; jj++) {
4093 int LP1 = 0;
4094 if ((iii == 1 && jj > 0))
4095 LP1 = 1;
4096 if (LP1 == 1 && jj == 5) {
4097 YY[k1] = ZZP[k1][iii][jj];
4098 }
4099 }
4100 }
4101 }
4102 }
4103 gr6 = new TGraph(NP, XX, YY);
4104 gr6->SetLineColor(6);
4105 gr6->SetMarkerColor(6);
4106 gr6->SetMarkerStyle(20);
4107 gr6->SetMarkerSize(1.2);
4108 gr6->SetTitle("ieta 21");
4109 gr6->SetFillStyle(0);
4110 gr6->Draw("P");
4111
4112
4113 TLegend leg(xx1, yy1, xx2, yy2, " HEP L1");
4114 leg.AddEntry(gr1);
4115 leg.AddEntry(gr2);
4116 leg.AddEntry(gr3);
4117 leg.AddEntry(gr4);
4118 leg.AddEntry(gr5);
4119 leg.AddEntry(gr6);
4120 leg.SetMargin(0.3);
4121 leg.SetBorderSize(1);
4122 leg.DrawClone("Same");
4123
4124 }
4125
4126
4127 if (Method == 16) {
4128 for (int k1 = 0; k1 < NP; k1++) {
4129 YY[k1] = 0.;
4130 }
4131 c1->Range(0, 0, 25, 18);
4132 c1->SetFillColor(0);
4133 c1->Clear();
4134 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
4135 pad1->Draw();
4136 pad1->cd();
4137 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
4138
4139 pad1->DrawFrame(x1, y1, x2, y2);
4140
4141 pad1->GetFrame()->SetFillColor(0);
4142 pad1->Draw();
4143 gPad->SetGridy();
4144 gPad->SetGridx();
4145 cout << "****Draw 16 Dependencies: P r vs t for L1 and eta 17 18 19 20 24 26 29 **** " << endl;
4146
4147
4148
4149
4150
4151 for (int k1 = 0; k1 < NP; k1++) {
4152 YY[k1] = 0.;
4153 if (LASERruns[k1] > "305874") {
4154 for (int iii = 2; iii < 3; iii++) {
4155 for (int jj = 1; jj <= 13; jj++) {
4156 int LP1 = 0;
4157 if ((iii == 2 && jj > 0))
4158 LP1 = 1;
4159 if (LP1 == 1 && jj == 13) {
4160 YY[k1] = ZZP[k1][iii][jj];
4161 }
4162 }
4163 }
4164 } else {
4165 for (int iii = 1; iii < 2; iii++) {
4166 for (int jj = 1; jj <= 13; jj++) {
4167 int LP1 = 0;
4168 if ((iii == 1 && jj > 0))
4169 LP1 = 1;
4170 if (LP1 == 1 && jj == 13) {
4171 YY[k1] = ZZP[k1][iii][jj];
4172 }
4173 }
4174 }
4175 }
4176 }
4177 gr1 = new TGraph(NP, XX, YY);
4178 gr1->SetLineColor(1);
4179 gr1->SetMarkerColor(1);
4180 gr1->SetMarkerStyle(20);
4181 gr1->SetMarkerSize(1.2);
4182 gr1->SetTitle("ieta 29");
4183 gr1->SetFillStyle(0);
4184 gr1->Draw("P");
4185
4186
4187 for (int k1 = 0; k1 < NP; k1++) {
4188 YY[k1] = 0.;
4189 if (LASERruns[k1] > "305874") {
4190 for (int iii = 2; iii < 3; iii++) {
4191 for (int jj = 1; jj <= 13; jj++) {
4192 int LP1 = 0;
4193 if ((iii == 2 && jj > 0))
4194 LP1 = 1;
4195 if (LP1 == 1 && jj == 10) {
4196 YY[k1] = ZZP[k1][iii][jj];
4197 }
4198 }
4199 }
4200 } else {
4201 for (int iii = 1; iii < 2; iii++) {
4202 for (int jj = 1; jj <= 13; jj++) {
4203 int LP1 = 0;
4204 if ((iii == 1 && jj > 0))
4205 LP1 = 1;
4206 if (LP1 == 1 && jj == 10) {
4207 YY[k1] = ZZP[k1][iii][jj];
4208 }
4209 }
4210 }
4211 }
4212 }
4213 gr2 = new TGraph(NP, XX, YY);
4214 gr2->SetLineColor(2);
4215 gr2->SetMarkerColor(2);
4216 gr2->SetMarkerStyle(20);
4217 gr2->SetMarkerSize(1.2);
4218 gr2->SetTitle("ieta 27");
4219 gr2->SetFillStyle(0);
4220 gr2->Draw("P");
4221
4222 for (int k1 = 0; k1 < NP; k1++) {
4223 YY[k1] = 0.;
4224 if (LASERruns[k1] > "305874") {
4225 for (int iii = 2; iii < 3; iii++) {
4226 for (int jj = 1; jj <= 13; jj++) {
4227 int LP1 = 0;
4228 if ((iii == 2 && jj > 0))
4229 LP1 = 1;
4230 if (LP1 == 1 && jj == 8) {
4231 YY[k1] = ZZP[k1][iii][jj];
4232 }
4233 }
4234 }
4235 } else {
4236 for (int iii = 1; iii < 2; iii++) {
4237 for (int jj = 1; jj <= 13; jj++) {
4238 int LP1 = 0;
4239 if ((iii == 1 && jj > 0))
4240 LP1 = 1;
4241 if (LP1 == 1 && jj == 8) {
4242 YY[k1] = ZZP[k1][iii][jj];
4243 }
4244 }
4245 }
4246 }
4247 }
4248 gr3 = new TGraph(NP, XX, YY);
4249 gr3->SetLineColor(3);
4250 gr3->SetMarkerColor(3);
4251 gr3->SetMarkerStyle(20);
4252 gr3->SetMarkerSize(1.2);
4253 gr3->SetTitle("ieta 25");
4254 gr3->SetFillStyle(0);
4255 gr3->Draw("P");
4256
4257
4258
4259
4260 for (int k1 = 0; k1 < NP; k1++) {
4261 YY[k1] = 0.;
4262 if (LASERruns[k1] > "305874") {
4263 for (int iii = 2; iii < 3; iii++) {
4264 for (int jj = 1; jj <= 13; jj++) {
4265 int LP1 = 0;
4266 if ((iii == 2 && jj > 0))
4267 LP1 = 1;
4268 if (LP1 == 1 && jj == 3) {
4269 YY[k1] = ZZP[k1][iii][jj];
4270 }
4271 }
4272 }
4273 } else {
4274 for (int iii = 1; iii < 2; iii++) {
4275 for (int jj = 1; jj <= 13; jj++) {
4276 int LP1 = 0;
4277 if ((iii == 1 && jj > 0))
4278 LP1 = 1;
4279 if (LP1 == 1 && jj == 3) {
4280 YY[k1] = ZZP[k1][iii][jj];
4281 }
4282 }
4283 }
4284 }
4285 }
4286 gr5 = new TGraph(NP, XX, YY);
4287 gr5->SetLineColor(5);
4288 gr5->SetMarkerColor(5);
4289 gr5->SetMarkerStyle(20);
4290 gr5->SetMarkerSize(1.2);
4291 gr5->SetTitle("ieta 22");
4292 gr5->SetFillStyle(0);
4293 gr5->Draw("P");
4294
4295
4296
4297
4298 for (int k1 = 0; k1 < NP; k1++) {
4299 YY[k1] = 0.;
4300 if (LASERruns[k1] > "305874") {
4301 for (int iii = 2; iii < 3; iii++) {
4302 for (int jj = 1; jj <= 13; jj++) {
4303 int LP1 = 0;
4304 if ((iii == 2 && jj > 0))
4305 LP1 = 1;
4306 if (LP1 == 1 && jj == 1) {
4307 YY[k1] = ZZP[k1][iii][jj];
4308 }
4309 }
4310 }
4311 } else {
4312 for (int iii = 1; iii < 2; iii++) {
4313 for (int jj = 1; jj <= 13; jj++) {
4314 int LP1 = 0;
4315 if ((iii == 1 && jj > 0))
4316 LP1 = 1;
4317 if (LP1 == 1 && jj == 1) {
4318 YY[k1] = ZZP[k1][iii][jj];
4319 }
4320 }
4321 }
4322 }
4323 }
4324 gr7 = new TGraph(NP, XX, YY);
4325 gr7->SetLineColor(7);
4326 gr7->SetMarkerColor(7);
4327 gr7->SetMarkerStyle(20);
4328 gr7->SetMarkerSize(1.2);
4329 gr7->SetTitle("ieta 21");
4330 gr7->SetFillStyle(0);
4331 gr7->Draw("P");
4332
4333
4334 TLegend leg(xx1, yy1, xx2, yy2, " HEP L1");
4335 leg.AddEntry(gr1);
4336 leg.AddEntry(gr2);
4337 leg.AddEntry(gr3);
4338 leg.AddEntry(gr4);
4339 leg.AddEntry(gr5);
4340 leg.SetMargin(0.3);
4341 leg.SetBorderSize(1);
4342 leg.DrawClone("Same");
4343 }
4344
4345
4346
4347 if (Method == 17) {
4348 for (int k1 = 0; k1 < NP; k1++) {
4349 YY[k1] = 0.;
4350 }
4351
4352
4353 cout << "**** Positive direction: take good runs eta 21 **** " << endl;
4354 c1->Range(0, 0, 25, 18);
4355 c1->SetFillColor(0);
4356 c1->Clear();
4357 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
4358 pad1->Draw();
4359 pad1->cd();
4360 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
4361
4362 pad1->DrawFrame(x1, y5, x2, y6);
4363
4364 pad1->GetFrame()->SetFillColor(0);
4365 pad1->Draw();
4366 gPad->SetGridy();
4367 gPad->SetGridx();
4368 cout << "****Draw 17 Dependencies: P r vs t for L7 and eta 21 22 23 25 27 28 **** " << endl;
4369
4370
4371
4372
4373
4374
4375 for (int k1 = 0; k1 < NP; k1++) {
4376 YY[k1] = 0.;
4377 if (LASERruns[k1] > "305874") {
4378 for (int iii = 3; iii < 6; iii++) {
4379 for (int jj = 1; jj < 13; jj++) {
4380 int LP7 = 0;
4381 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4382 LP7 = 1;
4383 if (LP7 == 1 && jj == 12) {
4384 YY[k1] = ZZP[k1][iii][jj];
4385 }
4386 }
4387 }
4388 } else {
4389 for (int iii = 2; iii < 4; iii++) {
4390 for (int jj = 1; jj < 13; jj++) {
4391 int LP7 = 0;
4392 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4393 LP7 = 1;
4394 if (LP7 == 1 && jj == 12) {
4395 YY[k1] = ZZP[k1][iii][jj];
4396 }
4397 }
4398 }
4399 }
4400 }
4401 gr1 = new TGraph(NP, XX, YY);
4402 gr1->SetLineColor(1);
4403 gr1->SetMarkerColor(1);
4404 gr1->SetMarkerStyle(20);
4405 gr1->SetTitle("ieta 28");
4406 gr1->SetFillStyle(0);
4407 gr1->SetMarkerSize(1.2);
4408 gr1->Draw("P");
4409
4410 for (int k1 = 0; k1 < NP; k1++) {
4411 YY[k1] = 0.;
4412 if (LASERruns[k1] > "305874") {
4413 for (int iii = 3; iii < 6; iii++) {
4414 for (int jj = 1; jj < 13; jj++) {
4415 int LP7 = 0;
4416 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4417 LP7 = 1;
4418 if (LP7 == 1 && jj == 11) {
4419 YY[k1] = ZZP[k1][iii][jj];
4420 }
4421 }
4422 }
4423 } else {
4424 for (int iii = 2; iii < 4; iii++) {
4425 for (int jj = 1; jj < 13; jj++) {
4426 int LP7 = 0;
4427 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4428 LP7 = 1;
4429 if (LP7 == 1 && jj == 11) {
4430 YY[k1] = ZZP[k1][iii][jj];
4431 }
4432 }
4433 }
4434 }
4435 }
4436 gr2 = new TGraph(NP, XX, YY);
4437 gr2->SetLineColor(2);
4438 gr2->SetMarkerColor(2);
4439 gr2->SetMarkerStyle(20);
4440 gr2->SetTitle("ieta 27");
4441 gr2->SetFillStyle(0);
4442 gr2->SetMarkerSize(1.2);
4443 gr2->Draw("P");
4444
4445 for (int k1 = 0; k1 < NP; k1++) {
4446 YY[k1] = 0.;
4447 if (LASERruns[k1] > "305874") {
4448 for (int iii = 3; iii < 6; iii++) {
4449 for (int jj = 1; jj < 13; jj++) {
4450 int LP7 = 0;
4451 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4452 LP7 = 1;
4453 if (LP7 == 1 && jj == 9) {
4454 YY[k1] = ZZP[k1][iii][jj];
4455 }
4456 }
4457 }
4458 } else {
4459 for (int iii = 2; iii < 4; iii++) {
4460 for (int jj = 1; jj < 13; jj++) {
4461 int LP7 = 0;
4462 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4463 LP7 = 1;
4464 if (LP7 == 1 && jj == 9) {
4465 YY[k1] = ZZP[k1][iii][jj];
4466 }
4467 }
4468 }
4469 }
4470 }
4471 gr3 = new TGraph(NP, XX, YY);
4472 gr3->SetLineColor(3);
4473 gr3->SetMarkerColor(3);
4474 gr3->SetMarkerStyle(20);
4475 gr3->SetTitle("ieta 25");
4476 gr3->SetFillStyle(0);
4477 gr3->SetMarkerSize(1.2);
4478 gr3->Draw("P");
4479
4480 for (int k1 = 0; k1 < NP; k1++) {
4481 YY[k1] = 0.;
4482 if (LASERruns[k1] > "305874") {
4483 for (int iii = 3; iii < 6; iii++) {
4484 for (int jj = 1; jj < 13; jj++) {
4485 int LP7 = 0;
4486 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4487 LP7 = 1;
4488 if (LP7 == 1 && jj == 7) {
4489 YY[k1] = ZZP[k1][iii][jj];
4490 }
4491 }
4492 }
4493 } else {
4494 for (int iii = 2; iii < 4; iii++) {
4495 for (int jj = 1; jj < 13; jj++) {
4496 int LP7 = 0;
4497 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4498 LP7 = 1;
4499 if (LP7 == 1 && jj == 7) {
4500 YY[k1] = ZZP[k1][iii][jj];
4501 }
4502 }
4503 }
4504 }
4505 }
4506 gr4 = new TGraph(NP, XX, YY);
4507 gr4->SetLineColor(4);
4508 gr4->SetMarkerColor(4);
4509 gr4->SetMarkerStyle(20);
4510 gr4->SetMarkerSize(1.2);
4511 gr4->SetTitle("ieta 23");
4512 gr4->SetFillStyle(0);
4513 gr4->SetMarkerSize(1.2);
4514 gr4->Draw("P");
4515
4516 for (int k1 = 0; k1 < NP; k1++) {
4517 YY[k1] = 0.;
4518 if (LASERruns[k1] > "305874") {
4519 for (int iii = 3; iii < 6; iii++) {
4520 for (int jj = 1; jj < 13; jj++) {
4521 int LP7 = 0;
4522 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4523 LP7 = 1;
4524 if (LP7 == 1 && jj == 6) {
4525 YY[k1] = ZZP[k1][iii][jj];
4526 }
4527 }
4528 }
4529 } else {
4530 for (int iii = 2; iii < 4; iii++) {
4531 for (int jj = 1; jj < 13; jj++) {
4532 int LP7 = 0;
4533 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4534 LP7 = 1;
4535 if (LP7 == 1 && jj == 6) {
4536 YY[k1] = ZZP[k1][iii][jj];
4537 }
4538 }
4539 }
4540 }
4541 }
4542 gr5 = new TGraph(NP, XX, YY);
4543 gr5->SetLineColor(5);
4544 gr5->SetMarkerColor(5);
4545 gr5->SetMarkerStyle(20);
4546 gr5->SetTitle("ieta 22");
4547 gr5->SetFillStyle(0);
4548 gr5->SetMarkerSize(1.2);
4549 gr5->Draw("P");
4550
4551 for (int k1 = 0; k1 < NP; k1++) {
4552 YY[k1] = 0.;
4553 if (LASERruns[k1] > "305874") {
4554 for (int iii = 3; iii < 6; iii++) {
4555 for (int jj = 1; jj < 13; jj++) {
4556 int LP7 = 0;
4557 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4558 LP7 = 1;
4559 if (LP7 == 1 && jj == 5) {
4560 YY[k1] = ZZP[k1][iii][jj];
4561 }
4562 }
4563 }
4564 } else {
4565 for (int iii = 2; iii < 4; iii++) {
4566 for (int jj = 1; jj < 13; jj++) {
4567 int LP7 = 0;
4568 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4569 LP7 = 1;
4570 if (LP7 == 1 && jj == 5) {
4571 YY[k1] = ZZP[k1][iii][jj];
4572 }
4573 }
4574 }
4575 }
4576
4577
4578
4579
4580 }
4581 gr6 = new TGraph(NP, XX, YY);
4582 gr6->SetLineColor(6);
4583 gr6->SetMarkerColor(6);
4584 gr6->SetMarkerStyle(20);
4585 gr6->SetMarkerSize(1.2);
4586 gr6->SetTitle("ieta 21");
4587 gr6->SetFillStyle(0);
4588 gr6->Draw("P");
4589
4590 TLegend leg(xx1, yy1, xx2, yy2, " HEP L7");
4591 leg.AddEntry(gr1);
4592 leg.AddEntry(gr2);
4593 leg.AddEntry(gr3);
4594 leg.AddEntry(gr4);
4595 leg.AddEntry(gr5);
4596 leg.AddEntry(gr6);
4597 leg.SetMargin(0.3);
4598 leg.SetBorderSize(1);
4599 leg.DrawClone("Same");
4600 }
4601
4602
4603 if (Method == 18) {
4604 for (int k1 = 0; k1 < NP; k1++) {
4605 YY[k1] = 0.;
4606 }
4607 c1->Range(0, 0, 25, 18);
4608 c1->SetFillColor(0);
4609 c1->Clear();
4610 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
4611 pad1->Draw();
4612 pad1->cd();
4613 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
4614
4615 pad1->DrawFrame(x1, y5, x2, y6);
4616
4617 pad1->GetFrame()->SetFillColor(0);
4618 pad1->Draw();
4619 gPad->SetGridy();
4620 gPad->SetGridx();
4621 cout << "****Draw 18 Dependencies: P r vs t for L7 and eta 17 18 19 20 24 26 29 **** " << endl;
4622
4623
4624
4625
4626
4627
4628 for (int k1 = 0; k1 < NP; k1++) {
4629 YY[k1] = 0.;
4630 if (LASERruns[k1] > "305874") {
4631 for (int iii = 3; iii < 6; iii++) {
4632 for (int jj = 1; jj < 13; jj++) {
4633 int LP7 = 0;
4634 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4635 LP7 = 1;
4636 if (LP7 == 1 && jj == 10) {
4637 YY[k1] = ZZP[k1][iii][jj];
4638 }
4639 }
4640 }
4641 } else {
4642 for (int iii = 2; iii < 4; iii++) {
4643 for (int jj = 1; jj < 13; jj++) {
4644 int LP7 = 0;
4645 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4646 LP7 = 1;
4647 if (LP7 == 1 && jj == 10) {
4648 YY[k1] = ZZP[k1][iii][jj];
4649 }
4650 }
4651 }
4652 }
4653 }
4654 gr1 = new TGraph(NP, XX, YY);
4655 gr1->SetLineColor(1);
4656 gr1->SetMarkerColor(1);
4657 gr1->SetMarkerStyle(20);
4658 gr1->SetTitle("ieta 26");
4659 gr1->SetFillStyle(0);
4660 gr1->SetMarkerSize(1.2);
4661 gr1->Draw("P");
4662
4663 for (int k1 = 0; k1 < NP; k1++) {
4664 YY[k1] = 0.;
4665 if (LASERruns[k1] > "305874") {
4666 for (int iii = 3; iii < 6; iii++) {
4667 for (int jj = 1; jj < 13; jj++) {
4668 int LP7 = 0;
4669 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4670 LP7 = 1;
4671 if (LP7 == 1 && jj == 8) {
4672 YY[k1] = ZZP[k1][iii][jj];
4673 }
4674 }
4675 }
4676 } else {
4677 for (int iii = 2; iii < 4; iii++) {
4678 for (int jj = 1; jj < 13; jj++) {
4679 int LP7 = 0;
4680 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4681 LP7 = 1;
4682 if (LP7 == 1 && jj == 8) {
4683 YY[k1] = ZZP[k1][iii][jj];
4684 }
4685 }
4686 }
4687 }
4688 }
4689 gr2 = new TGraph(NP, XX, YY);
4690 gr2->SetLineColor(2);
4691 gr2->SetMarkerColor(2);
4692 gr2->SetMarkerStyle(20);
4693 gr2->SetTitle("ieta 24");
4694 gr2->SetFillStyle(0);
4695 gr2->SetMarkerSize(1.2);
4696 gr2->Draw("P");
4697
4698 for (int k1 = 0; k1 < NP; k1++) {
4699 YY[k1] = 0.;
4700 if (LASERruns[k1] > "305874") {
4701 for (int iii = 3; iii < 6; iii++) {
4702 for (int jj = 1; jj < 13; jj++) {
4703 int LP7 = 0;
4704 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4705 LP7 = 1;
4706 if (LP7 == 1 && jj == 4) {
4707 YY[k1] = ZZP[k1][iii][jj];
4708 }
4709 }
4710 }
4711 } else {
4712 for (int iii = 2; iii < 4; iii++) {
4713 for (int jj = 1; jj < 13; jj++) {
4714 int LP7 = 0;
4715 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4716 LP7 = 1;
4717 if (LP7 == 1 && jj == 4) {
4718 YY[k1] = ZZP[k1][iii][jj];
4719 }
4720 }
4721 }
4722 }
4723 }
4724 gr3 = new TGraph(NP, XX, YY);
4725 gr3->SetLineColor(3);
4726 gr3->SetMarkerColor(3);
4727 gr3->SetMarkerStyle(20);
4728 gr3->SetTitle("ieta 20");
4729 gr3->SetFillStyle(0);
4730 gr3->SetMarkerSize(1.2);
4731 gr3->Draw("P");
4732
4733 for (int k1 = 0; k1 < NP; k1++) {
4734 YY[k1] = 0.;
4735 if (LASERruns[k1] > "305874") {
4736 for (int iii = 3; iii < 6; iii++) {
4737 for (int jj = 1; jj < 13; jj++) {
4738 int LP7 = 0;
4739 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4740 LP7 = 1;
4741 if (LP7 == 1 && jj == 3) {
4742 YY[k1] = ZZP[k1][iii][jj];
4743 }
4744 }
4745 }
4746 } else {
4747 for (int iii = 2; iii < 4; iii++) {
4748 for (int jj = 1; jj < 13; jj++) {
4749 int LP7 = 0;
4750 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4751 LP7 = 1;
4752 if (LP7 == 1 && jj == 3) {
4753 YY[k1] = ZZP[k1][iii][jj];
4754 }
4755 }
4756 }
4757 }
4758 }
4759 gr4 = new TGraph(NP, XX, YY);
4760 gr4->SetLineColor(4);
4761 gr4->SetMarkerColor(4);
4762 gr4->SetMarkerStyle(20);
4763 gr4->SetTitle("ieta 19");
4764 gr4->SetFillStyle(0);
4765 gr4->SetMarkerSize(1.2);
4766 gr4->Draw("P");
4767
4768 for (int k1 = 0; k1 < NP; k1++) {
4769 YY[k1] = 0.;
4770 if (LASERruns[k1] > "305874") {
4771 for (int iii = 3; iii < 6; iii++) {
4772 for (int jj = 1; jj < 13; jj++) {
4773 int LP7 = 0;
4774 if ((iii == 3 && jj == 1) || (iii == 4 && jj > 1 && jj < 10) || (iii == 5 && jj > 9 && jj < 13))
4775 LP7 = 1;
4776 if (LP7 == 1 && jj == 2) {
4777 YY[k1] = ZZP[k1][iii][jj];
4778 }
4779 }
4780 }
4781 } else {
4782 for (int iii = 2; iii < 4; iii++) {
4783 for (int jj = 1; jj < 13; jj++) {
4784 int LP7 = 0;
4785 if ((iii == 2 && jj > 0 && jj < 11) || (iii == 3 && jj > 10 && jj < 13))
4786 LP7 = 1;
4787 if (LP7 == 1 && jj == 2) {
4788 YY[k1] = ZZP[k1][iii][jj];
4789 }
4790 }
4791 }
4792 }
4793 }
4794 gr5 = new TGraph(NP, XX, YY);
4795 gr5->SetLineColor(5);
4796 gr5->SetMarkerColor(5);
4797 gr5->SetMarkerStyle(20);
4798 gr5->SetTitle("ieta 18");
4799 gr5->SetFillStyle(0);
4800 gr5->SetMarkerSize(1.2);
4801 gr5->Draw("P");
4802
4803
4804
4805
4806 TLegend leg(xx1, yy1, xx2, yy2, " HEP L7");
4807 leg.AddEntry(gr1);
4808 leg.AddEntry(gr2);
4809 leg.AddEntry(gr3);
4810 leg.AddEntry(gr4);
4811 leg.AddEntry(gr5);
4812 leg.SetMargin(0.3);
4813 leg.SetBorderSize(1);
4814 leg.DrawClone("Same");
4815 }
4816
4817
4818 if (Method == 21) {
4819 cout << "****Draw 21 Dependencies: P r vs t for Depth3 and some eta 27 28 *** " << endl;
4820
4821 gStyle->SetOptTitle(kFALSE);
4822 gStyle->SetPalette(kSolar);
4823 for (int k1 = 0; k1 < NP; k1++) {
4824 YY[k1] = 0.;
4825 }
4826 c1->Range(0, 0, 25, 18);
4827
4828 c1->SetFillColor(0);
4829 c1->Clear();
4830 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
4831 pad1->Draw();
4832 pad1->cd();
4833 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
4834
4835 pad1->DrawFrame(x1, y1, x2, y2);
4836
4837 pad1->GetFrame()->SetFillColor(0);
4838 pad1->Draw();
4839 gPad->SetGridy();
4840 gPad->SetGridx();
4841
4842
4843
4844 for (int k1 = 0; k1 < NP; k1++) {
4845 YY[k1] = ZZP[k1][3][11];
4846 YYE[k1] = 0.01 * YY[k1];
4847 }
4848
4849 ge1 = new TGraphErrors(NP, XX, YY, XXE, YYE);
4850 ge1->SetLineColor(1);
4851 ge1->SetMarkerColor(1);
4852 ge1->SetMarkerStyle(20);
4853 ge1->SetMarkerSize(1.2);
4854 ge1->SetTitle("ieta 27");
4855 ge1->SetName("gr1");
4856 ge1->SetFillStyle(0);
4857 ge1->Draw("P");
4858 TF1 *myfit = new TF1("myfit", "[0] + [1]*x**{[2]}", 1., 2.);
4859 myfit->SetParName(0, "c0");
4860 myfit->SetParName(1, "c1");
4861 myfit->SetParName(2, "c2");
4862 myfit->SetParameter(0, 0.01);
4863 myfit->SetParameter(1, 1.1);
4864 myfit->SetParameter(2, 0.12);
4865 gStyle->SetOptFit(1101);
4866 myfit->SetLineColor(1);
4867 ge1->Fit("myfit");
4868
4869 for (int k1 = 0; k1 < NP; k1++) {
4870 YY[k1] = ZZP[k1][3][12];
4871 YYE[k1] = 0.01 * YY[k1];
4872 }
4873
4874 ge2 = new TGraphErrors(NP, XX, YY, XXE, YYE);
4875 ge2->SetLineColor(2);
4876 ge2->SetMarkerColor(2);
4877 ge2->SetMarkerStyle(20);
4878 ge2->SetMarkerSize(1.2);
4879 ge2->SetTitle("ieta 28");
4880 ge2->SetName("ge2");
4881 ge2->SetFillStyle(0);
4882 ge2->Draw("P");
4883
4884
4885 TF1 *myfit2 = new TF1("myfit2", "[0]*sin(x) + [1]*exp(-[2]*x)", 0, 2);
4886 myfit2->SetParName(0, "c0");
4887 myfit2->SetParName(1, "c1");
4888 myfit2->SetParName(2, "c2");
4889 myfit2->SetParameter(0, 0.01);
4890 myfit2->SetParameter(1, 1.1);
4891 myfit2->SetParameter(2, 0.12);
4892 gStyle->SetOptFit(1101);
4893 myfit2->SetLineColor(2);
4894 ge2->Fit("myfit2");
4895
4896
4897 TLegend leg(xx1, yy1, xx2, yy2, " HEP Depth3");
4898
4899 leg.AddEntry(ge1);
4900 leg.AddEntry(ge2);
4901 leg.SetMargin(0.3);
4902 leg.SetBorderSize(1);
4903 leg.DrawClone("Same");
4904
4905
4906
4907
4908 }
4909
4910
4911
4912 if (Method == 22) {
4913 cout << "****Draw 22 Dependencies: P: draw <Rijk>_ev.run for all runs *** " << endl;
4914 gPad->SetGridy();
4915 gPad->SetGridx();
4916 gPad->SetLogy();
4917 gStyle->SetOptStat(1110000);
4918 rwithphi->SetMarkerStyle(20);
4919 rwithphi->SetMarkerSize(0.4);
4920 rwithphi->GetYaxis()->SetLabelSize(0.04);
4921 rwithphi->SetTitle("all runs & channels");
4922 rwithphi->SetXTitle("HEP <Rijk>");
4923 rwithphi->SetYTitle("N");
4924 rwithphi->SetMarkerColor(4);
4925 rwithphi->SetLineColor(4);
4926
4927 rwithphi->Draw("PL");
4928 }
4929
4930
4931 if (Method == 24) {
4932 cout << "****Draw 24 Dependencies: M: draw <Rijk>_ev.run for all runs *** " << endl;
4933 gPad->SetGridy();
4934 gPad->SetGridx();
4935 gPad->SetLogy();
4936 gStyle->SetOptStat(1110000);
4937 rwithphi->SetMarkerStyle(20);
4938 rwithphi->SetMarkerSize(0.4);
4939 rwithphi->GetYaxis()->SetLabelSize(0.04);
4940 rwithphi->SetTitle("all runs & channels");
4941 rwithphi->SetXTitle("HEM <Rijk>");
4942 rwithphi->SetYTitle("N");
4943 rwithphi->SetMarkerColor(4);
4944 rwithphi->SetLineColor(4);
4945
4946 rwithphi->Draw("PL");
4947 }
4948
4949
4950
4951 if (Method == 23) {
4952 cout << "****Draw 23 Dependencies: P:gains LED runs used *** " << endl;
4953 gPad->SetGridy();
4954 gPad->SetGridx();
4955 gPad->SetLogy();
4956 gStyle->SetOptStat(1110000);
4957 gainkoeff->SetMarkerStyle(20);
4958 gainkoeff->SetMarkerSize(0.4);
4959 gainkoeff->GetYaxis()->SetLabelSize(0.04);
4960 gainkoeff->SetTitle("all channels");
4961 gainkoeff->SetXTitle("gains");
4962 gainkoeff->SetYTitle("N");
4963 gainkoeff->SetMarkerColor(4);
4964 gainkoeff->SetLineColor(4);
4965
4966 gainkoeff->Draw("PL");
4967 }
4968
4969 if (Method == 25) {
4970 cout << "****Draw 25 R(Theta) Dependencies for last LASER runs for M&P *** " << endl;
4971
4972
4973 c1->Range(0, 0, 25, 18);
4974 c1->SetFillColor(0);
4975 c1->Clear();
4976 TPad *pad1 = new TPad("pad1", "This is pad1", 0.02, 0.02, 0.98, 0.98, 0);
4977 pad1->Draw();
4978 pad1->cd();
4979 pad1->Range(-0.255174, -19.25, 2.29657, -6.75);
4980
4981
4982
4983 pad1->DrawFrame(16., 0.70, 30., 1.05);
4984
4985
4986 pad1->GetFrame()->SetFillColor(0);
4987 pad1->Draw();
4988 gPad->SetGridy();
4989 gPad->SetGridx();
4990
4991
4992
4993
4994 gr1 = new TGraph(13, THETA, RRM);
4995 gr1->SetLineColor(2);
4996 gr1->SetMarkerColor(2);
4997 gr1->SetMarkerStyle(25);
4998 gr1->SetMarkerSize(1.5);
4999 gr1->SetTitle("HEM");
5000 gr1->SetFillStyle(0);
5001 gr1->Draw("PL");
5002 gr2 = new TGraph(13, THETA, RRP);
5003 gr2->SetLineColor(4);
5004 gr2->SetMarkerColor(4);
5005 gr2->SetMarkerStyle(25);
5006 gr2->SetMarkerSize(1.5);
5007 gr2->SetTitle("HEP");
5008 gr2->SetFillStyle(0);
5009 gr2->Draw("PL");
5010
5011 cout << "*******THETA7 0 = " << THETA7[0] << " 1 = " << THETA7[1] << " 2= " << THETA7[2] << endl;
5012 cout << "*******RRM7 0 = " << RRM7[0] << " 1 = " << RRM7[1] << " 2= " << RRM7[2] << endl;
5013
5014 cout << "**************************************************************************************** " << endl;
5015
5016 cout << "******THETA7 9 = " << THETA7[9] << " 10 = " << THETA7[10] << " 11= " << THETA7[11] << endl;
5017 cout << "*******RRM7 9 = " << RRM7[9] << " 10 = " << RRM7[10] << " 11= " << RRM7[11] << endl;
5018
5019 gr3 = new TGraph(11, THETA7, RRM7);
5020 gr3->SetLineColor(2);
5021 gr3->SetMarkerColor(2);
5022 gr3->SetMarkerStyle(26);
5023 gr3->SetMarkerSize(1.5);
5024 gr3->SetTitle("HEM");
5025 gr3->SetFillStyle(0);
5026 gr3->Draw("PL");
5027 gr4 = new TGraph(11, THETA7, RRP7);
5028 gr4->SetLineColor(4);
5029 gr4->SetMarkerColor(4);
5030 gr4->SetMarkerStyle(26);
5031 gr4->SetMarkerSize(1.5);
5032 gr4->SetTitle("HEP");
5033 gr4->SetFillStyle(0);
5034 gr4->Draw("PL");
5035
5036
5037 TLegend leg(xx1, yy1 + 0.1, xx2, yy2 + 0.1, " L1 & L7 ");
5038 leg.AddEntry(gr1);
5039 leg.AddEntry(gr2);
5040 leg.AddEntry(gr3);
5041 leg.AddEntry(gr4);
5042 leg.SetMargin(0.3);
5043 leg.SetBorderSize(1);
5044 leg.DrawClone("Same");
5045
5046
5047
5048
5049
5050
5051
5052
5053
5054
5055
5056
5057
5058
5059
5060
5061
5062
5063
5064 }
5065
5066
5067 if (Method == 26) {
5068 cout << "****Draw 26 Dependencies: M: draw <RRijk>_ev.run for all runs *** " << endl;
5069 gPad->SetGridy();
5070 gPad->SetGridx();
5071 gPad->SetLogy();
5072 gStyle->SetOptStat(1110000);
5073 rrwithphiM->SetMarkerStyle(20);
5074 rrwithphiM->SetMarkerSize(0.4);
5075 rrwithphiM->GetYaxis()->SetLabelSize(0.04);
5076 rrwithphiM->SetTitle("all runs & channels");
5077 rrwithphiM->SetXTitle("HEM <RRijk>");
5078 rrwithphiM->SetYTitle("N");
5079 rrwithphiM->SetMarkerColor(4);
5080 rrwithphiM->SetLineColor(4);
5081
5082 rrwithphiM->Draw("PL");
5083 }
5084
5085
5086 if (Method == 27) {
5087 cout << "****Draw 27 Dependencies: P: draw <RRijk>_ev.run for all runs *** " << endl;
5088 gPad->SetGridy();
5089 gPad->SetGridx();
5090 gPad->SetLogy();
5091 gStyle->SetOptStat(1110000);
5092 rrwithphiP->SetMarkerStyle(20);
5093 rrwithphiP->SetMarkerSize(0.4);
5094 rrwithphiP->GetYaxis()->SetLabelSize(0.04);
5095 rrwithphiP->SetTitle("all runs & channels");
5096 rrwithphiP->SetXTitle("HEP <RRijk>");
5097 rrwithphiP->SetYTitle("N");
5098 rrwithphiP->SetMarkerColor(4);
5099 rrwithphiP->SetLineColor(4);
5100
5101 rrwithphiP->Draw("PL");
5102 }
5103
5104
5105 if (Method == 28) {
5106 cout << "****Draw 28 Dependencies: M: Mean *** " << endl;
5107 gPad->SetGridy();
5108 gPad->SetGridx();
5109 gPad->SetLogy();
5110 gStyle->SetOptStat(1110000);
5111 c1->SetFillColor(0);
5112
5113 Ceff1M->SetMarkerStyle(20);
5114 Ceff1M->SetMarkerSize(0.4);
5115 Ceff1M->GetYaxis()->SetLabelSize(0.04);
5116 Ceff1M->SetTitle("M: all runs & channels");
5117 Ceff1M->SetXTitle("M: Mean <Aijk>");
5118 Ceff1M->SetYTitle("N");
5119 Ceff1M->SetMarkerColor(4);
5120 Ceff1M->SetLineColor(4);
5121 Ceff1M->SetMinimum(0.7);
5122
5123 Ceff1M->Draw("PL");
5124 }
5125
5126
5127 if (Method == 29) {
5128 cout << "****Draw 29 Dependencies: P:Mean *** " << endl;
5129 gPad->SetGridy();
5130 gPad->SetGridx();
5131 gPad->SetLogy();
5132 gStyle->SetOptStat(1110000);
5133 c1->SetFillColor(0);
5134
5135 Ceff1M->SetMarkerStyle(20);
5136 Ceff1M->SetMarkerSize(0.4);
5137 Ceff1M->GetYaxis()->SetLabelSize(0.04);
5138 Ceff1M->SetTitle("P: all runs & channels");
5139 Ceff1M->SetXTitle("P: Mean <Aijk>");
5140 Ceff1M->SetYTitle("N");
5141 Ceff1M->SetMarkerColor(4);
5142 Ceff1M->SetLineColor(4);
5143 Ceff1M->SetMinimum(0.7);
5144
5145 Ceff1M->Draw("PL");
5146 }
5147
5148
5149 if (Method == 30) {
5150 cout << "****Draw 30 Dependencies:M PHI *** " << endl;
5151 phidepdrawing30->Divide(phidependenceM, phidependenceM0, 1, 1, "B");
5152 phidepdrawing30->Sumw2();
5153 gPad->SetGridy();
5154 gPad->SetGridx();
5155 gPad->SetLogy(0);
5156 phidepdrawing30->SetMarkerStyle(20);
5157 phidepdrawing30->SetMarkerSize(1.4);
5158 phidepdrawing30->GetYaxis()->SetLabelSize(0.04);
5159 phidepdrawing30->SetXTitle("Phi");
5160 phidepdrawing30->SetYTitle("RR");
5161 phidepdrawing30->SetMarkerColor(2);
5162 phidepdrawing30->SetLineColor(2);
5163 phidepdrawing30->Draw("Error");
5164
5165 }
5166
5167
5168 if (Method == 31) {
5169 cout << "****Draw 31 Dependencies:P PHI *** " << endl;
5170 phidepdrawing31->Divide(phidependenceP, phidependenceP0, 1, 1, "B");
5171 phidepdrawing31->Sumw2();
5172 gPad->SetGridy();
5173 gPad->SetGridx();
5174 gPad->SetLogy(0);
5175 phidepdrawing31->SetMarkerStyle(20);
5176 phidepdrawing31->SetMarkerSize(1.4);
5177 phidepdrawing31->GetYaxis()->SetLabelSize(0.04);
5178 phidepdrawing31->SetXTitle("Phi");
5179 phidepdrawing31->SetYTitle("RR");
5180 phidepdrawing31->SetMarkerColor(4);
5181 phidepdrawing31->SetLineColor(4);
5182 phidepdrawing31->Draw("Error");
5183
5184 }
5185
5186
5187 if (Method == 32) {
5188 cout << "****Draw 32 Dependencies: M:RMS *** " << endl;
5189 gPad->SetGridy();
5190 gPad->SetGridx();
5191 gPad->SetLogy();
5192 gStyle->SetOptStat(1110000);
5193 c1->SetFillColor(0);
5194 Ceff1R->SetMarkerStyle(20);
5195 Ceff1R->SetMarkerSize(0.4);
5196 Ceff1R->GetYaxis()->SetLabelSize(0.04);
5197 Ceff1R->SetTitle("M: all runs & channels");
5198 Ceff1R->SetXTitle("M: RMS <Aijk>");
5199 Ceff1R->SetYTitle("N");
5200 Ceff1R->SetMarkerColor(4);
5201 Ceff1R->SetLineColor(4);
5202 Ceff1R->SetMinimum(0.7);
5203 Ceff1R->Draw("PL");
5204 }
5205
5206
5207 if (Method == 33) {
5208 cout << "****Draw 33 Dependencies: P:RMS *** " << endl;
5209 gPad->SetGridy();
5210 gPad->SetGridx();
5211 gPad->SetLogy();
5212 gStyle->SetOptStat(1110000);
5213 c1->SetFillColor(0);
5214 Ceff1R->SetMarkerStyle(20);
5215 Ceff1R->SetMarkerSize(0.4);
5216 Ceff1R->GetYaxis()->SetLabelSize(0.04);
5217 Ceff1R->SetTitle("P: all runs & channels");
5218 Ceff1R->SetXTitle("P: RMS <Aijk>");
5219 Ceff1R->SetYTitle("N");
5220 Ceff1R->SetMarkerColor(4);
5221 Ceff1R->SetLineColor(4);
5222 Ceff1R->SetMinimum(0.7);
5223 Ceff1R->Draw("PL");
5224 }
5225
5226
5227
5228 if (Method == 34) {
5229 cout << "****Draw 34 P:IntLumi fb-1 vs. Days *** " << endl;
5230 gPad->SetGridy();
5231 gPad->SetGridx();
5232 gPad->SetLogy(0);
5233 gStyle->SetOptStat(000000000);
5234 lumiplot->GetOption();
5235 lumiplot->SetMarkerStyle(20);
5236 lumiplot->SetMarkerSize(0.4);
5237 lumiplot->GetYaxis()->SetLabelSize(0.04);
5238 lumiplot->SetXTitle("days");
5239 lumiplot->SetYTitle("IntLumi fb-1");
5240 lumiplot->SetMarkerColor(4);
5241 lumiplot->SetLineColor(4);
5242 lumiplot->SetMinimum(0.);
5243 lumiplot->SetMaximum(170.);
5244
5245 lumiplot->Draw("hist PL");
5246 }
5247
5248
5249 if (Method == 35) {
5250 cout << "****Draw 35 P:IntLumi fb-1 vs. Days (up to 171) *** " << endl;
5251 gPad->SetGridy();
5252 gPad->SetGridx();
5253 gPad->SetLogy(0);
5254 gStyle->SetOptStat(000000000);
5255 lumiplot1->GetOption();
5256 lumiplot1->SetMarkerStyle(20);
5257 lumiplot1->SetMarkerSize(0.4);
5258 lumiplot1->GetYaxis()->SetLabelSize(0.04);
5259 lumiplot1->SetXTitle("days");
5260 lumiplot1->SetYTitle("IntLumi fb-1");
5261 lumiplot1->SetMarkerColor(4);
5262 lumiplot1->SetLineColor(4);
5263 lumiplot1->SetMinimum(0.);
5264 lumiplot1->SetMaximum(45.);
5265
5266 lumiplot1->Draw("hist PL");
5267 }
5268
5269
5270 if (Method == 36) {
5271 cout << "****Draw 36 P:IntLumi fb-1 vs. Days (2017 only) *** " << endl;
5272 gPad->SetGridy();
5273 gPad->SetGridx();
5274 gPad->SetLogy(0);
5275 gStyle->SetOptStat(000000000);
5276 lumiplot2->GetOption();
5277 lumiplot2->SetMarkerStyle(20);
5278 lumiplot2->SetMarkerSize(0.4);
5279 lumiplot2->GetYaxis()->SetLabelSize(0.04);
5280 lumiplot2->SetXTitle("days");
5281 lumiplot2->SetYTitle("IntLumi fb-1");
5282 lumiplot2->SetMarkerColor(4);
5283 lumiplot2->SetLineColor(4);
5284 lumiplot2->SetMinimum(0.);
5285 lumiplot2->SetMaximum(55.);
5286
5287 lumiplot2->Draw("hist PL");
5288 }
5289
5290
5291 if (Method == 37) {
5292 cout << "****Draw 37 P:IntLumi fb-1 vs. Days (2018 only) *** " << endl;
5293 gPad->SetGridy();
5294 gPad->SetGridx();
5295 gPad->SetLogy(0);
5296 gStyle->SetOptStat(000000000);
5297 lumiplot3->GetOption();
5298 lumiplot3->SetMarkerStyle(20);
5299 lumiplot3->SetMarkerSize(0.4);
5300 lumiplot3->GetYaxis()->SetLabelSize(0.04);
5301 lumiplot3->SetXTitle("days");
5302 lumiplot3->SetYTitle("IntLumi fb-1");
5303 lumiplot3->SetMarkerColor(4);
5304 lumiplot3->SetLineColor(4);
5305 lumiplot3->SetMinimum(0.);
5306 lumiplot3->SetMaximum(75.);
5307
5308 lumiplot3->Draw("hist PL");
5309 }
5310
5311
5312
5313
5314
5315
5316
5317
5318
5319
5320
5321
5322
5323
5324
5325 TString fdir = "amt2019raddam/";
5326
5327 TString plotName = fdir + MethodName[Method] + ".png";
5328 c1->SaveAs(plotName);
5329 cout << "****PLOT " << Method << " drawn **** " << endl;
5330
5331 }
5332 ++Method;
5333 }
5334
5335 }
5336
5337
5338
5339
5340 TString fdir = "amt2019raddam/";
5341
5342
5343 phidepdrawing30->Draw("Error");
5344
5345 phidepdrawing31->Draw("ErrorSame");
5346 TString plotName = fdir + "PLOT100.png";
5347 c1->SaveAs(plotName);
5348 cout << "****PLOT " << 100 << " drawn **** " << endl;
5349
5350 gStyle->SetOptStat(1110000);
5351 phiRRM->SetMarkerStyle(20);
5352 phiRRM->SetMarkerSize(1.4);
5353 phiRRM->GetYaxis()->SetLabelSize(0.04);
5354 phiRRM->SetXTitle("RR");
5355 phiRRM->SetYTitle("N");
5356 phiRRM->SetMarkerColor(2);
5357 phiRRM->SetLineColor(2);
5358 phiRRM->Draw("Error");
5359 phiRRP->SetMarkerStyle(20);
5360 phiRRP->SetMarkerSize(1.4);
5361 phiRRP->GetYaxis()->SetLabelSize(0.04);
5362 phiRRP->SetXTitle("RR");
5363 phiRRP->SetYTitle("N");
5364 phiRRP->SetMarkerColor(4);
5365 phiRRP->SetLineColor(4);
5366 phiRRP->Draw("ErrorSame");
5367 plotName = fdir + "PLOT101.png";
5368 c1->SaveAs(plotName);
5369 cout << "****PLOT " << 101 << " drawn **** " << endl;
5370
5371
5372
5373
5374
5375
5376
5377
5378
5379
5380
5381
5382
5383
5384
5385
5386 gSystem->Exit(0);
5387
5388 }
5389