Skip to content

Commit

Permalink
---
Browse files Browse the repository at this point in the history 
yaml
---
r: 46363
b: "refs/heads/CMSSW_7_1_X"
c: 678b08f
h: "refs/heads/CMSSW_7_1_X"
i:
  46361: 602bdd5
  46359: 32b275f
v: v3
  • Loading branch information
Chi Nhan Nguyen committed Jun 23, 2008
1 parent 5e786ba commit ba855e4
Show file tree
Hide file tree
Showing 4 changed files with 52 additions and 14 deletions.
2 changes: 1 addition & 1 deletion [refs]
View file
Original file line number Diff line number Diff line change
@@ -1,3 +1,3 @@
---
refs/heads/gh-pages: ac155dadd83efa75cad55c0508a57a2b9dd3d66c
"refs/heads/CMSSW_7_1_X": 4a12e59caefa990485fb79daa825cf39a51053d0
"refs/heads/CMSSW_7_1_X": 678b08fb0126191e66b436436d0069232d415841
22 changes: 11 additions & 11 deletions trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -282,7 +282,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
}
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;
Expand All @@ -292,7 +292,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
}
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;
Expand All @@ -302,7 +302,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
}
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;
Expand All @@ -312,7 +312,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
}
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;
Expand All @@ -322,7 +322,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
}
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;
Expand All @@ -332,7 +332,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
}
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;
Expand All @@ -350,7 +350,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
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;
Expand All @@ -367,7 +367,7 @@ void OHltTree::Loop(vector<int> * iCount, vector<int> * sPureCount, vector<int>
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;
Expand Down Expand Up @@ -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
Expand All @@ -614,15 +614,15 @@ 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
for (int i=0;i<NohPhot;i++) {
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++;
Expand Down
4 changes: 2 additions & 2 deletions trunk/HLTrigger/HLTanalyzers/test/rates/OHltTree.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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);
Expand Down
38 changes: 38 additions & 0 deletions trunk/HLTrigger/HLTanalyzers/test/rates/README
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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

0 comments on commit ba855e4

Please sign in to comment.