<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Dear all,<br>
I write you again concerning BDX DSTs, following-up my previous
email (you find it below).<br>
<br>
The first version of BDX DSXTs has been implemented, and is already
working for the data being measured in Catania now with the second
version of the prototype. The strategy that I followed - after some
discussion with other people - is the second one outlined in my
previous email.<br>
Basically:<br>
<br>
1) Each event is saved in a Root TTree as a TEvent class. <br>
1A) The concept of "event" is clear for the current version of the
prototype (1 event == 1 trigger), or for simulations (1 event == 1
primary particle). It will be maybe different if we will implement a
triggerless, pipeline DAQ system, but for the moment I neglect this
and focus on the two previous cases.<br>
<br>
2) The TEvent is, basically, a collection of information (it is very
similar to the event structure we have in HPS - I mutuated it from
there). You can find the relevant code in the BDX distribution,
src/libraries/EventBuilder/TEvent.h. The relevant members are:<br>
<br>
.... <br>
<br>
private:<br>
<br>
TEventHeader *m_eventHeader;<br>
vector<TClonesArray*> m_collections;<br>
....<br>
<br>
In other words, the TEvent contains:<br>
<br>
-> The Event Header (run number, event number, event type, event
time, trigger words)<br>
-> A vector of TClonesArray: each TClonesArray is a collection of
"objects" (or "hits"), of a given type. This means that the event
itself is "system-independent", the same structure is good for real
data (prototype) / simulations / .. , provided that the event
builder saves in the TEvent the relevant collections.<br>
<br>
I send you an example macro (and a datafile) showing how to open
this kind of DST, understand what is inside the TEvent, and read it.
<br>
Clearly, the are many other ways to read the DSTs. Here, I am using
a root TSelector, and a root macro. One can use a standalone C++
code, and use a normal loop instead of a TSelector.<br>
<br>
Pre-requirements:<br>
<br>
1) BDX software installed on a non - MAC machine (still need to
figure out how to build root dictionaries on a MAC)<br>
1A) The env var BDXRECO_ROOT pointing to the installation folder<br>
2) Root6 installed (jlab common software, version 2.0, is now using
root 6. This is why I am using it for DSTs).<br>
<br>
Files:<br>
1) A root file (out.root) containing the TTree with events (EventDST
tree). <br>
1A) The root file has been created with 4000 events from prototype
data. Collections saved are calibrated hits from calorimeter, int
veto, and Ext Veto.<br>
1B) There're also two other TTress (RunInfo and EventHeader). Ignore
them.<br>
2) Two files contaning the TSelector-derived class
(BDXDSTSelector.cc and BDXDSTSelector.h)<br>
3) A macro file to run the analysis (doAna.C)<br>
<br>
Instructions:<br>
1) Save all files to a folder of your choice.<br>
2) Make sure BDXRECO_ROOT env var is set, pointing to the folder
where the BDX software is installed<br>
3) Open root and issue these commands:<br>
* .L doAna.C<br>
* doAna("out.root")<br>
<br>
The macro is quite simple. It is checking that the event has a
collection named "CalorimeterHits", made of objects of class
"CalorimeterHit". If so, it gets this collection, and loops over the
CalorimeterHit objects. If it founds a hit with sector=0, X=1, Y=1,
it fills a histogram. <br>
<br>
<br>
<br>
What still has to be done (somehow in order of priorities):<br>
<br>
1) Have root dictionaries also built on MAC (need to access to a MAC
machine with root6 installed).<br>
2) Decide what to save in DSTs for Catania measurements. Currently,
for each event 3 collections are saved: Calibrated hits for
calorimeter, IntVeto, ExtVeto.<br>
2A) Document the Event structure for Catania measurements once
defined.<br>
3) Write the event builder for simulations (different collections
will be saved in this case), both for prototype simulations and full
detector simulations.<br>
<br>
<br>
<br>
Technical details:<br>
1) The DST idea here shown makes heavy use of ROOT reflectivity
features.<br>
2) The DST reader, as shown here, is proof compatible (provided one
takes care in setting to 0 all the pointers to objects that will be
crated in the SlaveBegin method)<br>
3) Although not easily (and absolutely not recommended), one can
interact with ROOT DSTs also from root command line. The following
set of instructions should produce the same histogram as in the
attached macro<br>
<br>
gSystem->Load("$BDXRECO_ROOT/lib/libbdxReco.so")<br>
TFile f("out.root")<br>
TH1D hSimple("hSimple","hSimple",400,-0.1,400)<br>
EventDST->Scan("m_collections.GetName()") -> From
this command I see the collection with calorimeter hits is always
the first (istance == 0)<br>
EventDST->Draw("((CalorimeterHit*)(m_collections[0].At(2)))->E>>hSimple")
-> The hit #2 is the one I am interested in<br>
<br>
Bests<br>
<br>
Andrea<br>
<br>
<br>
<br>
<br>
<blockquote
cite="mid:aefe7a9d-b4ed-9e67-3626-76a06412eacf@ge.infn.it"
type="cite">
<div class="moz-forward-container"> <br>
-------- Forwarded Message --------
<table class="moz-email-headers-table" border="0"
cellpadding="0" cellspacing="0">
<tbody>
<tr>
<th align="RIGHT" nowrap="nowrap" valign="BASELINE">Subject:
</th>
<td>[BDXlist] BDX software: event builder, DSTs. First
ideas and thoughts</td>
</tr>
<tr>
<th align="RIGHT" nowrap="nowrap" valign="BASELINE">Date:
</th>
<td>Wed, 13 Jul 2016 17:33:15 +0200</td>
</tr>
<tr>
<th align="RIGHT" nowrap="nowrap" valign="BASELINE">From:
</th>
<td>Andrea Celentano <a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:Andrea.Celentano@ge.infn.it"><Andrea.Celentano@ge.infn.it></a></td>
</tr>
<tr>
<th align="RIGHT" nowrap="nowrap" valign="BASELINE">To: </th>
<td><a moz-do-not-send="true"
class="moz-txt-link-abbreviated"
href="mailto:bdxlist@jlab.org">bdxlist@jlab.org</a></td>
</tr>
</tbody>
</table>
<br>
<br>
<pre>Dear all,
this (quite long) e-mail is a first tentative to discuss how to proceed
with the BDX software, and specifically how to define an "event", and
how to produce from "events" an "output" (DST) that can be analyzed by
people not directly involved in the JANA framework development.
The following points are well defined:
1) Whatever the DAQ system will be, we can always talk about an "event"
as a single "entity" that is written on file. Events contain multiple
information from different detector elements, and are independent.
Events should be self-consistent, and do not depend on other events.
1A) In a "fixed DAQ window" scheme, this mean that all the acquisition
windows should be long enough to contain data coming at different times,
but physically related. Example: muon+ with e+ emission, e+ impinging on
a crystal. We can expect to have a hit in a VETO (due to entering mu+) ~
2 us before the hit in the crystal. These two hits must be in the same
"event", i.e. in the same piece of information written to the file.
1B) In a sophisticated trigger-less system, as the one we foresee, the
concept of event is much more dynamic and versatile. There will be one
(or multiple) online trigger algorithms looking at the full data-stream
from the full detector, identifying proper combination of hits - where a
"combination" is defined by the physics - and write that combination as
an "event".
1C) In the MC, this kind of implementation is very natural: each "event"
is associated with a single primary particle. All the hits resulting in
any evolution of this particle - and of its daughters - are an event.
2) Reconstruction is performed using the JANA framework. By
"reconstruction" I mean the algorithms that, from the RAW evio files -
or whatever format it will be - permits to obtain calibrated information
for each hit in each sub-detector. Also, in the reconstruction
elaborated quantities are produced - basically clusters in the
calorimeter, both "single-module" clusters and "multi-module" clusters.
2A) In the reconstruction, each event - as discussed in the previous
point - is absolutely independent from the others. It is not foreseen at
all to use information from event A when reconstructing event B.
3) DSTs will be ROOT files, with a ROOT tree, where a C++ object - the
"event" is saved. In other words, "events" are C++ objects. This gives
maximum flexibility in defining what to write in an event.
Here is what I think we should write in an event - and how do to this.
A) General data:
RunN
EventN
EventType (real, MC, ...)
Absolute event time (for real data)
... others ...
This should go in an "event header", that is by itself a C++ object. The
event has a pointer to this object.
B) Trigger data
This will really depend on the type of DAQ we will use, but information
about which trigger selected the event should be there.
C) MC-truth data
For MC-only: which kind of primary event was simulated, resulting in
this event? Again, a C++ object, the event has a pointer to this (for MC
only)
D) Low-level calibrated data
For each sub-detector, all the hits, calibrated in energy and time. Note
that a given sub-detector element can have more than one hit. Example:
mu+ entering ext veto, int veto, and stopping in crystal. Mu+ then decay
to e+, e+ exits from the crystal and hit again int. veto, the same
element as before. All this information should be in the same "event".
Each hit is a C++ object. Hits are contained in proper C++ collections -
can be simple vectors, or more powerful and performant collections, such
as ROOT TClonesArrays.
E) Elaborated data
I think that here we want to start with calorimeter clusters.
* Should we do clustering in the full detector, or start with clusters
in each calorimeter module, and then combine these?
* What are the other kind of data we need here?
--------------------------------------------------------
Finally, more "technical" question, that I think it is worth to discuss.
The structure of the event can be something very specific, i.e. a class
with specific data in it, already predefined, and not changeable, like:
class TEvent{
...
...
vector<CalorimeterHit> theCalorimeterHits;
CalorimeterCluster theCluster;
...
...
}
or something more versatile, like:
class TEvent{
...
...
vector<TClonesArray*> theRawHitCollections;
vector<TClonesArray*> theReconstructedObjectsCollections;
vector<TObject*> theObjectsInThisEvent;
...
...
}
The advantage of the first approach is that it is easier to use in an
analysis (you need the calorimeter hits: you have them in a vector), and
is self-explicative at the code-level. However, the structure of an
event is fixed, and can't change... (example, if later you want to to
add another vector in the event, you need to modify the Jana code that
is producing the event, and the event class itself, being potentially
not backward-compatible).
The advantage of the second approach (that I prefer :) ) is that the
structure of an event is really versatile. Each specific event (thanks
to ROOT reflectivity) is self-explicative, and one can do something like
(please note this is just pseudo-code):
TEvent event; //an event, get it from the DST
event->ListRawHitCollections() //which are the raw hits in this event?
...
if (event->hasCollection("CalorimeterHit")){
... loop on the CalorimeterHits ...
}
For a different project, I already developed an analysis code like this,
that I'll be happy to show at one of the next BDX meetings. The
disadvantage of this approach is a longer learning curve to use DSTs.
Bests,
Andrea
_______________________________________________
BDXlist mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:BDXlist@jlab.org">BDXlist@jlab.org</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="https://mailman.jlab.org/mailman/listinfo/bdxlist">https://mailman.jlab.org/mailman/listinfo/bdxlist</a>
</pre>
</div>
</blockquote>
<br>
</body>
</html>