Particle Flow studies for MET Trigger with FTK Tracks

Steps:

1. D3PD production with eFlowRec package
2. Analysis with Marianna's script
3. Resolution plot
4. Turn-on curve

D3PD Production

D3PD data samples used for Trigger/FTK meeting on 23 January 2014:

• FTKConv (FTK refitted tracks) are used as track input to eFlowRec:
  • ZH-> nunubb, mu=60: user.giulini.mc12_14TeV.161849.Pythia8_AU2CTEQ6L1_ZH125_nunubb.recon.DESD_FTK.e2103_s1682_s1691_r4711_r4910.tO.eFlow.FTKconv.2/
  • ZH-> nunubb, mu=80: user.giulini.mc12_14TeV.161849.Pythia8_AU2CTEQ6L1_ZH125_nunubb.recon.DESD_FTK.e2103_s1682_s1691_r4712_r4911.tO.eFlow.FTKconv.2/
  • JZ1, mu=60: user.giulini.mc12_14TeV.147901.Pythia8_AU2CT10_jetjet_JZ1.recon.DESD_FTK.e1712_s1682_s1691_r4711_r4910.tO.eFlow.FTKconv.2/
  • JZ2, mu=60: user.giulini.mc12_14TeV.147901.Pythia8_AU2CT10_jetjet_JZ1.recon.DESD_FTK.e1712_s1682_s1691_r4711_r4910.tO.eFlow.FTKconv.2/
  • JZ3, mu=80: user.giulini.mc12_14TeV.147903.Pythia8_AU2CT10_jetjet_JZ3.recon.DESD_FTK.e1712_s1682_s1691_r4712_r4911.tO.eFlow.FTKconv.2/
• Trk (TrackParticleCandidate) are used as track input to eFlowRec:
  • ZH-> nunubb, mu=60: user.giulini.mc12_14TeV.161849.Pythia8_AU2CTEQ6L1_ZH125_nunubb.recon.DESD_FTK.e2103_s1682_s1691_r4711_r4910.tO.eFlow.Trk.2/
  • ZH-> nunubb, mu=80: user.giulini.mc12_14TeV.161849.Pythia8_AU2CTEQ6L1_ZH125_nunubb.recon.DESD_FTK.e2103_s1682_s1691_r4712_r4911.tO.eFlow.Trk.2/
  • JZ1, mu=60: user.giulini.mc12_14TeV.147901.Pythia8_AU2CT10_jetjet_JZ1.recon.DESD_FTK.e1712_s1682_s1691_r4711_r4910.tO.eFlow.Trk.2/
  • JZ2, mu=60: user.giulini.mc12_14TeV.147902.Pythia8_AU2CT10_jetjet_JZ2.recon.DESD_FTK.e1712_s1682_s1691_r4711_r4910.tO.eFlow.Trk.2/
  • JZ3, mu=80: user.giulini.mc12_14TeV.147903.Pythia8_AU2CT10_jetjet_JZ3.recon.DESD_FTK.e1712_s1682_s1691_r4712_r4911.tO.eFlow.Trk.2/
• NTUP_COMMON:
  • ZH-> nunubb, mu=60:user.giulini.mc12_14TeV.161849.Pythia8_AU2CTEQ6L1_ZH125_nunubb.recon.DESD_FTK.e2103_s1682_s1691_r4711_r4910.3
  • ZH-> nunubb, mu=80: user.giulini.mc12_14TeV.161849.Pythia8_AU2CTEQ6L1_ZH125_nunubb.recon.DESD_FTK.e2103_s1682_s1691_r4712_r4911.3/
  • JZ1, mu=60: user.giulini.mc12_14TeV.147901.Pythia8_AU2CT10_jetjet_JZ1.recon.DESD_FTK.e1712_s1682_s1691_r4711_r4910.3/
  • JZ2, mu=60: user.giulini.mc12_14TeV.147902.Pythia8_AU2CT10_jetjet_JZ2.recon.DESD_FTK.e1712_s1682_s1691_r4711_r4910.3a/
  • JZ3, mu=80: user.giulini.mc12_14TeV.147903.Pythia8_AU2CT10_jetjet_JZ3.recon.DESD_FTK.e1712_s1682_s1691_r4712_r4911.3

Analysis

Here are provided the information how to set up and run Marianna's code.

In folder $TestArea, you setup ROOT:

export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase
source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh
localSetupROOT 5.34.01-x86_64-slc5-gcc4.3

Check out fastjet and install it

cd $TestArea
curl -O http://fastjet.fr/repo/fastjet-3.0.6.tar.gz 
tar zxvf fastjet-3.0.6.tar.gz
cd fastjet-3.0.6/
./configure --prefix=$PWD/../fastjet-install
make 
make check
make install
cd ..

Now installing ApplyJet

svn co svn+ssh://$USER@svn.cern.ch/reps/atlasoff/Reconstruction/Jet/ApplyJetCalibration/tags/ApplyJetCalibration-00-03-19 ApplyJetCalibration
cd ApplyJetCalibration/cmt
make -f Makefile.Standalone
cp EflowJES_1.config $TestArea/ApplyJetCalibration/data/CalibrationFactors/EflowJES_1.config
cp EflowConfig_1.config $TestArea/ApplyJetCalibration/data/CalibrationConfigs/EflowConfig_1.config

Now for TrkMETUtil

# Go to https://svnweb.cern.ch/cern/wsvn/atlas-khoo/khoo/TrkMETUtil/doc/README_TrackMETUtil
# and download the package: atlas-khoo-TrkMETUtil.r152904.tar.gz

tar zxvf atlas-khoo-TrkMETUtil.r152904.tar.gz
cd TrkMETUtil*/cmt
make -f Makefile.Standalone

Now for MissingETUtility

svn co svn+ssh://$USER@svn.cern.ch/reps/atlasoff/Reconstruction/MissingETUtility/tags/MissingETUtility-01-02-05 MissingETUtility
cd MissingETUtility/cmt/
make -f Makefile.Standalone

Now copy in your $TestArea these files:

• EflowConfig_1.config
• EflowJES_1.config
• MET_tO_eFlow_2.C: To run eflow

• MET_tO_eFlow_2.h

• CommonTrkMETUtil.C: to run on NTUP_COMMON

• CommonTrkMETUtil.h

• RunFTKmetStandROOT.C: RunFTKmetStandROOT.C

cp EflowJES_1.config ApplyJetCalibration/data/CalibrationFactors/EflowJES_1.config
cp EflowConfig_1.config ApplyJetCalibration/data/CalibrationConfigs/EflowConfig_1.config
mkdir ftk_met
mv *.C ftk_met/
mv *.h ftk_met/

To run them, you have to create a file with a list of NTUP files.

And then adjust the flags:

doNTUPCOMMON should be 1 if you want to run over the NTUP_COMMON with CommonTrkMETUtil.C

do eflowcommon should be 1 if you want to run Marianna's eflow script. If you then want to use FTK tracks you should create the input list in "tO.eFlow.FTKconv.2.%s.txt" and have the flag doFTKeflowcommon in MET_tO_eflow_2.C uncommentd. Otherwise to run over offlline track pflow comment //#define FTKTRACKS and put the doFTKeflowcommon=0

To run RunFTKmetStandROOT.C

root -l RunFTKmetStandROOT.C

Possible problems

Depending on the ROOT source it is possible that in both the .h and .C

JetAnalysisCalib::

is not neededin the declaration, therefore uncomment //#define laptop in both the .C and .h

Known issues

In the DESD_FTK produced in September

• the FTK_TrackParticle does not have the charge information
• they don't have the necessary trigger information included for MET Yasu is asking for another reprocessing