ZeeJets Analysis with 2011 data

Contents:

Sych. Effort (Matteo et. al.): Feb 20

  • /castor/cern.ch/user/m/marone/ZJets/SynchTest/Nov08RunA44X.root ( entries: 5906)
  • Selection Steps:
    • 1) ask for event to pass trigger: "HLT_Ele17_----------------------------Ele8_---------------------------"
    • 2) if event has "pat::ElectronCollection" size < 2 .... skip it.
    • 3) pick the two leading pf-electrons
    • 4) check them for opposite charge requirement
    • 5) check the invariant mass: |InvMass-91|<20
    • 6) both electrons pass pt>20 and |eta|<2.4 (exclude: >1.446 && <1.566 )
    • 7) both pass ID conditions like:
More... Close
  if(e.isEB()) {
    sieie <0.01 &&
    fabs(delphi)<0.06 &&
    fabs(detain)<0.004 &&
    hoe <0.04 ;
  }
  else if(e.isEE()){
    sieie <0.03 &&
    fabs(delphi)<0.03 &&
    fabs(detain)<0.007&&
    hoe <0.15 ;
  }
    • 8) both the electrons pass det-based isolation as:
More... Close
  Double_t Aecal[2]    = {0.101, 0.046}; //   barrel/endcap
  Double_t Ahcal[2]    = {0.021 , 0.040 }; //   barrel/endcap
  enum detID{barrel=0, endcap=1};

    if(e.isEB()){
        IsoPass = (tiso/pt)< 0.09
        && (max(0. ,eiso - Aecal[barrel]*(lepIsoRho))/pt)<  0.07
        &&  (max(0.,hiso - Ahcal[barrel]*(lepIsoRho))/pt)< 0.1;
        }
  else if(e1.isEE()){
        IsoPass = (tiso/pt)< 0.04
        && (max(0. ,eiso - Aecal[barrel]*(lepIsoRho))/pt)<  0.05
        &&  (max(0.,hiso - Ahcal[barrel]*(lepIsoRho))/pt)< 0.025;
        }
    • 9) both pass conversion rejection cuts:
        ! (fabs(dcot) < 0.02 && fabs(dist) <=0.02 ) && nmshits <=0 

    • 10) PF-Jet Selection and Counting (AK5):

if(runOnMC == False):
        process.GlobalTag.globaltag = cms.string('GR_R_44_V13::All')# autoCond[ 'startup' ] )
else :
        process.GlobalTag.globaltag = cms.string('START44_V12::All')# autoCond[ 'startup' ] )


jetCorrections=('AK5PFchs', ['L1FastJet','L2Relative','L3Absolute','L2L3Residual']

jetID = jet->chargedHadronEnergyFraction() > 0 &&
        jet->chargedEmEnergyFraction()<0.99 &&
        jet->nConstituents()>1 &&
        jet->neutralHadronEnergyFraction()<0.99 &&
        jet->neutralEmEnergyFraction()<0.99 &&
        jet->chargedMultiplicity()>0.0;

     if(jetPt<30)continue;

     if(abs(jetEta)>2.4)continue;

     if(dR(electron,jet)<0.3))continue;

DATA = rfio:/castor/cern.ch/user/m/marone/ZJets/SynchTest/Nov08RunA44X.root (we start with 5906 data-events)

MC = rfio:/castor/cern.ch/user/m/marone/ZJets/SynchTest/DYJetsToLL_TuneZ2_PU_S6_START44.root (we start with 11378 mc-events)

Results:

Now what the results are:
STEP nEvents (data) nEvents (MC)
0 5906 11378
1 2833 1701
2-3 1666 1490
4 1567 1451
5 1362 1343
6 1144 1129
7 910 935
8 795 791
9 753 748
>=1jet 107 131
>=2jet 22 30
>=3jet 2 4
>=4jet 1 0
>=5jet 0 0

Note : No Filter , No JSON for data-file

FROM MATTEO:

1) HLT 2833
2) EleCollection.size() >= 2    2423  
3) WP80 -> 868
4) ele > 20 Gev -> 807
5) opposite charge requid: 802

They have ~6% more for Data.

How do we pick the electron member-functions:

  double q1 = e1.threeCharge();
  double pt1     = e1.pt();
  double eta1     = e1.superCluster()->eta();
  double phi1     = e1.phi();
  double tiso1    = e1.dr03TkSumPt();
  double eiso1    = e1.dr03EcalRecHitSumEt();
  double hiso1    = e1.dr03HcalTowerSumEt();
  double sieie1   = e1.sigmaIetaIeta();
  double delphi1  = e1.deltaPhiSuperClusterTrackAtVtx();
  double detain1  = e1.deltaEtaSuperClusterTrackAtVtx();
  double dcot1    = e1.convDcot();
  double dist1    = e1.convDist();
  double hoe1     = e1.hadronicOverEm();
  double nmshits1 = e1.gsfTrack().get()->trackerExpectedHitsInner().numberOfHits();


Sych. Effort (Kostas et. al.) Feb 7

(Results To be Loaded on Common TWIKI page after More Cross Checks.)

  • CMSSW_4_4_2, PF2PAT
  • /DoubleElectron/Run2011A-05Aug2011-v1/AOD (total events = 6717907)
  • Trigger Paths:
