void EnumDataKey(int nTabLevel, ISpDataKey *pSpDataKey) {
HRESULT hr = S_OK;
// enumerate subkeys recursively
for (ULONG k = 0; ; k++) {
LPWSTR key = nullptr;
hr = pSpDataKey->EnumKeys(k, &key);
if (SPERR_NO_MORE_ITEMS == hr) {
// done
break;
} else if (FAILED(hr)) {
ERR(L"ISpDataKey::EnumKeys failed: hr = 0x%08x", hr);
continue;
}
CoTaskMemFreeOnExit fKey(key);
Tab(nTabLevel);
LOG(L"%s", key);
ISpDataKey *pSubKey = nullptr;
hr = pSpDataKey->OpenKey(key, &pSubKey);
if (FAILED(hr)) {
ERR(L"ISpDataKey::OpenKeys failed: hr = 0x%08x", hr);
continue;
}
ReleaseOnExit rSubKey(pSubKey);
EnumDataKey(nTabLevel + 1, pSubKey);
}
// enumerate values
for (ULONG v = 0; ; v++) {
LPWSTR val = nullptr;
hr = pSpDataKey->EnumValues(v, &val);
if (SPERR_NO_MORE_ITEMS == hr) {
// done
break;
} else if (FAILED(hr)) {
ERR(L"ISpDataKey::EnumKeys failed: hr = 0x%08x", hr);
continue;
}
CoTaskMemFreeOnExit fVal(val);
// how do we know whether it's a string or a DWORD?
LPWSTR data = nullptr;
hr = pSpDataKey->GetStringValue(val, &data);
if (FAILED(hr)) {
ERR(L"ISpDataKey::GetStringValue failed: hr = 0x%08x", hr);
continue;
}
CoTaskMemFreeOnExit fData(data);
Tab(nTabLevel);
if (0 == wcscmp(val, L"")) {
LOG(L"(default) = %s", data);
} else {
LOG(L"%s = %s", val, data);
}
}
}
AP_MSG_HANDLER_METHOD(NetSimModule, NetSim_ContinueHttp)
{
int ok = 1;
Http* pHttp = 0;
if (http_.Get(pMsg->hClient, pHttp)) {
try {
String sIndexFile = Apollo::getModuleConfig(MODULE_NAME, "HTTP/IndexFile", "");
if (!sIndexFile) { throw ApException(LOG_CONTEXT, "No index file name"); }
xFile fIndex(sIndexFile);
if (!fIndex.Load()) { throw ApException(LOG_CONTEXT, "f.Load(%s) failed", _sz(sIndexFile)); }
String sIndex;
fIndex.GetData(sIndex);
List lDataFiles;
KeyValueBlob2List(sIndex, lDataFiles, "\r\n", " ", "");
Elem* e = lDataFiles.FindByName(pHttp->sUrl_);
if (!e) { throw ApException(LOG_CONTEXT, "No data file for url=%s", _sz(pHttp->sUrl_)); }
String sDataFilename = e->getString();
if (!sDataFilename) { throw ApException(LOG_CONTEXT, "No data file for url=%s", _sz(pHttp->sUrl_)); }
String sDataFilepath = Apollo::getModuleConfig(MODULE_NAME, "HTTP/DataFilesBase", "");
sDataFilepath += String::filenamePathSeparator();
sDataFilepath += sDataFilename;
xFile fData(sDataFilepath);
if (!fData.Load()) { throw ApException(LOG_CONTEXT, "f.Load(%s) failed", _sz(sDataFilepath)); }
String sHeaderFilename = sDataFilepath;
sHeaderFilename += ".http";
xFile fHeader(sHeaderFilename);
if (!fHeader.Load()) { throw ApException(LOG_CONTEXT, "f.Load(%s) failed", _sz(sHeaderFilename)); }
{
ApAsyncMessage<Msg_Net_HTTP_Connected> msg;
msg->hClient = pHttp->hClient_;
msg.Post();
}
Apollo::Sleep(50);
String sHeader;
fHeader.GetData(sHeader);
{
ApAsyncMessage<Msg_Net_HTTP_Header> msg;
msg->hClient = pHttp->hClient_;
msg->nStatus = 200;
List lHeader;
KeyValueLfBlob2List(sHeader, lHeader);
msg->kvHeader.fromList(lHeader);
msg.Post();
}
Apollo::Sleep(50);
Buffer sbData;
fData.GetData(sbData);
size_t nSize = 0;
while (nSize < sbData.Length()) {
int nFragment = sbData.Length() - nSize;
if (nFragment > 50000) {
nFragment = 50000;
}
{
ApAsyncMessage<Msg_Net_HTTP_DataIn> msg;
msg->hClient = pHttp->hClient_;
msg->sbData.SetData(sbData.Data() + nSize, nFragment);
msg.Post();
}
Apollo::Sleep(50);
nSize += nFragment;
}
} catch (ApException& ex) {
apLog_Error((LOG_CHANNEL, LOG_CONTEXT, "%s url=%s", _sz(ex.getText()), _sz(pHttp->sUrl_)));
ApAsyncMessage<Msg_Net_HTTP_Failed> msg;
msg->hClient = pHttp->hClient_;
msg->sMessage = ex.getText();
msg.Post();
ok = 0;
}
{
ApAsyncMessage<Msg_Net_HTTP_Closed> msg;
msg->hClient = pHttp->hClient_;
msg.Post();
}
if (http_.Unset(pMsg->hClient)) {
delete pHttp;
pHttp = 0;
}
}
pMsg->apStatus = ok ? ApMessage::Ok : ApMessage::Error;
}
void fakeStudyMu( TCut Cuts_ = "ptL1>0",
TString variable_ = "MtLeg1",
TString XTitle_ = "M_{T}(#mu,MET)",
TString Unities_ = "GeV",
Int_t nBins_ = 12, Float_t xMin_=0, Float_t xMax_=120,
Float_t magnifySgn_ = 10,
Float_t hltEff_ = 1.0,
Int_t enableHLTmatching_ = 1,
Int_t logy_ = 0 )
{
float Lumi = 24.86+159.15;
TCanvas *c1 = new TCanvas("c1","",5,30,650,600);
c1->SetGrid(0,0);
c1->SetFillStyle(4000);
c1->SetFillColor(10);
c1->SetTicky();
c1->SetObjectStat(0);
c1->SetLogy(logy_);
TLegend* leg = new TLegend(0.60,0.50,0.85,0.85,NULL,"brNDC");
leg->SetFillStyle(0);
leg->SetBorderSize(0);
leg->SetFillColor(10);
leg->SetTextSize(0.04);
leg->SetHeader("#mu+#tau_{had}");
THStack* aStack = new THStack("aStack",Form("CMS Preliminary 2011 #sqrt{s}=7 TeV L=%.0f pb^{-1}", Lumi ));
TH1F* hSiml = new TH1F();
TH1F* hSgn = new TH1F();
TH1F* hSgn1 = new TH1F();
TH1F* hSgn2 = new TH1F();
TH1F* hData = new TH1F();
// OpenNTuples
TString fDataName = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleRun2011-Mu_Open_MuTauStream.root";
TString fSignalNameVBF = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleVBFH115-Mu-powheg-PUS1_Open_MuTauStream.root";
TString fSignalNameGGH = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleGGFH115-Mu-powheg-PUS1_Open_MuTauStream.root";
TString fBackgroundNameDYTauTau = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleDYJets-Mu-50-madgraph-PUS1_Open_MuTauStream.root";
TString fBackgroundNameDYMuMu = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleDYJets-Mu-50-madgraph-PUS1_Open_MuTauStream.root";
TString fBackgroundNameWJets = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleWJets-Mu-madgraph-PUS1_Open_MuTauStream.root";
TString fBackgroundNameQCD = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleQCDmu_Open_MuTauStream.root";
TString fBackgroundNameTTbar = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleTTJets-Mu-madgraph-PUS1_Open_MuTauStream.root";
TString fBackgroundNameSingleTop = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleSingleTop-Mu_Open_MuTauStream.root";
TString fBackgroundNameDiBoson = "/data_CMS/cms/lbianchini/VbfJetsStudy/OpenNtuples/MuTauStream2011_iter2/Inclusive/nTupleDiBoson-Mu_Open_MuTauStream.root";
TFile *fData(0);
TFile *fSignalVBF(0);
TFile *fSignalGGH(0);
TFile *fBackgroundDYTauTau(0);
TFile *fBackgroundDYMuMu(0);
TFile *fBackgroundWJets(0);
TFile *fBackgroundQCD(0);
TFile *fBackgroundTTbar(0);
TFile *fBackgroundSingleTop(0);
TFile *fBackgroundDiBoson(0);
fData = TFile::Open( fDataName );
fSignalVBF = TFile::Open( fSignalNameVBF );
fSignalGGH = TFile::Open( fSignalNameGGH );
fBackgroundDYTauTau = TFile::Open( fBackgroundNameDYTauTau );
fBackgroundDYMuMu = TFile::Open( fBackgroundNameDYMuMu );
fBackgroundWJets = TFile::Open( fBackgroundNameWJets );
fBackgroundQCD = TFile::Open( fBackgroundNameQCD );
fBackgroundTTbar = TFile::Open( fBackgroundNameTTbar );
fBackgroundSingleTop= TFile::Open( fBackgroundNameSingleTop );
fBackgroundDiBoson = TFile::Open( fBackgroundNameDiBoson );
if(!fSignalVBF || !fBackgroundDYTauTau || !fBackgroundWJets || !fBackgroundQCD || !fBackgroundTTbar || !fBackgroundSingleTop ||
!fSignalGGH || !fBackgroundDYMuMu || !fBackgroundDiBoson || !fData){
std::cout << "ERROR: could not open files" << std::endl;
exit(1);
}
TString tree = "outTreePtOrd";
TTree *data = (TTree*)fData->Get(tree);
TTree *signalVBF = (TTree*)fSignalVBF->Get(tree);
TTree *signalGGH = (TTree*)fSignalGGH->Get(tree);
TFile* dummy1 = new TFile("dummy1.root","RECREATE");
TTree *backgroundDYTauTau = ((TTree*)fBackgroundDYTauTau->Get(tree))->CopyTree("isTauLegMatched>0.5");
TTree *backgroundDYMuMu = ((TTree*)fBackgroundDYMuMu->Get(tree))->CopyTree("isTauLegMatched<0.5");
cout <<backgroundDYTauTau->GetEntries() << " -- " << backgroundDYMuMu->GetEntries() << endl;
TTree *backgroundWJets = (TTree*)fBackgroundWJets->Get(tree);
TTree *backgroundQCD = (TTree*)fBackgroundQCD->Get(tree);
TTree *backgroundTTbar = (TTree*)fBackgroundTTbar->Get(tree);
TTree *backgroundSingleTop = (TTree*)fBackgroundSingleTop->Get(tree);
TTree *backgroundDiBoson = (TTree*)fBackgroundDiBoson->Get(tree);
// here I choose the order in the stack
std::vector<string> samples;
samples.push_back("ggH115");
samples.push_back("qqH115");
samples.push_back("DiBoson");
samples.push_back("SingleTop");
samples.push_back("TTbar");
samples.push_back("Wjets");
samples.push_back("QCD");
samples.push_back("ZfakeTau");
samples.push_back("Ztautau");
samples.push_back("Data");
// here I define the map between a sample name and its tree
std::map<std::string,TTree*> tMap;
tMap["Data"]=data;
tMap["ggH115"]=signalGGH;
tMap["qqH115"]=signalVBF;
tMap["Ztautau"]=backgroundDYTauTau;
tMap["ZfakeTau"]=backgroundDYMuMu;
tMap["Wjets"]=backgroundWJets;
tMap["QCD"]=backgroundQCD;
tMap["TTbar"]=backgroundTTbar;
tMap["SingleTop"]=backgroundSingleTop;
tMap["DiBoson"]=backgroundDiBoson;
Float_t pt1, pt2, Deta, Mjj, eta1, eta2;
Float_t diTauSVFitPt,diTauVisPt,diTauVisEta,diTauSVFitEta,diTauVisMass,diTauSVFitMass,ptL1,ptL2,etaL1,etaL2,diTauCharge,MtLeg1,numPV,combRelIsoLeg1DBeta,combRelIsoLeg1,pZetaCutVar, MEt, ptVeto, leadTrackPt, jetsBtagHE1, jetsBtagHE2;
Float_t sampleWeight,puWeight;
std::map<TString,Float_t> vMap;
for( unsigned iter=0; iter<samples.size(); iter++){
std::map<std::string,TTree*>::iterator it = tMap.find(samples[iter]);
TString h1Name = "h1_"+it->first;
TH1F* h1 = new TH1F( h1Name ,Form("CMS Preliminary 2011 #sqrt{s}=7 TeV L=%.0f pb^{-1}", Lumi) , nBins_ ,xMin_ , xMax_);
TFile* dummy = new TFile("dummy.root","RECREATE");
TCut Cuts = Cuts_;
TCut hlt = enableHLTmatching_ ? "( ((HLTmatch==2 && run<=163261) || (HLTmatch==1 && run>163261)) && ((HLTmu==1 && run<=163261) || (HLTx==1 && run>163261)) )" : "( ((HLTmu==1 && run<=163261) || (HLTx==1 && run>163261)) )";
if((it->first).find("Data")!=string::npos) Cuts = Cuts_ && hlt;
TTree* currentTree = (TTree*)(it->second)->CopyTree(Cuts);
Int_t counter = 0;
currentTree->SetBranchAddress( "pt1", &pt1 );
currentTree->SetBranchAddress( "pt2", &pt2 );
currentTree->SetBranchAddress( "Deta",&Deta );
currentTree->SetBranchAddress( "Mjj", &Mjj );
currentTree->SetBranchAddress( "diTauSVFitPt",&diTauSVFitPt);
currentTree->SetBranchAddress( "diTauSVFitEta",&diTauSVFitEta);
currentTree->SetBranchAddress( "diTauSVFitMass",&diTauSVFitMass);
currentTree->SetBranchAddress( "diTauVisPt",&diTauVisPt);
currentTree->SetBranchAddress( "diTauVisEta",&diTauVisEta);
currentTree->SetBranchAddress( "diTauVisMass",&diTauVisMass);
currentTree->SetBranchAddress( "ptL1", &ptL1 );
currentTree->SetBranchAddress( "ptL2", &ptL2 );
currentTree->SetBranchAddress( "etaL1", &etaL1 );
currentTree->SetBranchAddress( "etaL2", &etaL2 );
currentTree->SetBranchAddress( "combRelIsoLeg1DBeta",&combRelIsoLeg1DBeta);
currentTree->SetBranchAddress( "combRelIsoLeg1",&combRelIsoLeg1);
currentTree->SetBranchAddress( "diTauCharge",&diTauCharge);
currentTree->SetBranchAddress( "MtLeg1",&MtLeg1);
currentTree->SetBranchAddress( "pZetaCutVar",&pZetaCutVar);
currentTree->SetBranchAddress( "numPV",&numPV);
currentTree->SetBranchAddress( "sampleWeight",&sampleWeight);
currentTree->SetBranchAddress( "puWeight",&puWeight);
currentTree->SetBranchAddress( "ptVeto",&ptVeto);
currentTree->SetBranchAddress( "MEt",&MEt);
currentTree->SetBranchAddress( "leadTrackPt",&leadTrackPt);
currentTree->SetBranchAddress( "jetsBtagHE1",&jetsBtagHE1);
currentTree->SetBranchAddress( "jetsBtagHE2",&jetsBtagHE2);
for (Long64_t ievt=0; ievt<currentTree->GetEntries();ievt++) {
currentTree->GetEntry(ievt);
if (ievt%10000 == 0){
std::cout << (it->first) << " ---> processing event: " << ievt << " ..." <<std::endl;
}
vMap["diTauSVFitPt"] = diTauSVFitPt;
vMap["diTauSVFitEta"] = diTauSVFitEta;
vMap["diTauSVFitMass"]= diTauSVFitMass;
vMap["diTauVisPt"] = diTauVisPt;
vMap["diTauVisEta"] = diTauVisEta;
vMap["diTauVisMass"] = diTauVisMass;
vMap["ptL1"] = ptL1;
vMap["ptL2"] = ptL2;
vMap["etaL1"] = etaL1;
vMap["etaL2"] = etaL2;
vMap["diTauCharge"]= Float_t(diTauCharge);
vMap["MtLeg1"]= MtLeg1;
vMap["pZetaCutVar"] = pZetaCutVar;
vMap["numPV"]= numPV;
vMap["combRelIsoLeg1"]= combRelIsoLeg1;
vMap["combRelIsoLeg1DBeta"]= combRelIsoLeg1DBeta;
vMap["sampleWeight"]= sampleWeight;
vMap["puWeight"]= puWeight;
vMap["puWeight"]= puWeight;
vMap["jetsBtagHE1"]= jetsBtagHE1;
vMap["jetsBtagHE2"]= jetsBtagHE2;
vMap["leadTrackPt"]= leadTrackPt;
vMap["pt1"]= pt1;
vMap["pt2"]= pt2;
vMap["Deta"]= Deta;
vMap["Mjj"]= Mjj;
vMap["ptVeto"]= ptVeto;
if((it->first).find("Data")==string::npos) sampleWeight*=(Lumi/1000*hltEff_*puWeight);
counter++;
h1->Fill( vMap[variable_], sampleWeight);
}// end loop
if( (it->first).find("ZfakeTau")!=string::npos ) {
h1->SetFillColor(7);
leg->AddEntry(h1,"Z+jets, fake-#tau","F");
}
if( (it->first).find("Ztautau")!=string::npos ) {
h1->SetFillColor(kRed);
leg->AddEntry(h1,"Z+jets, genuine-#tau","F");
}
if( (it->first).find("TTbar")!=string::npos ) {
h1->SetFillColor(kBlue);
leg->AddEntry(h1,"t#bar{t}+jets","F");
}
if( (it->first).find("SingleTop")!=string::npos ) {
h1->SetFillColor(29);
leg->AddEntry(h1,"single-t","F");
}
if( (it->first).find("Wjets")!=string::npos ) {
h1->SetFillColor(kGreen);
leg->AddEntry(h1,"W+jets","F");
}
if( (it->first).find("DiBoson")!=string::npos ) {
h1->SetFillColor(38);
leg->AddEntry(h1,"Di-Boson","F");
}
if( (it->first).find("QCD")!=string::npos ) {
h1->SetFillColor(42);
leg->AddEntry(h1,"QCD-multijets","F");
}
if((it->first).find("qqH115")!=string::npos){
hSgn1 = (TH1F*)h1->Clone("hSgn1");
hSgn1->Scale(magnifySgn_);
hSgn1->SetLineWidth(2);
h1->SetFillColor(kBlack);
h1->SetFillStyle(3004);
h1->SetLineColor(kBlack);
leg->AddEntry(h1,Form("VBF H(115) X %.0f",magnifySgn_),"F");
}
if((it->first).find("ggH115")!=string::npos){
hSgn2 = (TH1F*)h1->Clone("hSgn2");
hSgn2->Scale(magnifySgn_);
hSgn2->SetLineWidth(2);
h1->SetFillColor(kBlack);
h1->SetFillStyle(3005);
h1->SetLineColor(kBlack);
leg->AddEntry(h1,Form("GGF H(115) X %.0f",magnifySgn_),"F");
}
if((it->first).find("Data")!=string::npos){
hData = (TH1F*)h1->Clone("hData");
hData->Sumw2();
hData->SetMarkerStyle(20);
hData->SetMarkerSize(1.2);
hData->SetMarkerColor(kBlack);
hData->SetLineColor(kBlack);
hData->SetXTitle(XTitle_+" ("+Unities_+")");
hData->SetYTitle(Form(" Events/(%.0f %s)", hData->GetBinWidth(1), Unities_.Data() ) );
hData->SetTitleSize(0.04,"X");
hData->SetTitleSize(0.05,"Y");
hData->SetTitleOffset(0.95,"Y");
leg->AddEntry(hData,"DATA","P");
}
if(iter==0) hSiml=(TH1F*)h1->Clone("hSiml");
else if((it->first).find("Data")==string::npos) hSiml->Add(h1);
if((it->first).find("Data")==string::npos) aStack->Add(h1);
if(VERBOSE) cout<<(it->first) << " ==> "
<< h1->Integral() << " +/- "
<< TMath::Sqrt(h1->GetEntries())*(h1->Integral()/h1->GetEntries())
<< endl;
}
// all signal summed together:
hSgn = (TH1F*)hSgn1->Clone("hSgn");
hSgn->Add(hSgn1,hSgn2,1,1);
if(VERBOSE) cout<< "S/sqrt(B) ==> "
<< hSgn->Integral()/ TMath::Sqrt(hSiml->Integral()) << " +/- "
<< (1./2)*TMath::Sqrt(hSiml->GetEntries())*(hSiml->GetSumOfWeights())/hSiml->Integral()*( hSgn->Integral()/ TMath::Sqrt(hSiml->Integral()) )
<< endl;
hData->SetXTitle(XTitle_+" ("+Unities_+")");
hData->SetYTitle(Form(" Events/(%.1f %s)", hData->GetBinWidth(1), Unities_.Data() ) );
hData->SetTitleSize(0.04,"X");
hData->SetTitleSize(0.05,"Y");
hData->SetTitleOffset(0.95,"Y");
hData->Draw("P");
aStack->Draw("HISTSAME");
hData->Draw("PSAME");
//hSgn1->Draw("HISTSAME");
//hSgn2->Draw("HISTSAME");
TH1F* hStack = (TH1F*)aStack->GetHistogram();
hStack->SetXTitle(XTitle_+" ("+Unities_+")");
hStack->SetYTitle(Form(" Events/(%.0f %s)", hStack->GetBinWidth(1), Unities_.Data() ) );
hStack->SetTitleSize(0.04,"X");
hStack->SetTitleSize(0.05,"Y");
hStack->SetTitleOffset(0.95,"Y");
//aStack->GetYaxis()->SetRangeUser(0.0,std::max(hData->GetMaximum(),hStack->GetMaximum())*1.1);
leg->Draw();
}
