The e/gamma Analysis Package.

• Purpose
• Package location
• Mailing list - Meetings
• Package description
• Production ntuples
• Running the package
• Definition of the variables used
• Cust applied
• Results
• Useful notes/remarks

The main purpose of the e/gamma Analysis Package is to provide a simple and robust way for calculating the electrons' (including posistrons) and gammas' selection efficiencies for each of the Trigger Levels of the ATLAS Detector. The electrons can come either from single electron simulated/reconstructed events, or from some interesting physics events (signal), such as:

The implementation is extended to include jets, as a background for rate calculations. The package is based on the code written by M. Diaz-Gomez, E. Moyse, V. Perez-Reale, C. Padilla, A. Wildauer. For more information about the original version, visit the initial web-page of the package.

The package is located at the ATLAS CVS Repository:

offline/Trigger/TrigAnalysis/TrigEgammaAnalysis/

The current, latest version of the package is:

There is also available the former TrigNtEgamma package. Although this package is not maintained anymore, it is useful for the oldest versions of the trigger e/gamma analysis.

The ATLAS mailing list for the package is: atlas-trig-egamma@cern.ch

Subscribers of the list are informed about the regular meetings (normally 1 hour before the PESA meetings) and about the latest modifications and tags of the package.

The agendas and the minutes of the meetings can be found at the CDS Agenda System.

• Structure
  

  The following diagram gives the structure of the Package. Click on it to enlarge.

  




• Steering file

  At the Steering File, which has been introduced and replaced the GUI, the user defines the options for running the package.The user can switch between the different data sets without having to build the code again.

  Reading the steering file is a case sensitive procedure. Therefore the boolean values should be YES or NO . The other options are:

  • Number_Of_Events= {ALL, an integer}.
  • Data_set= {e25_lumi05, e30_lumi10, jets17_lumi02, jets25_lumi10, e30_lumi10kineZ, e25_lumi02kineZ}.
  • Production= {6.0.4, 7.0.2, 8.0.2}.
  • Triggering= {e25, e30, jets}.



• Histograms

  The histograms are now saved in a file, for the different reconstruction algorithms (i.e Egamma_Histograms_ALGORITHM.root), whenever the user asks for producing them. The histograms are saved in different directories into the file, depending on the level they refer to. Then, it is easy for the user to apply any draw options on the histograms or even to compare histograms from different data sets and different reconstruction algorithms.

  For example, in the current version of the code, the following histograms are produced:
  • Level1Calo: The Pt distribution of the generated particles.
  • Level2Calo: The Et distribution after the Level2Calo.
  • Level2IDandIDCalo: The Eta-resolution of the selected track algorithm.
  • EventFilterCalo: The Et distribution after the EventFilterCalo.
  • EventFilterID: The distribution of the transverse impact parameter Ao .
  • EventFilterIDCalo: The distribution of the ratio E/Px at the end of the triggering.

  Tips:

  • To change the name of output histogram file, just edit the appropriate line at the Analysis.C file and build the code.
  • To add a new histogram: (i) Declare the histogram at the header file (ii) Initialize it into the constructor, at the source file (iii) Fill it wherever necessary (iv) Write it into the DrawHistos() method.

The productions refer to single electron (or positron) events with Pt = 25 GeV (Low Luminosity) and with Pt = 30 GeV (High Luminosity). The di-jet events are produced with 17 GeV jets (Low Luminosity) and 25 GeV jets (High Luminosity).

For a further and more detailed description about the production's features, visit the site at RAL

Find below the location of the ntuples.

• Ntuples produced with the 6.0.4 (DC1) ATLAS software release
• Low Luminosity

/castor/cern.ch/atlas/project/dc1/recon/dc1.002026.lumi02.recon.010/dc1.002026.lumi02.recon.010._00001-00010.test.e_minus_25_250.root (1 ntuple with electrons)

/castor/cern.ch/atlas/transfer/emoyse/dc1.002027.lumi02.recon.eplus25.root (1 ntuple with positrons)

• High Luminosity

/castor/cern.ch/atlas/project/dc1/recon/dc1.002021.lumi10.recon.010/dc1.002021.lumi10.recon.010._00001-00040.test.e_minus_30_250.root (1 ntuple with electrons)

/castor/cern.ch/atlas/transfer/emoyse/dc1.002023.lumi10.recon.eplus30.root (1 ntuple with positrons)

• Low Luminosity di-jet events

/castor/cern.ch/atlas/project/dc1/recon/dc1.002000.lumi02.recon.010/

• High Luminosity di-jet events

/castor/cern.ch/atlas/project/dc1/recon/dc1.002001.lumi10.recon.010/



• Ntuples produced with the 7.0.2 ATLAS software release
• Low Luminosity

/castor/cern.ch/user/b/baines/trigprod/rel702++2003-12-12/lumi02/002026/his/ (10 ntuples + 10 ntuples with the KineZ )

• High Luminosity

/castor/cern.ch/user/b/baines/trigprod/rel702++2003-12-12/lumi10/002021/his/ (40 ntuples + 40 ntuples with the KineZ)



• Ntuples produced with the 8.0.2 ATLAS software release
• Low Luminosity

/castor/cern.ch/user/b/baines/trigprod/8.0.2/2004-05-14/lumi02/002026/his/ (10 ntuples)

