Page History

...

1. Reconstruct data with GBL (to this purpose, you need to include the GblDriver and GblOutputDriver in your steering file).

   1. For Montecarlo data:

      1. run the readoout with the steering file: HPS2014ReadoutNoPileup.lcsim
      2. run the recon with the steering file: HPS2014OfflineNoPileUp.lcsim
   2. For real data use as steering file: /org/hps/steering/recon/EngineeringRun2015FullReconGbl2.lcsim
2. Check that at the end of reconstruction a out.gbl ascii file (or, named as you decided in the GblOutputDriver) is produced.
3. Remember that by default the geometry is taken from the database. If you want to force the use of your own geometry, you must provide it in the compact.xml file in a given detector. For MonteCarlo data, set the run number to zero during readout and reconstruction. This is done adding the flag -Drun=0 when running the readout. For real data, use the -DdisableSvtAlignmentConstants flag.. IMPORTANT: remember to re-compile hps-java before running each time you change the compact.xml file! (this is the most common error).
4. The out.gblfile is read by a python procedure. You must download with git the current version of the software from the github repository as descrived in the following. This will create a hps-gbl directory. After having configured your account and username for git usage, issue the following commands:
   • git clone https://github.com/perhansson/hps-gbl
   • cd hps-gbl
   • git submodule init
   • git submodule update
   • git pull
   (the second and third command need to be issued just upon installation, and they are needed since you will have to download something else from github -see step 1 of the next subsection- and some directories with needed utils are shared in the github repository).
5. remember to install the GBL software, if you already haven't it. In a directory parallel to hps-gbl download the GBL software using svn:
   • svn co http://svnsrv.desy.de/public/GeneralBrokenLines/tags/v01-16-0204 GeneralBrokenLines
   (or check the newest release, and get it). To compile it:
   • cd GeneralBrokenLines/cpp
   • mkdir build; cd build
   • cmake ../
   • make install
   • make doc (if you want it)
   Note: if you have installed the latest cmake version, it probably won't compile. You must prevent the compilation to search for C++11 support (the default for newest cmake). To do this, you have to set as compilation flag -std=c++0xadding it to the c++ compilation line. Either you do it in the cmake configuration files, or (quickest) you add by hand this flag at the end of the CXX_FLAGS line, in the two files:
   • GeneralBrokenLines/cpp/build/CMakeFiles/GBL.dir/flags.make
   • GeneralBrokenLines/cpp/build/examples/CMakeFiles/GBLpp.dir/flags.make
6. the gbl python procedure reads the out.gbl file and prepares the binary read by Millepede. You mustrun python from the hps-gbl directory. This is the shortest syntax (-h shows all possible options):
   • cd hps-gbl
   • python gbltst-hps.py [-h] out.gbl
   A heap of pdf files are produced containing plots of several quantities for top/bottom halves, with long names that should be self-explaining (but at the moment they are not and they are too long, this must be improved). You will also file a .root file containing the single root histograms, and a .ps file containing a summary of the plots ready to be printed.
   The file gbltst-hps.py contains the instructions to extract the useful information on tracks and hits from the ascii file and write the input file for Millepede. If you want to add/modify some of the output plots/histograms, you have to modify both the gbl_plots.py file (in which they have to be booked) and the gbltst-hps.py file, in which they have to be filled.
   Note: root must be compiled including the python support, otherwise python stops with an error complaining about root libraries missing. A good idea is to put in your profile and instruction to run automatically $ROOTSYS/bin/thisroot.(c)sh, which provides the correct root-python environment and libraries for your system.
   At the end of python run, you should also find a MilleBinaryISN.dat file (the name could slightly change), which is the input file to be read by Millepede.

...