#include "TROOT.h"

#include "TStyle.h"

#include "TCanvas.h"

#include "TH2F.h"

#include "TGraph.h"

#include "TPolyLine.h"

#include "TF1.h"

#include "TLegend.h"

#include "TPad.h"

#include "TLatex.h"

#include "TSpline.h"

#include "TFile.h"

using namespace std;

void myStyle();

void setTDRStyle();

void ComboPlotLQ1()

{

// **********************************************

// * Input parameters *

// **********************************************

static const int numbetas = 178;

Double_t beta_vals[178]={0.006,0.008,0.01,0.012,0.014,0.016,0.018,0.02,0.022,0.024,0.026,0.028,0.03,0.032,0.034,0.036,0.038,0.04,0.042,0.044,0.046,0.048,0.05,0.052,0.054,0.056,0.058,0.06,0.062,0.064,0.066,0.068,0.07,0.072,0.074,0.076,0.078,0.08,0.082,0.084,0.086,0.088,0.09,0.092,0.094,0.096,0.098,0.1,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.19,0.2,0.21,0.22,0.23,0.24,0.25,0.26,0.27,0.28,0.29,0.3,0.31,0.32,0.33,0.34,0.35,0.36,0.37,0.38,0.39,0.4,0.41,0.42,0.43,0.44,0.45,0.46,0.47,0.48,0.49,0.5,0.51,0.52,0.53,0.54,0.55,0.56,0.57,0.58,0.59,0.6,0.61,0.62,0.63,0.64,0.65,0.66,0.67,0.68,0.69,0.7,0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.78,0.79,0.8,0.81,0.82,0.83,0.84,0.85,0.86,0.87,0.88,0.89,0.9,0.902,0.904,0.906,0.908,0.91,0.912,0.914,0.916,0.918,0.92,0.922,0.924,0.926,0.928,0.93,0.932,0.934,0.936,0.938,0.94,0.942,0.944,0.946,0.948,0.95,0.952,0.954,0.956,0.958,0.96,0.962,0.964,0.966,0.968,0.97,0.972,0.974,0.976,0.978,0.98,0.982,0.984,0.986,0.988,0.99,0.992,0.994,0.996,0.998,0.9995};

Double_t m_expected_combo[178]={300.0,300.0,300.0,300.0,300.0,300.0,300.0,362.1,406.1,437.9,467.2,495.2,512.8,529.0,544.0,556.7,567.6,577.0,587.1,596.5,606.4,615.5,624.1,633.7,641.9,650.0,656.8,663.1,669.8,675.7,681.5,686.9,692.4,698.0,702.4,705.8,709.2,712.7,716.1,719.0,722.4,725.3,728.3,731.3,734.0,737.0,739.9,742.4,754.5,764.3,773.1,781.8,789.3,796.9,803.9,811.4,818.5,825.0,831.4,837.8,843.8,849.0,854.2,859.4,863.8,868.4,872.9,877.2,881.5,885.1,888.6,892.3,895.9,898.9,901.9,904.3,906.9,909.5,911.7,913.9,916.0,918.6,920.5,922.7,925.0,927.0,929.3,931.6,933.3,935.6,937.8,939.8,942.3,944.4,946.6,948.4,951.1,952.5,954.8,956.6,958.9,960.8,962.7,965.1,966.8,968.9,970.9,972.7,974.8,976.2,978.1,980.0,982.1,983.5,985.5,987.6,989.3,991.3,993.3,995.1,996.9,999.1,1000.9,1002.4,1004.0,1005.6,1007.3,1009.2,1009.6,1009.9,1010.6,1010.7,1011.0,1011.7,1011.7,1012.0,1012.4,1012.9,1013.1,1013.8,1014.1,1014.6,1014.8,1015.2,1015.6,1015.8,1016.3,1016.8,1017.1,1017.5,1018.1,1018.5,1019.1,1019.2,1019.6,1020.2,1020.4,1020.9,1021.5,1022.0,1022.2,1022.5,1023.0,1023.5,1024.1,1024.7,1024.8,1025.3,1025.8,1026.0,1026.8,1027.3,1027.6,1028.1,1028.5,1029.3,1029.5,1030.0};

Double_t m_observed_combo[178]={300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,354.8,372.6,391.9,418.1,509.8,514.2,518.4,522.4,526.1,529.7,533.2,536.6,539.7,542.8,545.7,548.5,551.4,554.4,557.2,560.0,562.6,565.2,567.6,570.0,572.3,574.6,576.7,578.8,580.9,582.9,584.9,586.9,588.8,590.6,592.4,594.1,595.7,597.3,598.9,607.3,615.0,622.1,628.6,634.4,639.8,644.7,649.2,769.7,774.0,778.0,781.6,785.1,788.4,791.7,794.9,797.9,800.8,803.5,806.2,808.6,810.9,813.0,815.1,817.0,819.0,820.9,822.8,824.5,826.2,827.9,829.6,831.2,833.1,834.7,836.3,837.9,839.5,841.0,842.5,843.9,845.3,846.7,848.1,849.5,860.3,875.3,887.1,900.8,904.9,909.1,913.2,917.3,922.1,926.8,931.8,937.1,942.7,948.7,951.2,952.8,954.4,956.3,957.8,959.2,960.7,962.5,963.8,965.4,967.1,968.7,970.1,971.7,973.5,975.1,976.9,978.6,980.5,982.3,984.1,984.5,984.9,985.3,985.6,986.0,986.4,986.8,987.3,987.7,988.1,988.2,988.7,989.1,989.5,989.9,990.3,990.7,991.2,991.6,992.0,992.5,992.9,993.4,993.8,994.3,994.7,995.2,995.7,996.1,996.6,997.1,997.5,998.0,998.5,999.0,999.5,1000.1,1000.9,1001.1,1001.8,1002.7,1003.3,1004.0,1004.7,1005.3,1005.9,1006.5,1007.2,1007.9,1008.4};

Double_t m_1sigma_combo[356]={300.0,300.0,300.0,300.0,354.0,411.8,457.4,494.2,518.9,540.8,559.6,573.9,588.3,601.1,614.8,627.8,640.5,651.3,660.8,670.6,678.3,687.3,694.9,701.7,706.9,712.1,716.5,721.4,725.8,730.0,734.4,738.5,742.4,746.4,750.1,753.3,756.3,759.3,762.2,765.0,767.8,770.3,773.2,775.8,778.4,781.0,783.2,785.5,796.8,808.5,819.4,829.9,839.8,848.8,857.2,865.0,872.3,879.8,885.6,891.7,897.8,902.4,906.2,910.1,913.5,916.9,920.0,923.1,926.4,929.5,932.3,935.3,938.3,940.9,943.6,946.2,949.0,952.2,955.1,957.8,960.7,963.5,966.2,968.9,971.7,974.1,976.7,979.3,981.5,983.8,986.1,988.4,990.8,993.1,995.2,997.2,999.7,1001.5,1002.9,1004.7,1006.8,1008.0,1009.8,1011.3,1013.0,1014.8,1016.6,1018.2,1020.0,1021.3,1022.8,1024.3,1025.9,1027.3,1028.9,1030.6,1031.9,1033.8,1035.1,1036.8,1038.2,1039.8,1041.5,1043.1,1044.7,1046.3,1048.0,1049.7,1050.4,1050.7,1051.0,1051.4,1051.7,1052.3,1052.3,1052.7,1053.0,1053.6,1053.7,1054.0,1054.3,1055.0,1055.0,1055.3,1055.7,1056.1,1056.4,1057.0,1057.3,1057.7,1058.1,1058.4,1058.7,1059.1,1059.4,1059.8,1060.2,1060.8,1061.1,1061.8,1061.9,1062.3,1062.7,1063.1,1063.7,1064.0,1064.4,1064.8,1065.4,1065.6,1066.2,1066.8,1067.1,1067.7,1067.7,1068.4,1068.8,1069.3,987.6,986.9,986.9,986.1,985.1,984.6,984.1,983.7,983.0,982.6,982.0,981.5,981.2,980.4,980.1,979.3,978.8,978.3,977.9,977.4,976.9,976.3,976.0,975.2,974.7,974.7,973.8,973.4,972.8,972.3,971.9,971.7,970.9,970.5,970.1,969.9,969.3,968.8,968.4,967.6,967.2,966.7,966.3,965.9,965.9,965.1,964.7,964.4,963.8,963.4,963.0,960.8,958.8,956.8,955.2,952.9,951.0,948.7,946.3,943.9,941.2,939.0,937.4,934.6,932.3,930.8,928.3,926.5,924.2,922.4,920.7,918.4,916.8,914.7,913.6,911.1,909.3,907.3,905.4,903.7,901.9,900.2,898.3,895.8,893.9,891.6,889.4,887.6,885.5,883.5,881.4,879.3,876.5,874.1,871.3,868.7,866.3,863.7,861.3,857.9,855.4,852.1,849.8,846.4,843.1,840.4,836.9,833.2,829.3,825.9,822.0,818.1,814.4,809.8,805.2,800.6,796.2,791.6,787.0,782.4,777.3,771.7,766.0,760.2,753.6,746.2,737.1,727.9,717.6,706.8,692.3,688.3,684.3,680.5,676.5,672.9,668.9,664.1,659.9,655.6,650.9,645.5,639.8,633.9,627.4,621.0,614.7,608.4,602.2,596.0,589.2,581.9,575.5,568.5,561.5,554.0,544.8,534.1,522.9,511.8,500.2,479.3,462.6,440.3,417.5,391.3,360.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0};

Double_t m_2sigma_combo[356]={300.0,300.0,300.0,412.7,469.9,512.5,542.8,564.9,583.8,602.5,620.6,637.4,653.2,665.9,677.8,689.1,700.0,706.8,713.8,720.2,726.0,732.4,737.9,743.8,748.7,753.9,757.6,761.7,765.7,769.3,773.2,776.8,780.3,783.7,786.8,790.2,793.4,796.6,799.4,803.0,806.5,809.6,812.9,815.9,819.2,822.4,825.0,827.9,842.1,855.1,866.3,877.4,886.7,895.7,902.9,908.6,913.8,918.5,922.9,927.5,931.6,935.7,939.7,943.6,947.4,951.2,955.4,959.6,963.6,967.6,971.0,974.9,978.6,982.1,985.4,988.6,992.1,995.0,998.3,1001.0,1003.2,1005.7,1008.0,1010.2,1012.6,1014.8,1017.0,1019.0,1021.2,1023.1,1025.2,1027.1,1029.0,1030.9,1032.8,1034.3,1036.5,1038.1,1039.7,1041.6,1043.5,1044.7,1046.5,1047.9,1050.0,1051.0,1052.5,1054.0,1055.7,1057.1,1058.8,1060.0,1061.5,1062.8,1064.3,1065.8,1067.0,1068.7,1070.0,1071.4,1073.0,1074.2,1075.9,1077.4,1078.8,1080.5,1082.0,1083.8,1084.2,1084.4,1084.8,1085.4,1085.7,1085.9,1086.4,1086.6,1087.0,1087.3,1087.6,1087.9,1088.4,1088.9,1089.0,1089.4,1090.1,1090.3,1090.7,1091.1,1091.3,1091.9,1092.2,1092.6,1092.7,1093.2,1093.5,1094.1,1094.4,1094.8,1095.1,1095.8,1096.0,1096.7,1097.0,1097.7,1098.0,1098.0,1098.7,1099.0,1099.7,1100.0,1100.3,1100.9,1101.6,1102.0,1102.0,1102.5,1103.0,1103.6,945.5,944.4,944.1,942.9,941.5,940.7,940.1,939.2,938.6,938.1,937.1,936.3,935.2,934.4,934.1,933.3,932.8,931.6,930.9,930.3,929.8,928.8,928.1,927.5,926.6,926.0,925.1,924.5,924.0,923.4,922.5,922.5,921.4,921.2,920.2,920.3,919.2,918.7,917.7,917.3,916.4,916.3,915.4,915.2,914.4,914.2,913.2,913.1,912.1,912.1,911.0,908.5,905.8,903.0,900.8,898.8,896.6,894.7,892.9,891.0,889.1,887.2,886.0,884.0,882.2,880.7,878.7,877.0,874.5,872.5,870.7,868.7,867.2,864.9,863.6,861.3,859.5,857.5,855.6,853.9,851.7,850.1,848.1,845.9,844.1,842.2,840.5,838.9,836.5,834.4,832.2,830.0,827.7,825.8,823.5,821.2,818.8,816.6,814.3,811.5,808.8,806.2,803.9,800.6,798.0,795.5,792.3,789.5,786.3,783.4,780.2,776.6,773.6,769.8,766.2,762.5,758.3,754.3,749.8,744.3,738.7,732.4,726.2,719.3,711.9,704.4,693.7,680.8,666.8,651.4,631.6,626.8,622.1,617.2,613.1,608.3,604.6,599.3,594.3,590.2,584.8,580.0,575.0,570.0,565.0,559.8,554.6,548.3,540.8,532.5,524.9,516.8,508.7,501.1,484.7,470.4,459.2,440.8,423.8,407.4,384.1,361.3,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0};

Double_t m_expected_lljj[178]={300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,301.6,306.2,311.1,315.4,338.5,361.9,385.8,407.9,426.6,445.3,462.0,476.9,492.1,508.0,524.6,541.1,555.3,567.1,578.3,589.9,600.5,612.8,624.7,637.1,648.0,659.6,671.8,682.7,694.0,704.3,713.7,723.1,732.2,741.7,750.3,758.6,767.1,775.2,783.2,791.4,800.0,805.2,810.5,815.7,821.0,826.1,831.0,836.4,841.6,846.2,851.5,857.6,862.9,868.0,873.5,878.9,884.0,888.8,893.4,899.0,904.7,911.6,918.0,924.1,929.8,935.7,942.0,947.1,951.9,955.8,959.6,963.5,967.0,970.6,974.1,977.8,981.3,985.1,988.0,991.6,994.7,998.0,1001.1,1004.0,1004.7,1004.9,1005.5,1006.3,1006.8,1007.5,1007.9,1008.5,1008.9,1009.2,1009.9,1010.6,1011.0,1011.7,1012.0,1012.7,1013.0,1013.8,1014.1,1014.8,1015.2,1015.9,1016.3,1017.0,1017.4,1018.0,1018.5,1019.0,1019.8,1020.0,1020.6,1021.1,1021.5,1022.0,1022.6,1023.0,1023.7,1024.3,1024.5,1025.3,1025.6,1026.2,1026.4,1027.3,1027.6,1028.1,1028.5,1029.1,1029.5,1030.0};

Double_t m_observed_lljj[178]={300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.6,304.3,307.9,311.4,314.8,332.0,348.7,365.7,382.5,399.0,410.4,420.9,431.0,440.6,450.1,458.0,465.2,472.1,478.8,485.6,492.2,498.5,524.9,551.8,556.8,561.8,566.6,571.2,575.7,580.1,584.5,588.8,592.9,596.9,601.7,609.0,616.3,623.6,630.8,637.8,644.7,713.9,719.8,725.6,731.2,736.9,742.5,748.0,756.6,766.6,777.0,787.3,798.0,802.9,806.5,809.9,813.5,816.9,820.2,823.5,826.8,830.1,833.4,836.7,840.1,843.2,846.2,849.2,864.7,884.0,951.4,954.0,956.6,959.2,961.7,964.2,966.7,969.1,971.5,974.0,976.3,978.7,981.0,983.3,985.6,986.0,986.5,986.9,987.4,987.8,988.3,988.7,989.2,989.6,990.1,990.5,991.0,991.4,991.8,992.3,992.7,993.2,993.6,994.1,994.5,994.9,995.4,995.8,996.2,996.7,997.1,997.5,998.0,998.4,998.8,999.2,999.7,1000.1,1000.6,1001.1,1001.6,1002.1,1002.6,1003.1,1003.6,1004.1,1004.6,1005.1,1005.6,1006.1,1006.6,1007.1,1007.5,1008.0,1008.4};

Double_t m_expected_lvjj[178]={300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,362.0,418.0,462.7,496.8,515.1,530.7,546.2,558.2,569.0,578.7,589.4,598.8,608.4,617.5,627.2,636.3,644.2,652.2,658.7,665.4,671.6,676.8,682.4,688.1,694.1,699.4,703.2,706.6,709.9,713.6,716.5,719.7,723.0,725.7,728.9,731.8,734.4,737.1,740.2,742.4,754.5,764.1,772.7,781.0,788.1,795.2,802.3,808.8,815.4,821.8,827.8,833.2,838.7,843.4,847.8,852.0,856.4,859.9,863.3,866.5,869.4,872.4,874.9,877.6,879.8,881.8,883.7,885.2,886.7,888.0,889.1,890.2,891.3,892.4,892.6,893.6,893.9,894.8,894.4,894.8,894.4,894.8,893.9,893.6,892.6,892.4,891.3,890.2,889.1,888.0,886.7,885.2,883.7,881.8,879.8,877.6,874.9,872.4,869.4,866.5,863.3,859.9,856.4,852.0,847.8,843.4,838.7,833.2,827.8,821.8,815.4,808.8,802.3,795.2,788.1,781.0,772.7,764.1,754.5,742.4,740.2,737.1,734.4,731.8,728.9,725.7,723.0,719.7,716.5,713.6,709.9,706.6,703.2,699.4,694.1,688.1,682.4,676.8,671.6,665.4,658.7,652.2,644.2,636.3,627.2,617.5,608.4,598.8,589.4,578.7,569.0,558.2,546.2,530.7,515.1,496.8,462.7,418.0,362.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0};

Double_t m_observed_lvjj[178]={300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,363.2,380.5,398.8,424.8,511.1,515.7,520.1,524.3,528.4,532.2,536.0,539.5,542.9,546.2,549.4,552.5,555.6,558.5,561.3,564.0,566.7,569.3,571.8,574.2,576.5,578.8,581.1,583.3,585.5,587.7,589.8,591.8,593.7,595.6,597.5,599.3,601.2,611.0,619.8,628.1,635.5,642.4,648.7,770.3,775.6,780.5,785.1,789.3,793.4,797.2,800.9,804.9,808.6,812.1,815.4,818.6,821.5,824.2,826.7,829.0,831.3,833.4,835.5,837.3,839.0,840.5,841.8,843.0,844.1,845.0,845.8,846.4,847.0,847.4,847.7,847.9,847.9,847.9,847.7,847.4,847.0,846.4,845.8,845.0,844.1,843.0,841.8,840.5,839.0,837.3,835.5,833.4,831.3,829.0,826.7,824.2,821.5,818.6,815.4,812.1,808.6,804.9,800.9,797.2,793.4,789.3,785.1,780.5,775.6,770.3,648.7,642.4,635.5,628.1,619.8,611.0,601.2,599.3,597.5,595.6,593.7,591.8,589.8,587.7,585.5,583.3,581.1,578.8,576.5,574.2,571.8,569.3,566.7,564.0,561.3,558.5,555.6,552.5,549.4,546.2,542.9,539.5,536.0,532.2,528.4,524.3,520.1,515.7,511.1,424.8,398.8,380.5,363.2,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0,300.0};

// INFORMATION PASTING OVER

Double_t beta_valsdub[numbetas*2];

int ind = numbetas;

for (unsigned int ii = 0; ii < numbetas*2; ++ii){

if (ii < numbetas) beta_valsdub[ii] = beta_vals[ii];

else {

ind = ind -1;

beta_valsdub[ii] = beta_vals[ind];

}

}

for (unsigned int ii = 0; ii < numbetas; ++ii){

if (beta_vals[ii] == 0.5) {

std::cout<<"## ---------------- B = 1/2 Limits ----------------------"<<std::endl;

std::cout<<"## lvjj Expected : "<<m_expected_lvjj[ii]<<std::endl;

std::cout<<"## lvjj Observed : "<<m_observed_lvjj[ii]<<std::endl;

std::cout<<"## combo Expected : "<<m_expected_combo[ii]<<std::endl;

std::cout<<"## combo Observed : "<<m_observed_combo[ii]<<std::endl;

std::cout<<"## "<<std::endl;

}

if (beta_vals[ii] > 0.999) {

std::cout<<"## ---------------- B = 1 Limits ----------------------"<<std::endl;

std::cout<<"## lljj Expected : "<<m_expected_lljj[ii]<<std::endl;

std::cout<<"## lljj Observed : "<<m_observed_lljj[ii]<<std::endl;

std::cout<<"## combo Expected : "<<m_expected_combo[ii]<<std::endl;

std::cout<<"## combo Observed : "<<m_observed_combo[ii]<<std::endl;

std::cout<<"## "<<std::endl;

}

}

// filename for the final plot (NB: changing the name extension changes the file format)

string fileName0 = "LQ1Combination.eps";

string fileName1 = "LQ1Combination.pdf";

string fileName2 = "LQ1Combination.png";

// axes labels for the final plot

string title = ";M_{LQ} [GeV];#beta";

// integrated luminosity

string sqrts = "#sqrt{s} = 8 TeV";

// region excluded by Tevatron limits (1 fb-1)

Double_t x_excl_D0[11] = {214.39,235.13,254.08,268.12,275.92,283.95,289.08,293.09,295.99,297.10,298.89};

Double_t y_excl_D0[11] = {0,0.10,0.20,0.30,0.40,0.50,0.60,0.70,0.80,0.90,1};

// region excluded by ATLAS limits (35 pb-1)

Double_t x_excl_ATLAS[13] = {200,221.91,243.07,254.12,265.17,287.83,304.68,318.54,331.46,342.88,353.75,363.11,375.47};

Double_t y_excl_ATLAS[13] = {0.095,0.10,0.126,0.151,0.20,0.30,0.40,0.50,0.60,0.70,0.80,0.90,1.0};

// region excluded by Tevatron limits (5.4 fb-1)

Double_t x_excl_D0_new[11] = {243.13,275.12,308.18,320.02,325.95,327.96,325.95,320.02,308.18,275.12,243.13};

Double_t y_excl_D0_new[11] = {0.05,0.10,0.20,0.30,0.40,0.50,0.60,0.70,0.80,0.90,0.95};

//Double_t AtlasComboMass[13] = {250,253.997194950912,278.541374474053,305.539971949509,341.374474053296,352.173913043478,362.973352033661,371.318373071529,385.063113604488,395.371669004208,413.534361851332,421.879382889201,250};

//Double_t AtlasComboBeta[13] ={.1,0.12280701754386,0.14619883040936,0.19493177387914,0.34892787524366,0.39766081871345,0.51851851851852,0.60038986354776,0.74269005847953,0.79727095516569,0.93567251461988,0.99805068226121,1.0};

// Double_t AtlasComboMass[34] = {250.88,284.14,00.0,332.10,354.57,357.86,361.14,380.60,401.43,417.33,437.61,459.81,463.92,477.07,493.79,511.05,521.19,541.75,553.81,565.86,581.49,596.29,613.55,629.45,641.50,647.53,653.56,662.61,667.27,674.67,679.87,683.16,250,250};

// Double_t AtlasComboBeta[34] ={0.04541595925297,0.04626485568761,0.04668930390492,0.04711375212224,0.04711375212224,0.04753820033956,0.04881154499151,0.05772495755518,0.06748726655348,0.07555178268251,0.08446519524618,0.09507640067912,0.09847198641766,0.1213921901528,0.14983022071307,0.17996604414261,0.19821731748727,0.26146010186757,0.29966044142615,0.36332767402377,0.4393039049236,0.51273344651952,0.58828522920204,0.65195246179966,0.69694397283531,0.74363327674024,0.78989813242784,0.87054329371817,0.90789473684211,0.94482173174873,0.97453310696095,1.00042444821732,1.00169779286927,0.04499151103565};

// Double_t CMSComboBeta[165]={ 0.034,0.036,0.038,0.04,0.042,0.044,0.046,0.048,0.05,0.052,0.054,0.056,0.058,0.06,0.062,0.064,0.066,0.068,0.07,0.072,0.074,0.076,0.078,0.08,0.082,0.084,0.086,0.088,0.09,0.092,0.094,0.096,0.098,0.1,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.19,0.2,0.21,0.22,0.23,0.24,0.25,0.26,0.27,0.28,0.29,0.3,0.31,0.32,0.33,0.34,0.35,0.36,0.37,0.38,0.39,0.4,0.41,0.42,0.43,0.44,0.45,0.46,0.47,0.48,0.49,0.5,0.51,0.52,0.53,0.54,0.56,0.57,0.58,0.59,0.6,0.61,0.62,0.63,0.64,0.65,0.66,0.67,0.68,0.69,0.7,0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.78,0.79,0.8,0.81,0.82,0.83,0.84,0.85,0.86,0.87,0.88,0.89,0.9,0.902,0.904,0.906,0.908,0.91,0.912,0.914,0.916,0.918,0.92,0.922,0.924,0.926,0.928,0.93,0.932,0.934,0.936,0.938,0.94,0.942,0.944,0.946,0.948,0.95,0.952,0.954,0.956,0.958,0.96,0.962,0.964,0.966,0.968,0.97,0.972,0.974,0.976,0.978,0.98,0.982,0.984,0.986,0.988,0.99,0.992,0.994,0.996,0.998,0.9995,1.0,0.034};

// Double_t CMSComboMass[165]={ 300.,345.9,356.6,363.4,369.9,376.1,382.1,387.8,393.3,398.7,404.2,409.7,414.9,419.3,424.1,429.0,433.8,438.5,443.0,447.2,451.8,457.2,462.4,467.4,471.9,476.6,481.2,485.6,490.0,494.3,498.5,500.7,502.5,504.3,512.4,520.0,527.0,533.3,539.1,545.0,550.1,561.5,573.1,583.2,592.5,600.2,602.4,604.7,607.0,609.1,611.1,613.1,615.0,617.0,618.7,620.5,622.5,624.0,625.7,627.6,628.9,630.8,632.5,634.0,635.6,637.2,638.8,640.3,641.9,643.5,645.2,646.8,648.1,649.8,653.8,657.2,661.5,665.8,673.7,678.0,682.3,686.5,691.2,695.5,699.9,704.8,710.2,714.8,719.6,725.0,730.0,734.8,740.9,746.8,750.5,754.0,757.9,760.9,764.6,767.6,771.4,774.0,777.9,781.0,784.2,787.1,791.0,794.1,797.3,800.8,804.0,807.6,811.2,811.8,812.2,813.3,814.0,814.4,815.0,815.5,816.2,817.2,817.5,818.4,819.3,819.7,820.5,821.0,821.4,822.8,823.2,823.5,824.5,825.4,825.7,826.3,827.3,828.0,828.5,829.2,829.9,830.8,831.0,831.8,832.8,832.7,834.2,834.0,834.7,835.8,836.1,837.0,838.0,838.4,838.6,839.3,840.5,841.2,841.6,841.6,842.1,843.5,844.0,300,300};

static const int n_exclude = 204;

Double_t beta_CMSATLAS[n_exclude]={0.018,0.018,0.02,0.022,0.024,0.026,0.028,0.03,0.032,0.034,0.036,0.038,0.04,0.042,0.044,0.046,0.048,0.05,0.052,0.054,0.056,0.058,0.06,0.062,0.064,0.066,0.068,0.07,0.072,0.074,0.076,0.078,0.08,0.082,0.084,0.086,0.088,0.09,0.092,0.094,0.096,0.098,0.1,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.19,0.2,0.21,0.22,0.23,0.24,0.26,0.27,0.28,0.29,0.3,0.31,0.32,0.33,0.34,0.35,0.36,0.37,0.38,0.39,0.4,0.41,0.42,0.43,0.44,0.45,0.46,0.47,0.48,0.49,0.5,0.51,0.52,0.53,0.54,0.55,0.56,0.57,0.58,0.59,0.6,0.61,0.62,0.63,0.64,0.65,0.66,0.67,0.68,0.69,0.7,0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.78,0.79,0.8,0.81,0.82,0.83,0.84,0.85,0.86,0.87,0.88,0.89,0.9,0.902,0.904,0.906,0.908,0.91,0.912,0.914,0.916,0.918,0.92,0.922,0.924,0.926,0.928,0.93,0.932,0.934,0.936,0.938,0.94,0.942,0.944,0.946,0.948,0.95,0.952,0.954,0.956,0.958,0.96,0.962,0.964,0.966,0.968,0.97,0.972,0.974,0.976,0.978,0.98,0.982,0.984,0.986,0.988,0.99,0.992,0.994,0.996,0.998,0.9995,0.9995, 0.9459430964276001, 0.9113315376333271, 0.880788245243444, 0.8298887187568689, 0.7789832314403973, 0.7321489909432641, 0.6934572440848077, 0.6506942504070129, 0.552954522593408, 0.5977471788509756, 0.5020579765217812, 0.4572384965296794, 0.4246416982404375, 0.39408648419076153, 0.3676025169604238, 0.3044564654514025, 0.23723022587819786, 0.34112153014503444, 0.26778842034282224, 0.18626214984781253, 0.20463342758895597, 0.16789087210666898, 0.15155521777519476, 0.13521956344372055, 0.11888390911224611, 0.10254527436582372, 0.09028086685368775, 0.07185892205842359, 0.08412332957055957, 0.06366576136562607, 0.05140135385349032};

Double_t mass_CMSATLAS [n_exclude]={300.,317.5,342.2,358.0,369.5,379.8,389.5,398.6,407.1,414.6,422.1,429.1,435.6,441.5,447.3,452.8,457.1,461.1,465.1,468.7,472.3,475.7,479.3,482.9,485.9,488.7,491.5,494.1,496.6,499.9,501.4,502.8,504.1,505.4,506.8,508.0,509.6,510.8,511.9,513.0,514.0,515.6,516.6,521.5,525.9,530.2,533.7,537.3,540.3,543.1,545.7,548.1,552.0,561.4,569.0,576.5,583.9,596.0,600.1,603.2,605.6,607.2,609.4,611.8,613.3,615.7,617.6,619.6,621.8,623.1,625.2,626.6,628.6,630.0,632.2,633.5,635.5,636.8,638.8,640.4,642.5,643.8,645.1,646.9,648.7,651.9,655.3,661.7,666.6,671.8,678.4,685.3,693.6,704.1,716.6,729.2,746.9,751.7,753.7,756.1,758.6,760.7,763.3,765.5,767.7,770.0,772.3,774.7,776.9,779.3,781.7,783.9,786.5,788.6,790.9,793.4,795.5,797.7,800.0,802.5,805.2,807.2,807.8,808.6,808.9,809.5,809.9,810.5,810.4,811.3,811.7,812.2,812.9,813.2,813.9,814.4,814.6,815.3,815.8,816.1,816.6,817.2,817.7,818.1,818.4,818.9,819.6,820.2,820.4,821.2,821.7,822.0,822.3,823.0,823.5,823.7,824.2,824.7,825.0,825.8,826.2,826.8,827.1,827.8,828.3,828.6,829.0,829.9,830.3,830.8,831.1,831.3, 659.1508052708638, 654.4655929721814, 652.1229868228403, 648.6090775988287, 645.0951683748169, 639.2386530014641, 633.3821376281112, 627.5256222547583, 621.6691068814055, 609.9560761346999, 613.4699853587115, 607.6134699853586, 593.5578330893118, 583.0161054172767, 574.8169838945827, 566.6178623718887, 550.2196193265007, 530.3074670571009, 559.5900439238652, 539.6778916544655, 499.8535871156661, 519.7657393850658, 479.94143484626636, 460.02928257686676, 440.11713030746705, 420.2049780380672, 399.121522693997, 379.2093704245973, 339.38506588579787, 359.2972181551976, 319.47291361639816, 299.56076134699845};

Double_t mass_CMS[116]={300.0000,319.8520,331.7728,339.8246,351.0135,359.7973,370.1496,379.8745,389.8086,399.8472,409.4675,420.0290,429.8584,439.4788,449.3082,460.0788,469.4900,478.3784,486.8485,494.2728,498.7693,501.8018,508.4942,515.2912,522.0882,527.5257,532.4405,536.7278,541.8517,545.7207,549.0669,551.5766,557.4324,564.9614,573.7452,580.5421,586.1889,592.5676,597.2732,600.0965,602.0833,603.2336,605.7432,608.1483,610.8671,613.0631,615.0499,617.0367,618.7098,621.1149,623.7291,625.9250,627.4936,630.2124,633.2449,635.5454,638.2642,640.7738,642.9698,645.8977,647.6754,649.4530,650.0804,653.4266,656.4591,658.1322,660.8510,665.2429,669.6348,673.9221,677.7912,682.0785,685.5293,689.3983,694.5222,699.1232,704.2471,709.3710,713.7629,719.3050,723.5923,728.2979,730.9122,733.5264,736.8726,741.8919,746.8066,750.7802,755.5904,759.5640,761.1326,764.5833,768.8707,772.1123,773.3671,777.8636,782.9875,786.4382,788.9479,793.4443,798.4636,801.7053,806.2017,810.2799,814.7764,817.1815,823.3510,827.6384,832.6577,836.2130,838.1998,840.9186,841.9643,842.2780,300.0000,300.0000};

Double_t beta_CMS[116]={0.0333,0.0337,0.0345,0.0350,0.0357,0.0376,0.0411,0.0445,0.0478,0.0516,0.0550,0.0593,0.0636,0.0681,0.0724,0.0767,0.0804,0.0841,0.0881,0.0914,0.0934,0.0963,0.1043,0.1135,0.1225,0.1303,0.1382,0.1453,0.1541,0.1613,0.1678,0.1710,0.1756,0.1824,0.1903,0.1966,0.2027,0.2098,0.2154,0.2197,0.2271,0.2328,0.2440,0.2547,0.2679,0.2781,0.2896,0.2992,0.3090,0.3218,0.3360,0.3485,0.3591,0.3764,0.3936,0.4092,0.4252,0.4411,0.4563,0.4722,0.4854,0.4974,0.5005,0.5085,0.5167,0.5215,0.5281,0.5378,0.5489,0.5597,0.5697,0.5786,0.5876,0.5956,0.6076,0.6185,0.6279,0.6383,0.6478,0.6586,0.6668,0.6757,0.6811,0.6874,0.6929,0.7014,0.7093,0.7205,0.7330,0.7441,0.7503,0.7589,0.7727,0.7815,0.7886,0.7994,0.8154,0.8281,0.8350,0.8470,0.8626,0.8724,0.8855,0.8974,0.9094,0.9200,0.9356,0.9498,0.9632,0.9746,0.9823,0.9897,0.9961,1.0000,1.0000,0.0333};

static const int n_ex_atlas = 33;

Double_t mass_ATLAS[n_ex_atlas] = {300.,319.472913616,359.297218155,339.385065886,379.209370425,399.121522694,420.204978038,440.117130307,460.029282577,479.941434846,519.765739385,499.853587116,539.677891654,559.590043924,530.307467057,550.219619327,566.617862372,574.816983895,583.016105417,593.557833089,607.613469985,613.469985359,609.956076135,621.669106881,627.525622255,633.382137628,639.238653001,645.095168375,648.609077599,652.122986823,654.465592972,659.150805271, 300.};

Double_t beta_ATLAS[n_ex_atlas] = {0.0514013538535,0.0636657613656,0.0841233295706,0.0718589220584,0.0902808668537,0.102545274366,0.118883909112,0.135219563444,0.151555217775,0.167890872107,0.204633427589,0.186262149848,0.267788420343,0.341121530145,0.237230225878,0.304456465451,0.36760251696,0.394086484191,0.42464169824,0.45723849653,0.502057976522,0.597747178851,0.552954522593,0.650694250407,0.693457244085,0.732148990943,0.77898323144,0.829888718757,0.880788245243,0.911331537633,0.945943096428,0.998881226739, 1.};

// turn on/off batch mode

gROOT->SetBatch(kTRUE);

// set ROOT style

// myStyle();

setTDRStyle();

// gStyle->SetPadLeftMargin(0.14);

gROOT->ForceStyle();

TH2F *bg = new TH2F("bg",title.c_str(), 20, 300., 1350., 100, 0., 1.);

bg->SetStats(kFALSE);

bg->GetXaxis()->CenterTitle();

bg->GetYaxis()->CenterTitle();

bg->SetTitleOffset(1.,"X");

bg->SetTitleOffset(1.05,"Y");

// bg->GetXaxis()->SetNdivisions(505);

TCanvas *c = new TCanvas("c","",800,800);

c->cd();

bg->Draw();

TGraph *gr_excl_D0 = new TGraph(11,x_excl_D0,y_excl_D0);

gr_excl_D0->SetLineWidth(3);

gr_excl_D0->SetLineColor(49);

gr_excl_D0->SetLineStyle(6);

//gr_excl_D0->Draw("C");

TPolyLine*gr_excl_ATLAS = new TPolyLine(n_ex_atlas,mass_ATLAS,beta_ATLAS,"");

//gr_excl_ATLAS->SetLineWidth(3);

gr_excl_ATLAS->SetLineColor(8);

gr_excl_ATLAS->SetFillColor(8);

gr_excl_ATLAS->SetFillStyle(3345);

TPolyLine*gr_excl_CMS = new TPolyLine(116,mass_CMS,beta_CMS,"");

//gr_excl_CMS->SetLineWidth(3);

gr_excl_CMS->SetLineColor(14);

gr_excl_CMS->SetFillColor(14);

gr_excl_CMS->SetFillStyle(3344);

TPolyLine*gr_excl_CMSborder = new TPolyLine(n_exclude,mass_CMSATLAS,beta_CMSATLAS,"");

//gr_excl_CMS->SetLineWidth(3);

gr_excl_CMSborder->SetLineColor(14);

gr_excl_CMSborder->SetFillColor(14);

gr_excl_CMSborder->SetFillStyle(3344);

TGraph *gr_excl_D0_new = new TGraph(11,x_excl_D0_new,y_excl_D0_new);

gr_excl_D0_new->SetLineWidth(3);

gr_excl_D0_new->SetLineColor(30);

gr_excl_D0_new->SetLineStyle(7);

//gr_excl_D0_new->Draw("C");

// ------------------------------------ LLJJ CURVES -------------------------------------//

TPolyLine *beta_vs_m_lljj_observed = new TPolyLine(numbetas, m_observed_lljj, beta_vals);

TPolyLine *beta_vs_m_lljj_expected = new TPolyLine(numbetas, m_expected_lljj, beta_vals);

//TPolyLine *beta_vs_m_lljj_observed_shade = new TPolyLine(numbetas*2, m_2sigma_lljj, beta_valsdub);

beta_vs_m_lljj_observed->SetLineWidth(2);

beta_vs_m_lljj_observed->SetLineColor(kBlue);

beta_vs_m_lljj_observed->SetLineStyle(1);

beta_vs_m_lljj_expected->SetLineWidth(2);

beta_vs_m_lljj_expected->SetLineColor(kBlue);

beta_vs_m_lljj_expected->SetLineStyle(2);

//beta_vs_m_lljj_observed_shade->SetLineWidth(3);

//beta_vs_m_lljj_observed_shade->SetLineColor(kGreen+2);

//beta_vs_m_lljj_observed_shade->SetLineStyle(4);

//beta_vs_m_lljj_observed_shade->SetFillStyle(3425);

//beta_vs_m_lljj_observed_shade->SetFillColor(kGreen+2);

//beta_vs_m_lljj_observed_shade->Draw("f");

// ------------------------------------ LVJJ CURVES -------------------------------------//

TPolyLine *beta_vs_m_lvjj_observed = new TPolyLine(numbetas, m_observed_lvjj, beta_vals);

TPolyLine *beta_vs_m_lvjj_expected = new TPolyLine(numbetas, m_expected_lvjj, beta_vals);

//TPolyLine *beta_vs_m_lvjj_observed_shade = new TPolyLine(numbetas*2, m_2sigma_lvjj, beta_valsdub);

beta_vs_m_lvjj_observed->SetLineWidth(2);

beta_vs_m_lvjj_observed->SetLineColor(kRed);

beta_vs_m_lvjj_observed->SetLineStyle(1);

beta_vs_m_lvjj_expected->SetLineWidth(2);

beta_vs_m_lvjj_expected->SetLineColor(kRed);

beta_vs_m_lvjj_expected->SetLineStyle(2);

//beta_vs_m_lvjj_observed_shade->SetFillStyle(3452);

//beta_vs_m_lvjj_observed_shade->SetFillColor(kRed);

//beta_vs_m_lvjj_observed_shade->Draw("f");

// ------------------------------------ COMBO CURVES -------------------------------------//

TPolyLine *beta_vs_m_comb_observed = new TPolyLine(numbetas, m_observed_combo, beta_vals);

TPolyLine *beta_vs_m_comb_expected = new TPolyLine(numbetas, m_expected_combo, beta_vals);

TGraph *beta_vs_m_comb_expected_shade2 = new TGraph(numbetas*2, m_2sigma_combo, beta_valsdub);

TGraph *beta_vs_m_comb_expected_shade1 = new TGraph(numbetas*2, m_1sigma_combo, beta_valsdub);

beta_vs_m_comb_observed->SetLineWidth(4);

beta_vs_m_comb_observed->SetLineColor(kBlack);

beta_vs_m_comb_expected->SetLineStyle(1);

//beta_vs_m_comb_expected_shade->SetLineWidth(1);

beta_vs_m_comb_expected_shade1->SetFillStyle(1001);

beta_vs_m_comb_expected_shade1->SetFillColor(kGreen);

beta_vs_m_comb_expected_shade2->SetFillStyle(1001);

beta_vs_m_comb_expected_shade2->SetFillColor(kYellow);

beta_vs_m_comb_expected->SetLineWidth(4);

beta_vs_m_comb_expected->SetLineColor(kBlack);

beta_vs_m_comb_expected->SetLineStyle(2);

// ------------------------------------ DRAW -------------------------------------//

beta_vs_m_comb_expected_shade2->Draw("f");

beta_vs_m_comb_expected_shade1->Draw("f");

gr_excl_ATLAS->Draw("L");

gr_excl_ATLAS->Draw("F");

gr_excl_CMSborder->Draw("L");

gr_excl_CMSborder->Draw("F");

beta_vs_m_lljj_expected->Draw("C");

beta_vs_m_lvjj_expected->Draw("C");

beta_vs_m_comb_expected->Draw("C");

beta_vs_m_lljj_observed->Draw("C");

beta_vs_m_lvjj_observed->Draw("C");

beta_vs_m_comb_observed->Draw("C");

gPad->RedrawAxis();

TLegend *legend = new TLegend(.63,.15,.93,.44);

legend->SetBorderSize(1);

legend->SetFillColor(0);

//legend->SetFillStyle(0);

legend->SetTextFont(42);

legend->SetMargin(0.15);

legend->SetHeader("95% CL limits");

legend->AddEntry(beta_vs_m_comb_observed,"CMS eejj + e#nujj (Obs.)","l");

legend->AddEntry(beta_vs_m_comb_expected,"CMS eejj + e#nujj (Exp.)","l");

legend->AddEntry(beta_vs_m_lvjj_observed,"CMS e#nujj (Obs.)","l");

legend->AddEntry(beta_vs_m_lvjj_expected,"CMS e#nujj (Exp.)","l");

legend->AddEntry(beta_vs_m_lljj_observed,"CMS eejj (Obs.)","l");

legend->AddEntry(beta_vs_m_lljj_expected,"CMS eejj (Exp.)","l");

legend->AddEntry(gr_excl_ATLAS,"ATLAS, 7 TeV, 1.03 fb^{-1}","f");

legend->AddEntry(gr_excl_CMSborder,"CMS, 7 TeV, 5.0 fb^{-1}","f");

//legend->AddEntry(gr_excl_D0,"D#oslash (Obs.), 1 fb^{-1}","l");

legend->Draw();

TLatex l1;

l1.SetTextAlign(12);

l1.SetTextFont(42);

l1.SetNDC();

l1.SetTextSize(0.04);

// l1.DrawLatex(0.7,0.41,"CMS 2012");

// l1.DrawLatex(0.7,0.36,sqrts.c_str());

// l1.DrawLatex(0.7,0.47,"CMS Preliminary");

// l1.DrawLatex(0.7,0.41,"Preliminary");

// l1.DrawLatex(0.7,0.36,sqrts.c_str());

// l1.DrawLatex(0.14,0.93,"CMS Work In Progress #sqrt{s} = 8 TeV 19.7 fb^{-1}");

//l1.DrawLatex(0.14,0.93,"CMS #it{Preliminary} #sqrt{s} = 8 TeV 19.7 fb^{-1}");

l1.DrawLatex(0.14,0.93,"CMS #sqrt{s} = 8 TeV 19.7 fb^{-1}");

TLatex l2;

l2.SetTextAlign(12);

l2.SetTextFont(42);

l2.SetNDC();

l2.SetTextSize(0.04);

l2.SetTextColor(kBlue);

l2.SetTextAngle(30);

l2.DrawLatex(0.24,0.35,"eejj");

l2.SetTextColor(kRed);

l2.SetTextAngle(-20);

l2.DrawLatex(0.25,0.84,"e#nujj");

l2.SetTextColor(kBlack);

l2.SetTextAngle(79);

l2.DrawLatex(0.65,0.50,"eejj + e#nujj");

// c->SetGridx();

// c->SetGridy();

c->SetLogy(0);

c->SaveAs(fileName0.c_str());

c->SaveAs(fileName1.c_str());

c->SaveAs(fileName2.c_str());

delete legend;

delete beta_vs_m_lljj_observed;

delete beta_vs_m_lvjj_observed;

delete beta_vs_m_comb_observed;

delete gr_excl_D0;

delete gr_excl_ATLAS;

delete gr_excl_D0_new;

delete bg;

delete c;

}

