Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 14:21:04

0001 void maketable_HLT(string model, string oldversion, string newversion, bool oldrelease) {
0002   
0003   //  TFile* file = new TFile(string("html/SUSYVal/histo_"+model+"/outputfile.root").c_str());
0004   TFile* file = new TFile(string("html/SUSYVal_"+oldversion+"_vs_"+newversion+"/"+model+"/outputfile.root").c_str());
0005   TH1D* histo1;
0006   TH1D* histo2;
0007   vector<string> tablelines;
0008   
0009   bool oldrelease = true;
0010 
0011   vector<double> e_eff;
0012   vector<string> e_line;
0013 
0014   vector<double> p_eff;
0015   vector<string> p_line;
0016 
0017   vector<double> jet_eff;
0018   vector<string> jet_line;
0019 
0020   vector<double > m_eff;
0021   vector<string> m_line;
0022 
0023   for(int iplot=0; iplot<8;iplot++) {
0024     char name[256];
0025     if(oldrelease == true)
0026       sprintf(name,"HltPaths_eff_%i",iplot);
0027     else
0028       sprintf(name,"HltPaths_eff_%i",iplot+8);
0029     histo = (TH1D*) file->Get(name);
0030 
0031     for(int iline=0; iline< histo->GetXaxis()->GetNbins(); iline++) {
0032       if(string(histo->GetXaxis()->GetBinLabel(iline+1)).find("Total") 
0033      != string::npos) 
0034     break;
0035       
0036       double eff = histo->GetBinContent(iline+1);
0037       double err = histo->GetBinError(iline+1);
0038       
0039       sprintf(name,"| %s | %f +/- %f |",histo->GetXaxis()->GetBinLabel(iline+1),eff, err);
0040 
0041       if(string(name).find("Electron") != string::npos) {
0042     e_eff.push_back(eff);
0043     e_line.push_back(string(name));
0044       }
0045 
0046       if(string(name).find("Photon") != string::npos) {
0047     p_eff.push_back(eff);
0048     p_line.push_back(string(name));
0049       }
0050 
0051       if(string(name).find("Jet") != string::npos ||
0052      string(name).find("jet") != string::npos) {
0053     jet_eff.push_back(eff);
0054     jet_line.push_back(string(name));
0055       }
0056 
0057       if(string(name).find("Mu") != string::npos) {
0058     m_eff.push_back(eff);
0059     m_line.push_back(string(name));
0060       }
0061     }
0062   }
0063   
0064   sort(e_eff.begin(),   e_eff.end());
0065   sort(p_eff.begin(),   p_eff.end());
0066   sort(jet_eff.begin(), jet_eff.end());
0067   sort(m_eff.begin(),   m_eff.end());
0068 
0069 
0070   vector<string> firstfour;
0071   for(int i=0; i<4; i++) {
0072     char eff_ch[256];
0073     sprintf(eff_ch,"%d",e_eff[i]);
0074     for(int j=0; j<int(e_line.size()); j++) {
0075       if(e_line[j].find(eff_ch) != string::npos) {
0076     firstfour.push_back(e_line[j]);
0077     e_line[j] = string("none");
0078     break;
0079       }
0080     }
0081   }
0082   
0083   for(int i=0; i<firstfour.size(); i++) 
0084     cout << firstfour[i] << endl;
0085 
0086   firstfour.clear();
0087 
0088   for(int i=0; i<4; i++) {
0089     char eff_ch[256];
0090     sprintf(eff_ch,"%d",p_eff[i]);
0091     for(int j=0; j<int(p_line.size()); j++) {
0092       if(p_line[j].find(eff_ch) != string::npos) {
0093     firstfour.push_back(p_line[j]);
0094     p_line[j] = "none";
0095     break;
0096       }
0097     }
0098   }
0099 
0100   for(int i=0; i<firstfour.size(); i++) 
0101     cout << firstfour[i] << endl;
0102 
0103   firstfour.clear();
0104 
0105   for(int i=0; i<4; i++) {
0106     char eff_ch[256];
0107     sprintf(eff_ch,"%d",jet_eff[i]);
0108     for(int j=0; j<int(jet_line.size()); j++) {
0109       if(jet_line[j].find(eff_ch) != string::npos) {
0110     firstfour.push_back(jet_line[j]);
0111     jet_line[j] = "none";
0112     break;
0113       }
0114     }
0115   }
0116   
0117   for(int i=0; i<firstfour.size(); i++) 
0118     cout << firstfour[i] << endl;
0119 
0120   firstfour.clear();
0121 
0122   for(int i=0; i<4; i++) {
0123     char eff_ch[256];
0124     sprintf(eff_ch,"%d",m_eff[i]);
0125     for(int j=0; j<int(m_line.size()); j++) {
0126       if(m_line[j].find(eff_ch) != string::npos) {
0127     firstfour.push_back(m_line[j]);
0128     m_line[j] = "none";
0129     break;
0130       }
0131     }
0132   }
0133   
0134   for(int i=0; i<firstfour.size(); i++) 
0135     cout << firstfour[i] << endl;
0136 }
0137 
0138 
0139 void maketable_HLT() {
0140   string oldversion = "183";
0141   string newversion = "200";
0142   cout << "---++ RELEASE: " << oldversion << endl;
0143   cout << "---++ LM1 " << endl;
0144   maketable_HLT("LM1",oldversion,newversion,true);
0145   cout << "---++ LM5 " << endl;
0146   maketable_HLT("LM5",oldversion,newversion,true);
0147   cout << "--- ++ LM9p " << endl;
0148   maketable_HLT("LM9p",oldversion,newversion,true);
0149   cout << "---++ GM1b " << endl;
0150   maketable_HLT("GM1b",oldversion,newversion,true);
0151   cout << "---++ RSgrav " << endl;
0152   maketable_HLT("RSgrav",oldversion,newversion,true);
0153   cout << "---++ Zprime " << endl;
0154   maketable_HLT("Zprime",oldversion,newversion,true);
0155   cout << endl;
0156   cout << endl;
0157   cout << endl;
0158   cout << endl;
0159   cout << endl;
0160   cout << "---++ " << newversion << endl;
0161   cout << "---++ LM1 " << endl;
0162   maketable_HLT("LM1",oldversion,newversion,false);
0163   cout << "---++ LM5 " << endl;
0164   maketable_HLT("LM5",oldversion,newversion,false);
0165   cout << "---++ LM9p " << endl;
0166   maketable_HLT("LM9p",oldversion,newversion,false);
0167   cout << "---++ GM1b " << endl;
0168   maketable_HLT("GM1b",oldversion,newversion,false);
0169   cout << "---++ RSgrav " << endl;
0170   maketable_HLT("RSgrav",oldversion,newversion,false);
0171   cout << "---++ Zprime " << endl;
0172   maketable_HLT("Zprime",oldversion,newversion,false);
0173 }
0174 
0175