EvtPsi3Sdecay Class Reference

#include <EvtPsi3Sdecay.hh>

List of all members.

Public Member Functions
bool	choseDecay ()
	EvtPsi3Sdecay (double ecms, EvtParticle *parent)
int	findMode ()
int	findPoints ()
double	polint (std::vector< double > points)
double	theProb (std::vector< double > myxs, int ich)
virtual	~EvtPsi3Sdecay ()
Private Attributes
std::vector< double >	d0bardst0
std::vector< double >	d0d0bar
std::vector< double >	d0dst0bar
EvtId *	Daughters
std::vector< double >	dpdm
std::vector< double >	dspdsm
std::vector< double >	dst0dst0bar
std::vector< double >	dstmdp
std::vector< double >	dstpdm
std::vector< double >	dstpdstm
double	Ecms
int	Ndaugs
int	nsize
int	theLocation
EvtParticle *	theParent
double	theXsection [50]
std::vector< double >	x

---

Constructor & Destructor Documentation

EvtPsi3Sdecay::EvtPsi3Sdecay (  double  ecms, EvtParticle *  parent )  [inline]

{ //initializer Ecms = ecms; theParent = parent; Ndaugs=parent->getNDaug(); nsize = 14; x.clear(); // open charm cross section, see PRD 80, 072001 double xx[14]={3.72968, 3.97, 3.99, 4.01, 4.015, 4.03, 4.06, 4.12, 4.14, 4.16, 4.17, 4.18, 4.2, 4.26}; // 14 energy points double y0[14]={0., 86, 133, 76, 10, 334, 410, 303, 177, 167, 177, 179, 180, 86}; // 0) D0D0bar cross section double y1[14]={0., 137, 90, 135, 38, 196, 480, 310, 200, 200, 182, 197, 181, 94}; // 1) D+D- double y2[14]={0., 1140, 1370, 1660, 1920, 1600, 1115, 700, 675, 626, 636, 605, 515, 540}; // 2)D0D*0bar double y3[14]={0., 1140, 1370, 1660, 1920, 1600, 1115, 700, 675, 626, 636, 605, 515, 540}; // 3)D0bar D*0 double y4[14]={0., 0, 0, 0, 213, 2000, 2290, 2550, 2443, 2566, 2363, 2173, 1830, 269}; // 4)D*0 D*0bar double y5[14]={0., 1115, 1375, 1650, 1851, 1650, 1085, 780, 688, 688, 642, 648, 535, 511}; // 5)D*+D- double y6[14]={0., 1115, 1375, 1650, 1851, 1650, 1085, 780, 688, 688, 642, 648, 535, 511}; // 6)D*-D+ double y7[14]={0., 0, 0, 0, 0, 1400, 2390, 2280, 2556, 2479, 2357, 2145, 1564, 237}; // 7)D*+D*- double y8[14]={0., 102, 133, 269, 250, 174, 51, 26, 25, 15, 34, 7, 15, 47}; // 8)Ds+ Ds- d0d0bar.clear(); dpdm.clear(); d0dst0bar.clear(); dst0dst0bar.clear(); d0bardst0.clear(); dstpdm.clear(); dstmdp.clear(); dstpdstm.clear(); dspdsm.clear(); for(int i=0;i<14;i++){ x.push_back(xx[i]); d0d0bar.push_back(y0[i]); dpdm.push_back(y1[i]); d0dst0bar.push_back(y2[i]); d0bardst0.push_back(y3[i]); dst0dst0bar.push_back(y4[i]); dstpdm.push_back( y5[i]); dstmdp.push_back( y6[i]); dstpdstm.push_back(y7[i]); dspdsm.push_back( y8[i]); } }

virtual EvtPsi3Sdecay::~EvtPsi3Sdecay (   )  [inline, virtual]

{}

---

Member Function Documentation

bool EvtPsi3Sdecay::choseDecay (   )

{ //determing accept or reject a generated decay // findPoints(); double d0d0bar_xs=polint(d0d0bar); double dpdm_xs = polint(dpdm); double d0dst0bar_xs = polint(d0dst0bar); double d0bardst0_xs = polint(d0bardst0); double dst0dst0bar_xs = polint(dst0dst0bar); double dstpdm_xs = polint(dstpdm); double dstmdp_xs = polint(dstmdp); double dstpdstm_xs = polint(dstpdstm); double dspdsm_xs = polint(dspdsm); int ich = findMode(); // std::cout<<"calculated XS "<< d0d0bar_xs<<" "<<dpdm_xs<<" "<<d0dst0bar_xs<<" "<<d0bardst0_xs<< " "<<dst0dst0bar_xs<<" "<<dstpdm_xs<< " "<<dstmdp_xs<<" "<<dstpdstm_xs<<" "<<dspdsm_xs<<std::endl; std::vector<double> myxs; myxs.clear(); myxs.push_back(d0d0bar_xs); myxs.push_back(dpdm_xs); myxs.push_back(d0dst0bar_xs); myxs.push_back(d0bardst0_xs); myxs.push_back(dst0dst0bar_xs); myxs.push_back(dstpdm_xs); myxs.push_back(dstmdp_xs); myxs.push_back(dstpdstm_xs); myxs.push_back(dspdsm_xs); double Prop0 = theProb(myxs,ich-1); double Prop1 = theProb(myxs,ich); double pm= EvtRandom::Flat(0.,1); bool flag = false; if( Prop0 < pm && pm<= Prop1 ) flag = true; //--- debuging if(flag) { // std::cout<<"findMode= "<<ich<<std::endl; // for(int i=0;i<myxs.size();i++){ std::cout<<"channel "<<i<<" myxs: "<<myxs[i]<<std::endl;} //std::cout<<"prop0,prop1= "<<Prop0<<" "<<Prop1<<std::endl; } //------------- return flag; }

int EvtPsi3Sdecay::findMode (   )

{ // mode index: 0) D0D0bar, 1)D+D-; 2)D0D*0bar , 3)D0bar D*0, 4)D*0 D*0 bar, 5)D*+D- 6)D*-D+ 7)D*+ D*- 8) Ds+ Ds- std::string son0=EvtPDL::name(theParent->getDaug(0)->getId()); std::string son1=EvtPDL::name(theParent->getDaug(1)->getId()); std::string son2; if(Ndaugs == 3){son2=EvtPDL::name(theParent->getDaug(2)->getId());} if(Ndaugs == 2 || Ndaugs ==3 && son2=="gammaFSR"){ if(son0 == "D0" && son1 == "anti-D0" || son1 == "D0" && son0 == "anti-D0") {return 0;} else if(son0 == "D+" && son1 == "D-" || son1 == "D+" && son0 == "D-" ) {return 1;} else if(son0 == "D0" && son1 == "anti-D*0" || son1 == "D0" && son0 == "anti-D*0") {return 2;} else if(son0 == "anti-D0" && son1 == "D*0" || son1 == "anti-D0" && son0 == "D*0") {return 3;} else if(son0 == "D*0" && son1 == "anti-D*0" || son1 == "D*0" && son0 == "anti-D*0") {return 4;} else if(son0 == "D*+" && son1 == "D-" || son1 == "D*+" && son0 == "D-") {return 5;} else if(son0 == "D*-" && son1 == "D+" || son1 == "D*-" && son0 == "D+") {return 6;} else if(son0 == "D*+" && son1 == "D*-" || son1 == "D*+" && son0 == "D*-") {return 7;} else if(son0 == "D_s+" && son1 == "D_s-" || son1 == "D_s+" && son0 == "D_s-") {return 8;} } else{ std::cout<<"Not open charm decay"<<std::endl; ::abort(); } }

int EvtPsi3Sdecay::findPoints (   )

{ if(Ecms < x[0] ){theLocation =0;} else if(Ecms >=x[nsize]) {theLocation = nsize;} for (int i=0;i<nsize-1;i++){ if( x[i] <= Ecms && x[i+1] > Ecms) {theLocation=i;} } return theLocation; }

double EvtPsi3Sdecay::polint (  std::vector< double >  points  )

{ // EvtPolInt thePolInt(x,vy,Ecms); // double xs = thePolInt.getvalue(); theLocation = findPoints(); double xs = (vy[theLocation +1 ]- vy[theLocation]) / (x[theLocation+1]-x[theLocation])*(Ecms - x[theLocation])+vy[theLocation]; if(xs <0 ) xs = 0; // for(int i=0;i<vy.size();i++){ std::cout<<"channel "<<i<<" xsection: "<<vy[i]<<std::endl;} // std::cout<<"Ecms, theLocation= "<<Ecms<<" "<<theLocation<<std::endl; return xs; }

double EvtPsi3Sdecay::theProb (  std::vector< double >  myxs, int  ich )

{ int Nchannels=myxs.size(); std::vector <double> thexs; thexs.clear(); double sumxs=0; for(int j=0;j<Nchannels;j++){ sumxs += myxs[j]; thexs.push_back(sumxs); } return thexs[ich] / sumxs ; }

---

Member Data Documentation

std::vector<double> EvtPsi3Sdecay::d0bardst0 [private]

std::vector<double> EvtPsi3Sdecay::d0d0bar [private]

std::vector<double> EvtPsi3Sdecay::d0dst0bar [private]

EvtId* EvtPsi3Sdecay::Daughters [private]

std::vector<double> EvtPsi3Sdecay::dpdm [private]

std::vector<double> EvtPsi3Sdecay::dspdsm [private]

std::vector<double> EvtPsi3Sdecay::dst0dst0bar [private]

std::vector<double> EvtPsi3Sdecay::dstmdp [private]

std::vector<double> EvtPsi3Sdecay::dstpdm [private]

std::vector<double> EvtPsi3Sdecay::dstpdstm [private]

double EvtPsi3Sdecay::Ecms [private]

int EvtPsi3Sdecay::Ndaugs [private]

int EvtPsi3Sdecay::nsize [private]

int EvtPsi3Sdecay::theLocation [private]

EvtParticle* EvtPsi3Sdecay::theParent [private]

double EvtPsi3Sdecay::theXsection[50] [private]

std::vector<double> EvtPsi3Sdecay::x [private]

---

The documentation for this class was generated from the following files:

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Generator/BesEvtGen/BesEvtGen-00-01-11/src/EvtGen/EvtGenModels/EvtPsi3Sdecay.hh
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Generator/BesEvtGen/BesEvtGen-00-01-11/src/EvtGen/EvtGenModels/EvtPsi3Sdecay.cc

---