allPathsData =  (
    "HLT_Ele17_CaloIdL_CaloIsoVL_Ele8_CaloIdL_CaloIsoVL_v1",
    "HLT_Ele17_CaloIdL_CaloIsoVL_Ele8_CaloIdL_CaloIsoVL_v2",
    "HLT_Ele17_CaloIdL_CaloIsoVL_Ele8_CaloIdL_CaloIsoVL_v3",
    "HLT_Ele17_CaloIdL_CaloIsoVL_Ele8_CaloIdL_CaloIsoVL_v4",
    "HLT_Ele17_CaloIdL_CaloIsoVL_Ele8_CaloIdL_CaloIsoVL_v5",
    "HLT_Ele17_CaloIdL_CaloIsoVL_Ele8_CaloIdL_CaloIsoVL_v6",
    "HLT_Ele17_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_Ele8_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_v6",
    "HLT_Ele17_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_Ele8_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_v7" ,
    "HLT_Ele17_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_Ele8_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_v8"
     )

 ==  =  Total :
| Delivered LS | Delivered(/pb) | Selected LS | Recorded(/pb) |
---------------------------------------------------------------
|        15472 |        470.389 |       12467 |       368.037 |

  • Global Tag: "GR_R_44_V13::All"
  • Data Processed with some Filters:

process.primaryVertexFilter = cms.EDFilter("VertexSelector",
   src = cms.InputTag("offlinePrimaryVertices"),
   cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.Rho <= 2"),
   filter = cms.bool(True),   # otherwise it won't filter the events, just produce an empty vertex collection.
)

process.noscraping = cms.EDFilter("FilterOutScraping",
applyfilter = cms.untracked.bool(True),
debugOn = cms.untracked.bool(False),
numtrack = cms.untracked.uint32(10),
thresh = cms.untracked.double(0.25)
)

# HBHENoiseFilter
process.load('CommonTools/RecoAlgos/HBHENoiseFilter_cfi')
process.HBHENoiseFilter.maxRBXEMF = cms.double(0.01)

  • Events left here: 4333051 (Finally checked, if I don't apply filters then I get: 4338163 (exact as ETH) )

Event Selection:

  • two leading pf-Electrons with pT >20GeV and |eta|<2.4 (exclude gap-region)
  • both electrons have opposite charge and |inv-mass-91|<20
  • WP80 Id cuts with combined relative isolation <0.15 and evaluated as follows:
  //Evaluate Isolation Cuts:
  double lepIsoRho;// = -9999999999999;
  edm::Handle<double> rhoLepIso;//= -999999999;
  const edm::InputTag eventrhoLepIso("kt6PFJetsForIsolation", "rho");
  iEvent.getByLabel(eventrhoLepIso, rhoLepIso);
  if( *rhoLepIso == *rhoLepIso) lepIsoRho = *rhoLepIso;
  else  lepIsoRho =  999999.9;

//  Double_t Atracker[2] = {0., 0.}; //   barrel/endcap
  Double_t Aecal[2]    = {0.101, 0.046}; //   barrel/endcap
  Double_t Ahcal[2]    = {0.021 , 0.040 }; //   barrel/endcap
  enum detID{barrel=0, endcap=1};

  double combRelRhoSubtractedIso1 = 0;
  double combRelRhoSubtractedIso2 = 0;

  if(e1.isEB()) combRelRhoSubtractedIso1 = (tiso1 + max(0. ,eiso1 - Aecal[barrel]*(lepIsoRho)) + max(0.,hiso1 - Ahcal[barrel]*(lepIsoRho)) )/pt1;
  else if(e1.isEE())combRelRhoSubtractedIso1 = (tiso1 + max(0. ,eiso1 - Aecal[endcap]*(lepIsoRho)) + max(0.,hiso1 - Ahcal[endcap]*(lepIsoRho)) )/pt1;

  if(e2.isEB())combRelRhoSubtractedIso2 = (tiso2 + max(0. ,eiso2 - Aecal[barrel]*(lepIsoRho)) + max(0.,hiso2 - Ahcal[barrel]*(lepIsoRho)) )/pt2;
  else if(e2.isEE())combRelRhoSubtractedIso2 = (tiso2 + max(0. ,eiso2 - Aecal[endcap]*(lepIsoRho)) + max(0.,hiso2 - Ahcal[endcap]*(lepIsoRho)) )/pt2;

  bool bothPassIso = combRelRhoSubtractedIso1<0.15 && combRelRhoSubtractedIso2<0.15;

  • pf-Jets with pt> 30 GeV and |eta|<2.5, dR(e,jet)>0.3

Event Counting table:

Jet-bin Events ETH
>=0Jet 75149 -
>=1Jet 11758 16510
>=2Jet 2468 3543
>=3Jet 471 676
>=4Jet 104 137
>=5Jet 17 26

(Note: Above event counting in my case have NoFilters applied+Isolationcut-pfBased=0.2-the default one+ ETH cuts)


-- LovedeepKaurSaini - 20-Feb-2012

Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2012-02-20 - LovedeepKaurSaini
 
    • 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-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback