Rate_14Ngp_ThrCut/0000775000372200007340000000000012640052611014134 5ustar suleimanitsopsRate_14Ngp_ThrCut/csN14_interpolation.C0000644000372200007340000000207612640046616020112 0ustar suleimanitsops#include "Math/Polynomial.h" #include "Math/Interpolator.h" #include #include #include "TCanvas.h" #include "TGraph.h" Double_t csN14_interpolation(Double_t k){ using namespace std; const Int_t n1 = 94; Double_t k1[n1], cs1[n1]; // ifstream txtfile1; char file1[200]; sprintf(file1,"csN14.dat"); txtfile1.open(file1); // 14N(g,p)13C // E(MeV) sigma(barn) // Claudio cross section file. // // cout << "Data file is " << file1 << endl; for ( Int_t i = 0; i < n1; ++i ){ txtfile1 >> k1[i] >> cs1[i]; cs1[i] = cs1[i] *1e+9; // in nbarn (1nb = 1e-33 cm2) // cout << "energy = " << k1[n1] << " Cross Section = " << cs1[n1] << endl; } //// // You can choose among the following methods: // CSPLINE, LINEAR, POLYNOMIAL, // CSPLINE_PERIODIC, AKIMA, AKIMA_PERIODIC ROOT::Math::Interpolator inter(n1, ROOT::Math::Interpolation::kLINEAR); inter.SetData(n1, k1, cs1); // cout << " k =" << 8.0 << " Interpolation Value = " << inter.Eval(8.0) << endl; return inter.Eval(k); //// } Rate_14Ngp_ThrCut/csN14_interpolation_C.so0000775000372200007340000007046712640051500020615 0ustar suleimanitsopsELF'4K4 (FFFVVlGWW$$Ptd>>>QtdGNUIS{mË2w? DC9   ( @!jD X (\M\%F//% P9<=>?ABCEIJLMNOQSUVXYZ[\]_`abegijlorsuvwx{:K.^a-jUu526'/%+faN@QjR.J#;**Fs*ERz :E+sfulR bx_ x롢i$U>guK#䔡 uDO f)Pt͌3BE3hs-,6]P5o9K?(<+]|Khf}_e ! hZ + 2 8" __gmon_start___init_fini__cxa_finalize_Jv_RegisterClasses_ZN10TCanvasImp4LockEv_ZN10TCanvasImp6UnlockEv_ZN10TCanvasImp8IsLockedEv_ZN10TCanvasImpD2Ev_ZTV10TCanvasImp_ZN10TCanvasImpD1Ev_ZN10TCanvasImp5CloseEv_ZN10TCanvasImp11ForceUpdateEv_ZN10TCanvasImp7IconifyEv_ZN10TCanvasImp10InitWindowEv_ZN10TCanvasImp4ShowEv_ZNK10TCanvasImp9HasEditorEv_ZNK10TCanvasImp10HasMenuBarEv_ZNK10TCanvasImp12HasStatusBarEv_ZNK10TCanvasImp10HasToolBarEv_ZNK10TCanvasImp11HasToolTipsEv_ZN10TCanvasImp17GetWindowGeometryERiS0_RjS1__ZN10TCanvasImp13SetStatusTextEPKci_ZN10TCanvasImp17SetWindowPositionEii_ZN10TCanvasImp13SetWindowSizeEjj_ZN10TCanvasImp14SetWindowTitleEPKc_ZN10TCanvasImp13SetCanvasSizeEjj_ZN10TCanvasImp11ShowMenuBarEb_ZN10TCanvasImp13ShowStatusBarEb_ZN10TCanvasImp11RaiseWindowEv_ZN10TCanvasImp12ReallyDeleteEv_ZN10TCanvasImp10ShowEditorEb_ZN10TCanvasImp11ShowToolBarEb_ZN10TCanvasImp12ShowToolTipsEbG__cpp_dllrevcsN14_interpolation_C_ACLiC_dict_ZN48G__Sizep2memfunccsN14_interpolation_C_ACLiC_dict13sizep2memfuncEvG__cpp_setup_inheritancecsN14_interpolation_C_ACLiC_dictG__cpp_setup_memvarcsN14_interpolation_C_ACLiC_dictG__cpp_setup_memfunccsN14_interpolation_C_ACLiC_dictG__cpp_reset_tagtablecsN14_interpolation_C_ACLiC_dictG__csN14_interpolation_C_ACLiC_dictLN_TVectorTlEdoublegRG__csN14_interpolation_C_ACLiC_dictLN_TVectorTlEfloatgRG__csN14_interpolation_C_ACLiC_dictLN_vectorlEROOTcLcLTSchemaHelpercOallocatorlEROOTcLcLTSchemaHelpergRsPgRG__csN14_interpolation_C_ACLiC_dictLN_reverse_iteratorlEvectorlEROOTcLcLTSchemaHelpercOallocatorlEROOTcLcLTSchemaHelpergRsPgRcLcLiteratorgRG__csN14_interpolation_C_ACLiC_dictLN_vectorlETVirtualArraymUcOallocatorlETVirtualArraymUgRsPgRG__csN14_interpolation_C_ACLiC_dictLN_reverse_iteratorlEvectorlETVirtualArraymUcOallocatorlETVirtualArraymUgRsPgRcLcLiteratorgRG__csN14_interpolation_C_ACLiC_dictLN_vectorlEcomplexlEdoublegRcOallocatorlEcomplexlEdoublegRsPgRsPgRG__csN14_interpolation_C_ACLiC_dictLN_reverse_iteratorlEvectorlEcomplexlEdoublegRcOallocatorlEcomplexlEdoublegRsPgRsPgRcLcLiteratorgRG__csN14_interpolation_C_ACLiC_dictLN_pairlEunsignedsPintcOintgRG__csN14_interpolation_C_ACLiC_dictLN_vectorlEpairlEunsignedsPintcOintgRcOallocatorlEpairlEunsignedsPintcOintgRsPgRsPgRG__csN14_interpolation_C_ACLiC_dictLN_reverse_iteratorlEvectorlEpairlEunsignedsPintcOintgRcOallocatorlEpairlEunsignedsPintcOintgRsPgRsPgRcLcLiteratorgRG__set_cpp_environmentcsN14_interpolation_C_ACLiC_dict_ZNK10TCanvasImp3IsAEv_ZN10TCanvasImp5ClassEv_ZN10TCanvasImpD0Ev_ZdlPv_ZN49G__cpp_setup_initcsN14_interpolation_C_ACLiC_dictD2EvG__remove_setup_func_ZN49G__cpp_setup_initcsN14_interpolation_C_ACLiC_dictD1EvG__cpp_setup_tagtablecsN14_interpolation_C_ACLiC_dictG__get_linked_tagnum_fwdG__cpp_setup_funccsN14_interpolation_C_ACLiC_dictG__lastifuncpositionG__defined_typenameG__memfunc_setupG__resetifuncpositionG__cpp_setup_globalcsN14_interpolation_C_ACLiC_dictG__resetplocalG__resetglobalenvG__cpp_setup_typetablecsN14_interpolation_C_ACLiC_dictG__search_typename2G__setnewtypeG__get_linked_tagnum_Z52G__get_sizep2memfunccsN14_interpolation_C_ACLiC_dictvG__setsizep2memfuncG__cpp_setupcsN14_interpolation_C_ACLiC_dictG__check_setup_versionG__getsizep2memfunc_ZN13TVersionCheckC1Ei_ZNSt8ios_base4InitC1Ev_ZNSt8ios_base4InitD1Ev__cxa_atexit_ZN4ROOT14RegisterModuleEvG__add_setup_funcG__call_setup_funcsG__cpp_setup_initializercsN14_interpolation_C_ACLiC_dict_Unwind_Resume_Z19csN14_interpolationd_ZNSt8ios_baseC2Ev_ZTVSt9basic_iosIcSt11char_traitsIcEE_ZTTSt14basic_ifstreamIcSt11char_traitsIcEE_ZNSt9basic_iosIcSt11char_traitsIcEE4initEPSt15basic_streambufIcS1_E_ZTVSt14basic_ifstreamIcSt11char_traitsIcEE_ZNSt13basic_filebufIcSt11char_traitsIcEEC1Ev_ZNSt13basic_filebufIcSt11char_traitsIcEE4openEPKcSt13_Ios_Openmode_ZNSt9basic_iosIcSt11char_traitsIcEE5clearESt12_Ios_Iostate_ZNSi10_M_extractIdEERSiRT__ZN4ROOT4Math12InterpolatorC1EjNS0_13Interpolation4TypeE_ZN4ROOT4Math12Interpolator7SetDataEjPKdS3__ZNK4ROOT4Math12Interpolator4EvalEd_ZN4ROOT4Math12InterpolatorD1Ev_ZTVSt13basic_filebufIcSt11char_traitsIcEE_ZNSt13basic_filebufIcSt11char_traitsIcEE5closeEv_ZNSt12__basic_fileIcED1Ev_ZTVSt15basic_streambufIcSt11char_traitsIcEE_ZNSt6localeD1Ev_ZNSt8ios_baseD2Ev_ZNSt14basic_ifstreamIcSt11char_traitsIcEED1Ev_ZNSt13basic_filebufIcSt11char_traitsIcEED1Ev_ZNSt9basic_iosIcSt11char_traitsIcEED2Ev_ZNSt15basic_streambufIcSt11char_traitsIcEED2Ev__gxx_personality_v0G__doubleG__letdouble_ZGVZN4ROOT4Math4Util7EvalLogEdE7epsilon_ZTI10TCanvasImp_ZN10TCanvasImp11ShowMembersER16TMemberInspector_ZN10TCanvasImp8StreamerER7TBuffer_ZTS10TCanvasImp_ZTVN10__cxxabiv117__class_type_infoElibCint.solibCore.solibMathMore.solibstdc++.so.6libgcc_s.so.1libc.so.6_edata__bss_start_endcsN14_interpolation_C.soGCC_3.0CXXABI_1.3GLIBCXX_3.4.9GLIBCXX_3.4GLIBC_2.1.3 P&y @ӯk)t)si W WX ZZZ$Z,Z4ZYf YuY5YnY6Y7dZ.,YQ0Yl4Y8YU%E ݝ$ƅDžƅADžDž DžP@DžBDžD$@ $J $B B D$$D$DžcsN1Dž4.dafDž tD$$9D$@ $R1  D$ $$T$ ^uuD$D$^4$m|$ D$^4$D$݅\$4$ ݝ4$B B $@$$$ @DžA$݅[^_]Ë PB$D$alj4$ $[<$拍lj $?‹B Dž $<$뼉Nj$먉Ǎ@$iҍ&US4EZ"P$PT$P T$PT$ PT$PT$PT$P T$P$T$ @(D$$$ED$d$\$4[]ÐU]ÐU]ÐU1]Ð1[!U勁PE]ÐU]ÐU]ÐU]ÐU1]ÐU]ÐU1]ÐU1]ÐU1]ÐU1]ÐU1]ÐUEE EE1]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐU]ÐUS× []ÐUSwE$[]ÐUSG4$p[]Ë $ÐUVS_t&vЋu[^]US[PY[csN14_interpolation_C_ACLiC_dictreverse_iteratorvector >vector >G__cpp_setupcsN14_interpolation_C_ACLiC_dict()vector >reverse_iterator >::iterator>vector >reverse_iterator >::iterator>vector,allocator > >reverse_iterator,allocator > >::iterator>vector,allocator > >reverse_iterator,allocator > >::iterator>Double_td - 'Double_t' 0 - kcsN14_interpolationDouble 8 bytesvectorreverse_iteratorvectorvector >vector >pairvector >TVectorTTVectorTTVectorTTVectorTpair10TCanvasImp(knN;/@`Ptt$TD4dD$ 4@D`Tdt @`$4D T@d`t ,:zPR|  ,AB A @AB A ` AD A LB O AB A AB A AB A AD A  AB A @AD A `AD A AD A |AD A lAD A \+AB g lAB A  \AB A @LAB A `<AB A ,AB A AB A  AB A AB A AB A  AB A @AB A `AB A  AG A  AG A |AB A lAB A \AB A  LLB   @AB AWA  dAB AWA  .AB AhA  &AB A`A  \AB AA  AB AA  "AB A\A ,<AB CAAA  l )AB AcA  ,kAB AeA (xAB Feu AzPLR|5 0$_AB C AAA A 4AB G}A XZJ;p2 Wit # 8:o@    Y "HoXoo`o WW&$6$F$V$f$v$$$$$$$$$%%&%6%F%V%f%v%%%%%%%%%&&&&6&F&V&f&v&&&&&&&&&''&'6'F'V'f'v'''r>c>c;cL;c;c;c(<c`<c>s<c<cGCC: (GNU) 4.4.7 20120313 (Red Hat 4.4.7-9)GCC: (GNU) 4.4.7 20120313 (Red Hat 4.4.7-11).symtab.strtab.shstrtab.note.gnu.build-id.gnu.hash.dynsym.dynstr.gnu.version.gnu.version_r.rel.dyn.rel.plt.init.text.fini.rodata.eh_frame_hdr.eh_frame.gcc_except_table.ctors.dtors.jcr.data.rel.ro.dynamic.got.got.plt.data.bss.comment$.oh8 @@ @ Ho``UoXXd Hm  " " v##0q$$|''8:8:T:T:h>>@@@@lFFNVF WGWG W G WGXHd Y I Z J\hZhJ0hJYJxP 4 h[@ `X "# $ ' 8: T: >@@FVWW WWX Y ZhZ VW(W5' KhZZlZh0( tWFW: 7 p2 OqZbrZq W~ WdZi( W9  Y 8" ( D( z"7" ( 9" > ! Im7" 9." (hZ/*" c8" @9" p) CY@8" x 3 9&" ue7!IFX }ȟ8 `!!M Y$hb#ֲ /c3"-5R7β^?rBwق1 l( Ԑ.(j8 'pT%,`2Tfcr?kbj-YmP' wR!Hx%i*B-@8M?*q恡`:IكxfG9hihʗiE tPZIPZ*i,ځD0Z[rZ)Az `hGDTi^ +ZKXVG!h vS.ĩF5˃&9*8z)L{k%ƚ YKBʫ(›٥> 1$QJU)c37$Ǵ6-F Jj4PRWMSg-n?Zcmu`-d]fgsͲu- p+H А s=7[~w*_1ijCgșҁJ ;qX?$3mqa>r`1 cCjz9dp0l2uOMkn(4G/Gs]kj.ӮFmRZk`uPKrPQ*F# lN:GoU ~+[.">cƓDf:$B8}%ݪ|J E@x/ry" DxJsXh,11Jg^!8Br1 C)X^_vH hc|bt^poC!dH]c֔PHvH~qr %4yn3:K-h,g0b)d36R<2 ӬfN~B*xѤJKAК(,sZfdgiTKBղ}T9X 324⚝vu=h9!qᙁ!UG,ӾD=PAtL?`zSVt{`|6iʼByowVw߆[0x[d"8ma|Bu]u[@ 2Bs|P%É ;BoR$+Wxa7w.Pja8'}bpRUL͗zE!8a$3}6|Vvf$x*,./؂33S1X0@ C ؀q}YWgj"muzIHl@{l7у>M}@wSnWmdGwo' &]_noxNȁqz{v{u~yQzx~BNx\jqwhP=GztH_`Khc03uJ؉hB(tpgdrp8~$7br|XN؊(i(Xp׌AX3 8U؈`hXq4XЈ]J!u'XVxmBrhWx츆^9ۖ0Xcpk.rHH(X0dp dCRrizHwiWeXODFK`6.9h؎2l#4rMؓ DBxDiFIn(4綔/bՔZ3Gb(Ty?6~Ggph6}d茣PZ< 2Nu)WmH5g]9Ojai{G _ɘ懐xyy "h2x0ɚٚfbK |8byaW @I"tkiٙ@X٘&FਝƝٍ E+i9九虒؋ўCv9׉yş6NyC92h>:qw :LI|*,e5-z$:>Y:":֠<>ɍ(IlP 38G{s(eY JJչ6*Oj5x{RE*pz|.Zs F:~@ШLW4ڧZЦة[Mꨗ: *&J lyAkŒo*ȗ(E4Dz3ny6sw9:7(aa!zfH.1vP}x4Gbw~;*{ )bIkx᭕!&A Pqʯ; Q7kGHQ&'[(*۲)lc\.˲fG6F)4\nJj_ N;&+t{ hc;xQٴ4@T#kG\AڰzP(tV8F@YȠY s[@ e[BGɥW)h-p¦3D^QBT!fgWٹ2`~q[| jPr{Zkxy10iSț/;++qc4ћ,c[-p4 k ˺[8龞ʽ˖=:[LrbIBk[Iuc,,w))<,Y\ *\e(oI8:lyivx;`SH&ÅL # Ť@FKqk<<jx[Z[3jRiyo<R +u À,\<܇Bw_Džhq_Ly{4//kPɓ( ɦ|:4w=A߂釯nxDܷ Ȳ~~8[2 0$[T[AO؜͍Zl7E&x3P$4 ,f\l62fHiB4 9I۳BŕЋ:d!|g,Ĩ{ Љ̘7þ< mX}}EYΰ1Ү)< @e: RӄqWx_4JTV-ǧlE yա0e'ѣLYiBĊa^"tMoF6q[ͦ- ֈ؅ҀMƟ ֈMLy-ؽjҔȇ[C-ڪ٬֮}r=١ՈٯMґHZٹà\ۂ]٠֩iԨۯ- <Ձܱי`_ݗ r*+o-Mڷ=\m<0lTA ުr݃-ڄ>fެ F+3LOAr;Sm0A$sakԄmmjcMˆ|CN%'m2y;h^9X~|be^~;MBFwmd}e,en=߭] Mv 9J܅nm]פ穮PWNM䢾ثݽڿo^cʍكͲ.ͮ3S]ڝ=e띎ڴ >>n꭮:^rn޸N~N>n gtݖN.Ά;]׎% }-$^%/Sd+.0/1/4?*q9PMd>&C_EoHI߶D/FOKU ;MON֗>&SOJ/eOfZ/j/hkS{ac?^luלWn/rsomXO.͇c-vxO??ޖ/ h/it(ʫ?4@^,UZ~blJcJTN3VOcaHBʶ^-D?!&LFn*=ʢnO r :B _;:3@WH4dȀXð)'0Wə E ։u` @+ {+[J+l:9 v,I ||<}K\lmumZ͎:SP!o_i'0})u> @B(HA4tX'A|._.i^9'P'>'$Pb*SIG:0Bԫ@*T"r,״M%=[ܢ$:mRv[kX]8ŵuxW6%#L91yWTǬUI;Ih;mzjԡIl;p5=ww˛?>ۻ?ۿ?`H`` .`>aNHa^ana~b"Hb&b*b.c2Hc6ވc:c>dBIdFdJ.SnZdMD KfHOf#fPf|x1ZP%(2Yśՠ\2pC@.ڒsbTgKM68PĢ::1dFum2d" !f9Ko2$61'ka.ƒ@U| (SV' +3Cc(pkR߂.j~41n^jL.VlH`&BVjn4\LjcZ'j%s@?CmqOY. ys  k14 U ɲ `&Gl@M=ď 1ҐdhdDR% M04 vo[3ݙ!|:h7GgU=wRK<ЙrE L[^FKsEA83XͣG9z`y5>=ř*iRMڥ)fJ;:[ՒQMj@^#a>Yc;˫Spx^Xz5k@f_v QRT-}cn]G}vvZaT^~qaiRY#'utnK!)oK_4uLa5k{{u&"1eVY!}X yDp,q"?`z b%ixӠq,HRn Y(/x7G+KNs '*r"(& K)Q%0pnĈn Kx/ kx? x$.Ox,n_ x4ox;J;Rate_14Ngp_ThrCut/g14Np_rate_ThrCut.C0000664000372200007340000001572712640051472017455 0ustar suleimanitsops// Double_t Schiff(Double_t Z, Double_t T, Double_t k){ // L. I. Schiff, Physical Review 83, 252 (1951). Double_t me = 0.511; // Electron mass (MeV) Double_t E = T + me; // Electron energy Double_t Ep = E-k; // Scattered electron energy (MeV) Double_t alpha = 1.0/137.036; // Fine-structure Constant Double_t re = 2.818e-13; // classical electron radius (cm) Double_t C = 111.0; Double_t M = 1.0 / ( (me*k/(2.0*E*Ep))**2 + (Z**1/3/C)**2 ); Double_t b = 2.0*E*Ep*Z**1/2 / (C*me*k); Double_t atan = TMath::ATan(b); // // Double_t schiff = 2.0*Z**2*re**2*alpha/k * ( ((E**2+Ep**2)/E**2 - 2.0*Ep/(3.0*E)) * (TMath::Log(M) + 1.0 - 2.0/b*atan) + Ep/E * ( 2.0/b**2*TMath::Log(1.0+b**2) + 4.0*(2.0-b**2)/(3.0*b**3)*atan - 8.0/(3.0*b**2)+2.0/9.0)); // cm2/MeV // cout << "E = " << E << " k = " << k << " Schiff = " << schiff <T) return 0.0; Int_t w = 100; Double_t Z = 29.0; // Cu Double_t ng = 0.0; Double_t DE = Delta/w; // for (Int_t i=0; i p + 13C // Rate versus electron beam kinetic energy gROOT->Reset(); // N14 Cross Section (nb) gROOT->ProcessLine(".L csN14_interpolation.C++"); // cout << " Cross Section (nb) = " << csN14_interpolation(8.0) << endl; //// #define MAX_nentries 950; TString pstit; // define style here // general parameters gStyle->SetOptDate(0); gStyle->SetOptTitle(0); gStyle->SetStatColor(10); gStyle->SetStatH(0.2); gStyle->SetStatW(0.3); gStyle->SetOptStat(0); // canvas parameters gStyle->SetFrameBorderMode(0); gStyle->SetFrameBorderSize(0); gStyle->SetFrameFillColor(0); gStyle->SetCanvasColor(0); // // pads parameters // gStyle->SetPadColor(39); gStyle->SetPadColor(0); gStyle->SetPadBorderMode(0); gStyle->SetPadBorderSize(0); gStyle->SetPadBottomMargin(0.18); gStyle->SetPadRightMargin(0.10); gStyle->SetPadLeftMargin(0.14); gStyle->SetLabelSize(0.04,"x"); gStyle->SetLabelSize(0.04,"y"); gStyle->SetTitleXSize(0.08); gStyle->SetPaperSize(10,12); gStyle->SetOptLogx(0); gStyle->SetOptLogy(1); // gStyle->SetNdivisions(505); gStyle->SetTitleYOffset(0.8); gStyle->SetTitleYSize(0.10); gROOT->ForceStyle(); //// Schiff Spectrum Double_t Nv = 6.02214179e+23; // Avogadro constant (/mol) Double_t L = 3.0; // cm, Chamber length Double_t rho = 0.846; // gm/cm^3, N2O density Double_t NMol = L*rho/44.0; // Number of N2O moles (/cm^2) Double_t NO = Nv*NMol; // Number of Oxygen nuclei (/cm^2) Double_t NN14 = 2.0*Nv*NMol; // Number of N nuclei (/cm^2) // Double_t NO16 = NO * 0.9976 / 1.0; // Number of Oxygen 16 nuclei Double_t Depletion = 1.0; // Depletion of Oxygen isotopes Double_t NO17 = NO * 0.00038 / Depletion; // Number of Oxygen 17 nuclei Double_t NO18 = NO * 0.00205 / Depletion; // Number of Oxygen 18 nuclei // Double_t BremSolidAngle = 6.0; // To find gammas hitting glass cell. // Factor to reduce Schiff (with 0.02 mm) to GEANT4 simulation (with 6 mm) // Double_t RadThickness = 0.002; // cm, 0.02 mm Radiator Double_t rhoCu = 8.94; // g/cm3 Cu Double_t MCu = 63.546; // g/mol Cu (A) Double_t Z = 29.0; // Cu Double_t fCu = rhoCu*RadThickness*Nv/MCu; // Number of Cu atoms/cm^2 Double_t Ib = 1.0; // Beam current (uA) Double_t e = 1.602e-19; // C Double_t Ne = 1.0e-6/e; // Number of electrons per 1 uC Double_t Q14Np = 7.551; // Q-Value of 14N(g,p)13C Double_t Delta = 0.001; // 0.01 MeV bins Double_t Neutron_KE[MAX_nentries]; Double_t weight[MAX_nentries]; Double_t T0[MAX_nentries]; // Electron beam kinetic energy (MeV) Double_t err1[MAX_nentries] = 0.0; Double_t cs14Ngp; Double_t rate14Ngp[MAX_nentries] = 0.0; Double_t ThrCut = 50.0; // Threshold energy cut of ion 13C in keV // if (!gROOT->GetClass("TGenPhaseSpace")) gSystem.Load("libPhysics"); Double_t M14N = 14.0 * 931.494061 + 2.8634; // MeV, mass of 14N nucleus TLorentzVector target(0.0, 0.0, 0.0, M14N); Double_t ProtonMass = 1.0*931.494061 + 7.289; // 938.272046 + me Double_t C13Mass = 13.0*931.494061 + 3.1250; Double_t masses[2] = {ProtonMass, C13Mass} ; // //// Generate Phase Space for (Int_t i=0; iSetGrid(); grO17gn->SetMinimum(1.0e-04); grO17gn->SetMaximum(1.0e+04); grO17gn->GetXaxis()->SetLimits(6.0,8.5); grO17gn->GetXaxis()->CenterTitle(); grO17gn->GetXaxis()->SetTitleOffset(1.2); grO17gn->GetXaxis()->SetTitleSize(0.06); TString tit = "Electron Kinetic Energy (MeV)"; grO17gn->GetXaxis()->SetTitle(tit.Data()); grO17gn->GetYaxis()->CenterTitle(); grO17gn->GetYaxis()->SetTitleOffset(1.0); grO17gn->GetYaxis()->SetTitleSize(0.06); TString tit ="Rate (Hz), 1 #muA"; grO17gn->GetYaxis()->SetTitle(tit.Data()); // grO17gn->SetLineColor(kCyan); grO17gn->SetLineWidth(2.0); grO17gn->SetLineStyle(1); grO17gn->Draw("A C"); // leg = new TLegend(0.20,0.75,0.45,0.85); leg->SetFillColor(0); leg->SetShadowColor(0); leg->AddEntry(grO17gn," ^{14}N(#gamma,p)^{13}C ","L"); leg->Draw(); // leg->SetHeader("I_{b} = 1 #muA"); leg->Draw(); // g1->cd(); pstit = "g14Np_rate_ThrCut"; pstit += ".gif"; g1->Print(pstit.Data()); //// } Rate_14Ngp_ThrCut/csN14_interpolation_C.d0000664000372200007340000000773612640051476020427 0ustar suleimanitsops # DO NOT DELETE ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/Polynomial.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/ParamFunction.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/IParamFunction.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/IFunction.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/IFunctionfwd.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/IParamFunctionfwd.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/Util.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/Interpolator.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Math/InterpolationTypes.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TCanvas.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TPad.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TVirtualPad.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TAttPad.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Rtypes.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/RConfig.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/RVersion.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/DllImport.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Rtypeinfo.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/snprintf.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/strlcpy.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TGenericClassInfo.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TSchemaHelper.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TVirtualX.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TNamed.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TObject.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TStorage.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TVersionCheck.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Riosfwd.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TBuffer.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TString.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TMathBase.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TAttLine.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TAttFill.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TAttText.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TAttMarker.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/GuiTypes.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/Buttons.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TQObject.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TAttCanvas.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TCanvasImp.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TGraph.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TVectorFfwd.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TVectorDfwd.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/TFitResultPtr.h ./csN14_interpolation_C.so: /usr/csite/pubtools/root/5.34.11/include/root/cintdictversion.h /usr/csite/pubtools/root/5.34.11/include/root/RVersion.h csN14_interpolation_C__ROOTBUILDVERSION= 5.34/11