get_xsec_stop.py

Go to the documentation of this file.

#!/usr/bin/env python
import os, sys, re
import glob
import string
from array import array 
from pprint import pprint
import ROOT

def stopCrossSection( stop_mass):
    if(stop_mass == 100): xsec = 1521.11
    elif(stop_mass == 105): xsec = 1233.18
    elif(stop_mass == 110): xsec = 1013.76
    elif(stop_mass == 115): xsec = 832.656
    elif(stop_mass == 120): xsec = 689.799
    elif(stop_mass == 125): xsec = 574.981
    elif(stop_mass == 130): xsec = 481.397
    elif(stop_mass == 135): xsec = 405.159
    elif(stop_mass == 140): xsec = 342.865
    elif(stop_mass == 145): xsec = 291.752
    elif(stop_mass == 150): xsec = 249.409
    elif(stop_mass == 155): xsec = 214.221
    elif(stop_mass == 160): xsec = 184.623
    elif(stop_mass == 165): xsec = 159.614
    elif(stop_mass == 170): xsec = 139.252
    elif(stop_mass == 175): xsec = 121.416
    elif(stop_mass == 180): xsec = 106.194
    elif(stop_mass == 185): xsec = 93.3347
    elif(stop_mass == 190): xsec = 82.2541
    elif(stop_mass == 195): xsec = 72.7397
    elif(stop_mass == 200): xsec = 64.5085
    elif(stop_mass == 205): xsec = 57.2279
    elif(stop_mass == 210): xsec = 50.9226
    elif(stop_mass == 215): xsec = 45.3761
    elif(stop_mass == 220): xsec = 40.5941
    elif(stop_mass == 225): xsec = 36.3818
    elif(stop_mass == 230): xsec = 32.6679
    elif(stop_mass == 235): xsec = 29.3155
    elif(stop_mass == 240): xsec = 26.4761
    elif(stop_mass == 245): xsec = 23.8853
    elif(stop_mass == 250): xsec = 21.5949
    elif(stop_mass == 255): xsec = 19.5614
    elif(stop_mass == 260): xsec = 17.6836
    elif(stop_mass == 265): xsec = 16.112 
    elif(stop_mass == 270): xsec = 14.6459
    elif(stop_mass == 275): xsec = 13.3231
    elif(stop_mass == 280): xsec = 12.1575
    elif(stop_mass == 285): xsec = 11.0925
    elif(stop_mass == 290): xsec = 10.1363
    elif(stop_mass == 295): xsec = 9.29002
    elif(stop_mass == 300): xsec = 8.51615
    elif(stop_mass == 305): xsec = 7.81428
    elif(stop_mass == 310): xsec = 7.17876
    elif(stop_mass == 315): xsec = 6.60266
    elif(stop_mass == 320): xsec = 6.08444
    elif(stop_mass == 325): xsec = 5.60471
    elif(stop_mass == 330): xsec = 5.17188
    elif(stop_mass == 335): xsec = 4.77871
    elif(stop_mass == 340): xsec = 4.41629
    elif(stop_mass == 345): xsec = 4.08881
    elif(stop_mass == 350): xsec = 3.78661
    elif(stop_mass == 355): xsec = 3.50911
    elif(stop_mass == 360): xsec = 3.25619
    elif(stop_mass == 365): xsec = 3.02472
    elif(stop_mass == 370): xsec = 2.8077 
    elif(stop_mass == 375): xsec = 2.61162
    elif(stop_mass == 380): xsec = 2.43031
    elif(stop_mass == 385): xsec = 2.26365
    elif(stop_mass == 390): xsec = 2.10786
    elif(stop_mass == 395): xsec = 1.9665 
    elif(stop_mass == 400): xsec = 1.83537
    elif(stop_mass == 405): xsec = 1.70927
    elif(stop_mass == 410): xsec = 1.60378
    elif(stop_mass == 415): xsec = 1.49798
    elif(stop_mass == 420): xsec = 1.39688
    elif(stop_mass == 425): xsec = 1.31169
    elif(stop_mass == 430): xsec = 1.22589
    elif(stop_mass == 435): xsec = 1.14553
    elif(stop_mass == 440): xsec = 1.07484
    elif(stop_mass == 445): xsec = 1.01019
    elif(stop_mass == 450): xsec = 0.948333
    elif(stop_mass == 455): xsec = 0.890847
    elif(stop_mass == 460): xsec = 0.836762
    elif(stop_mass == 465): xsec = 0.787221
    elif(stop_mass == 470): xsec = 0.740549
    elif(stop_mass == 475): xsec = 0.697075
    elif(stop_mass == 480): xsec = 0.655954
    elif(stop_mass == 485): xsec = 0.618562
    elif(stop_mass == 490): xsec = 0.582467
    elif(stop_mass == 495): xsec = 0.549524
    elif(stop_mass == 500): xsec = 0.51848 
    elif(stop_mass == 505): xsec = 0.489324
    elif(stop_mass == 510): xsec = 0.462439
    elif(stop_mass == 515): xsec = 0.436832
    elif(stop_mass == 520): xsec = 0.412828
    elif(stop_mass == 525): xsec = 0.390303
    elif(stop_mass == 530): xsec = 0.368755
    elif(stop_mass == 535): xsec = 0.348705
    elif(stop_mass == 540): xsec = 0.330157
    elif(stop_mass == 545): xsec = 0.312672
    elif(stop_mass == 550): xsec = 0.296128
    elif(stop_mass == 555): xsec = 0.280734
    elif(stop_mass == 560): xsec = 0.266138
    elif(stop_mass == 565): xsec = 0.251557
    elif(stop_mass == 570): xsec = 0.238537
    elif(stop_mass == 575): xsec = 0.226118
    elif(stop_mass == 580): xsec = 0.214557
    elif(stop_mass == 585): xsec = 0.203566
    elif(stop_mass == 590): xsec = 0.193079
    elif(stop_mass == 595): xsec = 0.183604
    elif(stop_mass == 600): xsec = 0.174599
    elif(stop_mass == 605): xsec = 0.166131
    elif(stop_mass == 610): xsec = 0.158242
    elif(stop_mass == 615): xsec = 0.150275
    elif(stop_mass == 620): xsec = 0.142787
    elif(stop_mass == 625): xsec = 0.136372
    elif(stop_mass == 630): xsec = 0.129886
    elif(stop_mass == 635): xsec = 0.123402
    elif(stop_mass == 640): xsec = 0.11795 
    elif(stop_mass == 645): xsec = 0.112008
    elif(stop_mass == 650): xsec = 0.107045
    elif(stop_mass == 655): xsec = 0.102081
    elif(stop_mass == 660): xsec = 0.09725 
    elif(stop_mass == 665): xsec = 0.0927515
    elif(stop_mass == 670): xsec = 0.0885084
    elif(stop_mass == 675): xsec = 0.0844877
    elif(stop_mass == 680): xsec = 0.0806192
    elif(stop_mass == 685): xsec = 0.0769099
    elif(stop_mass == 690): xsec = 0.0734901
    elif(stop_mass == 695): xsec = 0.0701805
    elif(stop_mass == 700): xsec = 0.0670476
    elif(stop_mass == 705): xsec = 0.0641426
    elif(stop_mass == 710): xsec = 0.0612942
    elif(stop_mass == 715): xsec = 0.0585678
    elif(stop_mass == 720): xsec = 0.0560753
    elif(stop_mass == 725): xsec = 0.0536438
    elif(stop_mass == 730): xsec = 0.0513219
    elif(stop_mass == 735): xsec = 0.0491001
    elif(stop_mass == 740): xsec = 0.0470801
    elif(stop_mass == 745): xsec = 0.045061 
    elif(stop_mass == 750): xsec = 0.0431418
    elif(stop_mass == 755): xsec = 0.0413447
    elif(stop_mass == 760): xsec = 0.0396264
    elif(stop_mass == 765): xsec = 0.0379036
    elif(stop_mass == 770): xsec = 0.0363856
    elif(stop_mass == 775): xsec = 0.0348796
    elif(stop_mass == 780): xsec = 0.0334669
    elif(stop_mass == 785): xsec = 0.0320548
    elif(stop_mass == 790): xsec = 0.0307373
    elif(stop_mass == 795): xsec = 0.0295348
    elif(stop_mass == 800): xsec = 0.0283338
    elif(stop_mass == 805): xsec = 0.0272206
    elif(stop_mass == 810): xsec = 0.0261233
    elif(stop_mass == 815): xsec = 0.0251107
    elif(stop_mass == 820): xsec = 0.0241099
    elif(stop_mass == 825): xsec = 0.0230866
    elif(stop_mass == 830): xsec = 0.0221834
    elif(stop_mass == 835): xsec = 0.0213766
    elif(stop_mass == 840): xsec = 0.0204715
    elif(stop_mass == 845): xsec = 0.0197653
    elif(stop_mass == 850): xsec = 0.0189612
    elif(stop_mass == 855): xsec = 0.0182516
    elif(stop_mass == 860): xsec = 0.0175509
    elif(stop_mass == 865): xsec = 0.0168336
    elif(stop_mass == 870): xsec = 0.0162314
    elif(stop_mass == 875): xsec = 0.015625 
    elif(stop_mass == 880): xsec = 0.0150143
    elif(stop_mass == 885): xsec = 0.0144112
    elif(stop_mass == 890): xsec = 0.0138979
    elif(stop_mass == 895): xsec = 0.0133962
    elif(stop_mass == 900): xsec = 0.0128895
    elif(stop_mass == 905): xsec = 0.0123843
    elif(stop_mass == 910): xsec = 0.0119837
    elif(stop_mass == 915): xsec = 0.0114713
    elif(stop_mass == 920): xsec = 0.0110688
    elif(stop_mass == 925): xsec = 0.0106631
    elif(stop_mass == 930): xsec = 0.0102629
    elif(stop_mass == 935): xsec = 0.0098874
    elif(stop_mass == 940): xsec = 0.00952142
    elif(stop_mass == 945): xsec = 0.00916636
    elif(stop_mass == 950): xsec = 0.00883465
    elif(stop_mass == 955): xsec = 0.00851073
    elif(stop_mass == 960): xsec = 0.00820884
    elif(stop_mass == 965): xsec = 0.00791403
    elif(stop_mass == 970): xsec = 0.00763112
    elif(stop_mass == 975): xsec = 0.00735655
    elif(stop_mass == 980): xsec = 0.00710317
    elif(stop_mass == 985): xsec = 0.00684867
    elif(stop_mass == 990): xsec = 0.00660695
    elif(stop_mass == 995): xsec = 0.00637546
    elif(stop_mass == 1000): xsec = 0.00615134
    elif(stop_mass == 1005): xsec = 0.00593765
    elif(stop_mass == 1010): xsec = 0.00572452
    elif(stop_mass == 1015): xsec = 0.00553094
    elif(stop_mass == 1020): xsec = 0.00533968
    elif(stop_mass == 1025): xsec = 0.00514619
    elif(stop_mass == 1030): xsec = 0.00497235
    elif(stop_mass == 1035): xsec = 0.00479906
    elif(stop_mass == 1040): xsec = 0.00463806
    elif(stop_mass == 1045): xsec = 0.00447537
    elif(stop_mass == 1050): xsec = 0.00432261
    elif(stop_mass == 1055): xsec = 0.00417983
    elif(stop_mass == 1060): xsec = 0.00403886
    elif(stop_mass == 1065): xsec = 0.0038962
    elif(stop_mass == 1070): xsec = 0.00376343
    elif(stop_mass == 1075): xsec = 0.00364174
    elif(stop_mass == 1080): xsec = 0.00352093
    elif(stop_mass == 1085): xsec = 0.00339813
    elif(stop_mass == 1090): xsec = 0.00328695
    elif(stop_mass == 1095): xsec = 0.00317628
    elif(stop_mass == 1100): xsec = 0.00307413
    elif(stop_mass == 1105): xsec = 0.00297377
    elif(stop_mass == 1110): xsec = 0.00287148
    elif(stop_mass == 1115): xsec = 0.00278078
    elif(stop_mass == 1120): xsec = 0.00268873
    elif(stop_mass == 1125): xsec = 0.00260821
    elif(stop_mass == 1130): xsec = 0.00251529
    elif(stop_mass == 1135): xsec = 0.00243484
    elif(stop_mass == 1140): xsec = 0.00236295
    elif(stop_mass == 1145): xsec = 0.00228192
    elif(stop_mass == 1150): xsec = 0.00221047
    elif(stop_mass == 1155): xsec = 0.00213907
    elif(stop_mass == 1160): xsec = 0.00206845
    elif(stop_mass == 1165): xsec = 0.0020063
    elif(stop_mass == 1170): xsec = 0.00194569
    elif(stop_mass == 1175): xsec = 0.0018741
    elif(stop_mass == 1180): xsec = 0.00182266
    elif(stop_mass == 1185): xsec = 0.00176211
    elif(stop_mass == 1190): xsec = 0.00170006
    elif(stop_mass == 1195): xsec = 0.00164968
    elif(stop_mass == 1200): xsec = 0.00159844
    elif(stop_mass == 1205): xsec = 0.0015472
    elif(stop_mass == 1210): xsec = 0.00149657
    elif(stop_mass == 1215): xsec = 0.00145544
    elif(stop_mass == 1220): xsec = 0.00140288
    elif(stop_mass == 1225): xsec = 0.00136155
    elif(stop_mass == 1230): xsec = 0.00131271
    elif(stop_mass == 1235): xsec = 0.0012717
    elif(stop_mass == 1240): xsec = 0.00123066
    elif(stop_mass == 1245): xsec = 0.00119994
    elif(stop_mass == 1250): xsec = 0.0011583
    elif(stop_mass == 1255): xsec = 0.00112694
    elif(stop_mass == 1260): xsec = 0.00108716
    elif(stop_mass == 1265): xsec = 0.00105517
    elif(stop_mass == 1270): xsec = 0.00102241
    elif(stop_mass == 1275): xsec = 0.00099129
    elif(stop_mass == 1280): xsec = 0.00096101
    elif(stop_mass == 1285): xsec = 0.00093239
    elif(stop_mass == 1290): xsec = 0.00090340
    elif(stop_mass == 1295): xsec = 0.00087695
    elif(stop_mass == 1300): xsec = 0.00085034
    elif(stop_mass == 1305): xsec = 0.00082443
    elif(stop_mass == 1310): xsec = 0.00079983
    elif(stop_mass == 1315): xsec = 0.00077522
    elif(stop_mass == 1320): xsec = 0.00075137
    elif(stop_mass == 1325): xsec = 0.00072891
    elif(stop_mass == 1330): xsec = 0.00070686
    elif(stop_mass == 1335): xsec = 0.00068537
    elif(stop_mass == 1340): xsec = 0.00066464
    elif(stop_mass == 1345): xsec = 0.00064480
    elif(stop_mass == 1350): xsec = 0.00062515
    elif(stop_mass == 1355): xsec = 0.00060680
    elif(stop_mass == 1360): xsec = 0.00058851
    elif(stop_mass == 1365): xsec = 0.00057050
    elif(stop_mass == 1370): xsec = 0.00055337
    elif(stop_mass == 1375): xsec = 0.00053664
    elif(stop_mass == 1380): xsec = 0.00052140
    elif(stop_mass == 1385): xsec = 0.00050500
    elif(stop_mass == 1390): xsec = 0.00049035
    elif(stop_mass == 1395): xsec = 0.00047616
    elif(stop_mass == 1400): xsec = 0.00046194
    elif(stop_mass == 1405): xsec = 0.00044817
    elif(stop_mass == 1410): xsec = 0.00043508
    elif(stop_mass == 1415): xsec = 0.00042296
    elif(stop_mass == 1420): xsec = 0.00041038
    elif(stop_mass == 1425): xsec = 0.00039810
    elif(stop_mass == 1430): xsec = 0.00038679
    elif(stop_mass == 1435): xsec = 0.00037572
    elif(stop_mass == 1440): xsec = 0.00036461
    elif(stop_mass == 1445): xsec = 0.00035396
    elif(stop_mass == 1450): xsec = 0.00034392
    elif(stop_mass == 1455): xsec = 0.00033388
    elif(stop_mass == 1460): xsec = 0.00032434
    elif(stop_mass == 1465): xsec = 0.0003153
    elif(stop_mass == 1470): xsec = 0.00030583
    elif(stop_mass == 1475): xsec = 0.00029681
    elif(stop_mass == 1480): xsec = 0.00028814
    elif(stop_mass == 1485): xsec = 0.00027971
    elif(stop_mass == 1490): xsec = 0.00027172
    elif(stop_mass == 1495): xsec = 0.00026427
    elif(stop_mass == 1500): xsec = 0.00025624
    elif(stop_mass == 1505): xsec = 0.00024885
    elif(stop_mass == 1510): xsec = 0.00024184
    elif(stop_mass == 1515): xsec = 0.00023443
    elif(stop_mass == 1520): xsec = 0.00022737
    elif(stop_mass == 1525): xsec = 0.00022104
    elif(stop_mass == 1530): xsec = 0.00021443
    elif(stop_mass == 1535): xsec = 0.00020809
    elif(stop_mass == 1540): xsec = 0.00020174
    elif(stop_mass == 1545): xsec = 0.00019639
    elif(stop_mass == 1550): xsec = 0.00019047
    elif(stop_mass == 1555): xsec = 0.00018518
    elif(stop_mass == 1560): xsec = 0.00017926
    elif(stop_mass == 1565): xsec = 0.00017402
    elif(stop_mass == 1570): xsec = 0.00016917
    elif(stop_mass == 1575): xsec = 0.00016386
    elif(stop_mass == 1580): xsec = 0.00015958
    elif(stop_mass == 1585): xsec = 0.00015471
    elif(stop_mass == 1590): xsec = 0.00015050
    elif(stop_mass == 1595): xsec = 0.00014562
    elif(stop_mass == 1600): xsec = 0.00014138
    elif(stop_mass == 1605): xsec = 0.00013713
    elif(stop_mass == 1610): xsec = 0.00013218
    elif(stop_mass == 1615): xsec = 0.00012792
    elif(stop_mass == 1620): xsec = 0.00012408
    elif(stop_mass == 1625): xsec = 0.00011982
    elif(stop_mass == 1630): xsec = 0.00011604
    elif(stop_mass == 1635): xsec = 0.00011276
    elif(stop_mass == 1640): xsec = 0.00010893
    elif(stop_mass == 1645): xsec = 0.00010574
    elif(stop_mass == 1650): xsec = 0.00010269
    elif(stop_mass == 1655): xsec = 0.00010011
    elif(stop_mass == 1660): xsec = 9.75763e-05
    elif(stop_mass == 1665): xsec = 9.52062e-05
    elif(stop_mass == 1670): xsec = 9.29857e-05
    elif(stop_mass == 1675): xsec = 9.08285e-05
    elif(stop_mass == 1680): xsec = 8.87433e-05
    elif(stop_mass == 1685): xsec = 8.66618e-05
    elif(stop_mass == 1690): xsec = 8.46535e-05
    elif(stop_mass == 1695): xsec = 8.27102e-05
    elif(stop_mass == 1700): xsec = 8.07774e-05
    elif(stop_mass == 1705): xsec = 7.8666e-05
    elif(stop_mass == 1710): xsec = 7.6572e-05
    elif(stop_mass == 1715): xsec = 7.45994e-05
    elif(stop_mass == 1720): xsec = 7.25199e-05
    elif(stop_mass == 1725): xsec = 7.05189e-05
    elif(stop_mass == 1730): xsec = 6.85712e-05
    elif(stop_mass == 1735): xsec = 6.67296e-05
    elif(stop_mass == 1740): xsec = 6.49184e-05
    elif(stop_mass == 1745): xsec = 6.30949e-05
    elif(stop_mass == 1750): xsec = 6.13637e-05
    elif(stop_mass == 1755): xsec = 5.97301e-05
    elif(stop_mass == 1760): xsec = 5.80751e-05
    elif(stop_mass == 1765): xsec = 5.65479e-05
    elif(stop_mass == 1770): xsec = 5.49998e-05
    elif(stop_mass == 1775): xsec = 5.35686e-05
    elif(stop_mass == 1780): xsec = 5.20828e-05
    elif(stop_mass == 1785): xsec = 5.07079e-05
    elif(stop_mass == 1790): xsec = 4.93948e-05
    elif(stop_mass == 1795): xsec = 4.80635e-05
    elif(stop_mass == 1800): xsec = 4.67492e-05
    elif(stop_mass == 1805): xsec = 4.55055e-05
    elif(stop_mass == 1810): xsec = 4.42835e-05
    elif(stop_mass == 1815): xsec = 4.30744e-05
    elif(stop_mass == 1820): xsec = 4.19954e-05
    elif(stop_mass == 1825): xsec = 4.08527e-05
    elif(stop_mass == 1830): xsec = 3.97561e-05
    elif(stop_mass == 1835): xsec = 3.87041e-05
    elif(stop_mass == 1840): xsec = 3.76008e-05
    elif(stop_mass == 1845): xsec = 3.66914e-05
    elif(stop_mass == 1850): xsec = 3.56995e-05
    elif(stop_mass == 1855): xsec = 3.47689e-05
    elif(stop_mass == 1860): xsec = 3.38528e-05
    elif(stop_mass == 1865): xsec = 3.29644e-05
    elif(stop_mass == 1870): xsec = 3.20679e-05
    elif(stop_mass == 1875): xsec = 3.12583e-05
    elif(stop_mass == 1880): xsec = 3.04342e-05
    elif(stop_mass == 1885): xsec = 2.96516e-05
    elif(stop_mass == 1890): xsec = 2.88952e-05
    elif(stop_mass == 1895): xsec = 2.81145e-05
    elif(stop_mass == 1900): xsec = 2.73974e-05
    elif(stop_mass == 1905): xsec = 2.66796e-05
    elif(stop_mass == 1910): xsec = 2.59941e-05
    elif(stop_mass == 1915): xsec = 2.52784e-05
    elif(stop_mass == 1920): xsec = 2.46598e-05
    elif(stop_mass == 1925): xsec = 2.39932e-05
    elif(stop_mass == 1930): xsec = 2.33737e-05
    elif(stop_mass == 1935): xsec = 2.27623e-05
    elif(stop_mass == 1940): xsec = 2.21454e-05
    elif(stop_mass == 1945): xsec = 2.15924e-05
    elif(stop_mass == 1950): xsec = 2.10232e-05
    elif(stop_mass == 1955): xsec = 2.05211e-05
    elif(stop_mass == 1960): xsec = 1.98996e-05
    elif(stop_mass == 1965): xsec = 1.9408e-05
    elif(stop_mass == 1970): xsec = 1.88974e-05
    elif(stop_mass == 1975): xsec = 1.84612e-05
    elif(stop_mass == 1980): xsec = 1.79562e-05
    elif(stop_mass == 1985): xsec = 1.75673e-05
    elif(stop_mass == 1990): xsec = 1.70612e-05
    elif(stop_mass == 1995): xsec = 1.66228e-05
    elif(stop_mass == 2000): xsec = 1.62355e-05
    else: xsec = 0.

    return xsec