[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Small Sample Rare Eta MC
Hi Guys,
I ran into Pawel at the grocery store tonight and he mentioned that
it might be useful to have this file in HBOOK/Ntuple format. So I wrote
a quick utility to convert it. The hbook version can be found here:
http://www.jlab.org/~davidl/PrimEx/primex_event.hbook
(Matt, feel free to copy this over to your web page so it is kept in the
same spot as the ROOT version.)
The source code for the conversion utility is in the Hall-D repository at:
https://halldsvn.jlab.org/repos/trunk/src/programs/Utilities/root2hbook
I hope this is useful.
Regards,
-David
Matthew Shepherd wrote:
>
> Hi Liping and Ashot,
>
> First of all, we may want to eventually keep traffic off of this
> halld-physics list. I imagine that once we start exchanging files
> there are going to be lots of technical questions. I'll send this
> first note to the whole list but maybe subsequent discussion should
> happen directly.
>
> I put a sample root file here:
>
> http://dustbunny.physics.indiana.edu/~mashephe/primex/
>
> Which comes from a subset of Lipings rare eta events.
>
> The file contains one tree called "event" (The dump of the fields is
> pasted below.) and a simple histogram of the total invariant mass of
> all photons in the calorimeter.
>
> The event tree is mainly two sets of arrays. One is indexed on n and
> contains the reconstructed four vectors of the n photons found in the
> detector. The other indexed on nGen contains the generated four
> vectors and type of the nGen generated pi^0, etas, or photons. These
> generated types correspond just to those you supplied in the ascii
> file and does not consider decays for example of pi^0 -> gamma gamma.
> So each of Liping's events will have 3 generated particles: 2 photons
> and a pi^0. Each of Ashot's will have 2 particles, just the two
> photons. (I do not record generator info for other particle types
> since it isn't of use for us at the moment.)
>
> Then genNeuE contains the total energy of all neutral particles and
> electrons and positrons. When generating the background I plan to
> preselect events by placing a loose cut on this variable. This way I
> do not spend huge CPU time to generate events that we will certainly
> toss out. This is important to be able to get enough statistics to
> get some feel for the background. As long as this cut is loose enough
> then I don't have to worry about this cut artificially reducing the
> background. A similar variable, neuETot, is the sum of all the
> reconstructed neutral energy (a good first variable to cut on in an
> analysis).
>
> I have not included raw hit information. This would increase file
> size dramatically. It could easily be done, but after our phone
> conversation today I'm not sure if you have plans to use it right away.
>
> -Matt
>
>
> ******************************************************************************
>
> *Tree :event : Event
> Reconstruction *
> *Entries : 10000 : Total = 1882594 bytes File Size =
> 1363534 *
> * : : Tree compression factor =
> 1.37 *
> ******************************************************************************
>
> *Br 0 :n :
> n/I *
> *Entries : 10000 : Total Size= 40677 bytes File Size
> = 4586 *
> *Baskets : 1 : Basket Size= 32000 bytes Compression=
> 6.98 *
> *............................................................................*
>
> *Br 1 :E :
> E[n]/F *
> *Entries : 10000 : Total Size= 228491 bytes File Size =
> 183329 *
> *Baskets : 8 : Basket Size= 32000 bytes Compression=
> 1.19 *
> *............................................................................*
>
> *Br 2 :px :
> px[n]/F *
> *Entries : 10000 : Total Size= 228505 bytes File Size =
> 187776 *
> *Baskets : 8 : Basket Size= 32000 bytes Compression=
> 1.16 *
> *............................................................................*
>
> *Br 3 :py :
> py[n]/F *
> *Entries : 10000 : Total Size= 228505 bytes File Size =
> 187842 *
> *Baskets : 8 : Basket Size= 32000 bytes Compression=
> 1.16 *
> *............................................................................*
>
> *Br 4 :pz :
> pz[n]/F *
> *Entries : 10000 : Total Size= 228505 bytes File Size =
> 184061 *
> *Baskets : 8 : Basket Size= 32000 bytes Compression=
> 1.18 *
> *............................................................................*
>
> *Br 5 :neuETot :
> neuETot/F *
> *Entries : 10000 : Total Size= 40719 bytes File Size =
> 26283 *
> *Baskets : 1 : Basket Size= 32000 bytes Compression=
> 1.22 *
> *............................................................................*
>
> *Br 6 :nGen :
> nGen/I *
> *Entries : 10000 : Total Size= 40698 bytes File Size
> = 263 *
> *Baskets : 1 : Basket Size= 32000 bytes Compression=
> 121.66 *
> *............................................................................*
>
> *Br 7 :genE :
> genE[nGen]/F *
> *Entries : 10000 : Total Size= 161217 bytes File Size =
> 118567 *
> *Baskets : 6 : Basket Size= 32000 bytes Compression=
> 1.30 *
> *............................................................................*
>
> *Br 8 :genPx :
> genPx[nGen]/F *
> *Entries : 10000 : Total Size= 161229 bytes File Size =
> 122566 *
> *Baskets : 6 : Basket Size= 32000 bytes Compression=
> 1.25 *
> *............................................................................*
>
> *Br 9 :genPy :
> genPy[nGen]/F *
> *Entries : 10000 : Total Size= 161229 bytes File Size =
> 122616 *
> *Baskets : 6 : Basket Size= 32000 bytes Compression=
> 1.25 *
> *............................................................................*
>
> *Br 10 :genPz :
> genPz[nGen]/F *
> *Entries : 10000 : Total Size= 161229 bytes File Size =
> 118466 *
> *Baskets : 6 : Basket Size= 32000 bytes Compression=
> 1.30 *
> *............................................................................*
>
> *Br 11 :type :
> type[nGen]/I *
> *Entries : 10000 : Total Size= 161210 bytes File Size =
> 14124 *
> *Baskets : 6 : Basket Size= 32000 bytes Compression=
> 10.88 *
> *............................................................................*
>
> *Br 12 :genNeuE :
> genNeuE/F *
> *Entries : 10000 : Total Size= 40719 bytes File Size =
> 20990 *
> *Baskets : 1 : Basket Size= 32000 bytes Compression=
> 1.52 *
> *............................................................................*
>
>
--
------------------------------------------------------------------------
David Lawrence Ph.D.
Staff Scientist Office: (757)269-5567 [[[ [ [ [
Jefferson Lab Pager: (757)584-5567 [ [ [ [ [ [
http://www.jlab.org/~davidl davidl@jlab.org [[[ [[ [[ [[[
------------------------------------------------------------------------