Pixels:

nlines = 512 ncol = 640

Loop to create image : for (int j=0;j<nlines;j++) {//nlines for (int i=0;i<ncol;i++) {//ncol } }

use cut to see only the aperture: if(i>cutX && i<630) if(j > cutY && j <440) where cutX = 290 and cutY = 100

if ( abs (z-before) > 200) filtering noisy pixels, where before = 25 ( background) and after the first iteration, before changes to z.

Histogram: h->SetBinContent(i-cutX,j-cutY,z); max: maximum intensity and maxX and maxY correspond to the coordinates.

Gaussian fit for pbars: TF2 * fit = new TF2( "f2" , "25+[0]*exp(- ( (cos([5])*cos([5])/(2*[2]*[2]) + sin([5])*sin([5])/(2*[4]*[4]))*pow(abs(x-[1]),2) + (-sin([5]*2)/(4*[2]*[2]) +sin([5]*2)/(4*[4]*[4]))*(x-[1])*(y-[3]) + (sin([5])*sin([5])/(2*[2]*[2]) + cos([5])*cos([5])/(2*[4]*[4]))*pow((abs(y-[3])),2)))",x1,x2,y1,y2); // ellypse for pbars - works

Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2013-09-09 - AndreaGutierrez
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback