TWiki> Main Web>TWikiUsers>StefanoManzoni>GammaGammaPlusMETNotes (2016-02-04, GiovanniMarchiori) EditAttachPDF

GammaGammaPlusMETNotes



References

Mailing lists:

Meetings

Supporting note

Run2 analysis

DC14 MC

  • signal: JIRA prodsys
    • fullsim (10k each):
      • mc14_13TeV.205027.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_150.merge.AOD.e3503_s1982_s2008_r5787_r5853
      • mc14_13TeV.205028.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_600.merge.AOD.e3503_s1982_s2008_r5787_r5853
      • mc14_13TeV.205029.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_1300.merge.AOD.e3503_s1982_s2008_r5787_r5853
    • fastsim (10k each):
      • mc14_13TeV.205027.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_150.merge.AOD.e3503_a266_a265_r5853
      • mc14_13TeV.205028.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_600.merge.AOD.e3503_a266_a265_r5853
      • mc14_13TeV.205029.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_1300.merge.AOD.e3503_a266_a265_r5853
  • bkg:

MC15 MC

Event selection

Derivations

  • There is no dedicated SUSY derivation for this analysis at the moment
  • We are investigating if we can use HIGG1D1 (currently only trigger cuts: "EF_g35_loose_g25_loose","EF_2g20vh_medium") or EXOT10 (trigger: HLT_g35_loose1_g25_loose1 || HLT_g35_medium1_g25_medium1 || HLT_g35_medium1_g25_medium1_L12EM15V || HLT_2g20_tight1 || HLT_2g50_loose; offline: 2e, 1e1g or 2g with pT>30 GeV)
  • We may write our own if these 2 are too big
  • To test these two follow the instructions below

Initial setup

  • rel19: 
    asetup rel_3,19.1.X.Y-VAL,AtlasDerivation,slc6,here
# or the cache you want: 19.1.4.6, 19.1.4.7
  • rel 20 
    asetup rel_3,20.1.X.Y-VAL,AtlasDerivation,slc6,here

Create work area 

setupWorkArea.py 
cd WorkArea/cmt 
cmt bro cmt config 
cmt bro gmake 
cd ../run

Run derivation code 

Reco_tf.py --inputAODFile xAOD_ggMET.pool.root --outputDAODFile test --reductionConf HIGG1D1,EXOT10

One signal sample is available on lxplus0069 in /tmp/gmarchio/mc14_13TeV.205027.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_150.merge.AOD.e3503_s1982_s2008_r5787_r5853_tid04845088_00

So for instance do ln -s /tmp/gmarchio/mc14_13TeV.205027.Herwigpp_UEEE4_CTEQ6L1_GGM_gl_bino_1400_150.merge.AOD.e3503_s1982_s2008_r5787_r5853_tid04845088_00/AOD.04845088._000005.pool.root.1 xAOD_ggMET.pool.root from the run directory before running. This file contains 5k events. The output is stored in /afs/cern.ch/work/g/gmarchio/public/ggplusMET/DAOD_EXOT10.root and /afs/cern.ch/work/g/gmarchio/public/ggplusMET/DAOD_HIGG1D1.root

Analysis code

Run1 analysis

Bkg MC sample

Event selection

Photon Selection

Pt requirement pt>75 GeV
Eta range fabs(ph_etas2)<1.37 OR 1.52<fabs(ph_etas2)<2.37
Tightness PhotonTightAR
Calorimeter isolation corrected topoEtcone40 < 4 GeV
Photon cleaning Reject if (ph_OQ&134217728!=0 AND (ph_reta>0.98 OR ph_rphi>1.0 OR (ph_OQ&67108864)!=0) )
Overlapping removal If the clusters of a photon and an electron are found within ∆R < 0.01. The object is interpreted as an electron and the photon is removed from the event
If the clusters of a photon and a jet are found within ∆R < 0.21. The object is interpreted as a photon and the jet is removed from the event
If a jet and a photon are found within 0.2 < ∆R < 0.4, the object is interpreted as jet and the photon is removed.

Electron Selection

There is no requirement to have any electrons in the signal selection, e.g. events with zero or any other number of electrons are selected.

Pt requirement pt>20 GeV
Eta range fabs(el_etas2)<1.37 OR 1.52<fabs(el_etas2)<2.37
Tightness ElectronMedium ++
Calorimeter isolation ptcone30/ET < 0.16
z0 requirement fabs(z0 sin(theta))<0.4 mm
Impact parameter significance fabs(d0/sigmad0)<5.0
Electron cleaning (el_OQ&1446)==0
Overlapping Removal If the clusters of a photon and an electron are found within ∆R < 0.01. The object is interpreted as an electron and the photon is removed from the event
  If a jet and an electron are found within ∆R < 0.2, the object is interpreted as an electron and the overlapping jet is removed from the event
  If a jet and an electron are found within 0.2 < ∆R < 0.4, the object is interpreted as jet and the electron is removed

Muon Selection

There is no requirements on the presence of muons in the final signal selection, e.g. events with zero or any other number of muons are selected.

For this analysis muons from the STACO muon algorithm chain are used, Combined or Segment-tagged.

Pt requirement pt>20 GeV
Eta range fabs(mu_eta)<2.4
Tightness Loose
B layer NOT expectBLayerHit OR nBLHits>0
Pixel hits nPixHits+nPixelDeadSensors>1
SCT hits nSCTHits+nSCTDeadSensors>6
Holes nPixHoles+nSCTHoles<3
TRT outliers N=nTRTOutliers+nTRTHits
if fabs(η) <1.9: require N>5 AND nTRTOutliers<0.9N
if fabs(η) >=1.9: if N>5 require nTRTOutliers<0.9N
Calorimeter isolation Etcone30/Et < 0.12
z0 requirement fabs(z0 sin(theta))<0.4 mm
Impact parameter significance fabs(d0/sigmad0)<3.0
Overlapping Removal If a jet and a muon are found within ∆R < 0.4, the object is interpreted as jet and the muon is removed

Jet selection

There is no requirement to have any number of jets in the final signal selection, e.g. events with zero or any other number of jets are selected.

Jets are reconstructed using the AntiKt4LCTopo jet algorithm.

Pt requirement pT > 40 GeV
Eta range fabs(eta)<2.8

Event Veto

  • If the event contains jets after overlap removal and at least one of the jets isBadLooseMinus

  • The event has at least one selected muon with a fabs(z0) > 1 mm or fabs(d0) > 0.2 mm, where the values are calculated with respect to the primary vertex

-- StefanoManzoni - 2014-12-11

Edit | Attach | Watch | Print version | History: r29 < r28 < r27 < r26 < r25 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r29 - 2016-02-04 - GiovanniMarchiori
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback