 |
 |
I response to the error message I obviously must be that the parameters in floating all parameters are not acceptable when it comes to fixed certain of those numbers and running submit.pl like so; ./submit.pl -data -all -allcut -depl -intpur 0.1 -parfile param.txt -scanmx. C.C. On Thu, 10 Aug 2006, Chukwudi Kweku Clarke wrote: > 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 >> >