void myStyle()

{

gStyle->Reset("Default");

gStyle->SetCanvasColor(0);

gStyle->SetPadColor(0);

gStyle->SetTitleFillColor(10);

gStyle->SetCanvasBorderMode(0);

gStyle->SetStatColor(0);

gStyle->SetPadBorderMode(0);

gStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame

gStyle->SetPadTickY(1);

gStyle->SetFrameBorderMode(0);

gStyle->SetPalette(1);

//gStyle->SetOptStat(kFALSE);

gStyle->SetOptStat(111110);

gStyle->SetOptFit(0);

gStyle->SetStatFont(42);

gStyle->SetPadLeftMargin(0.13);

gStyle->SetPadRightMargin(0.07);

// gStyle->SetTitleFont(42);

// gStyle->SetTitleFont(42, "XYZ");

// gStyle->SetLabelFont(42, "XYZ");

gStyle->SetStatY(.9);

}

void setTDRStyle() {

TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");

// For the canvas:

tdrStyle->SetCanvasBorderMode(0);

tdrStyle->SetCanvasColor(kWhite);

tdrStyle->SetCanvasDefH(600); //Height of canvas

tdrStyle->SetCanvasDefW(600); //Width of canvas

tdrStyle->SetCanvasDefX(0); //POsition on screen

tdrStyle->SetCanvasDefY(0);

// For the Pad:

tdrStyle->SetPadBorderMode(0);

// tdrStyle->SetPadBorderSize(Width_t size = 1);

tdrStyle->SetPadColor(kWhite);

tdrStyle->SetPadGridX(false);

tdrStyle->SetPadGridY(false);

tdrStyle->SetGridColor(0);

tdrStyle->SetGridStyle(3);

tdrStyle->SetGridWidth(1);

// For the frame:

tdrStyle->SetFrameBorderMode(0);

tdrStyle->SetFrameBorderSize(1);

tdrStyle->SetFrameFillColor(0);

tdrStyle->SetFrameFillStyle(0);

tdrStyle->SetFrameLineColor(1);

tdrStyle->SetFrameLineStyle(1);

tdrStyle->SetFrameLineWidth(1);

// For the histo:

tdrStyle->SetHistFillColor(63);

// tdrStyle->SetHistFillStyle(0);

tdrStyle->SetHistLineColor(1);

tdrStyle->SetHistLineStyle(0);

tdrStyle->SetHistLineWidth(1);

tdrStyle->SetMarkerStyle(20);

//For the fit/function:

tdrStyle->SetOptFit(1);

tdrStyle->SetFitFormat("5.4g");

tdrStyle->SetFuncColor(2);

tdrStyle->SetFuncStyle(1);

tdrStyle->SetFuncWidth(1);

//For the date:

tdrStyle->SetOptDate(0);

// For the statistics box:

tdrStyle->SetOptFile(0);

tdrStyle->SetOptStat(0); // To display the mean and RMS: SetOptStat("mr");

tdrStyle->SetStatColor(kWhite);

tdrStyle->SetStatFont(42);

tdrStyle->SetStatFontSize(0.025);

tdrStyle->SetStatTextColor(1);

tdrStyle->SetStatFormat("6.4g");

tdrStyle->SetStatBorderSize(1);

tdrStyle->SetStatH(0.1);

tdrStyle->SetStatW(0.15);

// Margins:

tdrStyle->SetPadTopMargin(0.1);

tdrStyle->SetPadBottomMargin(0.13);

tdrStyle->SetPadLeftMargin(0.13);

tdrStyle->SetPadRightMargin(0.05);

// For the Global title:

// tdrStyle->SetOptTitle(0);

tdrStyle->SetTitleFont(42);

tdrStyle->SetTitleColor(1);

tdrStyle->SetTitleTextColor(1);

tdrStyle->SetTitleFillColor(10);

tdrStyle->SetTitleFontSize(0.05);

// For the axis titles:

tdrStyle->SetTitleColor(1, "XYZ");

tdrStyle->SetTitleFont(42, "XYZ");

tdrStyle->SetTitleSize(0.06, "XYZ");

tdrStyle->SetTitleXOffset(0.9);

tdrStyle->SetTitleYOffset(1.05);

// For the axis labels:

tdrStyle->SetLabelColor(1, "XYZ");

tdrStyle->SetLabelFont(42, "XYZ");

tdrStyle->SetLabelOffset(0.007, "XYZ");

tdrStyle->SetLabelSize(0.04, "XYZ");

// For the axis:

tdrStyle->SetAxisColor(1, "XYZ");

tdrStyle->SetStripDecimals(kTRUE);

tdrStyle->SetTickLength(0.03, "XYZ");

tdrStyle->SetNdivisions(510, "XYZ");

tdrStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame

tdrStyle->SetPadTickY(1);

// Change for log plots:

tdrStyle->SetOptLogx(0);

tdrStyle->SetOptLogy(0);

tdrStyle->SetOptLogz(0);

tdrStyle->cd();

}