{
//=========Macro generated from canvas: c1/Cut values
//=========  (Wed Mar 11 23:07:56 2009) by ROOT version5.20/00
   TCanvas *c1 = new TCanvas("c1", "Cut values",5,69,600,600);
   gStyle->SetOptStat(0);
   gStyle->SetOptTitle(0);
   c1->Range(-0.68,-1113.108,3.32,5070.825);
   c1->SetFillColor(10);
   c1->SetBorderMode(0);
   c1->SetBorderSize(2);
   c1->SetTickx(1);
   c1->SetTicky(1);
   c1->SetLeftMargin(0.17);
   c1->SetRightMargin(0.08);
   c1->SetTopMargin(0.08);
   c1->SetBottomMargin(0.18);
   c1->SetFrameLineWidth(2);
   c1->SetFrameBorderMode(0);
   c1->SetFrameLineWidth(2);
   c1->SetFrameBorderMode(0);

   THStack hs("hs","stacked histograms");

   
   TH1 *rhonu = new TH1F("rhonu","1-prong events",100,0,3);
   rhonu->SetBinContent(5,1775.49);
   rhonu->SetBinContent(6,3141.987);
   rhonu->SetBinContent(7,1718.422);
   rhonu->SetBinContent(8,1176.823);
   rhonu->SetBinContent(9,967.7008);
   rhonu->SetBinContent(10,821.1399);
   rhonu->SetBinContent(11,728.2875);
   rhonu->SetBinContent(12,681.8793);
   rhonu->SetBinContent(13,622.2052);
   rhonu->SetBinContent(14,637.5435);
   rhonu->SetBinContent(15,649.352);
   rhonu->SetBinContent(16,645.0507);
   rhonu->SetBinContent(17,699.2615);
   rhonu->SetBinContent(18,802.9532);
   rhonu->SetBinContent(19,865.5247);
   rhonu->SetBinContent(20,1064.349);
   rhonu->SetBinContent(21,1335.469);
   rhonu->SetBinContent(22,1803.841);
   rhonu->SetBinContent(23,2419.604);
   rhonu->SetBinContent(24,3273.185);
   rhonu->SetBinContent(25,4079.212);
   rhonu->SetBinContent(26,4358.2);
   rhonu->SetBinContent(27,4163.096);
   rhonu->SetBinContent(28,3646.158);
   rhonu->SetBinContent(29,3106.5);
   rhonu->SetBinContent(30,2711.635);
   rhonu->SetBinContent(31,2393.017);
   rhonu->SetBinContent(32,2281.276);
   rhonu->SetBinContent(33,2076.671);
   rhonu->SetBinContent(34,1971.743);
   rhonu->SetBinContent(35,1936.022);
   rhonu->SetBinContent(36,1909.006);
   rhonu->SetBinContent(37,1903.206);
   rhonu->SetBinContent(38,1822.823);
   rhonu->SetBinContent(39,1761.115);
   rhonu->SetBinContent(40,1728.243);
   rhonu->SetBinContent(41,1626.401);
   rhonu->SetBinContent(42,1596.379);
   rhonu->SetBinContent(43,1505.096);
   rhonu->SetBinContent(44,1383.235);
   rhonu->SetBinContent(45,1325.118);
   rhonu->SetBinContent(46,1262.915);
   rhonu->SetBinContent(47,1205.165);
   rhonu->SetBinContent(48,1088.281);
   rhonu->SetBinContent(49,1055.634);
   rhonu->SetBinContent(50,1049.428);
   rhonu->SetBinContent(51,984.7456);
   rhonu->SetBinContent(52,909.6859);
   rhonu->SetBinContent(53,853.8308);
   rhonu->SetBinContent(54,859.6144);
   rhonu->SetBinContent(55,760.7219);
   rhonu->SetBinContent(56,784.5514);
   rhonu->SetBinContent(57,758.0751);
   rhonu->SetBinContent(58,722.5002);
   rhonu->SetBinContent(59,736.2607);
   rhonu->SetBinContent(60,648.1995);
   rhonu->SetBinContent(61,674.9014);
   rhonu->SetBinContent(62,628.7472);
   rhonu->SetBinContent(63,581.0217);
   rhonu->SetBinContent(64,577.3681);
   rhonu->SetBinContent(65,550.2404);
   rhonu->SetBinContent(66,577.8948);
   rhonu->SetBinContent(67,537.5273);
   rhonu->SetBinContent(68,533.1882);
   rhonu->SetBinContent(69,523.4176);
   rhonu->SetBinContent(70,533.2104);
   rhonu->SetBinContent(71,487.2624);
   rhonu->SetBinContent(72,497.385);
   rhonu->SetBinContent(73,463.9187);
   rhonu->SetBinContent(74,451.4818);
   rhonu->SetBinContent(75,463.7028);
   rhonu->SetBinContent(76,450.9517);
   rhonu->SetBinContent(77,435.9056);
   rhonu->SetBinContent(78,408.5019);
   rhonu->SetBinContent(79,384.4059);
   rhonu->SetBinContent(80,391.2084);
   rhonu->SetBinContent(81,376.5051);
   rhonu->SetBinContent(82,361.7858);
   rhonu->SetBinContent(83,355.5674);
   rhonu->SetBinContent(84,345.3557);
   rhonu->SetBinContent(85,379.9715);
   rhonu->SetBinContent(86,337.6803);
   rhonu->SetBinContent(87,355.4214);
   rhonu->SetBinContent(88,334.6867);
   rhonu->SetBinContent(89,341.2133);
   rhonu->SetBinContent(90,314.5396);
   rhonu->SetBinContent(91,310.0859);
   rhonu->SetBinContent(92,295.5701);
   rhonu->SetBinContent(93,292.8593);
   rhonu->SetBinContent(94,274.4832);
   rhonu->SetBinContent(95,264.5796);
   rhonu->SetBinContent(96,249.4636);
   rhonu->SetBinContent(97,226.85);
   rhonu->SetBinContent(98,198.6084);
   rhonu->SetBinContent(99,200.1035);
   rhonu->SetBinContent(100,170.6811);
   rhonu->SetBinContent(101,2504.777);
   rhonu->SetEntries(312083);
   rhonu->SetStats(0);

   Int_t ci;   // for color index setting
   ci = TColor::GetColor("#ff0000");
   rhonu->SetFillColor(1);

   ci = TColor::GetColor("#ff0000");
   rhonu->SetLineColor(1);
   rhonu->SetLineWidth(2);
   rhonu->SetMarkerStyle(20);
   rhonu->SetMarkerSize(0.7);
   rhonu->SetFillStyle(0);
   rhonu->GetXaxis()->SetTitle("mass(#pi^{#pm}+neutrals)/GeV");
   rhonu->GetXaxis()->SetNdivisions(506);
   rhonu->GetXaxis()->SetLabelFont(42);
   rhonu->GetXaxis()->SetLabelOffset(0.015);
   rhonu->GetXaxis()->SetLabelSize(0.06);
   rhonu->GetXaxis()->SetTitleSize(0.07);
   rhonu->GetXaxis()->SetTitleFont(42);
   rhonu->GetYaxis()->SetNdivisions(506);
   rhonu->GetYaxis()->SetLabelFont(42);
   rhonu->GetYaxis()->SetLabelOffset(0.015);
   rhonu->GetYaxis()->SetLabelSize(0.06);
   rhonu->GetYaxis()->SetTitleSize(0.07);
   rhonu->GetYaxis()->SetTitleOffset(1.1);
   rhonu->GetYaxis()->SetTitleFont(42);
   rhonu->GetZaxis()->SetLabelFont(42);
   rhonu->GetZaxis()->SetLabelOffset(0.015);
   rhonu->GetZaxis()->SetLabelSize(0.06);
   rhonu->GetZaxis()->SetTitleSize(0.07);
   rhonu->GetZaxis()->SetTitleOffset(1.1);
   rhonu->GetZaxis()->SetTitleFont(42);
   rhonu->Draw("");
   
   TH1 *pinu = new TH1F("pinu","pinu decay",100,0,3);
   pinu->SetBinContent(5,1773.649);
   pinu->SetBinContent(6,3135.003);
   pinu->SetBinContent(7,1706.055);
   pinu->SetBinContent(8,1166.297);
   pinu->SetBinContent(9,951.1977);
   pinu->SetBinContent(10,807.1381);
   pinu->SetBinContent(11,696.1543);
   pinu->SetBinContent(12,633.4525);
   pinu->SetBinContent(13,551.6505);
   pinu->SetBinContent(14,527.1134);
   pinu->SetBinContent(15,487.025);
   pinu->SetBinContent(16,460.4179);
   pinu->SetBinContent(17,438.2105);
   pinu->SetBinContent(18,445.686);
   pinu->SetBinContent(19,419.5264);
   pinu->SetBinContent(20,405.6292);
   pinu->SetBinContent(21,408.521);
   pinu->SetBinContent(22,415.8756);
   pinu->SetBinContent(23,422.45);
   pinu->SetBinContent(24,437.9404);
   pinu->SetBinContent(25,451.5994);
   pinu->SetBinContent(26,447.8793);
   pinu->SetBinContent(27,457.8082);
   pinu->SetBinContent(28,457.8873);
   pinu->SetBinContent(29,507.248);
   pinu->SetBinContent(30,569.1912);
   pinu->SetBinContent(31,647.9364);
   pinu->SetBinContent(32,726.4083);
   pinu->SetBinContent(33,807.7751);
   pinu->SetBinContent(34,873.4615);
   pinu->SetBinContent(35,964.9963);
   pinu->SetBinContent(36,1046.707);
   pinu->SetBinContent(37,1096.645);
   pinu->SetBinContent(38,1094.933);
   pinu->SetBinContent(39,1088.01);
   pinu->SetBinContent(40,1056.396);
   pinu->SetBinContent(41,1030.451);
   pinu->SetBinContent(42,1019.372);
   pinu->SetBinContent(43,972.8948);
   pinu->SetBinContent(44,900.4024);
   pinu->SetBinContent(45,851.2148);
   pinu->SetBinContent(46,809.0831);
   pinu->SetBinContent(47,760.7891);
   pinu->SetBinContent(48,693.3926);
   pinu->SetBinContent(49,647.0445);
   pinu->SetBinContent(50,648.0572);
   pinu->SetBinContent(51,590.812);
   pinu->SetBinContent(52,542.3213);
   pinu->SetBinContent(53,502.5187);
   pinu->SetBinContent(54,500.4206);
   pinu->SetBinContent(55,448.3078);
   pinu->SetBinContent(56,454.5009);
   pinu->SetBinContent(57,424.7928);
   pinu->SetBinContent(58,407.2769);
   pinu->SetBinContent(59,408.6448);
   pinu->SetBinContent(60,368.379);
   pinu->SetBinContent(61,369.3439);
   pinu->SetBinContent(62,349.2696);
   pinu->SetBinContent(63,320.539);
   pinu->SetBinContent(64,319.6725);
   pinu->SetBinContent(65,296.8555);
   pinu->SetBinContent(66,332.3219);
   pinu->SetBinContent(67,309.2258);
   pinu->SetBinContent(68,293.2114);
   pinu->SetBinContent(69,301.5376);
   pinu->SetBinContent(70,299.1092);
   pinu->SetBinContent(71,273.7531);
   pinu->SetBinContent(72,286.5549);
   pinu->SetBinContent(73,254.1962);
   pinu->SetBinContent(74,267.2871);
   pinu->SetBinContent(75,267.2553);
   pinu->SetBinContent(76,248.8762);
   pinu->SetBinContent(77,254.8471);
   pinu->SetBinContent(78,244.1466);
   pinu->SetBinContent(79,222.9995);
   pinu->SetBinContent(80,223.2058);
   pinu->SetBinContent(81,216.9144);
   pinu->SetBinContent(82,215.3241);
   pinu->SetBinContent(83,204.9315);
   pinu->SetBinContent(84,199.3098);
   pinu->SetBinContent(85,218.3334);
   pinu->SetBinContent(86,193.2661);
   pinu->SetBinContent(87,206.0075);
   pinu->SetBinContent(88,206.3885);
   pinu->SetBinContent(89,194.1612);
   pinu->SetBinContent(90,181.4863);
   pinu->SetBinContent(91,182.7655);
   pinu->SetBinContent(92,169.4304);
   pinu->SetBinContent(93,168.4559);
   pinu->SetBinContent(94,167.6909);
   pinu->SetBinContent(95,142.0872);
   pinu->SetBinContent(96,136.7703);
   pinu->SetBinContent(97,118.88);
   pinu->SetBinContent(98,114.5725);
   pinu->SetBinContent(99,116.5882);
   pinu->SetBinContent(100,108.1478);
   pinu->SetBinContent(101,1617.26);
   pinu->SetEntries(154052);
   pinu->SetStats(0);

   ci = TColor::GetColor("#0000ff");
   pinu->SetFillColor(ci);
   ci = TColor::GetColor("#0000ff");
   pinu->SetLineColor(ci);
   pinu->SetLineWidth(2);
   pinu->SetMarkerStyle(20);
   pinu->SetMarkerSize(0.7);
   pinu->GetXaxis()->SetNdivisions(506);
   pinu->GetXaxis()->SetLabelFont(42);
   pinu->GetXaxis()->SetLabelOffset(0.015);
   pinu->GetXaxis()->SetLabelSize(0.06);
   pinu->GetXaxis()->SetTitleSize(0.07);
   pinu->GetXaxis()->SetTitleFont(42);
   pinu->GetYaxis()->SetNdivisions(506);
   pinu->GetYaxis()->SetLabelFont(42);
   pinu->GetYaxis()->SetLabelOffset(0.015);
   pinu->GetYaxis()->SetLabelSize(0.06);
   pinu->GetYaxis()->SetTitleSize(0.07);
   pinu->GetYaxis()->SetTitleOffset(1.1);
   pinu->GetYaxis()->SetTitleFont(42);
   pinu->GetZaxis()->SetLabelFont(42);
   pinu->GetZaxis()->SetLabelOffset(0.015);
   pinu->GetZaxis()->SetLabelSize(0.06);
   pinu->GetZaxis()->SetTitleSize(0.07);
   pinu->GetZaxis()->SetTitleOffset(1.1);
   pinu->GetZaxis()->SetTitleFont(42);
//   pinu->Draw("same");
   
   TH1 *a1nu = new TH1F("a1nu","a1nu decay",100,0,3);
   a1nu->SetBinContent(5,1735.021);
   a1nu->SetBinContent(6,3004.702);
   a1nu->SetBinContent(7,1591.105);
   a1nu->SetBinContent(8,1042.979);
   a1nu->SetBinContent(9,840.9041);
   a1nu->SetBinContent(10,697.819);
   a1nu->SetBinContent(11,573.5668);
   a1nu->SetBinContent(12,511.9727);
   a1nu->SetBinContent(13,445.9406);
   a1nu->SetBinContent(14,408.5191);
   a1nu->SetBinContent(15,364.4025);
   a1nu->SetBinContent(16,329.2153);
   a1nu->SetBinContent(17,310.5885);
   a1nu->SetBinContent(18,300.1354);
   a1nu->SetBinContent(19,288.3904);
   a1nu->SetBinContent(20,270.5697);
   a1nu->SetBinContent(21,253.8728);
   a1nu->SetBinContent(22,280.4004);
   a1nu->SetBinContent(23,278.8452);
   a1nu->SetBinContent(24,308.0201);
   a1nu->SetBinContent(25,312.4704);
   a1nu->SetBinContent(26,297.5354);
   a1nu->SetBinContent(27,310.7974);
   a1nu->SetBinContent(28,319.5931);
   a1nu->SetBinContent(29,373.0551);
   a1nu->SetBinContent(30,434.5095);
   a1nu->SetBinContent(31,491.3391);
   a1nu->SetBinContent(32,574.3915);
   a1nu->SetBinContent(33,659.028);
   a1nu->SetBinContent(34,726.419);
   a1nu->SetBinContent(35,819.5474);
   a1nu->SetBinContent(36,888.7827);
   a1nu->SetBinContent(37,960.8137);
   a1nu->SetBinContent(38,946.5638);
   a1nu->SetBinContent(39,951.6208);
   a1nu->SetBinContent(40,925.1521);
   a1nu->SetBinContent(41,892.9532);
   a1nu->SetBinContent(42,884.0264);
   a1nu->SetBinContent(43,843.2476);
   a1nu->SetBinContent(44,747.554);
   a1nu->SetBinContent(45,720.0058);
   a1nu->SetBinContent(46,688.0731);
   a1nu->SetBinContent(47,626.4185);
   a1nu->SetBinContent(48,561.7645);
   a1nu->SetBinContent(49,542.4677);
   a1nu->SetBinContent(50,529.2819);
   a1nu->SetBinContent(51,483.9149);
   a1nu->SetBinContent(52,438.3731);
   a1nu->SetBinContent(53,393.2187);
   a1nu->SetBinContent(54,390.9174);
   a1nu->SetBinContent(55,350.4067);
   a1nu->SetBinContent(56,338.6364);
   a1nu->SetBinContent(57,312.858);
   a1nu->SetBinContent(58,301.9796);
   a1nu->SetBinContent(59,305.8426);
   a1nu->SetBinContent(60,273.6712);
   a1nu->SetBinContent(61,277.9406);
   a1nu->SetBinContent(62,252.0318);
   a1nu->SetBinContent(63,239.5568);
   a1nu->SetBinContent(64,245.042);
   a1nu->SetBinContent(65,213.1244);
   a1nu->SetBinContent(66,244.3849);
   a1nu->SetBinContent(67,220.4951);
   a1nu->SetBinContent(68,207.8519);
   a1nu->SetBinContent(69,223.7519);
   a1nu->SetBinContent(70,219.0666);
   a1nu->SetBinContent(71,201.1034);
   a1nu->SetBinContent(72,202.3794);
   a1nu->SetBinContent(73,164.9831);
   a1nu->SetBinContent(74,179.95);
   a1nu->SetBinContent(75,197.6624);
   a1nu->SetBinContent(76,175.2932);
   a1nu->SetBinContent(77,177.966);
   a1nu->SetBinContent(78,169.732);
   a1nu->SetBinContent(79,165.9482);
   a1nu->SetBinContent(80,156.7427);
   a1nu->SetBinContent(81,159.9012);
   a1nu->SetBinContent(82,156.0856);
   a1nu->SetBinContent(83,141.6014);
   a1nu->SetBinContent(84,130.5899);
   a1nu->SetBinContent(85,168.6177);
   a1nu->SetBinContent(86,144.3123);
   a1nu->SetBinContent(87,145.4169);
   a1nu->SetBinContent(88,148.1944);
   a1nu->SetBinContent(89,138.3669);
   a1nu->SetBinContent(90,127.5997);
   a1nu->SetBinContent(91,138.643);
   a1nu->SetBinContent(92,125.2698);
   a1nu->SetBinContent(93,116.5501);
   a1nu->SetBinContent(94,124.7143);
   a1nu->SetBinContent(95,105.7797);
   a1nu->SetBinContent(96,109.4619);
   a1nu->SetBinContent(97,93.03183);
   a1nu->SetBinContent(98,85.70242);
   a1nu->SetBinContent(99,96.95525);
   a1nu->SetBinContent(100,87.78473);
   a1nu->SetBinContent(101,1247.246);
   a1nu->SetEntries(123403);
   a1nu->SetStats(0);

   ci = TColor::GetColor("#ff00ff");
   a1nu->SetFillColor(ci);

   ci = TColor::GetColor("#ff00ff");
   a1nu->SetLineColor(ci);
   a1nu->SetLineWidth(2);
   a1nu->SetMarkerStyle(20);
   a1nu->SetMarkerSize(0.7);
   a1nu->GetXaxis()->SetNdivisions(506);
   a1nu->GetXaxis()->SetLabelFont(42);
   a1nu->GetXaxis()->SetLabelOffset(0.015);
   a1nu->GetXaxis()->SetLabelSize(0.06);
   a1nu->GetXaxis()->SetTitleSize(0.07);
   a1nu->GetXaxis()->SetTitleFont(42);
   a1nu->GetYaxis()->SetNdivisions(506);
   a1nu->GetYaxis()->SetLabelFont(42);
   a1nu->GetYaxis()->SetLabelOffset(0.015);
   a1nu->GetYaxis()->SetLabelSize(0.06);
   a1nu->GetYaxis()->SetTitleSize(0.07);
   a1nu->GetYaxis()->SetTitleOffset(1.1);
   a1nu->GetYaxis()->SetTitleFont(42);
   a1nu->GetZaxis()->SetLabelFont(42);
   a1nu->GetZaxis()->SetLabelOffset(0.015);
   a1nu->GetZaxis()->SetLabelSize(0.06);
   a1nu->GetZaxis()->SetTitleSize(0.07);
   a1nu->GetZaxis()->SetTitleOffset(1.1);
   a1nu->GetZaxis()->SetTitleFont(42);
//   a1nu->Draw("same");

   TH1 *piex = new TH1F("piex","piex decay",100,0,3);
   TH1 *piX = new TH1F("piX","piex decay",100,0,3);
   piex->Add(pinu,a1nu,1.0,-1.0);

    
   TH1 *enunu = new TH1F("enunu","enunu decay",100,0,3);
   enunu->SetBinContent(5,1735.021);
   enunu->SetBinContent(6,3004.702);
   enunu->SetBinContent(7,1591.07);
   enunu->SetBinContent(8,1042.979);
   enunu->SetBinContent(9,840.2788);
   enunu->SetBinContent(10,697.1588);
   enunu->SetBinContent(11,572.9415);
   enunu->SetBinContent(12,511.3474);
   enunu->SetBinContent(13,445.3153);
   enunu->SetBinContent(14,407.8239);
   enunu->SetBinContent(15,363.1518);
   enunu->SetBinContent(16,327.9297);
   enunu->SetBinContent(17,306.3858);
   enunu->SetBinContent(18,298.1896);
   enunu->SetBinContent(19,283.8734);
   enunu->SetBinContent(20,263.3102);
   enunu->SetBinContent(21,242.4771);
   enunu->SetBinContent(22,267.2684);
   enunu->SetBinContent(23,262.6913);
   enunu->SetBinContent(24,267.9353);
   enunu->SetBinContent(25,266.0624);
   enunu->SetBinContent(26,243.6171);
   enunu->SetBinContent(27,219.0226);
   enunu->SetBinContent(28,194.501);
   enunu->SetBinContent(29,194.5328);
   enunu->SetBinContent(30,188.7683);
   enunu->SetBinContent(31,151.0482);
   enunu->SetBinContent(32,168.6147);
   enunu->SetBinContent(33,161.8852);
   enunu->SetBinContent(34,167.3323);
   enunu->SetBinContent(35,154.0986);
   enunu->SetBinContent(36,172.1637);
   enunu->SetBinContent(37,174.5254);
   enunu->SetBinContent(38,163.5486);
   enunu->SetBinContent(39,164.7708);
   enunu->SetBinContent(40,170.672);
   enunu->SetBinContent(41,166.4691);
   enunu->SetBinContent(42,203.1484);
   enunu->SetBinContent(43,190.9908);
   enunu->SetBinContent(44,197.3838);
   enunu->SetBinContent(45,199.4345);
   enunu->SetBinContent(46,192.0384);
   enunu->SetBinContent(47,200.8279);
   enunu->SetBinContent(48,191.7623);
   enunu->SetBinContent(49,191.0354);
   enunu->SetBinContent(50,193.6001);
   enunu->SetBinContent(51,189.4324);
   enunu->SetBinContent(52,166.263);
   enunu->SetBinContent(53,175.186);
   enunu->SetBinContent(54,171.1577);
   enunu->SetBinContent(55,139.4082);
   enunu->SetBinContent(56,149.9691);
   enunu->SetBinContent(57,124.1272);
   enunu->SetBinContent(58,131.5168);
   enunu->SetBinContent(59,135.3418);
   enunu->SetBinContent(60,107.4144);
   enunu->SetBinContent(61,123.2161);
   enunu->SetBinContent(62,104.7417);
   enunu->SetBinContent(63,100.9833);
   enunu->SetBinContent(64,101.5071);
   enunu->SetBinContent(65,97.61543);
   enunu->SetBinContent(66,99.28506);
   enunu->SetBinContent(67,91.57796);
   enunu->SetBinContent(68,93.93961);
   enunu->SetBinContent(69,103.764);
   enunu->SetBinContent(70,94.73635);
   enunu->SetBinContent(71,87.85453);
   enunu->SetBinContent(72,88.13068);
   enunu->SetBinContent(73,76.15408);
   enunu->SetBinContent(74,69.45009);
   enunu->SetBinContent(75,86.63878);
   enunu->SetBinContent(76,75.4907);
   enunu->SetBinContent(77,80.87112);
   enunu->SetBinContent(78,79.69031);
   enunu->SetBinContent(79,83.61687);
   enunu->SetBinContent(80,71.59908);
   enunu->SetBinContent(81,73.05606);
   enunu->SetBinContent(82,80.6997);
   enunu->SetBinContent(83,68.85332);
   enunu->SetBinContent(84,60.55259);
   enunu->SetBinContent(85,77.53814);
   enunu->SetBinContent(86,69.65639);
   enunu->SetBinContent(87,64.44106);
   enunu->SetBinContent(88,72.22122);
   enunu->SetBinContent(89,75.76373);
   enunu->SetBinContent(90,62.98408);
   enunu->SetBinContent(91,72.95129);
   enunu->SetBinContent(92,56.94025);
   enunu->SetBinContent(93,56.10543);
   enunu->SetBinContent(94,61.8699);
   enunu->SetBinContent(95,55.61658);
   enunu->SetBinContent(96,52.97558);
   enunu->SetBinContent(97,38.45643);
   enunu->SetBinContent(98,35.57418);
   enunu->SetBinContent(99,47.06823);
   enunu->SetBinContent(100,37.55175);
   enunu->SetBinContent(101,609.3654);
   enunu->SetEntries(69137);
   enunu->SetStats(0);

   ci = TColor::GetColor("#9c9a00");

   TH1 *a1ex = new TH1F("a1ex","a1ex decay",100,0,3);
   TH1 *a1X = new TH1F("a1X","a1iX decay",100,0,3);
   a1ex->Add(a1nu,enunu,1.0,-1.0);



   piX->Add(enunu,piex);
   a1X->Add(piX,a1ex);



   ci = TColor::GetColor("#ff00ff");
   a1X->SetFillColor(2);

   ci = TColor::GetColor("#ff00ff");
   a1X->SetLineColor(2);
   a1X->SetLineWidth(2);
   a1X->SetMarkerStyle(20);
   a1X->SetMarkerSize(0.7);
   a1X->SetFillStyle(0);
   a1X->Draw("same");

   ci = TColor::GetColor("#0000ff");
   piX->SetFillColor(ci); 
   piX->SetFillStyle(0);
   ci = TColor::GetColor("#0000ff");
   piX->SetLineColor(ci);
   piX->SetLineWidth(2);
   piX->SetMarkerStyle(20);
   piX->SetMarkerSize(0.7);

   piX->Draw("same");


   ci = TColor::GetColor("#9c9a00");
   enunu->SetFillColor(6);
   enunu->SetLineColor(6);
   enunu->SetLineWidth(2);
   enunu->SetFillStyle(0);
   enunu->SetMarkerStyle(20);
   enunu->SetMarkerSize(0.7);
   enunu->GetXaxis()->SetNdivisions(506);
   enunu->GetXaxis()->SetLabelFont(42);
   enunu->GetXaxis()->SetLabelOffset(0.015);
   enunu->GetXaxis()->SetLabelSize(0.06);
   enunu->GetXaxis()->SetTitleSize(0.07);
   enunu->GetXaxis()->SetTitleFont(42);
   enunu->GetYaxis()->SetNdivisions(506);
   enunu->GetYaxis()->SetLabelFont(42);
   enunu->GetYaxis()->SetLabelOffset(0.015);
   enunu->GetYaxis()->SetLabelSize(0.06);
   enunu->GetYaxis()->SetTitleSize(0.07);
   enunu->GetYaxis()->SetTitleOffset(1.1);
   enunu->GetYaxis()->SetTitleFont(42);
   enunu->GetZaxis()->SetLabelFont(42);
   enunu->GetZaxis()->SetLabelOffset(0.015);
   enunu->GetZaxis()->SetLabelSize(0.06);
   enunu->GetZaxis()->SetTitleSize(0.07);
   enunu->GetZaxis()->SetTitleOffset(1.1);
   enunu->GetZaxis()->SetTitleFont(42);
   enunu->Draw("same");
   
   TH1 *mununu = new TH1F("mununu","mununu decay",100,0,3);
   mununu->SetBinContent(5,105.7416);
   mununu->SetBinContent(6,277.5852);
   mununu->SetBinContent(7,161.5296);
   mununu->SetBinContent(8,112.9694);
   mununu->SetBinContent(9,109.2809);
   mununu->SetBinContent(10,81.80757);
   mununu->SetBinContent(11,71.94191);
   mununu->SetBinContent(12,71.7705);
   mununu->SetBinContent(13,60.61926);
   mununu->SetBinContent(14,63.05075);
   mununu->SetBinContent(15,75.83673);
   mununu->SetBinContent(16,73.47504);
   mununu->SetBinContent(17,82.33128);
   mununu->SetBinContent(18,106.8526);
   mununu->SetBinContent(19,110.7728);
   mununu->SetBinContent(20,116.6485);
   mununu->SetBinContent(21,125.2349);
   mununu->SetBinContent(22,146.6933);
   mununu->SetBinContent(23,163.1649);
   mununu->SetBinContent(24,173.5862);
   mununu->SetBinContent(25,176.6461);
   mununu->SetBinContent(26,166.1202);
   mununu->SetBinContent(27,157.813);
   mununu->SetBinContent(28,134.3706);
   mununu->SetBinContent(29,141.2462);
   mununu->SetBinContent(30,138.6432);
   mununu->SetBinContent(31,113.3218);
   mununu->SetBinContent(32,129.0155);
   mununu->SetBinContent(33,126.5205);
   mununu->SetBinContent(34,131.4852);
   mununu->SetBinContent(35,114.5344);
   mununu->SetBinContent(36,143.227);
   mununu->SetBinContent(37,146.214);
   mununu->SetBinContent(38,140.1351);
   mununu->SetBinContent(39,143.719);
   mununu->SetBinContent(40,150.3851);
   mununu->SetBinContent(41,146.8076);
   mununu->SetBinContent(42,179.9793);
   mununu->SetBinContent(43,175.1829);
   mununu->SetBinContent(44,179.738);
   mununu->SetBinContent(45,179.1477);
   mununu->SetBinContent(46,176.5796);
   mununu->SetBinContent(47,183.7375);
   mununu->SetBinContent(48,180.192);
   mununu->SetBinContent(49,180.3698);
   mununu->SetBinContent(50,176.9605);
   mununu->SetBinContent(51,176.9606);
   mununu->SetBinContent(52,154.4165);
   mununu->SetBinContent(53,166.571);
   mununu->SetBinContent(54,160.6318);
   mununu->SetBinContent(55,133.5008);
   mununu->SetBinContent(56,140.4843);
   mununu->SetBinContent(57,114.7472);
   mununu->SetBinContent(58,120.6449);
   mununu->SetBinContent(59,124.8159);
   mununu->SetBinContent(60,99.21526);
   mununu->SetBinContent(61,107.9668);
   mununu->SetBinContent(62,93.5206);
   mununu->SetBinContent(63,87.6101);
   mununu->SetBinContent(64,91.12088);
   mununu->SetBinContent(65,90.21619);
   mununu->SetBinContent(66,90.49548);
   mununu->SetBinContent(67,84.59454);
   mununu->SetBinContent(68,84.48978);
   mununu->SetBinContent(69,92.09219);
   mununu->SetBinContent(70,84.83577);
   mununu->SetBinContent(71,81.46152);
   mununu->SetBinContent(72,78.9951);
   mununu->SetBinContent(73,72.26241);
   mununu->SetBinContent(74,60.00025);
   mununu->SetBinContent(75,83.37244);
   mununu->SetBinContent(76,70.90386);
   mununu->SetBinContent(77,75.24312);
   mununu->SetBinContent(78,72.77673);
   mununu->SetBinContent(79,71.87521);
   mununu->SetBinContent(80,63.60941);
   mununu->SetBinContent(81,65.10131);
   mununu->SetBinContent(82,74.37653);
   mununu->SetBinContent(83,60.30816);
   mununu->SetBinContent(84,53.18826);
   mununu->SetBinContent(85,63.46973);
   mununu->SetBinContent(86,59.96531);
   mununu->SetBinContent(87,57.14656);
   mununu->SetBinContent(88,60.93348);
   mununu->SetBinContent(89,68.98981);
   mununu->SetBinContent(90,58.88291);
   mununu->SetBinContent(91,63.39671);
   mununu->SetBinContent(92,50.72184);
   mununu->SetBinContent(93,48.67126);
   mununu->SetBinContent(94,54.40082);
   mununu->SetBinContent(95,49.01725);
   mununu->SetBinContent(96,49.32833);
   mununu->SetBinContent(97,37.17085);
   mununu->SetBinContent(98,34.87901);
   mununu->SetBinContent(99,45.78265);
   mununu->SetBinContent(100,35.19009);
   mununu->SetBinContent(101,592.7607);
   mununu->SetEntries(32104);
   mununu->SetStats(0);

   ci = TColor::GetColor("#009a00");
   mununu->SetFillColor(ci);

   ci = TColor::GetColor("#009a00");
   mununu->SetLineColor(ci);
   mununu->SetLineWidth(2);
   mununu->SetFillStyle(0);
  
   mununu->SetMarkerStyle(20);
   mununu->SetMarkerSize(0.7);
   mununu->GetXaxis()->SetNdivisions(506);
   mununu->GetXaxis()->SetLabelFont(42);
   mununu->GetXaxis()->SetLabelOffset(0.015);
   mununu->GetXaxis()->SetLabelSize(0.06);
   mununu->GetXaxis()->SetTitleSize(0.07);
   mununu->GetXaxis()->SetTitleFont(42);
   mununu->GetYaxis()->SetNdivisions(506);
   mununu->GetYaxis()->SetLabelFont(42);
   mununu->GetYaxis()->SetLabelOffset(0.015);
   mununu->GetYaxis()->SetLabelSize(0.06);
   mununu->GetYaxis()->SetTitleSize(0.07);
   mununu->GetYaxis()->SetTitleOffset(1.1);
   mununu->GetYaxis()->SetTitleFont(42);
   mununu->GetZaxis()->SetLabelFont(42);
   mununu->GetZaxis()->SetLabelOffset(0.015);
   mununu->GetZaxis()->SetLabelSize(0.06);
   mununu->GetZaxis()->SetTitleSize(0.07);
   mununu->GetZaxis()->SetTitleOffset(1.1);
   mununu->GetZaxis()->SetTitleFont(42);
   mununu->Draw("same");
   
   TH1 *other = new TH1F("other","other tau decay",100,0,3);
   other->SetBinContent(5,43.90666);
   other->SetBinContent(6,100.0374);
   other->SetBinContent(7,55.85781);
   other->SetBinContent(8,47.72847);
   other->SetBinContent(9,52.24231);
   other->SetBinContent(10,46.99839);
   other->SetBinContent(11,39.77374);
   other->SetBinContent(12,45.1605);
   other->SetBinContent(13,42.90359);
   other->SetBinContent(14,47.69676);
   other->SetBinContent(15,60.86364);
   other->SetBinContent(16,67.39312);
   other->SetBinContent(17,70.72928);
   other->SetBinContent(18,93.23498);
   other->SetBinContent(19,100.977);
   other->SetBinContent(20,110.0143);
   other->SetBinContent(21,117.3151);
   other->SetBinContent(22,139.3289);
   other->SetBinContent(23,159.9335);
   other->SetBinContent(24,169.6945);
   other->SetBinContent(25,170.323);
   other->SetBinContent(26,161.673);
   other->SetBinContent(27,149.163);
   other->SetBinContent(28,129.4029);
   other->SetBinContent(29,136.2435);
   other->SetBinContent(30,133.7104);
   other->SetBinContent(31,107.6938);
   other->SetBinContent(32,122.867);
   other->SetBinContent(33,120.9624);
   other->SetBinContent(34,125.2668);
   other->SetBinContent(35,109.5317);
   other->SetBinContent(36,138.71);
   other->SetBinContent(37,139.9606);
   other->SetBinContent(38,133.2914);
   other->SetBinContent(39,137.5355);
   other->SetBinContent(40,146.0078);
   other->SetBinContent(41,142.5001);
   other->SetBinContent(42,173.9704);
   other->SetBinContent(43,173.2021);
   other->SetBinContent(44,172.234);
   other->SetBinContent(45,173.4498);
   other->SetBinContent(46,171.0215);
   other->SetBinContent(47,180.4014);
   other->SetBinContent(48,172.342);
   other->SetBinContent(49,175.437);
   other->SetBinContent(50,173.6593);
   other->SetBinContent(51,173.6943);
   other->SetBinContent(52,150.7344);
   other->SetBinContent(53,161.7429);
   other->SetBinContent(54,155.9751);
   other->SetBinContent(55,129.0886);
   other->SetBinContent(56,133.7802);
   other->SetBinContent(57,111.0301);
   other->SetBinContent(58,116.3025);
   other->SetBinContent(59,117.4865);
   other->SetBinContent(60,94.14275);
   other->SetBinContent(61,102.9991);
   other->SetBinContent(62,85.77536);
   other->SetBinContent(63,80.83618);
   other->SetBinContent(64,86.11821);
   other->SetBinContent(65,85.31827);
   other->SetBinContent(66,85.94357);
   other->SetBinContent(67,81.84878);
   other->SetBinContent(68,76.63979);
   other->SetBinContent(69,84.83261);
   other->SetBinContent(70,78.09677);
   other->SetBinContent(71,76.49377);
   other->SetBinContent(72,73.61153);
   other->SetBinContent(73,68.44057);
   other->SetBinContent(74,52.56609);
   other->SetBinContent(75,75.41769);
   other->SetBinContent(76,66.56145);
   other->SetBinContent(77,69.12946);
   other->SetBinContent(78,66.73289);
   other->SetBinContent(79,68.53903);
   other->SetBinContent(80,59.09241);
   other->SetBinContent(81,58.29248);
   other->SetBinContent(82,71.80537);
   other->SetBinContent(83,55.72132);
   other->SetBinContent(84,47.56026);
   other->SetBinContent(85,57.18148);
   other->SetBinContent(86,56.62914);
   other->SetBinContent(87,53.21997);
   other->SetBinContent(88,52.24865);
   other->SetBinContent(89,59.78439);
   other->SetBinContent(90,52.21375);
   other->SetBinContent(91,58.46388);
   other->SetBinContent(92,46.34451);
   other->SetBinContent(93,42.31318);
   other->SetBinContent(94,47.4174);
   other->SetBinContent(95,45.43983);
   other->SetBinContent(96,48.1475);
   other->SetBinContent(97,37.17085);
   other->SetBinContent(98,34.87901);
   other->SetBinContent(99,45.19223);
   other->SetBinContent(100,34.59967);
   other->SetBinContent(101,588.4183);
   other->SetEntries(28949);
   other->SetStats(0);

   ci = TColor::GetColor("#00009c");
   other->SetFillColor(8);

   ci = TColor::GetColor("#00009c");
   other->SetLineColor(8);
   other->SetLineWidth(2);
   other->SetFillStyle(3013);
   other->SetMarkerStyle(20);
   other->SetMarkerSize(0.7);
   other->GetXaxis()->SetNdivisions(506);
   other->GetXaxis()->SetLabelFont(42);
   other->GetXaxis()->SetLabelOffset(0.015);
   other->GetXaxis()->SetLabelSize(0.06);
   other->GetXaxis()->SetTitleSize(0.07);
   other->GetXaxis()->SetTitleFont(42);
   other->GetYaxis()->SetNdivisions(506);
   other->GetYaxis()->SetLabelFont(42);
   other->GetYaxis()->SetLabelOffset(0.015);
   other->GetYaxis()->SetLabelSize(0.06);
   other->GetYaxis()->SetTitleSize(0.07);
   other->GetYaxis()->SetTitleOffset(1.1);
   other->GetYaxis()->SetTitleFont(42);
   other->GetZaxis()->SetLabelFont(42);
   other->GetZaxis()->SetLabelOffset(0.015);
   other->GetZaxis()->SetLabelSize(0.06);
   other->GetZaxis()->SetTitleSize(0.07);
   other->GetZaxis()->SetTitleOffset(1.1);
   other->GetZaxis()->SetTitleFont(42);
   other->Draw("same");
   
   TLegend *leg = new TLegend(0.6,0.6,0.9,0.9,NULL,"brNDC");
   leg->SetBorderSize(0);
   leg->SetLineColor(1);
   leg->SetLineStyle(1);
   leg->SetLineWidth(2);
   leg->SetFillColor(10);
   leg->SetFillStyle(1001);


   entry=leg->AddEntry("rhonu","#rho^{#pm}#nu","f");

   ci = TColor::GetColor("#ff0000");
   entry->SetFillColor(ci);
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#ff0000");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);


   entry=leg->AddEntry("a1X","a_{1}#nu","f");

   ci = TColor::GetColor("#ff00ff");
   entry->SetFillColor(ci);
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#ff00ff");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);

   TLegendEntry *entry=leg->AddEntry("piX","#pi^{#pm}#nu","f");
   ci = TColor::GetColor("#0000ff");
   entry->SetFillColor(ci);
   entry->SetFillStyle(1001);
   ci = TColor::GetColor("#0000ff");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);

   entry=leg->AddEntry("enunu","e#nu#nu","f");

   ci = TColor::GetColor("#9c9a00");
   entry->SetFillColor(ci);
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#9c9a00");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("mununu","#mu#nu#nu","f");

   ci = TColor::GetColor("#009a00");
   entry->SetFillColor(ci);
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#009a00");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("other","Other","f");

   ci = TColor::GetColor("#00009c");
   entry->SetFillColor(ci);
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#00009c");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   leg->Draw();
   

   tex = new TLatex(0.2,4100,"a)");
   tex->SetTextFont(42);
   tex->SetLineWidth(2);
   tex->Draw();

   c1->Modified();
   c1->cd();
   c1->SetSelected(c1);
   c1->SaveAs("090311-tau-invmass-1p-loi.eps");
   c1->SaveAs("090311-tau-invmass-1p-loi.gif");
}