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 {
83 std::cout <<
"Lambda0's id:" << p->
getId() <<
" p4 =" << p->
getP4Lab() << std::endl;
84
87 std::cout << "FromParent_Lambda: pi- p4[numstable] = " << p4[numstable] << std::endl;
88 numstable++;
89
92 std::cout << "FromParent_Lambda: p+ p4[numstable] = " << p4[numstable] << std::endl;
93 numstable++;
94 }
95 if ( p->
getId() == antiLambdaId )
96 {
97 std::cout <<
"anti-Lambda0's parent id:" << p->
getParent()->
getId()
99 std::cout <<
"anti-Lambda0's id:" << p->
getId() <<
" p4 =" << p->
getP4Lab() << std::endl;
100
103 std::cout << "FromParent_Lambda: pi+ p4[numstable] = " << p4[numstable] << std::endl;
104 numstable++;
105
108 std::cout << "FromParent_Lambda: anti-p- p4[numstable] = " << p4[numstable] << std::endl;
109 numstable++;
110 }
112 int ndaugFound = 0;
113 for ( int i = 0; i < numstable; i++ )
114 {
116 ndaugFound++;
117 }
118 if ( ndaugFound == 0 )
119 {
120 report(
ERROR,
"EvtGen" ) <<
"Phokhara has failed to do a decay ";
122 << endl;
123 assert( 0 );
124 }
125 if ( p->
getId() == LambdaId )
126 {
127
128
129
130
131
132 std::cout <<
"EvtFromParent_Lambda Lambda0's parent id:" << p->
getParent()->
getId()
134 std::cout <<
"EvtFromParent_Lambda SUMMARY: part p4" << p->
getP4Lab() << std::endl;
135 std::cout <<
"EvtFromParent_Lambda SUMMARY: Daug0 p4" << p->
getDaug( 0 )->
getP4Lab()
136 << std::endl;
137 std::cout <<
"EvtFromParent_Lambda SUMMARY: Daug1 p4" << p->
getDaug( 1 )->
getP4Lab()
138 << std::endl;
139 }
140 if ( p->
getId() == antiLambdaId )
141 {
142
143
144
145
146
147 std::cout <<
"EvtFromParent_Lambda antiLambda0's parent id:" << p->
getParent()->
getId()
149 std::cout <<
"EvtFromParent_Lambda SUMMARY: part p4" << p->
getP4Lab() << std::endl;
150 std::cout <<
"EvtFromParent_Lambda SUMMARY: Daug0 p4" << p->
getDaug( 0 )->
getP4Lab()
151 << std::endl;
152 std::cout <<
"EvtFromParent_Lambda SUMMARY: Daug1 p4" << p->
getDaug( 1 )->
getP4Lab()
153 << std::endl;
154 }
155 }
156 return;
157}
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]