Commit 67b677d8 authored by Justin Griffiths's avatar Justin Griffiths Committed by Graeme Stewart
Browse files

r21 tag (TauHistUtils-00-03-00)

	* make compatible with TauJet_v3
parent 5c6855b8
......@@ -33,8 +33,8 @@ void GeneralTauPlots::fill(const xAOD::TauJet& tau) {
m_oParamPlots.fill(tau);
m_tauCharge->Fill(tau.charge(), 1.);
m_tauNCoreTracks->Fill(tau.nTracks(), 1.);
m_tauNWideTracks->Fill(tau.nWideTracks(), 1.);
if( tau.hasDiscriminant(xAOD::TauJetParameters::BDTJetScore) )m_id_BDTJetScore->Fill(tau.discriminant(xAOD::TauJetParameters::BDTJetScore));
// m_tauNWideTracks->Fill(tau.nWideTracks(), 1.);
m_tauNWideTracks->Fill(tau.nTracks(xAOD::TauJetParameters::classifiedIsolation), 1.); //
if( tau.hasDiscriminant(xAOD::TauJetParameters::BDTJetScore) )m_id_BDTJetScore->Fill(tau.discriminant(xAOD::TauJetParameters::BDTJetScore));
if(tau.isTau(xAOD::TauJetParameters::JetBDTSigLoose)) m_pt_jetBDTloose->Fill( tau.pt()/1000,1);
if(tau.isTau(xAOD::TauJetParameters::JetBDTSigMedium)) m_pt_jetBDTmed ->Fill( tau.pt()/1000,1);
......
......@@ -23,7 +23,7 @@ void NewCorePlots::initializePlots(){
m_massTrkSys = Book1D("massTrkSys",m_sTauJetContainerName + "Tau massTrkSys ;massTrkSys ;# of Taus", 20, 0, 2500 );
m_trkWidth2 = Book1D("trkWidth2",m_sTauJetContainerName + "Tau trkWidth2 ;trkWidth2 ;# of Taus", 20, 0, 1.0 );
m_trFlightPathSig = Book1D("trFlightPathSig",m_sTauJetContainerName + "Tau trFlightPathSig ;trFlightPathSig ;# of Taus", 100, -10, 20 );
m_nPi0 = Book1D("NPio",m_sTauJetContainerName + "Tau NPio ;NPi0 ;# of Taus", 10, 0, 15 );
// m_nPi0 = Book1D("NPio",m_sTauJetContainerName + "Tau NPio ;NPi0 ;# of Taus", 10, 0, 15 );
m_EMRadius = Book1D("Seed_Calo_EMRadius",m_sTauJetContainerName + "Tau Seed Calo EMRadius ;EMRadius ;# of Taus", 10, -0.5,1 );
m_hadRadius = Book1D("Seed_Calo_hadRadius",m_sTauJetContainerName + "Tau Seed Calo hadRadius ;hadRadius ;# of Taus", 10, -0.5,2 );
m_isolFrac = Book1D("Seed_Calo_isolFrac",m_sTauJetContainerName + "Tau Seed Calo isolFrac ;isolFrac ;# of Taus", 40, -0.0,2 );
......@@ -37,9 +37,9 @@ void NewCorePlots::initializePlots(){
m_caloIsoCorrected = Book1D("Seed_Calo_caloIsoCorrected",m_sTauJetContainerName + "Tau Seed Calo caloIsoCorrected ;caloIsoCorrected ;# of Taus", 15, 0, 200 );
m_tauDRMax = Book1D("DRMax",m_sTauJetContainerName + " Tau DR Max track-seed; DRMax; # Taus",10,0.,0.5);
m_tauSflight = Book1D("Sflight",m_sTauJetContainerName + " Tau flight sign. ; Sflight; # Taus",100,-10.,20.);
m_mEflowTopo = Book1D("mEflowTopo",m_sTauJetContainerName + "Tau mEflowTopo; mEflowTopo; #Taus",30,0,300);
m_ptRatioEflowTopo = Book1D("ptRatioEflowTopo",m_sTauJetContainerName + "Tau ptRatioEflowTopo; ptRatioEflowTopo; # of Taus", 10, 0, 10);
m_nPi0Topo = Book1D("nPi0Topo",m_sTauJetContainerName + "Tau nPi0Topo; nPi0Topo; # of Taus", 100, -1, 1);
// m_mEflowTopo = Book1D("mEflowTopo",m_sTauJetContainerName + "Tau mEflowTopo; mEflowTopo; #Taus",30,0,300);
// m_ptRatioEflowTopo = Book1D("ptRatioEflowTopo",m_sTauJetContainerName + "Tau ptRatioEflowTopo; ptRatioEflowTopo; # of Taus", 10, 0, 10);
// m_nPi0Topo = Book1D("nPi0Topo",m_sTauJetContainerName + "Tau nPi0Topo; nPi0Topo; # of Taus", 100, -1, 1);
m_PSSFraction = Book1D("PSSFraction",m_sTauJetContainerName + "Tau PSSFraction; PSSFraction; # of Taus", 30, -15, 15);
m_ChPiEMEOverCaloEME = Book1D("ChPiEMEOverCaloEME",m_sTauJetContainerName + "Tau ChPiEMEOverCaloEME; ChPiEMEOverCaloEME; # of Taus", 30, -15, 15);
m_EMPOverTrkSysP = Book1D("EMPOverTrkSysP",m_sTauJetContainerName + "Tau EMPOverTrkSysP; EMPOverTrkSysP; # of Taus", 81, -1, 80);
......@@ -84,8 +84,8 @@ void NewCorePlots::fill(const xAOD::TauJet& tau) {
// test=tau.detail(xAOD::TauJetParameters::mEflow, avariable);
// if(test) m_mEflow->Fill(avariable,1.);
test=tau.detail(xAOD::TauJetParameters::nPi0, bvariable);
if(test) m_nPi0->Fill(bvariable,1.);
// test=tau.detail(xAOD::TauJetParameters::nPi0, bvariable);
// if(test) m_nPi0->Fill(bvariable,1.);
test=tau.detail(xAOD::TauJetParameters::EMRadius, avariable);
if(test) m_EMRadius->Fill(avariable,1.);
......@@ -120,14 +120,14 @@ void NewCorePlots::fill(const xAOD::TauJet& tau) {
test=tau.detail(xAOD::TauJetParameters::caloIsoCorrected, avariable);
if(test) m_caloIsoCorrected->Fill(avariable,1.);
test=tau.detail(xAOD::TauJetParameters::mEflowTopo, avariable);
if(test) m_mEflowTopo ->Fill(avariable/1000,1.);
// test=tau.detail(xAOD::TauJetParameters::mEflowTopo, avariable);
// if(test) m_mEflowTopo ->Fill(avariable/1000,1.);
test=tau.detail(xAOD::TauJetParameters::ptRatioEflowTopo, avariable);
if(test) m_ptRatioEflowTopo->Fill(avariable,1.);
// test=tau.detail(xAOD::TauJetParameters::ptRatioEflowTopo, avariable);
// if(test) m_ptRatioEflowTopo->Fill(avariable,1.);
test=tau.detail(xAOD::TauJetParameters::nPi0Topo, bvariable);
if(test) m_nPi0Topo->Fill(bvariable,1.);
// test=tau.detail(xAOD::TauJetParameters::nPi0Topo, bvariable);
// if(test) m_nPi0Topo->Fill(bvariable,1.);
test=tau.detail(xAOD::TauJetParameters::PSSFraction, avariable);
if(test) m_PSSFraction->Fill(avariable,1.);
......
......@@ -163,64 +163,64 @@ namespace Tau{
int panTauIntiger = 0;
float panTauFloat = 0;
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_isPanTauCandidate,panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_isPanTauCandidate,panTauIntiger);
m_pantau_CellBasedInput_isPanTauCandidate->Fill(panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_DecayMode,panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_DecayMode,panTauIntiger);
m_pantau_CellBasedInput_DecayMode->Fill(panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTValue_1p0n_vs_1p1n,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTValue_1p0n_vs_1p1n,panTauFloat);
m_pantau_CellBasedInput_BDTValue_1p0n_vs_1p1n->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTValue_1p1n_vs_1pXn,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTValue_1p1n_vs_1pXn,panTauFloat);
m_pantau_CellBasedInput_BDTValue_1p1n_vs_1pXn->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTValue_3p0n_vs_3pXn,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTValue_3p0n_vs_3pXn,panTauFloat);
m_pantau_CellBasedInput_BDTValue_3p0n_vs_3pXn->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Basic_NNeutralConsts,panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Basic_NNeutralConsts,panTauIntiger);
m_pantau_CellBasedInput_BDTVar_Basic_NNeutralConsts->Fill(panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Charged_JetMoment_EtDRxTotalEt,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Charged_JetMoment_EtDRxTotalEt,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Charged_JetMoment_EtDRxTotalEt->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Charged_StdDev_Et_WrtEtAllConsts,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Charged_StdDev_Et_WrtEtAllConsts->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Neutral_HLV_SumM,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Neutral_HLV_SumM,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Neutral_HLV_SumM->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Neutral_PID_BDTValues_BDTSort_1,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Neutral_PID_BDTValues_BDTSort_1->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Neutral_Shots_NPhotonsInSeed,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Neutral_Shots_NPhotonsInSeed,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Neutral_Shots_NPhotonsInSeed->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_CellBasedInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged,panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged,panTauFloat);
m_pantau_CellBasedInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_isPanTauCandidate,panTauIntiger);
m_pantau_eflowRecInput_isPanTauCandidate->Fill(panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_DecayMode,panTauIntiger);
m_pantau_eflowRecInput_DecayMode->Fill(panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTValue_1p0n_vs_1p1n,panTauFloat);
m_pantau_eflowRecInput_BDTValue_1p0n_vs_1p1n->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTValue_1p1n_vs_1pXn,panTauFloat);
m_pantau_eflowRecInput_BDTValue_1p1n_vs_1pXn->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTValue_3p0n_vs_3pXn,panTauFloat);
m_pantau_eflowRecInput_BDTValue_3p0n_vs_3pXn->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Basic_NPi0NeutConsts,panTauIntiger);
m_pantau_eflowRecInput_BDTVar_Basic_NPi0NeutConsts->Fill(panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Basic_NNeutralConsts,panTauIntiger);
m_pantau_eflowRecInput_BDTVar_Basic_NNeutralConsts->Fill(panTauIntiger);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Charged_Ratio_EtOverEtAllConsts,panTauFloat);
m_pantau_eflowRecInput_BDTVar_Charged_Ratio_EtOverEtAllConsts->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_HLV_SumM,panTauFloat);
m_pantau_eflowRecInput_BDTVar_Neutral_HLV_SumM->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_EtSort_1,panTauFloat);
m_pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_EtSort_1->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2,panTauFloat);
m_pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts,panTauFloat);
m_pantau_eflowRecInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_Mean_DRToLeading_WrtEtAllConsts,panTauFloat);
m_pantau_eflowRecInput_BDTVar_Neutral_Mean_DRToLeading_WrtEtAllConsts->Fill(panTauFloat);
thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged,panTauFloat);
m_pantau_eflowRecInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_isPanTauCandidate,panTauIntiger);
// m_pantau_eflowRecInput_isPanTauCandidate->Fill(panTauIntiger);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_DecayMode,panTauIntiger);
// m_pantau_eflowRecInput_DecayMode->Fill(panTauIntiger);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTValue_1p0n_vs_1p1n,panTauFloat);
// m_pantau_eflowRecInput_BDTValue_1p0n_vs_1p1n->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTValue_1p1n_vs_1pXn,panTauFloat);
// m_pantau_eflowRecInput_BDTValue_1p1n_vs_1pXn->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTValue_3p0n_vs_3pXn,panTauFloat);
// m_pantau_eflowRecInput_BDTValue_3p0n_vs_3pXn->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Basic_NPi0NeutConsts,panTauIntiger);
// m_pantau_eflowRecInput_BDTVar_Basic_NPi0NeutConsts->Fill(panTauIntiger);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Basic_NNeutralConsts,panTauIntiger);
// m_pantau_eflowRecInput_BDTVar_Basic_NNeutralConsts->Fill(panTauIntiger);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Charged_Ratio_EtOverEtAllConsts,panTauFloat);
// m_pantau_eflowRecInput_BDTVar_Charged_Ratio_EtOverEtAllConsts->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_HLV_SumM,panTauFloat);
// m_pantau_eflowRecInput_BDTVar_Neutral_HLV_SumM->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_EtSort_1,panTauFloat);
// m_pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_EtSort_1->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2,panTauFloat);
// m_pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts,panTauFloat);
// m_pantau_eflowRecInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Neutral_Mean_DRToLeading_WrtEtAllConsts,panTauFloat);
// m_pantau_eflowRecInput_BDTVar_Neutral_Mean_DRToLeading_WrtEtAllConsts->Fill(panTauFloat);
// thisTau.panTauDetail( xAOD::TauJetParameters::PanTauDetails::pantau_eflowRecInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged,panTauFloat);
// m_pantau_eflowRecInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged->Fill(panTauFloat);
m_oParamPlots.fill(thisTau);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment