Run154822/May reprocessing

Samples used for cut comparison :

  • merge.AOD : data10_7TeV.00154822.physics_L1Calo.merge.AOD.r1299_p161

Please see SusyObjectDefinition-2010-01.

Cut flow table

Cut number Cut information events for Run 154822
Cut 0 no cut 849966
Cut 1 apply GRL and L1 calo stream 789992
Cut 2 pass L1_J5 164389
Cut 3 jet definition for cut 3a, 3b, 3c : The series of cut3 is applied to all AntiKt4Topo jets with p_T > 10 GeV (EM-scale)
Please note : These cuts are applied at EM-scale (no JES applied yet), no eta cut, before any overlap removal
Cut 3a no jet with
   (n90<=5 && HECf>0.8)  
Cut 3b require all jets with calorimeter timing
Cut 3c require all jets with jet quality :
 abs(quality)<=0.8 || emf<=0.95 
Cut 4 Loop over all primary vertex and check there is >= 1 primary vertex with >4 tracks (no vertex type is required) 163712
Cut 5 Crack veto : use electron eta (4-momenta) 163675
Cut6-12 Jet definitions for cuts 6-12: consider AntiKt4Topo jets  
- at EM+JES scale  
- satisfying susy acceptance selection cuts : PT > 20 GeV, abs(η) < 2.5  
- after overlap removal (see link)  
Cut 6 Lepton veto (Veto 1 isolated lepton with p_T > 10 GeV) 163330
Cut 9 Leading jet pt(EM+JES) > 70 GeV 781
Cut 10 Jet selection 2 jets in eta region with p_T(EM+JES) > 30 GeV 653
Cut 11 Jet selection 3 jets in eta region with p_T(EM+JES) > 30 GeV 175
Cut 12 Jet selection 4 jets in eta region with p_T(EM+JES) > 30 GeV 40

Some Remarks

  • Use AntiKt4Topo (jet finding algo + constituent )
    • Use EM-scale (no jet level calibration factor) for cleaning cut.
    • Use EMJES, not H1(GCW+JES as we used for the previous MC comparison) for offline selection
  • Jet cleaning variables are calculated using EM scale energy, irrelevant to constituent- and jet- level calibration.
  • mET definition : use MET_Topo(EM scale) only in this cut flow at the moment
  • Use EMJES scale=1 for jets with EMJES=0 due to the bug (Although, you will get the identical numbers in the table even you use EMJES=0) --> This will be fixed in the next production and comparison

For AOD/ESD users

Please let me know if you find any problem.

May reprocessing (r1299_p161) + Tier0(f2xx) data (not frozen yet)

To recover EMJES=0 problem, you need

  • >=
  • JetCalibTools-00-00-70 (if <
  • JetUtils-01-01-27
  • AOD or ESD

Then, apply correct EMJES on ROOT (macro). See JetEtmissDataAnalysisRecommendationSummer2010#Jet_reconstruction_and_calibration

Some remarks:

  • Not rerun jets on AOD, otherwise jetQuality will be replaced by a jetQuality calculated with clusters (not cells).
  • If you want to use ESD, see April reprocessing part (below)

April reprocessing (r1239_p134) + Tier0(f249-f251) data

To recover EMJES=0 and jetQuality, you need

  • >=
  • JetCalibTools-00-00-70
  • JetUtils-01-01-27
  • ESD

Then, reconstruct jets on the fly:

from JetRec.JetGetters import *  
make_StandardJetGetter('AntiKt', 0.4, 'H1Topo', calibName='NumInv', numInvBase='EM', outputCollectionName='AntiKt4EMJESTopoJets')

Alternatively, you can apply correct EMJES on ROOT (macro). See JetEtmissDataAnalysisRecommendationSummer2010#Jet_reconstruction_and_calibration

and access jetQuality by

double jetQuality = JetCaloQualityUtils::jetQuality(jet); // fbadQ 

How to access variables to define bad/ugly jets

// bad jets
double jetQuality = jet->getJetQuality(); // fbadQ 
double time= jet->getJetTime(); // JetCaloQualityUtils::jetTime(jet) gives with sign (+/-) and more precision for ESD case
int n90cells = jet->getMoment("n90") // n90
double hecf = JetCaloQualityUtils::hecF(jet); // hecf
double emf = JetCaloHelper::jetEMFraction(jet);  // emf 
// ugly jets
double tileGap3F = JetCaloQualityUtils::tileGap3F(jet); // tileGap3F
double fcor = jet->getMoment("BCH_CORR_CELL");
// others
int SamplingMax;
double fracSamplingMax = JetCaloQualityUtils(jet, SamplingMax);

Major updates:
-- NaokoKanaya - 26-May-2010

%REVIEW% Never reviewed

Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r4 - 2010-06-01 - NaokoKanaya
    • 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-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback