void dumpStack(StackType *stk)
{
int i, j;
printf(" ------------ STACK -----------\n");
for (i=0; i<stk->numFrames; ++i) {
printf(" ---------- FRAME #%d: %s ----------\n", i, stk->frames[i].funcName);
for (j=0; j<stk->frames[i].numParms; ++j) {
printf(" -- Parm #%d:\n", j);
dumpVar(&(stk->frames[i].parms[j]));
}
}
printf(" -------- END OF STACK --------\n\n");
}
void QSpyWrk::newTrackHeader( const QString& name,
int trackno, int channel, int pitch,
int velocity, int port,
bool selected, bool muted, bool loop )
{
dump(0, trackno, QString::number(channel), "Track", name);
if (m_verbosity) {
dumpVar("pitch", pitch);
dumpVar("velocity",velocity);
dumpVar("port", port);
dumpVar("selected", selected);
dumpVar("muted", muted);
dumpVar("loop", loop);
}
}
void QSpyWrk::trackHeader( const QString& name1, const QString& name2,
int trackno, int channel, int pitch,
int velocity, int port,
bool selected, bool muted, bool loop )
{
dump(0, trackno, QString::number(channel), "Track", QString("name1='%2' name2='%3'").arg(name1).arg(name2));
if (m_verbosity) {
dumpVar("pitch", pitch);
dumpVar("velocity",velocity);
dumpVar("port", port);
dumpVar("selected", selected);
dumpVar("muted", muted);
dumpVar("loop", loop);
}
}
void QSpyWrk::globalVars()
{
dumpStr(0, "Global Vars", QString());
if (m_verbosity) {
dumpVar("Now", m_engine->getNow());
dumpVar("From", m_engine->getFrom());
dumpVar("Thru", m_engine->getThru());
dumpVar("KeySig", m_engine->getKeySig());
dumpVar("Clock", m_engine->getClock());
dumpVar("AutoSave", m_engine->getAutoSave());
dumpVar("PlayDelay", m_engine->getPlayDelay());
dumpVar("ZeroCtrls", m_engine->getZeroCtrls());
dumpVar("SendSPP", m_engine->getSendSPP());
dumpVar("SendCont", m_engine->getSendCont());
dumpVar("PatchSearch", m_engine->getPatchSearch());
dumpVar("AutoStop", m_engine->getAutoStop());
dumpVar("StopTime", m_engine->getStopTime());
dumpVar("AutoRewind", m_engine->getAutoRewind());
dumpVar("RewindTime", m_engine->getRewindTime());
dumpVar("MetroPlay", m_engine->getMetroPlay());
dumpVar("MetroRecord", m_engine->getMetroRecord());
dumpVar("MetroAccent", m_engine->getMetroAccent());
dumpVar("CountIn", m_engine->getCountIn());
dumpVar("ThruOn", m_engine->getThruOn());
dumpVar("AutoRestart", m_engine->getAutoRestart());
dumpVar("CurTempoOfs", m_engine->getCurTempoOfs());
dumpVar("TempoOfs1", m_engine->getTempoOfs1());
dumpVar("TempoOfs2", m_engine->getTempoOfs2());
dumpVar("TempoOfs3", m_engine->getTempoOfs3());
dumpVar("PunchEnabled", m_engine->getPunchEnabled());
dumpVar("PunchInTime", m_engine->getPunchInTime());
dumpVar("PunchOutTime", m_engine->getPunchOutTime());
dumpVar("EndAllTime", m_engine->getEndAllTime());
}
}
//_____________________________________________________________________________
void dumpCats(bool debug, TString fileName, TString dirName,
bool smearMassError) {
TFile* file = TFile::Open(fileName.Data());
TDirectory* theDir = (TDirectory*) file->FindObjectAny(dirName.Data());
TTree* theTree = (TTree*) theDir->Get("hPhotonTree");
UInt_t run, lumi, evt;
float rho, mass;
Int_t tth, vhLep, vhMet, vhHad, vbf, cat;
theTree->SetBranchAddress("run", &run);
theTree->SetBranchAddress("lumi",&lumi);
theTree->SetBranchAddress("evt", &evt);
theTree->SetBranchAddress("mass",&mass);
theTree->SetBranchAddress("rho",&rho);
theTree->SetBranchAddress("tthTag",&tth);
theTree->SetBranchAddress("VHLepTag",&vhLep);
theTree->SetBranchAddress("VHHadTag",&vhHad);
theTree->SetBranchAddress("vbfTag",&vbf);
float ph1e, ph1pt, ph1eerr, ph1eerrsmeared, teta1, phi1;
float ph2e, ph2pt, ph2eerr, ph2eerrsmeared, teta2, phi2;
theTree->SetBranchAddress("ph1.pt",&ph1pt);
theTree->SetBranchAddress("ph1.e",&ph1e);
theTree->SetBranchAddress("ph1.eerr",&ph1eerr);
theTree->SetBranchAddress("ph1.eerrsmeared",&ph1eerrsmeared);
theTree->SetBranchAddress("ph1.eta",&teta1);
theTree->SetBranchAddress("ph1.phi",&phi1);
theTree->SetBranchAddress("ph2.pt",&ph2pt);
theTree->SetBranchAddress("ph2.e",&ph2e);
theTree->SetBranchAddress("ph2.eerr",&ph2eerr);
theTree->SetBranchAddress("ph2.eerrsmeared",&ph2eerrsmeared);
theTree->SetBranchAddress("ph2.eta",&teta2);
theTree->SetBranchAddress("ph2.phi",&phi2);
Float_t ele1_pt, ele1_eta;
Float_t mu1_pt, mu1_eta;
theTree->SetBranchAddress("elePt" , &ele1_pt );
theTree->SetBranchAddress("eleEta" , &ele1_eta);
theTree->SetBranchAddress("muonPt" , &mu1_pt );
theTree->SetBranchAddress("muonEta", &mu1_eta );
float jet1pt, jet1eta, jet1phi;
float jet2pt, jet2eta, jet2phi;
theTree->SetBranchAddress("jet1pt",&jet1pt);
theTree->SetBranchAddress("jet1eta",&jet1eta);
theTree->SetBranchAddress("jet1phi",&jet1phi);
theTree->SetBranchAddress("jet2pt",&jet2pt);
theTree->SetBranchAddress("jet2eta",&jet2eta);
theTree->SetBranchAddress("jet2phi",&jet2phi);
float masserr, masserrwvtx, masserr_ns, masserrwvtx_ns, vtxprob, idmva_1, idmva_2;
theTree->SetBranchAddress("masserrsmeared",&masserr);
theTree->SetBranchAddress("masserrsmearedwrongvtx",&masserrwvtx);
theTree->SetBranchAddress("masserr",&masserr_ns);
theTree->SetBranchAddress("masserrwrongvtx",&masserrwvtx_ns);
theTree->SetBranchAddress("vtxprob",&vtxprob);
theTree->SetBranchAddress("ph1.idmva",&idmva_1);
theTree->SetBranchAddress("ph2.idmva",&idmva_2);
// MET tag stuff
float corrpfmet, corrpfmetphi, pfmet, pfmetphi;
theTree->SetBranchAddress("corrpfmet",&corrpfmet);
theTree->SetBranchAddress("corrpfmetphi",&corrpfmetphi);
theTree->SetBranchAddress("pfmet",&pfmet);
theTree->SetBranchAddress("pfmetphi",&pfmetphi);
float phigg, jetleadNoIDpt, jetleadNoIDphi, jetleadNoIDeta;
float ph1sceta, ph1scphi;
float ph2sceta, ph2scphi;
theTree->SetBranchAddress("phigg",&phigg);
theTree->SetBranchAddress("jetleadNoIDpt",&jetleadNoIDpt);
theTree->SetBranchAddress("jetleadNoIDphi",&jetleadNoIDphi);
theTree->SetBranchAddress("jetleadNoIDeta",&jetleadNoIDeta);
theTree->SetBranchAddress("ph1.sceta",&ph1sceta);
theTree->SetBranchAddress("ph1.scphi",&ph1scphi);
theTree->SetBranchAddress("ph2.sceta",&ph2sceta);
theTree->SetBranchAddress("ph2.scphi",&ph2scphi);
// Setup the diphoton BDT
Float_t rVtxSigmaMoM, wVtxSigmaMoM, cosDPhi;
Float_t pho1_ptOverM;
Float_t pho2_ptOverM;
Float_t diphoMVA;
TMVA::Reader* reader = new TMVA::Reader("Silent");
reader->AddVariable("masserrsmeared/mass" , &rVtxSigmaMoM);
reader->AddVariable("masserrsmearedwrongvtx/mass", &wVtxSigmaMoM);
reader->AddVariable("vtxprob" , &vtxprob );
reader->AddVariable("ph1.pt/mass" , &pho1_ptOverM);
reader->AddVariable("ph2.pt/mass" , &pho2_ptOverM);
reader->AddVariable("ph1.eta" , &teta1 );
reader->AddVariable("ph2.eta" , &teta2 );
reader->AddVariable("TMath::Cos(ph1.phi-ph2.phi)", &cosDPhi );
reader->AddVariable("ph1.idmva" , &idmva_1 );
reader->AddVariable("ph2.idmva" , &idmva_2 );
const char *diphotonWeights = (
"/home/veverka/cms/cmssw/031/CMSSW_5_3_10_patch1/src/MitPhysics/data/"
"HggBambu_SMDipho_Oct01_redqcdweightallsigevenbkg_BDTG.weights.xml"
);
reader->BookMVA("BDTG", diphotonWeights);
TRandom3 rng(0);
int eventCounter=0;
// Loop over the entries.
std::cout << "Looping over " << theTree->GetEntries() << " entries." << std::endl;
for (int i=0; i < theTree->GetEntries(); ++i) {
if (eventCounter > 9 && debug ) break;
if (debug) {
cout << "Processing entry " << i << " :" << endl
<< " mass: " << mass << endl
<< " ph1pt: " << ph1pt << endl
<< " ph2pt: " << ph2pt << endl
<< " idmva_1:" << idmva_1 << endl
<< " idmva_2:" << idmva_2 << endl;
}
theTree->GetEntry(i);
// MET category
vhMet = 0;
double dEtaJPh1 = ph1sceta - jetleadNoIDeta;
double dPhiJPh1 = TMath::ACos(TMath::Cos(ph1scphi - jetleadNoIDphi));
double dRJPh1 = TMath::Sqrt(TMath::Power(dEtaJPh1, 2) +
TMath::Power(dPhiJPh1, 2));
double dEtaJPh2 = ph2sceta - jetleadNoIDeta;
double dPhiJPh2 = TMath::ACos(TMath::Cos(ph2scphi - jetleadNoIDphi));
double dRJPh2 = TMath::Sqrt(TMath::Power(dEtaJPh2, 2) +
TMath::Power(dPhiJPh2, 2));
double dPhiMetGG = TMath::ACos(TMath::Cos(phigg - corrpfmetphi));
double dPhiMetJet = TMath::ACos(
TMath::Cos(TMath::Abs(jetleadNoIDphi - corrpfmetphi))
);
if (TMath::Abs(ph1sceta) < 1.4442 &&
TMath::Abs(ph2sceta) < 1.4442 &&
corrpfmet > 70. &&
ph1pt/mass > 45./120. &&
dPhiMetGG > 2.1 &&
(
jetleadNoIDpt < 50. ||
dRJPh1 < 0.5 ||
dRJPh2 < 0.5 ||
dPhiMetJet < 2.7
) &&
ph2pt > mass/4) {
vhMet = 1;
}
// Calculate needed variables for the diphoMVA
if (smearMassError) {
rVtxSigmaMoM = masserr / mass; // with smearing
wVtxSigmaMoM = masserrwvtx / mass; // with smearing
} else {
rVtxSigmaMoM = masserr_ns / mass; // no smearing
wVtxSigmaMoM = masserrwvtx_ns / mass; // no smearing
}
cosDPhi = TMath::Cos(phi1 - phi2);
pho1_ptOverM = ph1pt / mass;
pho2_ptOverM = ph2pt / mass;
diphoMVA = reader->EvaluateMVA("BDTG");
bool passPreselection = (mass > 100 &&
mass < 180 &&
ph1pt > mass/3 &&
ph2pt > mass/4 &&
idmva_1 > -0.2 &&
idmva_2 > -0.2);
if (passPreselection == false) {
if (debug) {
cout << " passPreselection: " << passPreselection << endl;
}
continue;
}
if (debug) {
cout << " ... passed preselection." << endl;
}
eventCounter++;
if (tth == 1) tth = 2;
else if (tth == 2) tth = 1;
if (vhHad == 2) vhHad = 1;
cat = kIncl0;
if (tth == 2) cat = kTTHLep;
else if (vhLep == 2) cat = kVHLepTight;
else if (vhLep == 1) cat = kVHLepLoose;
else if (vbf > 0) cat = kDijet0;
else if (vhMet == 1) cat = kVHMet;
else if (tth == 1) cat = kTTHHad;
else if (vhHad == 1) cat = kVHHad;
// if (cat == kIncl0 && diphoMVA < -0.4) continue;
// Event Variables
dumpVar("run" , run ); // 1
dumpVar("lumi" , lumi ); // 2
dumpVar("event" , evt ); // 3
dumpVar("cat" , cat );
dumpVar("tth" , tth );
dumpVar("vhLep" , vhLep );
dumpVar("vhMet" , vhMet );
dumpVar("vhHad" , vhHad );
dumpVar("vbf" , vbf );
// Leading Photon Variables
dumpVar("pho1_e" , ph1e ); // 10
dumpVar("pho1_eErr" , ph1eerr ); // 11
dumpVar("pho1_eta" , teta1 ); // 8
dumpVar("pho1_phi" , phi1 ); // 9
dumpVar("pho1_idMVA" , idmva_1 );
// Trailing Photon Variables
dumpVar("pho2_e" , ph2e ); // 36
dumpVar("pho2_eErr" , ph2eerr ); // 37
dumpVar("pho2_eta" , teta2 ); // 34
dumpVar("pho2_phi" , phi2 ); // 35
dumpVar("pho2_idMVA" , idmva_2 );
// Diphoton Variables
dumpVar("mass" , mass );
dumpVar("met" , corrpfmet );
dumpVar("met_phi" , corrpfmetphi );
dumpVar("uncorrMet" , pfmet );
dumpVar("uncorrMet_phi" , pfmetphi );
dumpVar("diphoMVA" , diphoMVA );
// Muon Variables
if (mu1_pt < 0) {
mu1_pt = -999;
mu1_eta = -999;
}
dumpVar("mu1_pt" , mu1_pt );
dumpVar("mu1_eta" , mu1_eta );
// Electron Variables
if (ele1_pt < 0) {
ele1_pt = -999;
ele1_eta = -999;
}
dumpVar("ele1_pt" , ele1_pt );
dumpVar("ele1_eta" , ele1_eta );
// Leading Jet Variables
if (jet1pt < 0) {
jet1pt = -999;
jet1eta = -999;
}
dumpVar("jet1_pt" , jet1pt ); // 69
dumpVar("jet1_eta" , jet1eta ); // 70
dumpVar("jet1_phi" , jet1phi ); // 70
// Trailing Jet Variables
if (jet2pt < 0) {
jet2pt = -999;
jet2eta = -999;
}
dumpVar("jet2_pt" , jet2pt ); // 72
dumpVar("jet2_eta" , jet2eta ); // 73
dumpVar("jet2_phi" , jet2phi ); // 70
std::cout << std::endl;
} // Loop over the tree entries.
return;
} // void dumpMvaInputs(bool debug, TString fileName)
//_____________________________________________________________________________
void dumpMvaInputs(bool debug, TString fileName) {
TFile* file = TFile::Open(fileName.Data());
const char *treeName = "PhotonTreeWriterPresel";
// const char *treeName = "PhotonTreeWriterPreselNoSmear";
TDirectory* theDir = (TDirectory*) file->FindObjectAny(treeName);
TTree* theTree = (TTree*) theDir->Get("hPhotonTree");
// open the MVA files, if requested
UInt_t run, lumi, evt;
UInt_t ph1index, ph1scindex;
UInt_t ph2index, ph2scindex;
float ph1pt, ph1sceta, ph1iso1, ph1iso2, ph1iso3, ph1cov, ph1hoe, ph1r9;
float ph2pt, ph2sceta, ph2iso1, ph2iso2, ph2iso3, ph2cov, ph2hoe, ph2r9;
float ph1e, ph2e;
float rho, mass;
float ph1ecalIso3,ph1ecalIso4,ph1trackIsoSel03,ph1trackIsoWorst04;
float ph2ecalIso3,ph2ecalIso4,ph2trackIsoSel03,ph2trackIsoWorst04;
float ph1eerr, ph1eerrsmeared;
float ph2eerr, ph2eerrsmeared;
UChar_t ph1hasconversion;
UChar_t ph2hasconversion;
Float_t scetawidth1, scphiwidth1;
Float_t scetawidth2, scphiwidth2;
// 2012 id mva
Float_t ph1_idmva_CoviEtaiPhi;
Float_t ph1_idmva_s4ratio;
Float_t ph1_idmva_GammaIso;
Float_t ph1_idmva_ChargedIso_selvtx;
Float_t ph1_idmva_ChargedIso_0p2_selvtx;
Float_t ph1_idmva_ChargedIso_worstvtx;
Float_t ph1_idmva_PsEffWidthSigmaRR;
Float_t ph2_idmva_CoviEtaiPhi;
Float_t ph2_idmva_s4ratio;
Float_t ph2_idmva_GammaIso;
Float_t ph2_idmva_ChargedIso_selvtx;
Float_t ph2_idmva_ChargedIso_0p2_selvtx;
Float_t ph2_idmva_ChargedIso_worstvtx;
Float_t ph2_idmva_PsEffWidthSigmaRR;
theTree->SetBranchAddress("run", &run);
theTree->SetBranchAddress("lumi",&lumi);
theTree->SetBranchAddress("evt", &evt);
theTree->SetBranchAddress("mass",&mass);
theTree->SetBranchAddress("rho",&rho);
theTree->SetBranchAddress("ph1.index",&ph1index);
theTree->SetBranchAddress("ph1.scindex",&ph1scindex);
theTree->SetBranchAddress("ph1.pt",&ph1pt);
theTree->SetBranchAddress("ph1.e",&ph1e);
theTree->SetBranchAddress("ph1.eerr",&ph1eerr);
theTree->SetBranchAddress("ph1.eerrsmeared",&ph1eerrsmeared);
theTree->SetBranchAddress("ph1.sceta",&ph1sceta);
theTree->SetBranchAddress("ph1.pfcic4_tIso1",&ph1iso1);
theTree->SetBranchAddress("ph1.pfcic4_tIso2",&ph1iso2);
theTree->SetBranchAddress("ph1.pfcic4_tIso3",&ph1iso3);
theTree->SetBranchAddress("ph1.pfcic4_covIEtaIEta",&ph1cov);
theTree->SetBranchAddress("ph1.pfcic4_HoE",&ph1hoe);
theTree->SetBranchAddress("ph1.pfcic4_R9",&ph1r9);
theTree->SetBranchAddress("ph1.pfcic4_ecalIso3",&ph1ecalIso3);
theTree->SetBranchAddress("ph1.pfcic4_ecalIso4",&ph1ecalIso4);
theTree->SetBranchAddress("ph1.pfcic4_trackIsoSel03",&ph1trackIsoSel03);
theTree->SetBranchAddress("ph1.pfcic4_trackIsoWorst04",&ph1trackIsoWorst04);
theTree->SetBranchAddress("ph1.hasconversion",&ph1hasconversion);
theTree->SetBranchAddress("ph1.scetawidth",&scetawidth1);
theTree->SetBranchAddress("ph1.scphiwidth",&scphiwidth1);
theTree->SetBranchAddress("ph1.idmva_CoviEtaiPhi",&ph1_idmva_CoviEtaiPhi);
theTree->SetBranchAddress("ph1.idmva_s4ratio",&ph1_idmva_s4ratio);
theTree->SetBranchAddress("ph1.idmva_GammaIso",&ph1_idmva_GammaIso);
theTree->SetBranchAddress("ph1.idmva_ChargedIso_selvtx",
&ph1_idmva_ChargedIso_selvtx);
theTree->SetBranchAddress("ph1.idmva_ChargedIso_0p2_selvtx",
&ph1_idmva_ChargedIso_0p2_selvtx);
theTree->SetBranchAddress("ph1.idmva_ChargedIso_worstvtx",
&ph1_idmva_ChargedIso_worstvtx);
theTree->SetBranchAddress("ph1.idmva_PsEffWidthSigmaRR",
&ph1_idmva_PsEffWidthSigmaRR);
theTree->SetBranchAddress("ph2.index",&ph2index);
theTree->SetBranchAddress("ph2.scindex",&ph2scindex);
theTree->SetBranchAddress("ph2.pt",&ph2pt);
theTree->SetBranchAddress("ph2.e",&ph2e);
theTree->SetBranchAddress("ph2.eerr",&ph2eerr);
theTree->SetBranchAddress("ph2.eerrsmeared",&ph2eerrsmeared);
theTree->SetBranchAddress("ph2.sceta",&ph2sceta);
theTree->SetBranchAddress("ph2.pfcic4_tIso1",&ph2iso1);
theTree->SetBranchAddress("ph2.pfcic4_tIso2",&ph2iso2);
theTree->SetBranchAddress("ph2.pfcic4_tIso3",&ph2iso3);
theTree->SetBranchAddress("ph2.pfcic4_covIEtaIEta",&ph2cov);
theTree->SetBranchAddress("ph2.pfcic4_HoE",&ph2hoe);
theTree->SetBranchAddress("ph2.pfcic4_R9",&ph2r9);
theTree->SetBranchAddress("ph2.pfcic4_ecalIso3",&ph2ecalIso3);
theTree->SetBranchAddress("ph2.pfcic4_ecalIso4",&ph2ecalIso4);
theTree->SetBranchAddress("ph2.pfcic4_trackIsoSel03",&ph2trackIsoSel03);
theTree->SetBranchAddress("ph2.pfcic4_trackIsoWorst04",&ph2trackIsoWorst04);
theTree->SetBranchAddress("ph2.hasconversion",&ph2hasconversion);
theTree->SetBranchAddress("ph2.scetawidth",&scetawidth2);
theTree->SetBranchAddress("ph2.scphiwidth",&scphiwidth2);
theTree->SetBranchAddress("ph2.idmva_CoviEtaiPhi",&ph2_idmva_CoviEtaiPhi);
theTree->SetBranchAddress("ph2.idmva_s4ratio",&ph2_idmva_s4ratio);
theTree->SetBranchAddress("ph2.idmva_GammaIso",&ph2_idmva_GammaIso);
theTree->SetBranchAddress("ph2.idmva_ChargedIso_selvtx",
&ph2_idmva_ChargedIso_selvtx);
theTree->SetBranchAddress("ph2.idmva_ChargedIso_0p2_selvtx",
&ph2_idmva_ChargedIso_0p2_selvtx);
theTree->SetBranchAddress("ph2.idmva_ChargedIso_worstvtx",
&ph2_idmva_ChargedIso_worstvtx);
theTree->SetBranchAddress("ph2.idmva_PsEffWidthSigmaRR",
&ph2_idmva_PsEffWidthSigmaRR);
float jet1pt, jet2pt, jet1eta, jet2eta, dijetmass, zeppenfeld, dphidijetgg;
theTree->SetBranchAddress("jet1pt",&jet1pt);
theTree->SetBranchAddress("jet2pt",&jet2pt);
theTree->SetBranchAddress("jet1eta",&jet1eta);
theTree->SetBranchAddress("jet2eta",&jet2eta);
theTree->SetBranchAddress("dijetmass",&dijetmass);
theTree->SetBranchAddress("zeppenfeld",&zeppenfeld);
theTree->SetBranchAddress("dphidijetgg",&dphidijetgg);
int numVtx;
theTree->SetBranchAddress("nVtx",&numVtx);
// presel vars
float ecalisodr03_1,ecalisodr03_2;
theTree->SetBranchAddress("ph1.ecalisodr03",&ecalisodr03_1);
theTree->SetBranchAddress("ph2.ecalisodr03",&ecalisodr03_2);
float hcalisodr03_1,hcalisodr03_2;
theTree->SetBranchAddress("ph1.hcalisodr03",&hcalisodr03_1);
theTree->SetBranchAddress("ph2.hcalisodr03",&hcalisodr03_2);
float trkisohollowdr03_1,trkisohollowdr03_2;
theTree->SetBranchAddress("ph1.trkisohollowdr03",&trkisohollowdr03_1);
theTree->SetBranchAddress("ph2.trkisohollowdr03",&trkisohollowdr03_2);
float hoveretower_1, hoveretower_2;
theTree->SetBranchAddress("ph1.hoveretower", &hoveretower_1);
theTree->SetBranchAddress("ph2.hoveretower", &hoveretower_2);
float sieie_1,sieie_2;
theTree->SetBranchAddress("ph1.sigietaieta",&sieie_1);
theTree->SetBranchAddress("ph2.sigietaieta",&sieie_2);
float idmva_ChargedIso_presel_1,idmva_ChargedIso_presel_2;
theTree->SetBranchAddress("ph1.idmva_ChargedIso_selvtx",&idmva_ChargedIso_presel_1);
theTree->SetBranchAddress("ph2.idmva_ChargedIso_selvtx",&idmva_ChargedIso_presel_2);
float idmva_ChargedIso_0p2_presel_1,idmva_ChargedIso_0p2_presel_2;
theTree->SetBranchAddress("ph1.idmva_ChargedIso_0p2_selvtx",&idmva_ChargedIso_0p2_presel_1);
theTree->SetBranchAddress("ph2.idmva_ChargedIso_0p2_selvtx",&idmva_ChargedIso_0p2_presel_2);
// float eleVeto_1, eleVeto_2;
// theTree->SetBranchAddress("ph1.eleVeto",&eleVeto_1);
// theTree->SetBranchAddress("ph2.eleVeto",&eleVeto_2);
float idmva_1, idmva_2;
theTree->SetBranchAddress("ph1.idmva",&idmva_1);
theTree->SetBranchAddress("ph2.idmva",&idmva_2);
float teta1, teta2;
theTree->SetBranchAddress("ph1.eta",&teta1);
theTree->SetBranchAddress("ph2.eta",&teta2);
float vtxprob, ptgg, phi1, phi2, masserr, masserrwvtx, masserr_ns, masserrwvtx_ns;
theTree->SetBranchAddress("vtxprob",&vtxprob);
theTree->SetBranchAddress("ptgg",&ptgg);
theTree->SetBranchAddress("ph1.phi",&phi1);
theTree->SetBranchAddress("ph2.phi",&phi2);
theTree->SetBranchAddress("masserrsmeared",&masserr);
theTree->SetBranchAddress("masserrsmearedwrongvtx",&masserrwvtx);
theTree->SetBranchAddress("masserr",&masserr_ns);
theTree->SetBranchAddress("masserrwrongvtx",&masserrwvtx_ns);
int vbfTag, leptonTag;
theTree->SetBranchAddress("vbfTag",&vbfTag);
theTree->SetBranchAddress("leptonTag",&leptonTag);
float corrpfmet, corrpfmetphi, pfmet, pfmetphi;
theTree->SetBranchAddress("corrpfmet",&corrpfmet);
theTree->SetBranchAddress("corrpfmetphi",&corrpfmetphi);
theTree->SetBranchAddress("pfmet",&pfmet);
theTree->SetBranchAddress("pfmetphi",&pfmetphi);
float mupt, mueta, mudr1, mudr2, mudz, mud0;
theTree->SetBranchAddress("muonPt",&mupt);
theTree->SetBranchAddress("muonEta",&mueta);
theTree->SetBranchAddress("muDR1",&mudr1);
theTree->SetBranchAddress("muDR2",&mudr2);
theTree->SetBranchAddress("muD0",&mud0);
theTree->SetBranchAddress("muDZ",&mudz);
float elept, eleeta, elesceta, eledr1, eledr2, eleMass1, eleMass2;
theTree->SetBranchAddress("elePt",&elept);
theTree->SetBranchAddress("eleEta",&eleeta);
theTree->SetBranchAddress("eleSCEta",&elesceta);
theTree->SetBranchAddress("eleDR1",&eledr1);
theTree->SetBranchAddress("eleDR2",&eledr2);
theTree->SetBranchAddress("eleMass1",&eleMass1);
theTree->SetBranchAddress("eleMass2",&eleMass2);
int elemisshits;
theTree->SetBranchAddress("eleNinnerHits",&elemisshits);
// vertext stuff
int vertexId1, vertexId2, vertexId3;
float vertexMva1, vertexMva2, vertexMva3;
//float ptbal, ptasym, sumpt2, p2conv;
Float_t ptbal, ptasym, sumpt2, p2conv;
int nleg1, nleg2;
theTree->SetBranchAddress("vtxInd1",&vertexId1);
theTree->SetBranchAddress("vtxInd2",&vertexId2);
theTree->SetBranchAddress("vtxInd3",&vertexId3);
theTree->SetBranchAddress("vtxMva1",&vertexMva1);
theTree->SetBranchAddress("vtxMva2",&vertexMva2);
theTree->SetBranchAddress("vtxMva3",&vertexMva3);
theTree->SetBranchAddress("vtxBestPtbal",&ptbal);
theTree->SetBranchAddress("vtxBestPtasym",&ptasym);
theTree->SetBranchAddress("vtxBestSumpt2",&sumpt2);
theTree->SetBranchAddress("vtxBestP2Conv",&p2conv);
theTree->SetBranchAddress("vtxNleg1",&nleg1);
theTree->SetBranchAddress("vtxNleg2",&nleg2);
float convVtxZ1, convVtxRes1, convVtxChi1;
float convVtxZ2, convVtxRes2, convVtxChi2;
int convVtxIdx1, convVtxIdx2, vtxNconv;
theTree->SetBranchAddress("vtxConv1Z",&convVtxZ1);
theTree->SetBranchAddress("vtxConv1DZ",&convVtxRes1);
theTree->SetBranchAddress("vtxConv1Prob",&convVtxChi1);
theTree->SetBranchAddress("vtxConvIdx1",&convVtxIdx1);
theTree->SetBranchAddress("vtxConv2Z",&convVtxZ2);
theTree->SetBranchAddress("vtxConv2DZ",&convVtxRes2);
theTree->SetBranchAddress("vtxConv2Prob",&convVtxChi2);
theTree->SetBranchAddress("vtxConvIdx2",&convVtxIdx2);
theTree->SetBranchAddress("vtxNconv",&vtxNconv);
float vtxz1, vtxz2, vtxz3;
theTree->SetBranchAddress("vtxMva1Z",&vtxz1);
theTree->SetBranchAddress("vtxMva2Z",&vtxz2);
theTree->SetBranchAddress("vtxMva3Z",&vtxz3);
float eleIdMva;
theTree->SetBranchAddress("eleIdMva",&eleIdMva);
// additional MET tag stuff.
float phigg, jetleadNoIDpt, jetleadNoIDphi, jetleadNoIDeta;
float ph1scphi, ph2scphi;
theTree->SetBranchAddress("phigg",&phigg);
theTree->SetBranchAddress("jetleadNoIDpt",&jetleadNoIDpt);
theTree->SetBranchAddress("jetleadNoIDphi",&jetleadNoIDphi);
theTree->SetBranchAddress("jetleadNoIDeta",&jetleadNoIDeta);
theTree->SetBranchAddress("ph1.scphi",&ph1scphi);
theTree->SetBranchAddress("ph2.scphi",&ph2scphi);
float ph1scrawe, ph1scpse, ph1e3x3, ph1e5x5, ph1e3x3seed, ph1e5x5seed;
float ph2scrawe, ph2scpse, ph2e3x3, ph2e5x5, ph2e3x3seed, ph2e5x5seed;
theTree->SetBranchAddress("ph1.scrawe",&ph1scrawe);
theTree->SetBranchAddress("ph1.scpse",&ph1scpse);
theTree->SetBranchAddress("ph1.e3x3",&ph1e3x3);
theTree->SetBranchAddress("ph1.e3x3seed",&ph1e3x3seed);
theTree->SetBranchAddress("ph1.e5x5",&ph1e5x5);
theTree->SetBranchAddress("ph1.e5x5seed",&ph1e5x5seed);
theTree->SetBranchAddress("ph2.scrawe",&ph2scrawe);
theTree->SetBranchAddress("ph2.scpse",&ph2scpse);
theTree->SetBranchAddress("ph2.e3x3",&ph2e3x3);
theTree->SetBranchAddress("ph2.e3x3seed",&ph2e3x3seed);
theTree->SetBranchAddress("ph2.e5x5",&ph2e5x5);
theTree->SetBranchAddress("ph2.e5x5seed",&ph2e5x5seed);
// Setup the diphoton BDT
Float_t rVtxSigmaMoM, wVtxSigmaMoM, cosDPhi;
Float_t pho1_ptOverM;
Float_t pho2_ptOverM;
Float_t diphoMVA;
TMVA::Reader* reader = new TMVA::Reader();
reader->AddVariable("masserrsmeared/mass" , &rVtxSigmaMoM);
reader->AddVariable("masserrsmearedwrongvtx/mass", &wVtxSigmaMoM);
reader->AddVariable("vtxprob" , &vtxprob );
reader->AddVariable("ph1.pt/mass" , &pho1_ptOverM);
reader->AddVariable("ph2.pt/mass" , &pho2_ptOverM);
reader->AddVariable("ph1.eta" , &teta1 );
reader->AddVariable("ph2.eta" , &teta2 );
reader->AddVariable("TMath::Cos(ph1.phi-ph2.phi)", &cosDPhi );
reader->AddVariable("ph1.idmva" , &idmva_1 );
reader->AddVariable("ph2.idmva" , &idmva_2 );
const char *diphotonWeights = (
"/home/veverka/cms/cmssw/031/CMSSW_5_3_10_patch1/src/MitPhysics/data/"
"HggBambu_SMDipho_Oct01_redqcdweightallsigevenbkg_BDTG.weights.xml"
);
reader->BookMVA("BDTG", diphotonWeights);
TRandom3 rng(0);
int eventCounter=0;
// Loop over the entries.
std::cout << "Looping over " << theTree->GetEntries() << " entries." << std::endl;
for (int i=0; i < theTree->GetEntries(); ++i) {
if (eventCounter > 9 && debug ) break;
theTree->GetEntry(i);
bool passPreselection = (mass > 100 &&
mass < 180 &&
ph1pt > mass/3 &&
ph2pt > mass/4);
if (passPreselection == false) continue;
eventCounter++;
// Calculate needed variables
// rVtxSigmaMoM = masserr_ns / mass; // no smearing
// wVtxSigmaMoM = masserrwvtx_ns / mass; // no smearing
rVtxSigmaMoM = masserr / mass; // with smearing
wVtxSigmaMoM = masserrwvtx / mass; // with smearing
cosDPhi = TMath::Cos(phi1 - phi2);
pho1_ptOverM = ph1pt / mass;
pho2_ptOverM = ph2pt / mass;
diphoMVA = reader->EvaluateMVA("BDTG");
// Event Variables
dumpVar("run" , run ); // 1
dumpVar("lumi" , lumi ); // 2
dumpVar("event" , evt ); // 3
dumpVar("rho" , rho ); // 4
// Leading Photon Variables
dumpVar("pho1_ind" , ph1index ); // 5
dumpVar("pho1_scInd" , /*ph1scindex*/ -999 ); // 6
dumpVar("pho1_pt" , ph1pt ); // 7
dumpVar("pho1_eta" , teta1 ); // 8
dumpVar("pho1_phi" , phi1 ); // 9
dumpVar("pho1_e" , ph1e ); // 10
dumpVar("pho1_eErr" , ph1eerr ); // 11
dumpVar("pho1_isConv" ,
(UInt_t) ph1hasconversion ); // 12
dumpVar("pho1_HoE" , hoveretower_1 ); // 13
dumpVar("pho1_hcalIso03" ,
hcalisodr03_1 - 0.005 * ph1pt ); // 14
dumpVar("pho1_trkIso03" ,
trkisohollowdr03_1 - 0.002 * ph1pt ); // 15
dumpVar("pho1_pfChargedIsoGood02",
ph1_idmva_ChargedIso_0p2_selvtx ); // 16
dumpVar("pho1_pfChargedIsoGood03",
ph1_idmva_ChargedIso_selvtx ); // 17
dumpVar("pho1_pfChargedIsoBad03" ,
ph1_idmva_ChargedIso_worstvtx ); // 18
dumpVar("pho1_pfPhotonIso03" , ph1_idmva_GammaIso ); // 19
// TODO: remove pho1_pfNeutralIso03, it's not used
// dumpVar("pho1_pfNeutralIso03" , -999 ); // 20
dumpVar("pho1_sieie" , sieie_1 ); // 21
dumpVar("pho1_cieip" , ph1_idmva_CoviEtaiPhi ); // 22
dumpVar("pho1_etaWidth" , scetawidth1 ); // 23
dumpVar("pho1_phiWidth" , scphiwidth1 ); // 24
dumpVar("pho1_r9" , ph1r9 ); // 25
// TODO: remove pho1_lambdaRatio, it's not used
// dumpVar("pho1_lambdaRatio" , -999 ); // 26
dumpVar("pho1_s4Ratio" , ph1_idmva_s4ratio ); // 27
dumpVar("pho1_scEta" , ph1sceta ); // 28
dumpVar("pho1_ESEffSigmaRR" ,
ph1_idmva_PsEffWidthSigmaRR ); // 29
dumpVar("pho1_ptOverM" , pho1_ptOverM ); // 30
dumpVar("pho1_scRawE" , ph1scrawe );
dumpVar("pho1_idMVA" , idmva_1 );
// Trailing Photon Variables
dumpVar("pho2_ind" , ph2index ); // 31
dumpVar("pho2_scInd" , /*ph2scindex*/ -999 ); // 32
dumpVar("pho2_pt" , ph2pt ); // 33
dumpVar("pho2_eta" , teta2 ); // 34
dumpVar("pho2_phi" , phi2 ); // 35
dumpVar("pho2_e" , ph2e ); // 36
dumpVar("pho2_eErr" , ph2eerr ); // 37
dumpVar("pho2_isConv" ,
(UInt_t) ph2hasconversion ); // 38
dumpVar("pho2_HoE" , hoveretower_2 ); // 39
dumpVar("pho2_hcalIso03" ,
hcalisodr03_1 - 0.005 * ph2pt ); // 40
dumpVar("pho2_trkIso03" ,
trkisohollowdr03_1 - 0.002 * ph2pt ); // 41
dumpVar("pho2_pfChargedIsoGood02",
ph2_idmva_ChargedIso_0p2_selvtx ); // 42
dumpVar("pho2_pfChargedIsoGood03",
ph2_idmva_ChargedIso_selvtx ); // 43
dumpVar("pho2_pfChargedIsoBad03" ,
ph2_idmva_ChargedIso_worstvtx ); // 44
dumpVar("pho2_pfPhotonIso03" , ph2_idmva_GammaIso ); // 45
// dumpVar("pho2_pfNeutralIso03" , -999 ); // 46
dumpVar("pho2_sieie" , sieie_2 ); // 47
dumpVar("pho2_cieip" , ph2_idmva_CoviEtaiPhi ); // 48
dumpVar("pho2_etaWidth" , scetawidth2 ); // 49
dumpVar("pho2_phiWidth" , scphiwidth2 ); // 50
dumpVar("pho2_r9" , ph2r9 ); // 51
// dumpVar("pho2_lambdaRatio" , -999 ); // 52
dumpVar("pho2_s4Ratio" , ph2_idmva_s4ratio ); // 53
dumpVar("pho2_scEta" , ph2sceta ); // 54
dumpVar("pho2_ESEffSigmaRR" ,
ph2_idmva_PsEffWidthSigmaRR ); // 55
dumpVar("pho2_ptOverM" , pho2_ptOverM ); // 56
dumpVar("pho2_scRawE" , ph2scrawe );
dumpVar("pho2_idMVA" , idmva_2 );
// Diphoton Variables
dumpVar("mass" , mass ); // 57
dumpVar("rVtxSigmaMoM" , rVtxSigmaMoM ); // 58
dumpVar("wVtxSigmaMoM" , wVtxSigmaMoM ); // 59
dumpVar("vtxProb" , vtxprob ); // 61
float logSPt2 = TMath::Log(sumpt2);
if (vertexId1 < 0) {
vertexId1 = ptbal = ptasym = logSPt2 = p2conv = vtxNconv = -999;
}
dumpVar("vtxIndex" , vertexId1 ); // 60
dumpVar("ptBal" , ptbal ); // 62
dumpVar("ptAsym" , ptasym ); // 63
dumpVar("logSPt2" , logSPt2 ); // 64
dumpVar("p2Conv" , p2conv ); // 65
dumpVar("nConv" , vtxNconv ); // 66
dumpVar("cosDPhi" , cosDPhi ); // 67
dumpVar("diphoMVA" , diphoMVA );
// Leading Jet Variables
if (jet1pt < 0) {
jet1pt = -999;
jet1eta = -999;
}
dumpVar("jet1_ind" , -999 ); // 68
dumpVar("jet1_pt" , jet1pt ); // 69
dumpVar("jet1_eta" , jet1eta ); // 70
// Trailing Jet Variables
if (jet2pt < 0) {
jet2pt = -999;
jet2eta = -999;
}
dumpVar("jet2_ind" , -999 ); // 71
dumpVar("jet2_pt" , jet2pt ); // 72
dumpVar("jet2_eta" , jet2eta ); // 73
// Dijet Variables
float dijet_dEta = abs(jet1eta - jet2eta);
if (jet1pt < 0 || jet2pt < 0) {
dijet_dEta = -999;
zeppenfeld = -999;
dphidijetgg = -999;
dijetmass = -999;
}
dumpVar("dijet_dEta" , dijet_dEta ); // 74
dumpVar("dijet_Zep" , zeppenfeld ); // 75
dumpVar("dijet_dPhi" , dphidijetgg ); // 76
dumpVar("dijet_mass" , dijetmass , false); // 77
std::cout << std::endl;
} // Loop over the tree entries.
return;
} // void dumpMvaInputs(bool debug, TString fileName)
