|
|
|
#include <MuonTrackValidator.h>
Public Member Functions | |
| void | analyze (const edm::Event &, const edm::EventSetup &) |
| Method called once per event. | |
| void | beginRun (edm::Run const &, edm::EventSetup const &) |
| Method called before the event loop. | |
| void | endRun (edm::Run const &, edm::EventSetup const &) |
| Method called at the end of the event loop. | |
| MuonTrackValidator (const edm::ParameterSet &pset) | |
| Constructor. | |
| virtual | ~MuonTrackValidator () |
| Destructor. | |
Private Member Functions | |
| void | getRecoMomentum (const reco::Track &track, double &pt, double &ptError, double &qoverp, double &qoverpError, double &lambda, double &lambdaError, double &phi, double &phiError) const |
| retrieval of reconstructed momentum components from reco::Track (== mean values for GSF) | |
| void | getRecoMomentum (const reco::GsfTrack &gsfTrack, double &pt, double &ptError, double &qoverp, double &qoverpError, double &lambda, double &lambdaError, double &phi, double &phiError) const |
| retrieval of reconstructed momentum components based on the mode of a reco::GsfTrack | |
Private Attributes | |
| edm::InputTag | associatormap |
| bool | BiDirectional_RecoToSim_association |
| std::vector< MonitorElement * > | chi2_vs_eta |
| std::vector< MonitorElement * > | chi2_vs_nhits |
| std::vector< MonitorElement * > | chi2_vs_phi |
| CosmicTrackingParticleSelector | cosmictpSelector |
| std::vector< MonitorElement * > | cotThetares_vs_eta |
| std::vector< MonitorElement * > | cotThetares_vs_pt |
| std::string | dirName_ |
| std::vector< MonitorElement * > | dxypull_vs_eta |
| std::vector< MonitorElement * > | dxyres_vs_eta |
| std::vector< MonitorElement * > | dxyres_vs_pt |
| std::vector< MonitorElement * > | dzpull_vs_eta |
| std::vector< MonitorElement * > | dzres_vs_eta |
| std::vector< MonitorElement * > | dzres_vs_pt |
| std::vector< MonitorElement * > | etares_vs_eta |
| std::vector< MonitorElement * > | h_assochi2 |
| std::vector< MonitorElement * > | h_assochi2_prob |
| std::vector< MonitorElement * > | h_chi2mean_vs_phi |
| std::vector< MonitorElement * > | h_chi2meanh |
| std::vector< MonitorElement * > | h_chi2meanhitsh |
| std::vector< MonitorElement * > | h_dxypulleta |
| std::vector< MonitorElement * > | h_dzpulleta |
| std::vector< MonitorElement * > | h_hits_phi |
| std::vector< MonitorElement * > | h_losthits |
| std::vector< MonitorElement * > | h_losthits_eta |
| std::vector< MonitorElement * > | h_nchi2 |
| std::vector< MonitorElement * > | h_nchi2_prob |
| std::vector< MonitorElement * > | h_phipulleta |
| std::vector< MonitorElement * > | h_phipullphi |
| std::vector< MonitorElement * > | h_ptpulleta |
| std::vector< MonitorElement * > | h_ptpullphi |
| std::vector< MonitorElement * > | h_ptshifteta |
| std::vector< MonitorElement * > | h_thetapulleta |
| std::vector< MonitorElement * > | h_thetapullphi |
| bool | MABH |
| double | maxPhi |
| double | minPhi |
| std::vector< MonitorElement * > | nhits_vs_phi |
| int | nintPhi |
| std::vector< MonitorElement * > | nlosthits_vs_eta |
| std::vector< MonitorElement * > | phimean_vs_eta_phi |
| std::vector< MonitorElement * > | phipull_vs_eta |
| std::vector< MonitorElement * > | phipull_vs_phi |
| std::vector< MonitorElement * > | phires_vs_eta |
| std::vector< MonitorElement * > | phires_vs_phi |
| std::vector< MonitorElement * > | phires_vs_pt |
| std::vector< MonitorElement * > | ptmean_vs_eta_phi |
| std::vector< MonitorElement * > | ptpull_vs_eta |
| std::vector< MonitorElement * > | ptpull_vs_phi |
| std::vector< MonitorElement * > | ptres_vs_eta |
| std::vector< MonitorElement * > | ptres_vs_phi |
| std::vector< MonitorElement * > | ptres_vs_pt |
| std::vector< MonitorElement * > | thetapull_vs_eta |
| std::vector< MonitorElement * > | thetapull_vs_phi |
| TrackingParticleSelector | tpSelector |
| bool | UseAssociators |
| bool | useGsf |
Class that produces histograms to validate Muon Track Reconstruction performances
Definition at line 17 of file MuonTrackValidator.h.
| MuonTrackValidator::MuonTrackValidator | ( | const edm::ParameterSet & | pset | ) | [inline] |
Constructor.
Definition at line 20 of file MuonTrackValidator.h.
References associatormap, MuonTrackValidatorBase::associators, BiDirectional_RecoToSim_association, cosmictpSelector, dirName_, edm::ParameterSet::dump(), edm::ParameterSet::getParameter(), MuonTrackValidatorBase::label, MABH, maxPhi, minPhi, nintPhi, tpSelector, UseAssociators, useGsf, MuonTrackValidatorBase::usemuon, and MuonTrackValidatorBase::usetracker.
:MuonTrackValidatorBase(pset){ dirName_ = pset.getParameter<std::string>("dirName"); associatormap = pset.getParameter< edm::InputTag >("associatormap"); UseAssociators = pset.getParameter< bool >("UseAssociators"); tpSelector = TrackingParticleSelector(pset.getParameter<double>("ptMinTP"), pset.getParameter<double>("minRapidityTP"), pset.getParameter<double>("maxRapidityTP"), pset.getParameter<double>("tipTP"), pset.getParameter<double>("lipTP"), pset.getParameter<int>("minHitTP"), pset.getParameter<bool>("signalOnlyTP"), pset.getParameter<bool>("chargedOnlyTP"), pset.getParameter<bool>("stableOnlyTP"), pset.getParameter<std::vector<int> >("pdgIdTP")); cosmictpSelector = CosmicTrackingParticleSelector(pset.getParameter<double>("ptMinTP"), pset.getParameter<double>("minRapidityTP"), pset.getParameter<double>("maxRapidityTP"), pset.getParameter<double>("tipTP"), pset.getParameter<double>("lipTP"), pset.getParameter<int>("minHitTP"), pset.getParameter<bool>("chargedOnlyTP"), pset.getParameter<std::vector<int> >("pdgIdTP")); minPhi = pset.getParameter<double>("minPhi"); maxPhi = pset.getParameter<double>("maxPhi"); nintPhi = pset.getParameter<int>("nintPhi"); useGsf = pset.getParameter<bool>("useGsf"); BiDirectional_RecoToSim_association = pset.getParameter<bool>("BiDirectional_RecoToSim_association"); // dump cfg parameters edm::LogVerbatim("MuonTrackValidator") << "constructing MuonTrackValidator: " << pset.dump(); MABH = false; if (!UseAssociators) { // flag MuonAssociatorByHits if (associators[0] == "MuonAssociationByHits") MABH = true; // reset string associators to the map label associators.clear(); associators.push_back(associatormap.label()); edm::LogVerbatim("MuonTrackValidator") << "--> associators reset to: " <<associators[0]; } // inform on which SimHits will be counted if (usetracker) edm::LogVerbatim("MuonTrackValidator") <<"\n usetracker = TRUE : Tracker SimHits WILL be counted"; else edm::LogVerbatim("MuonTrackValidator") <<"\n usetracker = FALSE : Tracker SimHits WILL NOT be counted"; if (usemuon) edm::LogVerbatim("MuonTrackValidator") <<" usemuon = TRUE : Muon SimHits WILL be counted"; else edm::LogVerbatim("MuonTrackValidator") <<" usemuon = FALSE : Muon SimHits WILL NOT be counted"<<std::endl; // loop over the reco::Track collections to validate: check for inconsistent input settings for (unsigned int www=0;www<label.size();www++) { std::string recoTracksLabel = label[www].label(); std::string recoTracksInstance = label[www].instance(); // tracks with hits only on tracker if (recoTracksLabel=="generalTracks" || (recoTracksLabel.find("cutsRecoTracks") != std::string::npos) || recoTracksLabel=="ctfWithMaterialTracksP5LHCNavigation" || recoTracksLabel=="hltL3TkTracksFromL2" || (recoTracksLabel=="hltL3Muons" && recoTracksInstance=="L2Seeded")) { if (usemuon) { edm::LogWarning("MuonTrackValidator") <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www] <<"\n with usemuon == true"<<"\n ---> please change to usemuon == false "; } if (!usetracker) { edm::LogWarning("MuonTrackValidator") <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www] <<"\n with usetracker == false"<<"\n ---> please change to usetracker == true "; } } // tracks with hits only on muon detectors else if (recoTracksLabel=="standAloneMuons" || recoTracksLabel=="standAloneSETMuons" || recoTracksLabel=="cosmicMuons" || recoTracksLabel=="hltL2Muons") { if (usetracker) { edm::LogWarning("MuonTrackValidator") <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www] <<"\n with usetracker == true"<<"\n ---> please change to usetracker == false "; } if (!usemuon) { edm::LogWarning("MuonTrackValidator") <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www] <<"\n with usemuon == false"<<"\n ---> please change to usemuon == true "; } } } // for (unsigned int www=0;www<label.size();www++) }
| virtual MuonTrackValidator::~MuonTrackValidator | ( | ) | [inline, virtual] |
| void MuonTrackValidator::analyze | ( | const edm::Event & | event, |
| const edm::EventSetup & | setup | ||
| ) | [virtual] |
Method called once per event.
Implements edm::EDAnalyzer.
Definition at line 231 of file MuonTrackValidator.cc.
References abs, asciidump::at, funct::cos(), edm::AssociationMap< Tag >::end(), f, edm::AssociationMap< Tag >::find(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), i, edm::Ref< C, T, F >::key(), edm::RefToBase< T >::key(), label, LogTrace, M_PI, min, ParticleBase::momentum(), DetId::Muon, edm::AssociationMap< Tag >::post_insert(), edm::Handle< T >::product(), dt_dqm_sourceclient_common_cff::reco, funct::sin(), mathSSE::sqrt(), tmp, DetId::Tracker, TrackingParticle::trackPSimHit(), ParticleBase::vertex(), and cms::Exception::what().
{
using namespace reco;
edm::LogInfo("MuonTrackValidator") << "\n====================================================" << "\n"
<< "Analyzing new event" << "\n"
<< "====================================================\n" << "\n";
edm::ESHandle<ParametersDefinerForTP> parametersDefinerTP;
setup.get<TrackAssociatorRecord>().get(parametersDefiner,parametersDefinerTP);
edm::Handle<TrackingParticleCollection> TPCollectionHeff ;
event.getByLabel(label_tp_effic,TPCollectionHeff);
const TrackingParticleCollection tPCeff = *(TPCollectionHeff.product());
edm::Handle<TrackingParticleCollection> TPCollectionHfake ;
event.getByLabel(label_tp_fake,TPCollectionHfake);
const TrackingParticleCollection tPCfake = *(TPCollectionHfake.product());
edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
event.getByLabel(bsSrc,recoBeamSpotHandle);
reco::BeamSpot bs = *recoBeamSpotHandle;
int w=0;
for (unsigned int ww=0;ww<associators.size();ww++){
for (unsigned int www=0;www<label.size();www++){
//
//get collections from the event
//
edm::Handle<View<Track> > trackCollection;
reco::RecoToSimCollection recSimColl;
reco::SimToRecoCollection simRecColl;
unsigned int trackCollectionSize = 0;
// if(!event.getByLabel(label[www], trackCollection)&&ignoremissingtkcollection_) continue;
if(!event.getByLabel(label[www], trackCollection)&&ignoremissingtkcollection_) {
recSimColl.post_insert();
simRecColl.post_insert();
}
else {
trackCollectionSize = trackCollection->size();
//associate tracks
if(UseAssociators){
edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
<< label[www].process()<<":"
<< label[www].label()<<":"
<< label[www].instance()<<" with "
<< associators[ww].c_str() <<"\n";
LogTrace("MuonTrackValidator") << "Calling associateRecoToSim method" << "\n";
recSimColl=associator[ww]->associateRecoToSim(trackCollection,
TPCollectionHfake,
&event);
LogTrace("MuonTrackValidator") << "Calling associateSimToReco method" << "\n";
simRecColl=associator[ww]->associateSimToReco(trackCollection,
TPCollectionHeff,
&event);
}
else{
edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
<< label[www].process()<<":"
<< label[www].label()<<":"
<< label[www].instance()<<" with "
<< associatormap.process()<<":"
<< associatormap.label()<<":"
<< associatormap.instance()<<"\n";
Handle<reco::SimToRecoCollection > simtorecoCollectionH;
event.getByLabel(associatormap,simtorecoCollectionH);
simRecColl= *(simtorecoCollectionH.product());
Handle<reco::RecoToSimCollection > recotosimCollectionH;
event.getByLabel(associatormap,recotosimCollectionH);
recSimColl= *(recotosimCollectionH.product());
}
}
//
//fill simulation histograms
//compute number of tracks per eta interval
//
edm::LogVerbatim("MuonTrackValidator") << "\n# of TrackingParticles: " << tPCeff.size() << "\n";
int ats = 0;
int st=0;
for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
bool TP_is_matched = false;
double quality = 0.;
bool Quality05 = false;
bool Quality075 = false;
TrackingParticleRef tpr(TPCollectionHeff, i);
TrackingParticle* tp=const_cast<TrackingParticle*>(tpr.get());
ParticleBase::Vector momentumTP;
ParticleBase::Point vertexTP;
double dxySim = 0;
double dzSim = 0;
//If the TrackingParticle is collison like, get the momentum and vertex at production state
if(parametersDefiner=="LhcParametersDefinerForTP")
{
if(! tpSelector(*tp)) continue;
momentumTP = tp->momentum();
vertexTP = tp->vertex();
//Calcualte the impact parameters w.r.t. PCA
ParticleBase::Vector momentum = parametersDefinerTP->momentum(event,setup,*tp);
ParticleBase::Point vertex = parametersDefinerTP->vertex(event,setup,*tp);
dxySim = (-vertex.x()*sin(momentum.phi())+vertex.y()*cos(momentum.phi()));
dzSim = vertex.z() - (vertex.x()*momentum.x()+vertex.y()*momentum.y())/sqrt(momentum.perp2()) * momentum.z()/sqrt(momentum.perp2());
}
//If the TrackingParticle is comics, get the momentum and vertex at PCA
if(parametersDefiner=="CosmicParametersDefinerForTP")
{
if(! cosmictpSelector(*tp,&bs,event,setup)) continue;
momentumTP = parametersDefinerTP->momentum(event,setup,*tp);
vertexTP = parametersDefinerTP->vertex(event,setup,*tp);
dxySim = (-vertexTP.x()*sin(momentumTP.phi())+vertexTP.y()*cos(momentumTP.phi()));
dzSim = vertexTP.z() - (vertexTP.x()*momentumTP.x()+vertexTP.y()*momentumTP.y())/sqrt(momentumTP.perp2()) * momentumTP.z()/sqrt(momentumTP.perp2());
}
edm::LogVerbatim("MuonTrackValidator") <<"--------------------Selected TrackingParticle #"<<tpr.key();
st++;
h_ptSIM[w]->Fill(sqrt(momentumTP.perp2()));
h_etaSIM[w]->Fill(momentumTP.eta());
h_vertposSIM[w]->Fill(sqrt(vertexTP.perp2()));
std::vector<std::pair<RefToBase<Track>, double> > rt;
if(simRecColl.find(tpr) != simRecColl.end()){
rt = (std::vector<std::pair<RefToBase<Track>, double> >) simRecColl[tpr];
if (rt.size()!=0) {
RefToBase<Track> assoc_recoTrack = rt.begin()->first;
edm::LogVerbatim("MuonTrackValidator")<<"-----------------------------associated Track #"<<assoc_recoTrack.key();
TP_is_matched = true;
ats++;
quality = rt.begin()->second;
edm::LogVerbatim("MuonTrackValidator") << "TrackingParticle #" <<tpr.key()
<< " with pt=" << sqrt(momentumTP.perp2())
<< " associated with quality:" << quality <<"\n";
if (MABH) {
if (quality > 0.75) {
Quality075 = true;
Quality05 = true;
}
else if (quality > 0.5) {
Quality05 = true;
}
}
}
}else{
edm::LogVerbatim("MuonTrackValidator")
<< "TrackingParticle #" << tpr.key()
<< " with pt,eta,phi: "
<< sqrt(momentumTP.perp2()) << " , "
<< momentumTP.eta() << " , "
<< momentumTP.phi() << " , "
<< " NOT associated to any reco::Track" << "\n";
}
for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
if (getEta(momentumTP.eta())>etaintervals[w][f]&&
getEta(momentumTP.eta())<etaintervals[w][f+1]) {
totSIMeta[w][f]++;
if (TP_is_matched) {
totASSeta[w][f]++;
if (MABH) {
if (Quality075) {
totASSeta_Quality075[w][f]++;
totASSeta_Quality05[w][f]++;
}
else if (Quality05) {
totASSeta_Quality05[w][f]++;
}
}
}
}
} // END for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
if (momentumTP.phi() > phiintervals[w][f]&&
momentumTP.phi() <phiintervals[w][f+1]) {
totSIM_phi[w][f]++;
if (TP_is_matched) {
totASS_phi[w][f]++;
if (MABH) {
if (Quality075) {
totASS_phi_Quality075[w][f]++;
totASS_phi_Quality05[w][f]++;
}
else if (Quality05) {
totASS_phi_Quality05[w][f]++;
}
}
}
}
} // END for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
if (getPt(sqrt(momentumTP.perp2()))>pTintervals[w][f]&&
getPt(sqrt(momentumTP.perp2()))<pTintervals[w][f+1]) {
totSIMpT[w][f]++;
if (TP_is_matched) {
totASSpT[w][f]++;
if (MABH) {
if (Quality075) {
totASSpT_Quality075[w][f]++;
totASSpT_Quality05[w][f]++;
}
else if (Quality05) {
totASSpT_Quality05[w][f]++;
}
}
}
}
} // END for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
if (dxySim>dxyintervals[w][f]&&
dxySim<dxyintervals[w][f+1]) {
totSIM_dxy[w][f]++;
if (TP_is_matched) {
totASS_dxy[w][f]++;
}
}
} // END for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
if (dzSim>dzintervals[w][f]&&
dzSim<dzintervals[w][f+1]) {
totSIM_dz[w][f]++;
if (TP_is_matched) {
totASS_dz[w][f]++;
}
}
} // END for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
for (unsigned int f=0; f<vertposintervals[w].size()-1; f++){
if (sqrt(vertexTP.perp2())>vertposintervals[w][f]&&
sqrt(vertexTP.perp2())<vertposintervals[w][f+1]) {
totSIM_vertpos[w][f]++;
if (TP_is_matched) {
totASS_vertpos[w][f]++;
}
}
} // END for (unsigned int f=0; f<vertposintervals[w].size()-1; f++){
for (unsigned int f=0; f<zposintervals[w].size()-1; f++){
if (vertexTP.z()>zposintervals[w][f]&&
vertexTP.z()<zposintervals[w][f+1]) {
totSIM_zpos[w][f]++;
if (TP_is_matched) {
totASS_zpos[w][f]++;
}
}
} // END for (unsigned int f=0; f<zposintervals[w].size()-1; f++){
std::vector<PSimHit> simhits;
if (usetracker && usemuon) {
simhits=tp->trackPSimHit();
}
else if (!usetracker && usemuon) {
simhits=tp->trackPSimHit(DetId::Muon);
}
else if (usetracker && !usemuon) {
simhits=tp->trackPSimHit(DetId::Tracker);
}
int tmp = std::min((int)(simhits.end()-simhits.begin()),int(maxHit-1));
edm::LogVerbatim("MuonTrackValidator") << "\t N simhits = "<< (int)(simhits.end()-simhits.begin())<<"\n";
totSIM_hit[w][tmp]++;
if (TP_is_matched) totASS_hit[w][tmp]++;
if (TP_is_matched)
{
RefToBase<Track> assoctrack = rt.begin()->first;
nrecHit_vs_nsimHit_sim2rec[w]->Fill( assoctrack->numberOfValidHits(),(int)(simhits.end()-simhits.begin() ));
}
} // End for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
if (st!=0) h_tracksSIM[w]->Fill(st);
//
//fill reconstructed track histograms
//
edm::LogVerbatim("MuonTrackValidator") << "\n# of reco::Tracks with "
<< label[www].process()<<":"
<< label[www].label()<<":"
<< label[www].instance()
<< ": " << trackCollectionSize << "\n";
int at=0;
int rT=0;
for(View<Track>::size_type i=0; i<trackCollectionSize; ++i){
bool Track_is_matched = false;
RefToBase<Track> track(trackCollection, i);
rT++;
std::vector<std::pair<TrackingParticleRef, double> > tp;
TrackingParticleRef tpr;
// new logic (bidirectional)
if (BiDirectional_RecoToSim_association) {
edm::LogVerbatim("MuonTrackValidator")<<"----------------------------------------Track #"<< track.key();
if(recSimColl.find(track) != recSimColl.end()) {
tp = recSimColl[track];
if (tp.size() != 0) {
tpr = tp.begin()->first;
// RtS and StR must associate the same pair !
if(simRecColl.find(tpr) != simRecColl.end()) {
std::vector<std::pair<RefToBase<Track>, double> > track_checkback = simRecColl[tpr];
RefToBase<Track> assoc_track_checkback;
assoc_track_checkback = track_checkback.begin()->first;
if ( assoc_track_checkback.key() == track.key() ) {
edm::LogVerbatim("MuonTrackValidator")<<"------------------associated TrackingParticle #"<<tpr.key();
Track_is_matched = true;
at++;
double Purity = tp.begin()->second;
double Quality = track_checkback.begin()->second;
edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
<< " associated with quality:" << Purity <<"\n";
if (MABH) h_PurityVsQuality[w]->Fill(Quality,Purity);
}
}
}
}
if (!Track_is_matched) edm::LogVerbatim("MuonTrackValidator")
<< "reco::Track #" << track.key() << " with pt=" << track->pt() << " NOT associated to any TrackingParticle" << "\n";
}
// old logic (bugged)
else {
if(recSimColl.find(track) != recSimColl.end()){
tp = recSimColl[track];
if (tp.size()!=0) {
Track_is_matched = true;
tpr = tp.begin()->first;
at++;
edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
<< " associated with quality:" << tp.begin()->second <<"\n";
}
} else {
edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
<< " NOT associated to any TrackingParticle" << "\n";
}
}
//Compute fake rate vs eta
for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
if (getEta(track->momentum().eta())>etaintervals[w][f]&&
getEta(track->momentum().eta())<etaintervals[w][f+1]) {
totRECeta[w][f]++;
if (Track_is_matched) {
totASS2eta[w][f]++;
}
}
} // End for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
if (track->momentum().phi()>phiintervals[w][f]&&
track->momentum().phi()<phiintervals[w][f+1]) {
totREC_phi[w][f]++;
if (Track_is_matched) {
totASS2_phi[w][f]++;
}
}
} // End for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
if (getPt(sqrt(track->momentum().perp2()))>pTintervals[w][f]&&
getPt(sqrt(track->momentum().perp2()))<pTintervals[w][f+1]) {
totRECpT[w][f]++;
if (Track_is_matched) {
totASS2pT[w][f]++;
}
}
} // End for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
if (track->dxy(bs.position())>dxyintervals[w][f]&&
track->dxy(bs.position())<dxyintervals[w][f+1]) {
totREC_dxy[w][f]++;
if (Track_is_matched) {
totASS2_dxy[w][f]++;
}
}
} // End for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
if (track->dz(bs.position())>dzintervals[w][f]&&
track->dz(bs.position())<dzintervals[w][f+1]) {
totREC_dz[w][f]++;
if (Track_is_matched) {
totASS2_dz[w][f]++;
}
}
} // End for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
int tmp = std::min((int)track->found(),int(maxHit-1));
totREC_hit[w][tmp]++;
if (Track_is_matched) totASS2_hit[w][tmp]++;
edm::LogVerbatim("MuonTrackValidator") << "\t N valid rechits = "<< (int)track->found() <<"\n";
//Fill other histos
try{
if (!Track_is_matched) continue;
if (associators[ww]=="TrackAssociatorByChi2"){
//association chi2
double assocChi2 = -tp.begin()->second;//in association map is stored -chi2
h_assochi2[www]->Fill(assocChi2);
h_assochi2_prob[www]->Fill(TMath::Prob((assocChi2)*5,5));
}
else if (associators[ww]=="TrackAssociatorByHits"){
double fraction = tp.begin()->second;
h_assocFraction[www]->Fill(fraction);
h_assocSharedHit[www]->Fill(fraction*track->numberOfValidHits());
}
//nchi2 and hits global distributions
h_nchi2[w]->Fill(track->normalizedChi2());
h_nchi2_prob[w]->Fill(TMath::Prob(track->chi2(),(int)track->ndof()));
h_hits[w]->Fill(track->numberOfValidHits());
h_losthits[w]->Fill(track->numberOfLostHits());
chi2_vs_nhits[w]->Fill(track->numberOfValidHits(),track->normalizedChi2());
h_charge[w]->Fill( track->charge() );
//Get tracking particle parameters at point of closest approach to the beamline
ParticleBase::Vector momentumTP = parametersDefinerTP->momentum(event,setup,*(tpr.get()));
ParticleBase::Point vertexTP = parametersDefinerTP->vertex(event,setup,*(tpr.get()));
double ptSim = sqrt(momentumTP.perp2());
double qoverpSim = tpr->charge()/sqrt(momentumTP.x()*momentumTP.x()+momentumTP.y()*momentumTP.y()+momentumTP.z()*momentumTP.z());
double thetaSim = momentumTP.theta();
double lambdaSim = M_PI/2-momentumTP.theta();
double phiSim = momentumTP.phi();
double dxySim = (-vertexTP.x()*sin(momentumTP.phi())+vertexTP.y()*cos(momentumTP.phi()));
double dzSim = vertexTP.z() - (vertexTP.x()*momentumTP.x()+vertexTP.y()*momentumTP.y())/sqrt(momentumTP.perp2()) * momentumTP.z()/sqrt(momentumTP.perp2());
TrackBase::ParameterVector rParameters = track->parameters();
double qoverpRec(0);
double qoverpErrorRec(0);
double ptRec(0);
double ptErrorRec(0);
double lambdaRec(0);
double lambdaErrorRec(0);
double phiRec(0);
double phiErrorRec(0);
//loop to decide whether to take gsfTrack (utilisation of mode-function) or common track
const GsfTrack* gsfTrack(0);
if(useGsf){
gsfTrack = dynamic_cast<const GsfTrack*>(&(*track));
if (gsfTrack==0) edm::LogInfo("MuonTrackValidator") << "Trying to access mode for a non-GsfTrack";
}
if (gsfTrack) {
// get values from mode
getRecoMomentum(*gsfTrack, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
}
else {
// get values from track (without mode)
getRecoMomentum(*track, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
}
double thetaRec = track->theta();
double ptError = ptErrorRec;
double ptres = ptRec - ptSim;
double etares = track->eta()-momentumTP.Eta();
double dxyRec = track->dxy(bs.position());
double dzRec = track->dz(bs.position());
// eta residue; pt, k, theta, phi, dxy, dz pulls
double qoverpPull=(qoverpRec-qoverpSim)/qoverpErrorRec;
double thetaPull=(lambdaRec-lambdaSim)/lambdaErrorRec;
double phiDiff = phiRec - phiSim;
if (abs(phiDiff) > M_PI) {
if (phiDiff >0.) phiDiff = phiDiff - 2.*M_PI;
else phiDiff = phiDiff + 2.*M_PI;
}
double phiPull=phiDiff/phiErrorRec;
double dxyPull=(dxyRec-dxySim)/track->dxyError();
double dzPull=(dzRec-dzSim)/track->dzError();
double contrib_Qoverp = ((qoverpRec-qoverpSim)/qoverpErrorRec)*
((qoverpRec-qoverpSim)/qoverpErrorRec)/5;
double contrib_dxy = ((dxyRec-dxySim)/track->dxyError())*((dxyRec-dxySim)/track->dxyError())/5;
double contrib_dz = ((dzRec-dzSim)/track->dzError())*((dzRec-dzSim)/track->dzError())/5;
double contrib_theta = ((lambdaRec-lambdaSim)/lambdaErrorRec)*
((lambdaRec-lambdaSim)/lambdaErrorRec)/5;
double contrib_phi = (phiDiff/phiErrorRec)*(phiDiff/phiErrorRec)/5;
LogTrace("MuonTrackValidator") << "assocChi2=" << tp.begin()->second << "\n"
<< "" << "\n"
<< "ptREC=" << ptRec << "\n"
<< "etaREC=" << track->eta() << "\n"
<< "qoverpREC=" << qoverpRec << "\n"
<< "dxyREC=" << dxyRec << "\n"
<< "dzREC=" << dzRec << "\n"
<< "thetaREC=" << track->theta() << "\n"
<< "phiREC=" << phiRec << "\n"
<< "" << "\n"
<< "qoverpError()=" << qoverpErrorRec << "\n"
<< "dxyError()=" << track->dxyError() << "\n"
<< "dzError()=" << track->dzError() << "\n"
<< "thetaError()=" << lambdaErrorRec << "\n"
<< "phiError()=" << phiErrorRec << "\n"
<< "" << "\n"
<< "ptSIM=" << ptSim << "\n"
<< "etaSIM=" << momentumTP.Eta() << "\n"
<< "qoverpSIM=" << qoverpSim << "\n"
<< "dxySIM=" << dxySim << "\n"
<< "dzSIM=" << dzSim << "\n"
<< "thetaSIM=" << M_PI/2-lambdaSim << "\n"
<< "phiSIM=" << phiSim << "\n"
<< "" << "\n"
<< "contrib_Qoverp=" << contrib_Qoverp << "\n"
<< "contrib_dxy=" << contrib_dxy << "\n"
<< "contrib_dz=" << contrib_dz << "\n"
<< "contrib_theta=" << contrib_theta << "\n"
<< "contrib_phi=" << contrib_phi << "\n"
<< "" << "\n"
<<"chi2PULL="<<contrib_Qoverp+contrib_dxy+contrib_dz+contrib_theta+contrib_phi<<"\n";
h_pullQoverp[w]->Fill(qoverpPull);
h_pullTheta[w]->Fill(thetaPull);
h_pullPhi[w]->Fill(phiPull);
h_pullDxy[w]->Fill(dxyPull);
h_pullDz[w]->Fill(dzPull);
h_pt[w]->Fill(ptres/ptError);
h_eta[w]->Fill(etares);
etares_vs_eta[w]->Fill(getEta(track->eta()),etares);
//chi2 and #hit vs eta: fill 2D histos
chi2_vs_eta[w]->Fill(getEta(track->eta()),track->normalizedChi2());
nhits_vs_eta[w]->Fill(getEta(track->eta()),track->numberOfValidHits());
nDThits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonDTHits());
nCSChits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonCSCHits());
nRPChits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonRPCHits());
nlosthits_vs_eta[w]->Fill(getEta(track->eta()),track->numberOfLostHits());
//resolution of track params: fill 2D histos
dxyres_vs_eta[w]->Fill(getEta(track->eta()),dxyRec-dxySim);
ptres_vs_eta[w]->Fill(getEta(track->eta()),(ptRec-ptSim)/ptRec);
dzres_vs_eta[w]->Fill(getEta(track->eta()),dzRec-dzSim);
phires_vs_eta[w]->Fill(getEta(track->eta()),phiDiff);
cotThetares_vs_eta[w]->Fill(getEta(track->eta()), cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
//same as before but vs pT
dxyres_vs_pt[w]->Fill(getPt(ptRec),dxyRec-dxySim);
ptres_vs_pt[w]->Fill(getPt(ptRec),(ptRec-ptSim)/ptRec);
dzres_vs_pt[w]->Fill(getPt(ptRec),dzRec-dzSim);
phires_vs_pt[w]->Fill(getPt(ptRec),phiDiff);
cotThetares_vs_pt[w]->Fill(getPt(ptRec), cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
//pulls of track params vs eta: fill 2D histos
dxypull_vs_eta[w]->Fill(getEta(track->eta()),dxyPull);
ptpull_vs_eta[w]->Fill(getEta(track->eta()),ptres/ptError);
dzpull_vs_eta[w]->Fill(getEta(track->eta()),dzPull);
phipull_vs_eta[w]->Fill(getEta(track->eta()),phiPull);
thetapull_vs_eta[w]->Fill(getEta(track->eta()),thetaPull);
//plots vs phi
nhits_vs_phi[w]->Fill(phiRec,track->numberOfValidHits());
chi2_vs_phi[w]->Fill(phiRec,track->normalizedChi2());
ptmean_vs_eta_phi[w]->Fill(phiRec,getEta(track->eta()),ptRec);
phimean_vs_eta_phi[w]->Fill(phiRec,getEta(track->eta()),phiRec);
ptres_vs_phi[w]->Fill(phiRec,(ptRec-ptSim)/ptRec);
phires_vs_phi[w]->Fill(phiRec,phiDiff);
ptpull_vs_phi[w]->Fill(phiRec,ptres/ptError);
phipull_vs_phi[w]->Fill(phiRec,phiPull);
thetapull_vs_phi[w]->Fill(phiRec,thetaPull);
std::vector<PSimHit> simhits;
if (usetracker && usemuon) {
simhits=tpr.get()->trackPSimHit();
}
else if (!usetracker && usemuon) {
simhits=tpr.get()->trackPSimHit(DetId::Muon);
}
else if (usetracker && !usemuon) {
simhits=tpr.get()->trackPSimHit(DetId::Tracker);
}
nrecHit_vs_nsimHit_rec2sim[w]->Fill(track->numberOfValidHits(), (int)(simhits.end()-simhits.begin() ));
} // End of try{
catch (cms::Exception e){
LogTrace("MuonTrackValidator") << "exception found: " << e.what() << "\n";
}
} // End of for(View<Track>::size_type i=0; i<trackCollectionSize; ++i){
if (at!=0) h_tracks[w]->Fill(at);
h_fakes[w]->Fill(rT-at);
edm::LogVerbatim("MuonTrackValidator") << "Total Simulated: " << st << "\n"
<< "Total Associated (simToReco): " << ats << "\n"
<< "Total Reconstructed: " << rT << "\n"
<< "Total Associated (recoToSim): " << at << "\n"
<< "Total Fakes: " << rT-at << "\n";
nrec_vs_nsim[w]->Fill(rT,st);
w++;
} // End of for (unsigned int www=0;www<label.size();www++){
} //END of for (unsigned int ww=0;ww<associators.size();ww++){
}
| void MuonTrackValidator::beginRun | ( | edm::Run const & | , |
| edm::EventSetup const & | setup | ||
| ) | [virtual] |
Method called before the event loop.
this are needed to calculate efficiency during tha harvesting for the automated validation
Reimplemented from edm::EDAnalyzer.
Definition at line 29 of file MuonTrackValidator.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), DQMStore::bookProfile2D(), DQMStore::cd(), dbe_, TrackerOfflineValidation_Dqm_cff::dirName, edm::EventSetup::get(), getTH1F(), getTH2F(), DQMStore::goUp(), edm::InputTag::instance(), j, label, edm::InputTag::label(), max(), min, nint(), edm::InputTag::process(), edm::ESHandle< T >::product(), linker::replace(), and DQMStore::setCurrentFolder().
{
// dbe_->showDirStructure();
int j=0;
for (unsigned int ww=0;ww<associators.size();ww++){
for (unsigned int www=0;www<label.size();www++){
dbe_->cd();
InputTag algo = label[www];
string dirName=dirName_;
if (algo.process()!="")
dirName+=algo.process()+"_";
if(algo.label()!="")
dirName+=algo.label()+"_";
if(algo.instance()!="")
dirName+=algo.instance()+"_";
if (dirName.find("Tracks")<dirName.length()){
dirName.replace(dirName.find("Tracks"),6,"");
}
string assoc= associators[ww];
if (assoc.find("Track")<assoc.length()){
assoc.replace(assoc.find("Track"),5,"");
}
dirName+=assoc;
std::replace(dirName.begin(), dirName.end(), ':', '_');
dbe_->setCurrentFolder(dirName.c_str());
setUpVectors();
dbe_->goUp();
string subDirName = dirName + "/simulation";
dbe_->setCurrentFolder(subDirName.c_str());
h_ptSIM.push_back( dbe_->book1D("ptSIM", "generated p_{t}", 5500, 0, 110 ) );
h_etaSIM.push_back( dbe_->book1D("etaSIM", "generated pseudorapidity", 500, -2.5, 2.5 ) );
h_tracksSIM.push_back( dbe_->book1D("tracksSIM","number of simulated tracks",200,-0.5,99.5) );
h_vertposSIM.push_back( dbe_->book1D("vertposSIM","Transverse position of sim vertices",100,0.,120.) );
dbe_->cd();
dbe_->setCurrentFolder(dirName.c_str());
h_tracks.push_back( dbe_->book1D("tracks","number of reconstructed tracks",200,-0.5,19.5) );
h_fakes.push_back( dbe_->book1D("fakes","number of fake reco tracks",20,-0.5,19.5) );
h_charge.push_back( dbe_->book1D("charge","charge",3,-1.5,1.5) );
h_hits.push_back( dbe_->book1D("hits", "number of hits per track", nintHit,minHit,maxHit ) );
h_losthits.push_back( dbe_->book1D("losthits", "number of lost hits per track", nintHit,minHit,maxHit) );
h_nchi2.push_back( dbe_->book1D("chi2", "normalized #chi^{2}", 200, 0, 20 ) );
h_nchi2_prob.push_back( dbe_->book1D("chi2_prob", "normalized #chi^{2} probability",100,0,1));
h_recoeta.push_back( dbe_->book1D("num_reco_eta","N of reco track vs eta",nint,min,max) );
h_assoceta.push_back( dbe_->book1D("num_assoc(simToReco)_eta","N of associated tracks (simToReco) vs eta",nint,min,max) );
h_assoc2eta.push_back( dbe_->book1D("num_assoc(recoToSim)_eta","N of associated (recoToSim) tracks vs eta",nint,min,max) );
h_simuleta.push_back( dbe_->book1D("num_simul_eta","N of simulated tracks vs eta",nint,min,max) );
h_recopT.push_back( dbe_->book1D("num_reco_pT","N of reco track vs pT",nintpT,minpT,maxpT) );
h_assocpT.push_back( dbe_->book1D("num_assoc(simToReco)_pT","N of associated tracks (simToReco) vs pT",nintpT,minpT,maxpT) );
h_assoc2pT.push_back( dbe_->book1D("num_assoc(recoToSim)_pT","N of associated (recoToSim) tracks vs pT",nintpT,minpT,maxpT) );
h_simulpT.push_back( dbe_->book1D("num_simul_pT","N of simulated tracks vs pT",nintpT,minpT,maxpT) );
//
h_recohit.push_back( dbe_->book1D("num_reco_hit","N of reco track vs hit",nintHit,minHit,maxHit) );
h_assochit.push_back( dbe_->book1D("num_assoc(simToReco)_hit","N of associated tracks (simToReco) vs hit",nintHit,minHit,maxHit) );
h_assoc2hit.push_back( dbe_->book1D("num_assoc(recoToSim)_hit","N of associated (recoToSim) tracks vs hit",nintHit,minHit,maxHit) );
h_simulhit.push_back( dbe_->book1D("num_simul_hit","N of simulated tracks vs hit",nintHit,minHit,maxHit) );
//
h_recophi.push_back( dbe_->book1D("num_reco_phi","N of reco track vs phi",nintPhi,minPhi,maxPhi) );
h_assocphi.push_back( dbe_->book1D("num_assoc(simToReco)_phi","N of associated tracks (simToReco) vs phi",nintPhi,minPhi,maxPhi) );
h_assoc2phi.push_back( dbe_->book1D("num_assoc(recoToSim)_phi","N of associated (recoToSim) tracks vs phi",nintPhi,minPhi,maxPhi) );
h_simulphi.push_back( dbe_->book1D("num_simul_phi","N of simulated tracks vs phi",nintPhi,minPhi,maxPhi) );
h_recodxy.push_back( dbe_->book1D("num_reco_dxy","N of reco track vs dxy",nintDxy,minDxy,maxDxy) );
h_assocdxy.push_back( dbe_->book1D("num_assoc(simToReco)_dxy","N of associated tracks (simToReco) vs dxy",nintDxy,minDxy,maxDxy) );
h_assoc2dxy.push_back( dbe_->book1D("num_assoc(recoToSim)_dxy","N of associated (recoToSim) tracks vs dxy",nintDxy,minDxy,maxDxy) );
h_simuldxy.push_back( dbe_->book1D("num_simul_dxy","N of simulated tracks vs dxy",nintDxy,minDxy,maxDxy) );
h_recodz.push_back( dbe_->book1D("num_reco_dz","N of reco track vs dz",nintDz,minDz,maxDz) );
h_assocdz.push_back( dbe_->book1D("num_assoc(simToReco)_dz","N of associated tracks (simToReco) vs dz",nintDz,minDz,maxDz) );
h_assoc2dz.push_back( dbe_->book1D("num_assoc(recoToSim)_dz","N of associated (recoToSim) tracks vs dz",nintDz,minDz,maxDz) );
h_simuldz.push_back( dbe_->book1D("num_simul_dz","N of simulated tracks vs dz",nintDz,minDz,maxDz) );
h_assocvertpos.push_back( dbe_->book1D("num_assoc(simToReco)_vertpos","N of associated tracks (simToReco) vs transverse vert position",nintVertpos,minVertpos,maxVertpos) );
h_simulvertpos.push_back( dbe_->book1D("num_simul_vertpos","N of simulated tracks vs transverse vert position",nintVertpos,minVertpos,maxVertpos) );
h_assoczpos.push_back( dbe_->book1D("num_assoc(simToReco)_zpos","N of associated tracks (simToReco) vs z vert position",nintZpos,minZpos,maxZpos) );
h_simulzpos.push_back( dbe_->book1D("num_simul_zpos","N of simulated tracks vs z vert position",nintZpos,minZpos,maxZpos) );
h_eta.push_back( dbe_->book1D("eta", "pseudorapidity residue", 1000, -0.1, 0.1 ) );
h_pt.push_back( dbe_->book1D("pullPt", "pull of p_{t}", 100, -10, 10 ) );
h_pullTheta.push_back( dbe_->book1D("pullTheta","pull of #theta parameter",250,-25,25) );
h_pullPhi.push_back( dbe_->book1D("pullPhi","pull of #phi parameter",250,-25,25) );
h_pullDxy.push_back( dbe_->book1D("pullDxy","pull of dxy parameter",250,-25,25) );
h_pullDz.push_back( dbe_->book1D("pullDz","pull of dz parameter",250,-25,25) );
h_pullQoverp.push_back( dbe_->book1D("pullQoverp","pull of qoverp parameter",250,-25,25) );
if (associators[ww]=="TrackAssociatorByChi2"){
h_assochi2.push_back( dbe_->book1D("assocChi2","track association #chi^{2}",1000000,0,100000) );
h_assochi2_prob.push_back(dbe_->book1D("assocChi2_prob","probability of association #chi^{2}",100,0,1));
} else if (associators[ww]=="TrackAssociatorByHits"){
h_assocFraction.push_back( dbe_->book1D("assocFraction","fraction of shared hits",200,0,2) );
h_assocSharedHit.push_back(dbe_->book1D("assocSharedHit","number of shared hits",20,0,20));
}
chi2_vs_nhits.push_back( dbe_->book2D("chi2_vs_nhits","#chi^{2} vs nhits",25,0,25,100,0,10) );
h_chi2meanhitsh.push_back( dbe_->bookProfile("chi2mean_vs_nhits","mean #chi^{2} vs nhits",25,0,25,100,0,10) );
etares_vs_eta.push_back( dbe_->book2D("etares_vs_eta","etaresidue vs eta",nint,min,max,200,-0.1,0.1) );
nrec_vs_nsim.push_back( dbe_->book2D("nrec_vs_nsim","nrec vs nsim",20,-0.5,19.5,20,-0.5,19.5) );
chi2_vs_eta.push_back( dbe_->book2D("chi2_vs_eta","chi2_vs_eta",nint,min,max, 200, 0, 20 ));
h_chi2meanh.push_back( dbe_->bookProfile("chi2mean","mean #chi^{2} vs #eta",nint,min,max, 200, 0, 20) );
chi2_vs_phi.push_back( dbe_->book2D("chi2_vs_phi","#chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20 ) );
h_chi2mean_vs_phi.push_back( dbe_->bookProfile("chi2mean_vs_phi","mean of #chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20) );
nhits_vs_eta.push_back( dbe_->book2D("nhits_vs_eta","nhits vs eta",nint,min,max,nintHit,minHit,maxHit) );
nDThits_vs_eta.push_back( dbe_->book2D("nDThits_vs_eta","# DT hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
nCSChits_vs_eta.push_back( dbe_->book2D("nCSChits_vs_eta","# CSC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
nRPChits_vs_eta.push_back( dbe_->book2D("nRPChits_vs_eta","# RPC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
h_DThits_eta.push_back( dbe_->bookProfile("DThits_eta","mean # DT hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
h_CSChits_eta.push_back( dbe_->bookProfile("CSChits_eta","mean # CSC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
h_RPChits_eta.push_back( dbe_->bookProfile("RPChits_eta","mean # RPC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
h_hits_eta.push_back( dbe_->bookProfile("hits_eta","mean #hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
nhits_vs_phi.push_back( dbe_->book2D("nhits_vs_phi","#hits vs #phi",nintPhi,minPhi,maxPhi,nintHit,minHit,maxHit) );
h_hits_phi.push_back( dbe_->bookProfile("hits_phi","mean #hits vs #phi",nintPhi,minPhi,maxPhi, nintHit,minHit,maxHit) );
nlosthits_vs_eta.push_back( dbe_->book2D("nlosthits_vs_eta","nlosthits vs eta",nint,min,max,nintHit,minHit,maxHit) );
h_losthits_eta.push_back( dbe_->bookProfile("losthits_eta","losthits_eta",nint,min,max,nintHit,minHit,maxHit) );
ptres_vs_eta.push_back(dbe_->book2D("ptres_vs_eta","ptres_vs_eta",nint,min,max, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
ptres_vs_phi.push_back( dbe_->book2D("ptres_vs_phi","p_{t} res vs #phi",nintPhi,minPhi,maxPhi, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
ptres_vs_pt.push_back(dbe_->book2D("ptres_vs_pt","ptres_vs_pt",nintpT,minpT,maxpT, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
cotThetares_vs_eta.push_back(dbe_->book2D("cotThetares_vs_eta","cotThetares_vs_eta",nint,min,max,cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
cotThetares_vs_pt.push_back(dbe_->book2D("cotThetares_vs_pt","cotThetares_vs_pt",nintpT,minpT,maxpT, cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
phires_vs_eta.push_back(dbe_->book2D("phires_vs_eta","phires_vs_eta",nint,min,max, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
phires_vs_pt.push_back(dbe_->book2D("phires_vs_pt","phires_vs_pt",nintpT,minpT,maxpT, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
phires_vs_phi.push_back(dbe_->book2D("phires_vs_phi","#phi res vs #phi",nintPhi,minPhi,maxPhi,phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
dxyres_vs_eta.push_back(dbe_->book2D("dxyres_vs_eta","dxyres_vs_eta",nint,min,max,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
dxyres_vs_pt.push_back( dbe_->book2D("dxyres_vs_pt","dxyres_vs_pt",nintpT,minpT,maxpT,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
dzres_vs_eta.push_back(dbe_->book2D("dzres_vs_eta","dzres_vs_eta",nint,min,max,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
dzres_vs_pt.push_back(dbe_->book2D("dzres_vs_pt","dzres_vs_pt",nintpT,minpT,maxpT,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
ptmean_vs_eta_phi.push_back(dbe_->bookProfile2D("ptmean_vs_eta_phi","mean p_{t} vs #eta and #phi",nintPhi,minPhi,maxPhi,nint,min,max,1000,0,1000));
phimean_vs_eta_phi.push_back(dbe_->bookProfile2D("phimean_vs_eta_phi","mean #phi vs #eta and #phi",nintPhi,minPhi,maxPhi,nint,min,max,nintPhi,minPhi,maxPhi));
//pulls of track params vs eta: to be used with fitslicesytool
dxypull_vs_eta.push_back(dbe_->book2D("dxypull_vs_eta","dxypull_vs_eta",nint,min,max,100,-10,10));
ptpull_vs_eta.push_back(dbe_->book2D("ptpull_vs_eta","ptpull_vs_eta",nint,min,max,100,-10,10));
dzpull_vs_eta.push_back(dbe_->book2D("dzpull_vs_eta","dzpull_vs_eta",nint,min,max,100,-10,10));
phipull_vs_eta.push_back(dbe_->book2D("phipull_vs_eta","phipull_vs_eta",nint,min,max,100,-10,10));
thetapull_vs_eta.push_back(dbe_->book2D("thetapull_vs_eta","thetapull_vs_eta",nint,min,max,100,-10,10));
//pulls of track params vs phi
ptpull_vs_phi.push_back(dbe_->book2D("ptpull_vs_phi","p_{t} pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
phipull_vs_phi.push_back(dbe_->book2D("phipull_vs_phi","#phi pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
thetapull_vs_phi.push_back(dbe_->book2D("thetapull_vs_phi","#theta pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
nrecHit_vs_nsimHit_sim2rec.push_back( dbe_->book2D("nrecHit_vs_nsimHit_sim2rec","nrecHit vs nsimHit (Sim2RecAssoc)",nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
nrecHit_vs_nsimHit_rec2sim.push_back( dbe_->book2D("nrecHit_vs_nsimHit_rec2sim","nrecHit vs nsimHit (Rec2simAssoc)",nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
if (MABH) {
h_PurityVsQuality.push_back( dbe_->book2D("PurityVsQuality","Purity vs Quality (MABH)",20,0.01,1.01,20,0.01,1.01) );
h_assoceta_Quality05.push_back( dbe_->book1D("num_assoc(simToReco)_eta_Q05","N of associated tracks (simToReco) vs eta (Quality>0.5)",nint,min,max) );
h_assoceta_Quality075.push_back( dbe_->book1D("num_assoc(simToReco)_eta_Q075","N of associated tracks (simToReco) vs eta (Quality>0.75)",nint,min,max) );
h_assocpT_Quality05.push_back( dbe_->book1D("num_assoc(simToReco)_pT_Q05","N of associated tracks (simToReco) vs pT (Quality>0.5)",nintpT,minpT,maxpT) );
h_assocpT_Quality075.push_back( dbe_->book1D("num_assoc(simToReco)_pT_Q075","N of associated tracks (simToReco) vs pT (Quality>0.75)",nintpT,minpT,maxpT) );
h_assocphi_Quality05.push_back( dbe_->book1D("num_assoc(simToReco)_phi_Q05","N of associated tracks (simToReco) vs phi (Quality>0.5)",nintPhi,minPhi,maxPhi) );
h_assocphi_Quality075.push_back( dbe_->book1D("num_assoc(simToReco)_phi_Q075","N of associated tracks (simToReco) vs phi (Quality>0.75)",nintPhi,minPhi,maxPhi) );
}
if(useLogPt){
BinLogX(dzres_vs_pt[j]->getTH2F());
BinLogX(dxyres_vs_pt[j]->getTH2F());
BinLogX(phires_vs_pt[j]->getTH2F());
BinLogX(cotThetares_vs_pt[j]->getTH2F());
BinLogX(ptres_vs_pt[j]->getTH2F());
BinLogX(h_recopT[j]->getTH1F());
BinLogX(h_assocpT[j]->getTH1F());
BinLogX(h_assoc2pT[j]->getTH1F());
BinLogX(h_simulpT[j]->getTH1F());
if (MABH) {
BinLogX(h_assocpT_Quality05[j]->getTH1F());
BinLogX(h_assocpT_Quality075[j]->getTH1F());
}
j++;
}
}
}
if (UseAssociators) {
edm::ESHandle<TrackAssociatorBase> theAssociator;
for (unsigned int w=0;w<associators.size();w++) {
setup.get<TrackAssociatorRecord>().get(associators[w],theAssociator);
associator.push_back( theAssociator.product() );
}
}
}
| void MuonTrackValidator::endRun | ( | edm::Run const & | , |
| edm::EventSetup const & | |||
| ) | [virtual] |
Method called at the end of the event loop.
Reimplemented from edm::EDAnalyzer.
Definition at line 854 of file MuonTrackValidator.cc.
References dbe_, label, dbtoconf::out, and DQMStore::save().
{
int w=0;
for (unsigned int ww=0;ww<associators.size();ww++){
for (unsigned int www=0;www<label.size();www++){
//chi2 and #hit vs eta: get mean from 2D histos
doProfileX(chi2_vs_eta[w],h_chi2meanh[w]);
doProfileX(nhits_vs_eta[w],h_hits_eta[w]);
doProfileX(nDThits_vs_eta[w],h_DThits_eta[w]);
doProfileX(nCSChits_vs_eta[w],h_CSChits_eta[w]);
doProfileX(nRPChits_vs_eta[w],h_RPChits_eta[w]);
doProfileX(nlosthits_vs_eta[w],h_losthits_eta[w]);
//vs phi
doProfileX(chi2_vs_nhits[w],h_chi2meanhitsh[w]);
doProfileX(chi2_vs_phi[w],h_chi2mean_vs_phi[w]);
doProfileX(nhits_vs_phi[w],h_hits_phi[w]);
fillPlotFromVector(h_recoeta[w],totRECeta[w]);
fillPlotFromVector(h_simuleta[w],totSIMeta[w]);
fillPlotFromVector(h_assoceta[w],totASSeta[w]);
fillPlotFromVector(h_assoc2eta[w],totASS2eta[w]);
fillPlotFromVector(h_recopT[w],totRECpT[w]);
fillPlotFromVector(h_simulpT[w],totSIMpT[w]);
fillPlotFromVector(h_assocpT[w],totASSpT[w]);
fillPlotFromVector(h_assoc2pT[w],totASS2pT[w]);
fillPlotFromVector(h_recohit[w],totREC_hit[w]);
fillPlotFromVector(h_simulhit[w],totSIM_hit[w]);
fillPlotFromVector(h_assochit[w],totASS_hit[w]);
fillPlotFromVector(h_assoc2hit[w],totASS2_hit[w]);
fillPlotFromVector(h_recophi[w],totREC_phi[w]);
fillPlotFromVector(h_simulphi[w],totSIM_phi[w]);
fillPlotFromVector(h_assocphi[w],totASS_phi[w]);
fillPlotFromVector(h_assoc2phi[w],totASS2_phi[w]);
fillPlotFromVector(h_recodxy[w],totREC_dxy[w]);
fillPlotFromVector(h_simuldxy[w],totSIM_dxy[w]);
fillPlotFromVector(h_assocdxy[w],totASS_dxy[w]);
fillPlotFromVector(h_assoc2dxy[w],totASS2_dxy[w]);
fillPlotFromVector(h_recodz[w],totREC_dz[w]);
fillPlotFromVector(h_simuldz[w],totSIM_dz[w]);
fillPlotFromVector(h_assocdz[w],totASS_dz[w]);
fillPlotFromVector(h_assoc2dz[w],totASS2_dz[w]);
fillPlotFromVector(h_simulvertpos[w],totSIM_vertpos[w]);
fillPlotFromVector(h_assocvertpos[w],totASS_vertpos[w]);
fillPlotFromVector(h_simulzpos[w],totSIM_zpos[w]);
fillPlotFromVector(h_assoczpos[w],totASS_zpos[w]);
if (MABH) {
fillPlotFromVector(h_assoceta_Quality05[w] ,totASSeta_Quality05[w]);
fillPlotFromVector(h_assoceta_Quality075[w],totASSeta_Quality075[w]);
fillPlotFromVector(h_assocpT_Quality05[w] ,totASSpT_Quality05[w]);
fillPlotFromVector(h_assocpT_Quality075[w],totASSpT_Quality075[w]);
fillPlotFromVector(h_assocphi_Quality05[w] ,totASS_phi_Quality05[w]);
fillPlotFromVector(h_assocphi_Quality075[w],totASS_phi_Quality075[w]);
}
w++;
}
}
if ( out.size() != 0 && dbe_ ) dbe_->save(out);
}
| void MuonTrackValidator::getRecoMomentum | ( | const reco::Track & | track, |
| double & | pt, | ||
| double & | ptError, | ||
| double & | qoverp, | ||
| double & | qoverpError, | ||
| double & | lambda, | ||
| double & | lambdaError, | ||
| double & | phi, | ||
| double & | phiError | ||
| ) | const [private] |
retrieval of reconstructed momentum components from reco::Track (== mean values for GSF)
Definition at line 927 of file MuonTrackValidator.cc.
References reco::TrackBase::lambda(), reco::TrackBase::lambdaError(), reco::TrackBase::phi(), reco::TrackBase::phiError(), reco::TrackBase::pt(), reco::TrackBase::ptError(), reco::TrackBase::qoverp(), and reco::TrackBase::qoverpError().
{
pt = track.pt();
ptError = track.ptError();
qoverp = track.qoverp();
qoverpError = track.qoverpError();
lambda = track.lambda();
lambdaError = track.lambdaError();
phi = track.phi();
phiError = track.phiError();
}
| void MuonTrackValidator::getRecoMomentum | ( | const reco::GsfTrack & | gsfTrack, |
| double & | pt, | ||
| double & | ptError, | ||
| double & | qoverp, | ||
| double & | qoverpError, | ||
| double & | lambda, | ||
| double & | lambdaError, | ||
| double & | phi, | ||
| double & | phiError | ||
| ) | const [private] |
retrieval of reconstructed momentum components based on the mode of a reco::GsfTrack
Definition at line 941 of file MuonTrackValidator.cc.
References reco::GsfTrack::lambdaMode(), reco::GsfTrack::lambdaModeError(), reco::GsfTrack::phiMode(), reco::GsfTrack::phiModeError(), reco::GsfTrack::ptMode(), reco::GsfTrack::ptModeError(), reco::GsfTrack::qoverpMode(), and reco::GsfTrack::qoverpModeError().
{
pt = gsfTrack.ptMode();
ptError = gsfTrack.ptModeError();
qoverp = gsfTrack.qoverpMode();
qoverpError = gsfTrack.qoverpModeError();
lambda = gsfTrack.lambdaMode();
lambdaError = gsfTrack.lambdaModeError();
phi = gsfTrack.phiMode();
phiError = gsfTrack.phiModeError();
}
Definition at line 139 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
bool MuonTrackValidator::BiDirectional_RecoToSim_association [private] |
Definition at line 150 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
std::vector<MonitorElement*> MuonTrackValidator::chi2_vs_eta [private] |
Definition at line 169 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::chi2_vs_nhits [private] |
Definition at line 158 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::chi2_vs_phi [private] |
Definition at line 160 of file MuonTrackValidator.h.
Definition at line 147 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
std::vector<MonitorElement*> MuonTrackValidator::cotThetares_vs_eta [private] |
Definition at line 175 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::cotThetares_vs_pt [private] |
Definition at line 176 of file MuonTrackValidator.h.
std::string MuonTrackValidator::dirName_ [private] |
Definition at line 138 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
std::vector<MonitorElement*> MuonTrackValidator::dxypull_vs_eta [private] |
Definition at line 179 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::dxyres_vs_eta [private] |
Definition at line 175 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::dxyres_vs_pt [private] |
Definition at line 176 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::dzpull_vs_eta [private] |
Definition at line 179 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::dzres_vs_eta [private] |
Definition at line 175 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::dzres_vs_pt [private] |
Definition at line 176 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::etares_vs_eta [private] |
Definition at line 158 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_assochi2 [private] |
Definition at line 166 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_assochi2_prob [private] |
Definition at line 166 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_chi2mean_vs_phi [private] |
Definition at line 172 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_chi2meanh [private] |
Definition at line 170 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_chi2meanhitsh [private] |
Definition at line 172 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_dxypulleta [private] |
Definition at line 181 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_dzpulleta [private] |
Definition at line 181 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_hits_phi [private] |
Definition at line 171 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_losthits [private] |
Definition at line 155 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_losthits_eta [private] |
Definition at line 170 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_nchi2 [private] |
Definition at line 155 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_nchi2_prob [private] |
Definition at line 155 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_phipulleta [private] |
Definition at line 181 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_phipullphi [private] |
Definition at line 182 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_ptpulleta [private] |
Definition at line 181 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_ptpullphi [private] |
Definition at line 182 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_ptshifteta [private] |
Definition at line 159 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_thetapulleta [private] |
Definition at line 181 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::h_thetapullphi [private] |
Definition at line 182 of file MuonTrackValidator.h.
bool MuonTrackValidator::MABH [private] |
Definition at line 152 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
double MuonTrackValidator::maxPhi [private] |
Reimplemented from MuonTrackValidatorBase.
Definition at line 141 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
double MuonTrackValidator::minPhi [private] |
Reimplemented from MuonTrackValidatorBase.
Definition at line 141 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
std::vector<MonitorElement*> MuonTrackValidator::nhits_vs_phi [private] |
Definition at line 160 of file MuonTrackValidator.h.
int MuonTrackValidator::nintPhi [private] |
Reimplemented from MuonTrackValidatorBase.
Definition at line 142 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
std::vector<MonitorElement*> MuonTrackValidator::nlosthits_vs_eta [private] |
Definition at line 169 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::phimean_vs_eta_phi [private] |
Definition at line 163 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::phipull_vs_eta [private] |
Definition at line 179 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::phipull_vs_phi [private] |
Definition at line 180 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::phires_vs_eta [private] |
Definition at line 175 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::phires_vs_phi [private] |
Definition at line 160 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::phires_vs_pt [private] |
Definition at line 176 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::ptmean_vs_eta_phi [private] |
Definition at line 163 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::ptpull_vs_eta [private] |
Definition at line 179 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::ptpull_vs_phi [private] |
Definition at line 180 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::ptres_vs_eta [private] |
Definition at line 175 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::ptres_vs_phi [private] |
Definition at line 160 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::ptres_vs_pt [private] |
Definition at line 176 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::thetapull_vs_eta [private] |
Definition at line 179 of file MuonTrackValidator.h.
std::vector<MonitorElement*> MuonTrackValidator::thetapull_vs_phi [private] |
Definition at line 180 of file MuonTrackValidator.h.
Definition at line 146 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
bool MuonTrackValidator::UseAssociators [private] |
Definition at line 140 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
bool MuonTrackValidator::useGsf [private] |
Definition at line 143 of file MuonTrackValidator.h.
Referenced by MuonTrackValidator().
1.7.3