(r10) AtlfastCorrectors10TevV15 < ATLAS

ATLAS Web>AtlfastCorrectors10TevV15 (revision 10) (raw view)~~EditAttach~~

---+ ATLFAST CORRECTORS VERSION 15.6.1.2 FOR 10TEV

%TOC%

---++ %MAKETEXT{"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.

---++++ %MAKETEXT{"Producing Efficiencies and Contaminations"}%
---+++++ %MAKETEXT{"Checking Out And Setting Up <nop>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
   * <nop>AnalysisSkeleton.cxx
   * <nop>AnalysisSkeleton.h
   * mytry.C
   * mytry_eta.C
   * <nop>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 <nop>AtlfastC.

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

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

   * IMPORTANT- You will have to included a bug fix in the <nop>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.

<pre>import JetRec.ParticleJetCompatibility</pre>

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.

   * [[http://ppewww.physics.gla.ac.uk/~mwright/ATLFAST-EFF-CON.html][Logfile for producing Atlfast efficiencies and contaminations]]

---++++ %MAKETEXT{"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.


---+++++%MAKETEXT{"<strong>Loose Electrons</strong>"}%

<!-- img src="%ATTACHURLPATH%/LooseElecElec.jpg" alt="LooseElecElec.jpg" width='596' height='572' / -->

Fig 1.1 <A HREF="%ATTACHURLPATH%/LooseElecElec.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecElec.jpg" height="200" /></A> Fig 1.2<A HREF="%ATTACHURLPATH%/LooseElecElecEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecElecEta.jpg" height="200" /></A> Fig 1.3 <A HREF="%ATTACHURLPATH%/LooseElecPhoton.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecPhoton.jpg" height="200" /></A> Fig 1.4 <A HREF="%ATTACHURLPATH%/LooseElecPhoEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecPhoEta.jpg" height="200" /></A> Fig1.5 <A HREF="%ATTACHURLPATH%/LooseElecJet.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecJet.jpg" height="200" /></A>

Fig 1.6 <A HREF="%ATTACHURLPATH%/LooseElecJetEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecJetEta.jpg" height="200" /></A> Fig1.7 <A HREF="%ATTACHURLPATH%/LooseElecMuidMuon.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecMuidMuon.jpg" height="200" /></A> Fig1.8 <A HREF="%ATTACHURLPATH%/LooseElecMuidMuonEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/LooseElecMuidMuonEta.jpg" height="200" /></A>

---+++++%MAKETEXT{"<strong>Medium Electrons</strong>"}%

Fig2.1<A HREF="%ATTACHURLPATH%/MedElecElec.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MedElecElec.jpg" height="200" /></A> Fig2.2<A HREF="%ATTACHURLPATH%/MedElecElecEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MedElecElecEta.jpg" height="200" /></A> Fig2.3 <A HREF="%ATTACHURLPATH%/MediumElecPho.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MediumElecPho.jpg" height="200" /></A> Fig2.4 <A HREF="%ATTACHURLPATH%/MediumElecPhoEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MediumElecPhoEta.jpg" height="200" /></A> Fig2.5 <A HREF="%ATTACHURLPATH%/MediumElecJet.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MediumElecJet.jpg" height="200" /></A>

Fig2.6<A HREF="%ATTACHURLPATH%/MediumElecJetEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MediumElecJetEta.jpg" height="200" /></A> Fig2.7 <A HREF="%ATTACHURLPATH%/MediumElecMuidMuon.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MediumElecMuidMuon.jpg" height="200" /></A> Fig2.8 <A HREF="%ATTACHURLPATH%/MediumElecMuidMuonEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/MediumElecMuidMuonEta.jpg" height="200" /></A>

---+++++%MAKETEXT{"<strong>Tight Electrons</strong>"}%

Fig3.1<A HREF="%ATTACHURLPATH%/TigElecElec.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TigElecElec.jpg" height="200" /></A> Fig3.2 <A HREF="%ATTACHURLPATH%/TigElecElecEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TigElecElecEta.jpg" height="200" /></A> Fig3.3<A HREF="%ATTACHURLPATH%/TigElecPho.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TigElecPho.jpg" height="200" /></A> Fig3.4<A HREF="%ATTACHURLPATH%/TightElecPhoEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TightElecPhoEta.jpg" height="200" /></A> Fig3.5 <A HREF="%ATTACHURLPATH%/TigElecJet.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TigElecJet.jpg" height="200" /></A>

Fig3.6 <A HREF="%ATTACHURLPATH%/TightElecJetEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TightElecJetEta.jpg" height="200" /></A> Fig3.7 <A HREF="%ATTACHURLPATH%/TigElecMuidMuon.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TigElecMuidMuon.jpg" height="200" /></A> Fig3.8 <A HREF="%ATTACHURLPATH%/TightElecMuidMuonEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/TightElecMuidMuonEta.jpg" height="200" /></A>


---++++ %MAKETEXT{"Photon Efficiencies and Contaminations using a H->gamgam dataset"}%

---+++++%MAKETEXT{"<strong>Photon Efficiency</strong>"}%
Fig4.1 <A HREF="%ATTACHURLPATH%/PhoPho.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoPho.jpg" height="200" /> Fig4.2 <A HREF="%ATTACHURLPATH%/PhoPhoEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoPhoEta.jpg" height="200" />


---+++++%MAKETEXT{"<strong>Photon to Electron Contamination</strong>"}%

Fig5.1 <A HREF="%ATTACHURLPATH%/PhoLooseElec.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoLooseElec.jpg" height="200" /> Fig5.2 <A HREF="%ATTACHURLPATH%/PhoLooseElecEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoLooseElecEta.jpg" height="200" /> Fig5.3 <A HREF="%ATTACHURLPATH%/PhoMediumElec.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoMediumElec.jpg" height="200" /> Fig5.3 <A HREF="%ATTACHURLPATH%/PhoMediumElecEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoMediumElecEta.jpg" height="200" /> Fig5.4 <A HREF="%ATTACHURLPATH%/PhoTightElec.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoTightElec.jpg" height="200" />

Fig5.5 <A HREF="%ATTACHURLPATH%/PhoTightElecEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoTightElecEta.jpg" height="200" />


---+++++%MAKETEXT{"<strong>Photon to C4 Jet Contamination</strong>"}%

Fig6.1 <A HREF="%ATTACHURLPATH%/PhoJet.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoJet.jpg" height="200" /> Fig6.2 <A HREF="%ATTACHURLPATH%/PhoJetEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoJetEta.jpg" height="200" />

---+++++%MAKETEXT{"<strong>Photon to Muid Muon Contamination</strong>"}%

Fig7.1 <A HREF="%ATTACHURLPATH%/PhoMuidMuon.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoMuidMuon.jpg" height="200" /> Fig7.2 <A HREF="%ATTACHURLPATH%/PhoMuidMuonEta.jpg" target="_blank"><img width="200" alt="" src="%ATTACHURLPATH%/PhoMuidMuonEta.jpg" height="200" />




---++ %MAKETEXT{"Timing Study"}%

A timing Study for Atlfast and <nop>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-&gt;ee EVGEN samples. The EVGEN datasets where then put into <nop>POOLtoAtlfasttoAOD.py job option that produces an Atlfast AOD. With the AOD we create a corrected AOD using <nop>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 <nop>PythiatoAtlfasttoAOD.py job option to produce an Atlfast AOD.

---++++ %MAKETEXT{"Running <nop>POOLtoAtlfasttoAOD.py"}%

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

---+++++ %MAKETEXT{"Changes to <nop>POOLtoAtlfasttoAOD.py"}%

   * Remember only EVGEN datasets can be read into <nop>POOLtoAtlfasttoAOD.py.

   * Input dataset to <nop>Atlfast_ReadPOOLFile.py, which is located in <nop>POOLtoAtlfasttoAOD.py (need to get_file <nop>Atlfast_ReadPOOLFile.py)

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

<pre># 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
</pre>

---++++ %MAKETEXT{"Running <nop>AtlfastC_AODtoAOD.py"}%

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

---+++++ %MAKETEXT{"Changes to <nop>AtlfastC_AODtoAOD.py"}%

   * Edit the input of <nop>AtlfastC_AODtoAOD.py to the output of <nop>POOLtoAtlfasttoAOD.py.

---++++ %MAKETEXT{"Timing Table"}%

Results of the Timing Study

<pre>Version   |  Sample   |  ATLFAST  |  ATLFASTC  | ATLFASTC(All ALG EFF ON)
------------------------------------------------------------------------
14.5.0        Z-&gt;ee
15.6.1.2      Z-&gt;ee
         
14.5.0        ttbar
15.6.1.2      ttbar
</pre>

---++++ %MAKETEXT{"Timing Samples"}%

Information on the EVGEN and AOD datasets.

---+++++ %MAKETEXT{"AMI Z->ee "}%

   * [[http://ami.in2p3.fr:8080/AMI/servlet/net.hep.atlas.Database.Bookkeeping.AMI.Servlet.Command?linkId=512][mc08.106050.PythiaZee_1Lepton.evgen.EVNT.e347]]

---+++++ %MAKETEXT{"AMI ttbar"}%

   * [[http://ami.in2p3.fr:8080/AMI/servlet/net.hep.atlas.Database.Bookkeeping.AMI.Servlet.Command?linkId=512][mc08.105206.AcerMCttbarHerwig.evgen.EVNT.e429]]

---+++++ %MAKETEXT{"Location of EVGEN datasets"}%

   * Z-&gt;ee

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

   * ttbar

/data/atlas04/mwright/ATLFAST/ttSamples/

---+++++ %MAKETEXT{"Location of Atlfast AODS"}%

   * Z-&gt;ee

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

   * ttbar

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

---++++ %MAKETEXT{"Running <nop>AtlFast<nop> And <nop>AtlFastC<nop> Logfiles"}%

   * [[http://ppewww.physics.gla.ac.uk/~mwright/Atlfast-PooltoAtlfasttoAOD.html][Running Atlfast with PooltoAtlfasttoAOD.py JO]]
   * [[http://ppewww.physics.gla.ac.uk/~mwright/AtlFastCrunning.html][Running AtlFastC with AtlfastC_AODtoAOD.py JO]]

-- Main.MichaelWright - 2010-03-03

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
c	ALLOutputCanvasEtaALL_electron_electron.C	r1	manage	11.2 K	2010-03-08 - 16:47	MichaelWright	LooseElecEffEta
jpg	LooseElecElec.jpg	r1	manage	17.9 K	2010-03-08 - 16:35	MichaelWright	LooseElecEff
jpg	LooseElecElecEta.jpg	r1	manage	15.5 K	2010-03-08 - 16:48	MichaelWright	LooseElecEffEta
jpg	LooseElecJet.jpg	r1	manage	16.0 K	2010-03-08 - 19:22	MichaelWright	LooseElecJet
jpg	LooseElecJetEta.jpg	r1	manage	12.8 K	2010-03-08 - 19:24	MichaelWright	LooseElecJetEta
jpg	LooseElecMuidMuon.jpg	r1	manage	15.7 K	2010-03-09 - 09:01	MichaelWright	LooseElecMuidMuon
jpg	LooseElecMuidMuonEta.jpg	r1	manage	12.4 K	2010-03-09 - 09:01	MichaelWright	LooseElecMuidMuonEta
jpg	LooseElecPhoEta.jpg	r1	manage	14.2 K	2010-03-08 - 18:52	MichaelWright	LooseElecPhoEta
jpg	LooseElecPhoton.jpg	r1	manage	14.6 K	2010-03-08 - 18:24	MichaelWright	LooseElecPho
jpg	MedElecElec.jpg	r1	manage	18.4 K	2010-03-08 - 16:50	MichaelWright	MedElecEff
jpg	MedElecElecEta.jpg	r1	manage	18.4 K	2010-03-08 - 16:50	MichaelWright	MedElecEffEta
jpg	MediumElecJet.jpg	r1	manage	14.3 K	2010-03-09 - 10:19	MichaelWright	MediumElecJet
jpg	MediumElecJetEta.jpg	r1	manage	13.9 K	2010-03-09 - 10:19	MichaelWright	MediumElecJetEta
jpg	MediumElecMuidMuonEta.jpg	r1	manage	12.2 K	2010-03-09 - 11:22	MichaelWright	MediumElecMuidMuonEta
jpg	MediumElecPho.jpg	r1	manage	13.6 K	2010-03-09 - 09:29	MichaelWright	MediumElecPho
jpg	MediumElecPhoEta.jpg	r1	manage	12.9 K	2010-03-09 - 09:30	MichaelWright	MediumElecPhoEta
jpg	PhoJet.jpg	r1	manage	17.0 K	2010-03-17 - 10:56	MichaelWright
jpg	PhoJetEta.jpg	r1	manage	16.8 K	2010-03-17 - 10:57	MichaelWright
jpg	PhoLooseElec.jpg	r1	manage	19.5 K	2010-03-17 - 10:34	MichaelWright
jpg	PhoMediumElec.jpg	r1	manage	16.7 K	2010-03-17 - 10:44	MichaelWright
jpg	PhoMediumElecEta.jpg	r1	manage	12.9 K	2010-03-17 - 10:44	MichaelWright
jpg	PhoMuidMuon.jpg	r1	manage	15.8 K	2010-03-17 - 10:59	MichaelWright
jpg	PhoMuidMuonEta.jpg	r1	manage	12.4 K	2010-03-17 - 11:00	MichaelWright
jpg	PhoPho.jpg	r1	manage	20.5 K	2010-03-17 - 10:28	MichaelWright
jpg	PhoPhoEta.jpg	r1	manage	16.8 K	2010-03-17 - 10:43	MichaelWright
jpg	PhoTightElec.jpg	r1	manage	16.2 K	2010-03-17 - 10:47	MichaelWright
jpg	PhoTightElecEta.jpg	r1	manage	12.8 K	2010-03-17 - 10:47	MichaelWright
jpg	TigElecElec.jpg	r1	manage	19.1 K	2010-03-09 - 11:55	MichaelWright	TightElecPho
jpg	TigElecJet.jpg	r1	manage	17.4 K	2010-03-09 - 12:40	MichaelWright	TightElecJet
jpg	TigElecMuidMuon.jpg	r1	manage	10.7 K	2010-03-09 - 12:59	MichaelWright	TightElecMuidMuon
jpg	TigElecPho.jpg	r1	manage	16.4 K	2010-03-09 - 11:58	MichaelWright	TightElecPho
jpg	TightElecJetEta.jpg	r1	manage	14.3 K	2010-03-09 - 12:40	MichaelWright	TightElecJetEta
jpg	TightElecMuidMuonEta.jpg	r1	manage	12.1 K	2010-03-09 - 13:00	MichaelWright	TightElecMuidMuonEta
jpg	TightElecPhoEta.jpg	r1	manage	14.5 K	2010-03-09 - 11:58	MichaelWright	TightElecPhoEta
jpg	MediumElecMuidMuon.jpg	r2 r1	manage	12.0 K	2010-03-09 - 11:21	MichaelWright	MediumElecMuidMuon
jpg	PhoLooseElecEta.jpg	r2 r1	manage	16.1 K	2010-03-17 - 10:36	MichaelWright
jpg	TigElecElecEta.jpg	r2 r1	manage	16.0 K	2010-03-09 - 11:55	MichaelWright	TightElecPhoEta

Topic revision: r10 - 2010-03-17 - MichaelWright

ATLAS