The following cuts were used:
TCut cut1 = Form("(B.ts.ps.a_p[%d]<50&&B.ts.ps.a_p[%d]<50)",n-1,n+1); TCut cut3 = Form("(B.sum.rt[%d]>100||B.sum.rt[%d]>100)",n1-1,n1); where n is preshower block number and n1=n%27 is coresponding total sum number.
The function used: formula = "exp([0]+[1]*x/[7])"; formula += "+[2]/sqrt(2*pi)/[4]*exp(-0.5*(x/[7]-[3])^2/[4]^2)"; formula += "+[5]*(x/[7])^1.3*exp(-([6]*x/[7])^2.3)"; TF1* fg = new TF1("fg",formula.Data(),150,4000);
The block on the edges are skiped. Block numbering starts from 0 to 53.
1 221.183 28 222.079 2 196.232 29 145.667 3 210.284 30 145.667 4 172.74 31 203.346 5 209.553 32 227.223 6 196.118 33 216.358 7 213.047 34 203.491 8 202.888 35 235.032 9 201.771 36 220.53 10 205.905 37 221.613 11 194.164 38 210.48 12 203.008 39 203.186 13 208.744 40 218.387 14 202.429 41 215.046 15 199.352 42 220.931 16 193.878 43 214.435 17 194.417 44 194.328 18 203.895 45 200.38 19 191.366 46 253.562 20 226.572 47 226.011 21 197.163 48 187.788 22 194.089 49 175.215 23 203.69 50 183.176 24 226.388 51 249.376 25 216.508 52 145.667
Scripts could be found in the directory:
adaq@adaql?:~/e02013/COSMICS/nerses/BBPreShowerCalib/