%CERTIFY%

TemplateH2012

Introduction

This page as it is now is only meant as a starting point and should keep evolving. It uses the winter2012 HSG6 settings as starting point;

Common selection

Object Selection

object cut cut value comments
electrons Order: LAr hole veto; E smearing; run CaloIso tool (Should we do the LAr hole stuff?)
  alg ElectronAODCollection  
  alg author==1 OR 3 D3PD : el_author (we give up on forward electrons)
  alg cut ElectronMedium++ *...mediumPP
  eT> ET > 25 GeV (author 1 and 3) transverse Energy = cluster E /cosh(track eta); D3PD: if el_nSCTHits[i] + el_nPixHits[i] < 4: ET = el_cl_pt[i] , else : ET = el_cl_E[i]/cosh(el_tracketa[i])
  eta abs(eta_cluster)< 0.0-1.37 OR 1.52-2.47 D3PD: if el_nSCTHits[i] > 3: eta = el_tracketa[i] , else: eta = c.el_cl_eta[i]
  track Isolation sum pT tracks/pT < 0.1 within cone ΔR=0.2. D3PD : el_ptcone20/pt < 0.1 , Not applied for forward electrons.
  calo Isolation sum ET tracks/ET < 0.14 within cone ΔR=0.3. D3PD : el_Etcone30/ET < 0.14, Use pT corrected and NPV corrected values: details. Not applied for forward electrons.
  OQ flag ((el_OQ&1446)==0) avoid dead OTX region; D3PD: el_OQ
muons  
  alg MuidMuonCollection  
  alg cut author==12 and medium D3PD : mu_muid_author and mu_muid_medium
  pT> 25 GeV D3PD : mu_muid_pt
  eta abs(eta)< 2.7 D3PD : mu_muid_eta
  b layer ! expectBLayerHit OR numberOfBLayerHits > 0 D3PD : ! mu_muid_expectBLayerHit OR mu_muid_nBLHits > 0
  hits Number of pixel hits + number of crossed dead pixel sensors > 1 D3PD : mu_muid_nPixHits + mu_muid_nPixelDeadSensors > 1
  hits Number of SCT hits + number of crossed dead SCT sensors > 5 D3PD : mu_muid_nSCTHits + mu_muid_nSCTDeadSensors > 5
  hits Number of pixel holes + number of SCT holes < 3 D3PD : mu_muid_nPixHoles + mu_muid_nSCTHoles < 3
  hits TRT: denote n = nTRTHits + nTRTOutliers:
for abs(eta)<1.9 require n>5 and nTRTOutliers/n < 0.9,
for abs(eta)>=1.9: if n>5 then require nTRToutliers/n < 0.9
D3PD : use mu_muid_nTRTHits , mu_muid_nTRTOutliers and mu_muid_nTRTHits
  track Isolation sum pT tracks/pT < 0.1 within cone ΔR=0.2. D3PD : mu_muid_ptcone20/pt < 0.1
  calo Isolation sum ET tracks/ET < 0.14 within cone ΔR=0.3. D3PD : mu_muid_etcone30/ET < 0.14 (for ET we use mu_muid_pt, notice small e in etcone30)
jets anti-kt4 pt > 20, abs(eta) < 3.2
overlap removal remove 0.4 jets which are closer than 0.2 to good electrons
   
b-tagging  
  alg MV1 use this script. Usage: double mv1Eval(double w_IP3D, double w_SV1, double w_JetFitterCombNN, double jet_pt, double jet_eta)
  cut weight > 0.601713 efficiency: 70%
MET / sumET  
  MET METRefFinal  

Overlap removals

Triggers

  • Lowest un-prescaled lepton triggers:
    • B-J: e20_medium
    • K: e22_medium
    • L-M: e22vh_medium1
    • B-I: mu18_MG
    • J-M: mu18_MG_medium
  • Met triggers:

Data Periods

Lepton and tau+met triggers: MC11b:

period Run number in MC11b Run numbers in data conditions triggers in data and MC fraction in MC11b integrated luminosity fraction in data
B-D 180164 178044-180481 All FEBs ok e20_medium, mu18 3.2% 173.61/pb 3.75%
E-H 183003 180614-184169 6 missing FEBs e20_medium, mu18 17.4% 948.67/pb 20.50%
I-K 186169 185353-187815 2 missing FEBs e22_medium, mu18_medium 25.8% 1112.91/pb 24.05%
L-M 189751 188921-191933 2 missing FEBs e22vh_medium1, mu18_medium 53.5% 2392.85/pb 51.70%

LAr hole

electron isolation

GRL

Use this GRL for the whole 2011 dataset! Or download the latest xml file from this repository. The tag to use is DetStatus-v36-pro10.

Important: when using the GRL, you must also veto the time windows around a LAr noise burst by removing events with larError==2!

More information on how to calculate the integrated luminosity for a given GRL

The total integrated luminosity for the GRL given above: 4.661/fb. More info here

MC Reweighting

MC weights

Pile-up Reweighting

Pile-up Reweighting

The most up-to-date recommendation from the pileup reweighting developers is now available here: [[https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/ExtendedPileupReweighting]ExtendedPileupReweighting]]

Pile-up reweighting 2.0 in mc11b and mc11c

Use these files as input for the tool:
  • The MC11c distributions file: to-do
  • Note that for mc11c is should be not necessary to do mu rescaling.
  • Data: To-do . This corresponds to the GRL linked above.

You need PileupReweighting-00-02-00 and these lines to initialize it

Root::TPileupReweighting* reweight = new Root::TPileupReweighting( "PileUpReweighting" )
reweight->SetUnrepresentedDataAction(1)
reweight->AddConfigFile("path to mc11b_mu.root")
reweight->SetDataScaleFactors(0.83)    //to scale the mu, only for mc11b!!!
reweight->AddLumiCalcFile("path to data mu file")
reweight->initialize()

//then, for each event, get the weight as:
double pileupWeight = reweight->GetCombinedWeight(physis->RunNumber, 105200, physics->averageIntPerXing)

Data and MC samples

Data information

Make sure not to use the Atlfast samples (tag _aXXX instead of _sYYYY).

CP recommendations

Edit | Attach | Watch | Print version | History: r8 < r7 < r6 < r5 < r4 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r8 - 2012-04-02 - unknown
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 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