Data analysis

Please find general instructions concerning the analysis with Vetra v5r1 here. This chapter focuses on specific tasks for which a few scripts have been defined to make the analysis more efficient. Two tasks have been considered so far, the slice test and the module power up. The following instructions cover both, please replace task with either "slicetest" or "modulepowerup".

Vetra

  1. Log on to the analysis farm plus on the online network.
  2. Only temporary: Setup the CMT environment by executing "source /group/velo/sw/scripts/setup_pre64bit.sh".
  3. Set up Vetra by executing "setenvVetra v5r1".
  4. Run the setup script in the Vetra installation group directory: "source $User_release_area/Vetra_v5r1/Velo/Vetra/v5r1/cmt/setup.sh"
  5. Go to the Vetra scripts directory: "cd $VETRAROOT/scripts". In here you find a couple of shell scripts to run different Vetra jobs:
    • taskDelayscan.sh: digitisation delay scan
    • taskCabletest.sh: cable test
    • taskSnapshot.sh: standard data taking without test pulses
    • taskSnapshotZS.sh: standard data taking with test pulses
  6. For the cable test and the standard data taking the scripts are executed with the following syntax:
    ./scriptName.sh filepath/filename.dat
    (e.g. ./slicetestCabletest.sh /daqarea/data/velo/slicetest/slicetest_slot_VL05_CT_tell1_207_cabletest_080208_1040.dat).
  7. Wait for Vetra to finish.
  8. The Root output file is written to the corresponding directory in /group/velo/task
    (e.g. /group/velo/slicetest/cabletest/slicetest_slot_VL05_CT_tell1_207_cabletest_080208_1040.root).

Root

If you start a new session instead of continuing with the same session that you used for VETRA, please repeat steps 1,3,4 (not 2) of the Vetra section to set up the environment properly. At the moment you have to start a separate session for Root.

Delay scan

  1. Change to the directory to /group/velo/task/delayscan
  2. Start root.
  3. Load the analysis macro via ".L ../macros/taskDelayscan.C"
  4. The macro has 4 arguments: taskDelayscan(rootFilename, tell1SerialNumber, print=true, quit=false)
    • The name of the root file has to be given in double quotes without the extension ".root".
    • The TELL1 serial number has to be provided as an integer.
    • The boolean print determines whether the plots are saved as ps/png files. The default is true.
    • The boolean quit determines whether Root is quit after the macro execution. This option is useful when the analysis is performed in batch mode via shell script. The default is false.
  5. The text file rootFilename_optimaldelay.txt contains the phase an cycle delays to be set for each PP-FPGA.

Cable test

  1. Change to the directory to /group/velo/task/cabletest
  2. Start root.
  3. Load the analysis macro via ".L ../macros/taskCabletest.C"
  4. The macro has 4 arguments: taskCabletest(rootFilename, tell1SerialNumber, print=true, quit=false)
    • The name of the root file has to be given in double quotes without the extension ".root".
    • The TELL1 serial number has to be provided as an integer.
    • The boolean print determines whether the plots are saved as ps/png files. The default is true.
    • The boolean quit determines whether Root is quit after the macro execution. This option is useful when the analysis is performed in batch mode via a shell script. The default is false.

Noise and Pedestals before and after CM suppression

  1. Change to the directory /group/velo/task/snapshot (no test pluses) or /group/velo/task/snapshotZS (test pulses)
  2. For this analysis two python scripts can be used:
  3. Invoke the script via "python -i ../macros/SensorSnapshot.py"
  4. Wait for PYTHON to start.
  5. Call the analysis macro at the command prompt (>>>) via "SensorSnapshot(rootFilename, [tell1SerialNumber])"
    • The name of the root file has to be given in double quotes with or without the extension ".root". Both versions work.
    • The TELL1 serial number will be automatically detected from the contents of the file if the root file contains data from only one TELL1. If more than one TELL1 is present please provide the TELL1 serial number as an integer.
    • If you want to analyze data from all the tell1s in one go then invoke the script "SensorSnapshotSara.py".

-- OlafBehrendt - 24 Apr 2008

Edit | Attach | Watch | Print version | History: r5 < r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r5 - 2008-07-02 - StefanoDeCapua
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    LHCb All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback