, view all tags
-- PeterWinter - 16 Feb 2009

This page will keep a description of how I use to create the final lifetime histograms from the output shrubs from ./mta....

Subdividing the data sets into groups of same magnet setting

First thing is to separate the runs into data sets according to Sara's classification scheme. One important step for that is to sort them into according groups of magnet settings (mainly +125A and -125A). That is done with ~/jobscripts/ExtractMagnetForDataSets.pl at NPL.

Creating the links

After the mta has produced the intermediate shrubs, it might be necessary to create appropriate links to the data sets. Since Brendan had the original trees in his directories on Abe in /cfs/projects/lqo/bkiburg/run11_4/treepass2/mta/, I created data set directories in /cfs/scratch/projects/lqo/pwinter/run11_4/treepass2/. Then, to create the links for each of the data sets, I have a macro in /u/ac/pwinter/jobscripts called MakeLinksForDatasets.pl. It will make the links according to the run numbers that are listed in the current $datasetfile = $HOME."/jobscripts/lists/datasets/".$dataset.".txt";

Merging the shrubs on the nodes

Transfer to NPL

Conversion to lifetime histograms

  1. In order to create histograms, I have a script on NPL at ~/jobscripts/MakeHistogramExe.pl. If you run it, it requires thre input parameters:
MakeHistogramExe.pl input_directory data_set_name prefix where
    • input_directory: Is the directory, where your shrub files are
    • data_set_name: It's the name of the current data_set
    • prefix: Is the beginning of your shrub file names before the run number, i.e. if your shrubs are Pixel_34536-36776.root, then prefix="Pixel_"
If you run this script, it will loop over all the root files with this prefix and create a HistogramCreation_$data_set_name.exe file. You might need to adapt the script since it now expects shrub file names with a range of run numbers, i.e. Pixel_34536-36776.root rather than Pixel_34536.root! You should also edit the output directory in this MakeHistogramExe.pl to match your needs.
  1. You can now launch the HistogramCreation_$data_set_name.exe files from the step before. It will simply call several times ROOT wth the macro CreateHistogramsFromTreeWithLoop.C on the according shrub files. The CreateHistogramsFromTreeWithLoop.C will mainly produce various TH3 histogams with the axis:
    • x-axis: Time dt (electron-muon)
    • y-axis: Some parameter of your choice like gondola number, nContEH, extraEL
    • z-axis: TPC subvolume number
Depending on what you want to do, you need to create other histograms with the y-axis being adapted to the parameter of interest.

(Condor job later?)

Fitting the data using the class TAllScans

Edit | Attach | Watch | Print version | History: r12 | r9 < r8 < r7 < r6 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r7 - 2009-04-16 - PeterWinter
  • Edit
  • Attach
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback