Go to the documentation of this file.00001 #include "FWCore/Framework/interface/ESHandle.h"
00002 #include "FWCore/Framework/interface/EDAnalyzer.h"
00003 #include "FWCore/Framework/interface/Event.h"
00004 #include "DataFormats/Common/interface/Handle.h"
00005 #include "FWCore/Framework/interface/EventSetup.h"
00006 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
00007
00008 #include "FWCore/ServiceRegistry/interface/Service.h"
00009 #include "CommonTools/UtilAlgos/interface/TFileService.h"
00010
00011 #include "FWCore/ParameterSet/interface/ParameterSet.h"
00012 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
00013 #include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
00014 #include "TrackingTools/KalmanUpdators/interface/KFUpdator.h"
00015 #include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
00016 #include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHitBuilder.h"
00017 #include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
00018 #include "TrackingTools/TrackFitters/interface/KFTrajectoryFitter.h"
00019 #include "TrackingTools/TrackFitters/interface/KFTrajectorySmoother.h"
00020 #include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
00021 #include "MagneticField/Engine/interface/MagneticField.h"
00022 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
00023 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
00024 #include "RecoTracker/SingleTrackPattern/interface/CosmicTrajectoryBuilder.h"
00025 #include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementError.h"
00026 #include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementVector.h"
00027 #include "RecoLocalTracker/ClusterParameterEstimator/interface/StripClusterParameterEstimator.h"
00028
00029
00030 #include "TROOT.h"
00031 #include "TFile.h"
00032 #include "TH1F.h"
00033 #include "TH2F.h"
00034 #include <vector>
00035 #include "TTree.h"
00036 #include <iostream>
00037 #include <cstdlib>
00038 #include <cstdio>
00039 #include "Riostream.h"
00040 #include "TRandom2.h"
00041
00042
00043 class HitEff : public edm::EDAnalyzer {
00044 public:
00045 explicit HitEff(const edm::ParameterSet& conf);
00046 double checkConsistency(StripClusterParameterEstimator::LocalValues parameters, double xx, double xerr);
00047 bool isDoubleSided(unsigned int iidd) const;
00048 bool check2DPartner(unsigned int iidd, std::vector<TrajectoryMeasurement> traj);
00049 virtual ~HitEff();
00050 unsigned int checkLayer(unsigned int iidd);
00051
00052 private:
00053 virtual void beginJob();
00054 virtual void endJob();
00055 virtual void analyze(const edm::Event& e, const edm::EventSetup& c);
00056
00057
00058
00059 edm::ParameterSet conf_;
00060
00061 TTree* traj;
00062 int events,EventTrackCKF;
00063
00064 unsigned int layers;
00065 bool DEBUG;
00066 unsigned int whatlayer;
00067
00068
00069
00070 float TrajGlbX, TrajGlbY, TrajGlbZ;
00071 float TrajLocX, TrajLocY, TrajLocErrX, TrajLocErrY, TrajLocAngleX, TrajLocAngleY;
00072 float ClusterLocX, ClusterLocY, ClusterLocErrX, ClusterLocErrY, ClusterStoN;
00073 float ResX, ResXSig;
00074 unsigned int ModIsBad; unsigned int Id; unsigned int SiStripQualBad; bool withinAcceptance;
00075 int nHits, nLostHits;
00076 float p, pT, chi2;
00077 unsigned int trajHitValid, run, event, bunchx;
00078 float timeDT, timeDTErr;
00079 int timeDTDOF;
00080 float timeECAL, dedx;
00081 int dedxNOM;
00082 int tquality;
00083 int istep;
00084 };
00085
00086
00087