egammaPID determination

ElIsEMMask is an Athena algorithm that will print the IsEM masks for the Athena release it is run with.

Usage

Set up Athena with your desired version and change to your TestArea, e.g.:

source cmthome/setup.sh -tag=14.5.0
cd $TestArea

The program package is in the SVN: https://svnsrv.desy.de/k5viewvc/atlas/Top/clange/Athena/ElIsEMMask/. Check it out to your TestArea:

svn co https://svnsrv.desy.de/desy/atlas/Top/clange/Athena/ElIsEMMask/trunk ElIsEMMask

If you use Athena release 14, do instead:

svn co https://svnsrv.desy.de/desy/atlas/Top/clange/Athena/ElIsEMMask/tags/ElIsEMMask-14-05-00 ElIsEMMask

Go to the cmt directory of the package and compile:

cd ElIsEMMask/cmt
cmt config
source setup.sh
cmt make

Run the program (no output files will be produced):

athena ../share/ElIsEMMask_jobOptions.py

All information will be printed within the program's INFO output.

Results for different Athena releases

In release 15 additional isolation variables have been introduced with the Iso appendix. To distinguish those more easily, the non-isolated variables a la Medium are now additionally available as MediumNonIso (but they are the same, i.e. Medium = MediumNonIso).

For details on the cuts refer to http://alxr.usatlas.bnl.gov/lxr/source/atlas/Reconstruction/egamma/egammaEvent/egammaEvent/egammaPIDdefs.h

Further details

// new definitions of cuts as in rel>=15.2.0
// @brief Loose electron selection
const unsigned int ElectronLoose = 
CALOMIDDLE_ELECTRON | HADLEAKETA_ELECTRON;
// @brief Medium electron selection
const unsigned int ElectronMedium = 
CALO_ELECTRON | TRACKINGNOBLAYER_ELECTRON | TRACKMATCHDETA_ELECTRON;
// @brief Medium electron selection with isolation requirement
const unsigned int ElectronMediumIso = 
CALO_ELECTRON | TRACKINGNOBLAYER_ELECTRON | TRACKMATCHDETA_ELECTRON | ISOLATION_ELECTRON;
// @brief obsolete - kept not crash clients
const unsigned int ElectronMediumNoIso = ElectronMedium;
// @brief Tight electron selection
const unsigned int ElectronTight = 
CALO_ELECTRON | TRACKING_ELECTRON | TRACKMATCH_ELECTRON | TRT_ELECTRON ;
// @brief obsolete - kept not to crash clients
const unsigned int ElectronTightTRTNoIso = ElectronTight;
// @brief obsolete - kept not to crash clients
const unsigned int ElectronTightNoIsolation = ElectronTight;
// @brief Tight electron selection with isolation requirement
const unsigned int ElectronTightIso = 
CALO_ELECTRON | TRACKING_ELECTRON | TRACKMATCH_ELECTRON | TRT_ELECTRON | ISOLATION_ELECTRON;

more detailed:

// @brief all cuts in calorimeter (except isolation)
const unsigned int CALO_ELECTRON = HADLEAKETA_ELECTRON | CALOSTRIPS_ELECTRON | CALOMIDDLE_ELECTRON ;

and now the highest detail

// @brief cuts of hadronic leakage
const unsigned int HADLEAKETA_ELECTRON = 
0x1 << ClusterEtaRange_Electron        | 
0x1 << ClusterHadronicLeakage_Electron;

// @brief cuts in strips (with ClusterStripsDEmaxs1)
const unsigned int CALOSTRIPS_ELECTRON =
0x1 << ClusterStripsEratio_Electron     |
0x1 << ClusterStripsDeltaEmax2_Electron |
0x1 << ClusterStripsDeltaE_Electron     |
0x1 << ClusterStripsWtot_Electron       |
0x1 << ClusterStripsFracm_Electron      |
0x1 << ClusterStripsWeta1c_Electron     |
0x1 << ClusterStripsDEmaxs1_Electron    ;

//@brief cuts in middle sampling
const unsigned int CALOMIDDLE_ELECTRON =    
0x1 << ClusterMiddleEnergy_Electron     |
0x1 << ClusterMiddleEratio37_Electron   |
0x1 << ClusterMiddleWidth_Electron     ;

// @brief Track quality cuts except b-layer for electrons
const unsigned int TRACKINGNOBLAYER_ELECTRON =     
0x1 << TrackPixel_Electron   |
0x1 << TrackSi_Electron      |
0x1 << TrackA0_Electron;

// @brief Track cluster matching in eta for electrons
const unsigned int TRACKMATCHDETA_ELECTRON =
0x1 << TrackMatchEta_Electron;

In the anti-electron note they require 2 or 3 of the following cuts to fail:

• Had/Em - The ratio of the hadronic calorimeter energy of the cluster to its electromagnetic energy

• Chi^{2 - The Chi}2 comparison of the CES shower profile in the r-z view with the shower profile extracted from test beam electrons

• L_shr - The lateral shower profile, measure of how well the lateral shower development matches that expected from the em shower
• CES Delta X - difference between track and CES in x direction
• CES Delta Z - difference between track and CES in z direction

These cuts are ID cuts and NOT kinematic cuts that might bias the distributions!

detailed masks for rel. 15.5.1

detailed masks for rel. 15.6.9.8

detailed masks for rel. 15.6.12.9.1