ATLFAST CORRECTORS VERSION 15.6.1.2 FOR 10TEV

Efficiencies and Contaminations

I have produced a revised set of efficiencies and contaminations for e, gamma and jets using athena version 15.6.1.2. Here I will present the new efficiencies and investigate the cause of efficiencies and contaminations changes. I will start by giving details on how to produce the efficiencies and contaminations.

Producing Efficiencies and Contaminations

Checking Out And Setting Up AnalysisSkeleton.py

Before we can produce the efficiencies and contaminations we must check out the effmisid package. Once you have set up CMT and CVS we can checkout effmisid package (consult logfile ATLFAST-EFF-CON.txt for a quick guide on how to do this).

When checking out effmisid remember to make a new directory to check effmisid into.

Once effmisid is checked out you will find the following files.

  • Analysis2.cxx
  • Analysis2.h
  • AnalysisSkeleton.cxx
  • AnalysisSkeleton.h
  • mytry.C
  • mytry_eta.C
  • AnalysisSkeleton_jobOptions.py

I will describe briefly the purpose of the files. Analysis2.cxx and Analysis2.h create histograms of the ntuples from the truth and reconstructed containers. mytry.C and mytry_eta.C plot the efficiencies and contaminations from the histograms. mytry.C and mytry_eta.C also produces the .dat files that are used in AtlfastC.

AnalysisSkeleton.cxx and AnalysisSkeleton.h are responsible for making the ntuple from a AOD dataset (e.g. Z->ee). You must copy these two files to src and UserAnalysis directories in the athena PhysicsAnalysis package you set up. Once These files are copied, compile CMT by using

CMT bro gmake

. We are now ready to run athena with our new AnalysisSkeleton. First, we must copy AnalysisSkeleton_jobOptions.py to our run directory.

  • IMPORTANT- You will have to included a bug fix in the AnalysisSkeleton_jobOptions.py, as the container name for Tower jets has change in version 15 from version 14. The line of code that you have to insert is shown below.

import JetRec.ParticleJetCompatibility

We can now run over a sample and start to produce efficiencies. I have included a logfile below that outlines the steps for producing efficiencies.

Electron Efficiencies and Contaminations using a Z->ee dataset

The efficiencies presented here are produced using version 15.6.1.2 of athena. I will compare these efficiencies with version 14.5.0 efficiencies produced by Catherine Wright and version 15 efficiencies from the e/gamma group.

Loose Electrons

Fig 1.1 Fig 1.2 Fig 1.3 Fig 1.4 Fig1.5

Fig 1.6 Fig1.7 Fig1.8

Medium Electrons

Fig2.1 Fig2.2 Fig2.3 Fig2.4 Fig2.5

Fig2.6 Fig2.7 Fig2.8

Tight Electrons

Fig3.1 Fig3.2 Fig3.3 Fig3.4 Fig3.5

Fig3.6 Fig3.7 Fig3.8

Photon Efficiencies and Contaminations using a H->gamgam dataset

Timing Study

A timing Study for Atlfast and AtlfastC was carried out in Version 14.5.0 and Version 15.6.1.2. The samples used in this study where ttbar and Z->ee EVGEN samples. The EVGEN datasets where then put into POOLtoAtlfasttoAOD.py job option that produces an Atlfast AOD. With the AOD we create a corrected AOD using AtlfastC_AODtoAOD.py job option, resulting in an AOD with applied efficiencies to replicate the full simulation AODs.

It must also be noted that we can use the PythiatoAtlfasttoAOD.py job option to produce an Atlfast AOD.

Running POOLtoAtlfasttoAOD.py

I will now discuss the steps that where taken to produce an Atlfast AOD (I have included a Logfile below).

Changes to POOLtoAtlfasttoAOD.py

  • Remember only EVGEN datasets can be read into POOLtoAtlfasttoAOD.py.

  • Input dataset to Atlfast_ReadPOOLFile.py, which is located in POOLtoAtlfasttoAOD.py (need to get_file Atlfast_ReadPOOLFile.py)

  • Add ChronoStat service at the end of POOLtoAtlfasttoAOD.py (This outputs the running time of Atlfast and all the Algorithms used).

# Use auditors
theApp.Dlls += [ "GaudiAud" ]
#
# write out a summary of the time spent
# write out a short message upon entering or leaving each algorithm
#
theAuditorSvc.Auditors  += [ "NameAuditor" ]
ChronoStatSvc = Service ( "ChronoStatSvc")
#ChronoStatSvc.NumberOfSkippedEventsForMemStat = 50
theAuditorSvc.Auditors  += [ "MemStatAuditor" ]
MemStatAuditor = theAuditorSvc.auditor( "MemStatAuditor" )
theApp.AuditAlgorithms=True 
MemStatAuditor.OutputLevel = INFO
#theApp.AuditServices=True

Running AtlfastC_AODtoAOD.py

I will now discuss the steps that where taken to produce an AtlfastC AOD (I have included a Logfile below).

Changes to AtlfastC_AODtoAOD.py

  • Edit the input of AtlfastC_AODtoAOD.py to the output of POOLtoAtlfasttoAOD.py.

