QIE10 Chip Test

In order to study the different variables from the QIE10 chip test at Fermilab, this twiki explains how to use the code for plotting the different variables.

1. Get the code for plotting your histograms

After you are logged in a lpc machine please follow this instructions:

- This code will use pyROOT. Create a working area and set a CMSSW environment (e.g. CMSSW_5_3_16)

mkdir YOURWORKINGAREA
cd YOURWORKINGAREA
### If you are using csh shell
setenv SCRAM_ARCH slc6_amd64_gcc481 
### If you are using Bash shell
export SCRAM_ARCH=slc6_amd64_gcc481 
cmsrel CMSSW_7_3_0_pre1
cd CMSSW_7_3_0_pre1/src
cmsenv

- Get the code from GitHub :

git clone git@github.com:danielguerrero/QIE10.git
cd QIE10

2. Make the TH1F and TH2F plots

The script runPlotter.py will generate the plots that you want.This plotter can help you to make plots with these features:

  • Plot the chip variables (TH1F )
    • Add the label to the axis and title
    • Number of bins
    • Stats box showing: Mean, RMS, Entries, Overflow and Underflow
    • Add a caption
  • Plot the chip variables (TH2F)
    • Compare two chip variables
    • Add the label to the axis and title
    • Number of bins
    • Add a caption
    • Generate a colz, scat plot by editing the code

The code will ask you what kind of plot you want like this:

Would you like to make TH1D plots (Yes/No)?
  Would you like to add the cut bars to the plots (Yes/No)? 
Would you like to make TH2D plots (Yes/No)? 

In order to get your plots and add the labels that you want, you need to create some text files and write the information in columns as it is explained below:

1. qie10_histograms.txt: "Histogramname"

2. qie10_labels.txt: "LabelinXaxis LabelinYaxis "

3. qie10_cuts.txt: "ValueLowCut ValueHiCut ". Note: If you want to add some bars of the cuts from the chip, you have to create this text file.

4. qie10_bins.txt:"Nbins"

In the case of TH2F plots you need also to create a text files:

5. qie10_2dhistograms.txt: "HistonameXaxis HistonameYaxis"

After that, you just have to run the script over the ntuple in this way (file.root is the ntuple you get from the chip):

python runPlotter.py file.root

The output are some folders with the plots (you can edit the name of the folder in the code):

  • qie10plots (TH1F)
  • qieplots2dscat (TH2F)
Here is some of examples of how your plots will look:

- TH1F variable plot:

Volt_Isetp.png

- TH2F scatter plot

Volt_InSig_VS_Volt_InRef_scat.png

- TH2F colz plot

Volt_InSig_VS_Volt_InRef_colz.png

3. Make Comparison Histogram

Let's say you want to compare chips based on cut on one of the variable (e.g. Volt_ref <0.39). The runPlottersplit.py can make the comparison with these features:

  • A plot which will compare the shape of the distributions
  • Add poisson errors
  • Make a hatched band comparing plot in aiming to see if there is a correlation between the two groups (e.g. Volt_ref <0.39 and Volt_ref >0.39)
  • Set labels and title
  • Add a caption
In order to get your plots and add the labels that you want, you need to create some text files and write the information in columns as it is explained below:

1. qie10_histograms.txt: "Histogramname"

2. qie10_labels.txt: "LabelinXaxis LabelinYaxis "

3. qie10_bins.txt:"Nbins"

After that, you just have to run the script over the ntuple in this way (file.root is the ntuple you get from the chip):

python runPlottersplit.py file.root

The output is folder with the plots (you can edit the name of the folder in the code): qie10splits

Here is an example of a plot using the plotter:

Volt_Isetp_split.png

Topic attachments
I Attachment History Action Size Date Who Comment
PDFpdf Volt_InSig_VS_Volt_InRef_colz.pdf r1 manage 14.8 K 2015-02-21 - 18:11 DanielGuerrero  
PNGpng Volt_InSig_VS_Volt_InRef_colz.png r1 manage 62.2 K 2015-02-21 - 18:57 DanielGuerrero  
PDFpdf Volt_InSig_VS_Volt_InRef_scat.pdf r1 manage 23.2 K 2015-02-21 - 18:11 DanielGuerrero  
PNGpng Volt_InSig_VS_Volt_InRef_scat.png r1 manage 68.3 K 2015-02-21 - 18:57 DanielGuerrero  
PDFpdf Volt_Isetp.pdf r1 manage 13.7 K 2015-02-21 - 18:11 DanielGuerrero  
PNGpng Volt_Isetp.png r1 manage 54.4 K 2015-02-21 - 18:57 DanielGuerrero  
PDFpdf Volt_Isetp_split.pdf r1 manage 17.9 K 2015-02-21 - 18:11 DanielGuerrero  
PNGpng Volt_Isetp_split.png r1 manage 62.9 K 2015-02-21 - 18:57 DanielGuerrero  
Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2015-02-21 - DanielGuerrero
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

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.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback