39 int main(
int narg,
char *argv[])
41 const char *fname =
"primex_event.root";
42 if(narg>1)fname = argv[1];
45 TFile *
f =
new TFile(fname);
46 if(!f || !f->IsOpen()){
47 cerr<<
"Unable to open file \""<<fname<<
"\"!"<<endl;
50 cout<<
"Opened input file \""<<fname<<
"\""<<endl;
53 TTree *tree = (TTree*)f->Get(
"event");
55 cerr<<
"Unable to find TTree \"event\" in ROOT file!"<<endl;
61 tree->SetBranchAddress(
"n", &event.
n);
62 tree->SetBranchAddress(
"E", &event.
E);
63 tree->SetBranchAddress(
"px", &event.
px);
64 tree->SetBranchAddress(
"py", &event.
py);
65 tree->SetBranchAddress(
"pz", &event.
pz);
66 tree->SetBranchAddress(
"neuETot", &event.
neuETot);
67 tree->SetBranchAddress(
"nGen", &event.
nGen);
68 tree->SetBranchAddress(
"genE", &event.
genE);
69 tree->SetBranchAddress(
"genPx", &event.
genPx);
70 tree->SetBranchAddress(
"genPy", &event.
genPy);
71 tree->SetBranchAddress(
"genPz", &event.
genPz);
72 tree->SetBranchAddress(
"type", &event.
type);
73 tree->SetBranchAddress(
"genNeuE", &event.
genNeuE);
86 sprintf(ntp_str,
"n[0,%d]:I,E(n):R,px(n),py(n),pz(n),neuETot,nGen[0,%d]:I,genE(nGen):R,genPx(nGen),genPy(nGen),genPz(nGen),type(nGen):I,genNeuE:R",
MAX_PARTS,
MAX_PARTS);
88 hbname(10,
"EVNT", &event.
n, ntp_str);
91 int Nevents = tree->GetEntries();
100 if(i%100 == 0)cout<<
" "<<i<<
" events processed \r";
104 cout<<Nevents<<
" events processed total."<<endl;
119 ******************************************************************************
120 *Tree :
event :
Event Reconstruction *
121 *Entries : 10000 : Total = 1882594 bytes File Size = 1363534 *
122 * : : Tree compression factor = 1.37 *
123 ******************************************************************************
125 *Entries : 10000 : Total Size= 40677 bytes File Size = 4586 *
126 *Baskets : 1 : Basket Size= 32000 bytes Compression= 6.98 *
127 *............................................................................*
129 *Entries : 10000 : Total Size= 228491 bytes File Size = 183329 *
130 *Baskets : 8 : Basket Size= 32000 bytes Compression= 1.19 *
131 *............................................................................*
132 *Br 2 :px : px[n]/
F *
133 *Entries : 10000 : Total Size= 228505 bytes File Size = 187776 *
134 *Baskets : 8 : Basket Size= 32000 bytes Compression= 1.16 *
135 *............................................................................*
136 *Br 3 :
py :
py[n]/
F *
137 *Entries : 10000 : Total Size= 228505 bytes File Size = 187842 *
138 *Baskets : 8 : Basket Size= 32000 bytes Compression= 1.16 *
139 *............................................................................*
140 *Br 4 :pz : pz[n]/
F *
141 *Entries : 10000 : Total Size= 228505 bytes File Size = 184061 *
142 *Baskets : 8 : Basket Size= 32000 bytes Compression= 1.18 *
143 *............................................................................*
144 *Br 5 :neuETot : neuETot/
F *
145 *Entries : 10000 : Total Size= 40719 bytes File Size = 26283 *
146 *Baskets : 1 : Basket Size= 32000 bytes Compression= 1.22 *
147 *............................................................................*
148 *Br 6 :nGen : nGen/
I *
149 *Entries : 10000 : Total Size= 40698 bytes File Size = 263 *
150 *Baskets : 1 : Basket Size= 32000 bytes Compression= 121.66 *
151 *............................................................................*
152 *Br 7 :genE : genE[nGen]/
F *
153 *Entries : 10000 : Total Size= 161217 bytes File Size = 118567 *
154 *Baskets : 6 : Basket Size= 32000 bytes Compression= 1.30 *
155 *............................................................................*
156 *Br 8 :genPx : genPx[nGen]/
F *
157 *Entries : 10000 : Total Size= 161229 bytes File Size = 122566 *
158 *Baskets : 6 : Basket Size= 32000 bytes Compression= 1.25 *
159 *............................................................................*
160 *Br 9 :genPy : genPy[nGen]/
F *
161 *Entries : 10000 : Total Size= 161229 bytes File Size = 122616 *
162 *Baskets : 6 : Basket Size= 32000 bytes Compression= 1.25 *
163 *............................................................................*
164 *Br 10 :genPz : genPz[nGen]/
F *
165 *Entries : 10000 : Total Size= 161229 bytes File Size = 118466 *
166 *Baskets : 6 : Basket Size= 32000 bytes Compression= 1.30 *
167 *............................................................................*
168 *Br 11 :type : type[nGen]/
I *
169 *Entries : 10000 : Total Size= 161210 bytes File Size = 14124 *
170 *Baskets : 6 : Basket Size= 32000 bytes Compression= 10.88 *
171 *............................................................................*
172 *Br 12 :genNeuE : genNeuE/
F *
173 *Entries : 10000 : Total Size= 40719 bytes File Size = 20990 *
174 *Baskets : 1 : Basket Size= 32000 bytes Compression= 1.52 *
175 *............................................................................*
void hbname(int id, const char *CHBLOK, void *VAR, const char *CHFORM)
void hropen(int lun, char *name, char *filename, char *status, int stor, int istat)
sprintf(text,"Post KinFit Cut")
void hbnt(int id, char *CHTITL, char *CHOPT)
void hrout(int num, int icycle, char *opt)
void hrend(char *filename)
int main(int argc, char *argv[])