SiStripRecHitsValid.h

Go to the documentation of this file.

#ifndef SiStripRecHitsValid_h
#define SiStripRecHitsValid_h

/* \class SiStripRecHitsValid
 *
 * Analyzer to validate RecHits in the Strip tracker
 *
 * \author Patrizia Azzi, INFN PD
 *
 * \version   1st version May 2006
 *
 ************************************************************/

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/GeometrySurface/interface/Plane.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2DCollection.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2DCollection.h"

#include <string>
#include <utility>

class DQMStore;
class MonitorElement;
class PSimHit;
class StripGeomDetUnit;

class SiStripRecHitsValid : public edm::EDAnalyzer {

 public:
  
  SiStripRecHitsValid(const edm::ParameterSet& conf);
  
  ~SiStripRecHitsValid();

 protected:

  virtual void analyze(const edm::Event& e, const edm::EventSetup& c);
  void beginJob();
  void beginRun( const edm::Run& r, const edm::EventSetup& c );
  void endJob();
  
 private: 
  //Back-End Interface
  std::pair<LocalPoint,LocalVector> projectHit( const PSimHit& hit, const StripGeomDetUnit* stripDet,
                        const BoundPlane& plane );
  std::vector<PSimHit> matched;
  std::string outputFile_;
  DQMStore* dbe_;
  MonitorElement*  meNumTotRphi;
  MonitorElement*  meNumTotSas;
  MonitorElement*  meNumTotMatched;
  MonitorElement*  meNumRphiTIB;
  MonitorElement*  meNumSasTIB;
  MonitorElement*  meNumMatchedTIB;
  MonitorElement*  meNumRphiTOB;
  MonitorElement*  meNumSasTOB;
  MonitorElement*  meNumMatchedTOB;
  MonitorElement*  meNumRphiTID;
  MonitorElement*  meNumSasTID;
  MonitorElement*  meNumMatchedTID;
  MonitorElement*  meNumRphiTEC;
  MonitorElement*  meNumSasTEC;
  MonitorElement*  meNumMatchedTEC;

  //TIB
  MonitorElement* meNstpRphiTIB[4];
  MonitorElement* meAdcRphiTIB[4];
  MonitorElement* mePosxRphiTIB[4];
  MonitorElement* meErrxRphiTIB[4];
  MonitorElement* meResRphiTIB[4];
  MonitorElement* mePullLFRphiTIB[4];
  MonitorElement* mePullMFRphiTIB[4];
  MonitorElement* meChi2RphiTIB[4];
  MonitorElement* meNstpSasTIB[4];
  MonitorElement* meAdcSasTIB[4];
  MonitorElement* mePosxSasTIB[4];
  MonitorElement* meErrxSasTIB[4];
  MonitorElement* meResSasTIB[4];
  MonitorElement* mePullLFSasTIB[4];
  MonitorElement* mePullMFSasTIB[4];
  MonitorElement* meChi2SasTIB[4];
  MonitorElement* mePosxMatchedTIB[2];
  MonitorElement* mePosyMatchedTIB[2];
  MonitorElement* meErrxMatchedTIB[2];
  MonitorElement* meErryMatchedTIB[2];
  MonitorElement* meResxMatchedTIB[2];
  MonitorElement* meResyMatchedTIB[2];
  MonitorElement* meChi2MatchedTIB[2];
  //TOB
  MonitorElement* meNstpRphiTOB[6];
  MonitorElement* meAdcRphiTOB[6];
  MonitorElement* mePosxRphiTOB[6];
  MonitorElement* meErrxRphiTOB[6];
  MonitorElement* meResRphiTOB[6];
  MonitorElement* mePullLFRphiTOB[6];
  MonitorElement* mePullMFRphiTOB[6];
  MonitorElement* meChi2RphiTOB[6];
  MonitorElement* meNstpSasTOB[2];
  MonitorElement* meAdcSasTOB[2];
  MonitorElement* mePosxSasTOB[2];
  MonitorElement* meErrxSasTOB[2];
  MonitorElement* meResSasTOB[2];
  MonitorElement* mePullLFSasTOB[2];
  MonitorElement* mePullMFSasTOB[2];
  MonitorElement* meChi2SasTOB[2];
  MonitorElement* mePosxMatchedTOB[2];
  MonitorElement* mePosyMatchedTOB[2];
  MonitorElement* meErrxMatchedTOB[2];
  MonitorElement* meErryMatchedTOB[2];
  MonitorElement* meResxMatchedTOB[2];
  MonitorElement* meResyMatchedTOB[2];
  MonitorElement* meChi2MatchedTOB[2];
  //TID
  MonitorElement* meNstpRphiTID[3];
  MonitorElement* meAdcRphiTID[3];
  MonitorElement* mePosxRphiTID[3];
  MonitorElement* meErrxRphiTID[3];
  MonitorElement* meResRphiTID[3];
  MonitorElement* mePullLFRphiTID[3];
  MonitorElement* mePullMFRphiTID[3];
  MonitorElement* meChi2RphiTID[3];
  MonitorElement* meNstpSasTID[2];
  MonitorElement* meAdcSasTID[2];
  MonitorElement* mePosxSasTID[2];
  MonitorElement* meErrxSasTID[2];
  MonitorElement* meResSasTID[2];
  MonitorElement* mePullLFSasTID[2];
  MonitorElement* mePullMFSasTID[2];
  MonitorElement* meChi2SasTID[2];

