|
|
|
#include <MuonTrackAnalyzer.h>
Analyzer of the Muon tracks
Analyzer of the StandAlone muon tracks
Definition at line 40 of file MuonTrackAnalyzer.h.
Definition at line 43 of file MuonTrackAnalyzer.h.
| MuonTrackAnalyzer::MuonTrackAnalyzer | ( | const edm::ParameterSet & | pset | ) |
Constructor.
Definition at line 49 of file MuonTrackAnalyzer.cc.
References dbe_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), MuonServiceProxy_cff::MuonServiceProxy, MuonUpdatorAtVertex_cff::MuonUpdatorAtVertex, cmsCodeRules::cppFunctionSkipper::operator, and dbtoconf::out.
{
// service parameters
ParameterSet serviceParameters = pset.getParameter<ParameterSet>("ServiceParameters");
// the services
theService = new MuonServiceProxy(serviceParameters);
theTracksLabel = pset.getParameter<InputTag>("Tracks");
doTracksAnalysis = pset.getUntrackedParameter<bool>("DoTracksAnalysis",true);
doSeedsAnalysis = pset.getUntrackedParameter<bool>("DoSeedsAnalysis",false);
if(doSeedsAnalysis){
theSeedsLabel = pset.getParameter<InputTag>("MuonSeed");
ParameterSet updatorPar = pset.getParameter<ParameterSet>("MuonUpdatorAtVertexParameters");
theSeedPropagatorName = updatorPar.getParameter<string>("Propagator");
theUpdator = new MuonUpdatorAtVertex(updatorPar,theService);
}
theCSCSimHitLabel = pset.getParameter<InputTag>("CSCSimHit");
theDTSimHitLabel = pset.getParameter<InputTag>("DTSimHit");
theRPCSimHitLabel = pset.getParameter<InputTag>("RPCSimHit");
theEtaRange = (EtaRange) pset.getParameter<int>("EtaRange");
// number of sim tracks
numberOfSimTracks=0;
// number of reco tracks
numberOfRecTracks=0;
dbe_ = edm::Service<DQMStore>().operator->();
out = pset.getUntrackedParameter<string>("rootFileName");
dirName_ = pset.getUntrackedParameter<std::string>("dirName");
}
| MuonTrackAnalyzer::~MuonTrackAnalyzer | ( | ) | [virtual] |
Destructor.
Definition at line 86 of file MuonTrackAnalyzer.cc.
{
if (theService) delete theService;
}
| void MuonTrackAnalyzer::analyze | ( | const edm::Event & | event, |
| const edm::EventSetup & | eventSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 190 of file MuonTrackAnalyzer.cc.
References LogDebug.
{
LogDebug("MuonTrackAnalyzer") << "Run: " << event.id().run() << " Event: " << event.id().event();
// Update the services
theService->update(eventSetup);
Handle<SimTrackContainer> simTracks;
event.getByLabel("g4SimHits",simTracks);
fillPlots(event,simTracks);
if(doTracksAnalysis)
tracksAnalysis(event,eventSetup,simTracks);
if(doSeedsAnalysis)
seedsAnalysis(event,eventSetup,simTracks);
}
| void MuonTrackAnalyzer::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 90 of file MuonTrackAnalyzer.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::cd(), dbe_, TrackerOfflineValidation_Dqm_cff::dirName, hChi2, edm::InputTag::instance(), edm::InputTag::label(), edm::InputTag::process(), linker::replace(), DQMStore::setCurrentFolder(), and DQMStore::showDirStructure().
{
dbe_->showDirStructure();
dbe_->cd();
InputTag algo = theTracksLabel;
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,"");
}
std::replace(dirName.begin(), dirName.end(), ':', '_');
dbe_->setCurrentFolder(dirName.c_str());
//dbe_->goUp();
std::string simName = dirName;
simName+="/SimTracks";
hSimTracks = new HTrackVariables(simName.c_str(),"SimTracks");
dbe_->cd();
dbe_->setCurrentFolder(dirName.c_str());
// Create the root file
//theFile = new TFile(theRootFileName.c_str(), "RECREATE");
if(doSeedsAnalysis){
dbe_->cd();
dbe_->setCurrentFolder(dirName.c_str());
hRecoSeedInner = new HTrack(dirName.c_str(),"RecoSeed","Inner");
hRecoSeedPCA = new HTrack(dirName.c_str(),"RecoSeed","PCA");
}
if(doTracksAnalysis){
dbe_->cd();
dbe_->setCurrentFolder(dirName.c_str());
hRecoTracksPCA = new HTrack(dirName.c_str(),"RecoTracks","PCA");
hRecoTracksInner = new HTrack(dirName.c_str(),"RecoTracks","Inner");
hRecoTracksOuter = new HTrack(dirName.c_str(),"RecoTracks","Outer");
dbe_->cd();
dbe_->setCurrentFolder(dirName.c_str());
// General Histos
hChi2 = dbe_->book1D("chi2","#chi^2",200,0,200);
hChi2VsEta = dbe_->book2D("chi2VsEta","#chi^2 VS #eta",120,-3.,3.,200,0,200);
hChi2Norm = dbe_->book1D("chi2Norm","Normalized #chi^2",400,0,100);
hChi2NormVsEta = dbe_->book2D("chi2NormVsEta","Normalized #chi^2 VS #eta",120,-3.,3.,400,0,100);
hHitsPerTrack = dbe_->book1D("HitsPerTrack","Number of hits per track",55,0,55);
hHitsPerTrackVsEta = dbe_->book2D("HitsPerTrackVsEta","Number of hits per track VS #eta",
120,-3.,3.,55,0,55);
hDof = dbe_->book1D("dof","Number of Degree of Freedom",55,0,55);
hDofVsEta = dbe_->book2D("dofVsEta","Number of Degree of Freedom VS #eta",120,-3.,3.,55,0,55);
hChi2Prob = dbe_->book1D("chi2Prob","#chi^2 probability",200,0,1);
hChi2ProbVsEta = dbe_->book2D("chi2ProbVsEta","#chi^2 probability VS #eta",120,-3.,3.,200,0,1);
hNumberOfTracks = dbe_->book1D("NumberOfTracks","Number of reconstructed tracks per event",200,0,200);
hNumberOfTracksVsEta = dbe_->book2D("NumberOfTracksVsEta",
"Number of reconstructed tracks per event VS #eta",
120,-3.,3.,10,0,10);
hChargeVsEta = dbe_->book2D("ChargeVsEta","Charge vs #eta gen",120,-3.,3.,4,-2.,2.);
hChargeVsPt = dbe_->book2D("ChargeVsPt","Charge vs P_{T} gen",250,0,200,4,-2.,2.);
hPtRecVsPtGen = dbe_->book2D("PtRecVsPtGen","P_{T} rec vs P_{T} gen",250,0,200,250,0,200);
hDeltaPtVsEta = dbe_->book2D("DeltaPtVsEta","#Delta P_{t} vs #eta gen",120,-3.,3.,500,-250.,250.);
hDeltaPt_In_Out_VsEta = dbe_->book2D("DeltaPt_In_Out_VsEta_","P^{in}_{t} - P^{out}_{t} vs #eta gen",120,-3.,3.,500,-250.,250.);
}
//theFile->cd();
}
| bool MuonTrackAnalyzer::checkMuonSimHitPresence | ( | const edm::Event & | event, |
| edm::Handle< edm::SimTrackContainer > | simTracks | ||
| ) | [private] |
Definition at line 432 of file MuonTrackAnalyzer.cc.
References abs.
{
// Get the SimHit collection from the event
Handle<PSimHitContainer> dtSimHits;
event.getByLabel(theDTSimHitLabel.instance(),theDTSimHitLabel.label(), dtSimHits);
Handle<PSimHitContainer> cscSimHits;
event.getByLabel(theCSCSimHitLabel.instance(),theCSCSimHitLabel.label(), cscSimHits);
Handle<PSimHitContainer> rpcSimHits;
event.getByLabel(theRPCSimHitLabel.instance(),theRPCSimHitLabel.label(), rpcSimHits);
map<unsigned int, vector<const PSimHit*> > mapOfMuonSimHits;
for(PSimHitContainer::const_iterator simhit = dtSimHits->begin();
simhit != dtSimHits->end(); ++simhit) {
if (abs(simhit->particleType()) != 13) continue;
mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
}
for(PSimHitContainer::const_iterator simhit = cscSimHits->begin();
simhit != cscSimHits->end(); ++simhit) {
if (abs(simhit->particleType()) != 13) continue;
mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
}
for(PSimHitContainer::const_iterator simhit = rpcSimHits->begin();
simhit != rpcSimHits->end(); ++simhit) {
if (abs(simhit->particleType()) != 13) continue;
mapOfMuonSimHits[simhit->trackId()].push_back(&*simhit);
}
bool presence = false;
for (SimTrackContainer::const_iterator simTrack = simTracks->begin();
simTrack != simTracks->end(); ++simTrack){
if (abs(simTrack->type()) != 13) continue;
map<unsigned int, vector<const PSimHit*> >::const_iterator mapIterator =
mapOfMuonSimHits.find(simTrack->trackId());
if (mapIterator != mapOfMuonSimHits.end())
presence = true;
}
return presence;
}
| void MuonTrackAnalyzer::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 171 of file MuonTrackAnalyzer.cc.
References dbe_, dbtoconf::out, and DQMStore::save().
{
LogInfo("MuonTrackAnalyzer")<< "Number of Sim tracks: " << numberOfSimTracks;
LogInfo("MuonTrackAnalyzer") << "Number of Reco tracks: " << numberOfRecTracks;
if(doTracksAnalysis){
double eff = hRecoTracksPCA->computeEfficiency(hSimTracks);
LogInfo("MuonTrackAnalyzer") <<" *Track Efficiency* = "<< eff <<"%";
}
if(doSeedsAnalysis){
double eff = hRecoSeedInner->computeEfficiency(hSimTracks);
LogInfo("MuonTrackAnalyzer")<<" *Seed Efficiency* = "<< eff <<"%";
}
if ( out.size() != 0 && dbe_ ) dbe_->save(out);
}
| void MuonTrackAnalyzer::fillPlots | ( | TrajectoryStateOnSurface & | recoTSOS, |
| SimTrack & | simState, | ||
| HTrack * | histo, | ||
| MuonPatternRecoDumper & | debug | ||
| ) | [private] |
Definition at line 345 of file MuonTrackAnalyzer.cc.
References HTrack::computeResolutionAndPull(), MuonPatternRecoDumper::dumpTSOS(), HTrack::Fill(), HTrack::FillDeltaR(), TrajectoryStateOnSurface::globalMomentum(), LogTrace, CoreSimTrack::momentum(), dt_dqm_sourceclient_common_cff::reco, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
LogTrace("MuonTrackAnalyzer") << debug.dumpTSOS(recoTSOS)<<endl;
histo->Fill(recoTSOS);
GlobalVector tsosVect = recoTSOS.globalMomentum();
math::XYZVectorD reco(tsosVect.x(), tsosVect.y(), tsosVect.z());
double deltaRVal = deltaR<double>(reco.eta(),reco.phi(),
simTrack.momentum().eta(),simTrack.momentum().phi());
histo->FillDeltaR(deltaRVal);
histo->computeResolutionAndPull(recoTSOS,simTrack);
}
| void MuonTrackAnalyzer::fillPlots | ( | const edm::Event & | event, |
| edm::Handle< edm::SimTrackContainer > & | simTracks | ||
| ) | [private] |
Definition at line 296 of file MuonTrackAnalyzer.cc.
References abs, LogTrace, and mathSSE::sqrt().
{
if(!checkMuonSimHitPresence(event,simTracks)) return;
// Loop over the Sim tracks
SimTrackContainer::const_iterator simTrack;
LogTrace("MuonTrackAnalyzer")<<"Simulated tracks: "<<simTracks->size()<<endl;
for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
if (abs((*simTrack).type()) == 13) {
if( !isInTheAcceptance( (*simTrack).momentum().eta()) ) continue; // FIXME!!
numberOfSimTracks++;
LogTrace("MuonTrackAnalyzer")<<"Simualted muon:"<<endl;
LogTrace("MuonTrackAnalyzer")<<"Sim pT: "<<sqrt((*simTrack).momentum().perp2())<<endl;
LogTrace("MuonTrackAnalyzer")<<"Sim Eta: "<<(*simTrack).momentum().eta()<<endl; // FIXME
hSimTracks->Fill((*simTrack).momentum().mag(),
sqrt((*simTrack).momentum().perp2()),
(*simTrack).momentum().eta(),
(*simTrack).momentum().phi(),
-(*simTrack).type()/ abs((*simTrack).type()) ); // Double FIXME
LogTrace("MuonTrackAnalyzer") << "hSimTracks filled" << endl;
}
LogTrace("MuonTrackAnalyzer") << endl;
}
| void MuonTrackAnalyzer::fillPlots | ( | reco::TransientTrack & | track, |
| SimTrack & | simTrack | ||
| ) | [private] |
Definition at line 327 of file MuonTrackAnalyzer.cc.
References reco::TransientTrack::chi2(), ChiSquaredProbability(), hChi2, LogTrace, CoreSimTrack::momentum(), reco::TransientTrack::ndof(), reco::TransientTrack::normalizedChi2(), and reco::TransientTrack::recHitsSize().
{
LogTrace("MuonTrackAnalyzer")<<"Analizer: New track, chi2: "<<track.chi2()<<" dof: "<<track.ndof()<<endl;
hChi2->Fill(track.chi2());
hDof->Fill(track.ndof());
hChi2Norm->Fill(track.normalizedChi2());
hHitsPerTrack->Fill(track.recHitsSize());
hChi2Prob->Fill( ChiSquaredProbability(track.chi2(),track.ndof()) );
hChi2VsEta->Fill(simTrack.momentum().eta(),track.chi2());
hChi2NormVsEta->Fill(simTrack.momentum().eta(),track.normalizedChi2());
hChi2ProbVsEta->Fill(simTrack.momentum().eta(),ChiSquaredProbability(track.chi2(),track.ndof()));
hHitsPerTrackVsEta->Fill(simTrack.momentum().eta(),track.recHitsSize());
hDofVsEta->Fill(simTrack.momentum().eta(),track.ndof());
}
| void MuonTrackAnalyzer::fillPlots | ( | FreeTrajectoryState & | recoFTS, |
| SimTrack & | simTrack, | ||
| HTrack * | histo, | ||
| MuonPatternRecoDumper & | debug | ||
| ) | [private] |
Definition at line 361 of file MuonTrackAnalyzer.cc.
References HTrack::computeResolutionAndPull(), MuonPatternRecoDumper::dumpFTS(), HTrack::Fill(), HTrack::FillDeltaR(), LogTrace, CoreSimTrack::momentum(), FreeTrajectoryState::momentum(), dt_dqm_sourceclient_common_cff::reco, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
LogTrace("MuonTrackAnalyzer") << debug.dumpFTS(recoFTS)<<endl;
histo->Fill(recoFTS);
GlobalVector ftsVect = recoFTS.momentum();
math::XYZVectorD reco(ftsVect.x(), ftsVect.y(), ftsVect.z());
double deltaRVal = deltaR<double>(reco.eta(),reco.phi(),
simTrack.momentum().eta(),simTrack.momentum().phi());
histo->FillDeltaR(deltaRVal);
histo->computeResolutionAndPull(recoFTS,simTrack);
}
| TrajectoryStateOnSurface MuonTrackAnalyzer::getSeedTSOS | ( | const TrajectorySeed & | seed | ) | [private] |
Definition at line 482 of file MuonTrackAnalyzer.cc.
References PTrajectoryStateOnDet::detId(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::isValid(), oppositeToMomentum, query::result, TrajectorySeed::startingState(), GeometricSearchDet::surface(), and trajectoryStateTransform::transientState().
{
// Get the Trajectory State on Det (persistent version of a TSOS) from the seed
PTrajectoryStateOnDet pTSOD = seed.startingState();
// Transform it in a TrajectoryStateOnSurface
DetId seedDetId(pTSOD.detId());
const GeomDet* gdet = theService->trackingGeometry()->idToDet( seedDetId );
TrajectoryStateOnSurface initialState = trajectoryStateTransform::transientState(pTSOD, &(gdet->surface()), &*theService->magneticField());
// Get the layer on which the seed relies
const DetLayer *initialLayer = theService->detLayerGeometry()->idToLayer( seedDetId );
PropagationDirection detLayerOrder = oppositeToMomentum;
// ask for compatible layers
vector<const DetLayer*> detLayers;
detLayers = initialLayer->compatibleLayers( *initialState.freeState(),detLayerOrder);
TrajectoryStateOnSurface result = initialState;
if(detLayers.size()){
const DetLayer* finalLayer = detLayers.back();
const TrajectoryStateOnSurface propagatedState = theService->propagator(theSeedPropagatorName)->propagate(initialState, finalLayer->surface());
if(propagatedState.isValid())
result = propagatedState;
}
return result;
}
| pair< SimTrack, double > MuonTrackAnalyzer::getSimTrack | ( | TrajectoryStateOnSurface & | tsos, |
| edm::Handle< edm::SimTrackContainer > | simTracks | ||
| ) | [private] |
Definition at line 376 of file MuonTrackAnalyzer.cc.
References abs, TrajectoryStateOnSurface::globalMomentum(), LogTrace, query::result, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
// // Loop over the Sim tracks
// SimTrackContainer::const_iterator simTrack;
// SimTrack result;
// int mu=0;
// for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
// if (abs((*simTrack).type()) == 13) {
// result = *simTrack;
// ++mu;
// }
// if(mu != 1) LogTrace("MuonTrackAnalyzer") << "WARNING!! more than 1 simulated muon!!" <<endl;
// return result;
// Loop over the Sim tracks
SimTrackContainer::const_iterator simTrack;
SimTrack result;
double bestDeltaR = 10e5;
for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack){
if (abs((*simTrack).type()) != 13) continue;
// double newDeltaR = tsos.globalMomentum().basicVector().deltaR(simTrack->momentum().vect());
GlobalVector tsosVect = tsos.globalMomentum();
math::XYZVectorD vect(tsosVect.x(), tsosVect.y(), tsosVect.z());
double newDeltaR = deltaR<double>(vect.eta(),vect.phi(),
simTrack->momentum().eta(),simTrack->momentum().phi());
if ( newDeltaR < bestDeltaR ) {
LogTrace("MuonTrackAnalyzer") << "Matching Track with DeltaR = " << newDeltaR<<endl;
bestDeltaR = newDeltaR;
result = *simTrack;
}
}
return pair<SimTrack,double>(result,bestDeltaR);
}
| bool MuonTrackAnalyzer::isInTheAcceptance | ( | double | eta | ) | [private] |
Definition at line 419 of file MuonTrackAnalyzer.cc.
References abs, cond::ecalcond::all, Reference_intrackfit_cff::barrel, Reference_intrackfit_cff::endcap, and LogTrace.
{
switch(theEtaRange){
case all:
return ( abs(eta) <= 2.4 ) ? true : false;
case barrel:
return ( abs(eta) < 1.1 ) ? true : false;
case endcap:
return ( abs(eta) >= 1.1 && abs(eta) <= 2.4 ) ? true : false;
default:
{LogTrace("MuonTrackAnalyzer")<<"No correct Eta range selected!! "<<endl; return false;}
}
}
| void MuonTrackAnalyzer::seedsAnalysis | ( | const edm::Event & | event, |
| const edm::EventSetup & | eventSetup, | ||
| edm::Handle< edm::SimTrackContainer > | simTracks | ||
| ) |
Definition at line 211 of file MuonTrackAnalyzer.cc.
References debug, and LogTrace.
{
MuonPatternRecoDumper debug;
// Get the RecTrack collection from the event
Handle<TrajectorySeedCollection> seeds;
event.getByLabel(theSeedsLabel, seeds);
LogTrace("MuonTrackAnalyzer")<<"Number of reconstructed seeds: " << seeds->size()<<endl;
for(TrajectorySeedCollection::const_iterator seed = seeds->begin();
seed != seeds->end(); ++seed){
TrajectoryStateOnSurface seedTSOS = getSeedTSOS(*seed);
pair<SimTrack,double> sim = getSimTrack(seedTSOS,simTracks);
fillPlots(seedTSOS, sim.first,
hRecoSeedInner, debug);
std::pair<bool,FreeTrajectoryState> propSeed =
theUpdator->propagateToNominalLine(seedTSOS);
if(propSeed.first)
fillPlots(propSeed.second, sim.first,
hRecoSeedPCA, debug);
else
LogTrace("MuonTrackAnalyzer")<<"Error in seed propagation"<<endl;
}
}
| void MuonTrackAnalyzer::tracksAnalysis | ( | const edm::Event & | event, |
| const edm::EventSetup & | eventSetup, | ||
| edm::Handle< edm::SimTrackContainer > | simTracks | ||
| ) |
Definition at line 241 of file MuonTrackAnalyzer.cc.
References TrajectoryStateOnSurface::charge(), debug, TrajectoryStateOnSurface::globalMomentum(), LogTrace, CoreSimTrack::momentum(), PV3DBase< T, PVType, FrameType >::perp(), mathSSE::sqrt(), lumiQTWidget::t, and testEve_cfg::tracks.
{
MuonPatternRecoDumper debug;
// Get the RecTrack collection from the event
Handle<reco::TrackCollection> tracks;
event.getByLabel(theTracksLabel, tracks);
LogTrace("MuonTrackAnalyzer")<<"Reconstructed tracks: " << tracks->size() << endl;
hNumberOfTracks->Fill(tracks->size());
if(tracks->size()) numberOfRecTracks++;
// Loop over the Rec tracks
for(reco::TrackCollection::const_iterator t = tracks->begin(); t != tracks->end(); ++t) {
reco::TransientTrack track(*t,&*theService->magneticField(),theService->trackingGeometry());
TrajectoryStateOnSurface outerTSOS = track.outermostMeasurementState();
TrajectoryStateOnSurface innerTSOS = track.innermostMeasurementState();
TrajectoryStateOnSurface pcaTSOS = track.impactPointState();
pair<SimTrack,double> sim = getSimTrack(pcaTSOS,simTracks);
SimTrack simTrack = sim.first;
hNumberOfTracksVsEta->Fill(simTrack.momentum().eta(), tracks->size());
fillPlots(track,simTrack);
LogTrace("MuonTrackAnalyzer") << "State at the outer surface: " << endl;
fillPlots(outerTSOS,simTrack,hRecoTracksOuter,debug);
LogTrace("MuonTrackAnalyzer") << "State at the inner surface: " << endl;
fillPlots(innerTSOS,simTrack,hRecoTracksInner,debug);
LogTrace("MuonTrackAnalyzer") << "State at PCA: " << endl;
fillPlots(pcaTSOS,simTrack,hRecoTracksPCA,debug);
double deltaPt_in_out = innerTSOS.globalMomentum().perp()-outerTSOS.globalMomentum().perp();
hDeltaPt_In_Out_VsEta->Fill(simTrack.momentum().eta(),deltaPt_in_out);
double deltaPt_pca_sim = pcaTSOS.globalMomentum().perp()-sqrt(simTrack.momentum().Perp2());
hDeltaPtVsEta->Fill(simTrack.momentum().eta(),deltaPt_pca_sim);
hChargeVsEta->Fill(simTrack.momentum().eta(),pcaTSOS.charge());
hChargeVsPt->Fill(sqrt(simTrack.momentum().perp2()),pcaTSOS.charge());
hPtRecVsPtGen->Fill(sqrt(simTrack.momentum().perp2()),pcaTSOS.globalMomentum().perp());
}
LogTrace("MuonTrackAnalyzer")<<"--------------------------------------------"<<endl;
}
DQMStore* MuonTrackAnalyzer::dbe_ [private] |
Definition at line 83 of file MuonTrackAnalyzer.h.
std::string MuonTrackAnalyzer::dirName_ [private] |
Definition at line 84 of file MuonTrackAnalyzer.h.
bool MuonTrackAnalyzer::doSeedsAnalysis [private] |
Definition at line 98 of file MuonTrackAnalyzer.h.
bool MuonTrackAnalyzer::doTracksAnalysis [private] |
Definition at line 97 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChargeVsEta [private] |
Definition at line 113 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChargeVsPt [private] |
Definition at line 114 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChi2 [private] |
Definition at line 105 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChi2Norm [private] |
Definition at line 106 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChi2NormVsEta [private] |
Definition at line 118 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChi2Prob [private] |
Definition at line 109 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChi2ProbVsEta [private] |
Definition at line 121 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hChi2VsEta [private] |
Definition at line 117 of file MuonTrackAnalyzer.h.
Definition at line 123 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hDeltaPtVsEta [private] |
Definition at line 122 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hDof [private] |
Definition at line 108 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hDofVsEta [private] |
Definition at line 120 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hHitsPerTrack [private] |
Definition at line 107 of file MuonTrackAnalyzer.h.
Definition at line 119 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hNumberOfTracks [private] |
Definition at line 111 of file MuonTrackAnalyzer.h.
Definition at line 112 of file MuonTrackAnalyzer.h.
MonitorElement* MuonTrackAnalyzer::hPtRecVsPtGen [private] |
Definition at line 115 of file MuonTrackAnalyzer.h.
HTrack* MuonTrackAnalyzer::hRecoSeedInner [private] |
Definition at line 127 of file MuonTrackAnalyzer.h.
HTrack* MuonTrackAnalyzer::hRecoSeedPCA [private] |
Definition at line 128 of file MuonTrackAnalyzer.h.
HTrack* MuonTrackAnalyzer::hRecoTracksInner [private] |
Definition at line 130 of file MuonTrackAnalyzer.h.
HTrack* MuonTrackAnalyzer::hRecoTracksOuter [private] |
Definition at line 131 of file MuonTrackAnalyzer.h.
HTrack* MuonTrackAnalyzer::hRecoTracksPCA [private] |
Definition at line 129 of file MuonTrackAnalyzer.h.
HTrackVariables* MuonTrackAnalyzer::hSimTracks [private] |
Definition at line 125 of file MuonTrackAnalyzer.h.
int MuonTrackAnalyzer::numberOfRecTracks [private] |
Definition at line 135 of file MuonTrackAnalyzer.h.
int MuonTrackAnalyzer::numberOfSimTracks [private] |
Definition at line 134 of file MuonTrackAnalyzer.h.
std::string MuonTrackAnalyzer::out [private] |
Definition at line 86 of file MuonTrackAnalyzer.h.
Definition at line 93 of file MuonTrackAnalyzer.h.
Definition at line 94 of file MuonTrackAnalyzer.h.
EtaRange MuonTrackAnalyzer::theEtaRange [private] |
Definition at line 89 of file MuonTrackAnalyzer.h.
Definition at line 95 of file MuonTrackAnalyzer.h.
std::string MuonTrackAnalyzer::theSeedPropagatorName [private] |
Definition at line 99 of file MuonTrackAnalyzer.h.
Definition at line 92 of file MuonTrackAnalyzer.h.
MuonServiceProxy* MuonTrackAnalyzer::theService [private] |
Definition at line 101 of file MuonTrackAnalyzer.h.
Definition at line 91 of file MuonTrackAnalyzer.h.
MuonUpdatorAtVertex* MuonTrackAnalyzer::theUpdator [private] |
Definition at line 102 of file MuonTrackAnalyzer.h.
1.7.3