From ba855e4d0142976c155b708a481d7564da675d40 Mon Sep 17 00:00:00 2001 From: Chi Nhan Nguyen Date: Mon, 23 Jun 2008 13:51:20 +0000 Subject: [PATCH] --- yaml --- r: 46363 b: "refs/heads/CMSSW_7_1_X" c: 678b08fb0126191e66b436436d0069232d415841 h: "refs/heads/CMSSW_7_1_X" i: 46361: 602bdd5a0cb17fc860ea313fb74e156fb348db6b 46359: 32b275faf5c4e8c4b4aa9c84e660015dc058cb33 v: v3 --- [refs] | 2 +- .../HLTanalyzers/test/rates/OHltTree.cpp | 22 +++++------ .../HLTanalyzers/test/rates/OHltTree.h | 4 +- .../HLTrigger/HLTanalyzers/test/rates/README | 38 +++++++++++++++++++ 4 files changed, 52 insertions(+), 14 deletions(-) diff --git a/[refs] b/[refs] index 7365f52a5c4c6..2c48560670713 100644 --- a/[refs] +++ b/[refs] @@ -1,3 +1,3 @@ --- refs/heads/gh-pages: ac155dadd83efa75cad55c0508a57a2b9dd3d66c -"refs/heads/CMSSW_7_1_X": 4a12e59caefa990485fb79daa825cf39a51053d0 +"refs/heads/CMSSW_7_1_X": 678b08fb0126191e66b436436d0069232d415841 diff --git a/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.cpp b/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.cpp index 9d765e463ff5f..2475e33dd8e5c 100644 --- a/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.cpp +++ b/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.cpp @@ -282,7 +282,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector } else if (trignames[it].CompareTo("ElectronTau") == 0) { if ( L1_IsoEG10_TauJet20==1 ) { // L1 Seed - if(OpenHltTauPassed(15.,5.,3.,1,0.,0)>=1 && OpenHlt1ElectronPassed(10,1,0.06,3.,0.,0.)>=1 ) { + if(OpenHltTauPassed(15.,5.,3.,1,0.,0)>=1 && OpenHlt1ElectronPassed(10.,1,0.06,3.)>=1 ) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -292,7 +292,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector } else if (trignames[it].CompareTo("ElectronTau_NoL1") == 0) { if ( L1_IsoEG10_Jet15==1 ) { // L1 Seed - if(OpenHltTauPassed(15.,5.,3.,1,0.,0)>=1 && OpenHlt1ElectronPassed(10,1,0.06,3.,0.,0.)>=1 ) { + if(OpenHltTauPassed(15.,5.,3.,1,0.,0)>=1 && OpenHlt1ElectronPassed(10.,1,0.06,3.)>=1 ) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -302,7 +302,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector } else if (trignames[it].CompareTo("ElectronTau_NoL2") == 0) { if ( L1_IsoEG10_TauJet20==1 ) { // L1 Seed - if(OpenHltTauPassed(0.,999.,3.,1,0.,0)>=1 && OpenHlt1ElectronPassed(10,1,0.06,3.,0.,0.)>=1 ) { + if(OpenHltTauPassed(0.,999.,3.,1,0.,0)>=1 && OpenHlt1ElectronPassed(10.,1,0.06,3.)>=1 ) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -312,7 +312,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector } else if (trignames[it].CompareTo("ElectronTau_NoL25") == 0) { if ( L1_IsoEG10_TauJet20==1 ) { // L1 Seed - if(OpenHltTauPassed(15.,5.,0.,0,0.,0)>=1 && OpenHlt1ElectronPassed(10,1,0.06,3.,0.,0.)>=1 ) { + if(OpenHltTauPassed(15.,5.,0.,0,0.,0)>=1 && OpenHlt1ElectronPassed(10.,1,0.06,3.)>=1 ) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -322,7 +322,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector } else if (trignames[it].CompareTo("ElectronTau_NoSi") == 0) { if ( L1_IsoEG10_TauJet20==1 ) { // L1 Seed - if(OpenHltTauPassed(15.,5.,0.,0,0.,0)>=1 && OpenHlt1ElectronPassed(10,1,0.06,3.,0.,0.)>=1 ) { + if(OpenHltTauPassed(15.,5.,0.,0,0.,0)>=1 && OpenHlt1ElectronPassed(10.,1,0.06,3.)>=1 ) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -332,7 +332,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector } else if (trignames[it].CompareTo("ElectronMET") == 0) { if ( L1_SingleIsoEG12==1 ) { // L1 Seed - if(OpenHlt1ElectronPassed(10,1,0.06,3.,0.,0.)>=1 && recoMetCal>=35.) { + if(OpenHlt1ElectronPassed(10.,1,0.06,3.)>=1 && recoMetCal>=35.) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -350,7 +350,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector else if (trignames[it].CompareTo("OpenHLT1Electron") == 0) { if ( L1_SingleIsoEG12==1 ) { // L1 Seed //PrintOhltVariables(3,electron); - if(OpenHlt1ElectronPassed(15,1,0.06,3.,1.5,2.45)>=1) { + if(OpenHlt1ElectronPassed(15.,1,0.06,3.)>=1) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -367,7 +367,7 @@ void OHltTree::Loop(vector * iCount, vector * sPureCount, vector else if (trignames[it].CompareTo("OpenHLT1Photon") == 0) { if ( L1_SingleIsoEG12==1 ) { // L1 Seed //PrintOhltVariables(3,photon); - if(OpenHlt1PhotonPassed(30,1,0,1.5,6.,4.)>=1) { + if(OpenHlt1PhotonPassed(30.,1,0,1.5,6.,4.)>=1) { triggerBitNoPrescale[it] = true; if ((iCountNoPrescale[it]) % map_pathHLTPrescl.find(trignames[it])->second == 0) { triggerBit[it] = true; @@ -596,7 +596,7 @@ int OHltTree::OpenHltTauPassed(float Et,float Eiso, float L25Tpt, int L25Tiso, f } -int OHltTree::OpenHlt1ElectronPassed(double Et, double L1iso, double Tiso, double Hiso, double eoverpBR, double eoverpEC) +int OHltTree::OpenHlt1ElectronPassed(float Et, int L1iso, float Tiso, float Hiso) { int rc = 0; // Loop over all oh electrons @@ -614,7 +614,7 @@ int OHltTree::OpenHlt1ElectronPassed(double Et, double L1iso, double Tiso, doubl return rc; } -int OHltTree::OpenHlt1PhotonPassed(double Et, double L1iso, double Tiso, double Eiso, double HisoBR, double HisoEC) +int OHltTree::OpenHlt1PhotonPassed(float Et, int L1iso, float Tiso, float Eiso, float HisoBR, float HisoEC) { int rc = 0; // Loop over all oh photons @@ -622,7 +622,7 @@ int OHltTree::OpenHlt1PhotonPassed(double Et, double L1iso, double Tiso, double if ( ohPhotEt[i] > Et) if ( ohPhotL1iso[i] >= L1iso ) if( ohPhotTiso[i]<=Tiso ) - if( ohPhotEiso[i] < Eiso /* && isEgammaL1_PhotonSuperclusterMatched(i)*/ ) { + if( ohPhotEiso[i] < Eiso ) { if( (TMath::Abs(ohPhotEta[i]) < 1.5 && ohPhotHiso[i] < HisoBR ) || (1.5 < TMath::Abs(ohPhotEta[i]) && TMath::Abs(ohPhotEta[i]) < 2.5 && ohPhotHiso[i] < HisoEC ) ) rc++; diff --git a/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.h b/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.h index bd76f3ed6487c..ee6dbc78a1228 100644 --- a/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.h +++ b/trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.h @@ -1082,8 +1082,8 @@ public : ,int n=-1,bool doMuonCut=false,bool doElecCut=false, double muonPt=-999., double muonDr=-999.); void PrintOhltVariables(int level, int type); int OpenHltTauPassed(float Et,float Eiso, float L25Tpt, int L25Tiso,float L3Tpt, int L3Tiso); - int OpenHlt1ElectronPassed(double Et,double L1iso,double Tiso,double Hiso,double eoverpBR,double eoverpEC); - int OpenHlt1PhotonPassed(double Et,double L1iso,double Tiso,double Eiso,double HisoBR,double HisoEC); + int OpenHlt1ElectronPassed(float Et,int L1iso,float Tiso,float Hiso); + int OpenHlt1PhotonPassed(float Et,int L1iso,float Tiso,float Eiso,float HisoBR,float HisoEC); int OpenHlt1MuonPassed(double ptl1,double ptl2,double ptl3,double dr,int iso); int OpenHlt2MuonPassed(double ptl1,double ptl2,double ptl3,double dr,int iso); int OpenHlt1JetPassed(double pt); diff --git a/trunk/HLTrigger/HLTanalyzers/test/rates/README b/trunk/HLTrigger/HLTanalyzers/test/rates/README index f359d5281beb1..53dcfb98b8431 100644 --- a/trunk/HLTrigger/HLTanalyzers/test/rates/README +++ b/trunk/HLTrigger/HLTanalyzers/test/rates/README @@ -155,3 +155,41 @@ README for OpenHLT rate calculating scripts * Note that in CMSSW releases 2_X and higher, the L2 and L3 muon pT thresholds are implemented as flat cuts. In earlier releases, the thresholds were implemented as 90% CL upper limits. OpenHLT examples of the old pT threshold implementation are found at the beginning of the OpenHlt1MuonPassed and OpenHlt2MuonPassed examples. + + +-------------------------------------- +3d. Example OpenHLT electron triggers +-------------------------------------- + + One example OpenHLT electron trigger (HLT1Electron) are available. + + The algorithm for the single electron triggers is implemented in the OHltTree::OpenHlt1ElectronPassed method. + This takes five arguments: + + float Et: Et threshold in GeV + int L1iso: L1 is isolated (1: yes, 0: no) + float Tiso: Track isolation in GeV + float Hiso: Hcal isolation in GeV + + In the example code some cuts are mandatory and cannot be varied. + They have the following default values: + ohEleNewSC[i]>0 : Is 1, if electron found by Super Cluster Algo, otherwise 0 + ohEleHiso[i]/ohEleEt[i] < 0.05 : Hadronic fraction smaller than 0.05 + ohElePixelSeeds[i]>0 : Pixel seed found or not (1 or 0) + +-------------------------------------- +3e. Example OpenHLT photon triggers +-------------------------------------- + + One example OpenHLT electron trigger (HLT1Electron) are available. + + The algorithm for the single electron triggers is implemented in the OHltTree::OpenHlt1ElectronPassed method. + This takes five arguments: + + float Et: Et threshold in GeV + int L1iso: L1 is isolated (1: yes, 0: no) + float Tiso: Track isolation in GeV + float Eiso: Ecal isolation in GeV + float HisoBR: Barrel Hcal isolation in GeV + float HisoEC: Endcap Hcal isolation in GeV +