Commit a5e8835d authored by Vakhtang Tsulaia's avatar Vakhtang Tsulaia
Browse files

Merge branch 'master-response-histos' into 'master'

Adding response histograms without shift to zero

See merge request atlas/athena!47472
parents c8fac2e2 29e82278
......@@ -37,6 +37,10 @@ protected:
TProfile* m_etres_eta;
TProfile* m_etres_pt;
TH1F* m_etres_noShift;
TProfile* m_etres_noShift_eta;
TProfile* m_etres_noShift_pt;
TH1F* m_deltaRclosest;
};
......
......@@ -46,6 +46,10 @@ int EfficiencyResponseHistos::buildHistos(){
m_etres_eta =bookHisto( hbuilder.build<TProfile>("erhResponseVsEta") );
m_etres_pt =bookHisto( hbuilder.build<TProfile>("erhResponseVsPt") );
m_etres_noShift = bookHisto( hbuilder.build<TH1F>("erhResponse_noShift") );
m_etres_noShift_eta =bookHisto( hbuilder.build<TProfile>("erhResponseVsEta_noShift") );
m_etres_noShift_pt =bookHisto( hbuilder.build<TProfile>("erhResponseVsPt_noShift") );
m_deltaRclosest = bookHisto( hbuilder.build<TH1F>("erhDeltaR") );
......@@ -85,11 +89,18 @@ int EfficiencyResponseHistos::fillHistosFromContainer(const xAOD::JetContainer &
if( dr < 0.3) {
double relDiff = -999;
if (refPt > 0.) relDiff = ( matched->pt()* toGeV - refPt )/refPt;
double response = -999;
if (refPt > 0.){
relDiff = ( matched->pt()* toGeV - refPt )/refPt;
response = (matched->pt()* toGeV)/refPt;
}
m_etres->Fill( relDiff, weight );
m_etres_eta->Fill( refjet->eta(), relDiff, weight);
m_etres_pt->Fill( refPt, relDiff, weight);
m_etres_noShift->Fill( response, weight );
m_etres_noShift_eta->Fill( refjet->eta(), response, weight);
m_etres_noShift_pt->Fill( refPt, response, weight);
}
}
......
......@@ -218,6 +218,10 @@ jhm.addTool( EfficiencyResponseHistos("effresponse",
hdef('erhResponseVsEta', "Jet p_{T} Response vs #eta;#eta of jet;#frac{p_{T}^{Jet} - p_{T}^{Truth}}{p_{T}^{Truth}}",50,-5,5 ),
hdef('erhResponseVsPt', "Jet p_{T} Response vs p_{T};p_{T}^{Truth} of jet;#frac{p_{T}^{Jet} - p_{T}^{Truth}}{p_{T}^{Truth}}",50,0,1000 ),
hdef('erhResponse_noShift', "Jet p_{T} Response;#frac{p_{T}^{Jet}}{p_{T}^{Truth}};Number of jets",50,0,2 ),
hdef('erhResponseVsEta_noShift', "Jet p_{T} Response vs #eta;#eta of jet;#frac{p_{T}^{Jet}}{p_{T}^{Truth}}",50,-5,5 ),
hdef('erhResponseVsPt_noShift', "Jet p_{T} Response vs p_{T};p_{T}^{Truth} of jet;#frac{p_{T}^{Jet}}{p_{T}^{Truth}}",50,0,1000 ),
hdef('erhDeltaR', "#DeltaR between Jet and closest Truth Jet;#DeltaR;Number of jets",50,0,4 ),
]
......
Markdown is supported
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