BOSS_7.0.2: /home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Generator/BesEvtGen/BesEvtGen-00-03-58/src/EvtGen/EvtGenModels/EvtChi2BB2.cc Source File

00001 //--------------------------------------------------------------------------
00002 //
00003 // Environment:
00004 //      This software is part of models developed at BES collaboration
00005 //      based on the EvtGen framework.  If you use all or part
00006 //      of it, please give an appropriate acknowledgement.
00007 //
00008 // Copyright Information: See EvtGen/BesCopyright
00009 //      Copyright (A) 2006      Ping Rong-Gang @IHEP
00010 //
00011 // Module: EvtChi2BB2.cc
00012 //
00013 // Description: Routine to decay Chi2 to B10 B10bar.
00014 //              
00015 // Modification history:
00016 //
00017 // Pang C.Y. and Ping R.-G.       April, 2007       Module created
00018 //
00019 //------------------------------------------------------------------------
00020 // 
00021 #include "EvtGenBase/EvtPatches.hh"
00022 #include <stdlib.h>
00023 #include "EvtGenBase/EvtParticle.hh"
00024 #include "EvtGenBase/EvtGenKine.hh"
00025 #include "EvtGenBase/EvtPDL.hh"
00026 #include "EvtGenBase/EvtVector4C.hh"
00027 #include "EvtGenBase/EvtVector4R.hh"
00028 #include "EvtGenBase/EvtTensor4C.hh"
00029 #include "EvtGenBase/EvtReport.hh"
00030 #include "EvtGenBase/EvtdFunction.hh"
00031 #include "EvtGenBase/EvtHelSys.hh"
00032 #include "EvtChi2BB2.hh"
00033 #include <string>
00034 using std::endl;
00035 using std::cout;
00036 EvtChi2BB2::~EvtChi2BB2() {}
00037 
00038 void EvtChi2BB2::getName(std::string& model_name){
00039 
00040   model_name="Chi2BB2";     
00041 
00042 }
00043 
00044 
00045 EvtDecayBase* EvtChi2BB2::clone(){
00046 
00047   return new EvtChi2BB2;
00048 
00049 }
00050 
00051 void EvtChi2BB2::init(){
00052 // check that there are 0 arguments
00053  checkNArg(5); 
00054  checkNDaug(2);
00055  checkSpinParent(EvtSpinType::TENSOR);
00056  checkSpinDaughter(0,EvtSpinType::RARITASCHWINGER);
00057  checkSpinDaughter(1,EvtSpinType::RARITASCHWINGER);
00058 
00059  }
00060 
00061 /*
00062 void EvtChi2BB2::initProbMax() {
00063 
00064   //Hard coded... should not be hard to calculate...
00065  setProbMax(100.0);
00066 
00067 }      
00068 */
00069 void EvtChi2BB2::decay( EvtParticle *p){
00070 
00071   p->initializePhaseSpace(getNDaug(),getDaugs());
00072 
00073   EvtParticle *v,*s1;
00074   EvtVector4R pv,ps,ppr;
00075   
00076   v =p->getDaug(0);  
00077   s1=p->getDaug(1); 
00078   pv=v->getP4();
00079   ps=s1->getP4();
00080   ppr=p->getP4();
00081 
00082 //  Put phase space results into the daughters.
00083   EvtHelSys angles(ppr,pv); //using helicity sys. angles
00084   double theta =angles.getHelAng(1); 
00085   double phi  =angles.getHelAng(2);
00086   double gamma=0;
00087   double F00,F01,F02,F03,F10,F11,F13,F20,F22,F23,F30,F31,F32,F33;
00088   F00=getArg(0);
00089   F01=getArg(1); 
00090   F02=getArg(2);
00091   F03=getArg(3);
00092   F11=getArg(4);
00093 
00094   F33=F00;  F32=F01;  F31=F02; F30=F03; F22=F11;
00095   F10=F01;  F20=F02;  F23=F10; F13=F20; 
00096 //cout<<"********F00,F01,F02,F03,F10,F11,F13,F20,F22,F23,F30,F31,F32,F33="<<F00<<F01<<F02<<F03<<F10<<F11<<F13<<F20<<F22<<F23<<F30<<F31<<F32<<F33<<endl;
00097 
00098 
00099 
00100  vertex(0,0,0,Djmn(2, 0, 0,phi,theta,gamma)*F00);//Chi2  helicity =0 corresponding index 0;
00101  vertex(0,0,1,Djmn(2, 0,-1,phi,theta,gamma)*F01);
00102  vertex(0,0,2,Djmn(2, 0, 2,phi,theta,gamma)*F02);
00103  vertex(0,0,3,Djmn(2, 0, 1,phi,theta,gamma)*F03);
00104  vertex(0,1,0,Djmn(2, 0, 1,phi,theta,gamma)*F10);
00105  vertex(0,1,1,Djmn(2, 0, 0,phi,theta,gamma)*F11);
00106  vertex(0,1,2,0.0);
00107  vertex(0,1,3,Djmn(2, 0, 2,phi,theta,gamma)*F13);
00108  vertex(0,2,0,Djmn(2, 0,-2,phi,theta,gamma)*F20);
00109  vertex(0,2,1,0.0);
00110  vertex(0,2,2,Djmn(2, 0, 0,phi,theta,gamma)*F22);
00111  vertex(0,2,3,Djmn(2, 0,-1,phi,theta,gamma)*F23);
00112  vertex(0,3,0,Djmn(2, 0,-1,phi,theta,gamma)*F30);
00113  vertex(0,3,1,Djmn(2, 0,-2,phi,theta,gamma)*F31);
00114  vertex(0,3,2,Djmn(2, 0, 1,phi,theta,gamma)*F32);
00115  vertex(0,3,3,Djmn(2, 0, 0,phi,theta,gamma)*F33);
00116 
00117  vertex(1,0,0,Djmn(2, 2, 0,phi,theta,gamma)*F00);//Chi2  helicity =2 corresponding index 1;
00118  vertex(1,0,1,Djmn(2, 2,-1,phi,theta,gamma)*F01);
00119  vertex(1,0,2,Djmn(2, 2, 2,phi,theta,gamma)*F02);
00120  vertex(1,0,3,Djmn(2, 2, 1,phi,theta,gamma)*F03);
00121  vertex(1,1,0,Djmn(2, 2, 1,phi,theta,gamma)*F10);
00122  vertex(1,1,1,Djmn(2, 2, 0,phi,theta,gamma)*F11);
00123  vertex(1,1,2,0.0);
00124  vertex(1,1,3,Djmn(2, 2, 2,phi,theta,gamma)*F13);
00125  vertex(1,2,0,Djmn(2, 2,-2,phi,theta,gamma)*F20);
00126  vertex(1,2,1,0.0);
00127  vertex(1,2,2,Djmn(2, 2, 0,phi,theta,gamma)*F22);
00128  vertex(1,2,3,Djmn(2, 2,-1,phi,theta,gamma)*F23);
00129  vertex(1,3,0,Djmn(2, 2,-1,phi,theta,gamma)*F30);
00130  vertex(1,3,1,Djmn(2, 2,-2,phi,theta,gamma)*F31);
00131  vertex(1,3,2,Djmn(2, 2, 1,phi,theta,gamma)*F32);
00132  vertex(1,3,3,Djmn(2, 2, 0,phi,theta,gamma)*F33);
00133 
00134  vertex(2,0,0,Djmn(2,-2, 0,phi,theta,gamma)*F00);//Chi2  helicity =-2 corresponding index 2;
00135  vertex(2,0,1,Djmn(2,-2,-1,phi,theta,gamma)*F01);
00136  vertex(2,0,2,Djmn(2,-2, 2,phi,theta,gamma)*F02);
00137  vertex(2,0,3,Djmn(2,-2, 1,phi,theta,gamma)*F03);
00138  vertex(2,1,0,Djmn(2,-2, 1,phi,theta,gamma)*F10);
00139  vertex(2,1,1,Djmn(2,-2, 0,phi,theta,gamma)*F11);
00140  vertex(2,1,2,0.0); 
00141  vertex(2,1,3,Djmn(2,-2, 2,phi,theta,gamma)*F13);
00142  vertex(2,2,0,Djmn(2,-2,-2,phi,theta,gamma)*F20);
00143  vertex(2,2,1,0.0);
00144  vertex(2,2,2,Djmn(2,-2, 0,phi,theta,gamma)*F22);
00145  vertex(2,2,3,Djmn(2,-2,-1,phi,theta,gamma)*F23);
00146  vertex(2,3,0,Djmn(2,-2,-1,phi,theta,gamma)*F30);
00147  vertex(2,3,1,Djmn(2,-2,-2,phi,theta,gamma)*F31);
00148  vertex(2,3,2,Djmn(2,-2, 1,phi,theta,gamma)*F32);
00149  vertex(2,3,3,Djmn(2,-2, 0,phi,theta,gamma)*F33);
00150  
00151  vertex(3,0,0,Djmn(2, 1, 0,phi,theta,gamma)*F00);//Chi2  helicity =1 corresponding index 3;
00152  vertex(3,0,1,Djmn(2, 1,-1,phi,theta,gamma)*F01);
00153  vertex(3,0,2,Djmn(2, 1, 2,phi,theta,gamma)*F02);
00154  vertex(3,0,3,Djmn(2, 1, 1,phi,theta,gamma)*F03);
00155  vertex(3,1,0,Djmn(2, 1, 1,phi,theta,gamma)*F10);
00156  vertex(3,1,1,Djmn(2, 1, 0,phi,theta,gamma)*F11);
00157  vertex(3,1,2,0.0); 
00158  vertex(3,1,3,Djmn(2, 1, 2,phi,theta,gamma)*F13);
00159  vertex(3,2,0,Djmn(2, 1,-2,phi,theta,gamma)*F20);
00160  vertex(3,2,1,0.0);
00161  vertex(3,2,2,Djmn(2, 1, 0,phi,theta,gamma)*F22);
00162  vertex(3,2,3,Djmn(2, 1,-1,phi,theta,gamma)*F23);
00163  vertex(3,3,0,Djmn(2, 1,-1,phi,theta,gamma)*F30);
00164  vertex(3,3,1,Djmn(2, 1,-2,phi,theta,gamma)*F31);
00165  vertex(3,3,2,Djmn(2, 1, 1,phi,theta,gamma)*F32);
00166  vertex(3,3,3,Djmn(2, 1, 0,phi,theta,gamma)*F33);
00167 
00168  vertex(4,0,0,Djmn(2,-1, 0,phi,theta,gamma)*F00);//Chi2  helicity =-1 corresponding index 4;
00169  vertex(4,0,1,Djmn(2,-1,-1,phi,theta,gamma)*F01);
00170  vertex(4,0,2,Djmn(2,-1, 2,phi,theta,gamma)*F02);
00171  vertex(4,0,3,Djmn(2,-1, 1,phi,theta,gamma)*F03);
00172  vertex(4,1,0,Djmn(2,-1, 1,phi,theta,gamma)*F10);
00173  vertex(4,1,1,Djmn(2,-1, 0,phi,theta,gamma)*F11);
00174  vertex(4,1,2,0.0);
00175  vertex(4,1,3,Djmn(2,-1, 2,phi,theta,gamma)*F13);
00176  vertex(4,2,0,Djmn(2,-1,-2,phi,theta,gamma)*F20);
00177  vertex(4,2,1,0.0);
00178  vertex(4,2,2,Djmn(2,-1, 0,phi,theta,gamma)*F22);
00179  vertex(4,2,3,Djmn(2,-1,-1,phi,theta,gamma)*F23);
00180  vertex(4,3,0,Djmn(2,-1,-1,phi,theta,gamma)*F30);
00181  vertex(4,3,1,Djmn(2,-1,-2,phi,theta,gamma)*F31);
00182  vertex(4,3,2,Djmn(2,-1, 1,phi,theta,gamma)*F32);
00183  vertex(4,3,3,Djmn(2,-1, 0,phi,theta,gamma)*F33);
00184   
00185  return ;
00186 
00187 }
00188 
00189 
00190 
00191