Giter Club home page Giter Club logo

ele_vid's Introduction

Ele_VID

Learning the Ele VID method and correlated information, follow https://twiki.cern.ch/twiki/bin/view/Main/VIDTutorial2017:

Getting the code //install release

scram project -n CMSSW-8-0-26-p1-VID CMSSW CMSSW_8_0_26_patch1
cd CMSSW-8-0-26-p1-VID/src  
cmsenv
voms-proxy-init -voms cms

copy needed cmssw code

git cms-init
git cms-addpkg RecoEgamma/ElectronIdentification
git cms-addpkg PhysicsTools/SelectorUtils
git clone https://github.com/ikrav/EgammaWork.git
cd EgammaWork/
git checkout ntupler_and_VID_demos_8.0.3
cd ../

and change the input rootfile in runElectrons_VID_CutBased_Summer16_80X_demo.py:

inputFilesMiniAOD = cms.untracked.vstring('file:/eos/cms/store/mc/RunIISummer16MiniAODv2/DY2JetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/80000/E0976CCD-33BE-E611-879D-E41D2D08DFE0.root'
    )

then cp the ElectronNtuplerVIDDemo.cc to the plugins directory,

cd $CMSSW_BASE/src
scram b

In order to check that everything is ok, run the demo ntuplizer:

cd EgammaWork/ElectronNtupler/test
cmsRun runElectrons_VID_CutBased_Summer16_80X_demo.py

the output is like following:

Begin processing the 64th record. Run 1, Event 381641, LumiSection 256 at 18-Jan-2021 07:53:45.361 CET
MENG test---->>>ele size---->>> 2ele pt---->>> 83.220009
IdBit:2

Meng test--->>>1019

DEBUG CutFlow, full info for cand with pt=83.220009:
    CutFlow name= cutBasedElectronID-Summer16-80X-V1-medium    decision is 0
 Index                               cut name              isMasked    value-cut-upon     pass?
   0                                              MinPtCut_0    0        83.220009          1
   1                              GsfEleSCEtaMultiRangeCut_0    0        0.245101          1
   2                                   GsfEleDEtaInSeedCut_0    0        0.004014          0
   3                                       GsfEleDPhiInCut_0    0        0.011771          1
   4                         GsfEleFull5x5SigmaIEtaIEtaCut_0    0        0.008891          1
   5                               GsfEleHadronicOverEMCut_0    0        0.015428          1
   6                        GsfEleEInverseMinusPInverseCut_0    0        0.017043          1
   7                                 GsfEleEffAreaPFIsoCut_0    0        0.021371          1
   8                               GsfEleConversionVetoCut_0    0        1.000000          1
   9                                  GsfEleMissingHitsCut_0    0        0.000000          1
DEBUG CutFlow, the result with cut GsfEleFull5x5SigmaIEtaIEtaCut_0 masked out
    CutFlow name= cutBasedElectronID-Summer16-80X-V1-medium    decision is 0
 Index                               cut name              isMasked    value-cut-upon     pass?
   0                                              MinPtCut_0    0        83.220009          1
   1                              GsfEleSCEtaMultiRangeCut_0    0        0.245101          1
   2                                   GsfEleDEtaInSeedCut_0    0        0.004014          0
   3                                       GsfEleDPhiInCut_0    0        0.011771          1
   4                         GsfEleFull5x5SigmaIEtaIEtaCut_0    1        0.008891          1
   5                               GsfEleHadronicOverEMCut_0    0        0.015428          1
   6                        GsfEleEInverseMinusPInverseCut_0    0        0.017043          1
   7                                 GsfEleEffAreaPFIsoCut_0    0        0.021371          1
   8                               GsfEleConversionVetoCut_0    0        1.000000          1
   9                                  GsfEleMissingHitsCut_0    0        0.000000          1
ele size---->>> 2ele pt---->>> 57.086735
IdBit:2

Meng test--->>>1003

DEBUG CutFlow, full info for cand with pt=57.086735:
    CutFlow name= cutBasedElectronID-Summer16-80X-V1-medium    decision is 0
 Index                               cut name              isMasked    value-cut-upon     pass?
   0                                              MinPtCut_0    0        57.086735          1
   1                              GsfEleSCEtaMultiRangeCut_0    0        0.735629          1
   2                                   GsfEleDEtaInSeedCut_0    0        0.005562          0
   3                                       GsfEleDPhiInCut_0    0        0.004975          1
   4                         GsfEleFull5x5SigmaIEtaIEtaCut_0    0        0.010015          0
   5                               GsfEleHadronicOverEMCut_0    0        0.000000          1
   6                        GsfEleEInverseMinusPInverseCut_0    0        0.003293          1
   7                                 GsfEleEffAreaPFIsoCut_0    0        0.000000          1
   8                               GsfEleConversionVetoCut_0    0        1.000000          1
   9                                  GsfEleMissingHitsCut_0    0        0.000000          1
DEBUG CutFlow, the result with cut GsfEleFull5x5SigmaIEtaIEtaCut_0 masked out
    CutFlow name= cutBasedElectronID-Summer16-80X-V1-medium    decision is 0
 Index                               cut name              isMasked    value-cut-upon     pass?
   0                                              MinPtCut_0    0        57.086735          1
   1                              GsfEleSCEtaMultiRangeCut_0    0        0.735629          1
   2                                   GsfEleDEtaInSeedCut_0    0        0.005562          0
   3                                       GsfEleDPhiInCut_0    0        0.004975          1
   4                         GsfEleFull5x5SigmaIEtaIEtaCut_0    1        0.010015          1
   5                               GsfEleHadronicOverEMCut_0    0        0.000000          1
   6                        GsfEleEInverseMinusPInverseCut_0    0        0.003293          1
   7                                 GsfEleEffAreaPFIsoCut_0    0        0.000000          1
   8                               GsfEleConversionVetoCut_0    0        1.000000          1
   9                                  GsfEleMissingHitsCut_0    0        0.000000          1

the illustraion of the output: ele size---->>> 2 means this event has two electrons

ele pt---->>> 83.220009 means the pt of the first electron

IdBit:2 each electron has one Id bit

CutFlow name= cutBasedElectronID-Summer16-80X-V1-medium means the id in file cutBasedElectronID_Summer16_80X_V1_cff.py which contains several ids in 80X

Meng test--->>>1019 the 1019 means the bits of the ten cutflow, i.e., 1111111011(index 0 is the rightmost bit)

ele_vid's People

Contributors

menglu21 avatar

Watchers

 avatar

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.