57 {
58
59
61
63
64 if ( p->
getId() != LambdaId && p->
getId() != antiLambdaId )
65 {
66 std::cout << "Parent particle is required to be Lambda0 or antiLambda0" << std::endl;
67 abort();
68 }
69 int more = 0;
70 int numstable = 0;
71 EvtId evtnumstable[100];
72 EvtVector4R p4[20];
73
74
75
76
78 {
79 if ( p->
getId() == LambdaId )
80 {
81
84 std::cout << "FromParent_antiLambda: pi- p4[numstable] = " << p4[numstable] << std::endl;
85 numstable++;
86
89 std::cout << "FromParent_antiLambda: p+ p4[numstable] = " << p4[numstable] << std::endl;
90 numstable++;
91 }
92 if ( p->
getId() == antiLambdaId )
93 {
94 std::cout <<
"anti-Lambda0's parent id:" << p->
getParent()->
getId()
96 std::cout <<
"anti-Lambda0 id:" << p->
getId() <<
" p4 =" << p->
getP4Lab() << std::endl;
97
100 std::cout << "FromParent_antiLambda: pi+ p4[numstable] = " << p4[numstable] << std::endl;
101 numstable++;
102
105 std::cout << "FromParent_antiLambda: anti-p- p4[numstable] = " << p4[numstable]
106 << std::endl;
107 numstable++;
108 }
110 int ndaugFound = 0;
111 for ( int i = 0; i < numstable; i++ )
112 {
114 ndaugFound++;
115 }
116 if ( ndaugFound == 0 )
117 {
118 report(
ERROR,
"EvtGen" ) <<
"Phokhara has failed to do a decay ";
120 << endl;
121 assert( 0 );
122 }
123 if ( p->
getId() == LambdaId )
124 {
125
126
127
128
129
130 std::cout <<
"EvtFromParent_antiLambda Lambda0's parent id:" << p->
getParent()->
getId()
132 std::cout <<
"EvtFromParent_antiLambda SUMMARY: part p4" << p->
getP4Lab() << std::endl;
133 std::cout <<
"EvtFromParent_antiLambda SUMMARY: Daug0 p4" << p->
getDaug( 0 )->
getP4Lab()
134 << std::endl;
135 std::cout <<
"EvtFromParent_antiLambda SUMMARY: Daug1 p4" << p->
getDaug( 1 )->
getP4Lab()
136 << std::endl;
137 }
138 if ( p->
getId() == antiLambdaId )
139 {
140
141
142
143
144
145 std::cout << "EvtFromParent_antiLambda antiLambda0's parent id:"
147 << std::endl;
148 std::cout <<
"EvtFromParent_antiLambda SUMMARY: part p4" << p->
getP4Lab() << std::endl;
149 std::cout <<
"EvtFromParent_antiLambda SUMMARY: Daug0 p4" << p->
getDaug( 0 )->
getP4Lab()
150 << std::endl;
151 std::cout <<
"EvtFromParent_antiLambda SUMMARY: Daug1 p4" << p->
getDaug( 1 )->
getP4Lab()
152 << std::endl;
153 }
154 }
155 return;
156}
ostream & report(Severity severity, const char *facility)
static EvtId evtIdFromStdHep(int stdhep)
static std::string name(EvtId i)
static EvtId getId(const std::string &name)
void makeDaughters(int ndaug, EvtId *id)
virtual void init(EvtId part_n, const EvtVector4R &p4)=0
EvtParticle * getParent()
static int _NextLevelDauNum
EvtParticle * getDaug(int i)
static EvtVector4R _NextLevelP4[20]