How to Run qSigEx

All the code and scripts to run the ensemble test presented in this note can be found via the Carleton WEB page [5]. The ensemble test needs to produce many template PDFs with 1 millions events each. Then one needs to generate 10,000 samples and preform 10,000 fits. We used the qSigEx software developed at Carleton because it allow a flexible environment and reliable interface. The maximum likelihood method of Eq. (4) is embedded in qSigEx and interfaced with MINUIT. The fitted parameters and their errors are stored in a ROOT TTree and are very easy to access. The qSigEx manual can be found on the Carleton web page [5] or on manhattan [6].

The main program used to generate the ensemble test is bias.cxx. It calls mcgen.cxx and datagen.cxx. The function mcgen.cxx generates the template PDFs; while datagen.cxx performs the ensemble test by generating many experiments. The Monte Carlo method is coded in generate.cxx where the acceptance-rejection method [3,4] is used to generate the PDFs.