#include "TTree.h" #include "TFile.h" #include "TH1D.h" #include "TMath.h" #include "RooDataSet.h" #include "RooRealVar.h" #include "RooConstVar.h" #include "RooHistFunc.h" #include "RooDataHist.h" #include "RooProduct.h" #include "RooPoisson.h" #include "RooRealIntegral.h" #include "RooProjectedPdf.h" #include "RooWorkspace.h" #include "RooNumIntConfig.h" #include "RooStats/ModelConfig.h" #include "RooStats/FeldmanCousins.h" using namespace RooFit; using namespace RooStats; void GetInterval() { //Constant Double_t N_avo = 6.022E23 ; //mol^-1 Double_t pot = 5.3 ; //10^19 pot Double_t mass = 1.179 ; //kton Double_t konst = 1.E9 ; //mol/kton //Constant RooConstVar L ("L" ,"Baseline" ,732. ); RooConstVar theta13 ("theta13" ,"#theta_{13}" ,9. ); RooConstVar theta23 ("theta23" ,"#theta_{23}" ,40. ); RooConstVar dm2_atm ("dm2_atm" ,"#Deltam^{2}_{atm}" ,2.43E-3 ); //Variables RooRealVar Nobs ("Nobs" ,"Nobs" ,50. ,0. ,100. ,"events" ); RooRealVar Enu ("Enu" ,"Enu" ,200. ,0. ,400. ,"GeV" ); RooRealVar theta14 ("theta14" ,"#theta_{14}" ,12.5 ,0. ,25. ,"degree" ); RooRealVar theta24 ("theta24" ,"#theta_{24}" ,12.5 ,0. ,25. ,"degree" ); RooRealVar delta ("delta" ,"#delta_{2} - #delta_{3}" ,TMath::Pi() ,0. ,2*TMath::Pi() ,"rad" ); //Data TFile fin("data.root","RECREATE"); TTree* t = new TTree("t","t"); Double_t events = 6; t->Branch("Nobs",&events,"Nobs/D"); t->Fill(); RooDataSet data("data","data",t,RooArgSet(Nobs)); fin.Close(); //Histogram: beam numu flux TH1D *h_numu_flux = new TH1D("h_numu_flux","h_numu_flux",200,0,400); h_numu_flux->SetBinContent(1,29265); h_numu_flux->SetBinContent(2,62355); h_numu_flux->SetBinContent(3,112605); h_numu_flux->SetBinContent(4,220925); h_numu_flux->SetBinContent(5,364350); h_numu_flux->SetBinContent(6,309900); h_numu_flux->SetBinContent(7,309400); h_numu_flux->SetBinContent(8,322725); h_numu_flux->SetBinContent(9,309500); h_numu_flux->SetBinContent(10,328175); h_numu_flux->SetBinContent(11,362900); h_numu_flux->SetBinContent(12,330850); h_numu_flux->SetBinContent(13,229800); h_numu_flux->SetBinContent(14,137700); h_numu_flux->SetBinContent(15,83725); h_numu_flux->SetBinContent(16,49567.5); h_numu_flux->SetBinContent(17,27890); h_numu_flux->SetBinContent(18,16057.5); h_numu_flux->SetBinContent(19,11025.5); h_numu_flux->SetBinContent(20,8391.75); h_numu_flux->SetBinContent(21,7220.75); h_numu_flux->SetBinContent(22,6671.5); h_numu_flux->SetBinContent(23,6012.75); h_numu_flux->SetBinContent(24,5558); h_numu_flux->SetBinContent(25,5079.25); h_numu_flux->SetBinContent(26,4968.25); h_numu_flux->SetBinContent(27,4754.25); h_numu_flux->SetBinContent(28,4356.5); h_numu_flux->SetBinContent(29,4049.25); h_numu_flux->SetBinContent(30,3872); h_numu_flux->SetBinContent(31,3791.25); h_numu_flux->SetBinContent(32,3586.75); h_numu_flux->SetBinContent(33,3449.75); h_numu_flux->SetBinContent(34,3208); h_numu_flux->SetBinContent(35,3090); h_numu_flux->SetBinContent(36,2661.25); h_numu_flux->SetBinContent(37,2550.5); h_numu_flux->SetBinContent(38,2515.75); h_numu_flux->SetBinContent(39,2309.75); h_numu_flux->SetBinContent(40,2183.5); h_numu_flux->SetBinContent(41,2086.5); h_numu_flux->SetBinContent(42,1945.5); h_numu_flux->SetBinContent(43,1839.75); h_numu_flux->SetBinContent(44,1758); h_numu_flux->SetBinContent(45,1694.75); h_numu_flux->SetBinContent(46,1472.5); h_numu_flux->SetBinContent(47,1476.75); h_numu_flux->SetBinContent(48,1344.25); h_numu_flux->SetBinContent(49,1321.5); h_numu_flux->SetBinContent(50,1131); h_numu_flux->SetBinContent(51,1081.93); h_numu_flux->SetBinContent(52,1077.78); h_numu_flux->SetBinContent(53,1016.5); h_numu_flux->SetBinContent(54,972.05); h_numu_flux->SetBinContent(55,937.9); h_numu_flux->SetBinContent(56,978.475); h_numu_flux->SetBinContent(57,946.05); h_numu_flux->SetBinContent(58,865.625); h_numu_flux->SetBinContent(59,768.3); h_numu_flux->SetBinContent(60,669.925); h_numu_flux->SetBinContent(61,695); h_numu_flux->SetBinContent(62,608.35); h_numu_flux->SetBinContent(63,651.7); h_numu_flux->SetBinContent(64,563.3); h_numu_flux->SetBinContent(65,558.175); h_numu_flux->SetBinContent(66,524.25); h_numu_flux->SetBinContent(67,594.175); h_numu_flux->SetBinContent(68,483.125); h_numu_flux->SetBinContent(69,484.425); h_numu_flux->SetBinContent(70,470.325); h_numu_flux->SetBinContent(71,450.675); h_numu_flux->SetBinContent(72,449.35); h_numu_flux->SetBinContent(73,417.425); h_numu_flux->SetBinContent(74,402.475); h_numu_flux->SetBinContent(75,336.625); h_numu_flux->SetBinContent(76,341.95); h_numu_flux->SetBinContent(77,293.275); h_numu_flux->SetBinContent(78,269.325); h_numu_flux->SetBinContent(79,282.525); h_numu_flux->SetBinContent(80,236.025); h_numu_flux->SetBinContent(81,236.4); h_numu_flux->SetBinContent(82,244.875); h_numu_flux->SetBinContent(83,210.325); h_numu_flux->SetBinContent(84,194.475); h_numu_flux->SetBinContent(85,182.175); h_numu_flux->SetBinContent(86,197.025); h_numu_flux->SetBinContent(87,161.5); h_numu_flux->SetBinContent(88,146.812); h_numu_flux->SetBinContent(89,152.375); h_numu_flux->SetBinContent(90,132.65); h_numu_flux->SetBinContent(91,112.47); h_numu_flux->SetBinContent(92,113.565); h_numu_flux->SetBinContent(93,97.0575); h_numu_flux->SetBinContent(94,122.25); h_numu_flux->SetBinContent(95,106.89); h_numu_flux->SetBinContent(96,88.0025); h_numu_flux->SetBinContent(97,70.295); h_numu_flux->SetBinContent(98,65.9275); h_numu_flux->SetBinContent(99,91.2825); h_numu_flux->SetBinContent(100,42.98); h_numu_flux->SetBinContent(101,55.1325); h_numu_flux->SetBinContent(102,54.785); h_numu_flux->SetBinContent(103,58.89); h_numu_flux->SetBinContent(104,50.31); h_numu_flux->SetBinContent(105,51.9925); h_numu_flux->SetBinContent(106,44.72); h_numu_flux->SetBinContent(107,47.415); h_numu_flux->SetBinContent(108,39.1); h_numu_flux->SetBinContent(109,34.0525); h_numu_flux->SetBinContent(110,39.765); h_numu_flux->SetBinContent(111,27.94); h_numu_flux->SetBinContent(112,46.995); h_numu_flux->SetBinContent(113,32.875); h_numu_flux->SetBinContent(114,34.7475); h_numu_flux->SetBinContent(115,36.1725); h_numu_flux->SetBinContent(116,34.97); h_numu_flux->SetBinContent(117,25.795); h_numu_flux->SetBinContent(118,20.9); h_numu_flux->SetBinContent(119,22.3117); h_numu_flux->SetBinContent(120,17.445); h_numu_flux->SetBinContent(121,22.93); h_numu_flux->SetBinContent(122,17.3705); h_numu_flux->SetBinContent(123,16.5703); h_numu_flux->SetBinContent(124,18.86); h_numu_flux->SetBinContent(125,19.4097); h_numu_flux->SetBinContent(126,17.3035); h_numu_flux->SetBinContent(127,10.3575); h_numu_flux->SetBinContent(128,11.6818); h_numu_flux->SetBinContent(129,11.953); h_numu_flux->SetBinContent(130,8.9455); h_numu_flux->SetBinContent(131,11.6255); h_numu_flux->SetBinContent(132,14.835); h_numu_flux->SetBinContent(133,14.1835); h_numu_flux->SetBinContent(134,6.4355); h_numu_flux->SetBinContent(135,6.0152); h_numu_flux->SetBinContent(136,7.386); h_numu_flux->SetBinContent(137,3.4675); h_numu_flux->SetBinContent(138,9.02775); h_numu_flux->SetBinContent(139,8.11025); h_numu_flux->SetBinContent(140,7.34445); h_numu_flux->SetBinContent(141,12.2634); h_numu_flux->SetBinContent(142,9.34225); h_numu_flux->SetBinContent(143,7.58575); h_numu_flux->SetBinContent(144,5.33437); h_numu_flux->SetBinContent(145,5.43255); h_numu_flux->SetBinContent(146,4.035); h_numu_flux->SetBinContent(147,4.8075); h_numu_flux->SetBinContent(148,1.3335); h_numu_flux->SetBinContent(149,2.58123); h_numu_flux->SetBinContent(150,2.50257); h_numu_flux->SetBinContent(151,2.06363); h_numu_flux->SetBinContent(152,3.34523); h_numu_flux->SetBinContent(153,4.31787); h_numu_flux->SetBinContent(154,1.89035); h_numu_flux->SetBinContent(155,2.868); h_numu_flux->SetBinContent(156,3.27398); h_numu_flux->SetBinContent(157,3.0241); h_numu_flux->SetBinContent(158,0.9665); h_numu_flux->SetBinContent(159,1.93845); h_numu_flux->SetBinContent(160,1.34748); h_numu_flux->SetBinContent(161,1.0465); h_numu_flux->SetBinContent(162,1.29); h_numu_flux->SetBinContent(163,1.42654); h_numu_flux->SetBinContent(164,1.7871); h_numu_flux->SetBinContent(165,3.246); h_numu_flux->SetBinContent(166,1.29222); h_numu_flux->SetBinContent(167,0.52688); h_numu_flux->SetBinContent(168,2.28467); h_numu_flux->SetBinContent(169,0.467973); h_numu_flux->SetBinContent(170,0.029395); h_numu_flux->SetBinContent(171,0.877957); h_numu_flux->SetBinContent(172,0.44119); h_numu_flux->SetBinContent(173,0.00665875); h_numu_flux->SetBinContent(174,0.0026765); h_numu_flux->SetBinContent(175,0.429446); h_numu_flux->SetBinContent(176,0.431054); h_numu_flux->SetBinContent(177,0.00486525); h_numu_flux->SetBinContent(178,0.0019575); h_numu_flux->SetBinContent(179,0.000787575); h_numu_flux->SetBinContent(180,0.000316875); h_numu_flux->SetBinContent(181,0.000127488); h_numu_flux->SetBinContent(182,5.12975e-05); h_numu_flux->SetBinContent(183,2.06425e-05); h_numu_flux->SetBinContent(184,8.3035e-06); h_numu_flux->SetBinContent(185,3.3415e-06); h_numu_flux->SetBinContent(186,1.34457e-06); h_numu_flux->SetBinContent(187,5.40875e-07); h_numu_flux->SetBinContent(188,2.17625e-07); h_numu_flux->SetBinContent(189,8.75625e-08); h_numu_flux->SetBinContent(190,3.52325e-08); h_numu_flux->SetBinContent(191,1.4174e-08); h_numu_flux->SetBinContent(192,5.70325e-09); h_numu_flux->SetBinContent(193,2.2945e-09); h_numu_flux->SetBinContent(194,9.2335e-10); h_numu_flux->SetBinContent(195,3.71475e-10); h_numu_flux->SetBinContent(196,1.49475e-10); h_numu_flux->SetBinContent(197,6.0135e-11); h_numu_flux->SetBinContent(198,2.41925e-11); h_numu_flux->SetBinContent(199,9.73325e-12); h_numu_flux->SetBinContent(200,4.1415e-12); //Histogram: beam nue flux TH1D *h_nue_flux = new TH1D("h_nue_flux","h_nue_flux",200,0,400); h_nue_flux->SetBinContent(1,354.5); h_nue_flux->SetBinContent(2,832.7); h_nue_flux->SetBinContent(3,1158); h_nue_flux->SetBinContent(4,1397); h_nue_flux->SetBinContent(5,1527); h_nue_flux->SetBinContent(6,1503); h_nue_flux->SetBinContent(7,1452); h_nue_flux->SetBinContent(8,1391); h_nue_flux->SetBinContent(9,1417); h_nue_flux->SetBinContent(10,1356); h_nue_flux->SetBinContent(11,1108); h_nue_flux->SetBinContent(12,944); h_nue_flux->SetBinContent(13,878.5); h_nue_flux->SetBinContent(14,866.5); h_nue_flux->SetBinContent(15,755.3); h_nue_flux->SetBinContent(16,658.6); h_nue_flux->SetBinContent(17,573); h_nue_flux->SetBinContent(18,506.8); h_nue_flux->SetBinContent(19,467); h_nue_flux->SetBinContent(20,431.5); h_nue_flux->SetBinContent(21,360.9); h_nue_flux->SetBinContent(22,324.1); h_nue_flux->SetBinContent(23,332.6); h_nue_flux->SetBinContent(24,369.4); h_nue_flux->SetBinContent(25,303); h_nue_flux->SetBinContent(26,228.4); h_nue_flux->SetBinContent(27,212.9); h_nue_flux->SetBinContent(28,195.6); h_nue_flux->SetBinContent(29,179); h_nue_flux->SetBinContent(30,159.7); h_nue_flux->SetBinContent(31,150.7); h_nue_flux->SetBinContent(32,148.2); h_nue_flux->SetBinContent(33,144.6); h_nue_flux->SetBinContent(34,117.6); h_nue_flux->SetBinContent(35,100.5); h_nue_flux->SetBinContent(36,94.73); h_nue_flux->SetBinContent(37,81.31); h_nue_flux->SetBinContent(38,74.51); h_nue_flux->SetBinContent(39,68.76); h_nue_flux->SetBinContent(40,61.92); h_nue_flux->SetBinContent(41,59.66); h_nue_flux->SetBinContent(42,58.27); h_nue_flux->SetBinContent(43,55.12); h_nue_flux->SetBinContent(44,47.34); h_nue_flux->SetBinContent(45,44.66); h_nue_flux->SetBinContent(46,37.26); h_nue_flux->SetBinContent(47,34.4); h_nue_flux->SetBinContent(48,34.02); h_nue_flux->SetBinContent(49,31.17); h_nue_flux->SetBinContent(50,34.3); h_nue_flux->SetBinContent(51,31.32); h_nue_flux->SetBinContent(52,25.4); h_nue_flux->SetBinContent(53,23.2); h_nue_flux->SetBinContent(54,19.27); h_nue_flux->SetBinContent(55,17.54); h_nue_flux->SetBinContent(56,16.8); h_nue_flux->SetBinContent(57,14.77); h_nue_flux->SetBinContent(58,14.22); h_nue_flux->SetBinContent(59,16.95); h_nue_flux->SetBinContent(60,15.39); h_nue_flux->SetBinContent(61,14.43); h_nue_flux->SetBinContent(62,14.25); h_nue_flux->SetBinContent(63,13.44); h_nue_flux->SetBinContent(64,8.715); h_nue_flux->SetBinContent(65,10.84); h_nue_flux->SetBinContent(66,8.536); h_nue_flux->SetBinContent(67,7.087); h_nue_flux->SetBinContent(68,5.618); h_nue_flux->SetBinContent(69,5.768); h_nue_flux->SetBinContent(70,7.695); h_nue_flux->SetBinContent(71,7.987); h_nue_flux->SetBinContent(72,6.508); h_nue_flux->SetBinContent(73,5.593); h_nue_flux->SetBinContent(74,3.607); h_nue_flux->SetBinContent(75,3.364); h_nue_flux->SetBinContent(76,4.575); h_nue_flux->SetBinContent(77,3.094); h_nue_flux->SetBinContent(78,3.387); h_nue_flux->SetBinContent(79,4); h_nue_flux->SetBinContent(80,3.812); h_nue_flux->SetBinContent(81,3.361); h_nue_flux->SetBinContent(82,1.737); h_nue_flux->SetBinContent(83,1.427); h_nue_flux->SetBinContent(84,1.649); h_nue_flux->SetBinContent(85,2.085); h_nue_flux->SetBinContent(86,1.577); h_nue_flux->SetBinContent(87,1.062); h_nue_flux->SetBinContent(88,0.9862); h_nue_flux->SetBinContent(89,1.225); h_nue_flux->SetBinContent(90,0.611); h_nue_flux->SetBinContent(91,1.464); h_nue_flux->SetBinContent(92,1.116); h_nue_flux->SetBinContent(93,1.301); h_nue_flux->SetBinContent(94,1.618); h_nue_flux->SetBinContent(95,0.9483); h_nue_flux->SetBinContent(96,0.09176); h_nue_flux->SetBinContent(97,0.4581); h_nue_flux->SetBinContent(98,0.7964); h_nue_flux->SetBinContent(99,0.354); h_nue_flux->SetBinContent(100,0.0253); h_nue_flux->SetBinContent(101,0.5083); h_nue_flux->SetBinContent(102,0.8367); h_nue_flux->SetBinContent(103,1.23); h_nue_flux->SetBinContent(104,0.7847); h_nue_flux->SetBinContent(105,0.1652); h_nue_flux->SetBinContent(106,0.02558); h_nue_flux->SetBinContent(107,0.2008); h_nue_flux->SetBinContent(108,0.5459); h_nue_flux->SetBinContent(109,0.4325); h_nue_flux->SetBinContent(110,0.2735); h_nue_flux->SetBinContent(111,0.2654); h_nue_flux->SetBinContent(112,0.4196); h_nue_flux->SetBinContent(113,0.4373); h_nue_flux->SetBinContent(114,0.1098); h_nue_flux->SetBinContent(115,0.0007707); h_nue_flux->SetBinContent(116,0.00685); h_nue_flux->SetBinContent(117,0.005641); h_nue_flux->SetBinContent(118,0.007542); h_nue_flux->SetBinContent(119,0.0005112); h_nue_flux->SetBinContent(120,0.002875); h_nue_flux->SetBinContent(121,0.002212); h_nue_flux->SetBinContent(122,0.08301); h_nue_flux->SetBinContent(123,0.4561); h_nue_flux->SetBinContent(124,0.249); h_nue_flux->SetBinContent(125,0.00179); h_nue_flux->SetBinContent(126,0.00234); h_nue_flux->SetBinContent(127,0.2371); h_nue_flux->SetBinContent(128,0.1614); h_nue_flux->SetBinContent(129,0.0009274); h_nue_flux->SetBinContent(130,0.0007276); h_nue_flux->SetBinContent(131,0.0005537); h_nue_flux->SetBinContent(132,0.0005096); h_nue_flux->SetBinContent(133,0.0004217); h_nue_flux->SetBinContent(134,0.0003407); h_nue_flux->SetBinContent(135,0.0004183); h_nue_flux->SetBinContent(136,0.0001039); h_nue_flux->SetBinContent(137,0.0002196); h_nue_flux->SetBinContent(138,0.0001604); h_nue_flux->SetBinContent(139,0.0001669); h_nue_flux->SetBinContent(140,7.514e-05); h_nue_flux->SetBinContent(141,9.849e-05); h_nue_flux->SetBinContent(142,9.063e-05); h_nue_flux->SetBinContent(143,9.573e-05); h_nue_flux->SetBinContent(144,6.719e-05); h_nue_flux->SetBinContent(145,3.302e-05); h_nue_flux->SetBinContent(146,1.969e-05); h_nue_flux->SetBinContent(147,2.113e-05); h_nue_flux->SetBinContent(148,1.975e-05); h_nue_flux->SetBinContent(149,1.024e-05); h_nue_flux->SetBinContent(150,1.331e-05); h_nue_flux->SetBinContent(151,8.607e-06); h_nue_flux->SetBinContent(152,1.617e-05); h_nue_flux->SetBinContent(153,1.274e-05); h_nue_flux->SetBinContent(154,3.436e-06); h_nue_flux->SetBinContent(155,6.105e-06); h_nue_flux->SetBinContent(156,3.696e-06); h_nue_flux->SetBinContent(157,1.03e-05); h_nue_flux->SetBinContent(158,1.956e-05); h_nue_flux->SetBinContent(159,5.567e-06); h_nue_flux->SetBinContent(160,3.618e-06); h_nue_flux->SetBinContent(161,2.9e-06); h_nue_flux->SetBinContent(162,8.627e-07); h_nue_flux->SetBinContent(163,2.506e-06); h_nue_flux->SetBinContent(164,4.056e-06); h_nue_flux->SetBinContent(165,3.818e-06); h_nue_flux->SetBinContent(166,2.397e-06); h_nue_flux->SetBinContent(167,4.321e-07); h_nue_flux->SetBinContent(168,1.444e-07); h_nue_flux->SetBinContent(169,5.81e-08); h_nue_flux->SetBinContent(170,2.338e-08); h_nue_flux->SetBinContent(171,9.405e-09); h_nue_flux->SetBinContent(172,3.784e-09); h_nue_flux->SetBinContent(173,1.523e-09); h_nue_flux->SetBinContent(174,6.126e-10); h_nue_flux->SetBinContent(175,2.465e-10); h_nue_flux->SetBinContent(176,9.917e-11); h_nue_flux->SetBinContent(177,3.99e-11); h_nue_flux->SetBinContent(178,1.605e-11); h_nue_flux->SetBinContent(179,6.459e-12); h_nue_flux->SetBinContent(180,2.599e-12); h_nue_flux->SetBinContent(181,1.046e-12); h_nue_flux->SetBinContent(182,4.207e-13); h_nue_flux->SetBinContent(183,1.693e-13); h_nue_flux->SetBinContent(184,6.811e-14); h_nue_flux->SetBinContent(185,2.74e-14); h_nue_flux->SetBinContent(186,1.103e-14); h_nue_flux->SetBinContent(187,4.436e-15); h_nue_flux->SetBinContent(188,1.785e-15); h_nue_flux->SetBinContent(189,7.181e-16); h_nue_flux->SetBinContent(190,2.889e-16); h_nue_flux->SetBinContent(191,1.163e-16); h_nue_flux->SetBinContent(192,4.677e-17); h_nue_flux->SetBinContent(193,1.882e-17); h_nue_flux->SetBinContent(194,7.572e-18); h_nue_flux->SetBinContent(195,3.047e-18); h_nue_flux->SetBinContent(196,1.226e-18); h_nue_flux->SetBinContent(197,4.932e-19); h_nue_flux->SetBinContent(198,1.984e-19); h_nue_flux->SetBinContent(199,7.984e-20); h_nue_flux->SetBinContent(200,3.397e-20); //Histogram: nu cross section TH1D *h_nuXsec = new TH1D("h_nuXsec","h_nuXsec",200,0,400); h_nuXsec->SetBinContent(1,1.88482e-38); h_nuXsec->SetBinContent(2,3.33504e-38); h_nuXsec->SetBinContent(3,4.78527e-38); h_nuXsec->SetBinContent(4,6.23549e-38); h_nuXsec->SetBinContent(5,7.68571e-38); h_nuXsec->SetBinContent(6,9.13593e-38); h_nuXsec->SetBinContent(7,1.05862e-37); h_nuXsec->SetBinContent(8,1.20364e-37); h_nuXsec->SetBinContent(9,1.34866e-37); h_nuXsec->SetBinContent(10,1.49368e-37); h_nuXsec->SetBinContent(11,1.6387e-37); h_nuXsec->SetBinContent(12,1.78373e-37); h_nuXsec->SetBinContent(13,1.92875e-37); h_nuXsec->SetBinContent(14,2.07377e-37); h_nuXsec->SetBinContent(15,2.21879e-37); h_nuXsec->SetBinContent(16,2.36382e-37); h_nuXsec->SetBinContent(17,2.50884e-37); h_nuXsec->SetBinContent(18,2.65386e-37); h_nuXsec->SetBinContent(19,2.79888e-37); h_nuXsec->SetBinContent(20,2.9439e-37); h_nuXsec->SetBinContent(21,3.08893e-37); h_nuXsec->SetBinContent(22,3.23395e-37); h_nuXsec->SetBinContent(23,3.37897e-37); h_nuXsec->SetBinContent(24,3.52399e-37); h_nuXsec->SetBinContent(25,3.66902e-37); h_nuXsec->SetBinContent(26,3.81404e-37); h_nuXsec->SetBinContent(27,3.95906e-37); h_nuXsec->SetBinContent(28,4.10408e-37); h_nuXsec->SetBinContent(29,4.24911e-37); h_nuXsec->SetBinContent(30,4.39413e-37); h_nuXsec->SetBinContent(31,4.53915e-37); h_nuXsec->SetBinContent(32,4.68417e-37); h_nuXsec->SetBinContent(33,4.82919e-37); h_nuXsec->SetBinContent(34,4.97422e-37); h_nuXsec->SetBinContent(35,5.11924e-37); h_nuXsec->SetBinContent(36,5.26426e-37); h_nuXsec->SetBinContent(37,5.40928e-37); h_nuXsec->SetBinContent(38,5.55431e-37); h_nuXsec->SetBinContent(39,5.69933e-37); h_nuXsec->SetBinContent(40,5.84435e-37); h_nuXsec->SetBinContent(41,5.98937e-37); h_nuXsec->SetBinContent(42,6.13439e-37); h_nuXsec->SetBinContent(43,6.27942e-37); h_nuXsec->SetBinContent(44,6.42444e-37); h_nuXsec->SetBinContent(45,6.56946e-37); h_nuXsec->SetBinContent(46,6.71448e-37); h_nuXsec->SetBinContent(47,6.85951e-37); h_nuXsec->SetBinContent(48,7.00453e-37); h_nuXsec->SetBinContent(49,7.14955e-37); h_nuXsec->SetBinContent(50,7.29457e-37); h_nuXsec->SetBinContent(51,7.43959e-37); h_nuXsec->SetBinContent(52,7.58462e-37); h_nuXsec->SetBinContent(53,7.72964e-37); h_nuXsec->SetBinContent(54,7.87466e-37); h_nuXsec->SetBinContent(55,8.01968e-37); h_nuXsec->SetBinContent(56,8.16471e-37); h_nuXsec->SetBinContent(57,8.30973e-37); h_nuXsec->SetBinContent(58,8.45475e-37); h_nuXsec->SetBinContent(59,8.59977e-37); h_nuXsec->SetBinContent(60,8.74479e-37); h_nuXsec->SetBinContent(61,8.88982e-37); h_nuXsec->SetBinContent(62,9.03484e-37); h_nuXsec->SetBinContent(63,9.17986e-37); h_nuXsec->SetBinContent(64,9.32488e-37); h_nuXsec->SetBinContent(65,9.46991e-37); h_nuXsec->SetBinContent(66,9.61493e-37); h_nuXsec->SetBinContent(67,9.75995e-37); h_nuXsec->SetBinContent(68,9.90497e-37); h_nuXsec->SetBinContent(69,1.005e-36); h_nuXsec->SetBinContent(70,1.0195e-36); h_nuXsec->SetBinContent(71,1.034e-36); h_nuXsec->SetBinContent(72,1.04851e-36); h_nuXsec->SetBinContent(73,1.06301e-36); h_nuXsec->SetBinContent(74,1.07751e-36); h_nuXsec->SetBinContent(75,1.09201e-36); h_nuXsec->SetBinContent(76,1.10652e-36); h_nuXsec->SetBinContent(77,1.12102e-36); h_nuXsec->SetBinContent(78,1.13552e-36); h_nuXsec->SetBinContent(79,1.15002e-36); h_nuXsec->SetBinContent(80,1.16452e-36); h_nuXsec->SetBinContent(81,1.17903e-36); h_nuXsec->SetBinContent(82,1.19353e-36); h_nuXsec->SetBinContent(83,1.20803e-36); h_nuXsec->SetBinContent(84,1.22253e-36); h_nuXsec->SetBinContent(85,1.23704e-36); h_nuXsec->SetBinContent(86,1.25154e-36); h_nuXsec->SetBinContent(87,1.26604e-36); h_nuXsec->SetBinContent(88,1.28054e-36); h_nuXsec->SetBinContent(89,1.29504e-36); h_nuXsec->SetBinContent(90,1.30955e-36); h_nuXsec->SetBinContent(91,1.32405e-36); h_nuXsec->SetBinContent(92,1.33855e-36); h_nuXsec->SetBinContent(93,1.35305e-36); h_nuXsec->SetBinContent(94,1.36756e-36); h_nuXsec->SetBinContent(95,1.38206e-36); h_nuXsec->SetBinContent(96,1.39656e-36); h_nuXsec->SetBinContent(97,1.41106e-36); h_nuXsec->SetBinContent(98,1.42556e-36); h_nuXsec->SetBinContent(99,1.44007e-36); h_nuXsec->SetBinContent(100,1.45457e-36); h_nuXsec->SetBinContent(101,1.46907e-36); h_nuXsec->SetBinContent(102,1.48357e-36); h_nuXsec->SetBinContent(103,1.49808e-36); h_nuXsec->SetBinContent(104,1.51258e-36); h_nuXsec->SetBinContent(105,1.52708e-36); h_nuXsec->SetBinContent(106,1.54158e-36); h_nuXsec->SetBinContent(107,1.55608e-36); h_nuXsec->SetBinContent(108,1.57059e-36); h_nuXsec->SetBinContent(109,1.58509e-36); h_nuXsec->SetBinContent(110,1.59959e-36); h_nuXsec->SetBinContent(111,1.61409e-36); h_nuXsec->SetBinContent(112,1.6286e-36); h_nuXsec->SetBinContent(113,1.6431e-36); h_nuXsec->SetBinContent(114,1.6576e-36); h_nuXsec->SetBinContent(115,1.6721e-36); h_nuXsec->SetBinContent(116,1.6866e-36); h_nuXsec->SetBinContent(117,1.70111e-36); h_nuXsec->SetBinContent(118,1.71561e-36); h_nuXsec->SetBinContent(119,1.73011e-36); h_nuXsec->SetBinContent(120,1.74461e-36); h_nuXsec->SetBinContent(121,1.75912e-36); h_nuXsec->SetBinContent(122,1.77362e-36); h_nuXsec->SetBinContent(123,1.78812e-36); h_nuXsec->SetBinContent(124,1.80262e-36); h_nuXsec->SetBinContent(125,1.81712e-36); h_nuXsec->SetBinContent(126,1.83163e-36); h_nuXsec->SetBinContent(127,1.84613e-36); h_nuXsec->SetBinContent(128,1.86063e-36); h_nuXsec->SetBinContent(129,1.87513e-36); h_nuXsec->SetBinContent(130,1.88964e-36); h_nuXsec->SetBinContent(131,1.90414e-36); h_nuXsec->SetBinContent(132,1.91864e-36); h_nuXsec->SetBinContent(133,1.93314e-36); h_nuXsec->SetBinContent(134,1.94764e-36); h_nuXsec->SetBinContent(135,1.96215e-36); h_nuXsec->SetBinContent(136,1.97665e-36); h_nuXsec->SetBinContent(137,1.99115e-36); h_nuXsec->SetBinContent(138,2.00565e-36); h_nuXsec->SetBinContent(139,2.02016e-36); h_nuXsec->SetBinContent(140,2.03466e-36); h_nuXsec->SetBinContent(141,2.04916e-36); h_nuXsec->SetBinContent(142,2.06366e-36); h_nuXsec->SetBinContent(143,2.07816e-36); h_nuXsec->SetBinContent(144,2.09267e-36); h_nuXsec->SetBinContent(145,2.10717e-36); h_nuXsec->SetBinContent(146,2.12167e-36); h_nuXsec->SetBinContent(147,2.13617e-36); h_nuXsec->SetBinContent(148,2.15068e-36); h_nuXsec->SetBinContent(149,2.16518e-36); h_nuXsec->SetBinContent(150,2.17968e-36); h_nuXsec->SetBinContent(151,2.19418e-36); h_nuXsec->SetBinContent(152,2.20868e-36); h_nuXsec->SetBinContent(153,2.22319e-36); h_nuXsec->SetBinContent(154,2.23769e-36); h_nuXsec->SetBinContent(155,2.25219e-36); h_nuXsec->SetBinContent(156,2.26669e-36); h_nuXsec->SetBinContent(157,2.2812e-36); h_nuXsec->SetBinContent(158,2.2957e-36); h_nuXsec->SetBinContent(159,2.3102e-36); h_nuXsec->SetBinContent(160,2.3247e-36); h_nuXsec->SetBinContent(161,2.3392e-36); h_nuXsec->SetBinContent(162,2.35371e-36); h_nuXsec->SetBinContent(163,2.36821e-36); h_nuXsec->SetBinContent(164,2.38271e-36); h_nuXsec->SetBinContent(165,2.39721e-36); h_nuXsec->SetBinContent(166,2.41172e-36); h_nuXsec->SetBinContent(167,2.42622e-36); h_nuXsec->SetBinContent(168,2.44072e-36); h_nuXsec->SetBinContent(169,2.45522e-36); h_nuXsec->SetBinContent(170,2.46972e-36); h_nuXsec->SetBinContent(171,2.48423e-36); h_nuXsec->SetBinContent(172,2.49873e-36); h_nuXsec->SetBinContent(173,2.51323e-36); h_nuXsec->SetBinContent(174,2.52773e-36); h_nuXsec->SetBinContent(175,2.54224e-36); h_nuXsec->SetBinContent(176,2.55674e-36); h_nuXsec->SetBinContent(177,2.57124e-36); h_nuXsec->SetBinContent(178,2.58574e-36); h_nuXsec->SetBinContent(179,2.60024e-36); h_nuXsec->SetBinContent(180,2.61475e-36); h_nuXsec->SetBinContent(181,2.62925e-36); h_nuXsec->SetBinContent(182,2.64375e-36); h_nuXsec->SetBinContent(183,2.65825e-36); h_nuXsec->SetBinContent(184,2.67276e-36); h_nuXsec->SetBinContent(185,2.68726e-36); h_nuXsec->SetBinContent(186,2.70176e-36); h_nuXsec->SetBinContent(187,2.71626e-36); h_nuXsec->SetBinContent(188,2.73076e-36); h_nuXsec->SetBinContent(189,2.74527e-36); h_nuXsec->SetBinContent(190,2.75977e-36); h_nuXsec->SetBinContent(191,2.77427e-36); h_nuXsec->SetBinContent(192,2.78877e-36); h_nuXsec->SetBinContent(193,2.80328e-36); h_nuXsec->SetBinContent(194,2.81778e-36); h_nuXsec->SetBinContent(195,2.83228e-36); h_nuXsec->SetBinContent(196,2.84678e-36); h_nuXsec->SetBinContent(197,2.86128e-36); h_nuXsec->SetBinContent(198,2.87579e-36); h_nuXsec->SetBinContent(199,2.89029e-36); h_nuXsec->SetBinContent(200,2.90479e-36); //Histogram: location efficiency TH1D *h_eff_loc = new TH1D("h_eff_loc","h_eff_loc",200,0,400); h_eff_loc->SetBinContent(1,0.092929); h_eff_loc->SetBinContent(2,0.239795); h_eff_loc->SetBinContent(3,0.346637); h_eff_loc->SetBinContent(4,0.424362); h_eff_loc->SetBinContent(5,0.480907); h_eff_loc->SetBinContent(6,0.522041); h_eff_loc->SetBinContent(7,0.551966); h_eff_loc->SetBinContent(8,0.573736); h_eff_loc->SetBinContent(9,0.589573); h_eff_loc->SetBinContent(10,0.601094); h_eff_loc->SetBinContent(11,0.609476); h_eff_loc->SetBinContent(12,0.615573); h_eff_loc->SetBinContent(13,0.620009); h_eff_loc->SetBinContent(14,0.623236); h_eff_loc->SetBinContent(15,0.625583); h_eff_loc->SetBinContent(16,0.627291); h_eff_loc->SetBinContent(17,0.628534); h_eff_loc->SetBinContent(18,0.629437); h_eff_loc->SetBinContent(19,0.630095); h_eff_loc->SetBinContent(20,0.630573); h_eff_loc->SetBinContent(21,0.630921); h_eff_loc->SetBinContent(22,0.631174); h_eff_loc->SetBinContent(23,0.631359); h_eff_loc->SetBinContent(24,0.631493); h_eff_loc->SetBinContent(25,0.63159); h_eff_loc->SetBinContent(26,0.631661); h_eff_loc->SetBinContent(27,0.631712); h_eff_loc->SetBinContent(28,0.63175); h_eff_loc->SetBinContent(29,0.631777); h_eff_loc->SetBinContent(30,0.631797); h_eff_loc->SetBinContent(31,0.631812); h_eff_loc->SetBinContent(32,0.631822); h_eff_loc->SetBinContent(33,0.63183); h_eff_loc->SetBinContent(34,0.631835); h_eff_loc->SetBinContent(35,0.631839); h_eff_loc->SetBinContent(36,0.631842); h_eff_loc->SetBinContent(37,0.631844); h_eff_loc->SetBinContent(38,0.631846); h_eff_loc->SetBinContent(39,0.631847); h_eff_loc->SetBinContent(40,0.631848); h_eff_loc->SetBinContent(41,0.631849); h_eff_loc->SetBinContent(42,0.631849); h_eff_loc->SetBinContent(43,0.631849); h_eff_loc->SetBinContent(44,0.63185); h_eff_loc->SetBinContent(45,0.63185); h_eff_loc->SetBinContent(46,0.63185); h_eff_loc->SetBinContent(47,0.63185); h_eff_loc->SetBinContent(48,0.63185); h_eff_loc->SetBinContent(49,0.63185); h_eff_loc->SetBinContent(50,0.63185); h_eff_loc->SetBinContent(51,0.63185); h_eff_loc->SetBinContent(52,0.63185); h_eff_loc->SetBinContent(53,0.63185); h_eff_loc->SetBinContent(54,0.63185); h_eff_loc->SetBinContent(55,0.63185); h_eff_loc->SetBinContent(56,0.63185); h_eff_loc->SetBinContent(57,0.63185); h_eff_loc->SetBinContent(58,0.63185); h_eff_loc->SetBinContent(59,0.63185); h_eff_loc->SetBinContent(60,0.63185); h_eff_loc->SetBinContent(61,0.63185); h_eff_loc->SetBinContent(62,0.63185); h_eff_loc->SetBinContent(63,0.63185); h_eff_loc->SetBinContent(64,0.63185); h_eff_loc->SetBinContent(65,0.63185); h_eff_loc->SetBinContent(66,0.63185); h_eff_loc->SetBinContent(67,0.63185); h_eff_loc->SetBinContent(68,0.63185); h_eff_loc->SetBinContent(69,0.63185); h_eff_loc->SetBinContent(70,0.63185); h_eff_loc->SetBinContent(71,0.63185); h_eff_loc->SetBinContent(72,0.63185); h_eff_loc->SetBinContent(73,0.63185); h_eff_loc->SetBinContent(74,0.63185); h_eff_loc->SetBinContent(75,0.63185); h_eff_loc->SetBinContent(76,0.63185); h_eff_loc->SetBinContent(77,0.63185); h_eff_loc->SetBinContent(78,0.63185); h_eff_loc->SetBinContent(79,0.63185); h_eff_loc->SetBinContent(80,0.63185); h_eff_loc->SetBinContent(81,0.63185); h_eff_loc->SetBinContent(82,0.63185); h_eff_loc->SetBinContent(83,0.63185); h_eff_loc->SetBinContent(84,0.63185); h_eff_loc->SetBinContent(85,0.63185); h_eff_loc->SetBinContent(86,0.63185); h_eff_loc->SetBinContent(87,0.63185); h_eff_loc->SetBinContent(88,0.63185); h_eff_loc->SetBinContent(89,0.63185); h_eff_loc->SetBinContent(90,0.63185); h_eff_loc->SetBinContent(91,0.63185); h_eff_loc->SetBinContent(92,0.63185); h_eff_loc->SetBinContent(93,0.63185); h_eff_loc->SetBinContent(94,0.63185); h_eff_loc->SetBinContent(95,0.63185); h_eff_loc->SetBinContent(96,0.63185); h_eff_loc->SetBinContent(97,0.63185); h_eff_loc->SetBinContent(98,0.63185); h_eff_loc->SetBinContent(99,0.63185); h_eff_loc->SetBinContent(100,0.63185); h_eff_loc->SetBinContent(101,0.63185); h_eff_loc->SetBinContent(102,0.63185); h_eff_loc->SetBinContent(103,0.63185); h_eff_loc->SetBinContent(104,0.63185); h_eff_loc->SetBinContent(105,0.63185); h_eff_loc->SetBinContent(106,0.63185); h_eff_loc->SetBinContent(107,0.63185); h_eff_loc->SetBinContent(108,0.63185); h_eff_loc->SetBinContent(109,0.63185); h_eff_loc->SetBinContent(110,0.63185); h_eff_loc->SetBinContent(111,0.63185); h_eff_loc->SetBinContent(112,0.63185); h_eff_loc->SetBinContent(113,0.63185); h_eff_loc->SetBinContent(114,0.63185); h_eff_loc->SetBinContent(115,0.63185); h_eff_loc->SetBinContent(116,0.63185); h_eff_loc->SetBinContent(117,0.63185); h_eff_loc->SetBinContent(118,0.63185); h_eff_loc->SetBinContent(119,0.63185); h_eff_loc->SetBinContent(120,0.63185); h_eff_loc->SetBinContent(121,0.63185); h_eff_loc->SetBinContent(122,0.63185); h_eff_loc->SetBinContent(123,0.63185); h_eff_loc->SetBinContent(124,0.63185); h_eff_loc->SetBinContent(125,0.63185); h_eff_loc->SetBinContent(126,0.63185); h_eff_loc->SetBinContent(127,0.63185); h_eff_loc->SetBinContent(128,0.63185); h_eff_loc->SetBinContent(129,0.63185); h_eff_loc->SetBinContent(130,0.63185); h_eff_loc->SetBinContent(131,0.63185); h_eff_loc->SetBinContent(132,0.63185); h_eff_loc->SetBinContent(133,0.63185); h_eff_loc->SetBinContent(134,0.63185); h_eff_loc->SetBinContent(135,0.63185); h_eff_loc->SetBinContent(136,0.63185); h_eff_loc->SetBinContent(137,0.63185); h_eff_loc->SetBinContent(138,0.63185); h_eff_loc->SetBinContent(139,0.63185); h_eff_loc->SetBinContent(140,0.63185); h_eff_loc->SetBinContent(141,0.63185); h_eff_loc->SetBinContent(142,0.63185); h_eff_loc->SetBinContent(143,0.63185); h_eff_loc->SetBinContent(144,0.63185); h_eff_loc->SetBinContent(145,0.63185); h_eff_loc->SetBinContent(146,0.63185); h_eff_loc->SetBinContent(147,0.63185); h_eff_loc->SetBinContent(148,0.63185); h_eff_loc->SetBinContent(149,0.63185); h_eff_loc->SetBinContent(150,0.63185); h_eff_loc->SetBinContent(151,0.63185); h_eff_loc->SetBinContent(152,0.63185); h_eff_loc->SetBinContent(153,0.63185); h_eff_loc->SetBinContent(154,0.63185); h_eff_loc->SetBinContent(155,0.63185); h_eff_loc->SetBinContent(156,0.63185); h_eff_loc->SetBinContent(157,0.63185); h_eff_loc->SetBinContent(158,0.63185); h_eff_loc->SetBinContent(159,0.63185); h_eff_loc->SetBinContent(160,0.63185); h_eff_loc->SetBinContent(161,0.63185); h_eff_loc->SetBinContent(162,0.63185); h_eff_loc->SetBinContent(163,0.63185); h_eff_loc->SetBinContent(164,0.63185); h_eff_loc->SetBinContent(165,0.63185); h_eff_loc->SetBinContent(166,0.63185); h_eff_loc->SetBinContent(167,0.63185); h_eff_loc->SetBinContent(168,0.63185); h_eff_loc->SetBinContent(169,0.63185); h_eff_loc->SetBinContent(170,0.63185); h_eff_loc->SetBinContent(171,0.63185); h_eff_loc->SetBinContent(172,0.63185); h_eff_loc->SetBinContent(173,0.63185); h_eff_loc->SetBinContent(174,0.63185); h_eff_loc->SetBinContent(175,0.63185); h_eff_loc->SetBinContent(176,0.63185); h_eff_loc->SetBinContent(177,0.63185); h_eff_loc->SetBinContent(178,0.63185); h_eff_loc->SetBinContent(179,0.63185); h_eff_loc->SetBinContent(180,0.63185); h_eff_loc->SetBinContent(181,0.63185); h_eff_loc->SetBinContent(182,0.63185); h_eff_loc->SetBinContent(183,0.63185); h_eff_loc->SetBinContent(184,0.63185); h_eff_loc->SetBinContent(185,0.63185); h_eff_loc->SetBinContent(186,0.63185); h_eff_loc->SetBinContent(187,0.63185); h_eff_loc->SetBinContent(188,0.63185); h_eff_loc->SetBinContent(189,0.63185); h_eff_loc->SetBinContent(190,0.63185); h_eff_loc->SetBinContent(191,0.63185); h_eff_loc->SetBinContent(192,0.63185); h_eff_loc->SetBinContent(193,0.63185); h_eff_loc->SetBinContent(194,0.63185); h_eff_loc->SetBinContent(195,0.63185); h_eff_loc->SetBinContent(196,0.63185); h_eff_loc->SetBinContent(197,0.63185); h_eff_loc->SetBinContent(198,0.63185); h_eff_loc->SetBinContent(199,0.63185); h_eff_loc->SetBinContent(200,0.63185); //Histogram: identification efficiency TH1D *h_eff_id = new TH1D("h_eff_id","h_eff_id",200,0,400); h_eff_id->SetBinContent(1,0.0768534); h_eff_id->SetBinContent(2,0.213037); h_eff_id->SetBinContent(3,0.328796); h_eff_id->SetBinContent(4,0.427195); h_eff_id->SetBinContent(5,0.510836); h_eff_id->SetBinContent(6,0.581933); h_eff_id->SetBinContent(7,0.642367); h_eff_id->SetBinContent(8,0.693737); h_eff_id->SetBinContent(9,0.737403); h_eff_id->SetBinContent(10,0.77452); h_eff_id->SetBinContent(11,0.806071); h_eff_id->SetBinContent(12,0.83289); h_eff_id->SetBinContent(13,0.855687); h_eff_id->SetBinContent(14,0.875064); h_eff_id->SetBinContent(15,0.891536); h_eff_id->SetBinContent(16,0.905537); h_eff_id->SetBinContent(17,0.917438); h_eff_id->SetBinContent(18,0.927555); h_eff_id->SetBinContent(19,0.936154); h_eff_id->SetBinContent(20,0.943463); h_eff_id->SetBinContent(21,0.949677); h_eff_id->SetBinContent(22,0.954958); h_eff_id->SetBinContent(23,0.959448); h_eff_id->SetBinContent(24,0.963264); h_eff_id->SetBinContent(25,0.966507); h_eff_id->SetBinContent(26,0.969265); h_eff_id->SetBinContent(27,0.971608); h_eff_id->SetBinContent(28,0.973601); h_eff_id->SetBinContent(29,0.975294); h_eff_id->SetBinContent(30,0.976734); h_eff_id->SetBinContent(31,0.977957); h_eff_id->SetBinContent(32,0.978997); h_eff_id->SetBinContent(33,0.979881); h_eff_id->SetBinContent(34,0.980633); h_eff_id->SetBinContent(35,0.981272); h_eff_id->SetBinContent(36,0.981815); h_eff_id->SetBinContent(37,0.982276); h_eff_id->SetBinContent(38,0.982669); h_eff_id->SetBinContent(39,0.983002); h_eff_id->SetBinContent(40,0.983286); h_eff_id->SetBinContent(41,0.983526); h_eff_id->SetBinContent(42,0.983731); h_eff_id->SetBinContent(43,0.983905); h_eff_id->SetBinContent(44,0.984053); h_eff_id->SetBinContent(45,0.984179); h_eff_id->SetBinContent(46,0.984286); h_eff_id->SetBinContent(47,0.984377); h_eff_id->SetBinContent(48,0.984454); h_eff_id->SetBinContent(49,0.98452); h_eff_id->SetBinContent(50,0.984576); h_eff_id->SetBinContent(51,0.984623); h_eff_id->SetBinContent(52,0.984664); h_eff_id->SetBinContent(53,0.984698); h_eff_id->SetBinContent(54,0.984727); h_eff_id->SetBinContent(55,0.984752); h_eff_id->SetBinContent(56,0.984773); h_eff_id->SetBinContent(57,0.984791); h_eff_id->SetBinContent(58,0.984806); h_eff_id->SetBinContent(59,0.984819); h_eff_id->SetBinContent(60,0.98483); h_eff_id->SetBinContent(61,0.984839); h_eff_id->SetBinContent(62,0.984847); h_eff_id->SetBinContent(63,0.984854); h_eff_id->SetBinContent(64,0.98486); h_eff_id->SetBinContent(65,0.984865); h_eff_id->SetBinContent(66,0.984869); h_eff_id->SetBinContent(67,0.984872); h_eff_id->SetBinContent(68,0.984875); h_eff_id->SetBinContent(69,0.984878); h_eff_id->SetBinContent(70,0.98488); h_eff_id->SetBinContent(71,0.984882); h_eff_id->SetBinContent(72,0.984883); h_eff_id->SetBinContent(73,0.984885); h_eff_id->SetBinContent(74,0.984886); h_eff_id->SetBinContent(75,0.984887); h_eff_id->SetBinContent(76,0.984888); h_eff_id->SetBinContent(77,0.984888); h_eff_id->SetBinContent(78,0.984889); h_eff_id->SetBinContent(79,0.984889); h_eff_id->SetBinContent(80,0.98489); h_eff_id->SetBinContent(81,0.98489); h_eff_id->SetBinContent(82,0.98489); h_eff_id->SetBinContent(83,0.984891); h_eff_id->SetBinContent(84,0.984891); h_eff_id->SetBinContent(85,0.984891); h_eff_id->SetBinContent(86,0.984891); h_eff_id->SetBinContent(87,0.984891); h_eff_id->SetBinContent(88,0.984892); h_eff_id->SetBinContent(89,0.984892); h_eff_id->SetBinContent(90,0.984892); h_eff_id->SetBinContent(91,0.984892); h_eff_id->SetBinContent(92,0.984892); h_eff_id->SetBinContent(93,0.984892); h_eff_id->SetBinContent(94,0.984892); h_eff_id->SetBinContent(95,0.984892); h_eff_id->SetBinContent(96,0.984892); h_eff_id->SetBinContent(97,0.984892); h_eff_id->SetBinContent(98,0.984892); h_eff_id->SetBinContent(99,0.984892); h_eff_id->SetBinContent(100,0.984892); h_eff_id->SetBinContent(101,0.984892); h_eff_id->SetBinContent(102,0.984892); h_eff_id->SetBinContent(103,0.984892); h_eff_id->SetBinContent(104,0.984892); h_eff_id->SetBinContent(105,0.984892); h_eff_id->SetBinContent(106,0.984892); h_eff_id->SetBinContent(107,0.984892); h_eff_id->SetBinContent(108,0.984892); h_eff_id->SetBinContent(109,0.984892); h_eff_id->SetBinContent(110,0.984892); h_eff_id->SetBinContent(111,0.984892); h_eff_id->SetBinContent(112,0.984892); h_eff_id->SetBinContent(113,0.984892); h_eff_id->SetBinContent(114,0.984892); h_eff_id->SetBinContent(115,0.984892); h_eff_id->SetBinContent(116,0.984892); h_eff_id->SetBinContent(117,0.984892); h_eff_id->SetBinContent(118,0.984892); h_eff_id->SetBinContent(119,0.984892); h_eff_id->SetBinContent(120,0.984892); h_eff_id->SetBinContent(121,0.984892); h_eff_id->SetBinContent(122,0.984892); h_eff_id->SetBinContent(123,0.984892); h_eff_id->SetBinContent(124,0.984892); h_eff_id->SetBinContent(125,0.984892); h_eff_id->SetBinContent(126,0.984892); h_eff_id->SetBinContent(127,0.984892); h_eff_id->SetBinContent(128,0.984892); h_eff_id->SetBinContent(129,0.984892); h_eff_id->SetBinContent(130,0.984892); h_eff_id->SetBinContent(131,0.984892); h_eff_id->SetBinContent(132,0.984892); h_eff_id->SetBinContent(133,0.984892); h_eff_id->SetBinContent(134,0.984892); h_eff_id->SetBinContent(135,0.984892); h_eff_id->SetBinContent(136,0.984892); h_eff_id->SetBinContent(137,0.984892); h_eff_id->SetBinContent(138,0.984892); h_eff_id->SetBinContent(139,0.984892); h_eff_id->SetBinContent(140,0.984892); h_eff_id->SetBinContent(141,0.984892); h_eff_id->SetBinContent(142,0.984892); h_eff_id->SetBinContent(143,0.984892); h_eff_id->SetBinContent(144,0.984892); h_eff_id->SetBinContent(145,0.984892); h_eff_id->SetBinContent(146,0.984892); h_eff_id->SetBinContent(147,0.984892); h_eff_id->SetBinContent(148,0.984892); h_eff_id->SetBinContent(149,0.984892); h_eff_id->SetBinContent(150,0.984892); h_eff_id->SetBinContent(151,0.984892); h_eff_id->SetBinContent(152,0.984892); h_eff_id->SetBinContent(153,0.984892); h_eff_id->SetBinContent(154,0.984892); h_eff_id->SetBinContent(155,0.984892); h_eff_id->SetBinContent(156,0.984892); h_eff_id->SetBinContent(157,0.984892); h_eff_id->SetBinContent(158,0.984892); h_eff_id->SetBinContent(159,0.984892); h_eff_id->SetBinContent(160,0.984892); h_eff_id->SetBinContent(161,0.984892); h_eff_id->SetBinContent(162,0.984892); h_eff_id->SetBinContent(163,0.984892); h_eff_id->SetBinContent(164,0.984892); h_eff_id->SetBinContent(165,0.984892); h_eff_id->SetBinContent(166,0.984892); h_eff_id->SetBinContent(167,0.984892); h_eff_id->SetBinContent(168,0.984892); h_eff_id->SetBinContent(169,0.984892); h_eff_id->SetBinContent(170,0.984892); h_eff_id->SetBinContent(171,0.984892); h_eff_id->SetBinContent(172,0.984892); h_eff_id->SetBinContent(173,0.984892); h_eff_id->SetBinContent(174,0.984892); h_eff_id->SetBinContent(175,0.984892); h_eff_id->SetBinContent(176,0.984892); h_eff_id->SetBinContent(177,0.984892); h_eff_id->SetBinContent(178,0.984892); h_eff_id->SetBinContent(179,0.984892); h_eff_id->SetBinContent(180,0.984892); h_eff_id->SetBinContent(181,0.984892); h_eff_id->SetBinContent(182,0.984892); h_eff_id->SetBinContent(183,0.984892); h_eff_id->SetBinContent(184,0.984892); h_eff_id->SetBinContent(185,0.984892); h_eff_id->SetBinContent(186,0.984892); h_eff_id->SetBinContent(187,0.984892); h_eff_id->SetBinContent(188,0.984892); h_eff_id->SetBinContent(189,0.984892); h_eff_id->SetBinContent(190,0.984892); h_eff_id->SetBinContent(191,0.984892); h_eff_id->SetBinContent(192,0.984892); h_eff_id->SetBinContent(193,0.984892); h_eff_id->SetBinContent(194,0.984892); h_eff_id->SetBinContent(195,0.984892); h_eff_id->SetBinContent(196,0.984892); h_eff_id->SetBinContent(197,0.984892); h_eff_id->SetBinContent(198,0.984892); h_eff_id->SetBinContent(199,0.984892); h_eff_id->SetBinContent(200,0.984892); //Histogram: energy cut efficiency TH1D *h_eff_ecut = new TH1D("h_eff_ecut","Energy Cut Selection Efficiency",200,0,400); h_eff_ecut->SetBinContent(1,0.990607); h_eff_ecut->SetBinContent(2,0.986792); h_eff_ecut->SetBinContent(3,0.981455); h_eff_ecut->SetBinContent(4,0.97402); h_eff_ecut->SetBinContent(5,0.963714); h_eff_ecut->SetBinContent(6,0.949532); h_eff_ecut->SetBinContent(7,0.930208); h_eff_ecut->SetBinContent(8,0.90423); h_eff_ecut->SetBinContent(9,0.869936); h_eff_ecut->SetBinContent(10,0.825728); h_eff_ecut->SetBinContent(11,0.770459); h_eff_ecut->SetBinContent(12,0.703946); h_eff_ecut->SetBinContent(13,0.627478); h_eff_ecut->SetBinContent(14,0.544053); h_eff_ecut->SetBinContent(15,0.45808); h_eff_ecut->SetBinContent(16,0.374533); h_eff_ecut->SetBinContent(17,0.297848); h_eff_ecut->SetBinContent(18,0.231064); h_eff_ecut->SetBinContent(19,0.175512); h_eff_ecut->SetBinContent(20,0.131039); h_eff_ecut->SetBinContent(21,0.0965164); h_eff_ecut->SetBinContent(22,0.0703523); h_eff_ecut->SetBinContent(23,0.0508815); h_eff_ecut->SetBinContent(24,0.0365873); h_eff_ecut->SetBinContent(25,0.026198); h_eff_ecut->SetBinContent(26,0.0187016); h_eff_ecut->SetBinContent(27,0.0133208); h_eff_ecut->SetBinContent(28,0.0094733); h_eff_ecut->SetBinContent(29,0.00672948); h_eff_ecut->SetBinContent(30,0.00477655); h_eff_ecut->SetBinContent(31,0.00338843); h_eff_ecut->SetBinContent(32,0.00240274); h_eff_ecut->SetBinContent(33,0.0017033); h_eff_ecut->SetBinContent(34,0.00120722); h_eff_ecut->SetBinContent(35,0.000855495); h_eff_ecut->SetBinContent(36,0.000606185); h_eff_ecut->SetBinContent(37,0.000429498); h_eff_ecut->SetBinContent(38,0.000304295); h_eff_ecut->SetBinContent(39,0.000215582); h_eff_ecut->SetBinContent(40,0.000152728); h_eff_ecut->SetBinContent(41,0.000108198); h_eff_ecut->SetBinContent(42,7.66497e-05); h_eff_ecut->SetBinContent(43,5.42999e-05); h_eff_ecut->SetBinContent(44,3.84667e-05); h_eff_ecut->SetBinContent(45,2.72502e-05); h_eff_ecut->SetBinContent(46,1.93042e-05); h_eff_ecut->SetBinContent(47,1.36752e-05); h_eff_ecut->SetBinContent(48,9.68756e-06); h_eff_ecut->SetBinContent(49,6.8627e-06); h_eff_ecut->SetBinContent(50,4.86155e-06); h_eff_ecut->SetBinContent(51,3.44393e-06); h_eff_ecut->SetBinContent(52,2.43969e-06); h_eff_ecut->SetBinContent(53,1.72828e-06); h_eff_ecut->SetBinContent(54,1.22432e-06); h_eff_ecut->SetBinContent(55,8.67307e-07); h_eff_ecut->SetBinContent(56,6.14401e-07); h_eff_ecut->SetBinContent(57,4.35242e-07); h_eff_ecut->SetBinContent(58,3.08326e-07); h_eff_ecut->SetBinContent(59,2.18419e-07); h_eff_ecut->SetBinContent(60,1.54728e-07); h_eff_ecut->SetBinContent(61,1.09609e-07); h_eff_ecut->SetBinContent(62,7.76475e-08); h_eff_ecut->SetBinContent(63,5.50055e-08); h_eff_ecut->SetBinContent(64,3.8966e-08); h_eff_ecut->SetBinContent(65,2.76035e-08); h_eff_ecut->SetBinContent(66,1.95544e-08); h_eff_ecut->SetBinContent(67,1.38523e-08); h_eff_ecut->SetBinContent(68,9.81301e-09); h_eff_ecut->SetBinContent(69,6.95155e-09); h_eff_ecut->SetBinContent(70,4.92449e-09); h_eff_ecut->SetBinContent(71,3.48851e-09); h_eff_ecut->SetBinContent(72,2.47126e-09); h_eff_ecut->SetBinContent(73,1.75065e-09); h_eff_ecut->SetBinContent(74,1.24016e-09); h_eff_ecut->SetBinContent(75,8.7853e-10); h_eff_ecut->SetBinContent(76,6.22352e-10); h_eff_ecut->SetBinContent(77,4.40875e-10); h_eff_ecut->SetBinContent(78,3.12316e-10); h_eff_ecut->SetBinContent(79,2.21245e-10); h_eff_ecut->SetBinContent(80,1.5673e-10); h_eff_ecut->SetBinContent(81,1.11028e-10); h_eff_ecut->SetBinContent(82,7.86522e-11); h_eff_ecut->SetBinContent(83,5.57173e-11); h_eff_ecut->SetBinContent(84,3.94702e-11); h_eff_ecut->SetBinContent(85,2.79607e-11); h_eff_ecut->SetBinContent(86,1.98074e-11); h_eff_ecut->SetBinContent(87,1.40316e-11); h_eff_ecut->SetBinContent(88,9.94e-12); h_eff_ecut->SetBinContent(89,7.0415e-12); h_eff_ecut->SetBinContent(90,4.98821e-12); h_eff_ecut->SetBinContent(91,3.53365e-12); h_eff_ecut->SetBinContent(92,2.50324e-12); h_eff_ecut->SetBinContent(93,1.7733e-12); h_eff_ecut->SetBinContent(94,1.25621e-12); h_eff_ecut->SetBinContent(95,8.89899e-13); h_eff_ecut->SetBinContent(96,6.30405e-13); h_eff_ecut->SetBinContent(97,4.4658e-13); h_eff_ecut->SetBinContent(98,3.16357e-13); h_eff_ecut->SetBinContent(99,2.24108e-13); h_eff_ecut->SetBinContent(100,1.58758e-13); h_eff_ecut->SetBinContent(101,1.12465e-13); h_eff_ecut->SetBinContent(102,7.967e-14); h_eff_ecut->SetBinContent(103,5.64383e-14); h_eff_ecut->SetBinContent(104,3.9981e-14); h_eff_ecut->SetBinContent(105,2.83226e-14); h_eff_ecut->SetBinContent(106,2.00637e-14); h_eff_ecut->SetBinContent(107,1.42132e-14); h_eff_ecut->SetBinContent(108,1.00686e-14); h_eff_ecut->SetBinContent(109,7.13262e-15); h_eff_ecut->SetBinContent(110,5.05276e-15); h_eff_ecut->SetBinContent(111,3.57938e-15); h_eff_ecut->SetBinContent(112,2.53563e-15); h_eff_ecut->SetBinContent(113,1.79625e-15); h_eff_ecut->SetBinContent(114,1.27246e-15); h_eff_ecut->SetBinContent(115,9.01414e-16); h_eff_ecut->SetBinContent(116,6.38562e-16); h_eff_ecut->SetBinContent(117,4.52358e-16); h_eff_ecut->SetBinContent(118,3.20451e-16); h_eff_ecut->SetBinContent(119,2.27008e-16); h_eff_ecut->SetBinContent(120,1.60813e-16); h_eff_ecut->SetBinContent(121,1.1392e-16); h_eff_ecut->SetBinContent(122,8.07009e-17); h_eff_ecut->SetBinContent(123,5.71686e-17); h_eff_ecut->SetBinContent(124,4.04983e-17); h_eff_ecut->SetBinContent(125,2.8689e-17); h_eff_ecut->SetBinContent(126,2.03233e-17); h_eff_ecut->SetBinContent(127,1.43971e-17); h_eff_ecut->SetBinContent(128,1.01989e-17); h_eff_ecut->SetBinContent(129,7.22492e-18); h_eff_ecut->SetBinContent(130,5.11814e-18); h_eff_ecut->SetBinContent(131,3.62569e-18); h_eff_ecut->SetBinContent(132,2.56845e-18); h_eff_ecut->SetBinContent(133,1.81949e-18); h_eff_ecut->SetBinContent(134,1.28893e-18); h_eff_ecut->SetBinContent(135,9.13078e-19); h_eff_ecut->SetBinContent(136,6.46826e-19); h_eff_ecut->SetBinContent(137,4.58212e-19); h_eff_ecut->SetBinContent(138,3.24598e-19); h_eff_ecut->SetBinContent(139,2.29945e-19); h_eff_ecut->SetBinContent(140,1.62894e-19); h_eff_ecut->SetBinContent(141,1.15394e-19); h_eff_ecut->SetBinContent(142,8.17452e-20); h_eff_ecut->SetBinContent(143,5.79084e-20); h_eff_ecut->SetBinContent(144,4.10224e-20); h_eff_ecut->SetBinContent(145,2.90603e-20); h_eff_ecut->SetBinContent(146,2.05863e-20); h_eff_ecut->SetBinContent(147,1.45834e-20); h_eff_ecut->SetBinContent(148,1.03309e-20); h_eff_ecut->SetBinContent(149,7.31841e-21); h_eff_ecut->SetBinContent(150,5.18437e-21); h_eff_ecut->SetBinContent(151,3.67261e-21); h_eff_ecut->SetBinContent(152,2.60168e-21); h_eff_ecut->SetBinContent(153,1.84303e-21); h_eff_ecut->SetBinContent(154,1.30561e-21); h_eff_ecut->SetBinContent(155,9.24893e-22); h_eff_ecut->SetBinContent(156,6.55195e-22); h_eff_ecut->SetBinContent(157,4.64141e-22); h_eff_ecut->SetBinContent(158,3.28798e-22); h_eff_ecut->SetBinContent(159,2.32921e-22); h_eff_ecut->SetBinContent(160,1.65001e-22); h_eff_ecut->SetBinContent(161,1.16887e-22); h_eff_ecut->SetBinContent(162,8.2803e-23); h_eff_ecut->SetBinContent(163,5.86577e-23); h_eff_ecut->SetBinContent(164,4.15532e-23); h_eff_ecut->SetBinContent(165,2.94363e-23); h_eff_ecut->SetBinContent(166,2.08527e-23); h_eff_ecut->SetBinContent(167,1.47721e-23); h_eff_ecut->SetBinContent(168,1.04646e-23); h_eff_ecut->SetBinContent(169,7.41311e-24); h_eff_ecut->SetBinContent(170,5.25145e-24); h_eff_ecut->SetBinContent(171,3.72013e-24); h_eff_ecut->SetBinContent(172,2.63535e-24); h_eff_ecut->SetBinContent(173,1.86688e-24); h_eff_ecut->SetBinContent(174,1.3225e-24); h_eff_ecut->SetBinContent(175,9.36862e-25); h_eff_ecut->SetBinContent(176,6.63674e-25); h_eff_ecut->SetBinContent(177,4.70147e-25); h_eff_ecut->SetBinContent(178,3.33053e-25); h_eff_ecut->SetBinContent(179,2.35935e-25); h_eff_ecut->SetBinContent(180,1.67137e-25); h_eff_ecut->SetBinContent(181,1.184e-25); h_eff_ecut->SetBinContent(182,8.38745e-26); h_eff_ecut->SetBinContent(183,5.94168e-26); h_eff_ecut->SetBinContent(184,4.20909e-26); h_eff_ecut->SetBinContent(185,2.98172e-26); h_eff_ecut->SetBinContent(186,2.11226e-26); h_eff_ecut->SetBinContent(187,1.49632e-26); h_eff_ecut->SetBinContent(188,1.06e-26); h_eff_ecut->SetBinContent(189,7.50903e-27); h_eff_ecut->SetBinContent(190,5.31941e-27); h_eff_ecut->SetBinContent(191,3.76827e-27); h_eff_ecut->SetBinContent(192,2.66945e-27); h_eff_ecut->SetBinContent(193,1.89104e-27); h_eff_ecut->SetBinContent(194,1.33961e-27); h_eff_ecut->SetBinContent(195,9.48985e-28); h_eff_ecut->SetBinContent(196,6.72262e-28); h_eff_ecut->SetBinContent(197,4.76231e-28); h_eff_ecut->SetBinContent(198,3.37362e-28); h_eff_ecut->SetBinContent(199,2.38988e-28); h_eff_ecut->SetBinContent(200,1.69299e-28); //Expected normalized (to p.o.t.) energy distribution of selected nue events TH1D* h_nue_int = new TH1D(*h_nue_flux); h_nue_int->Multiply(h_nuXsec); h_nue_int->Multiply(h_eff_loc); h_nue_int->Multiply(h_eff_id); h_nue_int->Multiply(h_eff_ecut); h_nue_int->Scale(N_avo * pot * mass * konst); //Expected normalized (to p.o.t.) energy distribution of selected numu events TH1D* h_numu_int = new TH1D(*h_numu_flux); h_numu_int->Multiply(h_nuXsec); h_numu_int->Multiply(h_eff_loc); h_numu_int->Multiply(h_eff_id); h_numu_int->Multiply(h_eff_ecut); h_numu_int->Scale(N_avo * pot * mass * konst); //RooAbsVar from energy distribution histogram RooDataHist h_numu ("h_numu" , "h_numu" , RooArgList(Enu) , h_numu_int ); RooDataHist h_nue ("h_nue" , "h_nue" , RooArgList(Enu) , h_nue_int ); RooHistFunc f_numu ("f_numu", "f_numu", RooArgSet(Enu), h_numu ); RooHistFunc f_nue ("f_nue" , "f_nue" , RooArgSet(Enu), h_nue ); //Oscillation probability formula RooArgList varlist("VariableList"); varlist.add(theta23); varlist.add(theta13); varlist.add(theta14); varlist.add(theta24); varlist.add(delta ); varlist.add(dm2_atm); varlist.add(L ); varlist.add(Enu ); TString formula = "4 * (TMath::Power(TMath::Sin(@0 * TMath::Pi()/180.),2) * TMath::Power(TMath::Sin(@1 * TMath::Pi()/180.),2) * ( 1 - TMath::Power(TMath::Sin(@1 * TMath::Pi()/180.),2) - TMath::Power(TMath::Sin(@2 * TMath::Pi()/180.),2) - TMath::Power(TMath::Sin(@3 * TMath::Pi()/180.),2) ) + TMath::Sin(@0 * TMath::Pi()/180.) * TMath::Sin(@1 * TMath::Pi()/180.) * TMath::Sin(@2 * TMath::Pi()/180.) * TMath::Sin(@3 * TMath::Pi()/180.) * TMath::Cos(@4) ) * TMath::Power(TMath::Sin( 1.27 * @5 * @6 / @7 ),2) + 2 * TMath::Sin(@0 * TMath::Pi()/180.) * TMath::Sin(@1 * TMath::Pi()/180.) * TMath::Sin(@2 * TMath::Pi()/180.) * TMath::Sin(@3 * TMath::Pi()/180.) * TMath::Sin(@4) * TMath::Sin( 1.27 * @5 * @6 / @7 ) + 2 * TMath::Power(TMath::Sin(@2 * TMath::Pi()/180.),2) * TMath::Power(TMath::Sin(@3 * TMath::Pi()/180.),2) "; RooFormulaVar prob_osc("prob_osc","prob_osc",formula.Data(),varlist); //Expected normalized (to p.o.t.) energy distribution of selected oscillated nue events RooProduct f_osc ("f_osc", "f_osc", RooArgSet(prob_osc, f_numu)); //Integrale RooRealIntegral n_beam ("n_beam" ,"n_beam" , f_nue , RooArgSet(Enu)); RooRealIntegral n_osc ("n_osc" ,"n_osc" , f_osc , RooArgSet(Enu)); //Integrator configuration RooNumIntConfig* integral_beam_config = n_beam.getIntegratorConfig(); RooNumIntConfig* integral_osc_config = n_osc.getIntegratorConfig(); integral_beam_config-> setEpsAbs(0.1); integral_beam_config-> setEpsRel(0.01); integral_osc_config-> setEpsAbs(0.1); integral_osc_config-> setEpsRel(0.01); //Number of expected events RooFormulaVar Nexp ("Nexp","Nexp","@0 + @1",RooArgList(n_beam,n_osc)); //Poisson RooPoisson poisson ("poisson", "poisson", Nobs, Nexp); //Projection RooProjectedPdf* projpdf = (RooProjectedPdf*) poisson.createProjection(RooArgSet(theta24)); //Integrator configuration RooNumIntConfig* integral_projpdf_config = projpdf->getIntegratorConfig(); integral_projpdf_config->setEpsAbs(0.1); integral_projpdf_config->setEpsRel(0.1); //Workspace RooWorkspace w("ws","ws"); w.import(data); w.import(*projpdf); //ModelConfig ModelConfig model("model","model",&w); model.SetPdf (*projpdf); model.SetObservables (Nobs); model.SetParametersOfInterest (RooArgSet(theta14,delta)); //Feldman & Cousins FeldmanCousins fc(data,model); fc.SetConfidenceLevel(0.9); fc.UseAdaptiveSampling(true); fc.FluctuateNumDataEntries(false); fc.SetNBins(10); //Interval ConfInterval* interval = fc.GetInterval(); //Save Result TFile f("results/result.root","RECREATE"); w. Write(); model. Write(); interval-> Write(); f.Close(); }