Micromegas LAPP reconstruction

00001 {
00002 //=========Macro generated from canvas: c0/Prototype m2
00003 //=========  (Tue Dec 13 10:32:22 2011) by ROOT version5.28/00
00004    TCanvas *c0 = new TCanvas("c0", "Prototype m2",11,37,700,500);
00005    gStyle->SetOptFit(1);
00006    gStyle->SetOptStat(0);
00007    c0->SetHighLightColor(2);
00008    c0->Range(874,-26.25,2134,236.25);
00009    c0->SetFillColor(0);
00010    c0->SetBorderMode(0);
00011    c0->SetBorderSize(2);
00012    c0->SetFrameFillStyle(0);
00013    c0->SetFrameBorderMode(0);
00014    c0->SetFrameFillStyle(0);
00015    c0->SetFrameBorderMode(0);
00016    
00017    TH1I *ped = new TH1I("ped","",512,0,4096);
00018    ped->SetBinContent(145,3);
00019    ped->SetBinContent(146,133);
00020    ped->SetBinContent(147,115);
00021    ped->SetBinContent(150,4);
00022    ped->SetBinContent(151,2);
00023    ped->SetBinContent(152,2);
00024    ped->SetBinContent(153,2);
00025    ped->SetBinContent(154,2);
00026    ped->SetBinContent(155,2);
00027    ped->SetBinContent(156,2);
00028    ped->SetBinContent(157,2);
00029    ped->SetBinContent(158,2);
00030    ped->SetBinContent(159,1);
00031    ped->SetBinContent(160,1);
00032    ped->SetBinContent(162,3);
00033    ped->SetBinContent(163,2);
00034    ped->SetBinContent(166,3);
00035    ped->SetBinContent(167,1);
00036    ped->SetBinContent(171,1);
00037    ped->SetBinContent(174,1);
00038    ped->SetBinContent(175,1);
00039    ped->SetBinContent(178,1);
00040    ped->SetBinContent(180,2);
00041    ped->SetBinContent(183,1);
00042    ped->SetBinContent(184,1);
00043    ped->SetBinContent(186,1);
00044    ped->SetBinContent(187,2);
00045    ped->SetBinContent(188,1);
00046    ped->SetBinContent(189,1);
00047    ped->SetBinContent(197,1);
00048    ped->SetBinContent(207,1);
00049    ped->SetBinContent(239,1);
00050    ped->SetBinContent(336,1);
00051    ped->SetBinContent(478,1);
00052    ped->SetMinimum(0);
00053    ped->SetMaximum(210);
00054    ped->SetEntries(300);
00055    ped->SetStats(0);
00056    
00057    TF1 *gaus = new TF1("gaus","gaus",1140,1200);
00058    gaus->SetFillColor(19);
00059    gaus->SetFillStyle(0);
00060 
00061    Int_t ci;   // for color index setting
00062    ci = TColor::GetColor("#0000cc");
00063    gaus->SetLineColor(ci);
00064    gaus->SetLineWidth(3);
00065    gaus->SetChisquare(4);
00066    gaus->SetNDF(1);
00067    gaus->SetParameter(0,202.8509);
00068    gaus->SetParError(0,21.07621);
00069    gaus->SetParLimits(0,0,0);
00070    gaus->SetParameter(1,1167.705);
00071    gaus->SetParError(1,0.2486121);
00072    gaus->SetParLimits(1,0,0);
00073    gaus->SetParameter(2,4.031796);
00074    gaus->SetParError(2,0.3123966);
00075    gaus->SetParLimits(2,0,54.36239);
00076    ped->GetListOfFunctions()->Add(gaus);
00077    ped->GetXaxis()->SetTitle("ADC count");
00078    ped->GetXaxis()->SetRange(126,251);
00079    ped->GetYaxis()->SetTitle("N");
00080    ped->Draw("");
00081    
00082    TH1I *ped = new TH1I("ped","",512,0,4096);
00083    ped->SetBinContent(145,3);
00084    ped->SetBinContent(146,133);
00085    ped->SetBinContent(147,115);
00086    ped->SetBinContent(150,4);
00087    ped->SetBinContent(151,2);
00088    ped->SetBinContent(152,2);
00089    ped->SetBinContent(153,2);
00090    ped->SetBinContent(154,2);
00091    ped->SetBinContent(155,2);
00092    ped->SetBinContent(156,2);
00093    ped->SetBinContent(157,2);
00094    ped->SetBinContent(158,2);
00095    ped->SetBinContent(159,1);
00096    ped->SetBinContent(160,1);
00097    ped->SetBinContent(162,3);
00098    ped->SetBinContent(163,2);
00099    ped->SetBinContent(166,3);
00100    ped->SetBinContent(167,1);
00101    ped->SetBinContent(171,1);
00102    ped->SetBinContent(174,1);
00103    ped->SetBinContent(175,1);
00104    ped->SetBinContent(178,1);
00105    ped->SetBinContent(180,2);
00106    ped->SetBinContent(183,1);
00107    ped->SetBinContent(184,1);
00108    ped->SetBinContent(186,1);
00109    ped->SetBinContent(187,2);
00110    ped->SetBinContent(188,1);
00111    ped->SetBinContent(189,1);
00112    ped->SetBinContent(197,1);
00113    ped->SetBinContent(207,1);
00114    ped->SetBinContent(239,1);
00115    ped->SetBinContent(336,1);
00116    ped->SetBinContent(478,1);
00117    ped->SetMinimum(0);
00118    ped->SetMaximum(210);
00119    ped->SetEntries(300);
00120    ped->SetStats(0);
00121    
00122    TF1 *gaus = new TF1("gaus","gaus",1140,1200);
00123    gaus->SetFillColor(19);
00124    gaus->SetFillStyle(0);
00125 
00126    ci = TColor::GetColor("#0000cc");
00127    gaus->SetLineColor(ci);
00128    gaus->SetLineWidth(3);
00129    gaus->SetChisquare(4);
00130    gaus->SetNDF(1);
00131    gaus->SetParameter(0,202.8509);
00132    gaus->SetParError(0,21.07621);
00133    gaus->SetParLimits(0,0,0);
00134    gaus->SetParameter(1,1167.705);
00135    gaus->SetParError(1,0.2486121);
00136    gaus->SetParLimits(1,0,0);
00137    gaus->SetParameter(2,4.031796);
00138    gaus->SetParError(2,0.3123966);
00139    gaus->SetParLimits(2,0,54.36239);
00140    ped->GetListOfFunctions()->Add(gaus);
00141    ped->GetXaxis()->SetTitle("ADC count");
00142    ped->GetXaxis()->SetRange(126,251);
00143    ped->GetYaxis()->SetTitle("N");
00144    ped->Draw("SAME");
00145    
00146    TH1I *all = new TH1I("all","",512,0,4096);
00147    all->SetBinContent(145,1);
00148    all->SetBinContent(146,11);
00149    all->SetBinContent(147,24);
00150    all->SetBinContent(148,68);
00151    all->SetBinContent(149,69);
00152    all->SetBinContent(150,89);
00153    all->SetBinContent(151,89);
00154    all->SetBinContent(152,134);
00155    all->SetBinContent(153,121);
00156    all->SetBinContent(154,131);
00157    all->SetBinContent(155,102);
00158    all->SetBinContent(156,125);
00159    all->SetBinContent(157,127);
00160    all->SetBinContent(158,124);
00161    all->SetBinContent(159,91);
00162    all->SetBinContent(160,113);
00163    all->SetBinContent(161,102);
00164    all->SetBinContent(162,82);
00165    all->SetBinContent(163,96);
00166    all->SetBinContent(164,90);
00167    all->SetBinContent(165,65);
00168    all->SetBinContent(166,59);
00169    all->SetBinContent(167,75);
00170    all->SetBinContent(168,55);
00171    all->SetBinContent(169,71);
00172    all->SetBinContent(170,57);
00173    all->SetBinContent(171,49);
00174    all->SetBinContent(172,51);
00175    all->SetBinContent(173,43);
00176    all->SetBinContent(174,48);
00177    all->SetBinContent(175,43);
00178    all->SetBinContent(176,21);
00179    all->SetBinContent(177,27);
00180    all->SetBinContent(178,33);
00181    all->SetBinContent(179,21);
00182    all->SetBinContent(180,23);
00183    all->SetBinContent(181,16);
00184    all->SetBinContent(182,20);
00185    all->SetBinContent(183,19);
00186    all->SetBinContent(184,22);
00187    all->SetBinContent(185,15);
00188    all->SetBinContent(186,15);
00189    all->SetBinContent(187,14);
00190    all->SetBinContent(188,14);
00191    all->SetBinContent(189,13);
00192    all->SetBinContent(190,15);
00193    all->SetBinContent(191,8);
00194    all->SetBinContent(192,8);
00195    all->SetBinContent(193,11);
00196    all->SetBinContent(194,9);
00197    all->SetBinContent(195,11);
00198    all->SetBinContent(196,7);
00199    all->SetBinContent(197,6);
00200    all->SetBinContent(198,10);
00201    all->SetBinContent(199,9);
00202    all->SetBinContent(200,7);
00203    all->SetBinContent(201,6);
00204    all->SetBinContent(202,3);
00205    all->SetBinContent(203,4);
00206    all->SetBinContent(204,8);
00207    all->SetBinContent(205,8);
00208    all->SetBinContent(206,4);
00209    all->SetBinContent(207,6);
00210    all->SetBinContent(208,5);
00211    all->SetBinContent(209,4);
00212    all->SetBinContent(210,6);
00213    all->SetBinContent(211,9);
00214    all->SetBinContent(212,4);
00215    all->SetBinContent(213,4);
00216    all->SetBinContent(214,2);
00217    all->SetBinContent(215,2);
00218    all->SetBinContent(217,1);
00219    all->SetBinContent(228,1);
00220    all->SetBinContent(230,1);
00221    all->SetBinContent(236,1);
00222    all->SetBinContent(239,1);
00223    all->SetBinContent(272,1);
00224    all->SetBinContent(277,1);
00225    all->SetBinContent(286,1);
00226    all->SetBinContent(299,1);
00227    all->SetBinContent(309,1);
00228    all->SetBinContent(311,1);
00229    all->SetBinContent(318,1);
00230    all->SetBinContent(325,1);
00231    all->SetBinContent(329,1);
00232    all->SetBinContent(336,1);
00233    all->SetBinContent(344,1);
00234    all->SetBinContent(345,1);
00235    all->SetBinContent(346,1);
00236    all->SetBinContent(356,1);
00237    all->SetBinContent(364,1);
00238    all->SetBinContent(370,1);
00239    all->SetBinContent(397,1);
00240    all->SetBinContent(399,1);
00241    all->SetBinContent(412,1);
00242    all->SetBinContent(419,1);
00243    all->SetBinContent(431,1);
00244    all->SetBinContent(432,1);
00245    all->SetBinContent(436,1);
00246    all->SetBinContent(437,1);
00247    all->SetBinContent(438,1);
00248    all->SetBinContent(439,1);
00249    all->SetBinContent(440,2);
00250    all->SetBinContent(441,1);
00251    all->SetBinContent(443,1);
00252    all->SetBinContent(446,1);
00253    all->SetBinContent(448,2);
00254    all->SetBinContent(451,1);
00255    all->SetBinContent(453,1);
00256    all->SetBinContent(460,3);
00257    all->SetBinContent(465,1);
00258    all->SetBinContent(467,1);
00259    all->SetBinContent(475,2);
00260    all->SetBinContent(477,1);
00261    all->SetBinContent(479,1);
00262    all->SetBinContent(480,1);
00263    all->SetEntries(2904);
00264    all->SetStats(0);
00265    
00266    TF1 *landau = new TF1("landau","landau",1100,1800);
00267    landau->SetFillColor(19);
00268    landau->SetFillStyle(0);
00269 
00270    ci = TColor::GetColor("#cc6633");
00271    landau->SetLineColor(ci);
00272    landau->SetLineWidth(3);
00273    landau->SetChisquare(109.5263);
00274    landau->SetNDF(69);
00275    landau->SetParameter(0,736.9045);
00276    landau->SetParError(0,19.44123);
00277    landau->SetParLimits(0,0,0);
00278    landau->SetParameter(1,1241.799);
00279    landau->SetParError(1,1.408543);
00280    landau->SetParLimits(1,0,0);
00281    landau->SetParameter(2,32.28482);
00282    landau->SetParError(2,0.6532016);
00283    landau->SetParLimits(2,0,1069.392);
00284    all->GetListOfFunctions()->Add(landau);
00285    all->Draw("sames");
00286    
00287    TH1I *all = new TH1I("all","",512,0,4096);
00288    all->SetBinContent(145,1);
00289    all->SetBinContent(146,11);
00290    all->SetBinContent(147,24);
00291    all->SetBinContent(148,68);
00292    all->SetBinContent(149,69);
00293    all->SetBinContent(150,89);
00294    all->SetBinContent(151,89);
00295    all->SetBinContent(152,134);
00296    all->SetBinContent(153,121);
00297    all->SetBinContent(154,131);
00298    all->SetBinContent(155,102);
00299    all->SetBinContent(156,125);
00300    all->SetBinContent(157,127);
00301    all->SetBinContent(158,124);
00302    all->SetBinContent(159,91);
00303    all->SetBinContent(160,113);
00304    all->SetBinContent(161,102);
00305    all->SetBinContent(162,82);
00306    all->SetBinContent(163,96);
00307    all->SetBinContent(164,90);
00308    all->SetBinContent(165,65);
00309    all->SetBinContent(166,59);
00310    all->SetBinContent(167,75);
00311    all->SetBinContent(168,55);
00312    all->SetBinContent(169,71);
00313    all->SetBinContent(170,57);
00314    all->SetBinContent(171,49);
00315    all->SetBinContent(172,51);
00316    all->SetBinContent(173,43);
00317    all->SetBinContent(174,48);
00318    all->SetBinContent(175,43);
00319    all->SetBinContent(176,21);
00320    all->SetBinContent(177,27);
00321    all->SetBinContent(178,33);
00322    all->SetBinContent(179,21);
00323    all->SetBinContent(180,23);
00324    all->SetBinContent(181,16);
00325    all->SetBinContent(182,20);
00326    all->SetBinContent(183,19);
00327    all->SetBinContent(184,22);
00328    all->SetBinContent(185,15);
00329    all->SetBinContent(186,15);
00330    all->SetBinContent(187,14);
00331    all->SetBinContent(188,14);
00332    all->SetBinContent(189,13);
00333    all->SetBinContent(190,15);
00334    all->SetBinContent(191,8);
00335    all->SetBinContent(192,8);
00336    all->SetBinContent(193,11);
00337    all->SetBinContent(194,9);
00338    all->SetBinContent(195,11);
00339    all->SetBinContent(196,7);
00340    all->SetBinContent(197,6);
00341    all->SetBinContent(198,10);
00342    all->SetBinContent(199,9);
00343    all->SetBinContent(200,7);
00344    all->SetBinContent(201,6);
00345    all->SetBinContent(202,3);
00346    all->SetBinContent(203,4);
00347    all->SetBinContent(204,8);
00348    all->SetBinContent(205,8);
00349    all->SetBinContent(206,4);
00350    all->SetBinContent(207,6);
00351    all->SetBinContent(208,5);
00352    all->SetBinContent(209,4);
00353    all->SetBinContent(210,6);
00354    all->SetBinContent(211,9);
00355    all->SetBinContent(212,4);
00356    all->SetBinContent(213,4);
00357    all->SetBinContent(214,2);
00358    all->SetBinContent(215,2);
00359    all->SetBinContent(217,1);
00360    all->SetBinContent(228,1);
00361    all->SetBinContent(230,1);
00362    all->SetBinContent(236,1);
00363    all->SetBinContent(239,1);
00364    all->SetBinContent(272,1);
00365    all->SetBinContent(277,1);
00366    all->SetBinContent(286,1);
00367    all->SetBinContent(299,1);
00368    all->SetBinContent(309,1);
00369    all->SetBinContent(311,1);
00370    all->SetBinContent(318,1);
00371    all->SetBinContent(325,1);
00372    all->SetBinContent(329,1);
00373    all->SetBinContent(336,1);
00374    all->SetBinContent(344,1);
00375    all->SetBinContent(345,1);
00376    all->SetBinContent(346,1);
00377    all->SetBinContent(356,1);
00378    all->SetBinContent(364,1);
00379    all->SetBinContent(370,1);
00380    all->SetBinContent(397,1);
00381    all->SetBinContent(399,1);
00382    all->SetBinContent(412,1);
00383    all->SetBinContent(419,1);
00384    all->SetBinContent(431,1);
00385    all->SetBinContent(432,1);
00386    all->SetBinContent(436,1);
00387    all->SetBinContent(437,1);
00388    all->SetBinContent(438,1);
00389    all->SetBinContent(439,1);
00390    all->SetBinContent(440,2);
00391    all->SetBinContent(441,1);
00392    all->SetBinContent(443,1);
00393    all->SetBinContent(446,1);
00394    all->SetBinContent(448,2);
00395    all->SetBinContent(451,1);
00396    all->SetBinContent(453,1);
00397    all->SetBinContent(460,3);
00398    all->SetBinContent(465,1);
00399    all->SetBinContent(467,1);
00400    all->SetBinContent(475,2);
00401    all->SetBinContent(477,1);
00402    all->SetBinContent(479,1);
00403    all->SetBinContent(480,1);
00404    all->SetEntries(2904);
00405    all->SetStats(0);
00406    
00407    TF1 *landau = new TF1("landau","landau",1100,1800);
00408    landau->SetFillColor(19);
00409    landau->SetFillStyle(0);
00410 
00411    ci = TColor::GetColor("#cc6633");
00412    landau->SetLineColor(ci);
00413    landau->SetLineWidth(3);
00414    landau->SetChisquare(109.5263);
00415    landau->SetNDF(69);
00416    landau->SetParameter(0,736.9045);
00417    landau->SetParError(0,19.44123);
00418    landau->SetParLimits(0,0,0);
00419    landau->SetParameter(1,1241.799);
00420    landau->SetParError(1,1.408543);
00421    landau->SetParLimits(1,0,0);
00422    landau->SetParameter(2,32.28482);
00423    landau->SetParError(2,0.6532016);
00424    landau->SetParLimits(2,0,1069.392);
00425    all->GetListOfFunctions()->Add(landau);
00426    all->Draw("SAME");
00427    
00428    TPaveStats *ptstats = new TPaveStats(0.75,0.8,0.9,0.9,"brNDC");
00429    ptstats->SetName("stats");
00430    ptstats->SetBorderSize(1);
00431    ptstats->SetFillColor(0);
00432    ptstats->SetTextAlign(12);
00433 
00434    ci = TColor::GetColor("#0000cc");
00435    ptstats->SetTextColor(ci);
00436    TText *text = ptstats->AddText("Constant = 202.9 ");
00437    text = ptstats->AddText("Mean     =  1168 ");
00438    text = ptstats->AddText("Sigma    = 4.032 ");
00439    ptstats->SetOptStat(0);
00440    ptstats->SetOptFit(10001);
00441    ptstats->Draw();
00442    
00443    ptstats = new TPaveStats(0.75,0.7,0.9,0.8,"brNDC");
00444    ptstats->SetName("stats");
00445    ptstats->SetBorderSize(1);
00446    ptstats->SetFillColor(0);
00447    ptstats->SetTextAlign(12);
00448 
00449    ci = TColor::GetColor("#cc6633");
00450    ptstats->SetTextColor(ci);
00451    text = ptstats->AddText("Constant = 736.9 ");
00452    text = ptstats->AddText("MPV      =  1242 ");
00453    text = ptstats->AddText("Sigma    = 32.28 ");
00454    ptstats->SetOptStat(0);
00455    ptstats->SetOptFit(10001);
00456    ptstats->Draw();
00457    c0->Modified();
00458    c0->cd();
00459    c0->SetSelected(c0);
00460 }
/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/Renaud/ROOT/analog.C
