00001 /*********** using the following decay cards **************** 00002 Decay J/psi 00003 1 gamma eta_c PHSP; 00004 Enddecay 00005 00006 Decay eta_c 00007 1 phi phi PHSP; 00008 Enddecay 00009 00010 Decay phi 00011 1 K+ K- PHSP; 00012 Enddecay 00013 00014 End 00015 */ 00016 //---------------------------------- 00017 EvtVector4R pjsi,petac,pphi1,pphi2,pkp1,pkp2; 00018 pjsi = part->getP4(); 00019 petac= part->getDaug(1)->getP4(); //etac momentum 00020 pphi1= part->getDaug(1)->getDaug(0)->getP4(); //phi1 momentum 00021 pphi2= part->getDaug(1)->getDaug(1)->getP4(); //phi2 momentum 00022 00023 pkp1 = part->getDaug(1)->getDaug(0)->getDaug(0)->getP4(); //K+ from phi1 00024 pkp2 = part->getDaug(1)->getDaug(1)->getDaug(0)->getP4(); //K+ from phi2 00025 00026 EvtHelSys angles0(pjsi,petac); 00027 double theta0 = angles0.getHelAng(1); 00028 double phi0 = angles0.getHelAng(2); 00029 00030 EvtHelSys angles2(pphi1,pkp1); 00031 double theta2 = angles2.getHelAng(1); 00032 double phi2 = angles2.getHelAng(2); 00033 00034 EvtHelSys angles3(pphi2,pkp2); 00035 double theta3 = angles3.getHelAng(1); 00036 double phi3 = angles3.getHelAng(2); 00037 00038 double amps=(3+cos(2*theta0)) * sin(theta2)*sin(theta2) * sin(theta3)*sin(theta3) * sin(phi2+phi3)*sin(phi2+phi3); 00039