Loading...
Searching...
No Matches
EvtConExc Class Reference
#include <EvtConExc.hh>
Inheritance diagram for EvtConExc:
Public Member Functions | |
| EvtConExc () | |
| virtual | ~EvtConExc () |
| void | getName (std::string &name) |
| EvtDecayBase * | clone () |
| void | initProbMax () |
| void | init () |
| void | init_Br_ee () |
| void | decay (EvtParticle *p) |
| void | init_mode (int mode) |
| double | gamHXSection (EvtParticle *p, double El, double Eh, int nmc=100000) |
| double | gamHXSection (double s, double El, double Eh, int nmc=100000) |
| double | gamHXSection (double El, double Eh) |
| double | gamHXSection (double El, double Eh, int mode) |
| double | gamHXSection_er (double El, double Eh) |
| void | findMaxXS (EvtParticle *p) |
| double | difgamXs (EvtParticle *p) |
| double | difgamXs (double mhds, double sintheta) |
| double | Mhad_sampling (double *x, double *y) |
| double | ISR_ang_integrate (double x, double theta) |
| double | ISR_ang_sampling (double x) |
| bool | hadron_angle_sampling (EvtVector4R ppi, EvtVector4R pcm) |
| void | SetP4 (EvtParticle *part, double mhdr, double xeng, double theta) |
| void | SetP4Rvalue (EvtParticle *part, double mhdr, double xeng, double theta) |
| bool | gam_sampling (EvtParticle *p) |
| bool | xs_sampling (double xs) |
| bool | xs_sampling (double xs, double xs1) |
| bool | baryon_sampling (EvtVector4R pcm, EvtVector4R pi) |
| bool | meson_sampling (EvtVector4R pcm, EvtVector4R pi) |
| bool | VP_sampling (EvtVector4R pcm, EvtVector4R pi) |
| bool | angularSampling (EvtParticle *part) |
| bool | photonSampling (EvtParticle *part) |
| double | baryonAng (double mx) |
| double | Rad1 (double s, double x) |
| double | Rad2 (double s, double x) |
| double | Rad2difXs (EvtParticle *p) |
| double | Rad2difXs (double s, double x) |
| double | Rad2difXs (double s, double x, EvtXsection *myxsection) |
| double | Rad1difXs (EvtParticle *p) |
| double | Ros_xs (double mx, double bree, EvtId pid) |
| double | Li2 (double x) |
| double | SoftPhoton_xs (double s, double b) |
| double | lgr (double *x, double *y, int n, double t) |
| bool | islgr (double *x, double *y, int n, double t) |
| double | LLr (double *x, double *y, int n, double t) |
| int | selectMode (std::vector< int > vmod, double mhds) |
| void | findMaxXS (double mhds) |
| bool | gam_sampling (double mhds, double sintheta) |
| void | resetResMass () |
| void | getResMass () |
| bool | checkdecay (EvtParticle *p) |
| double | sumExc (double mx) |
| void | showResMass () |
| int | getNdaugs () |
| EvtId * | getDaugId () |
| int | getSelectedMode () |
| double | narrowRXS (double mxL, double mxH) |
| double | selectMass () |
| double | addNarrowRXS (double mhi, double binwidth) |
| void | ReadVP () |
| double | getVP (double cms) |
| void | mk_VXS (double Esig, double Egamcut, double EgamH, int midx) |
| int | get_mode_index (int mode) |
| double | getObsXsection (double mhds, int mode) |
| double | Egam2Mhds (double Egam) |
| std::vector< EvtId > | get_mode (int mode) |
| void | writeDecayTabel () |
| void | checkEvtRatio () |
| int | getModeIndex (int m) |
| std::string | commandName () |
| void | command (std::string cmd) |
| std::vector< std::string > | split (std::string str, std::string pattern) |
| void | InitPars () |
| double | energySpread (double mu, double sigma) |
| void | calAF (double myecms) |
| void | OutStaISR () |
| double | trapezoid (double s, double a, double b, int n) |
| double | trapezoid (double s, double El, double Eh, int n, EvtXsection *myxs) |
| Public Member Functions inherited from EvtDecayIncoherent | |
| void | makeDecay (EvtParticle *p) |
| virtual | ~EvtDecayIncoherent () |
| void | setDaughterSpinDensity (int daughter) |
| int | isDaughterSpinDensitySet (int daughter) |
| Public Member Functions inherited from EvtDecayBase | |
| virtual void | getName (std::string &name)=0 |
| virtual void | decay (EvtParticle *p)=0 |
| virtual void | makeDecay (EvtParticle *p)=0 |
| virtual EvtDecayBase * | clone ()=0 |
| virtual void | init () |
| virtual void | initProbMax () |
| virtual std::string | commandName () |
| virtual void | command (std::string cmd) |
| double | getProbMax (double prob) |
| double | resetProbMax (double prob) |
| EvtDecayBase () | |
| virtual | ~EvtDecayBase () |
| virtual bool | matchingDecay (const EvtDecayBase &other) const |
| EvtId | getParentId () |
| double | getBranchingFraction () |
| void | disableCheckQ () |
| void | checkQ () |
| int | getNDaug () |
| EvtId * | getDaugs () |
| EvtId | getDaug (int i) |
| int | getNArg () |
| int | getPHOTOS () |
| void | setPHOTOS () |
| void | setVerbose () |
| void | setSummary () |
| double * | getArgs () |
| std::string * | getArgsStr () |
| double | getArg (int j) |
| std::string | getArgStr (int j) |
| std::string | getModelName () |
| int | getDSum () |
| int | summary () |
| int | verbose () |
| void | saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr) |
| void | printSummary () |
| void | setProbMax (double prbmx) |
| void | noProbMax () |
| void | checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1) |
| void | checkNDaug (int d1, int d2=-1) |
| void | checkSpinParent (EvtSpinType::spintype sp) |
| void | checkSpinDaughter (int d1, EvtSpinType::spintype sp) |
| virtual int | nRealDaughters () |
Static Public Attributes | |
| static EvtXsection * | myxsection |
| static EvtXsection * | staxsection |
| static double | _cms |
| static double | XS_max |
| static double | SetMthr =0 |
| static int | getMode |
| static int | conexcmode =-1 |
| static int | multimode =0 |
| static double | mlow =0 |
| static double | mup =0 |
Additional Inherited Members | |
| Static Public Member Functions inherited from EvtDecayBase | |
| static void | findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10]) |
| static void | findMass (EvtParticle *p) |
| static double | findMaxMass (EvtParticle *p) |
| Protected Member Functions inherited from EvtDecayBase | |
| bool | daugsDecayedByParentModel () |
| Protected Attributes inherited from EvtDecayBase | |
| bool | _daugsDecayedByParentModel |
Detailed Description
Definition at line 138 of file EvtConExc.hh.
Constructor & Destructor Documentation
◆ EvtConExc()
|
inline |
◆ ~EvtConExc()
|
virtual |
Definition at line 205 of file EvtConExc.cc.
205 {
206// if(_mode==74110)checkEvtRatio();
208 if(mydbg){
209 // xs->Write();
210 myfile->Write();
211 }
212
215 gamH->deleteTree();
216
217 //the deletion of commands is really uggly!
218
219 int i;
220 if (nconexcdecays==0) {
221 delete [] commands;
222 commands=0;
223 return;
224 }
225
226 for(i=0;i<nconexcdecays;i++){
227 if (conexcdecays[i]==this){
228 conexcdecays[i]=conexcdecays[nconexcdecays-1];
229 nconexcdecays--;
230 if (nconexcdecays==0) {
231 delete [] commands;
232 commands=0;
233 }
234 return;
235 }
236 }
237
238
239}
Member Function Documentation
◆ addNarrowRXS()
| double EvtConExc::addNarrowRXS | ( | double | mhi, |
| double | binwidth | ||
| ) |
Definition at line 3326 of file EvtConExc.cc.
3326 {
3327 double myxs = 0;
3328 for(int i=0;i<ISRXS.size();i++){ // std::cout<<"ISRM "<<ISRM[i]<<std::endl;
3329 if( (ISRM[i]-mhi)<binwidth && ISRFLAG[i]){ISRFLAG[i]=0; myxs = ISRXS[i];}
3330 }
3331 //std::cout<<mhi<<" "<<binwidth<<" "<<myxs<<std::endl;
3332 return myxs;
3333}
◆ angularSampling()
| bool EvtConExc::angularSampling | ( | EvtParticle * | part | ) |
Definition at line 3215 of file EvtConExc.cc.
3215 {
3216 bool tagp,tagk;
3217 tagk=0;
3218 tagp=0;
3219 int nds = par->getNDaug();
3220 for(int i=0;i<par->getNDaug();i++){
3222 EvtVector4R p4i = par->getDaug(i)->getP4Lab();
3225 /*
3226 if(fabs(pdgcode)==2212){//proton and anti-proton
3227 alpha = baryonAng(p4i.mass());
3228 cosp = cos(p4i.theta());
3229 tagp=1;
3230 }else if(fabs(pdgcode)==321||fabs(pdgcode)==211){
3231 alpha=1;
3232 cosk = cos(p4i.theta());
3233 tagk=1;
3234 }else continue;
3235 */
3239 double xratio = ang/angmax;
3241 //std::cout<<"pdgcode "<<pdgcode<<std::endl;
3242 //std::cout<<ang<<" "<<angmax<<" "<<xi<<" "<<xratio<<std::endl;
3243 if(xi>xratio) return false;
3244 }//loop over duaghters
3245 //if(tagp) std::cout<<" cosp "<<cosp<<std::endl;
3246 //if(tagk) std::cout<<" cosk "<<cosk<<std::endl;
3247
3248 return true;
3249}
Definition: EvtId.hh:27
Definition: EvtVector4R.hh:29
Referenced by decay().
◆ baryon_sampling()
| bool EvtConExc::baryon_sampling | ( | EvtVector4R | pcm, |
| EvtVector4R | pi | ||
| ) |
Definition at line 2529 of file EvtConExc.cc.
2529 {
2531 double theta=angles.getHelAng(1);
2532 double phi =angles.getHelAng(2);
2534
2535 // double costh = pcm.dot(pi)/pcm.d3mag()/pi.d3mag();
2536 //std::cout<<"costheta: "<<costheta<<", "<<costh<<std::endl;
2541 double myang;
2543 if(pm< ang/myang) {return true;}else{return false;}
2544}
Definition: EvtHelSys.hh:44
Referenced by hadron_angle_sampling().
◆ baryonAng()
| double EvtConExc::baryonAng | ( | double | mx | ) |
Definition at line 3251 of file EvtConExc.cc.
3251 {
3253 double u = 0.938/mx;
3254 u = u*u;
3255 double u2 = (1+u)*(1+u);
3256 double uu = u*(1+6*u);
3259}
Referenced by baryon_sampling().
◆ calAF()
| void EvtConExc::calAF | ( | double | myecms | ) |
Definition at line 3687 of file EvtConExc.cc.
3687 {
3688
3689
3691 double Esig = 0.0001; //sigularity engy
3694
3697
3698 int nint=50;
3699 int nmc= 1000;
3701 //--- sigularity point x from 0 to 0.0001GeV
3704 _xs0 += sig_xs;
3705
3706 int Nstp=598;
3707 double stp=(EgamH - Egamcut)/Nstp;
3708 for(int i=0;i<Nstp;i++){//points within Egam(max) ~ Egamcut=25MeV
3709 double Eg0=EgamH - i*stp;
3710 double Eg1=EgamH - (i+1)*stp;
3712
3713 AA[i]=(Eg1+Eg0)/2;
3716 double binwidth = mh2-mhi;
3717 if(i==0) AF[0]=xsi;
3718 if(i>=1) AF[i]=AF[i-1]+xsi;
3719 RadXS[i]=xsi/stp;
3720 }
3721 //add the interval 0~Esig
3722 AA[598]=Egamcut; AA[599]=0; //AA[i]: E_gamma
3723 AF[598]=AF[597];
3724 AF[599]=AF[598]+ _xs0;
3725 //--
3726 for(int i=0;i<600;i++){
3728 }
3729 //--debugging
3732 double fisr=AF[599]/bornXS;
3733 myFisr.push_back(fisr);
3734 //std::cout<<"fisr= "<<fisr<<std::endl;
3735}
double trapezoid(double s, double a, double b, int n)
Definition: EvtConExc.cc:3752
Referenced by decay().
◆ checkdecay()
| bool EvtConExc::checkdecay | ( | EvtParticle * | p | ) |
Definition at line 3165 of file EvtConExc.cc.
3165 {
3167 double massflag=1;
3168 for(int i=0;i<nson;i++){
3171 }
3172 if(massflag==0){
3173 std::cout<<"Zero mass!"<<std::endl;
3174 return 0;
3175 }else{return 1;}
3176}
Referenced by decay().
◆ checkEvtRatio()
| void EvtConExc::checkEvtRatio | ( | ) |
Definition at line 3553 of file EvtConExc.cc.
3553 {
3556 //
3558 double xscon=VXS[im][599];
3559 double xscon0=xscon;
3563 for(int i=0;i<vmode.size();i++){
3564 //std::cout<<"i: "<<xscon<<" "<<VXS[Vmode[i]][599]<<std::endl;
3566 if(vmode[i]==74110) continue;
3567 xscon -= VXS[themode ][599];
3568 }
3569 if(xscon<0) xscon=0;
3570 //debugging
3571 for(int i=0;i<vmode.size();i++){
3573 if(vmode[i]==74110) continue;
3575 }
3577
3578
3579}
std::ofstream ofstream
◆ clone()
|
virtual |
◆ command()
|
virtual |
Reimplemented from EvtDecayBase.
Definition at line 3595 of file EvtConExc.cc.
3595 {
3596 if (ncommand==lcommand){
3597 lcommand=10+2*lcommand;
3598 std::string* newcommands=new std::string[lcommand];
3599 int i;
3600 for(i=0;i<ncommand;i++){
3601 newcommands[i]=commands[i];
3602 }
3603 delete [] commands;
3604 commands=newcommands;
3605 }
3606 commands[ncommand]=cmd;
3607 ncommand++;
3608}
◆ commandName()
|
virtual |
Reimplemented from EvtDecayBase.
Definition at line 3589 of file EvtConExc.cc.
3589 {
3590
3591 return std::string("ConExcPar");
3592
3593}
◆ decay()
|
virtual |
Implements EvtDecayBase.
Definition at line 1763 of file EvtConExc.cc.
1763 {
1764 //std::cout<<"_cms= "<<_cms<<" mode= "<<_mode<<std::endl;
1768 abort();
1769 }
1770 //-- to read _cms from database or to consider beam energy spread
1772 bool beamflag=0;
1774 EvtConExc::_cms=mvpho;
1775 beamflag=1;
1776 }
1777
1778 /*
1779 if(beamEnergySpread!=0){
1780 cmsspread=sqrt(2.)*beamEnergySpread;
1781 _cms = energySpread(_cms,cmsspread);
1782 beamflag=1;
1783 std::cout<<"test_cms "<<_cms<<std::endl;
1784 }
1785 */
1786
1788
1789
1790 //debugging
1791 //double mvpho=EvtPDL::getMeanMass(myvpho);
1792 //std::cout<<"myvpho= "<<mvpho<<std::endl;
1793
1794 //for fill root tree
1795 EvtVector4R vgam,hd1,hd2,vhds;
1796
1797 //first for Rscan generator with _mode=74110
1798 if(_mode==74110){
1799 //preparation of mode sampling
1800 std::vector<int> vmod; vmod.clear();
1801 int mn[83]={0,1,2,3,4,5,6,7,8,9,10,11,13,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,32,33,34,35,36,37,38,40,41,44,45,46, // 12,14, 30,31,39,42 and 43 is removed
1802 50,51,
1803 59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,
1804 90,91,92,93,94,95,96,
1805 74100,74101,74102,74103,74104,74105,74110};
1806 int mn2[84]={0,1,2,3,4,5,6,7,8,9,10,11,12,13,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,32,33,34,35,36,37,38,40,41,44,45,46,//mode 14,30,31,42, 43 are removed
1807 50,51,
1808 59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,
1809 90,91,92,93,94,95,96,
1810 74100,74101,74102,74103,74104,74105,74110};
1811 int mn3[77]={0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,44,45,46,
1812 50,51,
1813 68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,
1814 90,91,92,93,94,95,96,
1815 74100,74101,74102,74110}; //remove 43, 74103,74104,74105, this is included in PHOKHARA
1817
1819 for(int i=0;i<83;i++){vmod.push_back(mn[i]);}
1820 }else{
1821 for(int i=0;i<84;i++){vmod.push_back(mn2[i]);}
1822 }
1823
1824 //for(int i=0;i<76;i++){vmod.push_back(mn3[i]);}
1825
1826 int mymode;
1828
1829 //std::cout<<"AF[598], AF[598] "<<AF[598]<<" "<<AF[599]<<std::endl;
1830 //--
1831 //-- where AF[599]: cross section with full region (0-Egamma^max GeV), AF[598]: cross section within Egam=(0.025-Egamma^max GeV)
1832 if(pm <_xs0/(_xs0 + _xs1)){//without ISR photon for mode=74110
1833
1834 mhds=_cms;
1835 mymode=selectMode(vmod,mhds);
1836 _selectedMode = mymode;
1837 std::cout<<"A selected mode is "<<mymode<<" with Mhds= "<<mhds<<std::endl; //debugging
1840 init_mode(mymode);
1841 resetResMass();
1842
1843 if(_ndaugs>1){//for e+e- -> exclusive decays
1844 checkA:
1845 p->makeDaughters(_ndaugs,daugs);
1846 double totMass=0;
1847 for(int i=0;i< _ndaugs;i++){
1850 di->setMass(xmi);
1851 totMass+=xmi;
1852 }
1854 p->initializePhaseSpace( _ndaugs , daugs);
1856 vhds = p->getP4();
1859 }else{// for e+e- -> vector resonace, add a very soft photon
1860 EvtId mydaugs[2];
1861 mydaugs[0]=daugs[0];
1862 EvtPDL::reSetMass(mydaugs[0],mhds-0.001);
1863 //EvtPDL::reSetWidth(mydaugs[0],0);
1864 mydaugs[1]=gamId; //ISR photon
1865 p->makeDaughters(2,mydaugs);
1866 checkB:
1867 double totMass=0;
1868 for(int i=0;i< 2;i++){
1871 di->setMass(xmi);
1872 totMass+=xmi;
1873 }
1874 //std::cout<<"checkB "<<totMass<<" p->mass() "<<p->mass()<<" "<<checkdecay(p)<<std::endl;
1876 p->initializePhaseSpace(2,mydaugs);
1877
1881 }
1882 //-----
1883 }else{// with ISR photon for mode=74110
1884Sampling_mhds:
1885 mhds=Mhad_sampling(MH,AF);
1886 //std::cout<<"SetMthr= "<<SetMthr<<std::endl;
1890
1891 if(mydbg) mass2=mhds;
1892
1893 //generating events
1894 ModeSelection:
1895 EvtGlobalSet::ConExcPythia=1;
1896 mymode=selectMode(vmod,mhds);
1897 if(mymode==-10) goto Sampling_mhds;
1899 if(mhds<2.3 && mymode==74110) {goto ModeSelection;} // E<2.3 GeV, fully exclusive production
1900 std::cout<<"A selected mode is "<<mymode<<" with Mhds= "<<mhds<<std::endl; //debugging
1901 _selectedMode = mymode;
1904 init_mode(mymode);
1907 EvtPDL::reSetWidth(myvpho,0);
1908
1909 //debugging
1910 //for(int i=0;i<_ndaugs+1;i++){std::cout<<_ndaugs<<" "<<EvtPDL::name(daugs[i])<<std::endl;}
1911
1912 //decay e+e- ->gamma_ISR gamma*
1913 EvtId mydaugs[2];
1914 if(_ndaugs>1) { //gamma* -> exclusive decays
1915 resetResMass();
1916 mydaugs[0]=myvpho;
1917 }else{// vector resonance decays
1918 resetResMass();
1919 mydaugs[0]=daugs[0];
1920 EvtPDL::reSetMass(mydaugs[0],mhds);
1921 //EvtPDL::reSetWidth(mydaugs[0],0);
1922 }
1923 mydaugs[1]=gamId; //ISR photon
1924 int maxflag=0;
1925 int trycount=0;
1926 ISRphotonAng_sampling:
1927 double totMass=0;
1928 p->makeDaughters(2,mydaugs);
1929 for(int i=0;i< 2;i++){
1932 di->setMass(xmi);
1933 totMass+=xmi;
1934 }
1936 //std::cout<<"Ndaugs= "<<_ndaugs<<" Mhds= "<<mhds<<std::endl;
1939 //set the ISR photon angular distribution
1941
1949 trycount++;
1950
1951 } // with and without ISR sampling block
1952 // finish event generation
1953 // for debugging
1954 if(mydbg){
1955 pgam[0]=vgam.get(0);
1956 pgam[1]=vgam.get(1);
1957 pgam[2]=vgam.get(2);
1958 pgam[3]=vgam.get(3);
1959
1960 phds[0]=vhds.get(0);
1961 phds[1]=vhds.get(1);
1962 phds[2]=vhds.get(2);
1963 phds[3]=vhds.get(3);
1965 selectmode = mymode;
1966 xs->Fill();
1967 //std::cout<<"mass1= "<<mass1<<" mass2= "<<mass2<<std::endl;
1968 }
1969 _modeFlag[1]= _selectedMode;
1970 p->setIntFlag(_modeFlag);
1971 return;
1972 }//end block if(_mode==74110)
1973
1974 // ======================================================
1975 // multi-exclusive mode
1976 //=======================================================
1977 if(_mode==-100){
1978 int mymode;
1980
1981 //std::cout<<"AF[598], AF[598] "<<AF[598]<<" "<<AF[599]<<std::endl;
1982 //--
1983 //-- where AF[599]: cross section with full region (0-Egamma^max GeV), AF[598]: cross section within Egam=(0.025-Egamma^max GeV)
1984 if(pm <_xs0/(_xs0 + _xs1)){//without ISR photon for mode=74110
1985
1986 mhds=_cms;
1987 mymode=selectMode(vmd,mhds);
1988 _selectedMode = mymode;
1989 std::cout<<"A selected mode is "<<mymode<<" with Mhds= "<<mhds<<std::endl; //debugging
1992 init_mode(mymode);
1993 resetResMass();
1994
1995 if(_ndaugs>1){//for e+e- -> exclusive decays
1996 checkAA:
1997 p->makeDaughters(_ndaugs,daugs);
1998 double totMass=0;
1999 for(int i=0;i< _ndaugs;i++){
2002 di->setMass(xmi);
2003 totMass+=xmi;
2004 }
2006 p->initializePhaseSpace( _ndaugs , daugs);
2008 vhds = p->getP4();
2009 bool byn_ang;
2012 //p->printTree();
2013 if(_ndaugs==2){//angular distribution for the hadrons only for two-body decays
2014 byn_ang=hadron_angle_sampling(pd1, pcm);
2015 if(!byn_ang) goto checkAA;
2016 }
2017
2018 }else{// for e+e- -> vector resonace, add a very soft photon
2019 EvtId mydaugs[2];
2020 mydaugs[0]=daugs[0];
2021 EvtPDL::reSetMass(mydaugs[0],mhds-0.001);
2022 //EvtPDL::reSetWidth(mydaugs[0],0);
2023 mydaugs[1]=gamId; //ISR photon
2024 p->makeDaughters(2,mydaugs);
2025 checkBA:
2026 double totMass=0;
2027 for(int i=0;i< 2;i++){
2030 di->setMass(xmi);
2031 totMass+=xmi;
2032 }
2033 //std::cout<<"checkB "<<totMass<<" p->mass() "<<p->mass()<<" "<<checkdecay(p)<<std::endl;
2035 p->initializePhaseSpace(2,mydaugs);
2036
2040 }
2041 //-----
2042 }else{// with ISR photon for mode=-100
2043Sampling_mhds_A:
2044 mhds=Mhad_sampling(MH,AF);
2045 //std::cout<<"SetMthr= "<<SetMthr<<std::endl;
2049
2050 if(mydbg) mass2=mhds;
2051
2052 //generating events
2053 ModeSelection_A:
2054 mymode=selectMode(vmd,mhds);
2055 if(mymode==-10) goto Sampling_mhds_A;
2057 std::cout<<"A selected mode is "<<mymode<<" with Mhds= "<<mhds<<std::endl; //debugging
2058 _selectedMode = mymode;
2061 init_mode(mymode);
2064 EvtPDL::reSetWidth(myvpho,0);
2065
2066 //debugging
2067 // for(int i=0;i<_ndaugs+1;i++){std::cout<<_ndaugs<<" "<<EvtPDL::name(daugs[i])<<std::endl;}
2068
2069 //decay e+e- ->gamma_ISR gamma*
2070 EvtId mydaugs[2];
2071 if(_ndaugs>1) { //gamma* -> exclusive decays
2072 resetResMass();
2073 mydaugs[0]=myvpho;
2074 }else{// vector resonance decays
2075 resetResMass();
2076 mydaugs[0]=daugs[0];
2077 EvtPDL::reSetMass(mydaugs[0],mhds);
2078 //EvtPDL::reSetWidth(mydaugs[0],0);
2079 }
2080 mydaugs[1]=gamId; //ISR photon
2081 int maxflag=0;
2082 int trycount=0;
2083 ISRphotonAng_sampling_A:
2084 double totMass=0;
2085 p->makeDaughters(2,mydaugs);
2086 for(int i=0;i< 2;i++){
2089 di->setMass(xmi);
2090 totMass+=xmi;
2091 }
2093 //std::cout<<"Ndaugs= "<<_ndaugs<<" Mhds= "<<mhds<<std::endl;
2096 //set the ISR photon angular distribution
2098
2099 //--debugging
2100 //p->printTree();
2101
2109 trycount++;
2110
2111 } // with and without ISR sampling block
2112 _modeFlag[0]=-100;
2113 _modeFlag[1]= _selectedMode;
2114 p->setIntFlag(_modeFlag);
2115 // finish event generation
2116 // for debugging
2117 if(mydbg){
2118 pgam[0]=vgam.get(0);
2119 pgam[1]=vgam.get(1);
2120 pgam[2]=vgam.get(2);
2121 pgam[3]=vgam.get(3);
2122
2123 phds[0]=vhds.get(0);
2124 phds[1]=vhds.get(1);
2125 phds[2]=vhds.get(2);
2126 phds[3]=vhds.get(3);
2128 selectmode = mymode;
2129 xs->Fill();
2130 //std::cout<<"mass1= "<<mass1<<" mass2= "<<mass2<<std::endl;
2131 }
2132 return;
2133 }//end block if(_mode==-100)
2134
2135 //=======================================================
2136 // e+e- -> gamma_ISR gamma*
2137 //=======================================================
2138 if(VISR){
2141
2142 if(pm <_xs0/(_xs0 + _xs1)){//with a very soft ISR photon
2144 mhds = _cms-0.0001;
2145 }else{
2146 mhds=Mhad_sampling(MH,AF);
2147 EvtPDL::reSetMass(gid,mhds);
2148 }
2149 //debugging
2153 p->makeDaughters(2,daugs);
2154 for(int i=0;i< 2;i++){
2156 //if(EvtPDL::name(di->getId())=="gamma*") di->setVectorSpinDensity();
2158 di->setMass(xmi);
2159 }
2160 p->initializePhaseSpace(2,daugs);
2161 SetP4(p,mhds,xeng,theta);
2162 return;
2163 }
2164
2165
2166 //========================================================
2167 //=== for other mode generation with _mode != 74110
2168 //========================================================
2169
2170 if((_xs0 + _xs1)==0) {std::cout<<"EvtConExc::zero total xsection"<<std::endl;::abort();}
2172 double xsratio = _xs0/(_xs0 + _xs1);
2173 int iflag=2; //flag = 0 for ee->hadrons, 1 for ee->gamma_ISR hadrons, 2: mix 0 and 1 case
2177 }
2178
2179 // std::cout<<"xsratio= "<<xsratio<<std::endl;
2180
2181 if(pm<xsratio){// no ISR photon
2182 masscheck:
2183 double summassx=0;
2184 p->makeDaughters(_ndaugs,daugs);
2185 for(int i=0;i< _ndaugs;i++){
2188 di->setMass(xmi);
2189 summassx += xmi;
2190 //std::cout<<"PID and mass "<<di->getId()<<" "<<xmi<<std::endl;
2191 }
2193 angular_hadron:
2194 p->initializePhaseSpace(_ndaugs,daugs);
2195 bool byn_ang;
2198 if(_ndaugs==2){//angular distribution for the hadrons only for two-body decays
2199 byn_ang=hadron_angle_sampling(pd1, pcm);
2200 if(!byn_ang) goto angular_hadron;
2201 }
2202
2203 //for histogram
2205 mhds = _cms;
2206 //std::cout<<"cosp, mhds "<<cosp<<" "<<mhds<<std::endl;
2207 //p->printTree();
2208
2209 }else{// sampling m_hadrons and decay e+e- ->gamma gamma*
2214 EvtPDL::reSetMass(gid,mhdr);
2215 int ndaugs =2;
2216 EvtId mydaugs[2];
2219
2220 //for histogram
2221 mhds = mhdr;
2222 costheta = cos(theta);
2223 //--
2224
2225 p->makeDaughters(2,mydaugs);
2226 for(int i=0;i< 2;i++){
2229 //if(EvtPDL::name(di->getId())=="vhdr") di->setVectorSpinDensity();
2230 di->setMass(xmi);
2231 }
2232 p->initializePhaseSpace(2,mydaugs);
2234 //p->printTree();
2235
2236 //----
2237 }//end of gamma gamma*, gamma*->hadrons generation
2238 // p->printTree();
2239
2240 //-----------------------------------------
2241 // End of decays for all topology
2242 //----------------------------------------
2243 //================================= event analysis
2244 if(_nevt ==0){
2245 std::cout<<"The decay chain: "<<std::endl;
2246 p->printTree();
2247 }
2248 _nevt++;
2249 //--- for debuggint to make root file
2250 if(mydbg){
2251 xs->Fill();
2252 }
2253
2254 //----
2255 return ;
2256}
void SetP4Rvalue(EvtParticle *part, double mhdr, double xeng, double theta)
Definition: EvtConExc.cc:2959
int selectMode(std::vector< int > vmod, double mhds)
Definition: EvtConExc.cc:3003
bool hadron_angle_sampling(EvtVector4R ppi, EvtVector4R pcm)
Definition: EvtConExc.cc:2258
void SetP4(EvtParticle *part, double mhdr, double xeng, double theta)
Definition: EvtConExc.cc:2935
Definition: EvtParticle.hh:42
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1071
Definition: EvtXsection.hh:132
◆ difgamXs() [1/2]
| double EvtConExc::difgamXs | ( | double | mhds, |
| double | sintheta | ||
| ) |
Definition at line 2992 of file EvtConExc.cc.
2992 {// dsigma/dm/dcos\theta* (Eq(1)@PRD73,012005)
2994 double sin2theta = sintheta*sintheta;
3000 return difxs;
3001}
◆ difgamXs() [2/2]
| double EvtConExc::difgamXs | ( | EvtParticle * | p | ) |
Definition at line 2476 of file EvtConExc.cc.
2476 {// dsigma/dm/dcos\theta* (Eq(1)@PRD73,012005)
2477 // where p -> hdarons(0 ~ _ndaugs-1) + gamma
2479 for(int i=1;i<_ndaugs;i++){
2481 }
2486 sin2theta = 1-sin2theta*sin2theta;
2487
2493 // std::cout<<"pgam"<<pgam<<" mhs,egam,sin2theta,cms,x,wsx"<< mhs<<", "<<egam<<", "<<sin2theta<<", "<<cms<<", "<<x<<", "<<wsx<<std::endl;
2494
2496 double difxs = 2*mhs/cms/cms * wsx *xs;
2498 return difxs;
2499
2500}
Referenced by findMaxXS(), gam_sampling(), and gamHXSection().
◆ Egam2Mhds()
| double EvtConExc::Egam2Mhds | ( | double | Egam | ) |
◆ energySpread()
| double EvtConExc::energySpread | ( | double | mu, |
| double | sigma | ||
| ) |
Definition at line 3669 of file EvtConExc.cc.
3669 {
3670 //mu: mean value in Gaussian
3671 //sigma: variance in Gaussian
3672 rloop:
3674 double ri=0;
3677 ri += pm;
3678 }
3680 double xsig= sigma*eta+mu;//smapling Gaussion value
3683 if(xsig<mx0 || xsig>mx1) goto rloop;
3684 return xsig;
3685}
◆ findMaxXS() [1/2]
| void EvtConExc::findMaxXS | ( | double | mhds | ) |
Definition at line 2977 of file EvtConExc.cc.
2977 {// dsigma/dm/dcos\theta* (Eq(1)@PRD73,012005)
2978 maxXS=-1;
2979 for(int i=0;i<90000;i++){
2988 if(difxs>maxXS)maxXS=difxs;
2989 }
2990}
◆ findMaxXS() [2/2]
| void EvtConExc::findMaxXS | ( | EvtParticle * | p | ) |
Definition at line 2431 of file EvtConExc.cc.
2431 {
2432 //std::cout<<"nmc= "<<nmc<<std::endl;
2436 double totxs = 0;
2437 maxXS=-1;
2438 _er1=0;
2439 Rad2Xs =0;
2440 int nmc = 50000;
2441 for(int ii=0;ii<nmc;ii++){// integarted the gamma hadrons xsection
2442 gamH=EvtParticleFactory::particleFactory(xvp, p_init);
2443 int gamdaugs = _ndaugs+1;
2444
2445 EvtId theDaugs[10];
2446 for(int i=0;i<_ndaugs;i++){
2447 theDaugs[i] = daugs[i];
2448 }
2449 theDaugs[_ndaugs]=gamId;
2450
2451 gamH->makeDaughters(gamdaugs,theDaugs);
2452 loop:
2453 double totm=0;
2454 for(int i=0;i<gamdaugs;i++){
2457 di->setMass(xmi);
2458 totm += xmi;
2459 }
2460
2461 //std::cout<<totm<<" "<<p_init.get(0)<<std::endl;
2462 if(totm >= p_init.get(0) ) goto loop;
2463
2468 double sintheta = sqrt(1-costheta*costheta);
2469 bool acut=(sintheta>0.04) && (sintheta<0.96); //set photon auglar cut 2^o
2470 if(thexs>maxXS && acut ) {maxXS=thexs;}
2471 //gamH->deleteTree();
2472 }
2473
2474}
*********Class see also m_nmax DOUBLE PRECISION m_MasPhot DOUBLE PRECISION m_phsu DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_r2 DOUBLE PRECISION m_WtMass INTEGER m_nmax INTEGER m_Nevgen INTEGER m_IsFSR INTEGER m_MarTot *COMMON c_KarFin $ !Output file $ !Event serial number $ !alpha QED at Thomson limit $ !minimum energy at CMS for remooval $ !infrared dimensionless $ !dummy photon IR regulator $ !crude photon multiplicity enhancement factor *EVENT $ !MC crude volume of PhhSpace *Sfactors $ !YFS formfactor IR part only $ !YFS formfactor non IR finite part $ !mass weight
Definition: KarFin.h:34
static EvtParticle * particleFactory(EvtSpinType::spintype spinType)
Definition: EvtParticleFactory.cc:46
◆ gam_sampling() [1/2]
| bool EvtConExc::gam_sampling | ( | double | mhds, |
| double | sintheta | ||
| ) |
Definition at line 2510 of file EvtConExc.cc.
2510 {
2513 double xsratio = xs/maxXS;
2514 if(pm<xsratio){return true;}else{return false;}
2515}
◆ gam_sampling() [2/2]
| bool EvtConExc::gam_sampling | ( | EvtParticle * | p | ) |
Definition at line 2503 of file EvtConExc.cc.
2503 {//photon angular distribution sampling
2506 double xsratio = xs/maxXS;
2507 if(pm<xsratio){return true;}else{return false;}
2508}
Referenced by decay().
◆ gamHXSection() [1/4]
| double EvtConExc::gamHXSection | ( | double | El, |
| double | Eh | ||
| ) |
Definition at line 2377 of file EvtConExc.cc.
2377 {// using Gaussian integration subroutine from Cern lib
2378 std::cout<<" "<<std::endl;
2382 double Rad2Xs;
2383 if(_mode!=-2) {Rad2Xs= dgauss_(&Rad2difXs,&El,&Eh,&eps);}else{Rad2Xs = dgauss_(&Rad2difXs2,&El,&Eh,&eps);}
2384 if(Rad2Xs==0){
2386 eps += 0.01;
2387 if(_mode!=-2) {Rad2Xs= dgauss_(&Rad2difXs,&El,&Eh,&eps);}else{Rad2Xs = dgauss_(&Rad2difXs2,&El,&Eh,&eps);}
2388 if(!Rad2Xs) return Rad2Xs;
2389 }
2390 }
2391 return Rad2Xs; // the leading second order correction
2392}
double dgauss_(double(*fun)(double *), double *, double *, double *)
EvtStreamInputIterator< typename Generator::result_type > iter(Generator gen, int N=0)
Definition: EvtStreamInputIterator.hh:100
◆ gamHXSection() [2/4]
| double EvtConExc::gamHXSection | ( | double | El, |
| double | Eh, | ||
| int | mode | ||
| ) |
Definition at line 2395 of file EvtConExc.cc.
2395 {// using Gaussian integration subroutine from Cern lib
2396 std::cout<<" "<<std::endl;
2400 double Rad2Xs;
2401 if(mode!=-2) {Rad2Xs= dgauss_(&Rad2difXs,&El,&Eh,&eps);}else{Rad2Xs = dgauss_(&Rad2difXs2,&El,&Eh,&eps);}
2402 if(Rad2Xs==0){
2404 eps += 0.01;
2405 if(mode!=-2) {Rad2Xs= dgauss_(&Rad2difXs,&El,&Eh,&eps);}else{Rad2Xs = dgauss_(&Rad2difXs2,&El,&Eh,&eps);}
2406 if(!Rad2Xs) return Rad2Xs;
2407 }
2408 }
2409 return Rad2Xs; // the leading second order correction
2410}
◆ gamHXSection() [3/4]
| double EvtConExc::gamHXSection | ( | double | s, |
| double | El, | ||
| double | Eh, | ||
| int | nmc = 100000 |
||
| ) |
Definition at line 2345 of file EvtConExc.cc.
2345 {//El, Eh : the lower and higher energy of radiative photons
2346 //--for narrow resonance psi(2S),J/psi, phi, which can not calculated with the MC integration
2347 //double mxL = sqrt( s-2*Eh*sqrt(s) ); //low mass
2348 //double mxH = sqrt( s-2*El*sqrt(s) ); //high mass
2349 //double sigv = narrowRXS(mxL,mxH);
2350 //--
2351
2352 double totxs = 0;
2353 maxXS=-1;
2354 _er1=0;
2355 Rad2Xs =0;
2358 for(int i=0;i<nmc;i++){//It is found that the narrow resonance can not included in this MC integartion
2362 _er1 += differ2; //stored when calculate Rad2Xs
2363 // std::cout<<"i= "<<i<<" Rad2Xs= "<<Rad2Xs<<std::endl;
2364 }
2365 _er1/=nmc;
2366 _er1*=(xH-xL);
2367 //std::cout<<"_er1= "<<_er1<<std::endl;
2368 Rad2Xs/=nmc; // the leading second order correction
2369 double thexs= Rad2Xs*(xH-xL); //xH-xL is the Jacobi factor
2370 //std::cout<<"thexs= "<<thexs<<std::endl;
2371 //if(sigv != -1) thexs += sigv; //add narrow resonance ISR cross section
2372 return thexs;
2373}
◆ gamHXSection() [4/4]
| double EvtConExc::gamHXSection | ( | EvtParticle * | p, |
| double | El, | ||
| double | Eh, | ||
| int | nmc = 100000 |
||
| ) |
Definition at line 2286 of file EvtConExc.cc.
2286 {//El, Eh : the lower and higher energy of radiative photons
2287 //std::cout<<"nmc= "<<nmc<<std::endl;
2291 double totxs = 0;
2292 maxXS=-1;
2293 _er1=0;
2294 Rad2Xs =0;
2295 for(int ii=0;ii<nmc;ii++){// integarted the gamma hadrons xsection
2296 gamH=EvtParticleFactory::particleFactory(xvp, p_init);
2297 int gamdaugs = _ndaugs+1;
2298
2299 EvtId theDaugs[10];
2300 for(int i=0;i<_ndaugs;i++){
2301 theDaugs[i] = daugs[i];
2302 }
2303 theDaugs[_ndaugs]=gamId;
2304
2305 gamH->makeDaughters(gamdaugs,theDaugs);
2306
2307 for(int i=0;i<gamdaugs;i++){
2310 di->setMass(xmi);
2311 }
2312 loop:
2314 //-- slection:theta_gamma > 1 degree
2318 double sintheta = sqrt( pmag*pmag - pz*pz)/pmag;
2319 double onedegree = 1./180.*3.1415926;
2320 //if(sintheta < onedegree) {goto loop;}
2321 if(pmag <El || pmag >Eh) {goto loop;}
2322 //--------
2323 // std::cout<<"pmag= "<<pmag<<std::endl;
2324
2326 Rad2Xs += Rad2difXs( gamH );
2327 if(thexs>maxXS) {maxXS=thexs;}
2331 totxs += rad1xs;
2332 _er1 += differ;
2333 gamH->deleteDaughters();
2334 } //for int_i block
2335 _er1/=nmc;
2336
2337 Rad2Xs/=nmc; // the leading second order correction
2338 totxs/=nmc; // first order correction XS
2339
2340 // return totxs; // first order correction XS
2341 return Rad2Xs; // the leading second order correction
2342}
void deleteDaughters(bool keepChannel=false)
Definition: EvtParticle.cc:540
Referenced by init().
◆ gamHXSection_er()
| double EvtConExc::gamHXSection_er | ( | double | El, |
| double | Eh | ||
| ) |
Definition at line 2413 of file EvtConExc.cc.
2413 {// using Gaussian integration subroutine from Cern lib
2414 std::cout<<" "<<std::endl;
2418 double Rad2Xs;
2419 if(_mode !=-2){Rad2Xs = dgauss_(&Rad2difXs_er,&El,&Eh,&eps);}else{ Rad2Xs = dgauss_(&Rad2difXs_er2,&El,&Eh,&eps);}
2420 if(Rad2Xs==0){
2422 eps += 0.01;
2423 if(_mode !=-2){Rad2Xs = dgauss_(&Rad2difXs_er,&El,&Eh,&eps);}else{ Rad2Xs = dgauss_(&Rad2difXs_er2,&El,&Eh,&eps);}
2424 if(!Rad2Xs) return Rad2Xs;;
2425 }
2426 }
2427 return Rad2Xs; // the leading second order correction
2428}
◆ get_mode()
| std::vector< EvtId > EvtConExc::get_mode | ( | int | mode | ) |
Definition at line 1246 of file EvtConExc.cc.
1246 {
1247 int _ndaugs;
1248 EvtId daugs[100];
1249 if(mode==0){
1250 _ndaugs =2;
1253 }else if(mode==1){
1254 _ndaugs =2;
1257 }else if(mode==2){
1258 _ndaugs =2;
1261 }else if(mode==3){
1262 _ndaugs =2;
1265 }else if(mode==4){
1266 _ndaugs =2;
1269 }else if(mode==5){
1270 _ndaugs =2;
1273 }else if(mode ==6){
1274 _ndaugs =2;
1277 }else if(mode==7){
1278 _ndaugs =3;
1282 }else if(mode==8){
1283 _ndaugs =3;
1287 }else if(mode==9){
1288 _ndaugs =3;
1292 }else if(mode==10){
1293 _ndaugs =3;
1297 }else if(mode==11){
1298 _ndaugs =3;
1302 }else if(mode==12){
1303 _ndaugs =4;
1308 }else if(mode==13){
1309 _ndaugs =4;
1314 }else if(mode==14){
1315 _ndaugs =4;
1320 }else if(mode==15){
1321 _ndaugs =4;
1326 }else if(mode==16){
1327 _ndaugs =4;
1332 }else if(mode==17){
1333 _ndaugs =5;
1339 }else if(mode==18){
1340 _ndaugs =5;
1346 }else if(mode==19){
1347 _ndaugs =5;
1353 }else if(mode==20){
1354 _ndaugs =5;
1360 }else if(mode==21){
1361 _ndaugs =6;
1368 }else if(mode==22){
1369 _ndaugs =6;
1376 }else if(mode==23){
1377 _ndaugs =2;
1380 }else if(mode==24){
1381 _ndaugs =2;
1384 }else if(mode==25){
1385 _ndaugs =2;
1388 }else if(mode==26){
1389 _ndaugs =2;
1392 }else if(mode==27){
1393 _ndaugs =2;
1396 }else if(mode==28){
1397 _ndaugs =3;
1401 }else if(mode==29){
1402 _ndaugs =3;
1406 }else if(mode==30){
1407 _ndaugs =3;
1411 }else if(mode==31){
1412 _ndaugs =3;
1416 }else if(mode==32){
1417 _ndaugs =3;
1421 }else if(mode==33){
1422 _ndaugs =3;
1426 }else if(mode==34){
1427 _ndaugs =3;
1431 }else if(mode==35){
1432 _ndaugs =3;
1436 }else if(mode==36){
1437 _ndaugs =2;
1440 }else if(mode==37){
1441 _ndaugs =3;
1445 }else if(mode==38){
1446 _ndaugs =3;
1450 }else if(mode==39){
1451 _ndaugs =2;
1454 }else if(mode==40){
1455 _ndaugs =3;
1459 }else if(mode==41){
1460 _ndaugs =3;
1464 }else if(mode==42){
1465 _ndaugs =3;
1469 }else if(mode==43){
1470 _ndaugs =4;
1475 }else if(mode==44){
1476 _ndaugs =2;
1479 }else if(mode==45){
1480 _ndaugs =2;
1483 }else if(mode==46){
1484 _ndaugs =2;
1487 }else if(mode==47){
1488 _ndaugs =2;
1491 }else if(mode==48){
1492 _ndaugs =2;
1495 }else if(mode==49){
1496 _ndaugs =3;
1500 }else if(mode==50){
1501 _ndaugs =3;
1505 }else if(mode==51){
1506 _ndaugs =3;
1510 }else if(mode==52){
1511 _ndaugs =2;
1514 }else if(mode==59){
1515 _ndaugs =3;
1519 }else if(mode==60){
1520 _ndaugs =3;
1524 }else if(mode==61){
1525 _ndaugs =3;
1529 }else if(mode==62){
1530 _ndaugs =3;
1534 }else if(mode==63){
1535 _ndaugs =3;
1539 }else if(mode==64){
1540 _ndaugs =3;
1544 }else if(mode==65){
1545 _ndaugs =3;
1549 }else if(mode==66){
1550 _ndaugs =3;
1554 }else if(mode==67){
1555 _ndaugs =2;
1558 }else if(mode==68){
1559 _ndaugs =2;
1562 }else if(mode==69){
1563 _ndaugs =2;
1566 }else if(mode==70){
1567 _ndaugs =2;
1570 }else if(mode==71){
1571 _ndaugs =2;
1574 }else if(mode==72){
1575 _ndaugs =2;
1578 }else if(mode==73){
1579 _ndaugs =2;
1582 }else if(mode==74){
1583 _ndaugs =2;
1586 }else if(mode==75){
1587 _ndaugs =3;
1591 }else if(mode==76){
1592 _ndaugs =3;
1596 }else if(mode==77){
1597 _ndaugs =3;
1601 }else if(mode==78){
1602 _ndaugs =3;
1606 }else if(mode==79){
1607 _ndaugs =3;
1611 }else if(mode==80){
1612 _ndaugs =2;
1615 }else if(mode==81){
1616 _ndaugs =3;
1620 }else if(mode==82){
1621 _ndaugs =3;
1625 }else if(mode==83){
1626 _ndaugs =3;
1630 }else if(mode==84){
1631 _ndaugs =3;
1635 }else if(mode==85){
1636 _ndaugs =2;
1639 }else if(mode==90){
1640 _ndaugs =3;
1644 }else if(mode==91){
1645 _ndaugs =3;
1649 }else if(mode==92){
1650 _ndaugs =3;
1654 }else if(mode==93){
1655 _ndaugs =2;
1658 }else if(mode==94){
1659 _ndaugs =2;
1662 }else if(mode==95){
1663 _ndaugs =2;
1666 }else if(mode==96){
1667 _ndaugs =2;
1670 }else if(mode==111){
1671 _ndaugs =2;
1674 }else if(mode==112){
1675 _ndaugs =2;
1678 }else if(mode==74100){
1679 _ndaugs =2;
1682 }else if(mode==74101){
1683 _ndaugs =2;
1686 }else if(mode==74102){
1687 _ndaugs =2;
1690 }else if(mode==74103){
1691 _ndaugs =2;
1694 }else if(mode==74104){
1695 _ndaugs =2;
1698 }else if(mode==74105){
1699 _ndaugs =2;
1702 }else if(mode==74106){
1703 _ndaugs =2;
1706 }else if(mode==74107){
1707 _ndaugs =2;
1710 }else if(mode==74110){
1711 _ndaugs =2;
1715 EvtPDL::reSetMass(myvpho,mhds*0.9); //for calculating maxXS, it will be resent when this mode is selected
1716 _modeFlag.clear();
1717 _modeFlag.push_back(74110); //R-value generator tag
1718 _modeFlag.push_back(0); //to contain the mode selected
1719 }else if(mode==-1){
1720 _ndaugs = nson;
1721 for(int i=0;i<nson;i++){ daugs[i]=son[i]; }
1723 }else if(mode==-2){
1724 _ndaugs = nson;
1725 for(int i=0;i<nson;i++){ daugs[i]=son[i]; }
1726 }else if(mode==-100){
1727 _ndaugs =1;
1729 _modeFlag.clear();
1730 _modeFlag.push_back(-100); //R-value generator tag
1731 _modeFlag.push_back(0); //to contain the mode selected
1732
1733 }else if(mode==10000){//use for check ISR
1734 _ndaugs =2;
1737 }else {
1738 std::cout<<"Bad mode_index number " <<mode<<", please refer to the manual."<<std::endl;
1739 ::abort();
1740 }
1741 // std::cout<<_ndaugs<<", "<<daugs[0]<<daugs[1]<<std::endl;
1742
1743 if(VISR){
1744 _ndaugs=2;
1747 }
1748
1749 std::vector<EvtId> theDaugs;
1750 for(int i=0;i<_ndaugs;i++){
1751 theDaugs.push_back(daugs[i]);
1752 }
1753 if(theDaugs.size()) {return theDaugs;} else {std::cout<<"EvtConExc::get_mode: zero size"<<std::endl;abort();}
1754}
Referenced by selectMode().
◆ get_mode_index()
| int EvtConExc::get_mode_index | ( | int | mode | ) |
Definition at line 3440 of file EvtConExc.cc.
3440 {
3441 int i=0;
3442 for(int j=0;i<vmode.size();j++){
3443 if(mode==vmode[j]) return j;
3444 }
3445 std::cout<<" EvtConExc::get_mode_index: no index is found in vmode"<<std::endl;
3446 abort();
3447}
Referenced by getObsXsection().
◆ getDaugId()
|
inline |
◆ getModeIndex()
| int EvtConExc::getModeIndex | ( | int | m | ) |
Definition at line 3581 of file EvtConExc.cc.
3581 {
3582 for (int i=0;i<vmode.size();i++){
3583 if(m==vmode[i]) return i;
3584 }
3585 std::cout<<"EvtConExc::getModeIndex: no index in vmode found "<<std::endl;
3586 abort();
3587}
Referenced by checkEvtRatio(), and writeDecayTabel().
◆ getName()
|
virtual |
Implements EvtDecayBase.
Definition at line 241 of file EvtConExc.cc.
241 {
242
243 model_name="ConExc";
244
245}
◆ getNdaugs()
|
inline |
◆ getObsXsection()
| double EvtConExc::getObsXsection | ( | double | mhds, |
| int | mode | ||
| ) |
Definition at line 3449 of file EvtConExc.cc.
3449 {
3450 double hwid=(AA[0]-AA[1])/2.;
3454 for(int i=0;i<600;i++){
3455 if( (AA[i]-hwid)<=Egam && Egam<(AA[i]+hwid) ) return VXS[idx][i];
3456 }
3457
3459 abort();
3460}
◆ getResMass()
| void EvtConExc::getResMass | ( | ) |
Definition at line 3118 of file EvtConExc.cc.
3118 {
3120 mjsi = EvtPDL::getMeanMass(myres);
3121 wjsi = EvtPDL::getWidth(myres);
3122
3124 mpsip= EvtPDL::getMeanMass(myres);
3125 wpsip= EvtPDL::getWidth(myres);
3126
3128 mpsipp= EvtPDL::getMeanMass(myres);
3129 wpsipp= EvtPDL::getWidth(myres);
3130
3132 mphi = EvtPDL::getMeanMass(myres);
3133 wphi = EvtPDL::getWidth(myres);
3134
3136 momega= EvtPDL::getMeanMass(myres);
3137 womega= EvtPDL::getWidth(myres);
3138
3140 mrho0 = EvtPDL::getMeanMass(myres);
3141 wrho0 = EvtPDL::getWidth(myres);
3142
3144 mrho3s= EvtPDL::getMeanMass(myres);
3145 wrho3s= EvtPDL::getWidth(myres);
3146
3147
3149 momega2s= EvtPDL::getMeanMass(myres);
3150 womega2s= EvtPDL::getWidth(myres);
3151
3152}
Referenced by init().
◆ getSelectedMode()
|
inline |
Definition at line 210 of file EvtConExc.hh.
210{return _selectedMode;}
◆ getVP()
| double EvtConExc::getVP | ( | double | cms | ) |
Definition at line 3360 of file EvtConExc.cc.
3360 {
3361 double xx,r1,i1;
3362 double x1,y1,y2;
3363 xx=0;
3364 int mytag=1;
3365 for(int i=0;i<4001;i++){
3366 x1=vpx[i];
3367 y1=vpr[i];
3368 y2=vpi[i];
3369 if(xx<=mx && mx <x1){xx=x1; r1=y1; i1=y2; mytag=0; break;}
3370 xx=x1; r1=y1; i1=y2;
3371 }
3372 if(mytag==1){std::cout<<"No vacuum polarization value found"<<std::endl;abort();}
3373 EvtComplex cvp(r1,i1);
3375 if(3.0933<mx && mx<3.1035) return 1.; //J/psi peak excluded
3376 if(3.6810<mx && mx<3.6913) return 1.; //psi(2S) peak excluded
3377 return thevp*thevp;
3378}
Definition: EvtComplex.hh:28
Referenced by init(), Rad2(), and SoftPhoton_xs().
◆ hadron_angle_sampling()
| bool EvtConExc::hadron_angle_sampling | ( | EvtVector4R | ppi, |
| EvtVector4R | pcm | ||
| ) |
Definition at line 2258 of file EvtConExc.cc.
2258 {
2259 EvtVector4R pbst=-1*pcm;
2262 if((_mode>=0 && _mode<=5) || _mode==44 || _mode==96){//ee->two baryon mode, VP, SP, mode=-2 is excluded
2264 return byn_ang;
2265 }else if(_mode==6 || _mode==45 || _mode==46 || _mode==70 || _mode==71|| _mode==93){// ee->PP (pi+pi-,..) mode
2267 return msn_ang;
2268 }else if(_mode==23 || _mode==24 || _mode==25 || _mode==26 || _mode==27 || _mode==36 || _mode==39 || _mode==47 || _mode==48 || _mode==52 || _mode==68 || _mode==69 || _mode==72 || _mode==73 || _mode==80 || _mode==94 || _mode==95){
2270 return msn_ang;
2271 }else if(_mode==-2){
2276 return msn_ang;
2277 }else if(type0==EvtSpinType::VECTOR &&type1==EvtSpinType::SCALAR || type1==EvtSpinType::VECTOR &&type0==EvtSpinType::SCALAR){
2279 return msn_ang;
2280 }
2281 }
2282 return true;
2283}
EvtDiracSpinor boostTo(const EvtDiracSpinor &sp, const EvtVector4R p4)
Definition: EvtDiracSpinor.cc:75
bool VP_sampling(EvtVector4R pcm, EvtVector4R pi)
Definition: EvtConExc.cc:2557
bool meson_sampling(EvtVector4R pcm, EvtVector4R pi)
Definition: EvtConExc.cc:2546
bool baryon_sampling(EvtVector4R pcm, EvtVector4R pi)
Definition: EvtConExc.cc:2529
Referenced by decay().
◆ init()
|
virtual |
Reimplemented from EvtDecayBase.
Definition at line 254 of file EvtConExc.cc.
254 {
255 mydbg=false;
256 myFisr.clear();
257 ReadVP();
258 VXS.resize(120);
259 for(int i=0;i<120;i++){
260 VXS[i].resize(600);
261 }
262 _mode = getArg(0);
263 for(int i=0;i<97;i++){
264 if(53<=i && i<=58) continue;
265 if(86<=i && i<=89) continue;
266 if(_mode==74110 ||_mode==-100) vmode.push_back(i);
267 }
268 if(_mode==74110||_mode==-100){
269 vmode.push_back(74100);
270 vmode.push_back(74101);
271 vmode.push_back(74102);
272 vmode.push_back(74103);
273 vmode.push_back(74104);
274 vmode.push_back(74105);
275 vmode.push_back(74110);
276 }
277
278 std::cout<<"==== Parameters set by users====="<<std::endl;
279 for(int i=0;i<ncommand;i++){
280 std::cout<<commands[i]<<std::endl;
281 }
282 std::cout<<"==================================="<<std::endl;
283 InitPars();
284 std::cout<<threshold<<" "<<beamEnergySpread<<std::endl;
285 //----------------
286 // check that there are 0 arguments
287 // checkNArg(1);
288 //Vector ISR process
289 VISR=0;
291 if(getDaugs()[0]==EvtPDL::getId("gamma")&&getDaugs()[1]==EvtPDL::getId("gamma*") || getDaugs()[0]==EvtPDL::getId("gamma*")&&getDaugs()[1]==EvtPDL::getId("gamma")) VISR=1;
293 }
294
295 // cmsspread=0.0015; //CMC energy spread
296 testflag=0;
297 getResMass();
299 radflag=0;
300 mydbg=false;
301 _mode = getArg(0);
302 pdgcode = getArg(1);
303 pid = EvtPDL::evtIdFromStdHep(pdgcode );
304 nson = getNArg()-2;
305 std::cout<<"The decay daughters are "<<std::endl;
306 for(int i=2;i<getNArg();i++){son[i-2]=EvtPDL::evtIdFromStdHep(getArg(i));std::cout<<EvtPDL::name(son[i-2])<<" ";}
307 std::cout<<std::endl;
309 radflag=0;mydbg=false;
310 _mode = getArg(0);
311 nson = getNArg()-1;
314 _mode = getArg(0);
315 multimode=1;
316 vmd.clear();
318 std::cout<<" multi-exclusive mode "<<std::endl;
320 std::cout<<std::endl;
321 if(vmd[vmd.size()-1]==74110){std::cout<<"74110 is not allowd for multimode simulation"<<std::endl;abort();}
325 }else{std::cout<<"ConExc:number of parameter should be 1,2 or 3, but you set "<<getNArg()<<std::endl;::abort(); }
326 //--- debugging
327 //std::cout<<_mode<<" "<<pdgcode<<" "<<nson<<" "<<EvtPDL::name(son[0])<<" "<<EvtPDL::name(son[1])<<std::endl;
328
330 init_mode(_mode);
331 XS_max=-1;
332 //-----debugging to make root file
333 mydbg = false; //designer need to remove this line
334 if(mydbg){
335 myfile = new TFile("mydbg.root","RECREATE");
336 xs=new TTree ("xs","momentum of rad. gamma and hadrons");
337 xs->Branch("imode",&imode,"imode/I");
338 xs->Branch("pgam",pgam,"pgam[4]/D");
339 xs->Branch("phds",phds,"phds[4]/D");
340 xs->Branch("ph1",ph1,"ph1[4]/D");
341 xs->Branch("ph2",ph2,"ph2[4]/D");
342 xs->Branch("mhds",&mhds,"mhds/D");
343 xs->Branch("mass1",&mass1,"mass1/D");
344 xs->Branch("mass2",&mass2,"mass2/D");
345 xs->Branch("costheta",&costheta,"costheta/D");
346 xs->Branch("cosp",&cosp,"cosp/D");
347 xs->Branch("cosk",&cosk,"cosk/D");
348 xs->Branch("sumxs",&sumxs,"sumxs/D");
349 xs->Branch("selectmode",&selectmode,"selectmode/D");
350 }
351 //--- prepare the output information
353 EvtPDL::reSetWidth(parentId,0.0);
355 std::cout<<"parent mass = "<<parentMass<<std::endl;
356
357
358 EvtVector4R p_init(parentMass,0.0,0.0,0.0);
360 theparent = p;
363 if(_mode==-100) {
368 }
369 double mth=0;
371 if(_mode ==-1){
375 if(mx<mth){std::cout<<"The CMS energy is lower than the threshold of the final states"<<std::endl;abort();}
376 }else if(_mode == -2){
379 _cms = mx;
381
382 std::cout<<"The specified mode is "<<_mode<<std::endl;
383 EvtConExc::getMode = _mode;
384 //std::cout<<"getXlw= "<<mth<<std::endl;
385
386 Egamcut = 0.001; //set photon energy cut according to the BESIII detector
387 double Esig = 0.0001; //sigularity energy
391
392 //for vaccum polarization
393 float fe,fst2,fder,ferrder,fdeg,ferrdeg;
394 fe=_cms;
395 //using the vacuum pol. value as given by http://www-com.physik.hu-berlin.de/~fjeger/software.html
397 if(3.0943<_cms && _cms<3.102) vph=1;// For J/psi, the vacuum factor is included in the resonance
398 std::cout<<"sum of leptonic, hadronic contributions(u,d,s,c,b) to vacuum polarization: "<<vph<<" @"<<fe<<"GeV"<<std::endl;
399
400 //caculate the xs for ISR to psiprime, J/psi and phi narrow resonance.
401 double totxsIRS=0;
402 init_Br_ee();
404 for(int jj=1;jj<_ndaugs;jj++){
405 mthrld += EvtPDL::getMeanMass(daugs[jj]);
406 }
407
408 ISRXS.clear();ISRM.clear();ISRFLAG.clear();
409 for(int ii=0;ii<3;ii++){
411 if(mx<mR || _mode !=74110) continue;
413 std::cout<<"The corss section for gamma "<<EvtPDL::name(ResId[ii]).c_str()<<" is: "<< narRxs<<"*Br (nb)."<<std::endl;
414 ISRXS.push_back(narRxs);
415 ISRM.push_back(mR);
416 ISRFLAG.push_back(1.);
417 ISRID.push_back(ResId[ii]);
418 totxsIRS += narRxs;
419 }
420 std::cout<<"==========================================================================="<<std::endl;
421
422 //-- calculated with MC integration method
427 int nmc=1000000;
429 _er0 = _er1; // stored when calculate _xs0
430 std::cout<<"_er0= "<<_er0<<std::endl;
432 std::cout<<"_xs1= "<<_xs0<<std::endl;
433 double xs1_err = _er1;
434
435
436 //--- sigularity point x from 0 to 0.0001GeV
439 _xs0 += sig_xs;
440 std::cout<<"_xs0= "<<_xs0<<std::endl;
441
442 //prepare the observed cross section with interpotation function divdif in CERNLIB
443 testflag=1;
444 int Nstp=598;
445 double stp=(EgamH - Egamcut)/Nstp;
446 for(int i=0;i<Nstp;i++){//points within Egam(max) ~ Egamcut=25MeV
447 double Eg0=EgamH - i*stp;
448 double Eg1=EgamH - (i+1)*stp;
449 int nmc=20000;
450 int nint=100;
451 //double xsi= gamHXSection(s,Eg1,Eg0,nmc);
453 AA[i]=(Eg1+Eg0)/2;
456 double binwidth = mh2-mhi;
457 //if(_mode==74110) xsi += addNarrowRXS(mhi,binwidth);
458 if(i==0) AF[0]=xsi;
459 if(i>=1) AF[i]=AF[i-1]+xsi;
460 RadXS[i]=xsi/stp;
461 }
462 //add the interval 0~Esig
463 AA[598]=Egamcut; AA[599]=0; //AA[i]: E_gamma
464 AF[598]=AF[597];
465 AF[599]=AF[598]+ _xs0;
466 RadXS[599]=_xs0;
467
468
469 //prepare observed cross section for each mode
471 for(int i=0;i<vmode.size();i++){
473 }
475 //after mk_VXS, restore the myxsection to _mode
476 if(_mode !=-100){
479 }
480 //debugging VXS
481 /*
482 double xtmp=0;
483 for(int i=0;i<vmode.size();i++){
484 xtmp+=VXS[i][599];
485 for(int j=0;j<600;j++){
486 std::cout<<VXS[i][j]<<" ";
487 }
488 std::cout<<std::endl;
489 }
490 */
491 //output observed cross section for the gamma Resonance mode
492 //std::cout<<"vmode.size======="<<vmode.size()<<std::endl;
493 for(int i=0;i<vmode.size();i++){
494 int md=vmode[i];
495 if(md<74100 || md>74106) continue;
496 std::string partname="";
497 if(md==74100) {partname="J/psi";}
498 else if(md==74101) {partname="psi(2S)";}
499 else if(md==74102) {partname="psi(3770)";}
500 else if(md==74103) {partname="phi";}
501 else if(md==74104) {partname="omega";}
502 else if(md==74105) {partname="rho0";}
503 else if(md==74106) {partname="rho(3S)0";}
504 else{;}
505 std::cout<<"The observed cross section for gamma "<<partname<<": "<<VXS[i][599]<<" nb"<<std::endl;
506 }
507 //--
508
509 for(int i=0;i<600;i++){
511 //std::cout<<i<<" Egamma "<<AA[i]<<" Mhadons "<<MH[i]<<std::endl;
512 }
513 //for debugging
514 //for(int i=0;i<600;i++){double s=_cms*_cms; double mhds=sqrt(s-2*_cms*AA[i]);std::cout<<"Mhds="<<mhds<<" Egam="<<AA[i]<<" "<<AF[i]<<std::endl;}
515 //std::cout<<"VXS[79][599]: "<<VXS[79][599]<<" VXS[79][598]= "<<VXS[79][598]<<std::endl;
516 std::cout<<"EgamH = "<<EgamH<<", obsvXS = "<<_xs0+_xs1<<", _xs1 = "<<_xs1 << ", _xs0 = "<<_xs0<<std::endl;
517 std::cout<<"EgamH = "<<EgamH<<", AF[599] = "<< AF[599]<< ", AF[598] = "<<AF[598]<<", _xs0 = "<<_xs0<<std::endl;
518 //double Etst=0.0241838;
519 //std::cout<<"Etst="<<Etst<<" "<<gamHXSection(s,Etst,EgamH,nmc)<<" "<< lgr(AA,AF,600,Etst)<<std::endl; abort();
520
521 //for information output
524 double obsvXS = _xs0 + _xs1; //gamHXSection(mth ,mx);
525 double obsvXS_er= _er0 + xs1_err;
526 double corr = obsvXS/bornXS;
527 double corr_er =corr*sqrt(bornXS_er*bornXS_er/bornXS/bornXS + obsvXS_er*obsvXS_er/obsvXS/obsvXS);
528
529
530 if(bornXS==0){std::cout<<"EvtConExc::init : The Born cross section at this point is zero!"<<std::endl;abort();}
531 if(obsvXS==0){std::cout<<"EvtConExc::init : The observed cross section at this point is zero!"<<std::endl;abort();}
532
533 //_xs0 += bornXS;// events with very soft photon (Egam<0.025) are taken as the born process
534 //_er0 = sqrt(_er0*_er0 + bornXS_er*bornXS_er);
535
536 std::cout<<""<<std::endl;
537 std::cout<<"========= Generation using cross section (Egamcut="<<Egamcut<<" GeV) =============="<<std::endl;
538 std::cout<<" sqrt(s)= "<<mx<< " GeV "<<std::endl;
539 if(_mode>=0 || _mode ==-2 || _mode ==-100){
540 std::cout<<"The generated born cross section (sigma_b): "<<_xs0<<"+/-"<<_er0<<" in Unit "<<_unit.c_str()<<std::endl;
541 std::cout<<"The generated radiative correction cross section (sigma_isr): "<<_xs1<<"+/-"<<xs1_err<<" in Unit "<<_unit.c_str()<<std::endl;
542 std::cout<<"The Born process is sampled according to the ratio sigma_b/(sigma_b+sigma_isr)"<<std::endl;
543 std::cout<<"The radiative correction factor f_ISR*f_vacuum= sigma_obs/sigma_born(s) = "<<corr<<"+/-"<< fabs(corr_er)<<","<<std::endl;
544 std::cout<<"which is calcualted with these quantities:"<<std::endl;
545 std::cout<<"the observed cross section is "<<obsvXS<<"+/-"<< obsvXS_er<<_unit.c_str()<<std::endl;
546 std::cout<<"and the Born cross section (s) is "<<bornXS<<"+/-"<<bornXS_er<<_unit.c_str()<<std::endl;
547 std::cout<<"and the vaccum polarziation factor (lep. + hadr.) 1/|1-Pi|^2= "<<vph<<std::endl;
548 std::cout<<"Within the range "<<myxsection->getXlw()<<" ~ "<<myxsection->getXup()<<" GeV, "<<myxsection->getMsg().c_str()<<std::endl;
549 std::cout<<"==========================================================================="<<std::endl;
550 }else if(_mode==-1){
551 std::cout<<"The generated born cross section (sigma_b): "<<_xs0<<" *Br_ee"<<" in Unit "<<_unit.c_str()<<std::endl;
552 std::cout<<"The generated radiative correction cross section (sigma_isr): "<<_xs1<<"*Br_ee in Unit "<<_unit.c_str()<<std::endl;
553 std::cout<<"The Born process is sampled according to the ratio sigma_b/(sigma_b+sigma_isr)"<<std::endl;
554 std::cout<<"The radiative correction factor sigma_obs/sigma_born(s) = (1+delta) = "<<corr<<"+/-"<< fabs(corr_er)<<std::endl;
555 std::cout<<"==========================================================================="<<std::endl;
556 }
557 std::cout<<std::endl;
558 std::cout<<std::endl;
559
560 findMaxXS(p);
562 //--debugging
563 //std::cout<<"maxXS= "<<maxXS<<std::endl;
564
565 if(_xs0 == 0 && _xs1==0){std::cout<<"EvtConExc::ini() has zero cross section"<<std::endl;abort();}
566
567 std::cout<<std::endl;
568 std::cout<<std::endl;
571 //--- for debugging
572 if(mydbg && _mode!=74110){
574 double xx[10000],yy[10000],xer[10000],yer[10000];
575 for(int i=0;i<nd;i++){
579 xer[i]=0.;
580 //std::cout<<"yy "<<yy[i]<<std::endl;
581 }
582 myth=new TH1F("myth","Exp.data",200,xx[0],xx[nd]);
583 for(int i=0;i<nd;i++){
584 myth->Fill(xx[i],yy[i]);
585 }
586 myth->SetError(yer);
587 myth->Write();
588
589 //std::cout<<"mass1= "<<mass1<<" mass2= "<<mass2<<std::endl;
590 }else if(mydbg && _mode==74110){
592 double xx[10000],yy[10000],xer[10000],yer[10000],ysum[10000],yysum[10000];
593 for(int i=0;i<nd;i++){
597 xer[i]=0.;
598 std::cout<<"yy "<<yy[i]<<std::endl;
599 }
600
601#include "sty.C"
602 double xhigh=5.0;
604 TGraphErrors *Gdt = new TGraphErrors(nd,xx,yy,xer,yer);
605
609 for(int i=0;i<600;i++){
612 if(xsi==0) continue;
613 myth->Fill(mx,xsi);
614 //std::cout<<mx<<" "<<xsi<<std::endl;
615 }
616
617 for(int i=0;i<600;i++){
619 ysum[i]=sumExc(mx);
620 if(ysum[i]==0) continue;
621 Xsum->Fill(mx,ysum[i]);
622 //std::cout<<mx<<" "<<ysum[i]<<std::endl;
623 }
624
625 for(int i=0;i<nd;i++){
626 yysum[i]=sumExc(xx[i]);
627 }
628
629 myth->SetError(yer);
630 myth->Write();
631 Xsum->Write();
632
633 TGraphErrors *Gsum = new TGraphErrors(nd,xx,yysum,xer,yer);
634 //draw graph
635 double ex[600]={600*0};
636 double ey[600],ta[600];
637 double exz[600]={600*0};
638 double eyz[600]={600*0};
639 for(int i=0;i<600;i++){
640 ey[i]=AF[i]*0.001;
641 }
642 TGraphErrors *gr0 = new TGraphErrors(600,MH,AF,ex,ey);
643 TGraphErrors *gr1 = new TGraphErrors(600,MH,RadXS,exz,eyz);
644 TPostScript *ps = new TPostScript("xsobs.ps",111);
645 TCanvas *c1 = new TCanvas("c1","TGraphs for data",200,10,600,400);
646 gPad-> SetLogy(1);
647 // c1->SetLogy(1);
648 gStyle->SetCanvasColor(0);
649 gStyle->SetStatBorderSize(1);
650 c1->Divide(1,1);
651
652 c1->Update();
653 ps->NewPage();
654 c1->Draw();
655 c1->cd(1);
656 c1->SetLogy(1);
657 gr0->SetMarkerStyle(10);
658 gr0->Draw("AP");
659 gr0->GetXaxis()->SetTitle("M_{hadrons} (GeV)");
660 gr0->GetYaxis()->SetTitle("observed cross section (nb)");
661 gr0->Write();
662
663 c1->Update();
664 ps->NewPage();
665 c1->Draw();
666 c1->cd(1);
667 c1->SetLogy(1);
668 gr1->SetMarkerStyle(10);
669 gr1->Draw("AP");
670 gr1->GetXaxis()->SetTitle("M_{hadrons} (GeV)");
671 gr1->GetYaxis()->SetTitle("Rad*BornXS");
672 gr1->Write();
673
674 //--------- imposing graphs to a pad
675 TMultiGraph *mg = new TMultiGraph();
676 mg->SetTitle("Exclusion graphs");
677 Gdt->SetMarkerStyle(2);
678 Gdt->SetMarkerSize(1);
679 Gsum->SetLineColor(2);
680 Gsum->SetLineWidth(1);
681 mg->Add(Gdt);
682 mg->Add(Gsum);
683
684 c1->Update();
685 ps->NewPage();
686 c1->Draw();
687 c1->cd(1);
688 gPad-> SetLogy(1);
689 mg->Draw("APL");
690 mg->GetXaxis()->SetTitle("M_{hadrons} (GeV)");
691 mg->GetYaxis()->SetTitle("observed cross section (nb)");
692 mg->Write();
693 //-------
694
695 c1->Update();
696 ps->NewPage();
697 c1->Draw();
698 c1->cd(1);
699 Gdt->SetMarkerStyle(2);
700 Gdt->Draw("AP");
701 Gdt->GetXaxis()->SetTitle("M_{hadrons} (GeV)");
702 Gdt->GetYaxis()->SetTitle("observed cross section (nb)");
703 Gdt->Write();
704
705 c1->Update();
706 ps->NewPage();
707 c1->Draw();
708 c1->cd(1);
709 Gsum->SetMarkerStyle(2);
710 Gsum->SetMarkerSize(1);
711 Gsum->Draw("AP");
712 Gsum->SetLineWidth(0);
713 Gsum->GetXaxis()->SetTitle("M_{hadrons} (GeV)");
714 Gsum->GetYaxis()->SetTitle("observed cross section (nb)");
715 Gsum->Write();
716
717 c1->Update();
718 ps->NewPage();
719 c1->Draw();
720 c1->cd(1);
721 // gPad-> SetLogx(1);
722 Gdt->SetMarkerStyle(2);
723 Gdt->SetMarkerSize(1);
724 Gdt->SetMaximum(8000.0);
725 Gsum->SetLineColor(2);
726 Gsum->SetLineWidth(1.5);
727 Gsum->Draw("AL"); //A: draw axis
728 Gdt->Draw("P"); // superpose to the Gsum, without A-option
729 Gsum->GetXaxis()->SetTitle("M_{hadrons} (GeV)");
730 Gsum->GetYaxis()->SetTitle("cross section (nb)");
731 Gsum->Write();
732
733 ps->Close();
734 } //if(mydbg)_block
735 //-------------------------
736
737}
double gamHXSection(EvtParticle *p, double El, double Eh, int nmc=100000)
Definition: EvtConExc.cc:2286
void mk_VXS(double Esig, double Egamcut, double EgamH, int midx)
Definition: EvtConExc.cc:3401
◆ init_Br_ee()
| void EvtConExc::init_Br_ee | ( | ) |
Definition at line 2686 of file EvtConExc.cc.
2686 {
2687 // double psipp_ee=9.6E-06;
2688 double psip_ee =7.73E-03;
2689 double jsi_ee =5.94*0.01;
2690 double phi_ee =2.954E-04;
2691 // double omega_ee=7.28E-05;
2692 // double rho_ee = 4.72E-05;
2699 BR_ee.clear();
2700 ResId.clear();
2701
2702 //BR_ee.push_back(rho_ee);
2703 //BR_ee.push_back(omega_ee);
2704 BR_ee.push_back(phi_ee);
2705 BR_ee.push_back(jsi_ee);
2706 BR_ee.push_back(psip_ee);
2707 //BR_ee.push_back(psipp_ee);
2708
2709 //ResId.push_back(rhoId);
2710 //ResId.push_back(omegaId);
2711 ResId.push_back(phiId);
2712 ResId.push_back(jsiId);
2713 ResId.push_back(pspId);
2714 //ResId.push_back(psppId);
2715
2716}
Referenced by init().
◆ init_mode()
| void EvtConExc::init_mode | ( | int | mode | ) |
Definition at line 741 of file EvtConExc.cc.
741 {
742 if(mode==0){
743 _ndaugs =2;
746 }else if(mode==1){
747 _ndaugs =2;
750 }else if(mode==2){
751 _ndaugs =2;
754 }else if(mode==3){
755 _ndaugs =2;
758 }else if(mode==4){
759 _ndaugs =2;
762 }else if(mode==5){
763 _ndaugs =2;
766 }else if(mode==6){
767 _ndaugs =2;
770 }else if(mode==7){
771 _ndaugs =3;
775 }else if(mode==8){
776 _ndaugs =3;
780 }else if(mode==9){
781 _ndaugs =3;
785 }else if(mode==10){
786 _ndaugs =3;
790 }else if(mode==11){
791 _ndaugs =3;
795 }else if(mode==12){
796 _ndaugs =4;
801 }else if(mode==13){
802 _ndaugs =4;
807 }else if(mode==14){
808 _ndaugs =4;
813 }else if(mode==15){
814 _ndaugs =4;
819 }else if(mode==16){
820 _ndaugs =4;
825 }else if(mode==17){
826 _ndaugs =5;
832 }else if(mode==18){
833 _ndaugs =5;
839 }else if(mode==19){
840 _ndaugs =5;
846 }else if(mode==20){
847 _ndaugs =5;
853 }else if(mode==21){
854 _ndaugs =6;
861 }else if(mode==22){
862 _ndaugs =6;
869 }else if(mode==23){
870 _ndaugs =2;
873 }else if(mode==24){
874 _ndaugs =2;
877 }else if(mode==25){
878 _ndaugs =2;
881 }else if(mode==26){
882 _ndaugs =2;
885 }else if(mode==27){
886 _ndaugs =2;
889 }else if(mode==28){
890 _ndaugs =3;
894 }else if(mode==29){
895 _ndaugs =3;
899 }else if(mode==30){
900 _ndaugs =3;
904 }else if(mode==31){
905 _ndaugs =3;
909 }else if(mode==32){
910 _ndaugs =3;
914 }else if(mode==33){
915 _ndaugs =3;
919 }else if(mode==34){
920 _ndaugs =3;
924 }else if(mode==35){
925 _ndaugs =3;
929 }else if(mode==36){
930 _ndaugs =2;
933 }else if(mode==37){
934 _ndaugs =3;
938 }else if(mode==38){
939 _ndaugs =3;
943 }else if(mode==39){
944 _ndaugs =2;
947 }else if(mode==40){
948 _ndaugs =3;
952 }else if(mode==41){
953 _ndaugs =3;
957 }else if(mode==42){
958 _ndaugs =3;
962 }else if(mode==43){
963 _ndaugs =4;
968 }else if(mode==44){
969 _ndaugs =2;
972 }else if(mode==45){
973 _ndaugs =2;
976 }else if(mode==46){
977 _ndaugs =2;
980 }else if(mode==47){
981 _ndaugs =2;
984 }else if(mode==48){
985 _ndaugs =2;
988 }else if(mode==49){
989 _ndaugs =3;
993 }else if(mode==50){
994 _ndaugs =3;
998 }else if(mode==51){
999 _ndaugs =3;
1003 }else if(mode==52){
1004 _ndaugs =2;
1007 }else if(mode==59){
1008 _ndaugs =3;
1012 }else if(mode==60){
1013 _ndaugs =3;
1017 }else if(mode==61){
1018 _ndaugs =3;
1022 }else if(mode==62){
1023 _ndaugs =3;
1027 }else if(mode==63){
1028 _ndaugs =3;
1032 }else if(mode==64){
1033 _ndaugs =3;
1037 }else if(mode==65){
1038 _ndaugs =3;
1042 }else if(mode==66){
1043 _ndaugs =3;
1047 }else if(mode==67){
1048 _ndaugs =2;
1051 }else if(mode==68){
1052 _ndaugs =2;
1055 }else if(mode==69){
1056 _ndaugs =2;
1059 }else if(mode==70){
1060 _ndaugs =2;
1063 }else if(mode==71){
1064 _ndaugs =2;
1067 }else if(mode==72){
1068 _ndaugs =2;
1071 }else if(mode==73){
1072 _ndaugs =2;
1075 }else if(mode==74){
1076 _ndaugs =2;
1079 }else if(mode==75){
1080 _ndaugs =3;
1084 }else if(mode==76){
1085 _ndaugs =3;
1089 }else if(mode==77){
1090 _ndaugs =3;
1094 }else if(mode==78){
1095 _ndaugs =3;
1099 }else if(mode==79){
1100 _ndaugs =3;
1104 }else if(mode==80){
1105 _ndaugs =2;
1108 }else if(mode==81){
1109 _ndaugs =3;
1113 }else if(mode==82){
1114 _ndaugs =3;
1118 }else if(mode==83){
1119 _ndaugs =3;
1123 }else if(mode==84){
1124 _ndaugs =3;
1128 }else if(mode==85){
1129 _ndaugs =2;
1132 }else if(mode==90){
1133 _ndaugs =3;
1137 }else if(mode==91){
1138 _ndaugs =3;
1142 }else if(mode==92){
1143 _ndaugs =3;
1147 }else if(mode==93){
1148 _ndaugs =2;
1151 }else if(mode==94){
1152 _ndaugs =2;
1155 }else if(mode==95){
1156 _ndaugs =2;
1159 }else if(mode==96){
1160 _ndaugs =2;
1163 }else if(mode==111){
1164 _ndaugs =2;
1167 }else if(mode==112){
1168 _ndaugs =2;
1171 }else if(mode==74100){
1172 _ndaugs =1;
1174 }else if(mode==74101){
1175 _ndaugs =1;
1177 }else if(mode==74102){
1178 _ndaugs =1;
1180 }else if(mode==74103){
1181 _ndaugs =1;
1183 }else if(mode==74104){
1184 _ndaugs =1;
1186 }else if(mode==74105){
1187 _ndaugs =1;
1189 }else if(mode==74106){
1190 _ndaugs =1;
1192 }else if(mode==74107){
1193 _ndaugs =1;
1195 }else if(mode==74110){
1196 _ndaugs =1;
1199 EvtPDL::reSetMass(myvpho,mhds*0.9); //for calculating maxXS, it will be resent when this mode is selected
1200 _modeFlag.clear();
1201 _modeFlag.push_back(74110); //R-value generator tag
1202 _modeFlag.push_back(0); //to contain the mode selected
1203 }else if(mode==-1){
1204 _ndaugs = nson;
1205 for(int i=0;i<nson;i++){ daugs[i]=son[i]; }
1207 }else if(mode==-2){
1208 _ndaugs = nson;
1209 for(int i=0;i<nson;i++){ daugs[i]=son[i]; }
1210 }else if(mode==-100){
1211 _ndaugs =1;
1213 _modeFlag.clear();
1214 _modeFlag.push_back(-100); //R-value generator tag
1215 _modeFlag.push_back(0); //to contain the mode selected
1216 return;
1217 }else if(mode==10000){//use for check ISR
1218 _ndaugs =2;
1221 }else {
1222 std::cout<<"Bad mode_index number " <<mode<<", please refer to the manual."<<std::endl;
1223 ::abort();
1224 }
1225 // std::cout<<_ndaugs<<", "<<daugs[0]<<daugs[1]<<std::endl;
1226
1227 if(VISR){
1228 _ndaugs=2;
1231 }
1232
1233 double fmth=0;
1234 for(int i=0;i<_ndaugs;i++){
1236 fmth +=xmi;
1237 }
1240 if(!(_mode==74110 || _mode==-100)){
1241 if(Mthr<fmth) {std::cout<<"Low end of cross section: ("<<Mthr<<" < (mass of final state)"<<fmth<<") GeV."<< std::endl; abort();}
1242 }
1243}
Referenced by decay(), EvtRexc::decay(), init(), and sumExc().
◆ InitPars()
| void EvtConExc::InitPars | ( | ) |
Definition at line 3651 of file EvtConExc.cc.
3651 {
3652 threshold=0;
3653 beamEnergySpread=0;
3654 std::string pattern="=";
3655 for(int i=0;i<ncommand;i++){
3656 std::vector<std::string> result=split(commands[i],pattern);
3659 else{
3660 std::cout<<"Your parameter in decay card \""<<result[0]<<"\" incorect"<<std::endl;
3661 std::cout<<"Possible words: threshold, beamEnergySpread"<<std::endl;
3662 abort();
3663 }
3664 }
3665
3666}
std::vector< std::string > split(std::string str, std::string pattern)
Definition: EvtConExc.cc:3630
Referenced by init().
◆ initProbMax()
|
virtual |
Reimplemented from EvtDecayBase.
Definition at line 1757 of file EvtConExc.cc.
◆ islgr()
| bool EvtConExc::islgr | ( | double * | x, |
| double * | y, | ||
| int | n, | ||
| double | t | ||
| ) |
Definition at line 2834 of file EvtConExc.cc.
2835{ int n0=-1;
2836 double z;
2839 }
2840 double xstotal=y[599];
2841 if(n0==-1) {return 0;}else{
2842 double p1=y[n0]/xstotal;
2843 double p2=y[n0+1]/xstotal;
2845 if(p1<pm && pm<=p2) {return 1;}else{return 0;}
2846 }
2847}
◆ ISR_ang_integrate()
| double EvtConExc::ISR_ang_integrate | ( | double | x, |
| double | theta | ||
| ) |
Definition at line 2893 of file EvtConExc.cc.
2893 {
2894 //see Eq(11) in arXif:hep-ph/9910523, return h(theta)/h(pi)
2897 double cos2=costheta*costheta;
2898 double sin2=1-cos2;
2902 double hpi=L-1;
2903 double hq1= meE2/2*costheta/(sin2+meE2*cos2);
2904 double hq2= -0.5*log((1+sqrt(1-meE2)*costheta)/(1-sqrt(1-meE2)*costheta));
2906 double hq=(L-1)/2.+hq1+hq2+hq3;
2907 hq /= hpi; //Eq (11) in arXif:hep-ph/9910523
2908 return hq;
2909}
Referenced by ISR_ang_sampling().
◆ ISR_ang_sampling()
| double EvtConExc::ISR_ang_sampling | ( | double | x | ) |
Definition at line 2911 of file EvtConExc.cc.
2911 {
2912 double xx[180],yy[180];
2913 double dgr2rad=1./180.*3.1415926;
2914 xx[0]=0;
2915 xx[1]=5*dgr2rad; //first bin is 5 degree
2916 int nc=2;
2917 for(int i=6;i<=175;i++){ //last bin is 5 degree
2918 xx[nc]=i*dgr2rad;
2919 nc++;
2920 }
2921 xx[nc]=180*dgr2rad;
2922 for(int j=0;j<=nc;j++){
2923 yy[j]=ISR_ang_integrate(x,xx[j]);
2924 }
2926 int mypt=1;
2927 for(int j=1;j<=nc;j++){
2928 if(yy[j-1]<pm && pm<=yy[j]){mypt=j;break;}
2929 }
2930 pm= EvtRandom::Flat(0.,1);
2931 double mytheta= xx[mypt-1]+ (xx[mypt]-xx[mypt-1])*pm;
2932 return mytheta; //theta in rad unit
2933}
double ISR_ang_integrate(double x, double theta)
Definition: EvtConExc.cc:2893
Referenced by decay().
◆ lgr()
| double EvtConExc::lgr | ( | double * | x, |
| double * | y, | ||
| int | n, | ||
| double | t | ||
| ) |
Definition at line 2819 of file EvtConExc.cc.
2820{ int n0=-1;
2821 double z;
2824 }
2825 if(n0==-1) {return 0.0;}else{
2826 double k=(y[n0]-y[n0+1])/(x[n0]-x[n0+1]);
2828 }
2829 //std::cout<<"t="<<t<<" "<<x[n0+1]<<" "<<x[n0]<<std::endl;
2830 //std::cout<<"t="<<t<<" "<<y[n0+1]<<" "<<y[n0]<<std::endl;
2831 return z;
2832}
◆ Li2()
| double EvtConExc::Li2 | ( | double | x | ) |
Definition at line 2812 of file EvtConExc.cc.
2812 {
2813// double li2= -x +x*x/4. - x*x*x/9; // wangwp: may be not correct!
2815 return li2;
2816}
Referenced by SoftPhoton_xs().
◆ LLr()
| double EvtConExc::LLr | ( | double * | x, |
| double * | y, | ||
| int | n, | ||
| double | t | ||
| ) |
Definition at line 2849 of file EvtConExc.cc.
2850{ int n0=-1;
2851 double z;
2855 }
2856 if(n0==-1) {return 0.0;}else{
2857 double k=(y[n0]-y[n0+1])/(x[n0]-x[n0+1]);
2859 }
2860 //std::cout<<"t="<<t<<" "<<x[n0+1]<<" "<<x[n0]<<std::endl;
2861 //std::cout<<"t="<<t<<" "<<y[n0+1]<<" "<<y[n0]<<std::endl;
2862 return z;
2863}
◆ meson_sampling()
| bool EvtConExc::meson_sampling | ( | EvtVector4R | pcm, |
| EvtVector4R | pi | ||
| ) |
Definition at line 2546 of file EvtConExc.cc.
2546 {
2548 double theta=angles.getHelAng(1);
2549 double phi =angles.getHelAng(2);
2551
2553 double ang = 1 - costheta*costheta;
2554 if(pm< ang/1.) {return true;}else{return false;}
2555}
Referenced by hadron_angle_sampling().
◆ Mhad_sampling()
| double EvtConExc::Mhad_sampling | ( | double * | x, |
| double * | y | ||
| ) |
Definition at line 2865 of file EvtConExc.cc.
2865 {//sample ISR hadrons, return Mhadron
2866 //x=MH: array contains the Mhad
2867 //y=AF: array contains the Mhad*Xsection
2868 //n: specify how many bins in the hadron sampling
2871 double xrt=1;
2872 int mybin=1;
2875 if((y[i]/xst)<pm && pm<=(y[i+1]/xst)){
2876 mybin=i+1;
2877 break;
2878 }
2879 }
2880 pm= EvtRandom::Flat(0.,1);
2881 if(pm>1){std::cout<<"random larger than 1: "<<pm<<std::endl;}
2883
2884 if((mhd - _cms)>0){std::cout<<"selected mhd larger than Ecms: "<<mhd<<" > "<<x[mybin] <<std::endl;abort();}
2885 if(mhd<_mhdL){
2886 std::cout<<"the sample mhassrons is less than the low end of XS"<<mhd<<" <"<<_mhdL<<std::endl;
2888 abort();
2889 }
2890 return mhd;
2891}
Referenced by decay().
◆ mk_VXS()
| void EvtConExc::mk_VXS | ( | double | Esig, |
| double | Egamcut, | ||
| double | EgamH, | ||
| int | midx | ||
| ) |
Definition at line 3401 of file EvtConExc.cc.
3401 {
3402 //the mode index is put into vmode
3403 //initial
3404 //midx: mode index in vmode[midx] and VXS[midx][bin]
3405 int mode=vmode[midx];
3406 double uscale;
3407
3413 for(int i=0;i<600;i++){VXS[midx][i]=0;}
3414 return;
3415 }
3416
3417 int nint=50;
3418 int nmc=800;
3419 //double myxs0 =uscale*gamHXSection(s,Esig,Egamcut,nmc);
3423 myxs0 += sig_xs;
3424
3425 int Nstp=598;
3426 double stp=(EgamH - Egamcut)/Nstp;
3427 for(int i=0;i<Nstp;i++){//points within Egam(max) ~ Egamcut=25MeV
3428 double Eg0=EgamH - i*stp;
3429 double Eg1=EgamH - (i+1)*stp;
3430 //double xsi= uscale*gamHXSection(s,Eg1,Eg0,nmc);
3432 if(i==0) VXS[midx][0]=xsi;
3433 if(i>=1) VXS[midx][i]=VXS[midx][i-1]+xsi;
3434 }
3435 VXS[midx][598]=VXS[midx][597];
3436 VXS[midx][599]=VXS[midx][598]+ myxs0;
3437 //std::cout<<"mode "<<mode<<" "<<uscale<<" "<<VXS[midx][599]<<std::endl;
3438}
Referenced by init().
◆ narrowRXS()
| double EvtConExc::narrowRXS | ( | double | mxL, |
| double | mxH | ||
| ) |
Definition at line 3292 of file EvtConExc.cc.
3292 {
3293 //br for ee
3295 double kev2Gev=0.000001;
3296 psippee=0.262*kev2Gev;
3297 psipee =2.36*kev2Gev;
3298 jsiee =5.55*kev2Gev;
3299 phiee=4.266*0.001*2.954*0.0001;
3300 omegaee=0.6*kev2Gev;
3301 rhoee=7.04*kev2Gev;
3303 double sigv=0;
3304 double mv=0;
3305 double widee=0;
3306 double xpi=12*3.1415926*3.1415926;
3307 if(mxL<=3.0957 && 3.0957<=mxH || mxL<=3.0983 && 3.0983<=mxH) {
3308 widee = jsiee;
3309 mv = 3.096916;
3310 }else if(mxL<=3.6847 && 3.6847<=mxH || mxL<=3.6873 && 3.6873<=mxH){
3311 widee = psipee;
3312 mv = 3.686109;
3313 }else if(mxL<=1.01906 && 1.01906<=mxH || mxL<= 1.01986 && 1.01986<=mxH){
3314 widee = phiee;
3315 mv = 1.01946;
3316 }else{return -1.0;}
3317
3321 double unic = 3.89379304*100000; //in unit nb
3322 return sigv*unic; //vaccum factor has included in the Rad2
3323}
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA all others are from PDG *Only resonances with known widths into electron pairs are sept ************************************************************************C Declarations xmb DOUBLE PRECISION xcb DOUBLE PRECISION xmDcs DOUBLE PRECISION xmBbc DOUBLE PRECISION rh2ee DOUBLE PRECISION omepi DOUBLE PRECISION om3ee DOUBLE PRECISION phiee
Definition: RRes.h:37
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA all others are from PDG *Only resonances with known widths into electron pairs are sept ************************************************************************C Declarations xmb DOUBLE PRECISION xcb DOUBLE PRECISION xmDcs DOUBLE PRECISION xmBbc DOUBLE PRECISION rhoee
Definition: RRes.h:34
◆ OutStaISR()
| void EvtConExc::OutStaISR | ( | ) |
Definition at line 3738 of file EvtConExc.cc.
3738 {
3739 int ntot=myFisr.size();
3740 double mymu=0;
3741 for(int i=0;i<ntot;i++){mymu += myFisr[i];}
3742 mymu /= ntot;
3743 double mysig=0;
3744 for(int i=0;i<ntot;i++){ mysig += (myFisr[i]-mymu)*(myFisr[i]-mymu);}
3745 mysig /=ntot;
3746 mysig = sqrt(mysig);
3747 std::cout<<"========= Iteration times over ISR factor: "<<ntot<<std::endl;
3748 std::cout<<" ISR factor * VP factor= observedXS/BornXS: "<<mymu<<" + "<<mysig<<std::endl;
3749}
Referenced by ~EvtConExc().
◆ photonSampling()
| bool EvtConExc::photonSampling | ( | EvtParticle * | part | ) |
Definition at line 3261 of file EvtConExc.cc.
3261 {
3262 bool tagp,tagk;
3263 tagk=0;
3264 tagp=0;
3265 int nds = par->getNDaug();
3266 for(int i=0;i<par->getNDaug();i++){
3268 EvtVector4R p4i = par->getDaug(i)->getP4Lab();
3271
3272 if(pdgcode==111 ||pdgcode==221 || pdgcode==331){//for photon
3273 alpha = 1;
3274 }else continue;
3275
3279 double xratio = ang/angmax;
3281 //std::cout<<"pdgcode "<<pdgcode<<std::endl;
3282 //std::cout<<ang<<" "<<angmax<<" "<<xi<<" "<<xratio<<std::endl;
3283 if(xi>xratio) return false;
3284 }//loop over duaghters
3285 //if(tagp) std::cout<<" cosp "<<cosp<<std::endl;
3286 //if(tagk) std::cout<<" cosk "<<cosk<<std::endl;
3287
3288 return true;
3289}
Referenced by decay().
◆ Rad1()
| double EvtConExc::Rad1 | ( | double | s, |
| double | x | ||
| ) |
Definition at line 2568 of file EvtConExc.cc.
2568 { //first order correction
2569 //radiator correction upto second order, see Int. J. of Moder.Phys. A
2570 // Spectroscopy at B-factory using hard emision photon, by M. Benayoun, et al
2575 double L = 2*log(sme);
2577 return wsx;
2578}
Referenced by Rad1difXs().
◆ Rad1difXs()
| double EvtConExc::Rad1difXs | ( | EvtParticle * | p | ) |
Definition at line 2658 of file EvtConExc.cc.
2658 {// // first order xsection
2659 // where p -> hdarons(0 ~ _ndaugs-1) + gamma
2661 for(int i=1;i<_ndaugs;i++){
2663 }
2664
2666 double mrg = cms - summass;
2668 double mhs = pm*mrg + summass;
2669
2673
2674 // std::cout<<"pgam"<<pgam<<" mhs,egam,sin2theta,cms,x,wsx"<< mhs<<", "<<egam<<", "<<sin2theta<<", "<<cms<<", "<<x<<", "<<wsx<<std::endl;
2675
2678 double difxs = 2*mhs/cms/cms * wsx *xs;
2679
2681 differ *= mrg; //Jacobi factor for variable m_hds
2682 difxs *= mrg;
2683 return difxs;
2684}
Referenced by gamHXSection().
◆ Rad2()
| double EvtConExc::Rad2 | ( | double | s, |
| double | x | ||
| ) |
Definition at line 2580 of file EvtConExc.cc.
2580 {
2581 //radiator correction upto second order, see Int. J. of Moder.Phys. A, hep-ph/9910523, Eq(8)
2582 // Spectroscopy at B-factory using hard emision photon, by M. Benayoun, et al
2586 double xi2 = 1.64493407;
2587 double xi3=1.2020569;
2589 double L = 2*log(sme);
2591 double delta2 = (9./8.-2*xi2)*L*L - (45./16.-5.5*xi2-3*xi3)*L-6./5.*xi2*xi2-4.5*xi3-6*xi2*log(2.)+3./8.*xi2+57./12.;
2596 wsx = wsx + beta*beta/8. *wsx2;
2598 // return wsx*vph; //getVP(mymx);//vph is vaccum polarization factor
2600}
Referenced by narrowRXS(), Rad2difXs(), Rad2difXs(), Rad2difXs2(), Rad2difXs_er(), Rad2difXs_er2(), and Ros_xs().
◆ Rad2difXs() [1/3]
| double EvtConExc::Rad2difXs | ( | double | s, |
| double | x | ||
| ) |
Definition at line 2635 of file EvtConExc.cc.
2635 {// leading second order xsection
2639 //if(testflag==1) std::cout<<"s= "<<s<<" mhs= "<<mhs<<" x= "<<x<<" xs= "<<xs<<std::endl;
2641 double difxs = wsx *xs;
2643 return difxs;
2644}
◆ Rad2difXs() [2/3]
| double EvtConExc::Rad2difXs | ( | double | s, |
| double | x, | ||
| EvtXsection * | myxsection | ||
| ) |
Definition at line 2646 of file EvtConExc.cc.
2646 {// leading second order xsection
2650 //if(testflag==1) std::cout<<"s= "<<s<<" mhs= "<<mhs<<" x= "<<x<<" xs= "<<xs<<std::endl;
2652 double difxs = wsx *xs;
2654 return difxs;
2655}
◆ Rad2difXs() [3/3]
| double EvtConExc::Rad2difXs | ( | EvtParticle * | p | ) |
Definition at line 2604 of file EvtConExc.cc.
2604 {// leading second order xsection
2605 // where p -> hdarons(0 ~ _ndaugs-1) + gamma
2608 for(int i=1;i<_ndaugs;i++){
2611 }
2612
2615 double mrg = cms - summass;
2617 double mhs = pm*mrg + summass;
2618
2621 //double mhs = sqrt( s*(1-x) );
2623
2624 // std::cout<<"random : "<<pm<<std::endl;
2625
2628 double difxs = 2*mhs/cms/cms * wsx *xs;
2630 differ *= mrg; //Jacobi factor for variable m_hds
2631 difxs *= mrg;
2632 return difxs;
2633}
Referenced by gamHXSection(), and trapezoid().
◆ ReadVP()
| void EvtConExc::ReadVP | ( | ) |
Definition at line 3381 of file EvtConExc.cc.
3381 {
3382 vpx.clear();vpr.clear();vpi.clear();
3383 int vpsize=4001;
3384 vpx.resize(vpsize);
3385 vpr.resize(vpsize);
3386 vpi.resize(vpsize);
3387 std::string locvp=getenv("BESEVTGENROOT");
3388 locvp +="/share/vp.dat";
3389 ifstream m_inputFile;
3390 m_inputFile.open(locvp.c_str());
3391 double xx,r1,i1;
3392 double x1,y1,y2;
3393 xx=0;
3394 int mytag=1;
3395 for(int i=0;i<4001;i++){
3396 m_inputFile>>vpx[i]>>vpr[i]>>vpi[i];
3397 }
3398
3399}
std::ifstream ifstream
Referenced by init().
◆ resetResMass()
| void EvtConExc::resetResMass | ( | ) |
Definition at line 3083 of file EvtConExc.cc.
3083 {
3084 // EvtGen base class resetwidth has bugs, it make the mass of particle changed.
3086 EvtPDL::reSetMass(myres,mjsi);
3087 //EvtPDL::reSetWidth(myres,wjsi);
3088
3090 EvtPDL::reSetMass(myres,mpsip);
3091 //EvtPDL::reSetWidth(myres,wpsip);
3092
3094 EvtPDL::reSetMass(myres,mpsipp);
3095 //EvtPDL::reSetWidth(myres,wpsipp);
3096
3098 EvtPDL::reSetMass(myres,mphi);
3099 //EvtPDL::reSetWidth(myres,wphi);
3100
3102 EvtPDL::reSetMass(myres,momega);
3103 //EvtPDL::reSetWidth(myres,womega);
3104
3106 EvtPDL::reSetMass(myres,mrho0);
3107 //EvtPDL::reSetWidth(myres,wrho0);
3108
3110 EvtPDL::reSetMass(myres,mrho3s);
3111 //EvtPDL::reSetWidth(myres,wrho3s);
3112
3114 EvtPDL::reSetMass(myres,momega2s);
3115 //EvtPDL::reSetWidth(myres,womega2s);
3116}
Referenced by decay().
◆ Ros_xs()
| double EvtConExc::Ros_xs | ( | double | mx, |
| double | bree, | ||
| EvtId | pid | ||
| ) |
Definition at line 2718 of file EvtConExc.cc.
2718 {// in unit nb
2726 double nbar = 3.89379304*100000; // ! GeV-2 = 3.89379304*10^5 nbar
2727 return sigma*nbar;
2728}
Referenced by init().
◆ selectMass()
| double EvtConExc::selectMass | ( | ) |
Definition at line 3335 of file EvtConExc.cc.
3335 {
3336 double pm,mhdL,mhds;
3337 pm = EvtRandom::Flat(0.,1);
3338 mhdL =_mhdL;
3340 std::vector<double> sxs;
3341 for(int i=0;i<ISRID.size();i++){
3342 double ri=ISRXS[i]/AF[598];
3343 sxs.push_back(ri);
3344 }
3345 for(int j=0;j<sxs.size();j++){
3346 if(j>0) sxs[j] += sxs[j-1];
3347 }
3348 pm = EvtRandom::Flat(0.,1);
3349 if(pm<sxs[0]) {
3351 }
3352 for(int j=1;j<sxs.size();j++){
3353 double x0 = sxs[j-1];
3354 double x1 = sxs[j];
3356 }
3357 return mhds;
3358}
◆ selectMode()
| int EvtConExc::selectMode | ( | std::vector< int > | vmod, |
| double | mhds | ||
| ) |
Definition at line 3003 of file EvtConExc.cc.
3003 {
3004 //first get xsection for each mode in vmod array
3005 //std::cout<<"vmod.size, mhds "<<vmod.size()<<" "<<mhds<<std::endl;
3006 std::vector<int>excmod;
3007 excmod.push_back(0);
3008 excmod.push_back(1);
3009 excmod.push_back(2);
3010 excmod.push_back(6);
3011 excmod.push_back(7);
3012 excmod.push_back(12);
3013 excmod.push_back(13);
3014 excmod.push_back(45);
3015 excmod.push_back(46);
3016 double uscale = 1;
3017 std::vector<double> vxs;vxs.clear();
3018 for (int i=0;i<vmod.size();i++){
3019 int imod = vmod[i];
3023 //if(uscale!=1) std::cout<<"mode="<<imod<<" uscale= "<<uscale<<std::endl;
3024
3026
3027 if(i==0) {vxs.push_back(myxs);}
3028 else if(imod==74110){//for continuum process
3029 double xcon = myxs - vxs[i-1];
3030 if(xcon<0) {xcon=vxs[i-1];}else{xcon=myxs;}
3031 if(mhds<2.0) xcon=vxs[i-1]; //continuum cut: _cms energy less than 1.1, sampling phi,rho0 and omega resonance, see parp(2)=1.08 at Pythia.F
3032 vxs.push_back(xcon);
3033 }else{
3034 vxs.push_back(myxs+vxs[i-1]);
3035 }
3036 }//for_i_block
3037 //degugging
3038 // for(int i=0;i<vxs.size();i++){std::cout<<"Mhds="<<mhds<<" Mode="<<vmod[i]<<" xs_i= "<<vxs[i]<<" totalxs "<< vxs[vxs.size()-1]<<std::endl;}
3039
3040
3041 double totxs = vxs[vxs.size()-1];
3042 //if(totxs==0){std::cout<<"totalXs=0, vxs.size()= "<<vxs.size()<<std::endl;}
3043 int icount=0;
3044 int themode;
3045 mode_iter:
3047 if(pm <= vxs[0]/totxs) {
3048 themode= vmod[0];
3049 std::vector<EvtId> theDaug=get_mode(themode);
3050 EvtId daugs[100];
3051 for(int i=0;i<theDaug.size();i++){daugs[i]=theDaug[i];}
3053 if(_mode!=74110){return themode;}
3054 else if(exmode==-1){ return themode;}else{goto mycount;}
3055 }
3056
3057 for(int j=1;j<vxs.size();j++){
3058 double x0 = vxs[j-1]/totxs;
3059 double x1 = vxs[j]/totxs; //std::cout<<"j,x0,x1 "<<j<<" "<<x0<<" "<<x1<<std::endl;
3060 if(x0<pm && pm<=x1){
3061 themode=vmod[j];
3062 std::vector<EvtId> theDaug=get_mode(themode);
3063 EvtId daugs[100];
3064 for(int i=0;i<theDaug.size();i++){daugs[i]=theDaug[i];}
3066 if(_mode!=74110){return themode;}
3067 else if(exmode==-1){ return themode;} else{ break;}
3068 }
3069 }
3070
3071 mycount:
3072 icount++;
3073 if(icount<20000 ) goto mode_iter;
3074 //debugging
3075 //for(int i=0;i<vxs.size();i++){std::cout<<"Random="<<pm<<" Mode,Mhad "<<vmod[i]<<", "<<mhds<<" xs_i "<<vxs[i]<<" xsi/totalxs "<<vxs[i]/totxs<<std::endl;}
3076 std::cout<<"EvtConExc::selectMode can not find a mode with 20000 iteration for Mhadrons="<<mhds<<std::endl;
3077 return -10;
3078 //abort();
3079
3080}
static int inChannelList(EvtId parent, int ndaug, EvtId *daugs)
Definition: EvtDecayTable.cc:795
Referenced by decay().
◆ SetP4()
| void EvtConExc::SetP4 | ( | EvtParticle * | part, |
| double | mhdr, | ||
| double | xeng, | ||
| double | theta | ||
| ) |
Definition at line 2935 of file EvtConExc.cc.
2935 { //set the gamma and gamma* momentum according sampled results
2937 double phi=2*3.1415926*pm;
2939 double hdrE = sqrt(mhdr*mhdr+gamE*gamE);
2943 EvtVector4R p4[2];
2946
2947 for(int i=0;i<2;i++){
2949 p->getDaug(i)->init(myid,p4[i]);
2953 p->getDaug(i)->init(myid,PG);
2954 }
2955 }
2956 }
2957}
Referenced by decay().
◆ SetP4Rvalue()
| void EvtConExc::SetP4Rvalue | ( | EvtParticle * | part, |
| double | mhdr, | ||
| double | xeng, | ||
| double | theta | ||
| ) |
Definition at line 2959 of file EvtConExc.cc.
2959 { //set the gamma and gamma* momentum according sampled results
2961 double phi=2*3.1415926*pm;
2963 double hdrE = sqrt(mhdr*mhdr+gamE*gamE);
2967 EvtVector4R p4[2];
2970 for(int i=0;i<2;i++){
2972 p->getDaug(i)->init(myid,p4[i]);
2973 }
2974}
Referenced by decay().
◆ showResMass()
| void EvtConExc::showResMass | ( | ) |
Definition at line 3154 of file EvtConExc.cc.
3154 {
3155 std::cout<<"J/psi: "<<mjsi<<" "<<wjsi<<std::endl;
3156 std::cout<<"psipr: "<<mpsip<<" "<<wpsip<<std::endl;
3157 std::cout<<"psipp: "<<mpsipp<<" "<<wpsipp<<std::endl;
3158 std::cout<<"phi : "<<mphi<<" "<<wphi<<std::endl;
3159 std::cout<<"omega: "<<momega<<" "<<womega<<std::endl;
3160 std::cout<<"rho0 : "<<mrho0<<" "<<wrho0<<std::endl;
3161 std::cout<<"rho(3S)0: "<<mrho3s<<" "<<wrho3s<<std::endl;
3162 std::cout<<"omega(2S): "<<momega2s<<" "<<womega2s<<std::endl;
3163}
◆ SoftPhoton_xs()
| double EvtConExc::SoftPhoton_xs | ( | double | s, |
| double | b | ||
| ) |
Definition at line 2787 of file EvtConExc.cc.
2787 {
2791 double xi2 = 1.64493407;
2792 double xi3=1.2020569;
2794 double L = 2*log(sme);
2796 double delta2 = (9./8.-2*xi2)*L*L - (45./16.-5.5*xi2-3*xi3)*L-6./5.*xi2*xi2-4.5*xi3-6*xi2*log(2.)+3./8.*xi2+57./12.;
2799
2800 double beta2 = beta*beta;
2801 double b2 = b*b;
2802
2803 double xs=(-32*b*beta + 8*pow(b,2)*beta - 10*b*pow(beta,2) + pow(b,2)*pow(beta,2) + 32*pow(b,beta)*Delta -
2804 6*(3 - 4*b + pow(b,2))*pow(beta,2)*log(1 - b) - 32*b*pow(beta,2)*log(b) + 8*pow(b,2)*pow(beta,2)*log(b) +
2805 16*pow(beta,2)*Li2(b))/32. ;
2807 //return xs*vph; //getVP(mymx); //vph :the vaccum polarzation factor
2809
2810}
◆ split()
| std::vector< std::string > EvtConExc::split | ( | std::string | str, |
| std::string | pattern | ||
| ) |
Definition at line 3630 of file EvtConExc.cc.
3631{
3632 std::string::size_type pos;
3633 std::vector<std::string> result;
3634 str+=pattern;
3635 int size=str.size();
3636
3637 for(int i=0; i<size; i++)
3638 {
3639 pos=str.find(pattern,i);
3640 if(pos<size)
3641 {
3642 std::string s=str.substr(i,pos-i);
3643 result.push_back(s);
3644 i=pos+pattern.size()-1;
3645 }
3646 }
3647 return result;
3648}
Referenced by InitPars().
◆ sumExc()
| double EvtConExc::sumExc | ( | double | mx | ) |
Definition at line 3179 of file EvtConExc.cc.
3179 {//summation all cross section of exclusive decays
3180 std::vector<int> vmod; vmod.clear();
3181 int mn[83]={0,1,2,3,4,5,6,7,8,9,10,11,13,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,32,33,34,35,36,37,38,40,41,44,45,46, // 12,14, 30,31,39,42 and 43 is removed
3182 50,51,
3183 59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,
3184 90,91,92,93,94,95,96,
3185 74100,74101,74102,74103,74104,74105,74110};
3186 int mn2[84]={0,1,2,3,4,5,6,7,8,9,10,11,12,13,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,32,33,34,35,36,37,38,40,41,44,45,46,//mode 14,30,31,42, 43 are removed
3187 50,51,
3188 59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,
3189 90,91,92,93,94,95,96,
3190 74100,74101,74102,74103,74104,74105,74110};
3191
3193 for(int i=0;i<83;i++){vmod.push_back(mn[i]);}
3194 }else{
3195 for(int i=0;i<84;i++){vmod.push_back(mn2[i]);}
3196 }
3197
3198 // for(int i=0;i<79;i++){vmod.push_back(mn2[i]);}
3199
3200 double xsum=0;
3201 double uscale = 1;
3202 for(int i=0;i<vmod.size();i++){
3203 int mymode = vmod[i];
3204 if(mymode ==74110) continue;
3207 init_mode(mymode);
3210 //if(mx==3.0967) {std::cout<<vmod[i]<<" "<<uscale<<" "<<xsum<<std::endl;}
3211 }
3212 return xsum;
3213}
Referenced by init().
◆ trapezoid() [1/2]
| double EvtConExc::trapezoid | ( | double | s, |
| double | a, | ||
| double | b, | ||
| int | n | ||
| ) |
Definition at line 3752 of file EvtConExc.cc.
3753{
3756 double sum = 0.0;
3759 {
3760 sum += (gaps/2.0) * (Rad2difXs(s,xL + i*gaps, staxsection) + Rad2difXs(s,xL + (i+1)*gaps,staxsection) );
3761 }
3762 return sum;
3763}
◆ trapezoid() [2/2]
| double EvtConExc::trapezoid | ( | double | s, |
| double | El, | ||
| double | Eh, | ||
| int | n, | ||
| EvtXsection * | myxs | ||
| ) |
Definition at line 3765 of file EvtConExc.cc.
3766{
3769 double sum = 0.0;
3772 {
3774 }
3775 return sum;
3776}
◆ VP_sampling()
| bool EvtConExc::VP_sampling | ( | EvtVector4R | pcm, |
| EvtVector4R | pi | ||
| ) |
Definition at line 2557 of file EvtConExc.cc.
2557 {
2559 double theta=angles.getHelAng(1);
2560 double phi =angles.getHelAng(2);
2562
2564 double ang = 1 + costheta*costheta;
2565 if(pm< ang/2.) {return true;}else{return false;}
2566}
Referenced by hadron_angle_sampling().
◆ writeDecayTabel()
| void EvtConExc::writeDecayTabel | ( | ) |
Definition at line 3468 of file EvtConExc.cc.
3468 {
3471 ob.open("obsxs.dat");
3472 //
3474
3475 double xscon=VXS[im][599];
3476 //std::cout<<"tsize "<<tsize<<" "<<xscon<<" "<<VXS[70][599]<<std::endl;
3477 std::vector<int> Vmode;
3478 Vmode.push_back(6);//1:pi+pi-
3479 Vmode.push_back(13);//2:pi+pi-2pi0
3480 Vmode.push_back(12);//3:2(pi+pi-)
3481 Vmode.push_back(0);//4:ppbar
3482 Vmode.push_back(1);//5:nnbar
3483 Vmode.push_back(45);//6:K+K-
3484 Vmode.push_back(46);//7:K0K0bar
3485 Vmode.push_back(7);//8:pi+pi-pi0
3486 Vmode.push_back(2);//9:Lambda anti-Lambda
3487 Vmode.push_back(37);//10: pi+pi- eta
3488 std::vector<std::string> vmdl;
3489 vmdl.push_back("PHOKHARA_pipi");
3490 vmdl.push_back("PHOKHARA_pi0pi0pipi");
3491 vmdl.push_back("PHOKHARA_4pi");
3492 vmdl.push_back("PHOKHARA_ppbar");
3493 vmdl.push_back("PHOKHARA_nnbar");
3494 vmdl.push_back("PHOKHARA_KK");
3495 vmdl.push_back("PHOKHARA_K0K0");
3496 vmdl.push_back("PHOKHARA_pipipi0");
3497 vmdl.push_back("PHOKHARA_LLB");
3498 vmdl.push_back("PHOKHARA_pipieta");
3499
3500 ob<<"mode_index observed cross /nb"<<std::endl;
3501 for(int i=0;i<vmode.size();i++){
3503//cout<<"vmode["<<i<<"] = "<<vmode[i]<<", VXS["<<getModeIndex(vmode[i])<<"][599] = "<<XS[getModeIndex(vmode[i])][599]<<std::endl; // wangwp
3504 }
3505
3506 for(int i=0;i<Vmode.size();i++){
3507 //std::cout<<"i: "<<xscon<<" "<<VXS[Vmode[i]][599]<<std::endl;
3508 xscon -= VXS[ getModeIndex(Vmode[i]) ][599];
3509 }
3510 //debugging
3511 for(int i=0;i<Vmode.size();i++){
3512 // std::cout<<Vmode[i]<<"-th mode: "<<VXS[getModeIndex(Vmode[i])][599]<<" "<<std::endl;
3513 }
3514
3544 for(int i=0;i<Vmode.size();i++){
3546 }
3549}
Referenced by init().
◆ xs_sampling() [1/2]
| bool EvtConExc::xs_sampling | ( | double | xs | ) |
Definition at line 2517 of file EvtConExc.cc.
2517 {
2519 //std::cout<<"Rdm= "<<pm<<std::endl;
2521}
◆ xs_sampling() [2/2]
| bool EvtConExc::xs_sampling | ( | double | xs, |
| double | xs1 | ||
| ) |
Definition at line 2523 of file EvtConExc.cc.
2523 {
2525 // std::cout<<"Rdm= "<<pm<<std::endl;
2526 if(pm <xs/(xs0*1.1)) {return true;} else {return false;}
2527}
Member Data Documentation
◆ _cms
|
static |
Definition at line 190 of file EvtConExc.hh.
Referenced by baryon_sampling(), calAF(), checkEvtRatio(), decay(), difgamXs(), Egam2Mhds(), findMaxXS(), getObsXsection(), init(), ISR_ang_integrate(), Mhad_sampling(), mk_VXS(), narrowRXS(), Rad2difXs(), Rad2difXs2(), Rad2difXs_er(), Rad2difXs_er2(), selectMass(), SetP4(), SetP4Rvalue(), sumExc(), and writeDecayTabel().
◆ conexcmode
|
static |
Definition at line 211 of file EvtConExc.hh.
Referenced by decay(), and EvtRexc::decay().
◆ getMode
|
static |
Definition at line 194 of file EvtConExc.hh.
Referenced by init().
◆ mlow
|
static |
Definition at line 234 of file EvtConExc.hh.
Referenced by EvtDecay::energySpread(), init(), and mk_VXS().
◆ multimode
|
static |
Definition at line 211 of file EvtConExc.hh.
Referenced by EvtRexc::decay(), and init().
◆ mup
|
static |
Definition at line 234 of file EvtConExc.hh.
Referenced by EvtDecay::energySpread(), and init().
◆ myxsection
|
static |
Definition at line 189 of file EvtConExc.hh.
Referenced by decay(), difgamXs(), findMaxXS(), init(), init_mode(), mk_VXS(), Rad1difXs(), Rad2difXs(), Rad2difXs(), Rad2difXs2(), Rad2difXs_er(), Rad2difXs_er2(), selectMode(), sumExc(), and ~EvtConExc().
◆ SetMthr
|
static |
Definition at line 192 of file EvtConExc.hh.
Referenced by decay(), init(), and EvtDecay::initialize().
◆ staxsection
|
static |
Definition at line 189 of file EvtConExc.hh.
Referenced by calAF(), energySpread(), init(), trapezoid(), and ~EvtConExc().
◆ XS_max
|
static |
Definition at line 191 of file EvtConExc.hh.
Referenced by init(), Rad1difXs(), Rad2difXs(), Rad2difXs(), Rad2difXs2(), and xs_sampling().
The documentation for this class was generated from the following files:
