FDC軟件包在分波分析中的應(yīng)用

FDC軟件包在分波分析中的應(yīng)用

簡(jiǎn)單介紹FDC系統(tǒng)強(qiáng)相互作用的有效拉氏量的構(gòu)造規(guī)則FDC-PWA分波分析版本的簡(jiǎn)介簡(jiǎn)單演示總結(jié)BESIII物理分析講習(xí)班,2008年6月25日北京王建雄中國(guó)科學(xué)院高能物理研究所,理論室

簡(jiǎn)單介紹FeynmanDiagramCalculation(FDC),1993年開(kāi)始.PROGRESSINFDCPROJECT.Nucl.Instrum.Meth.A534:241-245,2004

FDC的主頁(yè)：所內(nèi)版：8/fdc/index.html(可以看到最新的結(jié)果)所外版：/lunwen/wjx/public_html/index.html

(更新不定期)

FDC-SM-and-Many-ExtensionsFDC-NRQCDFDC-MSSMFDC-EMTFDC-PWAFDC-LOOP

FDC-PWA1998開(kāi)始,鄒冰松,祝玉燦,..Thereisacommonkernel.TheSourceiswritteninREDUCE,RLISP.分波分析在BES物理分析中的重要性確定各種量子數(shù),質(zhì)量,寬度復(fù)雜結(jié)構(gòu)的譜、共振態(tài)、奇特態(tài)、膠球,以前應(yīng)用分波分析做理論計(jì)算，物理分析，編寫(xiě)Fortran程序都是手做的，就算不包括3/2、2、5/2,7/2自旋態(tài)，也得花幾個(gè)月的時(shí)間，如果包括了3/2、2、5/2,7/2自旋態(tài)，計(jì)算量和表達(dá)式會(huì)大幾十倍。另外人工做這些工作怎樣檢查錯(cuò)誤也很困難。而FDC系統(tǒng)，其中各個(gè)部份，各種核心算法在這么多年的發(fā)展和應(yīng)用中，已經(jīng)進(jìn)行了很多次的檢驗(yàn)。（例如手算與FDC算對(duì)比，F(xiàn)DC與KEK的Grace對(duì)比等）。2.FDC-PWA:FDC在BEPC-II物理分析中的應(yīng)用在FDC系統(tǒng)中實(shí)現(xiàn)了自動(dòng)化構(gòu)造分波分析有效拉氏量模型的部分（與鄒冰松，祝玉燦合作），以及產(chǎn)生包括3/5,5/2,7/2及整數(shù)自旋態(tài)的物理過(guò)程的費(fèi)曼圖，產(chǎn)生所有Fortran源程序以及最大似然法的擬合程序的部分，在其中結(jié)合最大似然法擬合程序的算法，對(duì)于矩陣元平方的算法在結(jié)構(gòu)上進(jìn)行了兩次大的優(yōu)化,使得產(chǎn)生的程序在實(shí)驗(yàn)分析擬合中計(jì)算速度兩次都有10倍的提高2.1強(qiáng)相互作用的有效拉氏量的構(gòu)造規(guī)則LorentzInvarianceC-parityconservationP-parityconservationCPconservationH=HABCDecayDecayABCDThecoefficientsofalltheindepenttermsareconstantThecoefficientsofalltheindepenttermsdependontwovaribales每個(gè)系數(shù)都有獨(dú)立的復(fù)數(shù)相因子.每個(gè)傳播子都有off-shell貢獻(xiàn)的壓低因子對(duì)于一個(gè)粒子列表,給出所有可能的有效相互作用頂角.ABC(offmassshell)ForNNscatting,T-channelverticesmustbeconsidered:TheCoefficientsofindependenttermsdependonthe:2

t=(PA

–

PB)

Thepropagatorforresonance:SuppersionFactorforinteractionverticeswhentheresonanceisoffmassshell:

Barrier,

ThereisanoptiontochoosethesuppersionfactorOptiontoapplyIsopinsymmetry:

GlobelOptionsinModelGeneration(infile:addvertices)%just_v_type:='ok;

Itwillgenerateallthetypeoftheinteractionverticesasreqiredandwithouttherealverticesifthelinewasuncommented.

%direct_model:='ok;

Thelineisusedtochosethetwodifferentkindsofeffectivelagrangian.Thedefaultischoosedtousethepure-orbital-angular-momentumcovarianttensorsandthecovariantspinwavefunctionforcontructionoftheeffectivelagrangian.Touncommentitmeanstochooseanaiveindependentlorentzsetofthecombinationofmomentum,spinwavefunction.

Thedifferenceofthesetwomodelsisfromtheoff-shellpartoftheresonaceparticles.%final_particle_list:='("p""\pi""K""eta""\pi^0""K^0");

Thislineisusedtoputtheparticlesintothelistofexternalparticles.ThelistisusedintheFeynmandiagramgenerationprocedurebyFDCandalltheparticlesinthelistwillnotappearedintheresonaceplace.

%keep_sm_part_in_tex:='ok;

TheoptionisusedtochoseiftheinteractionverticesintheStandardmodelarekeptinthe"model.tex"fileornot.Thedefaultiskeep_sm_part_in_tex:='nil;anditmeanstheverticesintheSMaredroped.

http:='ok;

Thisoptionisusedtochosethetypeof"model.tex"filewhichisgeneratedbyFDC.http:='ok;meansthe"model.tex"isforhttpversion.Thedefaultishttp:='nil;anditmeansthe"model.tex"isforPSversion.n_particles:='3;

Theisthelimittothenumberofparticlesineachvertex.Thesmallestnumberis3andthedefaultis3.Itcanbechangedinto4orbiggeranditmeansthatallpossibletheoneparticlesdecaysintothreeparticlesverticeswillbegeneratedifthenumberissetto4.

n_total_spin:='6;

Thisnumberisthelimittothesumofthespinofeachparticlesinaneffectiveinteractionvertex.Thedefaultis6.itmeansthattheFDCwilluse6forthe"n_totla_spin"ifitisnotgiven.

2.2FDC-PWA分波分析版本的簡(jiǎn)介有效拉氏量模型的自動(dòng)化產(chǎn)生費(fèi)曼圖的產(chǎn)生>diag矩陣元及平方的Fortran程序的產(chǎn)生>amp相空間程序及Likelihoodfitting的Fortran程序的產(chǎn)生>kine進(jìn)入Fortran目錄和編譯>cdfort>make開(kāi)始進(jìn)行Likelihood擬合>fit

FDC3.1,1999,開(kāi)始在BES的分波分析中應(yīng)用,模型不是完全自動(dòng)化產(chǎn)生的,并且沒(méi)有包括末態(tài)是純波色子過(guò)程FDC3.32003開(kāi)始應(yīng)用,模型產(chǎn)生完全的自動(dòng)化,包括末態(tài)是純波色子過(guò)程FDC-PWAApplication:

>gmodel>doall>cdfort>make>fit

FDCSystemNumericalresults&RequiredplotsinHbookformatPhysicalmodelevaluationFullFeynmanRulesPhysicalparametersLatexversionofthenmodelInputforphysicalmodelNumericalintegration,fittingTreatmentofkinematicsAmplitudemanipulationFeynmandiagramgenerationInputforphysicalprocessChoosephysicalmodelandprocessesFIG.1:FDCsystemflowchartTopreparephenomenologicalmodelListofallthenon-elementaryparticlesandtheirquantumnumbersThegeneratedphysicalmodelforFDCsystemLatexversionofthegeneartedphysicalmodelListofalltheparticlesListofallthepropagatorsListofalltheinteractionverticesConstructallthepossibleinteractionverticesfromalltheparticlesbyapplyingthefollowingconditions:H?=H,Lorentzinvariance,CPinvariance,Pinvariance,Cinvariance,Isospininvariance,Baryonnumberconservation,...StandardModelinputwithoutQCDFIG.3:SystemflowchartforphysicalprocessPhysicalProcessUsersshouldpreparetwofiles:pdata1.dat–experimenteventsdatafilepdata1.mc-phasespacemontecarloeventfileControlflagandparametersfilesgeneratedbyFDCwhichcanbechangedlaterbyusers:flag.inp,amptable.inp,fpara.inp,reson.inpOutput:,pep.res,mplot.hbook,dplot.hbookInputforaphysicalprocess:physicalmodelcanbechosenManyoptions,histograms,scatterplotscanbedemanded.GenerateFeynmanDiagramManipulateamplitudesforeachdiagramandgenerateFORTRAN77sourceforcalculationofamplitudesandtheirsquareFORTRAN77sourcetodolikelihoodfittingforallthefreeparametersthatwereintroducedinphysicalmodelFindandproperlytreatalltheresonnance,t-channelsingularities,...andgenerateFORTRAN77sourceforphasespaceintegralCompileFORTRAN77programsandrun'fit'forlikelihoodfitting系統(tǒng)的檢驗(yàn):1)FDC系統(tǒng)內(nèi)部的檢驗(yàn),各種核心程序:三種矩陣元,矩陣元平方算法的相互檢驗(yàn),