1 / 19

Energy Collection on the Space

Energy Collection on the Space. number of events. number of events. Generated P.E. Generated P.E. PhotoElectron Distribution. Linear Fit. number of events. number of events. Reconstructed Energy. Generated Energy. Energy Reconstruction. Generated Energy 4.544MeV

lynna
Télécharger la présentation

Energy Collection on the Space

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Energy Collection on the Space

  2. number of events number of events Generated P.E. Generated P.E. PhotoElectron Distribution Linear Fit number of events number of events Reconstructed Energy Generated Energy Energy Reconstruction Generated Energy 4.544MeV (for confirmation)

  3. Vertex Reconstruction Blue = true Red = reconstructed Overlapped each position # of events σ of dx = 104.6±0.9 σ of dy = 105.5±0.9 σ of dz = 119.5±1.0 Generated_X – reconstructed_X Generated_Y – reconstructed_Y Generated_Z – reconstructed_Z

  4. Reconstructed Radius # of events R_xy space R_xyz space Reconstructed Radius

  5. Surface Drawing # of events absolute Z axis R_xy space

  6. Processing chisqr2.C... ********** ** 1 **SET ERR 1 ********** PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 p0 1.00000e+00 1.00000e-01 -2.00000e+00 2.00000e+00 2 p1 0.00000e+00 1.00000e-01 -2.00000e+00 2.00000e+00 3 p2 0.00000e+00 1.00000e-01 -2.00000e+00 2.00000e+00 4 p3 0.00000e+00 1.00000e-01 -2.00000e+00 2.00000e+00 5 p4 0.00000e+00 1.00000e-01 -2.00000e+00 2.00000e+00 6 p5 0.00000e+00 1.00000e-01 -2.00000e+00 2.00000e+00 ********** ** 2 **MIGRAD 5e+04 ********** FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4. START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-04 FCN=794.106 FROM MIGRAD STATUS=INITIATE 24 CALLS 25 TOTAL EDM= unknown STRATEGY= 1 NO ERROR MATRIX EXT PARAMETER CURRENT GUESS STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 p0 1.00000e+00 1.00000e-01 5.77994e-02 -4.02331e+04 2 p1 0.00000e+00 1.00000e-01 5.00209e-02 -6.13908e+07 3 p2 0.00000e+00 1.00000e-01 5.00209e-02 -9.24939e+10 4 p3 0.00000e+00 1.00000e-01 5.00209e-02 -4.95488e+07 5 p4 0.00000e+00 1.00000e-01 5.00209e-02 -6.84853e+10 6 p5 0.00000e+00 1.00000e-01 5.00209e-02 -6.25011e+10

  7. MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. COVARIANCE MATRIX CALCULATED SUCCESSFULLY FCN=8.64062 FROM MIGRAD STATUS=CONVERGED 228 CALLS 229 TOTAL EDM=2.04309e-08 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 p0 1.00544e+00 8.23636e-03 1.54324e-06 7.35599e-02 2 p1 -1.91320e-05 1.34728e-05 1.30409e-09 4.99075e+01 3 p2 3.79252e-08 6.32446e-09 8.91741e-13 7.15699e+04 4 p3 -1.26356e-05 1.27494e-05 1.25196e-09 2.61855e+01 5 p4 1.47322e-08 5.80395e-09 8.38881e-13 -9.66419e+03 6 p5 2.76868e-09 5.83280e-09 1.22430e-12 4.86168e+04 EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 6 ERR DEF=1 6.784e-05 -7.793e-08 2.352e-11 -7.388e-08 2.132e-11 2.685e-11 -7.793e-08 1.815e-10 -7.559e-14 2.680e-11 3.860e-17 -3.102e-14 2.352e-11 -7.559e-14 4.000e-17 3.933e-17 -3.349e-20 1.459e-22 -7.388e-08 2.680e-11 3.933e-17 1.625e-10 -6.540e-14 -2.949e-14 2.132e-11 3.860e-17 -3.349e-20 -6.540e-14 3.369e-17 3.700e-20 2.685e-11 -3.102e-14 1.459e-22 -2.949e-14 3.700e-20 3.402e-17 PARAMETER CORRELATION COEFFICIENTS NO. GLOBAL 1 2 3 4 5 6 1 0.97691 1.000 -0.702 0.452 -0.704 0.446 0.559 2 0.99181 -0.702 1.000 -0.887 0.156 0.000 -0.395 3 0.98255 0.452 -0.887 1.000 0.000 -0.001 0.000 4 0.99157 -0.704 0.156 0.000 1.000 -0.884 -0.396 5 0.98165 0.446 0.000 -0.001 -0.884 1.000 0.001 6 0.96090 0.559 -0.395 0.000 -0.396 0.001 1.000

  8. Result Reconstructed Ratio(gen/rec) absolute Z axis R_xy space R_xy space • FitFCN = • 1.00544 • +0.523 * pow(10,-5) * rxy[entry] + 0.791 * pow(10,-5) * abs(int(RadialSum[entry].Z())) • (0.264*pow(10,-8) * pow(rxy[entry],2) + 0.679 * pow(10,-8) * pow(abs(int(RadialSum[entry].Z())),2)) • 0.361 * pow(10,-9) * rxy[entry] * abs(int(RadialSum[entry].Z())) Z axis

  9. Supplement

  10. Energy Reconstruction by Space- The ratio between generated and reconstructed on the ZvsR relation gen/rec ratio absolute Z-axis reconstructed radius R verse Z-axis R verse gen/rec Z verse gen/rec

  11. Motivation of The Curve Analysis Quadratic Equation Inverse Quadratic Equation

  12. Variable relation not applied plot gen/rec ratio Used all Reconstructed radius R verse Z R verse E Z verse E Used relation Calculated Energy[array] for calculating RvsEnergy = 0.42 * pow(10,-9) * pow(R[array],2) + 0.392* pow(10,-5) * R[array] + 0.99 Calculated Energy[array] for calculating R,ZvsEnergy R&Z relation -› 0.000116 * pow(R[array],2) + 0.3665 * R[array] + 1200 Z&Energy relation -› 0.563 * pow(10,-9) * pow(Z[array],2) + 0.105* pow(10,-5) * Z[array] + 0.9954

  13. 2.4m Con’t – result by R & Energy(using ) 3.1m Target Applied reslut 2.0m Reconstructed gen/rec ratio Box opt Zoom out reconstructed radius Binning: Plot1,2: X-axis:50,0,3500 Y-axis:50,0.98,1.02 Plot3: X-axis:50,0,3500 Y-axis:25,0.5,1.5

  14. Con’t – result by R,Z & Energy (using ) Applied reslut Reconstructed gen/rec ratio Box opt reconstructed radius Binning: Plot1,2: X-axis:50,0,3500 Y-axis:50,0.995,1.005 Plot3: X-axis:50,0,3500 Y-axis:25,0.5,1.5

  15. Con’t – result by R & Energy(using ) Applied reslut Reconstructed gen/rec ratio Box opt Zoom out reconstructed radius Binning: Plot1,2: X-axis:50,0,3500 Y-axis:50,0.997,1.003 Plot3: X-axis:50,0,3500 Y-axis:25,0.5,1.5

  16. Con’t – result by R,Z & Energy (using ) Applied reslut Reconstructed gen/rec ratio Box opt reconstructed radius Binning: Plot1,2: X-axis:50,0,3500 Y-axis:50,0.994,1.006 Plot3: X-axis:50,0,3500 Y-axis:25,0.5,1.5

  17. Cesium – 10 events Gamma X,Y,Z view

  18. Co60 – 1 event what's the meaning of this Warning? ### Run 0 start. Warning: G4MaterialPropertyVector::GetProperty ==> attempt to Retrieve Property below range 188 G4double PMmin = MPV.front()->GetPhotonMomentum(); 189 G4double minProp = MPV.front()->GetProperty(); 197 if (aPhotonMomentum < PMmin) 198 { 199 G4cout << "\nWarning: G4MaterialPropertyVector::GetProperty"; 200 G4cout << "\n==> attempt to Retrieve Property below range" 201 << G4endl; 202 return minProp; 203 } Now I am tracing this code

More Related