June 09, 2004

White noise check in a flat background::[Misc]

The white noise can be a flat or other distribution, here I simulate a flat white noise + a gaussian signal.

Using the scan method and unbined maximun likelihood fit to find out the max significance on a pure white noise.
</p> <p>#!/bin/sh</p> <p>NUMOFDATA=1000<br /> TIMES=10<br /> COUNT=0<br /> STEP=5</p> <p><br /> echo 10000 >log.data<br /> echo 10000 >log.img</p> <p>GenShell() {</p> <p><br /> echo "#!/bin/sh"<br /> echo <br /> echo "root <<EOF"<br /> echo <br /> echo <br /> echo 'gSystem->Load("libRooFitCore.so") ;'<br /> echo 'gSystem->Load("libRooFitModels.so") ;'<br /> echo 'RooRealVar mass("mass","mass",0 , 100 , "Gev/c^{2}") ; '<br /> echo 'RooRealVar mean("#mu","#mu",0, 5,"Gev/c^{2}") ;'<br /> echo 'RooRealVar sigma("#sigma","#sigma",3, "Gev/c^{2}") ;'<br /> echo 'RooRealVar gausFrac("Gaussian Fraction","Gaussian Fraction",0 , 1 ) ;'<br /> echo 'RooGaussian g("g", "Gaussian", mass, mean , sigma );'<br /> echo 'RooPolynomial p("p", "Polynominal", mass);'<br /> echo 'RooAddPdf model("model", "Gaussian + Polynominal",RooArgList(g,p),gausFrac);'<br /> echo 'gausFrac.setVal(0.0);'</p> <p><br /> echo <br /> echo </p> <p>echo 10000 >log.img</p> <p>while [ $COUNT -lt $TIMES ]<br /> do</p> <p>echo "RooDataSet *data=model.generate(mass, $NUMOFDATA);"<br /> echo 'data.write("generate.data");'<br /> echo <br /> echo<br /> MASS=0<br /> MAS2=0<br /> while [ $MAS2 -le 95 ]<br /> do <br /> MAS2=`expr $MASS + $STEP`<br /> #echo "mean.setVal( $MASS, $MAS2 );"<br /> echo "mean.setFitRange( $MASS, $MAS2 );"<br /> echo 'RooFitResult *fs = model.fitTo(*data, "rq");'<br /> echo 'fs->Print("v");'</p> <p>echo 'RooPlot *frame=mass.frame();'<br /> echo 'data->plotOn(frame);'<br /> echo 'model.paramOn(frame, data);'<br /> echo 'model.plotOn(frame);'<br /> echo 'model.plotOn(frame,Components(p), LineStyle(kDashed));'<br /> echo 'frame->Draw();'<br /> echo 'TCanvas *c = new TCanvas;'<br /> echo 'TH1F *h = new TH1F( );'<br /> echo 'frame->Draw();'<br /> echo 'TImage *img = TImage::Create();'<br /> echo 'img->FromPad(c);'<br /> echo 'img->WriteImage("canvas.png");'<br /> echo 'delete h;'<br /> echo 'delete img;'<br /> echo 'delete c;'<br /> echo '.!./move.sh img'<br /> MASS=`expr $MASS + $STEP`<br /> echo<br /> echo<br /> done<br /> echo</p> <p>echo '.!./move.sh data'<br /> COUNT=`expr $COUNT + 1`<br /> echo<br /> echo</p> <p>done</p> <p>echo '.q'<br /> echo <br /> echo <br /> echo "EOF"</p> <p>}</p> <p>GenShell > ./poloynamial.root</p> <p>chmod 755 poloynamial.root</p> <p>./poloynamial.root > output.data</p> <p>mv 100* data<br />

Trackback

Comments

Post a comment