| 1 |
|
#include <iostream> |
| 2 |
|
#include <fstream> |
| 3 |
|
#include <sstream> |
| 4 |
+ |
#include <cmath> |
| 5 |
+ |
#include <cstdlib> |
| 6 |
|
|
| 7 |
|
#include "Math/VectorUtil.h" |
| 8 |
|
|
| 54 |
|
lTitle << "p_bDiscriminator_" << i << "_" << fs; |
| 55 |
|
std::ostringstream lAxis; |
| 56 |
|
lAxis << ";" << btagAlgos.at(i) << ";N_{entries}" ; |
| 57 |
< |
p_bDiscriminator[i][fs] = aDir.make<TH1F>(lTitle.str().c_str(),lAxis.str().c_str(),4480,-101,11);//100,min[i],max[i]); |
| 57 |
> |
p_bDiscriminator[i][fs] = aDir.make<TH1F>(lTitle.str().c_str(),lAxis.str().c_str(),5280,-101,31);//100,min[i],max[i]); |
| 58 |
|
} |
| 59 |
|
} |
| 60 |
|
|
| 232 |
|
if (isMatched) p_bDiscriminator[3][lFlav]->Fill(aJet.btagVars().bProba); |
| 233 |
|
p_bDiscriminator[4][0]->Fill(aJet.btagVars().probability); |
| 234 |
|
if (isMatched) p_bDiscriminator[4][lFlav]->Fill(aJet.btagVars().probability); |
| 235 |
< |
p_bDiscriminator[5][0]->Fill(aJet.btagVars().sSV); |
| 236 |
< |
if (isMatched) p_bDiscriminator[5][lFlav]->Fill(aJet.btagVars().sSV); |
| 237 |
< |
p_bDiscriminator[6][0]->Fill(aJet.btagVars().softElectronByPt); |
| 238 |
< |
if (isMatched) p_bDiscriminator[6][lFlav]->Fill(aJet.btagVars().softElectronByPt); |
| 239 |
< |
p_bDiscriminator[7][0]->Fill(aJet.btagVars().softElectronByIP3d); |
| 240 |
< |
if (isMatched) p_bDiscriminator[7][lFlav]->Fill(aJet.btagVars().softElectronByIP3d); |
| 241 |
< |
p_bDiscriminator[8][0]->Fill(aJet.btagVars().softMuon); |
| 242 |
< |
if (isMatched) p_bDiscriminator[8][lFlav]->Fill(aJet.btagVars().softMuon); |
| 243 |
< |
p_bDiscriminator[9][0]->Fill(aJet.btagVars().softMuonByPt); |
| 244 |
< |
if (isMatched) p_bDiscriminator[9][lFlav]->Fill(aJet.btagVars().softMuonByPt); |
| 245 |
< |
p_bDiscriminator[10][0]->Fill(aJet.btagVars().softMuonByIP3d); |
| 246 |
< |
if (isMatched) p_bDiscriminator[10][lFlav]->Fill(aJet.btagVars().softMuonByIP3d); |
| 247 |
< |
p_bDiscriminator[11][0]->Fill(aJet.btagVars().tCHE); |
| 248 |
< |
if (isMatched) p_bDiscriminator[11][lFlav]->Fill(aJet.btagVars().tCHE); |
| 249 |
< |
p_bDiscriminator[12][0]->Fill(aJet.btagVars().tCHP); |
| 250 |
< |
if (isMatched) p_bDiscriminator[12][lFlav]->Fill(aJet.btagVars().tCHP); |
| 235 |
> |
p_bDiscriminator[5][0]->Fill(aJet.btagVars().sSVHE); |
| 236 |
> |
if (isMatched) p_bDiscriminator[5][lFlav]->Fill(aJet.btagVars().sSVHE); |
| 237 |
> |
p_bDiscriminator[6][0]->Fill(aJet.btagVars().sSVHP); |
| 238 |
> |
if (isMatched) p_bDiscriminator[6][lFlav]->Fill(aJet.btagVars().sSVHP); |
| 239 |
> |
p_bDiscriminator[7][0]->Fill(aJet.btagVars().softElectronByPt); |
| 240 |
> |
if (isMatched) p_bDiscriminator[7][lFlav]->Fill(aJet.btagVars().softElectronByPt); |
| 241 |
> |
p_bDiscriminator[8][0]->Fill(aJet.btagVars().softElectronByIP3d); |
| 242 |
> |
if (isMatched) p_bDiscriminator[8][lFlav]->Fill(aJet.btagVars().softElectronByIP3d); |
| 243 |
> |
p_bDiscriminator[9][0]->Fill(aJet.btagVars().softMuon); |
| 244 |
> |
if (isMatched) p_bDiscriminator[9][lFlav]->Fill(aJet.btagVars().softMuon); |
| 245 |
> |
p_bDiscriminator[10][0]->Fill(aJet.btagVars().softMuonByPt); |
| 246 |
> |
if (isMatched) p_bDiscriminator[10][lFlav]->Fill(aJet.btagVars().softMuonByPt); |
| 247 |
> |
p_bDiscriminator[11][0]->Fill(aJet.btagVars().softMuonByIP3d); |
| 248 |
> |
if (isMatched) p_bDiscriminator[11][lFlav]->Fill(aJet.btagVars().softMuonByIP3d); |
| 249 |
> |
p_bDiscriminator[12][0]->Fill(aJet.btagVars().tCHE); |
| 250 |
> |
if (isMatched) p_bDiscriminator[12][lFlav]->Fill(aJet.btagVars().tCHE); |
| 251 |
> |
p_bDiscriminator[13][0]->Fill(aJet.btagVars().tCHP); |
| 252 |
> |
if (isMatched) p_bDiscriminator[13][lFlav]->Fill(aJet.btagVars().tCHP); |
| 253 |
|
|
| 254 |
|
if (aJet.genJetVars().valid){ |
| 255 |
|
p_pTgenjet[0]->Fill(aJet.genJetVars().pT); |
| 306 |
|
p_n90->Fill(aJet.caloVars().n90); |
| 307 |
|
p_n60->Fill(aJet.caloVars().n60); |
| 308 |
|
} |
| 309 |
< |
|
| 309 |
> |
//JPT + PF |
| 310 |
> |
if (aJet.jetType() == 1 || aJet.jetType() == 2){ |
| 311 |
> |
p_chargedHadronEnergy->Fill(aJet.jptPfVars().chargedHadronEnergy); |
| 312 |
> |
p_chargedHadronEnergyFraction->Fill(aJet.jptPfVars().chargedHadronEnergyFraction); |
| 313 |
> |
p_neutralHadronEnergy->Fill(aJet.jptPfVars().neutralHadronEnergy); |
| 314 |
> |
p_neutralHadronEnergyFraction->Fill(aJet.jptPfVars().neutralHadronEnergyFraction); |
| 315 |
> |
p_chargedEmEnergy->Fill(aJet.jptPfVars().chargedEmEnergy); |
| 316 |
> |
p_chargedEmEnergyFraction->Fill(aJet.jptPfVars().chargedEmEnergyFraction); |
| 317 |
> |
p_neutralEmEnergy->Fill(aJet.jptPfVars().neutralEmEnergy); |
| 318 |
> |
p_neutralEmEnergyFraction->Fill(aJet.jptPfVars().neutralEmEnergyFraction); |
| 319 |
> |
p_chargedMultiplicity->Fill(aJet.jptPfVars().chargedMultiplicity); |
| 320 |
> |
p_muonMultiplicity->Fill(aJet.jptPfVars().muonMultiplicity); |
| 321 |
> |
} |
| 322 |
> |
|
| 323 |
|
// //================== PF Jet specific information ==================== |
| 324 |
|
if (aJet.jetType() == 1){ |
| 308 |
– |
p_chargedHadronEnergy->Fill(aJet.pfVars().chargedHadronEnergy); |
| 309 |
– |
p_chargedHadronEnergyFraction->Fill(aJet.pfVars().chargedHadronEnergyFraction); |
| 310 |
– |
p_neutralHadronEnergy->Fill(aJet.pfVars().neutralHadronEnergy); |
| 311 |
– |
p_neutralHadronEnergyFraction->Fill(aJet.pfVars().neutralHadronEnergyFraction); |
| 312 |
– |
p_chargedEmEnergy->Fill(aJet.pfVars().chargedEmEnergy); |
| 313 |
– |
p_chargedEmEnergyFraction->Fill(aJet.pfVars().chargedEmEnergyFraction); |
| 325 |
|
p_chargedMuEnergy->Fill(aJet.pfVars().chargedMuEnergy); |
| 326 |
|
p_chargedMuEnergyFraction->Fill(aJet.pfVars().chargedMuEnergyFraction); |
| 316 |
– |
p_neutralEmEnergy->Fill(aJet.pfVars().neutralEmEnergy); |
| 317 |
– |
p_neutralEmEnergyFraction->Fill(aJet.pfVars().neutralEmEnergyFraction); |
| 318 |
– |
p_chargedMultiplicity->Fill(aJet.pfVars().chargedMultiplicity); |
| 327 |
|
p_neutralMultiplicity->Fill(aJet.pfVars().neutralMultiplicity); |
| 320 |
– |
p_muonMultiplicity->Fill(aJet.pfVars().muonMultiplicity); |
| 328 |
|
} |
| 329 |
|
|
| 330 |
|
|
