| 1 |
msegala |
1.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
|
