| 1 |
//#ifndef TrackingAnalysis_Cosmics_HitEff_h
|
| 2 |
//#define TrackingAnalysis_Cosmics_HitEff_h
|
| 3 |
|
| 4 |
#include "FWCore/Framework/interface/ESHandle.h"
|
| 5 |
#include "FWCore/Framework/interface/EDAnalyzer.h"
|
| 6 |
#include "FWCore/Framework/interface/Event.h"
|
| 7 |
#include "DataFormats/Common/interface/Handle.h"
|
| 8 |
#include "FWCore/Framework/interface/EventSetup.h"
|
| 9 |
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
|
| 10 |
|
| 11 |
#include "FWCore/ServiceRegistry/interface/Service.h"
|
| 12 |
#include "CommonTools/UtilAlgos/interface/TFileService.h"
|
| 13 |
|
| 14 |
|
| 15 |
#include "DataFormats/Common/interface/EDProduct.h"
|
| 16 |
#include "FWCore/ParameterSet/interface/ParameterSet.h"
|
| 17 |
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
|
| 18 |
#include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
|
| 19 |
#include "TrackingTools/KalmanUpdators/interface/KFUpdator.h"
|
| 20 |
#include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
|
| 21 |
#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHitBuilder.h"
|
| 22 |
#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
|
| 23 |
#include "TrackingTools/TrackFitters/interface/KFTrajectoryFitter.h"
|
| 24 |
#include "TrackingTools/TrackFitters/interface/KFTrajectorySmoother.h"
|
| 25 |
#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
|
| 26 |
#include "MagneticField/Engine/interface/MagneticField.h"
|
| 27 |
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
|
| 28 |
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
|
| 29 |
#include "RecoTracker/SingleTrackPattern/interface/CosmicTrajectoryBuilder.h"
|
| 30 |
#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementError.h"
|
| 31 |
#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementVector.h"
|
| 32 |
#include "TrackingAnalysis/Cosmics/interface/TrackLocalAngle.h"
|
| 33 |
|
| 34 |
#include "TROOT.h"
|
| 35 |
#include "TFile.h"
|
| 36 |
#include "TH1F.h"
|
| 37 |
#include "TH2F.h"
|
| 38 |
#include <vector>
|
| 39 |
#include "TTree.h"
|
| 40 |
#include <iostream>
|
| 41 |
#include "stdlib.h"
|
| 42 |
#include "stdio.h"
|
| 43 |
#include "Riostream.h"
|
| 44 |
#include "TRandom2.h"
|
| 45 |
|
| 46 |
|
| 47 |
class HitEff : public edm::EDAnalyzer {
|
| 48 |
public:
|
| 49 |
explicit HitEff(const edm::ParameterSet& conf);
|
| 50 |
double checkConsistency(StripClusterParameterEstimator::LocalValues parameters, double xx, double xerr);
|
| 51 |
double checkConsistency(const SiStripRecHit2D* rechit, double xx, double xerr);
|
| 52 |
bool isDoubleSided(uint iidd) const;
|
| 53 |
bool check2DPartner(uint iidd, std::vector<TrajectoryMeasurement> traj);
|
| 54 |
virtual ~HitEff();
|
| 55 |
|
| 56 |
private:
|
| 57 |
//virtual void beginJob(const edm::EventSetup& c);
|
| 58 |
virtual void endJob();
|
| 59 |
virtual void analyze(const edm::Event& e, const edm::EventSetup& c);
|
| 60 |
|
| 61 |
// ----------member data ---------------------------
|
| 62 |
|
| 63 |
edm::ParameterSet conf_;
|
| 64 |
|
| 65 |
TTree* traj;
|
| 66 |
int events,EventTrackCKF;
|
| 67 |
|
| 68 |
TrackLocalAngle* theAngleFinder;
|
| 69 |
uint layers;
|
| 70 |
bool DEBUG;
|
| 71 |
uint whatlayer;
|
| 72 |
|
| 73 |
// Tree declarations
|
| 74 |
// Trajectory positions for modules included in the study
|
| 75 |
float TrajGlbX, TrajGlbY, TrajGlbZ;
|
| 76 |
float TrajLocX, TrajLocY, TrajLocErrX, TrajLocErrY, TrajLocAngleX, TrajLocAngleY;
|
| 77 |
float ClusterLocX, ClusterLocY, ClusterLocErrX, ClusterLocErrY, ClusterStoN;
|
| 78 |
float ResX, ResXSig;
|
| 79 |
uint ModIsBad; uint Id; uint SiStripQualBad; bool withinAcceptance;
|
| 80 |
uint run; uint event;
|
| 81 |
float timeDT, timeDTErr;
|
| 82 |
int timeDTDOF;
|
| 83 |
float timeECAL, dedx;
|
| 84 |
int dedxNOM;
|
| 85 |
};
|
| 86 |
|
| 87 |
|
| 88 |
//#endif
|
