GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » Is there a fast way of histogramming RhoCandLists?
Is there a fast way of histogramming RhoCandLists? [message #17245] Fri, 05 September 2014 12:50 Go to next message
MartinJGaluska is currently offline  MartinJGaluska
Messages: 203
Registered: March 2010
Location: Germany
first-grade participant
From: *physik.uni-giessen.de
Dear all,

I would like to histogram multiple kinematical variables of multiple RhoCandLists within an analysis task. Instead of defining a lot of histograms or rerunning the task with slight modifications, I was wondering if there already exists a short cut.

I was thinking about a class which would only need the RhoCandList and a prefix as input and would take care of creating (using a vector for instance), filling and saving (to the output root file) the corresponding histograms itself.

Kind regards,
Martin
Re: Is there a fast way of histogramming RhoCandLists? [message #17249 is a reply to message #17245] Fri, 05 September 2014 22:18 Go to previous messageGo to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *adsl.alicedsl.de
Hi Martin,
Do you mean RhoTuple? see more details here:https:// panda-wiki.gsi.de/foswiki/bin/view/Computing/PandaRootRhoTutorial#A_4._N _45tuple_Analysis_using_RhoTuple

Best regards,
Lu
Re: Is there a fast way of histogramming RhoCandLists? [message #17250 is a reply to message #17245] Fri, 05 September 2014 22:36 Go to previous message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *pool.mediaWays.net
There is even a faster way: The "QA Tools"

See how it is applied in macro/scrut/PndScrutAnaTask:
https://subversion.gsi.de/trac/fairroot/browser/pandaroot/trunk/macro/sc rut/PndScrutAnaTask.cxx

        RhoTuple* ntp1 = new RhoTuple("ntp1", "jpsi analysis");
        if (ntp1) ntp1->GetInternalTree()->SetDirectory(gDirectory);
        PndRhoTupleQA qa(fAnalysis,fIni.P());
...
               // dump information about composite candidate tree recursively (see PndTools/AnalysisTools/PndRhoTupleQA)
                qa.qaComp("j", jpsi[j], ntp1);


You'll get many ntuples which follow s sort of naming code. First your prefix, here "j", and at the end the observable in play, such as px, py, pz, p, e, .... In between you'll find "d0" for the first daughter in your composite. This cascades down, so e.g. the momentum of the second dughter of the first daughter is labelled with "jd0d1p".

See also https://subversion.gsi.de/trac/fairroot/browser/pandaroot/trunk/PndTools /AnalysisTools/PndRhoTupleQA.h for all the automated possibilities.

Cheers
Ralf
Previous Topic: MC Truth match with PHOTOS photons
Next Topic: vertex fitting
Goto Forum:
  


Current Time: Thu Mar 28 11:46:47 CET 2024

Total time taken to generate the page: 0.00809 seconds