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