  MonitorElement* mePosxMatchedTID[2];
  MonitorElement* mePosyMatchedTID[2];
  MonitorElement* meErrxMatchedTID[2];
  MonitorElement* meErryMatchedTID[2];
  MonitorElement* meResxMatchedTID[2];
  MonitorElement* meResyMatchedTID[2];
  MonitorElement* meChi2MatchedTID[2];
  //TEC
  MonitorElement* meNstpRphiTEC[7];
  MonitorElement* meAdcRphiTEC[7];
  MonitorElement* mePosxRphiTEC[7];
  MonitorElement* meErrxRphiTEC[7];
  MonitorElement* meResRphiTEC[7];
  MonitorElement* mePullLFRphiTEC[7];
  MonitorElement* mePullMFRphiTEC[7];
  MonitorElement* meChi2RphiTEC[7];
  MonitorElement* meNstpSasTEC[5];
  MonitorElement* meAdcSasTEC[5];
  MonitorElement* mePosxSasTEC[5];
  MonitorElement* meErrxSasTEC[5];
  MonitorElement* meResSasTEC[5];
  MonitorElement* mePullLFSasTEC[5];
  MonitorElement* mePullMFSasTEC[5];
  MonitorElement* meChi2SasTEC[5];
  MonitorElement* mePosxMatchedTEC[5];
  MonitorElement* mePosyMatchedTEC[5];
  MonitorElement* meErrxMatchedTEC[5];
  MonitorElement* meErryMatchedTEC[5];
  MonitorElement* meResxMatchedTEC[5];
  MonitorElement* meResyMatchedTEC[5];
  MonitorElement* meChi2MatchedTEC[5];

  static constexpr int MAXHIT = 1000;
  float rechitrphix[MAXHIT];
  float rechitrphierrx[MAXHIT];
  float rechitrphiy[MAXHIT];
  float rechitrphiz[MAXHIT];
  float rechitrphiphi[MAXHIT];
  float rechitrphires[MAXHIT];
  float rechitrphipullMF[MAXHIT];
  int   clusizrphi[MAXHIT];
  float cluchgrphi[MAXHIT];
  float rechitsasx[MAXHIT];
  float rechitsaserrx[MAXHIT];
  float rechitsasy[MAXHIT];
  float rechitsasz[MAXHIT];
  float rechitsasphi[MAXHIT];
  float rechitsasres[MAXHIT];
  float rechitsaspullMF[MAXHIT];
  int   clusizsas[MAXHIT];
  float cluchgsas[MAXHIT];
  float chi2rphi[MAXHIT];
  float chi2sas[MAXHIT];
  float chi2matched[MAXHIT];
  float rechitmatchedx[MAXHIT];
  float rechitmatchedy[MAXHIT];
  float rechitmatchedz[MAXHIT];
  float rechitmatchederrxx[MAXHIT];
  float rechitmatchederrxy[MAXHIT];
  float rechitmatchederryy[MAXHIT];
  float rechitmatchedphi[MAXHIT];
  float rechitmatchedresx[MAXHIT];
  float rechitmatchedresy[MAXHIT];
  float rechitmatchedchi2[MAXHIT];

  edm::ParameterSet conf_;
  //const StripTopology* topol;

  edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> siStripMatchedRecHit2DCollectionToken_;
  edm::EDGetTokenT<SiStripRecHit2DCollection> siStripRecHit2DCollection_rphi_Token_, siStripRecHit2DCollection_stereo_Token_;
};

#endif