 |
 |
Hi Antonio, I seem to be having problems when it come to getting the correct parameters for data depleted. I'm using Concezio's VVF-tag (cb-060801). I have been able to get numbers for the MC enriched and depleted, which you can have a look at on my webpage (http://www.slac.stanford.edu/~clarke1/protected/web/IP.html). The problem is, if I use the same instructions for data as I do for MC the following error is observed: RooRealVar::setRange(mes) new range named 'bmesint' created with bounds [5.27,5.292] USING CUTS mes> 5.22 && mx > 0.000000 && mx < 1.550000 && q2 > -999.000000 && q2 < 300.000000 && lepYaSe>0 && ksele != 0 && intpur > 0.100000 Number of entries in dataset: 1244 scalefactor= 0.284083 WARNING ar NOT FOUND in param.txt Error: illegal pointer to class object argpar 0x0 791 FILE:fittest.C LINE:557 *** Interpreter error recovered *** What I have done is ./submit.pl -data -bkg -lepcut -depl -dumppar -intpur 0.1 -mxl -30 -mxh 30 putting these restant parameters into a .txt file then ./submit.pl -data -sig -lepcut -depl -dumppar -intpur 0.1 -mxl -30 -mxh 30 putting these restant parameters into the same .txt file BTW the above both produce the same results, i.e., the flags, -sig and -bkg in conjuction with -data have no effect. Then using the .txt file like so ./submit.pl -data -all -allcut -depl -intpur 0.1 -parfile param.txt -scanmx and then I get the above error. Curiously, I have been able to get results for data depleted using the same paramters file (paramIP.txt) that I used for the MC studies and I was also able to get results using the parameters_DATA.txt file (having set fixed the values you mentioned below). Cheers, C.C. On Wed, 26 Jul 2006, Antonio Petrella wrote: > Hi Chukwudi > > so the goal is to compute the correction factor for different integrated > lumiosity cut: the formula is > > (S/P)_MC enriched > ------------------ * (S/P)_data depleted. > (S/P)_MC depleted > > To compute S/P you can use the macro in fittest directory. > > Let's take MC enriched as example and let's suppose the integrated > purity cut is 0.1 > > The logical steps are: > > - find bkg pdf parameter fitting only background MC on entire sample. > Let all parameters floating... the command would be: > > ./submit.pl -MC -bkg -lepcut -enr -dumppar -intpur 0.1 -mxl -30 -mxh 30 > > (with -dumppar you will have a dump of the results of the fit on a txt > file). > > - control that the fit is good ie. the error matrix is ok and fitted > values are all ok. > > - find signal pdf parameters fitting only signal MC on entire sample. > Let all paraemters floating. The command is: > > ./submit.pl -MC -sig -lepcut -enr -dumppar -intpur 0.1 -mxl -30 -mxh 30 > > - control that the fit is good ie. the error matrix is ok and fitted > values are all ok. > > Now take the two txt files with dumped parameters and create a new one > by copying the signal stuff from the signal dumped parameter file and > background related stuff from the background dumped parameter file. > You will use this new file as input parameters for scan in bin of mx. > Let's assume that you call this file "parameterforscan.txt" > set to fixed (1) in this file the following parameters: > > Thosig_r, Thosig_xc, sigma_r2, Thosig_n, Thosig_alpha, mean_cb, > sigma_cb, alpha_cb, n_cb, cbendpoint. > > Now you are ready to run the scan (allcut now): > > ./submit.pl -MC -all -allcut -enr -intpur 0.1 -parfile > parameterforscan.txt -mxscan > > - verify that the fits are good (look the plots and also the output of > the fits). > - look into the txt file with results: you will have also the S/P +- > error printed in this file (look also README, section 5- OUTPUT FILE). > > So you have S/P for MC enriched. Do the same procedure for MC depleted > and data depleted. At the end you can compute the correction factor. > > One more comment: instead of fitting MC all compontents, you can take > the P values doing a scan on background only > (./submit.pl -MC -bkg -allcut -enr -intpur 0.1 -parfile > parameterforscan.txt -mxscan), > and the S values counting events in the signal dataset. > > You can find this number on the shell output of the fit when running > > ./submit.pl -MC -sig -allcut -enr -intpur 0.1 -parfile > parameterforsca.txt -mxscan). In this case you avoid fitting all MC. > > I hope everything is clear. > > Let me know, if you have problems I can read emails from time to time... > > Bye > Antonio >