TWiki> Main Web>TWikiUsers>ChengguangZhu>ZhuUseROOT (revision 24)EditAttachPDF
HistgramDrawOptions

Can't use pointer returned from Root Class Object

A Root class pointer defined in a function, but when it is returned out of the function, it can't be used any more

solution: Becuase the memory of the Root pointers are all managed by Root, when the function finished, the memory will be released by Root ( c++ don't do this) use pointer->SetDirectory(0) to let Root not to manage the momery of this pointer. or TH1::AddDirectory(kFALSE) to let Root not to manage memoty fo any pointers of TH1.

search memory problem

http://www.veryword.com/2009/08/valgrind_tutorial.html

valgrind --tool=memcheck program

to find which lib produce the following error.

our example is 1. we copy a .h file from the A package to B package. 2. the .h define some class and objects 3. if load A and B lib, the object defined in .h are free-ed 2 times.

problem:

*** glibc detected *** /afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/bin/root.exe: double free or corruption (!prev): 0x0000000006b91a40 ***
======= Backtrace: =========
/lib64/libc.so.6[0x2aff5d850d7f]
/lib64/libc.so.6(cfree+0x4b)[0x2aff5d8511db]
/lib64/libc.so.6(exit+0xe5)[0x2aff5d813515]
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/lib/libCore.so(_ZN11TUnixSystem4ExitEib+0x85)[0x2aff5b4dc9c5]
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/lib/libCore.so(_ZN12TApplication11ProcessLineEPKcbPi+0x2b6)[0x2aff5b3fd5e6]
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/lib/libRint.so(_ZN5TRint15HandleTermInputEv+0x1a6)[0x2aff5c9d7a76]
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/lib/libCore.so(_ZN11TUnixSystem16CheckDescriptorsEv+0x174)[0x2aff5b4df9f4]
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/lib/libCore.so(_ZN11TUnixSystem16DispatchOneEventEb+0xe0)[0x2aff5b4e0010]
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/lib/libCore.so(_ZN7TSystem9InnerLoopEv+0x16)[0x2aff5b45cfb6]
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.28.00g/x86_64-slc5-gcc43-opt/root/lib/libCore.so(_ZN7TSystem3RunEv+0x7c)[0x2aff5b45f07c]

Problem in using TTree::SetBranchAddress

When use TTree::SetBranchAddress to link variables to some of the branches (not all branches), some time, the link is not working as expected.

for example, SetBranchAddress("branch_a",&var_a) , but after tree->GetEntry(i), var_a may in fact get the value from branch_b, In this case branch_b is not linked to any variable.

To solve the prolbem:

to switch off all other not-used branches with SetBranchStatus( "branch_b", 0);

TChain read only the first file in the chain

When use TChain to chain several files and use "maked class" to read the events, it happened that only the first file is read once and once:
for( int i=0; i<tchain.GetEntries(); i++){
  tchain.GetEntry(i);
  ...
  ...
}

The correct should be 
for( int i=0; i<tchain.GetEntries(); i++){
  tchain.GetTree()->GetEntry(i);
  ...
  ...
}

ERRORs and reason/solutions

  • error message
Error: Symbol #include is not defined in current scope  (tmpfile):1:
Error: Symbol exception is not defined in current scope  (tmpfile):1:
Syntax Error: #include <exception> (tmpfile):1:
Error: Symbol G__exception is not defined in current scope  (tmpfile):1:
Error: type G__exception not defined FILE:(tmpfile) LINE:1
*** Interpreter error recovered ***

Reason and Solution: Memory leakage problem. When the leakage is getting big enough it produces this error. In the case pointer is defined and memory is allocated but didn't deleted after using.

  • "[ROOT] Load Error: Failed to load Dynamic link library"
 This usually happens when the compiler version used to build ROOT is not the
same as the compiler available in the path (i.e. bash> g++ -v).

  • Statement too long
reloading /atlas/czhu/myAnalysisCode/my_method/./method2_a.cxx  0
Limitation: Statement too long method2_a.cxx:197:
cint: Security mode 0x7:0x2 Error: Too many '}' method2_a.cxx:208: 

reason: has the following statement, just in the file method2_a.cxx (not inside any function)
for(int i=0;i<10;i++)pa_mceffi[i]=0.;
  • Unexpected EOF, need to use +P or -p option
reloading /atlas/czhu/myAnalysisCode/notes/notes_ele25/./T7CSCnote.cxx  0
Error: Unexpected EOF G__fignorestream():3 T7CSCnote.cxx:4159:
Advice: You may need to use +P or -p option

reason:      ""  or {}   don't  exist in pairs

  • error message as below:
root [1] CheckObjects(1000015)
open new data...$HOME/analysisResults/5205_12605_jcali_atlfast/ANA.ntuple.root
open the 3nd data...$HOME/analysisResults/5205_12605_jcali_fullsim_truthHardProc/ANA.ntuple.root

 *** Break *** segmentation violation
 Generating stack trace...
In happens, when program runs switch{case break; .....}. A very long double vector is defined in one of the (case break), which make all the (case, break) after the definition can't be run, but the (case break) before this definition can be run correctly.

  • error incomplete type used in nested name specifier

g++ gives this message if you've forward-declared a type, like this

class MyClass;

and then you try and access one of its members, like maybe:

MyClass::doSomething()

g++ is complaining that it hasn't seen the body of class MyClass yet, so it has no way to know what MyClass::doSomething is.

(In C++ jargon, an "incomplete type" is a type that's been forward-declared but not yet defined.)

You can often fix this by adding the right #include at the top of the source file.

When reading a root file produced with root5.26/00d with root5.26/00, this error happen

The 5.26/00d is none year later update of 5.26/00

  • error in using proofError message:

Info in <TProofLite::MarkBad>:  events  |>...................| 0.00 %
 +++ Message from master at lxplus206.cern.ch : marking 0.0-lxplus206.cern.ch-1285753103-11416:-1 (0.0) as bad
 +++ Reason: undefined message in TProof::CollectInputFrom(...)


 +++ Most likely your code crashed
 +++ Please check the session logs for error messages either using
 +++ the 'Show logs' button or executing
 +++
 +++ root [] TProof::Mgr("lxplus206.cern.ch")->GetSessionLogs()->Display("*")


Error in <TPacketizerAdaptive::SplitPerHost>: The input list contains no elements

Solution: because m_goodRunListXmlFile doesn't exist.

1888    TDOMParser xmlparser;
1889 #if ROOT_VERSION_CODE >= ROOT_VERSION( 5, 17, 4 )
1890    xmlparser.SetReplaceEntities( kTRUE );
1891 #endif
1892    Int_t parseError = xmlparser.ParseFile( m_goodRunListXmlFile.c_str() );
1893    if( parseError ) {
1894 
1895       cout << "Loading of xml document \"" << m_goodRunListXmlFile
1896             << "\" failed";

Example of using TGraphErrors

    double sys_bt1    [5]={50, 55,  60,  65,  70 };     
    double sys_A_bt   [5]={0.718301,  0.861586,   0.855908,  0.907561,  0.858515  };
    double sys_A_bt_e [5]={0.107668,  0.0962074, 0.0879854,  0.0841657, 0.0797709 };

    TGraphErrors* sysAbt =new TGraphErrors(5, sys_bt1, sys_A_bt, 0, sys_A_bt_e );    
    sysAbt->SetTitle("A vs b tag efficiency");
    sysAbt->SetLineWidth(2);
    sysAbt->SetMarkerSize(1.1);
    sysAbt->GetXaxis()->SetTitle(" b tag efficiency (%)");
    sysAbt->GetYaxis()->SetTitle("A");    
    sysAbt->Draw("AP");
    sysAbt->Fit("pol1");

Example to creat 2D histgram of variable bins

  • This example produce 2 2D histgram, divide one by the other and draw.
const  int    mptbins = 8;
const  int    metabins = 7;  
const Float_t mxbin[mptbins+1]={20.,25,30,35,40,50,70,120.,300};
const Float_t mybin[metabins+1]={-2.5,-1.35,-1.15,-0.1,0.1,1.15,1.35,2.5};

plot_Efficiency2D(ntuple_tt_3_d, "m_reco_effi_pteta_d", weight_3*data*semimu*"selmuon_true==1 && mu_true_found_rec>0 ",  weight_3*semimu*"selmuon_true==1",  "ptmuon_true", mptbins, mxbin, 0.001, "etamuon_true",metabins, mybin,1);

void plot_Efficiency2D(TTree* tree,char* hist, TCut cutup, TCut cutdw,  char* varntuplex, int binx, const Float_t* xbin ,float scalex=1, char* varntupley, int biny, const Float_t* ybin, float scaley=1){
   tmpupref=new TH2F("tmpupref", "tmpupref", binx,xbin, biny,ybin);
   GetHistFromNtuple2D(tree, cutup, "tmpupref",varntuplex,scalex, varntupley,scaley);
   TH2F* tmpdwref=new TH2F("tmpdwref", "tmpdwref", binx,xbin, biny,ybin);   
   GetHistFromNtuple2D(tree, cutdw, "tmpdwref",varntuplex, scalex, varntupley,scaley);

   tmpdwref->Sumw2();
   tmpupref->Divide(tmpdwref);        
   tmpupref->Clone(hist);
   tmpupref->DrawCopy("lego")
   tmpupref->Delete();
   tmpdwref->Delete();
}
void GetHistFromNtuple2D( TTree* tree, TCut cut, char* histname, char* varntuplex,float scalex=1, char* varntupley, float scaley=1){
      char var[1000];
      sprintf( var , "(%s*%f):(%s*%f)>>%s", varntupley,scaley,varntuplex,scalex,histname);
      tree->Draw(var, cut,"goff");  
}

Options for drawing histgram

  • use gStyle->SetOptFit
 The type of information about fit parameters printed in the histogram
 statistics box can be selected via the parameter mode.
  The parameter mode can be = pcev  (default = 0111)
    p = 1;  print Probability
    c = 1;  print Chisquare/Number of degress of freedom
    e = 1;  print errors (if e=1, v must be 1)
    v = 1;  print name/values of parameters
  Example: gStyle->SetOptFit(1011);
           print fit probability, parameter names/values and errors.
    When "v"=1 is specified, only the non-fixed parameters are shown.
    When "v"=2 all parameters are shown.

  Note: gStyle->SetOptFit(1) means "default value", so it is equivalent to
        gStyle->SetOptFit(111)
  • Use gStyle->SetOptStat
 The type of information printed in the histogram statistics box
  can be selected via the parameter mode.
  The parameter mode can be = ksiourmen  (default = 000001111)
    k = 1;  kurtosis printed
    k = 2;  kurtosis and kurtosis error printed
    s = 1;  skewness printed
    s = 2;  skewness and skewness error printed
    i = 1;  integral of bins printed
    o = 1;  number of overflows printed
    u = 1;  number of underflows printed
    r = 1;  rms printed
    r = 2;  rms and rms error printed
    m = 1;  mean value printed
    m = 2;  mean and mean error values printed
    e = 1;  number of entries printed
    n = 1;  name of histogram is printed
  Example: gStyle->SetOptStat(11);
           print only name of histogram and number of entries.
           gStyle->SetOptStat(1101);  displays the name of histogram, mean value and RMS.
  WARNING: never call SetOptStat(000111); but SetOptStat(1111), 0001111 will
          be taken as an octal number !!
  WARNING: SetOptStat(1) is taken as SetOptStat(1111) (for back compatibility
           with older versions. If you want to print only the name of the histogram
           call SetOptStat(1000000001).
  NOTE that in case of 2-D histograms, when selecting just underflow (10000)
        or overflow (100000), the stats box will show all combinations
        of underflow/overflows and not just one single number!

use roofit and FFT with root


  do not have root access of the machine to make a private install area for FFTW, e.g. /home/myself/fftw

  1.  www.fftw.org for fft downloading
  2.  Untar fftw-XXX.tar.gz in /tmp, cd into the untarred directory
       ./configure --prefix=/home/myself/fftw CFLAGS=fPIC CXXFLAGS=fPIC
       ' followed by 'make install'.
       Substitute /home/myself/fftw with a directory of your choice. This
       procedure will install FFTW in the location designated by you

   3. in ROOT directory 
      Run 
      configure --enable-fftw3 --with-fftw3-incdir=/home/myself/fftw/include --with-fftw3-libdir=/home/myself/fftw/lib  --enable-roofit
      Then run 'make'

xmanager can't farword plot to terminal

solution run xmanager-braodcast once to fix it.

-- ChengguangZhu - 29 Nov 2007

Edit | Attach | Watch | Print version | History: r28 | r26 < r25 < r24 < r23 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r24 - 2014-02-20 - ChengguangZhu
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2021 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