cg Version

x :x = 1 NNXxi (1)i=1 xi :s 2 = 1N ; 1NX(xi ; x) 2 ' x 2 ; x 2 (x 2 = 1 Ni=1NXi=1x 2 i ): (2) x s 2 x = 1 N s2 ' 1 N [x2 ; x 2 ] (3) sx ' [ 1 N (x2 ; x 2 )] 1=2 (4)4.9 : Step 9! -----------------------------------! Calculate average and its deviation! -----------------------------------! ---------------! Peak efficiency! ---------------avpe = pefs/ncountpef2s=pef2s/ncountsigpe=dsqrt((pef2s-avpe*avpe)/ncount)avpe = avpe*100.0sigpe = sigpe*100.0write(6,350) avpe,sigpe350 FORMAT(' Peak efficiency =',G11.4,'+-',G9.2,' %')! ----------------! Total efficiency! ----------------avte = tefs/ncounttef2s = tef2s/ncountsigte = dsqrt((tef2s-avte*avte)/ncount)avte = avte*100.0sigte = sigte*100.0write(6,360) avte,sigte360 FORMAT(' Total efficiency =',G11.4,'+-',G9.2,' %')! --------------------------! Pulse height distribution! --------------------------write(6,370)370 FORMAT(/' Pulse height distribution ')do ie=1,50elow=deltae*(ie-1)eup=deltae*ieavph = phs(ie)/ncountph2s(ie)=ph2s(ie)/ncount13