1 / 6

Emittance Calculation in G4MICE – Status Chris Rogers, Imperial College February 4, 2004

Emittance Calculation in G4MICE – Status Chris Rogers, Imperial College February 4, 2004. Class CovMatrix. Class to calculate 6D covariance matrix and various emittances.

maxine-castillo
Download Presentation

Emittance Calculation in G4MICE – Status Chris Rogers, Imperial College February 4, 2004

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. Emittance Calculation in G4MICE – StatusChris Rogers, Imperial CollegeFebruary 4, 2004

  2. Class CovMatrix Class to calculate 6D covariance matrix and various emittances. All calculations follow Holzer's emittance definitions using the covariance matrix c.f. Holzer 03 - "Figure of merit for muon cooling - an algorithm for particle counting in coupled phase planes". I haven't put it on CVS but will do so in next ~couple of days. Word of warning - don't trust the z plane! At the moment Recon isn't producing sensible data for z. WIP.

  3. Public Methods AddEvent(Event *) Adds Event to the covariance matrix getCovMatrix() Calculates sigma_ij and returns the Covariance Matrix. getNDEmittance(axis, dimension) Calculates the N-Dimensional emittance. Define an axis 0=x,1=y,2=z if appropriate to the dimension (2, 4 or 6). getNormNDEmittance(axis, dimension) Returns ND Emittance/<pz> getNDSPEmittance(axis, dimension, * Event, NormNDEmittance, * inverseCovarianceMatrix) Returns the single particle emittance for Event. Dimension/Axis as above. Returns -1000 for bad events.

  4. Output Output looks like this: MC Cov Matrices In 1266 1.933 53.3 -3.034 634.5 -4.278 1.933 0.05114 2.542 0.00796 18.17 0.2395 53.3 2.542 1149 2.125 1397 -1.634 -3.034 0.00796 2.125 0.01769 3.802 0.1253 634.5 18.17 1397 3.802 3.437e+05 602 -4.278 0.2395 -1.634 0.1253 602 6.738 <snip> MC Emittances 6D 4Dxy 2Dx 2Dy 2Dz In 9.232e+04 83.47 24.54 12.49 4391 Out 5.85e+07 8865 184.4 270.5 2.721e+04 <snip> Single Particle Emittances MC In MC Out Recon In Recon Out 0.2268 16.08 0.02243 1.011 1.711 118.1 0.2686 2.765 etc etc

  5. Output Cont Use mice2root to draw 2D phase space and superimpose ellipses over them from Analysis to get things like: Looks sensible (Downstream MC data runmode normal 10,000 events)

  6. Summary Analysis package now calculates 2D, 4D, 6D multi-particle emittances. Also calculates 2D, 4D, 6D single particle Emittances. Results look sensible but no proper testing. Next: Look at other ways to calculate emittance Calculate other stuff from these figures – figures of merit etc Rigorous testing? Documentation?

More Related