1 / 17

Energy Correction for Dead Cells in ECAL: Parameterization and Simulation Results

The study focuses on developing a formula (F(x, L, E)) to correct the missing energy of dead cells in ECAL based on distance from the center of gravity, layer number, and energy of the shower. An alternative correction method is proposed along with simulation results for dead cell energy deposition. Parameterization analysis reveals correlations between parameters and layer number or energy. Comparison with previous correction methods shows improved results.

muthi
Download Presentation

Energy Correction for Dead Cells in ECAL: Parameterization and Simulation Results

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. ECAL Energy Correction for Dead Cells_August Jie Feng, Corrine Goy, Li Tao LAPP, Annecy 1

  2. Principle ● We assume that the missing energy of the dead cell is a function of the distance of the center of gravity from the dead cell(x), the Layer number(L) and the Energy of the shower(E), as shown in the following expression: ● MissingEnergy= F(x, L, E) ● Our purpose is to find the formula of F(x, L, E) 2

  3. Current Method ● Currently, the correction of the ECAL dead cell energy is: b=1/2(a+c) ● ● Where b is the energy of the dead cell and a and c are the energies of the neighboring cells. ● This formula works well except at the cell of the center of gravity. 3

  4. Current Method If the center of gravity of 1D Cluster happens to be at the dead cell, the missing energy will be about 75% of the layer, as shown in the Figure below. The correction formula above can correct 12% of the energy so about 60% of that will still be missing. 4

  5. Simulation of a dead cell lThe TProfile of deposited energy without dead cells can be found on the left, while that with a dead cell can be found on the right. We simulate that the dead cell is located at the cell[34] of the Layer[7].

  6. lSubtract the profile with a dead cell from that without any dead cell and we get the distribution of the missing energy, as shown in the figure below.

  7. A Direct Correction(Path 2) I study the Radial shower profile in the article by Grindhammer and Peters (http://arxiv.org/pdf/hep-ex/0001020v1.pdf , formula (23)) : 2 C 2 T 2 2 rR  rR     ( ) f r (1 ) p p 2 2 2 ) C R 2 2 2 ) T R ( ( r r I integrate it along y and get:     2 2 2 C 2 2 2 T 2 x 2 x x  y R  x  y R     ( ) f x (1 ) p p dy 2 2 2 2 ) C R 2 2 2 2 ) T R ( ( y y  R R x x R R x x     2 2 ln( 1) ln( 1) x R x R C x C T x T C T 2 2 2 C 2 T R R  2 2   ) 2(1    2 ( p )( ) p    2 C 2 2 C 2 3/2 ) 2 T 2 2 T 2 3/2 ) ( ( R x R R x R 7

  8. A Direct Correction(Path 2) ● I use the above function to fit the Missing Energy and get: ● It seems to be better than the one fitted by the Double Gaussian. 8

  9. Parameterization ● For the function: R R x x R R x x     2 2 ln( 1) ln( 1) x R x R C x C T x T C T 2 2 2 C 2 T R R  2 2     ( ) f x ( ) ( ) N N 1 2    2 C 2 2 C 2 3/2 ) 2 T 2 2 T 2 3/2 ) ( ( R x R R x R core tail ● Where x is the distance between the dead cell the and shower center, which is computed by the cell ratio method, 2p is the normalization of the core, R_C is the RMS of the core, 2(1-p) is the normalization of the tail and R_T is the RMS of the tail. 9

  10. Parameterization ● Example of the Cell_34 and Layer_7: 10

  11. Parameterization ● Idea: ● These 4 parameters is correlated with the Layer number L and the shower energy E. Thus,       ( , ) ( , ) ( , ) ( , ) T R L E ( , 0( ), 1( )...) N L p E p E R L p E p E N L p E p E   N R N R N L E R L E N L E 1 1 1 ( , 0( ), 1( )...) ( , 0( ), 1( )...) ( , 0( ), 1( )...) T R L p E p E c c c 2 2 2 T 11

  12. Parameterization The fitting parameters change while the Layer number changes. For the 100GeV Test Beam electron, we get: 12

  13. Parameterization ● Results:       ( , ) ( , ) ( , ) N L E 0( ) 0( ) 0( ) p ( , 1( ), 2( )) N R N N L E R L E p p E GammaDist L p E p E L E GammaDist L p E p   E 1 1   1( ) p E c c   ( , 1( ), 2( ))   E 2 2 2 ( , ) R L E 0( ) E 1( ) p E 2( ) E R p L p L T T 13

  14. Parameterization I also analyze those parameters with different Energies. 20GeV positrons, 80GeV positrons, 100GeV electrons, 120GeV electrons and 180GeV electrons. The graph below is the parameter N_1(L) for different energies. 14

  15. Parameterization Comparing the 5 curves, we can analyze the relation between the parameter and the Energy. For example, the graph below shows the analysis of p0. Thus, we can roughly get:   0 645.87 4.842 2 1.861  (  ) 142.76 p p p E E 1 2 15

  16. Conclusion As far as I can do, the deposited energy of a dead cell should be a function of distance from the Shower Center x, the Layer number and the energy of the particle E. Thus, R R x x R R x x     2 2 ln( 1) ln( 1) x R x R C x C T x T C T 2 2 2 C 2 T R R  2 2     ( , , ) f x L E ( ) ( ) N N 1 2    2 C 2 2 C 2 3/2 ) 2 T 2 2 T 2 3/2 ) ( ( R x R R x R where       ( , ) 645.87 R L E N L E N N L E 1 1  ( 142.76) ( ,4.842,1.8612) E  GammaDist L L  0.481 0.0412  ( , ) ( 28.211 ( , ) T R L E R N c c , ) 2 2  ( 62.7 8.186 0.8681  87 ) ( ,8.5554,1.737) 0.0357 4  E GammaDist L    2 R L L 16 T

  17. Comparison with the old correction The TProfile with the old correction(ADC[i]=(ADC[i-1]+ADC[i+1])/2) is shown in the left, while that with the new correction is shown in the right. 17

More Related