]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TEvtGen/EvtGenBase/EvtIdSet.cxx
L1phase shift corrected
[u/mrichter/AliRoot.git] / TEvtGen / EvtGenBase / EvtIdSet.cxx
CommitLineData
da0e9ce3 1//--------------------------------------------------------------------------
2//
3// Environment:
4// This software is part of the EvtGen package developed jointly
5// for the BaBar and CLEO collaborations. If you use all or part
6// of it, please give an appropriate acknowledgement.
7//
8// Copyright Information: See EvtGen/COPYRIGHT
9// Copyright (C) 1998 Caltech, UCSB
10//
11// Module: EvtId.cc
12//
13// Description: Class for particle Id used in EvtGen.
14//
15// Modification history:
16//
17// DJL Jan 4, 2000 Module created
18//
19//------------------------------------------------------------------------
20//
21#include "EvtGenBase/EvtPatches.hh"
22
23#include "EvtGenBase/EvtIdSet.hh"
24#include "EvtGenBase/EvtPDL.hh"
25#include <iostream>
26#include <string>
27
28EvtIdSet::EvtIdSet(const EvtId name1) {
29 _numInList=1;
30 _list=new EvtId[_numInList];
31
32 _list[0]=name1;
33}
34
35EvtIdSet::EvtIdSet(const std::string name1){
36 _numInList=1;
37 _list=new EvtId[_numInList];
38
39 _list[0]=EvtPDL::getId(name1);
40}
41
42
43EvtIdSet::EvtIdSet(const EvtId name1,
44 const EvtId name2){
45 _numInList=2;
46 _list=new EvtId[_numInList];
47
48 _list[0]=name1;
49 _list[1]=name2;
50}
51
52
53EvtIdSet::EvtIdSet(const std::string name1,
54 const std::string name2){
55 _numInList=2;
56 _list=new EvtId[_numInList];
57
58 _list[0]=EvtPDL::getId(name1);
59 _list[1]=EvtPDL::getId(name2);
60}
61
62
63EvtIdSet::EvtIdSet(const EvtId name1,
64 const EvtId name2,
65 const EvtId name3){
66 _numInList=3;
67 _list=new EvtId[_numInList];
68
69 _list[0]=name1;
70 _list[1]=name2;
71 _list[2]=name3;
72
73}
74
75
76EvtIdSet::EvtIdSet(const std::string name1,
77 const std::string name2,
78 const std::string name3){
79 _numInList=3;
80 _list=new EvtId[_numInList];
81
82 _list[0]=EvtPDL::getId(name1);
83 _list[1]=EvtPDL::getId(name2);
84 _list[2]=EvtPDL::getId(name3);
85}
86
87
88EvtIdSet::EvtIdSet(const EvtId name1,
89 const EvtId name2,
90 const EvtId name3,
91 const EvtId name4){
92 _numInList=4;
93 _list=new EvtId[_numInList];
94
95 _list[0]=name1;
96 _list[1]=name2;
97 _list[2]=name3;
98 _list[3]=name4;
99}
100
101
102EvtIdSet::EvtIdSet(const std::string name1,
103 const std::string name2,
104 const std::string name3,
105 const std::string name4){
106 _numInList=4;
107 _list=new EvtId[_numInList];
108
109 _list[0]=EvtPDL::getId(name1);
110 _list[1]=EvtPDL::getId(name2);
111 _list[2]=EvtPDL::getId(name3);
112 _list[3]=EvtPDL::getId(name4);
113}
114
115
116EvtIdSet::EvtIdSet(const EvtId name1,
117 const EvtId name2,
118 const EvtId name3,
119 const EvtId name4,
120 const EvtId name5){
121 _numInList=5;
122 _list=new EvtId[_numInList];
123
124 _list[0]=name1;
125 _list[1]=name2;
126 _list[2]=name3;
127 _list[3]=name4;
128 _list[4]=name5;
129}
130
131
132EvtIdSet::EvtIdSet(const std::string name1,
133 const std::string name2,
134 const std::string name3,
135 const std::string name4,
136 const std::string name5){
137 _numInList=5;
138 _list=new EvtId[_numInList];
139
140 _list[0]=EvtPDL::getId(name1);
141 _list[1]=EvtPDL::getId(name2);
142 _list[2]=EvtPDL::getId(name3);
143 _list[3]=EvtPDL::getId(name4);
144 _list[4]=EvtPDL::getId(name5);
145}
146
147
148EvtIdSet::EvtIdSet(const EvtId name1,
149 const EvtId name2,
150 const EvtId name3,
151 const EvtId name4,
152 const EvtId name5,
153 const EvtId name6){
154 _numInList=6;
155 _list=new EvtId[_numInList];
156
157 _list[0]=name1;
158 _list[1]=name2;
159 _list[2]=name3;
160 _list[3]=name4;
161 _list[4]=name5;
162 _list[5]=name6;
163}
164
165
166EvtIdSet::EvtIdSet(const std::string name1,
167 const std::string name2,
168 const std::string name3,
169 const std::string name4,
170 const std::string name5,
171 const std::string name6){
172 _numInList=6;
173 _list=new EvtId[_numInList];
174
175 _list[0]=EvtPDL::getId(name1);
176 _list[1]=EvtPDL::getId(name2);
177 _list[2]=EvtPDL::getId(name3);
178 _list[3]=EvtPDL::getId(name4);
179 _list[4]=EvtPDL::getId(name5);
180 _list[5]=EvtPDL::getId(name6);
181}
182
183
184EvtIdSet::EvtIdSet(const EvtId name1,
185 const EvtId name2,
186 const EvtId name3,
187 const EvtId name4,
188 const EvtId name5,
189 const EvtId name6,
190 const EvtId name7){
191 _numInList=7;
192 _list=new EvtId[_numInList];
193
194 _list[0]=name1;
195 _list[1]=name2;
196 _list[2]=name3;
197 _list[3]=name4;
198 _list[4]=name5;
199 _list[5]=name6;
200 _list[6]=name7;
201}
202
203
204EvtIdSet::EvtIdSet(const std::string name1,
205 const std::string name2,
206 const std::string name3,
207 const std::string name4,
208 const std::string name5,
209 const std::string name6,
210 const std::string name7){
211 _numInList=7;
212 _list=new EvtId[_numInList];
213
214 _list[0]=EvtPDL::getId(name1);
215 _list[1]=EvtPDL::getId(name2);
216 _list[2]=EvtPDL::getId(name3);
217 _list[3]=EvtPDL::getId(name4);
218 _list[4]=EvtPDL::getId(name5);
219 _list[5]=EvtPDL::getId(name6);
220 _list[6]=EvtPDL::getId(name7);
221}
222
223
224EvtIdSet::EvtIdSet(const EvtId name1,
225 const EvtId name2,
226 const EvtId name3,
227 const EvtId name4,
228 const EvtId name5,
229 const EvtId name6,
230 const EvtId name7,
231 const EvtId name8){
232 _numInList=8;
233 _list=new EvtId[_numInList];
234
235 _list[0]=name1;
236 _list[1]=name2;
237 _list[2]=name3;
238 _list[3]=name4;
239 _list[4]=name5;
240 _list[5]=name6;
241 _list[6]=name7;
242 _list[7]=name8;
243}
244
245
246EvtIdSet::EvtIdSet(const std::string name1,
247 const std::string name2,
248 const std::string name3,
249 const std::string name4,
250 const std::string name5,
251 const std::string name6,
252 const std::string name7,
253 const std::string name8){
254 _numInList=8;
255 _list=new EvtId[_numInList];
256
257 _list[0]=EvtPDL::getId(name1);
258 _list[1]=EvtPDL::getId(name2);
259 _list[2]=EvtPDL::getId(name3);
260 _list[3]=EvtPDL::getId(name4);
261 _list[4]=EvtPDL::getId(name5);
262 _list[5]=EvtPDL::getId(name6);
263 _list[6]=EvtPDL::getId(name7);
264 _list[7]=EvtPDL::getId(name8);
265}
266
267
268EvtIdSet::EvtIdSet(const EvtId name1,
269 const EvtId name2,
270 const EvtId name3,
271 const EvtId name4,
272 const EvtId name5,
273 const EvtId name6,
274 const EvtId name7,
275 const EvtId name8,
276 const EvtId name9){
277 _numInList=9;
278 _list=new EvtId[_numInList];
279
280 _list[0]=name1;
281 _list[1]=name2;
282 _list[2]=name3;
283 _list[3]=name4;
284 _list[4]=name5;
285 _list[5]=name6;
286 _list[6]=name7;
287 _list[7]=name8;
288 _list[8]=name9;
289}
290
291
292EvtIdSet::EvtIdSet(const std::string name1,
293 const std::string name2,
294 const std::string name3,
295 const std::string name4,
296 const std::string name5,
297 const std::string name6,
298 const std::string name7,
299 const std::string name8,
300 const std::string name9){
301 _numInList=9;
302 _list=new EvtId[_numInList];
303
304 _list[0]=EvtPDL::getId(name1);
305 _list[1]=EvtPDL::getId(name2);
306 _list[2]=EvtPDL::getId(name3);
307 _list[3]=EvtPDL::getId(name4);
308 _list[4]=EvtPDL::getId(name5);
309 _list[5]=EvtPDL::getId(name6);
310 _list[6]=EvtPDL::getId(name7);
311 _list[7]=EvtPDL::getId(name8);
312 _list[8]=EvtPDL::getId(name9);
313}
314
315
316EvtIdSet::EvtIdSet(const EvtId name1,
317 const EvtId name2,
318 const EvtId name3,
319 const EvtId name4,
320 const EvtId name5,
321 const EvtId name6,
322 const EvtId name7,
323 const EvtId name8,
324 const EvtId name9,
325 const EvtId name10){
326 _numInList=10;
327 _list=new EvtId[_numInList];
328
329 _list[0]=name1;
330 _list[1]=name2;
331 _list[2]=name3;
332 _list[3]=name4;
333 _list[4]=name5;
334 _list[5]=name6;
335 _list[6]=name7;
336 _list[7]=name8;
337 _list[8]=name9;
338 _list[9]=name10;
339}
340
341
342EvtIdSet::EvtIdSet(const std::string name1,
343 const std::string name2,
344 const std::string name3,
345 const std::string name4,
346 const std::string name5,
347 const std::string name6,
348 const std::string name7,
349 const std::string name8,
350 const std::string name9,
351 const std::string name10){
352 _numInList=10;
353 _list=new EvtId[_numInList];
354
355 _list[0]=EvtPDL::getId(name1);
356 _list[1]=EvtPDL::getId(name2);
357 _list[2]=EvtPDL::getId(name3);
358 _list[3]=EvtPDL::getId(name4);
359 _list[4]=EvtPDL::getId(name5);
360 _list[5]=EvtPDL::getId(name6);
361 _list[6]=EvtPDL::getId(name7);
362 _list[7]=EvtPDL::getId(name8);
363 _list[8]=EvtPDL::getId(name9);
364 _list[9]=EvtPDL::getId(name10);
365}
366
367
368EvtIdSet::EvtIdSet(const EvtId name1,
369 const EvtId name2,
370 const EvtId name3,
371 const EvtId name4,
372 const EvtId name5,
373 const EvtId name6,
374 const EvtId name7,
375 const EvtId name8,
376 const EvtId name9,
377 const EvtId name10,
378 const EvtId name11){
379 _numInList=11;
380 _list=new EvtId[_numInList];
381
382 _list[0]=name1;
383 _list[1]=name2;
384 _list[2]=name3;
385 _list[3]=name4;
386 _list[4]=name5;
387 _list[5]=name6;
388 _list[6]=name7;
389 _list[7]=name8;
390 _list[8]=name9;
391 _list[9]=name10;
392 _list[10]=name11;
393}
394
395
396EvtIdSet::EvtIdSet(const std::string name1,
397 const std::string name2,
398 const std::string name3,
399 const std::string name4,
400 const std::string name5,
401 const std::string name6,
402 const std::string name7,
403 const std::string name8,
404 const std::string name9,
405 const std::string name10,
406 const std::string name11){
407 _numInList=11;
408 _list=new EvtId[_numInList];
409
410 _list[0]=EvtPDL::getId(name1);
411 _list[1]=EvtPDL::getId(name2);
412 _list[2]=EvtPDL::getId(name3);
413 _list[3]=EvtPDL::getId(name4);
414 _list[4]=EvtPDL::getId(name5);
415 _list[5]=EvtPDL::getId(name6);
416 _list[6]=EvtPDL::getId(name7);
417 _list[7]=EvtPDL::getId(name8);
418 _list[8]=EvtPDL::getId(name9);
419 _list[9]=EvtPDL::getId(name10);
420 _list[10]=EvtPDL::getId(name11);
421}
422
423
424EvtIdSet::EvtIdSet(const EvtId name1,
425 const EvtId name2,
426 const EvtId name3,
427 const EvtId name4,
428 const EvtId name5,
429 const EvtId name6,
430 const EvtId name7,
431 const EvtId name8,
432 const EvtId name9,
433 const EvtId name10,
434 const EvtId name11,
435 const EvtId name12){
436 _numInList=12;
437 _list=new EvtId[_numInList];
438
439 _list[0]=name1;
440 _list[1]=name2;
441 _list[2]=name3;
442 _list[3]=name4;
443 _list[4]=name5;
444 _list[5]=name6;
445 _list[6]=name7;
446 _list[7]=name8;
447 _list[8]=name9;
448 _list[9]=name10;
449 _list[10]=name11;
450 _list[11]=name12;
451}
452
453
454EvtIdSet::EvtIdSet(const std::string name1,
455 const std::string name2,
456 const std::string name3,
457 const std::string name4,
458 const std::string name5,
459 const std::string name6,
460 const std::string name7,
461 const std::string name8,
462 const std::string name9,
463 const std::string name10,
464 const std::string name11,
465 const std::string name12){
466 _numInList=12;
467 _list=new EvtId[_numInList];
468
469 _list[0]=EvtPDL::getId(name1);
470 _list[1]=EvtPDL::getId(name2);
471 _list[2]=EvtPDL::getId(name3);
472 _list[3]=EvtPDL::getId(name4);
473 _list[4]=EvtPDL::getId(name5);
474 _list[5]=EvtPDL::getId(name6);
475 _list[6]=EvtPDL::getId(name7);
476 _list[7]=EvtPDL::getId(name8);
477 _list[8]=EvtPDL::getId(name9);
478 _list[9]=EvtPDL::getId(name10);
479 _list[10]=EvtPDL::getId(name11);
480 _list[11]=EvtPDL::getId(name12);
481}
482
483
484EvtIdSet::EvtIdSet(const EvtIdSet& set1){
485
486 _numInList=set1.sizeOfSet();
487 _list=new EvtId[_numInList];
488 int i;
489 for (i=0;i<_numInList;i++){
490 _list[i]=set1.getElem(i);
491 }
492
493}
494EvtIdSet::EvtIdSet(const EvtIdSet& set1, const EvtIdSet& set2){
495
496 _numInList=set1.sizeOfSet();
497 _list=new EvtId[_numInList];
498 int i;
499 for (i=0;i<_numInList;i++){
500 _list[i]=set1.getElem(i);
501 }
502 //then just append the second list.
503 this->append(set2);
504
505}
506
507int EvtIdSet::contains(const EvtId id){
508
509 int i;
510 for (i=0;i<_numInList;i++){
511 if ( _list[i] == id ) return 1;
512 }
513
514 return 0;
515}
516
517int EvtIdSet::contains(const std::string nm){
518
519 int i;
520 for (i=0;i<_numInList;i++){
521 if ( _list[i] == EvtPDL::getId(nm) ) return 1;
522 }
523
524 return 0;
525}
526
527
528void EvtIdSet::append(const EvtIdSet set1){
529
530 int combLen=_numInList+set1.sizeOfSet();
531 int uniqueLen=0;
532 EvtId *combSet;
533
534 combSet=new EvtId[combLen];
535
536 int i;
537 for (i=0;i<combLen;i++){
538 if ( i>=_numInList ) {
539
540 //check that there are no overlaps between lists
541 int j;
542 int isUnique=1;
543 for (j=0;j<_numInList;j++){
544 if ( _list[j]==set1.getElem(i-_numInList) ) {
545 isUnique=0;
546 }
547 }
548 if ( isUnique==1 ) {
549 combSet[uniqueLen]=set1.getElem(i-_numInList);
550 uniqueLen+=1;
551 }
552 }
553 else{
554 combSet[uniqueLen]=_list[i];
555 uniqueLen+=1;
556 }
557
558 delete _list;
559 _list=new EvtId[uniqueLen];
560
561 _numInList=uniqueLen;
562 for (i=0;i<_numInList;i++){
563 _list[i]=combSet[i];
564 }
565
566 delete combSet;
567
568 }
569}
570
571int EvtIdSet::sizeOfSet() const { return _numInList;}
572
573EvtId EvtIdSet::getElem(const int i) const { return _list[i];}
574
575
576
577
578
579
580
581