Line: 1 to 1 | ||||||||
---|---|---|---|---|---|---|---|---|
Top_Statistics | ||||||||
Line: 12 to 12 | ||||||||
| ||||||||
Changed: | ||||||||
< < |
example_xThis program calculates the cross-section and vtb. This section offers some addtional explaination and notes on this program. | |||||||
> > |
Program ChangesThis section offers some addtional explaination and notes on these programs. | |||||||
What you might need to modify within the program:
| ||||||||
Changed: | ||||||||
< < |
| |||||||
> > |
How to run the programsThe cross-section program is the most straightforward. To generate a cross-section (pretty fast):./example_x -nEns 10000Where -nEns specifies the number of ensembles (a higher number takes longer to run but is more accurate) To find vtb, just pass an arguement: ./example_x -vtb -nEns 10000The significance program is a little more complicated because it involves multiple steps. ./example_LLR_x -xs 323. -out atlasLLR_st_323.root -nEns 5000 # figure out the mean LLR value from the above run (it is in the print out), and then feed it to the following as -checkLLR ./example_LLR_x -xs 0. -out atlasLLR_st_0.root -nEns 1000000 -checkLLR -16.6482The checkLLR part isn't strictly nessesary, but you do need to generate ensembles at xs = 0 and xs = signal cross-section. The number of ensembles for the xs = signal cross-section case usually needs to be much larger than the xs = 0 case. The macro you will need to determine the actual significance (llr_plot.C, see below) will require both a xs = 0 and a xs = signal cross-section file to run. If you think you will have a large significance, you will need lots of ensembles, which will take a long time to run. In this case, you may want to generate lots of smaller sets of ensembles that use different seeds, like example_LLR_x -xs 323. -out atlasLLR_dt_323_623.root -nEns 1000000 -seed 623Once the output files are written, you can add them together (with hadd) and run them in a program to determine the significance. This program is misc/macros/llr_plot.C I have a slightly modified version of this macro, as well as a script to write lots of ensemble generation jobs for the grid. Please let me know if you would like either of these. Otherwise, I will put them on CVS after I have tried them out again. | |||||||