• At CERN
  1. Check out the package:

     cmt co -r TrigEgammaAnalysis-xx-xx-xx Trigger/TrigAnalysis/TrigEgammaAnalysis

  2. cd Trigger/TrigAnalysis/TrigEgammaAnalysis/TrigEgammaAnalysis-xx-xx-xx/Root
  3. source setup.sh
  4. make
  5. ./main.exe
• Locally
  1. At your home area, create a directory for the package:

     mkdir /home/EgammaAnalysis/

  2. Create a directory where you can put all the production ntuples for the release you want.

     For example: /home/EgammaAnalysis/ProductionNtuples/8.0.2/

  3. Go to the CVS Repository and, having selected the preferable version (tag) of the package, download the tarball. Save it at /home/EgammaAnalysis/ directory.
  4. Unpack the package:

  tar xzvf TrigEgammaAnalysis.tar.gz

  5. cd TrigEgammaAnalysis/Root
  6. Edit the Analysis.C file to give the correct path of the ntuples.

  For the above example, you should write: Path8 (Low or High) = /home/EgammaAnalysis/ProductionNtuples/8.0.2/

  Tip: A search for TString Path will give you the paths for all the releases used.

  7. Finally, compile... :

     make

  8. ...and run:

  ./main.exe

  Tip: As quite often the castor server is down, it is strongly advised to run the package locally. You will only have to transfer the ntuples once. Keep the Analysis.C file so that, for the future package updates, you won't have to fix the paths from the scratch.

• LEVEL 1:

  Variable's Name	Description
  L1em_nroi	Number of EM RoI defined by Level1
  L1em_core
  L1em_emclus	Energy in the central four cells
  L1em_emisol	Energy in the EM layer isolation ring
  L1em_hdisol	Energy deposited in hadronic ring
  L1em_hdcore	Energy deposited in hadronic core
  L1em_eta	Eta of the EM cluster defined by LVL1
  L1em_phi	Phi of the EM cluster defined by LVL1
  Etagen	Generated Eta of the particle
  Phigen	Generated Phi of the particle
  Ptgen	Generated Transverse Momentum of the particle

• LEVEL 2 Calo:

  Variable's Name	Description
  T2canclus	Number of LVL2 e/m clusters
  T2caeta	Weighted eta position of the e/m cluster, calculated in sampling. 2 in a 3x7 cell window
  T2caphi	Weighted phi position of the e/m cluster, calculated in sampling 2 in a 3x7 cell window
  T2caeme	Total energy deposited in LAr EM
  T2cahade	Total Energy deposited in LAr HEC and TileCal
  T2carcore	E3x7/E7x7 in samp. 2 with Enxm is the energy deposited in a window of nxm cells around the refined LVL1 position
  T2caeratio	(E1st - E2nd) / (E1st + E2nd) in sampling 1, where E1st and E2nd are the energies of the two highest maxima

• LEVEL 2 ID and IDCalo:

  Variable's Name	Description
  T2idntracks	Number of Tracks in the Inner Detector
  T2idalgo	Algorithm that created the track (0:Offline IDScan, 1:SiTrack, 2:Online IDScan, 3:TrtLUT, 4:TrtXk)
  T2idpt	The transverse momentum of the reconstructed track
  T2ideta	The eta of the reconstructed track
  T2idd0	R of closest approach to origin
  T2idz0	Z at closest approach to origin
  T2idphi0	Phi at closest approach to origin

• EVENT FILTER Calo:

  Variable's Name	Description
  egnc	Number of reconstructed cluster
  eget	Transverse Energy of the cluster
  egf1	Fraction of energy found in EM Sampling 1
  egeta	Eta of the cluster
  egphi	Phi of the cluster
  ege277	Uncorrected energy in 7x7 cells in EM Sampling 2
  egisem	e-ID flag. 0: electron, >=1: jet
  ege237	Uncorrected energy in 3x7 cells in EM Sampling 2
  egetha1	Et leakage into Sampling 1 of Hadronic Calorimeter
  egemins1	Energy of strip with minimum between maximum 1 and 2
  ege2tsts1	2nd maximum in strips
  egweta2	Corrected width in 3x5 cells in the 2nd Sampling
  egfracs1	Energy in strips outside core [( (E+-7) - (E+-3) ) / (E+-7)]
  egwtots1	Total width in EM Sampling 1 in 20 strips
  egweta1	Corrected width in 3 strips in the 1st Sampling
  et37

• EVENT FILTER ID:

  Variable's Name	Description
  egtrkmatchx	link to best matched track in ntuple (for xKalman tracks)
  Pattern	Bit pattern of SCT and Pixel hits associated
  Nsihits	Number of SCT and Pixel holes
  A0vert	Impact Parameter
  Sharehits	Number of Hits Shared with another reconstructed track

• EVENT FILTER IDCalo:

Variable's Name	Description
egeoverpx	Energy match of cluster with track, E(calo)/p(track)
egdeta1x	Delta eta of track extrapolated to Calo S1
egdeta2x	Delta eta of track extrapolated to Calo S2
Nstrawhits	Number of straws associated with the track
Ntrhits	Number of Transition Radiation hits associated with the track

The cuts which apply at each Trigger Level are of major importance in optimizing the selection efficiencies. It is on the schedule to perform a detail study for optimizing the cuts. Untill then, it is suggested to keep the cuts applied until now. While running the code, one must make sure that the following values HAVE BEEN USED in order for the results to be comparable.

The values-together with the variables (for both luminosities)-can also be found in the /doc directory at the CVS-repository.

The cuts listed below refer to single electron events.