| 1 |
#include <iostream>
|
| 2 |
#include <fstream>
|
| 3 |
|
| 4 |
#include "DataFormats/TauReco/interface/PFTau.h"
|
| 5 |
|
| 6 |
|
| 7 |
#include "UserCode/HbbAnalysis/interface/HistosTaus.hh"
|
| 8 |
|
| 9 |
namespace HbbAnalysis {//namespace
|
| 10 |
|
| 11 |
void HistosTaus::Initialise(TFileDirectory & aDir, std::string aName){
|
| 12 |
|
| 13 |
CreateHistos(aName,aDir);
|
| 14 |
|
| 15 |
p_leadTrk_pT = aDir.make<TH1F>("p_leadTrk_pT",";p_{T} of leading track (GeV);n_{entries}/GeV",100,0,100);
|
| 16 |
p_leadTrk_eta = aDir.make<TH1F>("p_leadTrk_eta",";#eta of leading track;n_{entries}",48,-2.4,2.4);
|
| 17 |
p_leadTrk_phi = aDir.make<TH1F>("p_leadTrk_phi",";#phi of leading track;n_{entries}",64,-3.2,3.2);
|
| 18 |
|
| 19 |
p_leadTrk_matchDist = aDir.make<TH1F>("p_leadTrk_matchDist",";#DeltaR(#tau-leading track);n_{entries}",100,-0.5,0.5);
|
| 20 |
|
| 21 |
p_leadTrk_IPxy = aDir.make<TH1F>("p_leadTrk_IPxy",";dxy(leading track-IP);n_{entries}",100,-0.1,0.1);
|
| 22 |
p_leadTrk_IPz = aDir.make<TH1F>("p_leadTrk_IPz",";dz(leading track-IP);n_{entries}",100,-1,1);
|
| 23 |
|
| 24 |
}
|
| 25 |
|
| 26 |
void HistosTaus::FillHistograms(const pat::Tau & aTau, edm::Handle<std::vector<reco::Vertex> > aRecoVertices){//FillHistograms
|
| 27 |
|
| 28 |
if ( aTau.leadTrack().isAvailable() && aTau.leadTrack().isNonnull() ) {
|
| 29 |
p_leadTrk_pT->Fill(aTau.leadTrack()->pt());
|
| 30 |
p_leadTrk_eta->Fill(aTau.leadTrack()->eta());
|
| 31 |
p_leadTrk_phi->Fill(aTau.leadTrack()->phi());
|
| 32 |
|
| 33 |
p_leadTrk_matchDist->Fill(reco::deltaR(aTau.leadTrack()->momentum(), aTau.p4()));
|
| 34 |
|
| 35 |
if ( aRecoVertices->size() >= 1 ) {
|
| 36 |
const reco::Vertex& thePrimaryEventVertex = (*aRecoVertices->begin());
|
| 37 |
p_leadTrk_IPxy->Fill(aTau.leadTrack()->dxy(thePrimaryEventVertex.position()));
|
| 38 |
p_leadTrk_IPz->Fill(aTau.leadTrack()->dz(thePrimaryEventVertex.position()));
|
| 39 |
}
|
| 40 |
}
|
| 41 |
|
| 42 |
FillBaseHistograms(aTau.pt(),aTau.eta(),aTau.phi());
|
| 43 |
|
| 44 |
}//FillHistograms
|
| 45 |
|
| 46 |
}//namespace
|
| 47 |
|
| 48 |
|
