Acceptance Challenge for Gen5
Last update 01/17/2005
Sample |
First file | Last file | All files |
| ttop0z | Section 1 | Section 40 | All |
| ttopel | Section 1 | Section 213 | All |
| ttopvh | Section 1 | Section 40 | All |
| ttophl | Section 1 | Section 250 | All |
STAGE 0, INITIAL: Clasify events by jet counting
int level = 4;
// Scale + Time + Relative + Multiple interactions
int nvertices = nZV12;
//Must read in and decode z vertex block above
int coneSize=0;
// 0=0.4, 1=0.7 and 2=1.0
int version=5;
// 5 = Monte Carlo(5.3.1) , 2=484 REMAKE data 4=4111 REMAKE data 0= 4111
MC
if(IsData()) version = 5;
int syscode=0;
// 0=correction, other values for systematic studies
int run = event->runNumber;
int imode = 0;
if(IsData()) imode = 1;
// imode is a flag for data=1, MC=0
JetEnergyCorrections b("JetCorrections",
"JetCorrections",level, nvertices, coneSize, version, syscode, run, imode);
STAGE 1, OBSV: MC Z position fiducial
fabs(pr->mcz) <= 60.0;
STAGE 2, TRIGGER:
bool fired_cem
= true;
bool fired_pem
= true;
//no pem trigger check coded yet
bool fired_cmup
= true;
bool fired_cmx
= true;
if(IsData() && GetTrigger() ){
fired_cem = mytool->testTriggerBitL3(event->runNumber,mytrigInfo,0);
fired_cmup= mytool->testTriggerBitL3(event->runNumber,mytrigInfo,1);
fired_cmx = mytool->testTriggerBitL3(event->runNumber,mytrigInfo,2);
if(event->runNumber
< 150145) fired_cmx = false;
//CMX trigger only good after this point
}
STAGE 3, >=1TLP:
at least one tight lepton
STAGE 4, MET:
corrected MET > 20
if(corrEt >= 8.0 && fabs(jt->EtaDetector) <=2.4){
corr_metx += (1.0-scale)*jt->Et*cos(jt->Phi);
//kohei way
corr_mety += (1.0-scale)*jt->Et*sin(jt->Phi);
//kohei way
}
STAGE 5, ISO: lepton isolation
ele->Isol < 0.1;
STAGE 6, !DILEP: dilepton veto
(smry->fnTightLepton+ smry->fnLooseLepton) >= 2 ;
STAGE 7, !Z: z veto
smry->fTopEventClass
& 0x8) == 0;
STAGE 8, !CONV: conversion and cosmic veto
ele->Conversion ==1
smry->fTopEventClass
& 0x20) == 0;
STAGE 9, !DIFFZ: different interaction veto
fabs(TightLepton->Z0
- zvertex) > 5.0;
STAGE 10, !QCD: QCD veto
if(nTightJet>=1){
TLorentzVector jlead(TightJet[0]->Px,TightJet[0]->Py, TightJet[0]->Pz,TightJet[0]->En);
double dphi = calcDeltaPhi(nu,jlead);
if(corrMetMag < 30.0 && ( dphi < 0.5 || dphi > 2.5)) isQCD
= true;
}
STAGE 11, SI GOOD: good run with silicon
bool
goodrun_si = true;
if(IsData() && KeepGoodRuns()) {
goodrun_nosi = mytool->good_run_nosi(event->runNumber);
goodrun_si = mytool->good_run_si(event->runNumber);
}
STAGE 12, TAGGABLE: events with taggable jets
TightJet->secvNtrk>=2;
STAGE 13, >=1+BTAG: events with, at least, one positive tagged jet
TightJet->secvTag== 1
STAGE 14, >=2+BTAG: events with, at least, two positive tagged jet
Two
jets with TightJet->secvTag== 1;
STAGE 15, >=1-BTAG: events with, at least, one negative tagged jet
TightJet->secvTag==-1;
STAGE 16, TAGGABLE JP+: events with positive JP taggable jets
TightJet->jpbNGdRPTrkPos>=2;
STAGE 17, >=1+ JPBTAG: events with, at least, one positive JP tagged jet
TightJet->jpbRPpos<0.01;
STAGE 18, >=2+ JPBTAG: events with, at least, two positive JP tagged jets
Two
jets with TightJet->jpbRPpos<0.01;
STAGE 19, TAGGABLE JP-: events with negative JP taggable jets
TightJet->jpbNGdRPTrkNeg>=2;
STAGE 20, >=1- JPBTAG: events with, at least, one negative JP tagged jet
TightJet->jpbRPneg<0.01;