16 int kqa,kq1,kq2,kq3,kqj,irt,kqx,kqn;
20 int ichg[109]={-1,2,-1,2,-1,2,-1,2,0,0,-3,0,-3,0,-3,0,
21 -3,0,0,0,0,0,0,3,0,0,0,0,0,0,3,0,3,6,0,0,3,6,0,0,-1,2,-1,2,-1,2,0,0,0,0,
22 -3,0,-3,0,-3,0,0,0,0,0,-1,2,-1,2,-1,2,0,0,0,0,
23 -3,0,-3,0,-3,0,3,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
29 kqn=kqa/1000000000%10;
39 if(kqa==0 || kqa >= 10000000) {
41 if(kqn==1) {hepchg=0;}
44 else if(kqa<=100) {hepchg = ichg[kqa-1];}
46 else if(kqa==100 || kqa==101) {hepchg = -3;}
48 else if(kqa==102 || kqa==104) {hepchg = -6;}
50 else if(kqj == 0) {hepchg = 0;}
52 else if(kqx>0 && irt<100)
55 if(kqa==1000017 || kqa==1000018) {hepchg = 0;}
56 if(kqa==1000034 || kqa==1000052) {hepchg = 0;}
57 if(kqa==1000053 || kqa==1000054) {hepchg = 0;}
58 if(kqa==5100061 || kqa==5100062) {hepchg = 6;}
64 hepchg = ichg[kq2-1]-ichg[kq1-1];
66 if((kq2==3) || (kq2==5)) {hepchg = ichg[kq1-1]-ichg[kq2-1];}
70 hepchg = ichg[kq3-1] + ichg[kq2-1];
75 hepchg = ichg[kq3-1]+ichg[kq2-1]+ichg[kq1-1];
79 if(Id<0 && hepchg!=0) {hepchg = -1*hepchg;}
89 o << t.
type() <<
", ";
std::ostream & operator<<(std::ostream &out, const ALILine &li)
float charge() const
charge
int chargeValue(const int &) const
int type() const
particle type (HEP PDT convension)
const math::XYZTLorentzVectorD & momentum() const
particle info...