diff --git a/.gitignore b/.gitignore index a78f54b..521c740 100644 --- a/.gitignore +++ b/.gitignore @@ -7,6 +7,7 @@ EventBuilder* *.txt *.err *.seq +*.png Mapper AnasenMS @@ -26,4 +27,3 @@ qqq_gains_det3.dat qqq_relative_gains.dat Armory/CorrelateQQQ.h QQQStage2.C -qqqIvERun395.png diff --git a/ELoss/Eloss_17F b/ELoss/Eloss_17F new file mode 100644 index 0000000..e516a41 --- /dev/null +++ b/ELoss/Eloss_17F @@ -0,0 +1,268 @@ +0.017002095 0.00676442 +0.03400419 0.0135288 +0.085010476 0.032227 +0.136016762 0.0489265 +0.170020952 0.0592235 +0.255031428 0.0823584 +0.340041904 0.102374 +0.42505238 0.120086 +0.510062856 0.135943 +0.595073332 0.150461 +0.680083808 0.164011 +0.765094284 0.176885 +0.85010476 0.189282 +0.935115236 0.201176 +1.020125712 0.212518 +1.190146664 0.233962 +1.360167616 0.254045 +1.530188568 0.27299 +1.70020952 0.291072 +2.040251424 0.325807 +2.380293328 0.359536 +2.720335232 0.392589 +3.060377136 0.425159 +3.40041904 0.457286 +4.2505238 0.536817 +5.10062856 0.615601 +5.95073332 0.693052 +6.80083808 0.769499 +7.65094284 0.845324 +8.5010476 0.920793 +10.20125712 1.07377 +11.90146664 1.23229 +13.60167616 1.39414 +15.30188568 1.55782 +17.0020952 1.72459 +21.252619 2.15811 +25.5031428 2.61883 +29.7536666 3.10864 +34.0041904 3.62853 +38.2547142 4.17854 +42.505238 4.75843 +51.0062856 6.00863 +59.5073332 7.38089 +68.0083808 8.87668 +76.5094284 10.4957 +85.010476 12.2377 +93.5115236 14.1011 +102.0125712 16.0845 +110.5136188 18.1875 +119.0146664 20.4097 +127.515714 22.7505 +136.0167616 25.2096 +144.5178092 27.7867 +153.0188568 30.4814 +161.5199044 33.2935 +170.020952 36.2195 +187.0230472 42.214 +204.0251424 48.6246 +221.0272376 55.449 +238.0293328 62.6845 +255.031428 70.3284 +272.0335232 78.3781 +289.0356184 86.8308 +306.0377136 95.684 +323.0398088 104.935 +340.041904 114.581 +357.0439992 124.619 +374.0460944 135.048 +391.0481896 145.863 +408.0502848 157.064 +425.05238 168.647 +459.0565704 192.95 +493.0607608 218.753 +527.0649512 246.038 +561.0691416 274.787 +595.073332 304.983 +646.0796176 352.95 +697.0859032 404.079 +748.0921888 458.315 +799.0984744 515.608 +850.10476 575.909 +935.115236 682.969 +1020.125712 798.048 +1105.136188 920.949 +1190.146664 1051.49 +1275.15714 1189.49 +1360.167616 1334.78 +1445.178092 1487.21 +1530.188568 1646.62 +1615.199044 1812.86 +1700.20952 1985.8 +1785.219996 2165.29 +1870.230472 2351.2 +1955.240948 2543.42 +2040.251424 2741.82 +2125.2619 2946.27 +2210.272376 3156.67 +2295.282852 3372.91 +2380.293328 3594.88 +2465.303804 3822.47 +2550.31428 4055.59 +2635.324756 4294.14 +2720.335232 4538.01 +2805.345708 4787.13 +2890.356184 5041.4 +2975.36666 5300.73 +3060.377136 5565.04 +3145.387612 5834.24 +3230.398088 6108.26 +3315.408564 6387.01 +3400.41904 6670.42 +3570.439992 7250.92 +3740.460944 7849.17 +3910.481896 8464.62 +4080.502848 9096.74 +4250.5238 9745.02 +4420.544752 10409 +4590.565704 11088.2 +4760.586656 11782.1 +4930.607608 12490.4 +5100.62856 13212.7 +5270.649512 13948.5 +5440.670464 14697.4 +5610.691416 15459.2 +5780.712368 16233.4 +5950.73332 17019.7 +6120.754272 17817.8 +6290.775224 18627.3 +6460.796176 19448 +6630.817128 20279.6 +6800.83808 21121.8 +6970.859032 21974.2 +7140.879984 22836.7 +7310.900936 23709 +7480.921888 24590.9 +7650.94284 25482 +7820.963792 26382.2 +7990.984744 27291.2 +8161.005696 28208.8 +8331.026648 29134.9 +8501.0476 30069.1 +8841.089504 31961.3 +9181.131408 33884 +9521.173312 35835.8 +9861.215216 37815.1 +10201.25712 39820.8 +10541.29902 41851.6 +10881.34093 43906.4 +11221.38283 45984.1 +11561.42474 48083.6 +11901.46664 50203.9 +12241.50854 52344.2 +12581.55045 54503.4 +12921.59235 56680.8 +13261.63426 58875.5 +13601.67616 61086.7 +13941.71806 63313.8 +14281.75997 65556.1 +14621.80187 67812.7 +14961.84378 70083.2 +15301.88568 72366.9 +15641.92758 74663.1 +15981.96949 76971.4 +16322.01139 79291.2 +16662.0533 81622 +17002.0952 83963.3 +17342.1371 86314.6 +17682.17901 88675.4 +18022.22091 91045.4 +18362.26282 93424 +18702.30472 95811 +19042.34662 98205.9 +19382.38853 100608 +19722.43043 103018 +20062.47234 105434 +20402.51424 107857 +20742.55614 110286 +21082.59805 112721 +21422.63995 115162 +21762.68186 117608 +22102.72376 120059 +22442.76566 122514 +22782.80757 124975 +23122.84947 127439 +23462.89138 129908 +23802.93328 132380 +24142.97518 134856 +24483.01709 137335 +24823.05899 139817 +25163.1009 142303 +25503.1428 144791 +25843.1847 147281 +26183.22661 149774 +26523.26851 152270 +26863.31042 154767 +27203.35232 157266 +27543.39422 159767 +27883.43613 162270 +28223.47803 164774 +28563.51994 167279 +28903.56184 169786 +29753.6666 176056 +30603.77136 182331 +31453.87612 188609 +32303.98088 194888 +33154.08564 201166 +34004.1904 207442 +35704.39992 219984 +37404.60944 232503 +39104.81896 244992 +40805.02848 257444 +42505.238 269853 +44205.44752 282214 +45905.65704 294523 +47605.86656 306775 +49306.07608 318968 +51006.2856 331099 +54406.70464 355166 +57807.12368 378963 +61207.54272 402479 +64607.96176 425709 +68008.3808 448649 +71408.79984 471298 +74809.21888 493655 +78209.63792 515723 +81610.05696 537503 +85010.476 558999 +89260.9998 585475 +93511.5236 611522 +102012.5712 662361 +110513.6188 711594 +119014.6664 759299 +127515.714 805555 +136016.7616 850436 +144517.8092 894016 +153018.8568 936363 +161519.9044 977541 +170020.952 1.02E+06 +187023.0472 1.09E+06 +204025.1424 1.17E+06 +221027.2376 1.24E+06 +238029.3328 1.30E+06 +255031.428 1.37E+06 +272033.5232 1.43E+06 +289035.6184 1.49E+06 +306037.7136 1.54E+06 +323039.8088 1.60E+06 +340041.904 1.65E+06 +357043.9992 1.70E+06 +374046.0944 1.75E+06 +391048.1896 1.80E+06 +408050.2848 1.84E+06 +425052.38 1.89E+06 +442054.4752 1.93E+06 +459056.5704 1.98E+06 +476058.6656 2.02E+06 +493060.7608 2.06E+06 +510062.856 2.10E+06 +544067.0464 2.17E+06 +578071.2368 2.24E+06 +612075.4272 2.31E+06 +646079.6176 2.38E+06 +680083.808 2.44E+06 +714087.9984 2.50E+06 +748092.1888 2.56E+06 +782096.3792 2.61E+06 +816100.5696 2.67E+06 +850104.76 2.72E+06 diff --git a/ELoss/Eloss_27Al b/ELoss/Eloss_27Al new file mode 100644 index 0000000..9e9f489 --- /dev/null +++ b/ELoss/Eloss_27Al @@ -0,0 +1,268 @@ +0.026981539 1087.74 +0.053963078 2175.49 +0.134907695 5438.8 +0.215852312 8557.19 +0.26981539 10540.7 +0.404723085 15079.2 +0.53963078 19011.9 +0.674538475 22467.9 +0.80944617 25531.3 +0.944353865 28303.5 +1.07926156 30859.2 +1.214169255 33253.7 +1.34907695 35529.1 +1.483984645 37685.6 +1.61889234 39719.2 +1.88870773 43505.2 +2.15852312 46980.1 +2.42833851 50199.6 +2.6981539 53228.6 +3.23778468 58941.2 +3.77741546 64372.8 +4.31704624 69606.4 +4.85667702 74683.1 +5.3963078 79618 +6.74538475 91580.1 +8.0944617 103093 +9.44353865 114105 +10.7926156 124746 +12.14169255 135119 +13.4907695 145269 +16.1889234 165367 +18.8870773 185754 +21.5852312 206324 +24.2833851 226926 +26.981539 247688 +33.72692375 300998 +40.4723085 357029 +47.21769325 416359 +53.963078 478956 +60.70846275 544719 +67.4538475 613499 +80.944617 759812 +94.4353865 917732 +107.926156 1.09E+06 +121.4169255 1.27E+06 +134.907695 1.46E+06 +148.3984645 1.67E+06 +161.889234 1.88E+06 +175.3800035 2.11E+06 +188.870773 2.35E+06 +202.3615425 2.60E+06 +215.852312 2.86E+06 +229.3430815 3.13E+06 +242.833851 3.41E+06 +256.3246205 3.71E+06 +269.81539 4.01E+06 +296.796929 4.64E+06 +323.778468 5.30E+06 +350.760007 6.00E+06 +377.741546 6.74E+06 +404.723085 7.53E+06 +431.704624 8.35E+06 +458.686163 9.21E+06 +485.667702 1.01E+07 +512.649241 1.10E+07 +539.63078 1.20E+07 +566.612319 1.30E+07 +593.593858 1.41E+07 +620.575397 1.52E+07 +647.556936 1.63E+07 +674.538475 1.75E+07 +728.501553 1.99E+07 +782.464631 2.25E+07 +836.427709 2.52E+07 +890.390787 2.81E+07 +944.353865 3.11E+07 +1025.298482 3.59E+07 +1106.243099 4.10E+07 +1187.187716 4.64E+07 +1268.132333 5.21E+07 +1349.07695 5.81E+07 +1483.984645 6.88E+07 +1618.89234 8.02E+07 +1753.800035 9.24E+07 +1888.70773 1.05E+08 +2023.615425 1.19E+08 +2158.52312 1.34E+08 +2293.430815 1.49E+08 +2428.33851 1.65E+08 +2563.246205 1.81E+08 +2698.1539 1.98E+08 +2833.061595 2.16E+08 +2967.96929 2.35E+08 +3102.876985 2.54E+08 +3237.78468 2.73E+08 +3372.692375 2.94E+08 +3507.60007 3.15E+08 +3642.507765 3.36E+08 +3777.41546 3.58E+08 +3912.323155 3.81E+08 +4047.23085 4.04E+08 +4182.138545 4.27E+08 +4317.04624 4.52E+08 +4451.953935 4.76E+08 +4586.86163 5.02E+08 +4721.769325 5.27E+08 +4856.67702 5.54E+08 +4991.584715 5.80E+08 +5126.49241 6.08E+08 +5261.400105 6.35E+08 +5396.3078 6.64E+08 +5666.12319 7.21E+08 +5935.93858 7.81E+08 +6205.75397 8.42E+08 +6475.56936 9.04E+08 +6745.38475 9.69E+08 +7015.20014 1.03E+09 +7285.01553 1.10E+09 +7554.83092 1.17E+09 +7824.64631 1.24E+09 +8094.4617 1.31E+09 +8364.27709 1.39E+09 +8634.09248 1.46E+09 +8903.90787 1.54E+09 +9173.72326 1.61E+09 +9443.53865 1.69E+09 +9713.35404 1.77E+09 +9983.16943 1.85E+09 +10252.98482 1.93E+09 +10522.80021 2.02E+09 +10792.6156 2.10E+09 +11062.43099 2.18E+09 +11332.24638 2.27E+09 +11602.06177 2.36E+09 +11871.87716 2.44E+09 +12141.69255 2.53E+09 +12411.50794 2.62E+09 +12681.32333 2.71E+09 +12951.13872 2.80E+09 +13220.95411 2.89E+09 +13490.7695 2.99E+09 +14030.40028 3.18E+09 +14570.03106 3.37E+09 +15109.66184 3.56E+09 +15649.29262 3.76E+09 +16188.9234 3.96E+09 +16728.55418 4.16E+09 +17268.18496 4.36E+09 +17807.81574 4.57E+09 +18347.44652 4.78E+09 +18887.0773 4.99E+09 +19426.70808 5.20E+09 +19966.33886 5.41E+09 +20505.96964 5.63E+09 +21045.60042 5.85E+09 +21585.2312 6.07E+09 +22124.86198 6.29E+09 +22664.49276 6.51E+09 +23204.12354 6.74E+09 +23743.75432 6.96E+09 +24283.3851 7.19E+09 +24823.01588 7.42E+09 +25362.64666 7.65E+09 +25902.27744 7.88E+09 +26441.90822 8.11E+09 +26981.539 8.34E+09 +27521.16978 8.57E+09 +28060.80056 8.81E+09 +28600.43134 9.04E+09 +29140.06212 9.28E+09 +29679.6929 9.52E+09 +30219.32368 9.75E+09 +30758.95446 9.99E+09 +31298.58524 1.02E+10 +31838.21602 1.05E+10 +32377.8468 1.07E+10 +32917.47758 1.10E+10 +33457.10836 1.12E+10 +33996.73914 1.14E+10 +34536.36992 1.17E+10 +35076.0007 1.19E+10 +35615.63148 1.22E+10 +36155.26226 1.24E+10 +36694.89304 1.27E+10 +37234.52382 1.29E+10 +37774.1546 1.31E+10 +38313.78538 1.34E+10 +38853.41616 1.36E+10 +39393.04694 1.39E+10 +39932.67772 1.41E+10 +40472.3085 1.44E+10 +41011.93928 1.46E+10 +41551.57006 1.49E+10 +42091.20084 1.51E+10 +42630.83162 1.54E+10 +43170.4624 1.56E+10 +43710.09318 1.59E+10 +44249.72396 1.61E+10 +44789.35474 1.64E+10 +45328.98552 1.66E+10 +45868.6163 1.69E+10 +47217.69325 1.75E+10 +48566.7702 1.81E+10 +49915.84715 1.87E+10 +51264.9241 1.94E+10 +52614.00105 2.00E+10 +53963.078 2.06E+10 +56661.2319 2.19E+10 +59359.3858 2.31E+10 +62057.5397 2.43E+10 +64755.6936 2.56E+10 +67453.8475 2.68E+10 +70152.0014 2.80E+10 +72850.1553 2.93E+10 +75548.3092 3.05E+10 +78246.4631 3.17E+10 +80944.617 3.29E+10 +86340.9248 3.53E+10 +91737.2326 3.76E+10 +97133.5404 4.00E+10 +102529.8482 4.23E+10 +107926.156 4.46E+10 +113322.4638 4.68E+10 +118718.7716 4.90E+10 +124115.0794 5.12E+10 +129511.3872 5.34E+10 +134907.695 5.55E+10 +141653.0798 5.82E+10 +148398.4645 6.07E+10 +161889.234 6.58E+10 +175380.0035 7.07E+10 +188870.773 7.54E+10 +202361.5425 8.00E+10 +215852.312 8.45E+10 +229343.0815 8.88E+10 +242833.851 9.30E+10 +256324.6205 9.71E+10 +269815.39 1.01E+11 +296796.929 1.09E+11 +323778.468 1.16E+11 +350760.007 1.23E+11 +377741.546 1.30E+11 +404723.085 1.36E+11 +431704.624 1.42E+11 +458686.163 1.48E+11 +485667.702 1.53E+11 +512649.241 1.59E+11 +539630.78 1.64E+11 +566612.319 1.69E+11 +593593.858 1.74E+11 +620575.397 1.79E+11 +647556.936 1.83E+11 +674538.475 1.88E+11 +701520.014 1.92E+11 +728501.553 1.96E+11 +755483.092 2.00E+11 +782464.631 2.04E+11 +809446.17 2.08E+11 +863409.248 2.16E+11 +917372.326 2.23E+11 +971335.404 2.29E+11 +1025298.482 2.36E+11 +1079261.56 2.42E+11 +1133224.638 2.48E+11 +1187187.716 2.54E+11 +1241150.794 2.59E+11 +1295113.872 2.65E+11 +1349076.95 2.70E+11 \ No newline at end of file diff --git a/ELoss/Eloss_alpha b/ELoss/Eloss_alpha new file mode 100644 index 0000000..1575be9 --- /dev/null +++ b/ELoss/Eloss_alpha @@ -0,0 +1,345 @@ +0.053103695 11959.9 +0.116582129 19201.1 +0.180060563 24466.2 +0.243538998 28696 +0.307017432 32395.8 +0.370495866 35815.8 +0.4339739 39049.4 +0.497451133 42211.1 +0.560932369 45315.6 +0.624409603 48392.9 +0.687886836 51473.4 +0.751368072 54567.7 +0.814845305 57693.8 +0.878322539 60883.5 +0.941803775 64073.2 +1.005281008 67274.6 +1.068758242 70625.6 +1.132235475 73976.6 +1.195716711 77327.7 +1.259193944 80859 +1.322671178 84405.9 +1.386152414 87952.9 +1.449629647 91666.6 +1.513106881 95430.9 +1.576588117 99195.1 +1.64006535 103102 +1.703542583 107097 +1.767019817 111093 +1.830501053 115200 +1.893978286 119436 +1.95745552 123672 +2.020936756 128023 +2.084413989 132632 +2.147891222 137241 +2.211372458 141849 +2.274849692 146458 +2.338326925 151067 +2.401804159 155677 +2.465285395 160793 +2.528762628 165908 +2.592239861 171023 +2.655721097 176139 +2.719198331 181254 +2.782675564 186369 +2.8461568 191841 +2.909634034 197468 +2.973111267 203094 +3.0365885 208721 +3.100069736 214347 +3.16354697 219973 +3.227024203 225800 +3.290505439 231938 +3.353982673 238075 +3.417459906 244212 +3.480941142 250349 +3.544418375 256486 +3.607895609 262668 +3.671372842 269313 +3.734854078 275958 +3.798331312 282603 +3.861808545 289248 +3.925289781 295893 +3.988767014 302538 +4.05223224 309870 +4.115713476 317394 +4.179194712 324917 +4.242675948 332441 +4.306157184 339965 +4.36963842 347488 +4.433119656 355012 +4.496600892 362535 +4.560082128 370059 +4.623563364 377583 +4.6870446 385106 +4.750525836 392630 +4.813967046 400154 +4.877448282 407677 +4.940929518 415201 +5.004410754 422747 +5.06789199 431506 +5.131373226 440265 +5.194854462 449023 +5.258335698 457782 +5.321816934 466541 +5.38529817 475299 +5.448779406 484058 +5.512260642 492817 +5.575741878 501575 +5.639223114 510334 +5.70270435 519093 +5.76614556 527851 +5.829626796 536610 +5.893108032 545369 +5.956589268 554127 +6.020070504 563195 +6.08355174 573167 +6.147032976 583139 +6.210514212 593110 +6.273995448 603082 +6.337476684 613054 +6.40095792 623026 +6.464439156 632997 +6.527920392 642969 +6.591401628 652941 +6.654842838 662913 +6.718324074 672885 +6.78180531 682856 +6.845286546 692828 +6.908767782 702800 +6.972249018 712772 +7.035730254 723331 +7.09921149 734498 +7.162692726 745666 +7.226173962 756833 +7.289655198 768000 +7.353136434 779168 +7.41661767 790335 +7.480098906 801502 +7.543540116 812670 +7.607021352 823837 +7.670502588 835004 +7.733983824 846172 +7.79746506 857339 +7.860946296 868507 +7.924427532 879674 +7.987908768 890841 +8.051390004 902868 +8.11487124 915217 +8.178352476 927566 +8.241833712 939915 +8.305314948 952263 +8.368796184 964612 +8.43227742 976961 +8.49571863 989310 +8.559199866 1.00E+06 +8.622681102 1.01E+06 +8.686162338 1.03E+06 +8.749643574 1.04E+06 +8.81312481 1.05E+06 +8.876606046 1.06E+06 +8.940087282 1.08E+06 +9.003568518 1.09E+06 +9.067049754 1.10E+06 +9.13053099 1.12E+06 +9.194012226 1.13E+06 +9.257493462 1.14E+06 +9.320974698 1.16E+06 +9.384415908 1.17E+06 +9.447897144 1.18E+06 +9.51137838 1.20E+06 +9.574859616 1.21E+06 +9.638340852 1.22E+06 +9.701822088 1.24E+06 +9.765303324 1.25E+06 +9.82878456 1.26E+06 +9.892265796 1.28E+06 +9.955747032 1.29E+06 +10.01922827 1.30E+06 +10.0827095 1.32E+06 +10.14619074 1.34E+06 +10.20967198 1.35E+06 +10.27311319 1.37E+06 +10.33659442 1.38E+06 +10.40007566 1.40E+06 +10.46355689 1.41E+06 +10.52703813 1.43E+06 +10.59051937 1.44E+06 +10.6540006 1.46E+06 +10.71748184 1.47E+06 +10.78096307 1.49E+06 +10.84444431 1.50E+06 +10.90792555 1.52E+06 +10.97140678 1.53E+06 +11.03488802 1.55E+06 +11.09836925 1.56E+06 +11.16185049 1.58E+06 +11.2252917 1.59E+06 +11.28877294 1.61E+06 +11.35225417 1.62E+06 +11.41573541 1.64E+06 +11.47921664 1.66E+06 +11.54269788 1.67E+06 +11.60617912 1.69E+06 +11.66966035 1.70E+06 +11.73314159 1.72E+06 +11.79662282 1.73E+06 +11.86010406 1.75E+06 +11.9235853 1.76E+06 +11.98706653 1.78E+06 +12.05054777 1.79E+06 +12.11398898 1.81E+06 +12.17747021 1.83E+06 +12.24095145 1.85E+06 +12.30443269 1.86E+06 +12.36791392 1.88E+06 +12.43139516 1.90E+06 +12.49487639 1.92E+06 +12.55835763 1.93E+06 +12.62183887 1.95E+06 +12.6853201 1.97E+06 +12.74880134 1.99E+06 +12.81228257 2.00E+06 +12.87576381 2.02E+06 +12.93924505 2.04E+06 +13.00268626 2.06E+06 +13.06616749 2.07E+06 +13.12964873 2.09E+06 +13.19312996 2.11E+06 +13.2566112 2.13E+06 +13.32009244 2.14E+06 +13.38357367 2.16E+06 +13.44705491 2.18E+06 +13.51053614 2.20E+06 +13.57401738 2.22E+06 +13.63749862 2.23E+06 +13.70097985 2.25E+06 +13.76446109 2.27E+06 +13.82794232 2.29E+06 +13.89142356 2.30E+06 +13.95486477 2.32E+06 +14.01834601 2.34E+06 +14.08182724 2.36E+06 +14.14530848 2.38E+06 +14.20878971 2.40E+06 +14.27227095 2.42E+06 +14.33575219 2.44E+06 +14.39923342 2.46E+06 +14.46271466 2.48E+06 +14.52619589 2.50E+06 +14.58967713 2.52E+06 +14.65315837 2.54E+06 +14.7166396 2.56E+06 +14.78012084 2.58E+06 +14.84356205 2.60E+06 +14.90704328 2.62E+06 +14.97052452 2.63E+06 +15.03400576 2.65E+06 +15.09748699 2.67E+06 +15.16096823 2.69E+06 +15.22444946 2.71E+06 +15.2879307 2.73E+06 +15.35141194 2.75E+06 +15.41489317 2.77E+06 +15.47837441 2.79E+06 +15.54185564 2.81E+06 +15.60533688 2.83E+06 +15.66881812 2.85E+06 +15.73225933 2.87E+06 +15.79574056 2.89E+06 +15.8592218 2.91E+06 +15.92270303 2.93E+06 +15.98618427 2.95E+06 +16.04966551 2.97E+06 +16.11314674 2.99E+06 +16.17662798 3.02E+06 +16.24010921 3.04E+06 +16.30359045 3.06E+06 +16.36707169 3.08E+06 +16.43055292 3.10E+06 +16.49403416 3.13E+06 +16.55751539 3.15E+06 +16.6209566 3.17E+06 +16.68443784 3.19E+06 +16.74791908 3.21E+06 +16.81140031 3.24E+06 +16.87488155 3.26E+06 +16.93836278 3.28E+06 +17.00184402 3.30E+06 +17.06532526 3.32E+06 +17.12880649 3.35E+06 +17.19228773 3.37E+06 +17.25576896 3.39E+06 +17.3192502 3.41E+06 +17.38273144 3.43E+06 +17.44621267 3.46E+06 +17.50969391 3.48E+06 +17.57313512 3.50E+06 +17.63661635 3.52E+06 +17.70009759 3.54E+06 +17.76357883 3.57E+06 +17.82706006 3.59E+06 +17.8905413 3.61E+06 +17.95402253 3.63E+06 +18.01750377 3.65E+06 +18.08098501 3.68E+06 +18.14446624 3.70E+06 +18.20794748 3.73E+06 +18.27142871 3.75E+06 +18.33490995 3.78E+06 +18.39839119 3.80E+06 +18.4618324 3.82E+06 +18.52531363 3.85E+06 +18.58879487 3.87E+06 +18.6522761 3.90E+06 +18.71575734 3.92E+06 +18.77923858 3.95E+06 +18.84271981 3.97E+06 +18.90620105 3.99E+06 +18.96968228 4.02E+06 +19.03316352 4.04E+06 +19.09664476 4.07E+06 +19.16012599 4.09E+06 +19.22360723 4.11E+06 +19.28708846 4.14E+06 +19.35052967 4.16E+06 +19.41401091 4.19E+06 +19.47749215 4.21E+06 +19.54097338 4.24E+06 +19.60445462 4.26E+06 +19.66793585 4.28E+06 +19.73141709 4.31E+06 +19.79489833 4.33E+06 +19.85837956 4.36E+06 +19.9218608 4.38E+06 +19.98534203 4.41E+06 +20.04882327 4.43E+06 +20.11230451 4.46E+06 +20.17578574 4.48E+06 +20.23926698 4.51E+06 +20.30270819 4.54E+06 +20.36618942 4.56E+06 +20.42967066 4.59E+06 +20.4931519 4.62E+06 +20.55663313 4.64E+06 +20.62011437 4.67E+06 +20.6835956 4.69E+06 +20.74707684 4.72E+06 +20.81055808 4.75E+06 +20.87403931 4.77E+06 +20.93752055 4.80E+06 +21.00100178 4.83E+06 +21.06448302 4.85E+06 +21.12796426 4.88E+06 +21.19140547 4.91E+06 +21.2548867 4.93E+06 +21.31836794 4.96E+06 +21.38184917 4.98E+06 +21.44533041 5.01E+06 +21.50881165 5.04E+06 +21.57229288 5.06E+06 +21.63577412 5.09E+06 +21.69925535 5.12E+06 +21.76273659 5.14E+06 +21.82621783 5.17E+06 +21.88969906 5.20E+06 diff --git a/ELoss/Eloss_p b/ELoss/Eloss_p new file mode 100644 index 0000000..3e227da --- /dev/null +++ b/ELoss/Eloss_p @@ -0,0 +1,268 @@ +0.00100784 0.00455783 +0.00201568 0.00911567 +0.0050392 0.0197891 +0.00806272 0.0284591 +0.0100784 0.0336328 +0.0151176 0.0449321 +0.0201568 0.0545806 +0.025196 0.0631952 +0.0302352 0.071111 +0.0352744 0.0785363 +0.0403136 0.0856119 +0.0453528 0.0924397 +0.050392 0.0990962 +0.0554312 0.105641 +0.0604704 0.112122 +0.0705488 0.125041 +0.0806272 0.138091 +0.0907056 0.151435 +0.100784 0.16519 +0.1209408 0.194217 +0.1410976 0.225513 +0.1612544 0.25917 +0.1814112 0.295144 +0.201568 0.333333 +0.25196 0.4377 +0.302352 0.55345 +0.352744 0.679377 +0.403136 0.814707 +0.453528 0.958947 +0.50392 1.11177 +0.604704 1.44231 +0.705488 1.80513 +0.806272 2.19946 +0.907056 2.62475 +1.00784 3.08057 +1.2598 4.35139 +1.51176 5.80613 +1.76372 7.44108 +2.01568 9.25315 +2.26764 11.2397 +2.5196 13.3984 +3.02352 18.2246 +3.52744 23.7179 +4.03136 29.8671 +4.53528 36.6625 +5.0392 44.0961 +5.54312 52.1607 +6.04704 60.8501 +6.55096 70.1587 +7.05488 80.0814 +7.5588 90.6136 +8.06272 101.751 +8.56664 113.49 +9.07056 125.827 +9.57448 138.759 +10.0784 152.267 +11.08624 180.081 +12.09408 209.992 +13.10192 241.974 +14.10976 276.006 +15.1176 312.066 +16.12544 350.134 +17.13328 390.189 +18.14112 432.214 +19.14896 476.19 +20.1568 522.101 +21.16464 569.93 +22.17248 619.661 +23.18032 671.279 +24.18816 724.77 +25.196 780.118 +27.21168 896.334 +29.22736 1019.82 +31.24304 1150.48 +33.25872 1288.2 +35.2744 1432.91 +38.29792 1662.88 +41.32144 1908.07 +44.34496 2168.23 +47.36848 2443.1 +50.392 2732.43 +55.4312 3246.18 +60.4704 3798.47 +65.5096 4388.33 +70.5488 5014.87 +75.588 5677.25 +80.6272 6374.65 +85.6664 7106.3 +90.7056 7871.46 +95.7448 8669.43 +100.784 9499.52 +105.8232 10361.1 +110.8624 11253.5 +115.9016 12176.2 +120.9408 13128.5 +125.98 14109.9 +131.0192 15119.8 +136.0584 16157.8 +141.0976 17223.3 +146.1368 18315.7 +151.176 19434.7 +156.2152 20579.8 +161.2544 21750.4 +166.2936 22946.2 +171.3328 24166.7 +176.372 25411.6 +181.4112 26680.3 +186.4504 27972.5 +191.4896 29287.8 +196.5288 30625.8 +201.568 31986.2 +211.6464 34772.7 +221.7248 37644.4 +231.8032 40598.6 +241.8816 43632.8 +251.96 46744.7 +262.0384 49931.8 +272.1168 53192 +282.1952 56523.1 +292.2736 59923.1 +302.352 63390 +312.4304 66921.8 +322.5088 70516.9 +332.5872 74173.4 +342.6656 77889.7 +352.744 81664 +362.8224 85494.9 +372.9008 89380.9 +382.9792 93320.4 +393.0576 97312.1 +403.136 101355 +413.2144 105447 +423.2928 109587 +433.3712 113774 +443.4496 118007 +453.528 122284 +463.6064 126605 +473.6848 130969 +483.7632 135373 +493.8416 139818 +503.92 144303 +524.0768 153386 +544.2336 162615 +564.3904 171984 +584.5472 181485 +604.704 191113 +624.8608 200861 +645.0176 210724 +665.1744 220697 +685.3312 230775 +705.488 240953 +725.6448 251226 +745.8016 261591 +765.9584 272043 +786.1152 282578 +806.272 293192 +826.4288 303883 +846.5856 314646 +866.7424 325478 +886.8992 336377 +907.056 347338 +927.2128 358361 +947.3696 369441 +967.5264 380577 +987.6832 391765 +1007.84 403003 +1027.9968 414290 +1048.1536 425622 +1068.3104 436998 +1088.4672 448416 +1108.624 459874 +1128.7808 471370 +1148.9376 482902 +1169.0944 494468 +1189.2512 506067 +1209.408 517698 +1229.5648 529358 +1249.7216 541047 +1269.8784 552762 +1290.0352 564502 +1310.192 576267 +1330.3488 588055 +1350.5056 599864 +1370.6624 611694 +1390.8192 623543 +1410.976 635410 +1431.1328 647295 +1451.2896 659195 +1471.4464 671111 +1491.6032 683041 +1511.76 694984 +1531.9168 706940 +1552.0736 718907 +1572.2304 730885 +1592.3872 742872 +1612.544 754869 +1632.7008 766874 +1652.8576 778887 +1673.0144 790907 +1693.1712 802933 +1713.328 814964 +1763.72 845063 +1814.112 875184 +1864.504 905318 +1914.896 935457 +1965.288 965594 +2015.68 995721 +2116.464 1.06E+06 +2217.248 1.12E+06 +2318.032 1.18E+06 +2418.816 1.24E+06 +2519.6 1.30E+06 +2620.384 1.35E+06 +2721.168 1.41E+06 +2821.952 1.47E+06 +2922.736 1.53E+06 +3023.52 1.59E+06 +3225.088 1.70E+06 +3426.656 1.82E+06 +3628.224 1.93E+06 +3829.792 2.04E+06 +4031.36 2.15E+06 +4232.928 2.26E+06 +4434.496 2.37E+06 +4636.064 2.48E+06 +4837.632 2.58E+06 +5039.2 2.68E+06 +5291.16 2.81E+06 +5543.12 2.94E+06 +6047.04 3.18E+06 +6550.96 3.42E+06 +7054.88 3.64E+06 +7558.8 3.87E+06 +8062.72 4.08E+06 +8566.64 4.29E+06 +9070.56 4.49E+06 +9574.48 4.69E+06 +10078.4 4.88E+06 +11086.24 5.25E+06 +12094.08 5.61E+06 +13101.92 5.94E+06 +14109.76 6.26E+06 +15117.6 6.57E+06 +16125.44 6.86E+06 +17133.28 7.14E+06 +18141.12 7.41E+06 +19148.96 7.67E+06 +20156.8 7.93E+06 +21164.64 8.17E+06 +22172.48 8.41E+06 +23180.32 8.63E+06 +24188.16 8.86E+06 +25196 9.07E+06 +26203.84 9.28E+06 +27211.68 9.48E+06 +28219.52 9.68E+06 +29227.36 9.87E+06 +30235.2 1.01E+07 +32250.88 1.04E+07 +34266.56 1.08E+07 +36282.24 1.11E+07 +38297.92 1.14E+07 +40313.6 1.17E+07 +42329.28 1.20E+07 +44344.96 1.23E+07 +46360.64 1.25E+07 +48376.32 1.28E+07 +50392 1.30E+07 diff --git a/MakeVertex.C b/MakeVertex.C index 85b5aa8..cb0cb5e 100755 --- a/MakeVertex.C +++ b/MakeVertex.C @@ -45,6 +45,9 @@ public: const int MAX_QQQ = 4; const int MAX_RING = 16; const int MAX_WEDGE = 16; +const double qqqpos=100.0; +const double vertexpos=14.2; +const double pcrad=37.0; double qqqGain[MAX_QQQ][MAX_RING][MAX_WEDGE] = {{{0}}}; bool qqqGainValid[MAX_QQQ][MAX_RING][MAX_WEDGE] = {{{false}}}; double qqqCalib[MAX_QQQ][MAX_RING][MAX_WEDGE] = {{{0}}}; @@ -248,53 +251,53 @@ Bool_t MakeVertex::Process(Long64_t entry) pc.CalIndex(); std::vector sx3Events; - if(sx3.multi>1) { - std::array Fsx3; - //std::cout << "-----" << std::endl; - for(int i=0; i=12) continue; - int id = sx3.id[i]; - if(sx3.ch[i]>=8) { - int sx3ch=sx3.ch[i]-8; - sx3ch=(sx3ch+3)%4; - if(sx3ch==0 || sx3ch==3) continue; - float value=sx3.e[i]; - int gch = sx3.id[i]*4+(sx3.ch[i]-8); - Fsx3.at(id).fillevent("BACK",sx3ch,value); - Fsx3.at(id).ts = static_cast(sx3.t[i]); - plotter->Fill2D("sx3backs_raw",100,0,100,800,0,4096,gch,sx3.e[i]); - } else { - int sx3ch=sx3.ch[i]/2; - double value=sx3.e[i]; - if(sx3.ch[i]%2==0) { - Fsx3.at(id).fillevent("FRONT_L",sx3ch,value*sx3RightGain[id][sx3ch]); - } else { - Fsx3.at(id).fillevent("FRONT_R",sx3ch,value); - } - } - } - for(int id=0; id<12; id++) { - Fsx3.at(id).validate(); - auto det = Fsx3.at(id); - bool no_charge_sharing_strict = det.valid_front_chans.size()==1 && det.valid_back_chans.size()==1; - if(det.valid) { - //std::cout << det.frontEL << " " << det.frontEL*sx3RightGain[id][det.stripF] << std::endl; - plotter->Fill2D("be_vs_x_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF)+"_b"+std::to_string(det.stripB),200,-1,1,800,0,8192, - det.frontX,det.backE,"evsx"); - //std::cout << sx3BackGain[id][det.stripF][det.stripB] << " " << sx3FrontGain[id][det.stripF] << std::endl; - plotter->Fill2D("matched_be_vs_x_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF),200,-30,30,800,0,8192, - det.frontX*sx3FrontGain[id][det.stripF]+sx3FrontOffset[id][det.stripF],det.backE*sx3BackGain[id][det.stripF][det.stripB],"evsx_matched"); - //plotter->Fill2D("fe_vs_x_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF)+"_"+std::to_string(det.stripB),200,-1,1,800,0,4096,det.frontX,det.backE,"evsx"); - plotter->Fill2D("l_vs_r_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF),800,0,4096,800,0,4096,det.frontEL,det.frontER,"l_vs_r"); - } - if(det.valid && (id ==9 || id==7 || id == 1 || id==3) && det.stripF!=DEFAULT_NULL && det.stripB!=DEFAULT_NULL) { - double z = det.frontX*sx3FrontGain[id][det.stripF]+sx3FrontOffset[id][det.stripF]; - double backE = det.backE*sx3BackGain[id][det.stripF][det.stripB]; - Event sx3ev(TVector3(0,0,z),backE,-1,det.ts,-1,det.stripB+4*id,det.stripF+4*id); - sx3Events.push_back(sx3ev); - } - } - } + // if(sx3.multi>1) { + // std::array Fsx3; + // //std::cout << "-----" << std::endl; + // for(int i=0; i=12) continue; + // int id = sx3.id[i]; + // if(sx3.ch[i]>=8) { + // int sx3ch=sx3.ch[i]-8; + // sx3ch=(sx3ch+3)%4; + // if(sx3ch==0 || sx3ch==3) continue; + // float value=sx3.e[i]; + // int gch = sx3.id[i]*4+(sx3.ch[i]-8); + // Fsx3.at(id).fillevent("BACK",sx3ch,value); + // Fsx3.at(id).ts = static_cast(sx3.t[i]); + // plotter->Fill2D("sx3backs_raw",100,0,100,800,0,4096,gch,sx3.e[i]); + // } else { + // int sx3ch=sx3.ch[i]/2; + // double value=sx3.e[i]; + // if(sx3.ch[i]%2==0) { + // Fsx3.at(id).fillevent("FRONT_L",sx3ch,value*sx3RightGain[id][sx3ch]); + // } else { + // Fsx3.at(id).fillevent("FRONT_R",sx3ch,value); + // } + // } + // } + // for(int id=0; id<12; id++) { + // Fsx3.at(id).validate(); + // auto det = Fsx3.at(id); + // bool no_charge_sharing_strict = det.valid_front_chans.size()==1 && det.valid_back_chans.size()==1; + // if(det.valid) { + // //std::cout << det.frontEL << " " << det.frontEL*sx3RightGain[id][det.stripF] << std::endl; + // plotter->Fill2D("be_vs_x_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF)+"_b"+std::to_string(det.stripB),200,-1,1,800,0,8192, + // det.frontX,det.backE,"evsx"); + // //std::cout << sx3BackGain[id][det.stripF][det.stripB] << " " << sx3FrontGain[id][det.stripF] << std::endl; + // plotter->Fill2D("matched_be_vs_x_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF)+"_"+std::to_string(id*4+det.stripF),200,-30,30,800,0,8192, + // det.frontX*sx3FrontGain[id][det.stripF]+sx3FrontOffset[id][det.stripF],det.backE*sx3BackGain[id][det.stripF][det.stripB],"evsx_matched"); + // //plotter->Fill2D("fe_vs_x_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF)+"_"+std::to_string(det.stripB),200,-1,1,800,0,4096,det.frontX,det.backE,"evsx"); + // plotter->Fill2D("l_vs_r_sx3_id_"+std::to_string(id)+"_f"+std::to_string(det.stripF),800,0,4096,800,0,4096,det.frontEL,det.frontER,"l_vs_r"); + // } + // if(det.valid && (id ==9 || id==7 || id == 1 || id==3) && det.stripF!=DEFAULT_NULL && det.stripB!=DEFAULT_NULL) { + // double z = det.frontX*sx3FrontGain[id][det.stripF]+sx3FrontOffset[id][det.stripF]; + // double backE = det.backE*sx3BackGain[id][det.stripF][det.stripB]; + // Event sx3ev(TVector3(0,0,z),backE,-1,det.ts,-1,det.stripB+4*id,det.stripF+4*id); + // sx3Events.push_back(sx3ev); + // } + // } + // } //return kTRUE; // QQQ Processing @@ -403,8 +406,8 @@ Bool_t MakeVertex::Process(Long64_t entry) double rho = 50. + (50. / 16.) * (chRing + 0.5); //"?" //z used to be 75+30+23=128 //we found a 12mm shift towards the vertex later --> 116 - Event qqqevent(TVector3(rho*TMath::Cos(theta),rho*TMath::Sin(theta),116), eRingMeV, eWedgeMeV, tRing, tWedge,chRing+qqq.id[i]*16, chWedge+qqq.id[i]*16); - Event qqqeventr(TVector3(rho*TMath::Cos(theta),rho*TMath::Sin(theta),116), eRing, eWedge, tRing, tWedge,chRing+qqq.id[i]*16, chWedge+qqq.id[i]*16); + Event qqqevent(TVector3(rho*TMath::Cos(theta),rho*TMath::Sin(theta),qqqpos), eRingMeV, eWedgeMeV, tRing, tWedge,chRing+qqq.id[i]*16, chWedge+qqq.id[i]*16); + Event qqqeventr(TVector3(rho*TMath::Cos(theta),rho*TMath::Sin(theta),qqqpos), eRing, eWedge, tRing, tWedge,chRing+qqq.id[i]*16, chWedge+qqq.id[i]*16); QQQ_Events.push_back(qqqevent); QQQ_Events_Raw.push_back(qqqeventr); plotter->Fill2D("QQQCartesianPlot", 200, -100, 100, 200, -100, 100, rho * TMath::Cos(theta), rho * TMath::Sin(theta), "hCalQQQ"); @@ -452,7 +455,7 @@ Bool_t MakeVertex::Process(Long64_t entry) double rho = 50. + (50. / 16.) * (chRing + 0.5); //"?" double x = rho * TMath::Cos(theta); double y = rho * TMath::Sin(theta); - hitPos.SetXYZ(x, y, (23 + 75 + 30)); + hitPos.SetXYZ(x, y, (qqqpos)); qqqenergy = eRingMeV; qqqtimestamp = tRing; HitNonZero = true; @@ -597,8 +600,8 @@ Bool_t MakeVertex::Process(Long64_t entry) plotter->Fill2D("dE_E_Cathodesx3B",400,0,10,800,0,10000,sx3event.Energy1*0.001,pcevent.Energy2); double sx3z = sx3event.pos.Z()+(75.0/2.0)-3.0; //w.r.t target origin at 90 for run12 double sx3rho = 88.0;//approximate barrel radius - double sx3theta = TMath::ATan2(sx3rho,sx3z-90); - double pczguess = 40.0/TMath::Tan(sx3theta) + 90.0; + double sx3theta = TMath::ATan2(sx3rho,sx3z-vertexpos); + double pczguess = pcrad/TMath::Tan(sx3theta) + vertexpos; plotter->Fill2D("pcz_vs_sx3pczguess",300,0,200,150,0,200,pczguess,pcevent.pos.Z()); plotter->Fill2D("pcz_vs_sx3pczguess"+std::to_string(sx3event.ch2),300,0,200,150,0,200,pczguess,pcevent.pos.Z()); plotter->Fill2D("pcz_vs_sx3z",300,0,200,150,0,200,sx3z,pcevent.pos.Z()); @@ -611,13 +614,13 @@ Bool_t MakeVertex::Process(Long64_t entry) plotter->Fill1D("dt_pcC_qqqW",640,-2000,2000,qqqevent.Time2 - pcevent.Time2); plotter->Fill2D("dE_E_AnodeQQQR",400,0,10,800,0,40000,qqqevent.Energy1,pcevent.Energy1); plotter->Fill2D("dE_E_CathodeQQQR",400,0,10,800,0,10000,qqqevent.Energy2,pcevent.Energy2); - double sinTheta = TMath::Sin((qqqevent.pos - TVector3(0,0,90)).Theta())/TMath::Sin((TVector3(51.5,0,128.) - TVector3(0,0,90)).Theta()); + double sinTheta = TMath::Sin((qqqevent.pos - TVector3(0,0,vertexpos)).Theta())/TMath::Sin((TVector3(51.5,0,qqqpos) - TVector3(0,0,vertexpos)).Theta()); plotter->Fill2D("dE2_E_AnodeQQQR",400,0,10,800,0,40000,qqqevent.Energy1,pcevent.Energy1*sinTheta); plotter->Fill2D("dE2_E_CathodeQQQR",400,0,10,800,0,10000,qqqevent.Energy2,pcevent.Energy2*sinTheta); if(qqqevent.pos.Phi() <= pcevent.pos.Phi()+TMath::Pi()/4. && qqqevent.pos.Phi() >= pcevent.pos.Phi()-TMath::Pi()/4.) { plotter->Fill1D("PCZ",800,-200,200,pcevent.pos.Z(),"phicut"); - double pcz_guess = 40.0/TMath::Tan((qqqevent.pos-TVector3(0,0,90)).Theta()) + 90; //this is ideally kept to be all QQQ+userinput for calibration of pcz + double pcz_guess = pcrad/TMath::Tan((qqqevent.pos-TVector3(0,0,vertexpos)).Theta()) + vertexpos; //this is ideally kept to be all QQQ+userinput for calibration of pcz plotter->Fill2D("pczguess_vs_pc",300,0,200,150,0,200,pcz_guess,pcevent.pos.Z(),"phicut"); plotter->Fill2D("pczguess_vs_pc_phi="+std::to_string(qqqevent.pos.Phi()*180./M_PI),300,0,200,150,0,200,pcz_guess,pcevent.pos.Z(),"phicut"); //plotter->Fill1D("PCZ",800,-200,200,pcevent.pos.Z(),"phicut"); diff --git a/MakeVertex.C:Zone.Identifier b/MakeVertex.C:Zone.Identifier deleted file mode 100644 index d6c1ec6..0000000 Binary files a/MakeVertex.C:Zone.Identifier and /dev/null differ diff --git a/MakeVertex.h:Zone.Identifier b/MakeVertex.h:Zone.Identifier deleted file mode 100644 index d6c1ec6..0000000 Binary files a/MakeVertex.h:Zone.Identifier and /dev/null differ