Timing Table

Results of the Timing Study

Version   |  Sample   |  ATLFAST  |  ATLFASTC  | ATLFASTC(All ALG EFF ON)
------------------------------------------------------------------------
14.5.0        Z->ee
15.6.1.2      Z->ee
         
14.5.0        ttbar
15.6.1.2      ttbar

Timing Samples

Information on the EVGEN and AOD datasets.

AMI Z->ee

AMI ttbar

Location of EVGEN datasets

  • Z->ee

/data/atlas07/mwright/V14-Zee-EVGEN-MIKE/

  • ttbar

/data/atlas04/mwright/ATLFAST/ttSamples/

Location of Atlfast AODS

  • Z->ee

/data/atlas07/mwright/ATLFAST/Timings/Zee/AtlFast/

  • ttbar

/data/atlas04/mwright/ATLFAST/ttSamples/ttAtlfastOutput

Running AtlFast And AtlFastC Logfiles

-- MichaelWright - 2010-03-03

Topic attachments
I Attachment History Action Size DateSorted ascending Who Comment
Unknown file formatC ALLOutputCanvasEtaALL_electron_electron.C r1 manage 11.2 K 2010-03-08 - 16:47 MichaelWright LooseElecEffEta
JPEGjpg LooseElecElec.jpg r1 manage 17.9 K 2010-03-08 - 16:35 MichaelWright LooseElecEff
JPEGjpg LooseElecElecEta.jpg r1 manage 15.5 K 2010-03-08 - 16:48 MichaelWright LooseElecEffEta
JPEGjpg LooseElecJet.jpg r1 manage 16.0 K 2010-03-08 - 19:22 MichaelWright LooseElecJet
JPEGjpg LooseElecJetEta.jpg r1 manage 12.8 K 2010-03-08 - 19:24 MichaelWright LooseElecJetEta
JPEGjpg LooseElecPhoEta.jpg r1 manage 14.2 K 2010-03-08 - 18:52 MichaelWright LooseElecPhoEta
JPEGjpg LooseElecPhoton.jpg r1 manage 14.6 K 2010-03-08 - 18:24 MichaelWright LooseElecPho
JPEGjpg MedElecElec.jpg r1 manage 18.4 K 2010-03-08 - 16:50 MichaelWright MedElecEff
JPEGjpg MedElecElecEta.jpg r1 manage 18.4 K 2010-03-08 - 16:50 MichaelWright MedElecEffEta
JPEGjpg LooseElecMuidMuon.jpg r1 manage 15.7 K 2010-03-09 - 09:01 MichaelWright LooseElecMuidMuon
JPEGjpg LooseElecMuidMuonEta.jpg r1 manage 12.4 K 2010-03-09 - 09:01 MichaelWright LooseElecMuidMuonEta
JPEGjpg MediumElecJet.jpg r1 manage 14.3 K 2010-03-09 - 10:19 MichaelWright MediumElecJet
JPEGjpg MediumElecJetEta.jpg r1 manage 13.9 K 2010-03-09 - 10:19 MichaelWright MediumElecJetEta
JPEGjpg MediumElecMuidMuon.jpg r2 r1 manage 12.0 K 2010-03-09 - 11:21 MichaelWright MediumElecMuidMuon
JPEGjpg MediumElecMuidMuonEta.jpg r1 manage 12.2 K 2010-03-09 - 11:22 MichaelWright MediumElecMuidMuonEta
JPEGjpg MediumElecPho.jpg r1 manage 13.6 K 2010-03-09 - 09:29 MichaelWright MediumElecPho
JPEGjpg MediumElecPhoEta.jpg r1 manage 12.9 K 2010-03-09 - 09:30 MichaelWright MediumElecPhoEta
JPEGjpg TigElecElec.jpg r1 manage 19.1 K 2010-03-09 - 11:55 MichaelWright TightElecPho
JPEGjpg TigElecElecEta.jpg r2 r1 manage 16.0 K 2010-03-09 - 11:55 MichaelWright TightElecPhoEta
JPEGjpg TigElecJet.jpg r1 manage 17.4 K 2010-03-09 - 12:40 MichaelWright TightElecJet
JPEGjpg TigElecMuidMuon.jpg r1 manage 10.7 K 2010-03-09 - 12:59 MichaelWright TightElecMuidMuon
JPEGjpg TigElecPho.jpg r1 manage 16.4 K 2010-03-09 - 11:58 MichaelWright TightElecPho
JPEGjpg TightElecJetEta.jpg r1 manage 14.3 K 2010-03-09 - 12:40 MichaelWright TightElecJetEta
JPEGjpg TightElecMuidMuonEta.jpg r1 manage 12.1 K 2010-03-09 - 13:00 MichaelWright TightElecMuidMuonEta
JPEGjpg TightElecPhoEta.jpg r1 manage 14.5 K 2010-03-09 - 11:58 MichaelWright TightElecPhoEta
Edit | Attach | Print version | History: r25 | r11 < r10 < r9 < r8 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r9 - 2010-03-17 - MichaelWright
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback