{ //=========Macro generated from canvas: c1/c1 //========= (Tue Mar 3 18:38:43 2009) by ROOT version5.14/00 TCanvas *c1 = new TCanvas("c1", "c1",14,34,600,600); gStyle->SetOptStat(0); c1->Range(-10.2,-267.568,49.8,1218.92); c1->SetFillColor(10); c1->SetBorderMode(0); c1->SetBorderSize(2); c1->SetTickx(); c1->SetTicky(); 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); TH1 *a = new TH1F("","",45,0,45); a->SetMinimum(0); a->SetMaximum(1100); a->SetEntries(40732); a->SetDirectory(0); a->SetStats(0); a->SetLineStyle(9); a->SetMarkerStyle(20); a->SetMarkerSize(0.7); a->GetXaxis()->SetTitle("E_{#pi}/GeV"); a->GetXaxis()->SetNdivisions(506); a->GetXaxis()->SetLabelFont(42); a->GetXaxis()->SetLabelOffset(0.015); a->GetXaxis()->SetLabelSize(0.06); a->GetXaxis()->SetTitleSize(0.07); a->GetXaxis()->SetTitleFont(42); a->GetYaxis()->SetTitle(""); a->GetYaxis()->SetNdivisions(506); a->GetYaxis()->SetLabelFont(42); a->GetYaxis()->SetLabelOffset(0.015); a->GetYaxis()->SetLabelSize(0.06); a->GetYaxis()->SetTitleSize(0.07); a->GetYaxis()->SetTitleOffset(1.1); a->GetYaxis()->SetTitleFont(42); a->GetZaxis()->SetLabelFont(42); a->GetZaxis()->SetLabelOffset(0.015); a->GetZaxis()->SetLabelSize(0.06); a->GetZaxis()->SetTitleSize(0.07); a->GetZaxis()->SetTitleOffset(1.1); a->GetZaxis()->SetTitleFont(42); a->Draw(""); TGraphErrors *gre = new TGraphErrors(45); gre->SetName("Graph"); gre->SetTitle("Graph"); gre->SetFillColor(1); gre->SetLineWidth(2); gre->SetMarkerStyle(24); gre->SetMarkerSize(0.7); gre->SetPoint(0,0.5,-1); gre->SetPointError(0,0,-1); gre->SetPoint(1,1.5,-1); gre->SetPointError(1,0,-1); gre->SetPoint(2,2.5,835.383); gre->SetPointError(2,0,169.785); gre->SetPoint(3,3.5,811.232); gre->SetPointError(3,0,83.5794); gre->SetPoint(4,4.5,820.985); gre->SetPointError(4,0,60.2182); gre->SetPoint(5,5.5,909.631); gre->SetPointError(5,0,57.5381); gre->SetPoint(6,6.5,895.972); gre->SetPointError(6,0,54.0836); gre->SetPoint(7,7.5,805.348); gre->SetPointError(7,0,48.992); gre->SetPoint(8,8.5,713.819); gre->SetPointError(8,0,44.3519); gre->SetPoint(9,9.5,747.495); gre->SetPointError(9,0,45.323); gre->SetPoint(10,10.5,621.392); gre->SetPointError(10,0,40.4441); gre->SetPoint(11,11.5,642.512); gre->SetPointError(11,0,41.1759); gre->SetPoint(12,12.5,574.956); gre->SetPointError(12,0,38.0994); gre->SetPoint(13,13.5,711.834); gre->SetPointError(13,0,43.3758); gre->SetPoint(14,14.5,632.053); gre->SetPointError(14,0,39.9242); gre->SetPoint(15,15.5,571.546); gre->SetPointError(15,0,37.3128); gre->SetPoint(16,16.5,532.478); gre->SetPointError(16,0,35.6571); gre->SetPoint(17,17.5,487.067); gre->SetPointError(17,0,33.7854); gre->SetPoint(18,18.5,497.454); gre->SetPointError(18,0,34.3219); gre->SetPoint(19,19.5,461.82); gre->SetPointError(19,0,32.9216); gre->SetPoint(20,20.5,442.087); gre->SetPointError(20,0,32.269); gre->SetPoint(21,21.5,404.592); gre->SetPointError(21,0,30.8758); gre->SetPoint(22,22.5,406.751); gre->SetPointError(22,0,31.3749); gre->SetPoint(23,23.5,318.151); gre->SetPointError(23,0,27.3149); gre->SetPoint(24,24.5,317.223); gre->SetPointError(24,0,27.544); gre->SetPoint(25,25.5,344.974); gre->SetPointError(25,0,29.2281); gre->SetPoint(26,26.5,270.047); gre->SetPointError(26,0,25.4218); gre->SetPoint(27,27.5,238.374); gre->SetPointError(27,0,23.8486); gre->SetPoint(28,28.5,257.966); gre->SetPointError(28,0,25.4101); gre->SetPoint(29,29.5,210.854); gre->SetPointError(29,0,22.9047); gre->SetPoint(30,30.5,186.016); gre->SetPointError(30,0,21.7022); gre->SetPoint(31,31.5,179.516); gre->SetPointError(31,0,21.7412); gre->SetPoint(32,32.5,213.506); gre->SetPointError(32,0,24.8319); gre->SetPoint(33,33.5,204.497); gre->SetPointError(33,0,24.8603); gre->SetPoint(34,34.5,98.5616); gre->SetPointError(34,0,16.074); gre->SetPoint(35,35.5,109.477); gre->SetPointError(35,0,17.4417); gre->SetPoint(36,36.5,120.173); gre->SetPointError(36,0,18.9229); gre->SetPoint(37,37.5,50.1321); gre->SetPointError(37,0,11.124); gre->SetPoint(38,38.5,51.2873); gre->SetPointError(38,0,11.707); gre->SetPoint(39,39.5,36.0016); gre->SetPointError(39,0,9.62167); gre->SetPoint(40,40.5,18.6281); gre->SetPointError(40,0,6.31938); gre->SetPoint(41,41.5,0); gre->SetPointError(41,0,0); gre->SetPoint(42,42.5,-1); gre->SetPointError(42,0,-1); gre->SetPoint(43,43.5,-1); gre->SetPointError(43,0,-1); gre->SetPoint(44,44.5,-1); gre->SetPointError(44,0,-1); TF1 *fit = new TF1("*fit",2.1,40,8); //The original function : Energy_Spectrum had originally been created by: //TF1 *fit = new TF1("fit",Energy_Spectrum,2.1,40,8); fit->SetRange(2.1,40); fit->SetName("fit"); fit->SetTitle("Energy_Spectrum"); fit->SetSavedPoint(0,729.901); fit->SetSavedPoint(1,869.854); fit->SetSavedPoint(2,901.456); fit->SetSavedPoint(3,890.329); fit->SetSavedPoint(4,879.495); fit->SetSavedPoint(5,868.881); fit->SetSavedPoint(6,858.434); fit->SetSavedPoint(7,848.114); fit->SetSavedPoint(8,837.893); fit->SetSavedPoint(9,827.747); fit->SetSavedPoint(10,817.659); fit->SetSavedPoint(11,807.615); fit->SetSavedPoint(12,797.605); fit->SetSavedPoint(13,787.619); fit->SetSavedPoint(14,777.65); fit->SetSavedPoint(15,767.694); fit->SetSavedPoint(16,757.745); fit->SetSavedPoint(17,747.8); fit->SetSavedPoint(18,737.857); fit->SetSavedPoint(19,727.912); fit->SetSavedPoint(20,717.966); fit->SetSavedPoint(21,708.017); fit->SetSavedPoint(22,698.063); fit->SetSavedPoint(23,688.105); fit->SetSavedPoint(24,678.143); fit->SetSavedPoint(25,668.177); fit->SetSavedPoint(26,658.207); fit->SetSavedPoint(27,648.233); fit->SetSavedPoint(28,638.258); fit->SetSavedPoint(29,628.28); fit->SetSavedPoint(30,618.302); fit->SetSavedPoint(31,608.325); fit->SetSavedPoint(32,598.35); fit->SetSavedPoint(33,588.377); fit->SetSavedPoint(34,578.41); fit->SetSavedPoint(35,568.448); fit->SetSavedPoint(36,558.494); fit->SetSavedPoint(37,548.549); fit->SetSavedPoint(38,538.615); fit->SetSavedPoint(39,528.693); fit->SetSavedPoint(40,518.786); fit->SetSavedPoint(41,508.895); fit->SetSavedPoint(42,499.022); fit->SetSavedPoint(43,489.169); fit->SetSavedPoint(44,479.337); fit->SetSavedPoint(45,469.529); fit->SetSavedPoint(46,459.747); fit->SetSavedPoint(47,449.992); fit->SetSavedPoint(48,440.267); fit->SetSavedPoint(49,430.573); fit->SetSavedPoint(50,420.913); fit->SetSavedPoint(51,411.288); fit->SetSavedPoint(52,401.701); fit->SetSavedPoint(53,392.155); fit->SetSavedPoint(54,382.65); fit->SetSavedPoint(55,373.189); fit->SetSavedPoint(56,363.774); fit->SetSavedPoint(57,354.407); fit->SetSavedPoint(58,345.091); fit->SetSavedPoint(59,335.827); fit->SetSavedPoint(60,326.619); fit->SetSavedPoint(61,317.467); fit->SetSavedPoint(62,308.374); fit->SetSavedPoint(63,299.342); fit->SetSavedPoint(64,290.375); fit->SetSavedPoint(65,281.472); fit->SetSavedPoint(66,272.638); fit->SetSavedPoint(67,263.874); fit->SetSavedPoint(68,255.182); fit->SetSavedPoint(69,246.565); fit->SetSavedPoint(70,238.024); fit->SetSavedPoint(71,229.563); fit->SetSavedPoint(72,221.183); fit->SetSavedPoint(73,212.886); fit->SetSavedPoint(74,204.675); fit->SetSavedPoint(75,196.552); fit->SetSavedPoint(76,188.52); fit->SetSavedPoint(77,180.58); fit->SetSavedPoint(78,172.735); fit->SetSavedPoint(79,164.987); fit->SetSavedPoint(80,157.339); fit->SetSavedPoint(81,149.792); fit->SetSavedPoint(82,142.349); fit->SetSavedPoint(83,135.012); fit->SetSavedPoint(84,127.784); fit->SetSavedPoint(85,120.666); fit->SetSavedPoint(86,113.661); fit->SetSavedPoint(87,106.772); fit->SetSavedPoint(88,100.001); fit->SetSavedPoint(89,93.3492); fit->SetSavedPoint(90,86.8198); fit->SetSavedPoint(91,80.4149); fit->SetSavedPoint(92,74.1368); fit->SetSavedPoint(93,67.9878); fit->SetSavedPoint(94,61.9702); fit->SetSavedPoint(95,56.0863); fit->SetSavedPoint(96,50.3384); fit->SetSavedPoint(97,44.7289); fit->SetSavedPoint(98,39.2599); fit->SetSavedPoint(99,33.9339); fit->SetSavedPoint(100,28.7531); fit->SetSavedPoint(101,2.1); fit->SetSavedPoint(102,40); fit->SetFillColor(10); fit->SetFillStyle(0); fit->SetMarkerStyle(20); fit->SetMarkerSize(0.7); fit->SetLineColor(2); fit->SetLineWidth(2); fit->SetChisquare(48.9295); fit->GetXaxis()->SetNdivisions(506); fit->GetXaxis()->SetLabelFont(42); fit->GetXaxis()->SetLabelOffset(0.015); fit->GetXaxis()->SetLabelSize(0.06); fit->GetXaxis()->SetTitleSize(0.07); fit->GetXaxis()->SetTitleFont(42); fit->GetYaxis()->SetNdivisions(506); fit->GetYaxis()->SetLabelFont(42); fit->GetYaxis()->SetLabelOffset(0.015); fit->GetYaxis()->SetLabelSize(0.06); fit->GetYaxis()->SetTitleSize(0.07); fit->GetYaxis()->SetTitleOffset(1.1); fit->GetYaxis()->SetTitleFont(42); fit->SetParameter(0,17503.7); fit->SetParError(0,297.892); fit->SetParLimits(0,0,0); fit->SetParameter(1,0.940505); fit->SetParError(1,0.0648315); fit->SetParLimits(1,-1,1); fit->SetParameter(2,1.05343); fit->SetParError(2,0); fit->SetParLimits(2,1.05343,1.05343); fit->SetParameter(3,-0.0886778); fit->SetParError(3,0); fit->SetParLimits(3,-0.0886778,-0.0886778); fit->SetParameter(4,0.034432); fit->SetParError(4,0); fit->SetParLimits(4,0.034432,0.034432); fit->SetParameter(5,1.06497); fit->SetParError(5,0); fit->SetParLimits(5,1.06497,1.06497); fit->SetParameter(6,0.00313619); fit->SetParError(6,0); fit->SetParLimits(6,0.00313619,0.00313619); fit->SetParameter(7,-0.000364934); fit->SetParError(7,0); fit->SetParLimits(7,-0.000364934,-0.000364934); gre->GetListOfFunctions()->Add(fit); gre->Draw("p same"); TF1 *fit = new TF1("*fit",2.1,40,8); //The original function : Energy_Spectrum had originally been created by: //TF1 *fit = new TF1("fit",Energy_Spectrum,2.1,40,8); fit->SetRange(2.1,40); fit->SetName("fit"); fit->SetTitle("Energy_Spectrum"); fit->SetSavedPoint(0,729.901); fit->SetSavedPoint(1,869.854); fit->SetSavedPoint(2,901.456); fit->SetSavedPoint(3,890.329); fit->SetSavedPoint(4,879.495); fit->SetSavedPoint(5,868.881); fit->SetSavedPoint(6,858.434); fit->SetSavedPoint(7,848.114); fit->SetSavedPoint(8,837.893); fit->SetSavedPoint(9,827.747); fit->SetSavedPoint(10,817.659); fit->SetSavedPoint(11,807.615); fit->SetSavedPoint(12,797.605); fit->SetSavedPoint(13,787.619); fit->SetSavedPoint(14,777.65); fit->SetSavedPoint(15,767.694); fit->SetSavedPoint(16,757.745); fit->SetSavedPoint(17,747.8); fit->SetSavedPoint(18,737.857); fit->SetSavedPoint(19,727.912); fit->SetSavedPoint(20,717.966); fit->SetSavedPoint(21,708.017); fit->SetSavedPoint(22,698.063); fit->SetSavedPoint(23,688.105); fit->SetSavedPoint(24,678.143); fit->SetSavedPoint(25,668.177); fit->SetSavedPoint(26,658.207); fit->SetSavedPoint(27,648.233); fit->SetSavedPoint(28,638.258); fit->SetSavedPoint(29,628.28); fit->SetSavedPoint(30,618.302); fit->SetSavedPoint(31,608.325); fit->SetSavedPoint(32,598.35); fit->SetSavedPoint(33,588.377); fit->SetSavedPoint(34,578.41); fit->SetSavedPoint(35,568.448); fit->SetSavedPoint(36,558.494); fit->SetSavedPoint(37,548.549); fit->SetSavedPoint(38,538.615); fit->SetSavedPoint(39,528.693); fit->SetSavedPoint(40,518.786); fit->SetSavedPoint(41,508.895); fit->SetSavedPoint(42,499.022); fit->SetSavedPoint(43,489.169); fit->SetSavedPoint(44,479.337); fit->SetSavedPoint(45,469.529); fit->SetSavedPoint(46,459.747); fit->SetSavedPoint(47,449.992); fit->SetSavedPoint(48,440.267); fit->SetSavedPoint(49,430.573); fit->SetSavedPoint(50,420.913); fit->SetSavedPoint(51,411.288); fit->SetSavedPoint(52,401.701); fit->SetSavedPoint(53,392.155); fit->SetSavedPoint(54,382.65); fit->SetSavedPoint(55,373.189); fit->SetSavedPoint(56,363.774); fit->SetSavedPoint(57,354.407); fit->SetSavedPoint(58,345.091); fit->SetSavedPoint(59,335.827); fit->SetSavedPoint(60,326.619); fit->SetSavedPoint(61,317.467); fit->SetSavedPoint(62,308.374); fit->SetSavedPoint(63,299.342); fit->SetSavedPoint(64,290.375); fit->SetSavedPoint(65,281.472); fit->SetSavedPoint(66,272.638); fit->SetSavedPoint(67,263.874); fit->SetSavedPoint(68,255.182); fit->SetSavedPoint(69,246.565); fit->SetSavedPoint(70,238.024); fit->SetSavedPoint(71,229.563); fit->SetSavedPoint(72,221.183); fit->SetSavedPoint(73,212.886); fit->SetSavedPoint(74,204.675); fit->SetSavedPoint(75,196.552); fit->SetSavedPoint(76,188.52); fit->SetSavedPoint(77,180.58); fit->SetSavedPoint(78,172.735); fit->SetSavedPoint(79,164.987); fit->SetSavedPoint(80,157.339); fit->SetSavedPoint(81,149.792); fit->SetSavedPoint(82,142.349); fit->SetSavedPoint(83,135.012); fit->SetSavedPoint(84,127.784); fit->SetSavedPoint(85,120.666); fit->SetSavedPoint(86,113.661); fit->SetSavedPoint(87,106.772); fit->SetSavedPoint(88,100.001); fit->SetSavedPoint(89,93.3492); fit->SetSavedPoint(90,86.8198); fit->SetSavedPoint(91,80.4149); fit->SetSavedPoint(92,74.1368); fit->SetSavedPoint(93,67.9878); fit->SetSavedPoint(94,61.9702); fit->SetSavedPoint(95,56.0863); fit->SetSavedPoint(96,50.3384); fit->SetSavedPoint(97,44.7289); fit->SetSavedPoint(98,39.2599); fit->SetSavedPoint(99,33.9339); fit->SetSavedPoint(100,28.7531); fit->SetSavedPoint(101,2.1); fit->SetSavedPoint(102,40); fit->SetFillColor(10); fit->SetFillStyle(0); fit->SetMarkerStyle(20); fit->SetMarkerSize(0.7); fit->SetLineColor(2); fit->SetLineWidth(2); fit->SetChisquare(48.9295); fit->GetXaxis()->SetNdivisions(506); fit->GetXaxis()->SetLabelFont(42); fit->GetXaxis()->SetLabelOffset(0.015); fit->GetXaxis()->SetLabelSize(0.06); fit->GetXaxis()->SetTitleSize(0.07); fit->GetXaxis()->SetTitleFont(42); fit->GetYaxis()->SetNdivisions(506); fit->GetYaxis()->SetLabelFont(42); fit->GetYaxis()->SetLabelOffset(0.015); fit->GetYaxis()->SetLabelSize(0.06); fit->GetYaxis()->SetTitleSize(0.07); fit->GetYaxis()->SetTitleOffset(1.1); fit->GetYaxis()->SetTitleFont(42); fit->SetParameter(0,17503.7); fit->SetParError(0,297.892); fit->SetParLimits(0,0,0); fit->SetParameter(1,0.940505); fit->SetParError(1,0.0648315); fit->SetParLimits(1,-1,1); fit->SetParameter(2,1.05343); fit->SetParError(2,0); fit->SetParLimits(2,1.05343,1.05343); fit->SetParameter(3,-0.0886778); fit->SetParError(3,0); fit->SetParLimits(3,-0.0886778,-0.0886778); fit->SetParameter(4,0.034432); fit->SetParError(4,0); fit->SetParLimits(4,0.034432,0.034432); fit->SetParameter(5,1.06497); fit->SetParError(5,0); fit->SetParLimits(5,1.06497,1.06497); fit->SetParameter(6,0.00313619); fit->SetParError(6,0); fit->SetParLimits(6,0.00313619,0.00313619); fit->SetParameter(7,-0.000364934); fit->SetParError(7,0); fit->SetParLimits(7,-0.000364934,-0.000364934); fit->Draw("same"); THStack *hs = new THStack(); hs->SetName("hs"); hs->SetTitle("Energy Spectrum"); TH1 * b = new TH1F("","",45,0,45); b->SetDirectory(0); b->SetStats(0); b->SetLineWidth(2); b->SetMarkerStyle(20); b->SetMarkerSize(0.7); b->GetXaxis()->SetNdivisions(506); b->GetXaxis()->SetLabelFont(42); b->GetXaxis()->SetLabelOffset(0.015); b->GetXaxis()->SetLabelSize(0.06); b->GetXaxis()->SetTitleSize(0.07); b->GetXaxis()->SetTitleFont(42); b->GetYaxis()->SetNdivisions(506); b->GetYaxis()->SetLabelFont(42); b->GetYaxis()->SetLabelOffset(0.015); b->GetYaxis()->SetLabelSize(0.06); b->GetYaxis()->SetTitleSize(0.07); b->GetYaxis()->SetTitleOffset(1.1); b->GetYaxis()->SetTitleFont(42); b->GetZaxis()->SetLabelFont(42); b->GetZaxis()->SetLabelOffset(0.015); b->GetZaxis()->SetLabelSize(0.06); b->GetZaxis()->SetTitleSize(0.07); b->GetZaxis()->SetTitleOffset(1.1); b->GetZaxis()->SetTitleFont(42); // hs->Add(b,""); TH1 *c = new TH1F("","",45,0,45); c->SetDirectory(0); c->SetStats(0); c->SetLineColor(2); c->SetLineWidth(2); c->SetMarkerStyle(20); c->SetMarkerSize(0.7); c->GetXaxis()->SetNdivisions(506); c->GetXaxis()->SetLabelFont(42); c->GetXaxis()->SetLabelOffset(0.015); c->GetXaxis()->SetLabelSize(0.06); c->GetXaxis()->SetTitleSize(0.07); c->GetXaxis()->SetTitleFont(42); c->GetYaxis()->SetNdivisions(506); c->GetYaxis()->SetLabelFont(42); c->GetYaxis()->SetLabelOffset(0.015); c->GetYaxis()->SetLabelSize(0.06); c->GetYaxis()->SetTitleSize(0.07); c->GetYaxis()->SetTitleOffset(1.1); c->GetYaxis()->SetTitleFont(42); c->GetZaxis()->SetLabelFont(42); c->GetZaxis()->SetLabelOffset(0.015); c->GetZaxis()->SetLabelSize(0.06); c->GetZaxis()->SetTitleSize(0.07); c->GetZaxis()->SetTitleOffset(1.1); c->GetZaxis()->SetTitleFont(42); // hs->Add(c,""); TH1 *d = new TH1F("","",45,0,45); d->SetDirectory(0); d->SetStats(0); d->SetLineColor(2); d->SetLineWidth(2); d->SetMarkerStyle(20); d->SetMarkerSize(0.7); d->GetXaxis()->SetNdivisions(506); d->GetXaxis()->SetLabelFont(42); d->GetXaxis()->SetLabelOffset(0.015); d->GetXaxis()->SetLabelSize(0.06); d->GetXaxis()->SetTitleSize(0.07); d->GetXaxis()->SetTitleFont(42); d->GetYaxis()->SetNdivisions(506); d->GetYaxis()->SetLabelFont(42); d->GetYaxis()->SetLabelOffset(0.015); d->GetYaxis()->SetLabelSize(0.06); d->GetYaxis()->SetTitleSize(0.07); d->GetYaxis()->SetTitleOffset(1.1); d->GetYaxis()->SetTitleFont(42); d->GetZaxis()->SetLabelFont(42); d->GetZaxis()->SetLabelOffset(0.015); d->GetZaxis()->SetLabelSize(0.06); d->GetZaxis()->SetTitleSize(0.07); d->GetZaxis()->SetTitleOffset(1.1); d->GetZaxis()->SetTitleFont(42); //hs->Add(d,""); TH1 *e = new TH1F("","",45,0,45); e->SetDirectory(0); e->SetStats(0); e->SetLineColor(2); e->SetLineWidth(2); e->SetMarkerStyle(20); e->SetMarkerSize(0.7); e->GetXaxis()->SetNdivisions(506); e->GetXaxis()->SetLabelFont(42); e->GetXaxis()->SetLabelOffset(0.015); e->GetXaxis()->SetLabelSize(0.06); e->GetXaxis()->SetTitleSize(0.07); e->GetXaxis()->SetTitleFont(42); e->GetYaxis()->SetNdivisions(506); e->GetYaxis()->SetLabelFont(42); e->GetYaxis()->SetLabelOffset(0.015); e->GetYaxis()->SetLabelSize(0.06); e->GetYaxis()->SetTitleSize(0.07); e->GetYaxis()->SetTitleOffset(1.1); e->GetYaxis()->SetTitleFont(42); e->GetZaxis()->SetLabelFont(42); e->GetZaxis()->SetLabelOffset(0.015); e->GetZaxis()->SetLabelSize(0.06); e->GetZaxis()->SetTitleSize(0.07); e->GetZaxis()->SetTitleOffset(1.1); e->GetZaxis()->SetTitleFont(42); // hs->Add(e,""); TH1 *f = new TH1F("","",45,0,45); f->SetBinContent(3,7.80783); f->SetBinContent(4,14.8459); f->SetBinContent(5,30.2941); f->SetBinContent(6,40.7615); f->SetBinContent(7,41.6118); f->SetBinContent(8,32.4833); f->SetBinContent(9,25.3291); f->SetBinContent(10,20.632); f->SetBinContent(11,20.7074); f->SetBinContent(12,20.327); f->SetBinContent(13,16.2816); f->SetBinContent(14,28.502); f->SetBinContent(15,20.0982); f->SetBinContent(16,14.9938); f->SetBinContent(17,10.6923); f->SetBinContent(18,11.5888); f->SetBinContent(19,11.9786); f->SetBinContent(20,9.14828); f->SetBinContent(21,5.66825); f->SetBinContent(22,3.58671); f->SetBinContent(23,9.35276); f->SetBinContent(24,2.98135); f->SetBinContent(25,5.55054); f->SetBinContent(26,3.34129); f->SetBinContent(27,2.10083); f->SetBinContent(28,1.15474); f->SetBinContent(29,1.56067); f->SetBinContent(30,0.585712); f->SetBinContent(31,3.32561); f->SetBinContent(32,1.8837); f->SetBinContent(33,4.6588); f->SetBinContent(34,15.9151); f->SetBinContent(35,3.09237); f->SetBinContent(36,1.60375); f->SetBinContent(37,0.600229); f->SetBinContent(38,1.34623); f->SetBinContent(39,1.32078); f->SetBinContent(40,1.33161); f->SetBinContent(41,1.03705); f->SetBinContent(42,0.495724); f->SetBinContent(43,0.309259); f->SetEntries(27923); f->SetDirectory(0); f->SetStats(0); f->SetFillColor(2); f->SetLineWidth(2); f->SetMarkerStyle(20); f->SetMarkerSize(0.7); f->GetXaxis()->SetNdivisions(506); f->GetXaxis()->SetLabelFont(42); f->GetXaxis()->SetLabelOffset(0.015); f->GetXaxis()->SetLabelSize(0.06); f->GetXaxis()->SetTitleSize(0.07); f->GetXaxis()->SetTitleFont(42); f->GetYaxis()->SetNdivisions(506); f->GetYaxis()->SetLabelFont(42); f->GetYaxis()->SetLabelOffset(0.015); f->GetYaxis()->SetLabelSize(0.06); f->GetYaxis()->SetTitleSize(0.07); f->GetYaxis()->SetTitleOffset(1.1); f->GetYaxis()->SetTitleFont(42); f->GetZaxis()->SetLabelFont(42); f->GetZaxis()->SetLabelOffset(0.015); f->GetZaxis()->SetLabelSize(0.06); f->GetZaxis()->SetTitleSize(0.07); f->GetZaxis()->SetTitleOffset(1.1); f->GetZaxis()->SetTitleFont(42); hs->Add(f,""); TH1 *g = new TH1F("","",45,0,45); g->SetBinContent(3,6.05209); g->SetBinContent(4,35.6424); g->SetBinContent(5,69.2899); g->SetBinContent(6,93.597); g->SetBinContent(7,104.193); g->SetBinContent(8,108.416); g->SetBinContent(9,107.511); g->SetBinContent(10,119.831); g->SetBinContent(11,98.9356); g->SetBinContent(12,103.739); g->SetBinContent(13,98.354); g->SetBinContent(14,111.214); g->SetBinContent(15,108.112); g->SetBinContent(16,103.861); g->SetBinContent(17,101.697); g->SetBinContent(18,92.5149); g->SetBinContent(19,94.0081); g->SetBinContent(20,89.6898); g->SetBinContent(21,88.592); g->SetBinContent(22,82.2635); g->SetBinContent(23,75.4745); g->SetBinContent(24,63.9704); g->SetBinContent(25,60.4377); g->SetBinContent(26,67.2868); g->SetBinContent(27,53.9182); g->SetBinContent(28,48.2042); g->SetBinContent(29,50.418); g->SetBinContent(30,41.5915); g->SetBinContent(31,33.2); g->SetBinContent(32,32.3282); g->SetBinContent(33,34.1196); g->SetBinContent(34,20.3952); g->SetBinContent(35,15.5164); g->SetBinContent(36,18.6249); g->SetBinContent(37,21.2909); g->SetBinContent(38,8.86635); g->SetBinContent(39,8.88208); g->SetBinContent(40,6.18048); g->SetBinContent(41,3.47314); g->SetBinContent(42,0.835862); g->SetEntries(5746); g->SetDirectory(0); g->SetStats(0); g->SetFillColor(5); g->SetLineWidth(2); g->SetMarkerStyle(20); g->SetMarkerSize(0.7); g->GetXaxis()->SetNdivisions(506); g->GetXaxis()->SetLabelFont(42); g->GetXaxis()->SetLabelOffset(0.015); g->GetXaxis()->SetLabelSize(0.06); g->GetXaxis()->SetTitleSize(0.07); g->GetXaxis()->SetTitleFont(42); g->GetYaxis()->SetNdivisions(506); g->GetYaxis()->SetLabelFont(42); g->GetYaxis()->SetLabelOffset(0.015); g->GetYaxis()->SetLabelSize(0.06); g->GetYaxis()->SetTitleSize(0.07); g->GetYaxis()->SetTitleOffset(1.1); g->GetYaxis()->SetTitleFont(42); g->GetZaxis()->SetLabelFont(42); g->GetZaxis()->SetLabelOffset(0.015); g->GetZaxis()->SetLabelSize(0.06); g->GetZaxis()->SetTitleSize(0.07); g->GetZaxis()->SetTitleOffset(1.1); g->GetZaxis()->SetTitleFont(42); hs->Add(g,""); hs->Draw("same"); TLegend *leg = new TLegend(0.5,0.62,0.89,0.89,NULL,"brNDC"); leg->SetBorderSize(0); leg->SetLineColor(1); leg->SetLineStyle(1); leg->SetLineWidth(2); leg->SetTextSize(0.035); leg->SetFillColor(10); leg->SetFillStyle(1001); TLegendEntry *entry=leg->AddEntry("g","Signal","F"); entry->SetFillColor(5); entry->SetFillStyle(1001); entry->SetLineColor(1); entry->SetLineStyle(1); entry->SetLineWidth(2); entry->SetMarkerColor(1); entry->SetMarkerStyle(21); entry->SetMarkerSize(1); entry->SetTextAlign(12); entry->SetTextColor(1); entry=leg->AddEntry("b","Background","F"); entry->SetFillColor(2); entry->SetFillStyle(1001); entry->SetLineColor(1); entry->SetLineStyle(1); entry->SetLineWidth(2); entry->SetMarkerColor(1); entry->SetMarkerStyle(21); entry->SetMarkerSize(1); entry->SetTextAlign(12); entry->SetTextColor(1); entry=leg->AddEntry("a","(Signal-Back.)/Effic.","PE"); entry->SetLineColor(1); entry->SetLineStyle(1); entry->SetLineWidth(1); entry->SetMarkerColor(1); entry->SetMarkerStyle(24); entry->SetMarkerSize(0.7); entry->SetTextAlign(12); entry->SetTextColor(1); entry=leg->AddEntry("f","Fit","L"); entry->SetLineColor(2); entry->SetLineStyle(1); entry->SetLineWidth(3); entry->SetMarkerColor(1); entry->SetMarkerStyle(0); entry->SetMarkerSize(0); entry->SetTextAlign(12); entry->SetTextColor(1); leg->Draw(); gre = new TGraphErrors(45); gre->SetName("Graph"); gre->SetTitle("Graph"); gre->SetFillColor(1); gre->SetLineWidth(2); gre->SetMarkerStyle(24); gre->SetMarkerSize(0.7); gre->SetPoint(0,0.5,-1); gre->SetPointError(0,0,-1); gre->SetPoint(1,1.5,-1); gre->SetPointError(1,0,-1); gre->SetPoint(2,2.5,835.383); gre->SetPointError(2,0,169.785); gre->SetPoint(3,3.5,811.232); gre->SetPointError(3,0,83.5794); gre->SetPoint(4,4.5,820.985); gre->SetPointError(4,0,60.2182); gre->SetPoint(5,5.5,909.631); gre->SetPointError(5,0,57.5381); gre->SetPoint(6,6.5,895.972); gre->SetPointError(6,0,54.0836); gre->SetPoint(7,7.5,805.348); gre->SetPointError(7,0,48.992); gre->SetPoint(8,8.5,713.819); gre->SetPointError(8,0,44.3519); gre->SetPoint(9,9.5,747.495); gre->SetPointError(9,0,45.323); gre->SetPoint(10,10.5,621.392); gre->SetPointError(10,0,40.4441); gre->SetPoint(11,11.5,642.512); gre->SetPointError(11,0,41.1759); gre->SetPoint(12,12.5,574.956); gre->SetPointError(12,0,38.0994); gre->SetPoint(13,13.5,711.834); gre->SetPointError(13,0,43.3758); gre->SetPoint(14,14.5,632.053); gre->SetPointError(14,0,39.9242); gre->SetPoint(15,15.5,571.546); gre->SetPointError(15,0,37.3128); gre->SetPoint(16,16.5,532.478); gre->SetPointError(16,0,35.6571); gre->SetPoint(17,17.5,487.067); gre->SetPointError(17,0,33.7854); gre->SetPoint(18,18.5,497.454); gre->SetPointError(18,0,34.3219); gre->SetPoint(19,19.5,461.82); gre->SetPointError(19,0,32.9216); gre->SetPoint(20,20.5,442.087); gre->SetPointError(20,0,32.269); gre->SetPoint(21,21.5,404.592); gre->SetPointError(21,0,30.8758); gre->SetPoint(22,22.5,406.751); gre->SetPointError(22,0,31.3749); gre->SetPoint(23,23.5,318.151); gre->SetPointError(23,0,27.3149); gre->SetPoint(24,24.5,317.223); gre->SetPointError(24,0,27.544); gre->SetPoint(25,25.5,344.974); gre->SetPointError(25,0,29.2281); gre->SetPoint(26,26.5,270.047); gre->SetPointError(26,0,25.4218); gre->SetPoint(27,27.5,238.374); gre->SetPointError(27,0,23.8486); gre->SetPoint(28,28.5,257.966); gre->SetPointError(28,0,25.4101); gre->SetPoint(29,29.5,210.854); gre->SetPointError(29,0,22.9047); gre->SetPoint(30,30.5,186.016); gre->SetPointError(30,0,21.7022); gre->SetPoint(31,31.5,179.516); gre->SetPointError(31,0,21.7412); gre->SetPoint(32,32.5,213.506); gre->SetPointError(32,0,24.8319); gre->SetPoint(33,33.5,204.497); gre->SetPointError(33,0,24.8603); gre->SetPoint(34,34.5,98.5616); gre->SetPointError(34,0,16.074); gre->SetPoint(35,35.5,109.477); gre->SetPointError(35,0,17.4417); gre->SetPoint(36,36.5,120.173); gre->SetPointError(36,0,18.9229); gre->SetPoint(37,37.5,50.1321); gre->SetPointError(37,0,11.124); gre->SetPoint(38,38.5,51.2873); gre->SetPointError(38,0,11.707); gre->SetPoint(39,39.5,36.0016); gre->SetPointError(39,0,9.62167); gre->SetPoint(40,40.5,18.6281); gre->SetPointError(40,0,6.31938); gre->SetPoint(41,41.5,0); gre->SetPointError(41,0,0); gre->SetPoint(42,42.5,-1); gre->SetPointError(42,0,-1); gre->SetPoint(43,43.5,-1); gre->SetPointError(43,0,-1); gre->SetPoint(44,44.5,-1); gre->SetPointError(44,0,-1); TF1 *fit = new TF1("*fit",2.1,40,8); //The original function : Energy_Spectrum had originally been created by: //TF1 *fit = new TF1("fit",Energy_Spectrum,2.1,40,8); fit->SetRange(2.1,40); fit->SetName("fit"); fit->SetTitle("Energy_Spectrum"); fit->SetSavedPoint(0,729.901); fit->SetSavedPoint(1,869.854); fit->SetSavedPoint(2,901.456); fit->SetSavedPoint(3,890.329); fit->SetSavedPoint(4,879.495); fit->SetSavedPoint(5,868.881); fit->SetSavedPoint(6,858.434); fit->SetSavedPoint(7,848.114); fit->SetSavedPoint(8,837.893); fit->SetSavedPoint(9,827.747); fit->SetSavedPoint(10,817.659); fit->SetSavedPoint(11,807.615); fit->SetSavedPoint(12,797.605); fit->SetSavedPoint(13,787.619); fit->SetSavedPoint(14,777.65); fit->SetSavedPoint(15,767.694); fit->SetSavedPoint(16,757.745); fit->SetSavedPoint(17,747.8); fit->SetSavedPoint(18,737.857); fit->SetSavedPoint(19,727.912); fit->SetSavedPoint(20,717.966); fit->SetSavedPoint(21,708.017); fit->SetSavedPoint(22,698.063); fit->SetSavedPoint(23,688.105); fit->SetSavedPoint(24,678.143); fit->SetSavedPoint(25,668.177); fit->SetSavedPoint(26,658.207); fit->SetSavedPoint(27,648.233); fit->SetSavedPoint(28,638.258); fit->SetSavedPoint(29,628.28); fit->SetSavedPoint(30,618.302); fit->SetSavedPoint(31,608.325); fit->SetSavedPoint(32,598.35); fit->SetSavedPoint(33,588.377); fit->SetSavedPoint(34,578.41); fit->SetSavedPoint(35,568.448); fit->SetSavedPoint(36,558.494); fit->SetSavedPoint(37,548.549); fit->SetSavedPoint(38,538.615); fit->SetSavedPoint(39,528.693); fit->SetSavedPoint(40,518.786); fit->SetSavedPoint(41,508.895); fit->SetSavedPoint(42,499.022); fit->SetSavedPoint(43,489.169); fit->SetSavedPoint(44,479.337); fit->SetSavedPoint(45,469.529); fit->SetSavedPoint(46,459.747); fit->SetSavedPoint(47,449.992); fit->SetSavedPoint(48,440.267); fit->SetSavedPoint(49,430.573); fit->SetSavedPoint(50,420.913); fit->SetSavedPoint(51,411.288); fit->SetSavedPoint(52,401.701); fit->SetSavedPoint(53,392.155); fit->SetSavedPoint(54,382.65); fit->SetSavedPoint(55,373.189); fit->SetSavedPoint(56,363.774); fit->SetSavedPoint(57,354.407); fit->SetSavedPoint(58,345.091); fit->SetSavedPoint(59,335.827); fit->SetSavedPoint(60,326.619); fit->SetSavedPoint(61,317.467); fit->SetSavedPoint(62,308.374); fit->SetSavedPoint(63,299.342); fit->SetSavedPoint(64,290.375); fit->SetSavedPoint(65,281.472); fit->SetSavedPoint(66,272.638); fit->SetSavedPoint(67,263.874); fit->SetSavedPoint(68,255.182); fit->SetSavedPoint(69,246.565); fit->SetSavedPoint(70,238.024); fit->SetSavedPoint(71,229.563); fit->SetSavedPoint(72,221.183); fit->SetSavedPoint(73,212.886); fit->SetSavedPoint(74,204.675); fit->SetSavedPoint(75,196.552); fit->SetSavedPoint(76,188.52); fit->SetSavedPoint(77,180.58); fit->SetSavedPoint(78,172.735); fit->SetSavedPoint(79,164.987); fit->SetSavedPoint(80,157.339); fit->SetSavedPoint(81,149.792); fit->SetSavedPoint(82,142.349); fit->SetSavedPoint(83,135.012); fit->SetSavedPoint(84,127.784); fit->SetSavedPoint(85,120.666); fit->SetSavedPoint(86,113.661); fit->SetSavedPoint(87,106.772); fit->SetSavedPoint(88,100.001); fit->SetSavedPoint(89,93.3492); fit->SetSavedPoint(90,86.8198); fit->SetSavedPoint(91,80.4149); fit->SetSavedPoint(92,74.1368); fit->SetSavedPoint(93,67.9878); fit->SetSavedPoint(94,61.9702); fit->SetSavedPoint(95,56.0863); fit->SetSavedPoint(96,50.3384); fit->SetSavedPoint(97,44.7289); fit->SetSavedPoint(98,39.2599); fit->SetSavedPoint(99,33.9339); fit->SetSavedPoint(100,28.7531); fit->SetSavedPoint(101,2.1); fit->SetSavedPoint(102,40); fit->SetFillColor(10); fit->SetFillStyle(0); fit->SetMarkerStyle(20); fit->SetMarkerSize(0.7); fit->SetLineColor(2); fit->SetLineWidth(2); fit->SetChisquare(48.9295); fit->GetXaxis()->SetNdivisions(506); fit->GetXaxis()->SetLabelFont(42); fit->GetXaxis()->SetLabelOffset(0.015); fit->GetXaxis()->SetLabelSize(0.06); fit->GetXaxis()->SetTitleSize(0.07); fit->GetXaxis()->SetTitleFont(42); fit->GetYaxis()->SetNdivisions(506); fit->GetYaxis()->SetLabelFont(42); fit->GetYaxis()->SetLabelOffset(0.015); fit->GetYaxis()->SetLabelSize(0.06); fit->GetYaxis()->SetTitleSize(0.07); fit->GetYaxis()->SetTitleOffset(1.1); fit->GetYaxis()->SetTitleFont(42); fit->SetParameter(0,17503.7); fit->SetParError(0,297.892); fit->SetParLimits(0,0,0); fit->SetParameter(1,0.940505); fit->SetParError(1,0.0648315); fit->SetParLimits(1,-1,1); fit->SetParameter(2,1.05343); fit->SetParError(2,0); fit->SetParLimits(2,1.05343,1.05343); fit->SetParameter(3,-0.0886778); fit->SetParError(3,0); fit->SetParLimits(3,-0.0886778,-0.0886778); fit->SetParameter(4,0.034432); fit->SetParError(4,0); fit->SetParLimits(4,0.034432,0.034432); fit->SetParameter(5,1.06497); fit->SetParError(5,0); fit->SetParLimits(5,1.06497,1.06497); fit->SetParameter(6,0.00313619); fit->SetParError(6,0); fit->SetParLimits(6,0.00313619,0.00313619); fit->SetParameter(7,-0.000364934); fit->SetParError(7,0); fit->SetParLimits(7,-0.000364934,-0.000364934); gre->GetListOfFunctions()->Add(fit); gre->Draw("p same"); tex = new TLatex(40,1000,"b)"); tex->SetTextFont(42); tex->SetLineWidth(2); tex->Draw(); c1->Modified(); c1->cd(); c1->SaveAs("Stau_Spectrum.gif"); c1->SetSelected(c1); }