#include <SiStripRecHitsValid.h>

Inheritance diagram for SiStripRecHitsValid:

Public Member Functions
	SiStripRecHitsValid (const edm::ParameterSet &conf)

	~SiStripRecHitsValid ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Protected Member Functions
virtual void	analyze (const edm::Event &e, const edm::EventSetup &c)

void	beginJob ()

void	endJob ()

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Member Functions
std::pair< LocalPoint, LocalVector >	projectHit (const PSimHit &hit, const StripGeomDetUnit *stripDet, const BoundPlane &plane)

Private Attributes
float	chi2matched [MAXHIT]

float	chi2rphi [MAXHIT]

float	chi2sas [MAXHIT]

float	cluchgrphi [MAXHIT]

float	cluchgsas [MAXHIT]

int	clusizrphi [MAXHIT]

int	clusizsas [MAXHIT]

edm::ParameterSet	conf_

DQMStore *	dbe_

std::vector< PSimHit >	matched

edm::InputTag	matchedRecHits_

MonitorElement *	meAdcRphiTEC [7]

MonitorElement *	meAdcRphiTIB [4]

MonitorElement *	meAdcRphiTID [3]

MonitorElement *	meAdcRphiTOB [6]

MonitorElement *	meAdcSasTEC [5]

MonitorElement *	meAdcSasTIB [4]

MonitorElement *	meAdcSasTID [2]

MonitorElement *	meAdcSasTOB [2]

MonitorElement *	meChi2MatchedTEC [5]

MonitorElement *	meChi2MatchedTIB [2]

MonitorElement *	meChi2MatchedTID [2]

MonitorElement *	meChi2MatchedTOB [2]

MonitorElement *	meChi2RphiTEC [7]

MonitorElement *	meChi2RphiTIB [4]

MonitorElement *	meChi2RphiTID [3]

MonitorElement *	meChi2RphiTOB [6]

MonitorElement *	meChi2SasTEC [5]

MonitorElement *	meChi2SasTIB [4]

MonitorElement *	meChi2SasTID [2]

MonitorElement *	meChi2SasTOB [2]

MonitorElement *	meErrxMatchedTEC [5]

MonitorElement *	meErrxMatchedTIB [2]

MonitorElement *	meErrxMatchedTID [2]

MonitorElement *	meErrxMatchedTOB [2]

MonitorElement *	meErrxRphiTEC [7]

MonitorElement *	meErrxRphiTIB [4]

MonitorElement *	meErrxRphiTID [3]

MonitorElement *	meErrxRphiTOB [6]

MonitorElement *	meErrxSasTEC [5]

MonitorElement *	meErrxSasTIB [4]

MonitorElement *	meErrxSasTID [2]

MonitorElement *	meErrxSasTOB [2]

MonitorElement *	meErryMatchedTEC [5]

MonitorElement *	meErryMatchedTIB [2]

MonitorElement *	meErryMatchedTID [2]

MonitorElement *	meErryMatchedTOB [2]

MonitorElement *	meNstpRphiTEC [7]

MonitorElement *	meNstpRphiTIB [4]

MonitorElement *	meNstpRphiTID [3]

MonitorElement *	meNstpRphiTOB [6]

MonitorElement *	meNstpSasTEC [5]

MonitorElement *	meNstpSasTIB [4]

MonitorElement *	meNstpSasTID [2]

MonitorElement *	meNstpSasTOB [2]

MonitorElement *	meNumMatchedTEC

MonitorElement *	meNumMatchedTIB

MonitorElement *	meNumMatchedTID

MonitorElement *	meNumMatchedTOB

MonitorElement *	meNumRphiTEC

MonitorElement *	meNumRphiTIB

MonitorElement *	meNumRphiTID

MonitorElement *	meNumRphiTOB

MonitorElement *	meNumSasTEC

MonitorElement *	meNumSasTIB

MonitorElement *	meNumSasTID

MonitorElement *	meNumSasTOB

MonitorElement *	meNumTotMatched

MonitorElement *	meNumTotRphi

MonitorElement *	meNumTotSas

MonitorElement *	mePosxMatchedTEC [5]

MonitorElement *	mePosxMatchedTIB [2]

MonitorElement *	mePosxMatchedTID [2]

MonitorElement *	mePosxMatchedTOB [2]

MonitorElement *	mePosxRphiTEC [7]

MonitorElement *	mePosxRphiTIB [4]

MonitorElement *	mePosxRphiTID [3]

MonitorElement *	mePosxRphiTOB [6]

MonitorElement *	mePosxSasTEC [5]

MonitorElement *	mePosxSasTIB [4]

MonitorElement *	mePosxSasTID [2]

MonitorElement *	mePosxSasTOB [2]

MonitorElement *	mePosyMatchedTEC [5]

MonitorElement *	mePosyMatchedTIB [2]

MonitorElement *	mePosyMatchedTID [2]

MonitorElement *	mePosyMatchedTOB [2]

MonitorElement *	mePullLFRphiTEC [7]

MonitorElement *	mePullLFRphiTIB [4]

MonitorElement *	mePullLFRphiTID [3]

MonitorElement *	mePullLFRphiTOB [6]

MonitorElement *	mePullLFSasTEC [5]

MonitorElement *	mePullLFSasTIB [4]

MonitorElement *	mePullLFSasTID [2]

MonitorElement *	mePullLFSasTOB [2]

MonitorElement *	mePullMFRphiTEC [7]

MonitorElement *	mePullMFRphiTIB [4]

MonitorElement *	mePullMFRphiTID [3]

MonitorElement *	mePullMFRphiTOB [6]

MonitorElement *	mePullMFSasTEC [5]

MonitorElement *	mePullMFSasTIB [4]

MonitorElement *	mePullMFSasTID [2]

MonitorElement *	mePullMFSasTOB [2]

MonitorElement *	meResRphiTEC [7]

MonitorElement *	meResRphiTIB [4]

MonitorElement *	meResRphiTID [3]

MonitorElement *	meResRphiTOB [6]

MonitorElement *	meResSasTEC [5]

MonitorElement *	meResSasTIB [4]

MonitorElement *	meResSasTID [2]

MonitorElement *	meResSasTOB [2]

MonitorElement *	meResxMatchedTEC [5]

MonitorElement *	meResxMatchedTIB [2]

MonitorElement *	meResxMatchedTID [2]

MonitorElement *	meResxMatchedTOB [2]

MonitorElement *	meResyMatchedTEC [5]

MonitorElement *	meResyMatchedTIB [2]

MonitorElement *	meResyMatchedTID [2]

MonitorElement *	meResyMatchedTOB [2]

std::string	outputFile_

float	rechitmatchedchi2 [MAXHIT]

float	rechitmatchederrxx [MAXHIT]

float	rechitmatchederrxy [MAXHIT]

float	rechitmatchederryy [MAXHIT]

float	rechitmatchedphi [MAXHIT]

float	rechitmatchedresx [MAXHIT]

float	rechitmatchedresy [MAXHIT]

float	rechitmatchedx [MAXHIT]

float	rechitmatchedy [MAXHIT]

float	rechitmatchedz [MAXHIT]

float	rechitrphierrx [MAXHIT]

float	rechitrphiphi [MAXHIT]

float	rechitrphipullMF [MAXHIT]

float	rechitrphires [MAXHIT]

float	rechitrphix [MAXHIT]

float	rechitrphiy [MAXHIT]

float	rechitrphiz [MAXHIT]

float	rechitsaserrx [MAXHIT]

float	rechitsasphi [MAXHIT]

float	rechitsaspullMF [MAXHIT]

float	rechitsasres [MAXHIT]

float	rechitsasx [MAXHIT]

float	rechitsasy [MAXHIT]

float	rechitsasz [MAXHIT]

edm::InputTag	rphiRecHits_

edm::InputTag	stereoRecHits_

Static Private Attributes
static const int	MAXHIT = 1000

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 51 of file SiStripRecHitsValid.h.

Constructor & Destructor Documentation

SiStripRecHitsValid::SiStripRecHitsValid ( const edm::ParameterSet & conf )

Definition at line 52 of file SiStripRecHitsValid.cc.

References dbe_, edm::ParameterSet::getUntrackedParameter(), trackerHits::histo, i, meAdcRphiTEC, meAdcRphiTIB, meAdcRphiTID, meAdcRphiTOB, meAdcSasTEC, meAdcSasTIB, meAdcSasTID, meAdcSasTOB, meChi2MatchedTEC, meChi2MatchedTIB, meChi2MatchedTID, meChi2MatchedTOB, meChi2RphiTEC, meChi2RphiTIB, meChi2RphiTID, meChi2RphiTOB, meChi2SasTEC, meChi2SasTIB, meChi2SasTID, meChi2SasTOB, meErrxMatchedTEC, meErrxMatchedTIB, meErrxMatchedTID, meErrxMatchedTOB, meErrxRphiTEC, meErrxRphiTIB, meErrxRphiTID, meErrxRphiTOB, meErrxSasTEC, meErrxSasTIB, meErrxSasTID, meErrxSasTOB, meErryMatchedTEC, meErryMatchedTIB, meErryMatchedTID, meErryMatchedTOB, meNstpRphiTEC, meNstpRphiTIB, meNstpRphiTID, meNstpRphiTOB, meNstpSasTEC, meNstpSasTIB, meNstpSasTID, meNstpSasTOB, meNumMatchedTEC, meNumMatchedTIB, meNumMatchedTID, meNumMatchedTOB, meNumRphiTEC, meNumRphiTIB, meNumRphiTID, meNumRphiTOB, meNumSasTEC, meNumSasTIB, meNumSasTID, meNumSasTOB, meNumTotMatched, meNumTotRphi, meNumTotSas, mePosxMatchedTEC, mePosxMatchedTIB, mePosxMatchedTID, mePosxMatchedTOB, mePosxRphiTEC, mePosxRphiTIB, mePosxRphiTID, mePosxRphiTOB, mePosxSasTEC, mePosxSasTIB, mePosxSasTID, mePosxSasTOB, mePosyMatchedTEC, mePosyMatchedTIB, mePosyMatchedTID, mePosyMatchedTOB, mePullLFRphiTEC, mePullLFRphiTIB, mePullLFRphiTID, mePullLFRphiTOB, mePullLFSasTEC, mePullLFSasTIB, mePullLFSasTID, mePullLFSasTOB, mePullMFRphiTEC, mePullMFRphiTIB, mePullMFRphiTID, mePullMFRphiTOB, mePullMFSasTEC, mePullMFSasTIB, mePullMFSasTID, mePullMFSasTOB, meResRphiTEC, meResRphiTIB, meResRphiTID, meResRphiTOB, meResSasTEC, meResSasTIB, meResSasTID, meResSasTOB, meResxMatchedTEC, meResxMatchedTIB, meResxMatchedTID, meResxMatchedTOB, meResyMatchedTEC, meResyMatchedTIB, meResyMatchedTID, meResyMatchedTOB, cmsCodeRules.cppFunctionSkipper::operator, and outputFile_.

                                                                :
   dbe_(0),  
   conf_(ps),
   matchedRecHits_( ps.getParameter<edm::InputTag>("matchedRecHits") ),
   rphiRecHits_( ps.getParameter<edm::InputTag>("rphiRecHits") ),
   stereoRecHits_( ps.getParameter<edm::InputTag>("stereoRecHits") ) {
 
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "sistriprechitshisto.root");
   dbe_ = Service<DQMStore>().operator->();
   //dbe_->showDirStructure();
   dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/SISTRIP");
 
   meNumTotRphi = dbe_->book1D("NumTotRphi","Num of RecHits rphi",100, 0, 10000);
   meNumTotSas = dbe_->book1D("NumTotSas","Num of RecHits sas",100, 0, 10000);
   meNumTotMatched = dbe_->book1D("NumTotMatched","Num of RecHits rmatched",100, 0, 10000);
 
   dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TIB");
   meNumRphiTIB = dbe_->book1D("NumRphiTIB","Num of RecHits rphi", 100, 0, 1000.);
   meNumSasTIB = dbe_->book1D("NumSasTIB","Num of RecHits sas", 100, 0, 1000.);
   meNumMatchedTIB = dbe_->book1D("NumMatchedTIB","Num of RecHits matched", 100, 0, 1000.);
 
   //one histo per Layer rphi hits
   for(int i = 0 ;i<4 ; i++) {
     Char_t histo[200];
     sprintf(histo,"Nstp_rphi_layer%dtib",i+1);
     meNstpRphiTIB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
     sprintf(histo,"Adc_rphi_layer%dtib",i+1);
     meAdcRphiTIB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
     sprintf(histo,"Posx_rphi_layer%dtib",i+1);
     mePosxRphiTIB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
     sprintf(histo,"Errx_rphi_layer%dtib",i+1);
     meErrxRphiTIB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.01); //<error>~20micron  
     sprintf(histo,"Res_rphi_layer%dtib",i+1);
     meResRphiTIB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);  
     sprintf(histo,"Pull_LF_rphi_layer%dtib",i+1);
     mePullLFRphiTIB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Pull_MF_rphi_layer%dtib",i+1);
     mePullMFRphiTIB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Chi2_rphi_layer%dtib",i+1);
     meChi2RphiTIB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
   }
 
   //one histo per Layer stereo and matched hits
   for(int i = 0 ;i<2 ; i++) {
     Char_t histo[200];
     sprintf(histo,"Nstp_sas_layer%dtib",i+1);
     meNstpSasTIB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
     sprintf(histo,"Adc_sas_layer%dtib",i+1);
     meAdcSasTIB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
     sprintf(histo,"Posx_sas_layer%dtib",i+1);
     mePosxSasTIB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
     sprintf(histo,"Errx_sas_layer%dtib",i+1);
     meErrxSasTIB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.01);  
     sprintf(histo,"Res_sas_layer%dtib",i+1);
     meResSasTIB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);  
     sprintf(histo,"Pull_LF_sas_layer%dtib",i+1);
     mePullLFSasTIB[i] = dbe_->book1D(histo,"Pull",100,-4.,4.);  
     sprintf(histo,"Pull_MF_sas_layer%dtib",i+1);
     mePullMFSasTIB[i] = dbe_->book1D(histo,"Pull",100,-4.,4.);  
     sprintf(histo,"Chi2_sas_layer%dtib",i+1);
     meChi2SasTIB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
 
     sprintf(histo,"Posx_matched_layer%dtib",i+1);
     mePosxMatchedTIB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);  
     sprintf(histo,"Posy_matched_layer%dtib",i+1);
     mePosyMatchedTIB[i] = dbe_->book1D(histo,"RecHit y coord.",100,-6.0, +6.0);  
     sprintf(histo,"Errx_matched_layer%dtib",i+1);
     meErrxMatchedTIB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.01);  
     sprintf(histo,"Erry_matched_layer%dtib",i+1);
     meErryMatchedTIB[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.05);  
     sprintf(histo,"Resx_matched_layer%dtib",i+1);
     meResxMatchedTIB[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.02, +0.02);  
     sprintf(histo,"Resy_matched_layer%dtib",i+1);
     meResyMatchedTIB[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);  
     sprintf(histo,"Chi2_matched_layer%dtib",i+1);
     meChi2MatchedTIB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
   }
 
   dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TOB");
   meNumRphiTOB = dbe_->book1D("NumRphiTOB","Num of RecHits rphi", 100, 0, 1000.);
   meNumSasTOB = dbe_->book1D("NumSasTOB","Num of RecHits sas", 100, 0, 1000.);
   meNumMatchedTOB = dbe_->book1D("NumMatchedTOB","Num of RecHits matched", 100, 0, 1000.);
   //one histo per Layer rphi hits
   for(int i = 0 ;i<6 ; i++) {
     Char_t histo[200];
     sprintf(histo,"Nstp_rphi_layer%dtob",i+1);
     meNstpRphiTOB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
     sprintf(histo,"Adc_rphi_layer%dtob",i+1);
     meAdcRphiTOB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
     sprintf(histo,"Posx_rphi_layer%dtob",i+1);
     mePosxRphiTOB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
     sprintf(histo,"Errx_rphi_layer%dtob",i+1);
     meErrxRphiTOB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.01);  
     sprintf(histo,"Res_rphi_layer%dtob",i+1);
     meResRphiTOB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);  
     sprintf(histo,"Pull_LF_rphi_layer%dtob",i+1);
     mePullLFRphiTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Pull_MF_rphi_layer%dtob",i+1);
     mePullMFRphiTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Chi2_rphi_layer%dtob",i+1);
     meChi2RphiTOB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
   }
 
   //one histo per Layer stereo and matched hits
   for(int i = 0 ;i<2 ; i++) {
     Char_t histo[200];
     sprintf(histo,"Nstp_sas_layer%dtob",i+1);
     meNstpSasTOB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
     sprintf(histo,"Adc_sas_layer%dtob",i+1);
     meAdcSasTOB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
     sprintf(histo,"Posx_sas_layer%dtob",i+1);
     mePosxSasTOB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
     sprintf(histo,"Errx_sas_layer%dtob",i+1);
     meErrxSasTOB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.01);  
     sprintf(histo,"Res_sas_layer%dtob",i+1);
     meResSasTOB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);  
     sprintf(histo,"Pull_LF_sas_layer%dtob",i+1);
     mePullLFSasTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Pull_MF_sas_layer%dtob",i+1);
     mePullMFSasTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Chi2_sas_layer%dtob",i+1);
     meChi2SasTOB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
 
     sprintf(histo,"Posx_matched_layer%dtob",i+1);
     mePosxMatchedTOB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);  
     sprintf(histo,"Posy_matched_layer%dtob",i+1);
     mePosyMatchedTOB[i] = dbe_->book1D(histo,"RecHit y coord.",100,-10.0, +10.0);  
     sprintf(histo,"Errx_matched_layer%dtob",i+1);
     meErrxMatchedTOB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.01);  
     sprintf(histo,"Erry_matched_layer%dtob",i+1);
     meErryMatchedTOB[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.1);  
     sprintf(histo,"Resx_matched_layer%dtob",i+1);
     meResxMatchedTOB[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.02, +0.02);  
     sprintf(histo,"Resy_matched_layer%dtob",i+1);
     meResyMatchedTOB[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);  
     sprintf(histo,"Chi2_matched_layer%dtob",i+1);
     meChi2MatchedTOB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
   }
 
   dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TID");
   meNumRphiTID = dbe_->book1D("NumRphiTID","Num of RecHits rphi", 100, 0, 1000.);
   meNumSasTID = dbe_->book1D("NumSasTID","Num of RecHits sas", 100, 0, 1000.);
   meNumMatchedTID = dbe_->book1D("NumMatchedTID","Num of RecHits matched", 100, 0, 1000.);
 
   //one histo per Ring rphi hits: 3 rings, 6 disks, 2 inner rings are glued 
   for(int i = 0 ;i<3 ; i++) {
     Char_t histo[200];
     sprintf(histo,"Nstp_rphi_layer%dtid",i+1);
     meNstpRphiTID[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
     sprintf(histo,"Adc_rphi_layer%dtid",i+1);
     meAdcRphiTID[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
     sprintf(histo,"Posx_rphi_layer%dtid",i+1);
     mePosxRphiTID[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
     sprintf(histo,"Errx_rphi_layer%dtid",i+1);
     meErrxRphiTID[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.5);  
     sprintf(histo,"Res_rphi_layer%dtid",i+1);
     meResRphiTID[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);  
     sprintf(histo,"Pull_LF_rphi_layer%dtid",i+1);
     mePullLFRphiTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Pull_MF_rphi_layer%dtid",i+1);
     mePullMFRphiTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Chi2_rphi_layer%dtid",i+1);
     meChi2RphiTID[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
   }
 
   //one histo per Ring stereo and matched hits
   for(int i = 0 ;i<2 ; i++) {
     Char_t histo[200];
     sprintf(histo,"Nstp_sas_layer%dtid",i+1);
     meNstpSasTID[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
     sprintf(histo,"Adc_sas_layer%dtid",i+1);
     meAdcSasTID[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
     sprintf(histo,"Posx_sas_layer%dtid",i+1);
     mePosxSasTID[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
     sprintf(histo,"Errx_sas_layer%dtid",i+1);
     meErrxSasTID[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.5);  
     sprintf(histo,"Res_sas_layer%dtid",i+1);
     meResSasTID[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);  
     sprintf(histo,"Pull_LF_sas_layer%dtid",i+1);
     mePullLFSasTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Pull_MF_sas_layer%dtid",i+1);
     mePullMFSasTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Chi2_sas_layer%dtid",i+1);
     meChi2SasTID[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
 
     sprintf(histo,"Posx_matched_layer%dtid",i+1);
     mePosxMatchedTID[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);  
     sprintf(histo,"Posy_matched_layer%dtid",i+1);
     mePosyMatchedTID[i] = dbe_->book1D(histo,"RecHit y coord.",100,-6.0, +6.0);  
     sprintf(histo,"Errx_matched_layer%dtid",i+1);
     meErrxMatchedTID[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.05);  
     sprintf(histo,"Erry_matched_layer%dtid",i+1);
     meErryMatchedTID[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.1);  
     sprintf(histo,"Resx_matched_layer%dtid",i+1);
     meResxMatchedTID[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.2, +0.2);  
     sprintf(histo,"Resy_matched_layer%dtid",i+1);
     meResyMatchedTID[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);  
     sprintf(histo,"Chi2_matched_layer%dtid",i+1);
     meChi2MatchedTID[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
   }
 
   dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TEC");
   meNumRphiTEC = dbe_->book1D("NumRphiTEC","Num of RecHits rphi", 100, 0, 1000.);
   meNumSasTEC = dbe_->book1D("NumSasTEC","Num of RecHits sas", 100, 0, 1000.);
   meNumMatchedTEC = dbe_->book1D("NumMatchedTEC","Num of RecHits matched", 100, 0, 1000.);
 
   //one histo per Ring rphi hits: 7 rings, 18 disks. Innermost 3 rings are same as TID above.  
   for(int i = 0 ;i<7 ; i++) {
     Char_t histo[200];
     sprintf(histo,"Nstp_rphi_layer%dtec",i+1);
     meNstpRphiTEC[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
     sprintf(histo,"Adc_rphi_layer%dtec",i+1);
     meAdcRphiTEC[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
     sprintf(histo,"Posx_rphi_layer%dtec",i+1);
     mePosxRphiTEC[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
     sprintf(histo,"Errx_rphi_layer%dtec",i+1);
     meErrxRphiTEC[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.5);  
     sprintf(histo,"Res_rphi_layer%dtec",i+1);
     meResRphiTEC[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);  
     sprintf(histo,"Pull_LF_rphi_layer%dtec",i+1);
     mePullLFRphiTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Pull_MF_rphi_layer%dtec",i+1);
     mePullMFRphiTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);  
     sprintf(histo,"Chi2_rphi_layer%dtec",i+1);
     meChi2RphiTEC[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
   }
 
   //one histo per Layer stereo and matched hits: rings 1,2,5 are double sided
   for(int i = 0 ;i<5 ; i++) {
     if(i == 0 || i == 1 || i == 4) {
       Char_t histo[200];
       sprintf(histo,"Nstp_sas_layer%dtec",i+1);
       meNstpSasTEC[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);  
       sprintf(histo,"Adc_sas_layer%dtec",i+1);
       meAdcSasTEC[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);  
       sprintf(histo,"Posx_sas_layer%dtec",i+1);
       mePosxSasTEC[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);  
       sprintf(histo,"Errx_sas_layer%dtec",i+1);
       meErrxSasTEC[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.5);  
       sprintf(histo,"Res_sas_layer%dtec",i+1);
       meResSasTEC[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);  
       sprintf(histo,"Pull_LF_sas_layer%dtec",i+1);
       mePullLFSasTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
       sprintf(histo,"Pull_MF_sas_layer%dtec",i+1);
       mePullMFSasTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
       sprintf(histo,"Chi2_sas_layer%dtec",i+1);
       meChi2SasTEC[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
       
       sprintf(histo,"Posx_matched_layer%dtec",i+1);
       mePosxMatchedTEC[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);  
       sprintf(histo,"Posy_matched_layer%dtec",i+1);
       mePosyMatchedTEC[i] = dbe_->book1D(histo,"RecHit y coord.",100,-8.0, +8.0);  
       sprintf(histo,"Errx_matched_layer%dtec",i+1);
       meErrxMatchedTEC[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.05);  
       sprintf(histo,"Erry_matched_layer%dtec",i+1);
       meErryMatchedTEC[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.1);  
       sprintf(histo,"Resx_matched_layer%dtec",i+1);
       meResxMatchedTEC[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.2, +0.2);  
       sprintf(histo,"Resy_matched_layer%dtec",i+1);
       meResyMatchedTEC[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);  
       sprintf(histo,"Chi2_matched_layer%dtec",i+1);
       meChi2MatchedTEC[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);  
     }
   }
 }

SiStripRecHitsValid::~SiStripRecHitsValid ( )

Definition at line 319 of file SiStripRecHitsValid.cc.

                                          {
  //if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
 }

Member Function Documentation

void SiStripRecHitsValid::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 332 of file SiStripRecHitsValid.cc.

References TrackerHitAssociator::associateHit(), chi2matched, chi2rphi, chi2sas, cluchgrphi, cluchgsas, clusizrphi, clusizsas, SiStripRecHit2D::cluster(), conf_, cond::rpcobgas::detid, error, edm::EventID::event(), MonitorElement::Fill(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), helper::getRange(), i, edm::EventBase::id(), TrackerGeometry::idToDet(), TrackerGeometry::idToDetUnit(), j, gen::k, prof2calltree::l, TOBDetId::layer(), TIBDetId::layer(), BaseSiTrackerRecHit2DLocalPos::localPosition(), PSimHit::localPosition(), BaseSiTrackerRecHit2DLocalPos::localPositionError(), m, matched, matchedRecHits_, MAXHIT, meAdcRphiTEC, meAdcRphiTIB, meAdcRphiTID, meAdcRphiTOB, meAdcSasTEC, meAdcSasTIB, meAdcSasTID, meAdcSasTOB, Topology::measurementError(), Topology::measurementPosition(), meChi2MatchedTEC, meChi2MatchedTIB, meChi2MatchedTID, meChi2MatchedTOB, meChi2RphiTEC, meChi2RphiTIB, meChi2RphiTID, meChi2RphiTOB, meChi2SasTEC, meChi2SasTIB, meChi2SasTID, meChi2SasTOB, meErrxMatchedTEC, meErrxMatchedTIB, meErrxMatchedTID, meErrxMatchedTOB, meErrxRphiTEC, meErrxRphiTIB, meErrxRphiTID, meErrxRphiTOB, meErrxSasTEC, meErrxSasTIB, meErrxSasTID, meErrxSasTOB, meErryMatchedTEC, meErryMatchedTIB, meErryMatchedTID, meErryMatchedTOB, meNstpRphiTEC, meNstpRphiTIB, meNstpRphiTID, meNstpRphiTOB, meNstpSasTEC, meNstpSasTIB, meNstpSasTID, meNstpSasTOB, meNumMatchedTEC, meNumMatchedTIB, meNumMatchedTID, meNumMatchedTOB, meNumRphiTEC, meNumRphiTIB, meNumRphiTID, meNumRphiTOB, meNumSasTEC, meNumSasTIB, meNumSasTID, meNumSasTOB, meNumTotMatched, meNumTotRphi, meNumTotSas, mePosxMatchedTEC, mePosxMatchedTIB, mePosxMatchedTID, mePosxMatchedTOB, mePosxRphiTEC, mePosxRphiTIB, mePosxRphiTID, mePosxRphiTOB, mePosxSasTEC, mePosxSasTIB, mePosxSasTID, mePosxSasTOB, mePosyMatchedTEC, mePosyMatchedTIB, mePosyMatchedTID, mePosyMatchedTOB, mePullLFRphiTEC, mePullLFRphiTIB, mePullLFRphiTID, mePullLFRphiTOB, mePullLFSasTEC, mePullLFSasTIB, mePullLFSasTID, mePullLFSasTOB, mePullMFRphiTEC, mePullMFRphiTIB, mePullMFRphiTID, mePullMFRphiTOB, mePullMFSasTEC, mePullMFSasTIB, mePullMFSasTID, mePullMFSasTOB, meResRphiTEC, meResRphiTIB, meResRphiTID, meResRphiTOB, meResSasTEC, meResSasTIB, meResSasTID, meResSasTOB, meResxMatchedTEC, meResxMatchedTIB, meResxMatchedTID, meResxMatchedTOB, meResyMatchedTEC, meResyMatchedTIB, meResyMatchedTID, meResyMatchedTOB, RecHit2DLocalPos::parameters(), RecHit2DLocalPos::parametersError(), position, projectHit(), csvReporter::r, dttmaxenums::R, rechitmatchederrxx, rechitmatchederrxy, rechitmatchederryy, rechitmatchedresx, rechitmatchedresy, rechitmatchedx, rechitmatchedy, rechitmatchedz, rechitrphierrx, rechitrphipullMF, rechitrphires, rechitrphix, rechitrphiy, rechitrphiz, rechitsaserrx, rechitsaspullMF, rechitsasres, rechitsasx, rechitsasy, rechitsasz, TIDDetId::ring(), TECDetId::ring(), rphiRecHits_, edm::EventID::run(), mathSSE::sqrt(), SiStripDetId::stereo(), GluedGeomDet::stereoDet(), stereoRecHits_, DetId::subdetId(), GeomDet::surface(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, StripGeomDetUnit::topology(), patCandidatesForDimuonsSequences_cff::tracker, MeasurementError::uu(), PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

                                                                             {
 
   LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();  
   //cout  << " Run = " << e.id().run() << " Event = " << e.id().event() << endl;  
   
   //--- get RecHits
   
   std::string rechitProducer = "SiStripRecHits2D";
   
   // Step A: Get Inputs 
   edm::Handle<SiStripMatchedRecHit2DCollection> rechitsmatched;
   edm::Handle<SiStripRecHit2DCollection> rechitsrphi;
   edm::Handle<SiStripRecHit2DCollection> rechitsstereo;
   e.getByLabel(matchedRecHits_, rechitsmatched);
   e.getByLabel(rphiRecHits_, rechitsrphi);
   e.getByLabel(stereoRecHits_, rechitsstereo);
 
   int numrechitrphi   =0;
   int numrechitsas    =0;
   int numrechitmatched=0;
 
   int totTibnumrechitrphi=0;
   int totTibnumrechitsas=0;
   int totTibnumrechitmatched=0;
   int totTobnumrechitrphi=0;
   int totTobnumrechitsas=0;
   int totTobnumrechitmatched=0;
   int totTidnumrechitrphi=0;
   int totTidnumrechitsas=0;
   int totTidnumrechitmatched=0;
   int totTecnumrechitrphi=0;
   int totTecnumrechitsas=0;
   int totTecnumrechitmatched=0;
   int totrechitrphi =0;
   int totrechitsas =0;
   int totrechitmatched =0;
   
   TrackerHitAssociator associate(e, conf_);
   
   edm::ESHandle<TrackerGeometry> pDD;
   es.get<TrackerDigiGeometryRecord> ().get (pDD);
   const TrackerGeometry &tracker(*pDD);
   // FIXME: this using of vector<DetId> is suboptimal, but I don't want to re-write the full class now
   std::vector<DetId> IDs; 
   IDs.reserve(rechitsrphi->size() + rechitsmatched->size() + rechitsstereo->size());
   std::transform(rechitsrphi->begin(), rechitsrphi->end(), std::back_inserter(IDs), helper::GetDetId() );
   std::transform(rechitsstereo->begin(), rechitsstereo->end(), std::back_inserter(IDs), helper::GetDetId() );
   std::transform(rechitsmatched->begin(), rechitsmatched->end(), std::back_inserter(IDs), helper::GetDetId() );
   // loop over detunits
   //  for(TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++){
   for(std::vector<DetId>::const_iterator it = IDs.begin(); it != IDs.end(); ++it ){//loop on rphi detector with hits
     uint32_t myid=((*it).rawId());       
     DetId detid = ((*it));
     
     // initialize here
     for(int i=0; i<MAXHIT; i++){
       rechitrphix[i] =0;
       rechitrphierrx[i] =0;
       rechitrphiy[i] =0;
       rechitrphiz[i] =0;
       rechitsasx[i] =0;
       rechitsaserrx[i] =0;
       rechitsasy[i] =0;
       rechitsasz[i] =0;
       clusizrphi[i] =0;
       clusizsas[i] =0;
       cluchgrphi[i] =0;
       cluchgsas[i] =0;
       rechitrphires[i]=-999.;
       rechitsasres[i]=-999.;
       rechitrphipullMF[i]=-999.;
       rechitsaspullMF[i]=-999.;
       chi2rphi[i] =0;
       chi2sas[i]=0;
       rechitmatchedx[i] =0;
       rechitmatchedy[i] =0;
       rechitmatchedz[i] =0;
       rechitmatchederrxx[i] =0;
       rechitmatchederrxy[i] =0;
       rechitmatchederryy[i] =0;
       rechitmatchedresx[i]=-999;
       rechitmatchedresy[i]=-999;
       chi2matched[i]=0;
     }
     
     numrechitrphi =0;
     //loop over rechits-rphi in the same subdetector
     std::pair<SiStripRecHit2DCollection::DetSet::const_iterator,SiStripRecHit2DCollection::DetSet::const_iterator> rechitrphiRange = helper::getRange(*rechitsrphi, detid);
     SiStripRecHit2DCollection::DetSet::const_iterator rechitrphiRangeIteratorBegin = rechitrphiRange.first;
     SiStripRecHit2DCollection::DetSet::const_iterator rechitrphiRangeIteratorEnd   = rechitrphiRange.second;
     SiStripRecHit2DCollection::DetSet::const_iterator iterrphi=rechitrphiRangeIteratorBegin;
     
     numrechitrphi = rechitrphiRangeIteratorEnd - rechitrphiRangeIteratorBegin;   
          
 
     int i=0;
     int i2=0;
 
     if(numrechitrphi > 0 ){
       totrechitrphi+=numrechitrphi;
       for(iterrphi=rechitrphiRangeIteratorBegin; iterrphi!=rechitrphiRangeIteratorEnd;++iterrphi){
     const GeomDetUnit *  det = tracker.idToDetUnit(detid);
     const StripGeomDetUnit * stripdet=(const StripGeomDetUnit*)(det);
     const StripTopology &topol=(StripTopology&)stripdet->topology();
     SiStripRecHit2D const rechit=*iterrphi;
     LocalPoint position=rechit.localPosition();
     LocalError error=rechit.localPositionError();
     MeasurementPoint Mposition;
     MeasurementError Merror;
     Mposition = topol.measurementPosition(position);
     Merror = topol.measurementError(position,error);
     SiStripRecHit2D::ClusterRef clust=rechit.cluster();
     int clusiz=0;
     int totcharge=0;
     clusiz = clust->amplitudes().size();
     const std::vector<uint8_t> amplitudes=clust->amplitudes();
       for(size_t ia=0; ia<amplitudes.size();ia++){
         totcharge+=amplitudes[ia];
       }
     rechitrphix[i] = position.x();
     rechitrphiy[i] = position.y();
     rechitrphiz[i] = position.z();
     rechitrphierrx[i] = error.xx();
     clusizrphi[i] = clusiz;
     cluchgrphi[i] = totcharge;
 
     matched.clear();
     matched = associate.associateHit(rechit);
     float mindist = 999999;
     float dist = 999999;
     PSimHit closest;
     if(!matched.empty()){
       for(vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++){
         dist = fabs(rechitrphix[i] - (*m).localPosition().x());
         if(dist<mindist){
           mindist = dist;
           closest = (*m);
         }
         rechitrphires[i] = rechitrphix[i] - closest.localPosition().x();
       }  
       rechitrphipullMF[i] = (Mposition.x() - (topol.measurementPosition(closest.localPosition())).x())/sqrt(Merror.uu());
 
       //chi2test compare rechit errors with the simhit position ( using null matrix for the simhit). 
       //Can spot problems in the geometry better than a simple residual. (thanks to BorisM)
       AlgebraicVector rhparameters = rechit.parameters();
       AlgebraicVector shparameters(2);
       shparameters[0] = closest.localPosition().x();
       shparameters[1] = closest.localPosition().y();
       AlgebraicVector r(rhparameters - shparameters);
       AlgebraicSymMatrix R  = rechit.parametersError();
       int ierr; 
       R.invert(ierr); // if (ierr != 0) throw exception;
       double est = R.similarity(r);
       //      std::cout << " ====== Chi2 test rphi hits ====== " << std::endl;
       //      std::cout << "RecHit param. = " << rhparameters << std::endl;
       //      std::cout << "RecHit errors = " << R << std::endl;
       //      std::cout << "SimHit param. = " << shparameters << std::endl;
       //      std::cout << " chi2  = " << est << std::endl;
       //      std::cout << "DEBUG BORIS,filling chi2rphi[i],i: " << i << std::endl;
       chi2rphi[i2] = est;
       i2++;
     }
     i++;
       }
     }
     
     //loop over rechits-sas in the same subdetector
     int j=0;
     int j2=0;
     numrechitsas=0;
     std::pair<SiStripRecHit2DCollection::DetSet::const_iterator,SiStripRecHit2DCollection::DetSet::const_iterator> rechitsasRange = helper::getRange(*rechitsstereo, detid);
     SiStripRecHit2DCollection::DetSet::const_iterator rechitsasRangeIteratorBegin = rechitsasRange.first;
     SiStripRecHit2DCollection::DetSet::const_iterator rechitsasRangeIteratorEnd   = rechitsasRange.second;
     SiStripRecHit2DCollection::DetSet::const_iterator itersas=rechitsasRangeIteratorBegin;
     numrechitsas = rechitsasRangeIteratorEnd - rechitsasRangeIteratorBegin;   
     if(numrechitsas > 0){
       totrechitsas+=numrechitsas;
       for(itersas=rechitsasRangeIteratorBegin; itersas!=rechitsasRangeIteratorEnd;++itersas){
     const GeomDetUnit *  det = tracker.idToDetUnit(detid);
     const StripGeomDetUnit * stripdet=(const StripGeomDetUnit*)(det);
     const StripTopology &topol=(StripTopology&)stripdet->topology();
     SiStripRecHit2D const rechit=*itersas;
     LocalPoint position=rechit.localPosition();
     LocalError error=rechit.localPositionError();
     MeasurementPoint Mposition;
     MeasurementError Merror;
     Mposition = topol.measurementPosition(position);
     Merror = topol.measurementError(position,error);
     SiStripRecHit2D::ClusterRef clust=rechit.cluster(); int clusiz=0;
     int totcharge=0;
     clusiz = clust->amplitudes().size();
     const std::vector<uint8_t> amplitudes=clust->amplitudes();
     for(size_t ia=0; ia<amplitudes.size();ia++){
       totcharge+=amplitudes[ia];
     }
     
     rechitsasx[j] = position.x();
     rechitsasy[j] = position.y();
     rechitsasz[j] = position.z();
     rechitsaserrx[j] = error.xx();
     clusizsas[j] = clusiz;
     cluchgsas[j] = totcharge;
     
     float mindist = 999999;
     float dist = 999999;
     PSimHit closest;
     matched.clear();
     matched = associate.associateHit(rechit);
     if(!matched.empty()){
       for(vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++){
         dist = fabs(rechitsasx[j] - (*m).localPosition().x());
         if(dist<mindist){
           mindist = dist;
           closest = (*m);
         }
         rechitsasres[j] = rechitsasx[j] - closest.localPosition().x();
       }  
       rechitsaspullMF[j] = (Mposition.x() - (topol.measurementPosition(closest.localPosition())).x())/sqrt(Merror.uu());
       //chi2test compare rechit errors with the simhit position ( using null matrix for the simhit). 
       //Can spot problems in the geometry better than a simple residual. (thanks to BorisM)
       AlgebraicVector rhparameters = rechit.parameters();
       AlgebraicVector shparameters(2);
       shparameters[0] = closest.localPosition().x();
       shparameters[1] = closest.localPosition().y();
       AlgebraicVector r(rhparameters - shparameters);
       AlgebraicSymMatrix R  = rechit.parametersError();
       int ierr; 
       R.invert(ierr); // if (ierr != 0) throw exception;
       double est = R.similarity(r);
       //      std::cout << " ====== Chi2 test sas hits ====== " << std::endl;
       //      std::cout << "RecHit param. = " << rhparameters << std::endl;
       //      std::cout << "RecHit errors = " << R << std::endl;
       //      std::cout << "SimHit param. = " << shparameters << std::endl;
       //      std::cout << " chi2  = " << est << std::endl;
       //      std::cout << "DEBUG BORIS,filling chi2sas[j],j: " << i << std::endl;
       chi2sas[j2] = est;
       j2++;
     }
     j++;
       }
     }
     
     //now matched hits
 
     int k=0;
     int k2=0;
     
     //loop over rechits-matched in the same subdetector
     numrechitmatched=0;
     std::pair<SiStripMatchedRecHit2DCollection::DetSet::const_iterator,SiStripMatchedRecHit2DCollection::DetSet::const_iterator> rechitmatchedRange = helper::getRange(*rechitsmatched, detid);
     SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorBegin = rechitmatchedRange.first;
     SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorEnd   = rechitmatchedRange.second;
     SiStripMatchedRecHit2DCollection::DetSet::const_iterator itermatched=rechitmatchedRangeIteratorBegin;
     numrechitmatched = rechitmatchedRangeIteratorEnd - rechitmatchedRangeIteratorBegin;   
     if(numrechitmatched > 0){
       totrechitmatched +=numrechitmatched;
 
       for(itermatched=rechitmatchedRangeIteratorBegin; itermatched!=rechitmatchedRangeIteratorEnd;++itermatched){
     SiStripMatchedRecHit2D const rechit=*itermatched;
     LocalPoint position=rechit.localPosition();
     LocalError error=rechit.localPositionError();
     
     float mindist = 999999;
     float distx = 999999;
     float disty = 999999;
     float dist = 999999;
     std::pair<LocalPoint,LocalVector> closestPair;
     matched.clear();
     //  const SiStripRecHit2D *mono = rechit.monoHit();
     //  const SiStripRecHit2D *st = rechit.stereoHit();
     //  LocalPoint monopos = mono->localPosition();
     //  LocalPoint stpos   = st->localPosition();
     
     rechitmatchedx[k] = position.x();
     rechitmatchedy[k] = position.y();
     rechitmatchedz[k] = position.z();
     rechitmatchederrxx[k] = error.xx();
     rechitmatchederrxy[k] = error.xy();
     rechitmatchederryy[k] = error.yy();
 
     //  std::cout << " before association " << std::endl;
     matched = associate.associateHit(rechit);
     //std::cout << " after association size = " << matched.size() << std::endl;
 
     if(!matched.empty()){
       //project simhit;
       const GluedGeomDet* gluedDet = (const GluedGeomDet*)tracker.idToDet(rechit.geographicalId());
       const StripGeomDetUnit* partnerstripdet =(StripGeomDetUnit*) gluedDet->stereoDet();
       std::pair<LocalPoint,LocalVector> hitPair;
 
       //std::cout << " RECHIT position = " << position << std::endl;      
       
       for(vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++){
         //project simhit;
         hitPair= projectHit((*m),partnerstripdet,gluedDet->surface());
         distx = fabs(rechitmatchedx[k] - hitPair.first.x());
         disty = fabs(rechitmatchedy[k] - hitPair.first.y());
         dist = sqrt(distx*distx+disty*disty);
         // std::cout << " Simhit position x = " << hitPair.first.x() 
         //      << " y = " << hitPair.first.y() << " dist = " << dist << std::endl;   
         if(dist<mindist){
           mindist = dist;
           closestPair = hitPair;
         }
       }
       //std::cout << " Closest position x = " << closestPair.first.x() 
       //      << " y = " << closestPair.first.y() << " dist = " << dist << std::endl;     
       rechitmatchedresx[k] = rechitmatchedx[k] - closestPair.first.x();
       rechitmatchedresy[k] = rechitmatchedy[k] - closestPair.first.y();
 
       //chi2test compare rechit errors with the simhit position ( using null matrix for the simhit). 
       //Can spot problems in the geometry better than a simple residual. (thanks to BorisM)
 
       AlgebraicVector rhparameters = rechit.parameters();
       LocalPoint sh = closestPair.first;
       AlgebraicVector shparameters(2);
       shparameters[0] = sh.x();
       shparameters[1] = sh.y();
       AlgebraicVector r(rhparameters - shparameters);
       AlgebraicSymMatrix R  = rechit.parametersError();
       int ierr; 
       R.invert(ierr); // if (ierr != 0) throw exception;
       double est = R.similarity(r);
       //      std::cout << " ====== Chi2 test matched ====== " << std::endl;
       //      std::cout << "RecHit param. = " << rhparameters << std::endl;
       //      std::cout << "RecHit errors = " << R << std::endl;
       //      std::cout << "SimHit param. = " << shparameters << std::endl;
       //      std::cout << " chi2  = " << est << std::endl;
       chi2matched[k2] = est;
       k2++;
     }
     
     k++;
       }
     }
     
     //for each detid
     if(numrechitrphi>0 || numrechitsas>0 || numrechitmatched){
       if (detid.subdetId() == int(StripSubdetector::TIB)){
     TIBDetId tibid(myid);
     int Tibnumrechitrphi    = numrechitrphi;
     int Tibnumrechitsas     = numrechitsas;
     int Tibnumrechitmatched = numrechitmatched;
 
     int Tibnumrechitrphi2    = i2;
     int Tibnumrechitsas2     = j2;
     int Tibnumrechitmatched2 = k2;
 
     totTibnumrechitrphi +=numrechitrphi;
     totTibnumrechitsas +=numrechitsas;
     totTibnumrechitmatched +=numrechitmatched;
 
     int ilay = tibid.layer() - 1; //for histogram filling
     
     if(tibid.stereo()==0){
       for(int k = 0; k<Tibnumrechitrphi; k++){
         meNstpRphiTIB[ilay]->Fill(clusizrphi[k]);
         meAdcRphiTIB[ilay]->Fill(cluchgrphi[k]);
         mePosxRphiTIB[ilay]->Fill(rechitrphix[k]);
         meErrxRphiTIB[ilay]->Fill(sqrt(rechitrphierrx[k]));
         meResRphiTIB[ilay]->Fill(rechitrphires[k]);
         mePullLFRphiTIB[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
         mePullMFRphiTIB[ilay]->Fill(rechitrphipullMF[k]);
 
       }
 
       for(int k = 0; k<Tibnumrechitrphi2; k++){
         meChi2RphiTIB[ilay]->Fill(chi2rphi[k]);
       }
 
     } else  if(tibid.stereo()==1){
       for(int kk = 0; kk < Tibnumrechitsas; kk++)       
         {
           meNstpSasTIB[ilay]->Fill(clusizsas[kk]);
           meAdcSasTIB[ilay]->Fill(cluchgsas[kk]);
           mePosxSasTIB[ilay]->Fill(rechitsasx[kk]);
           meErrxSasTIB[ilay]->Fill(sqrt(rechitsaserrx[kk]));
           meResSasTIB[ilay]->Fill(rechitsasres[kk]);
           mePullLFSasTIB[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
           mePullMFSasTIB[ilay]->Fill(rechitsaspullMF[kk]);
         }     
       for(int l = 0; l<Tibnumrechitsas2; l++){
         meChi2SasTIB[ilay]->Fill(chi2sas[l]);
       }
 
     }
     if(Tibnumrechitmatched>0){
       for(int kkk = 0; kkk<Tibnumrechitmatched; kkk++)
         {
           mePosxMatchedTIB[ilay]->Fill(rechitmatchedx[kkk]);
           mePosyMatchedTIB[ilay]->Fill(rechitmatchedy[kkk]);
           meErrxMatchedTIB[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
           meErryMatchedTIB[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
           meResxMatchedTIB[ilay]->Fill(rechitmatchedresx[kkk]);
           meResyMatchedTIB[ilay]->Fill(rechitmatchedresy[kkk]);
         }     
       for(int l = 0; l<Tibnumrechitmatched2; l++){
         meChi2MatchedTIB[ilay]->Fill(chi2matched[l]);
       }
 
     }
       }
 
 
       if (detid.subdetId() == int(StripSubdetector::TOB)){
     TOBDetId tobid(myid);
     int Tobnumrechitrphi    = numrechitrphi;
     int Tobnumrechitsas     = numrechitsas;
     int Tobnumrechitmatched = numrechitmatched;
     totTobnumrechitrphi +=numrechitrphi;
     totTobnumrechitsas +=numrechitsas;
     totTobnumrechitmatched +=numrechitmatched;
 
     int Tobnumrechitrphi2    = i2;
     int Tobnumrechitsas2     = j2;
     int Tobnumrechitmatched2 = k2;
 
     int ilay = tobid.layer() - 1; //for histogram filling
     
     if(tobid.stereo()==0){
       for(int k = 0; k<Tobnumrechitrphi; k++){
         meNstpRphiTOB[ilay]->Fill(clusizrphi[k]);
         meAdcRphiTOB[ilay]->Fill(cluchgrphi[k]);
         mePosxRphiTOB[ilay]->Fill(rechitrphix[k]);
         meErrxRphiTOB[ilay]->Fill(sqrt(rechitrphierrx[k]));
         meResRphiTOB[ilay]->Fill(rechitrphires[k]);
         mePullLFRphiTOB[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
         mePullMFRphiTOB[ilay]->Fill(rechitrphipullMF[k]);
       }
       for(int l = 0; l<Tobnumrechitrphi2; l++){
         meChi2RphiTOB[ilay]->Fill(chi2rphi[l]);
       }
     } else  if(tobid.stereo()==1){
       for(int kk = 0; kk < Tobnumrechitsas; kk++)       
         {
           meNstpSasTOB[ilay]->Fill(clusizsas[kk]);
           meAdcSasTOB[ilay]->Fill(cluchgsas[kk]);
           mePosxSasTOB[ilay]->Fill(rechitsasx[kk]);
           meErrxSasTOB[ilay]->Fill(sqrt(rechitsaserrx[kk]));
           meResSasTOB[ilay]->Fill(rechitsasres[kk]);
           mePullLFSasTOB[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
           mePullMFSasTOB[ilay]->Fill(rechitsaspullMF[kk]);
         }
       for(int l = 0; l<Tobnumrechitsas2; l++){
         meChi2SasTOB[ilay]->Fill(chi2sas[l]);
       }
       
     }
     if(Tobnumrechitmatched>0){
       for(int kkk = 0; kkk<Tobnumrechitmatched; kkk++)
         {
           mePosxMatchedTOB[ilay]->Fill(rechitmatchedx[kkk]);
           mePosyMatchedTOB[ilay]->Fill(rechitmatchedy[kkk]);
           meErrxMatchedTOB[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
           meErryMatchedTOB[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
           meResxMatchedTOB[ilay]->Fill(rechitmatchedresx[kkk]);
           meResyMatchedTOB[ilay]->Fill(rechitmatchedresy[kkk]);
         }     
       for(int l = 0; l<Tobnumrechitmatched2; l++){
         meChi2MatchedTOB[ilay]->Fill(chi2matched[l]);
       }
 
     }
       }
       if (detid.subdetId() == int(StripSubdetector::TID)){
     TIDDetId tidid(myid);
     int Tidnumrechitrphi    = numrechitrphi;
     int Tidnumrechitsas     = numrechitsas;
     int Tidnumrechitmatched = numrechitmatched;
     totTidnumrechitrphi +=numrechitrphi;
     totTidnumrechitsas +=numrechitsas;
     totTidnumrechitmatched +=numrechitmatched;
 
     int Tidnumrechitrphi2    = i2;
     int Tidnumrechitsas2     = j2;
     int Tidnumrechitmatched2 = k2;
 
     int ilay = tidid.ring() - 1; //for histogram filling
     
     if(tidid.stereo()==0){
       for(int k = 0; k<Tidnumrechitrphi; k++){
         meNstpRphiTID[ilay]->Fill(clusizrphi[k]);
         meAdcRphiTID[ilay]->Fill(cluchgrphi[k]);
         mePosxRphiTID[ilay]->Fill(rechitrphix[k]);
         meErrxRphiTID[ilay]->Fill(sqrt(rechitrphierrx[k]));
         meResRphiTID[ilay]->Fill(rechitrphires[k]);
         mePullLFRphiTID[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
         mePullMFRphiTID[ilay]->Fill(rechitrphipullMF[k]);
       }
       for(int l = 0; l<Tidnumrechitrphi2; l++){
         meChi2RphiTID[ilay]->Fill(chi2rphi[l]);
       }
 
     } else  if(tidid.stereo()==1){
       for(int kk = 0; kk < Tidnumrechitsas; kk++)       
         {
           meNstpSasTID[ilay]->Fill(clusizsas[kk]);
           meAdcSasTID[ilay]->Fill(cluchgsas[kk]);
           mePosxSasTID[ilay]->Fill(rechitsasx[kk]);
           meErrxSasTID[ilay]->Fill(sqrt(rechitsaserrx[kk]));
           meResSasTID[ilay]->Fill(rechitsasres[kk]);
           mePullLFSasTID[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
           mePullMFSasTID[ilay]->Fill(rechitsaspullMF[kk]);
         }     
       for(int l = 0; l<Tidnumrechitsas2; l++){
         meChi2SasTID[ilay]->Fill(chi2sas[l]);
       }
 
     }
     if(Tidnumrechitmatched>0){
       for(int kkk = 0; kkk<Tidnumrechitmatched; kkk++)
         {
           mePosxMatchedTID[ilay]->Fill(rechitmatchedx[kkk]);
           mePosyMatchedTID[ilay]->Fill(rechitmatchedy[kkk]);
           meErrxMatchedTID[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
           meErryMatchedTID[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
           meResxMatchedTID[ilay]->Fill(rechitmatchedresx[kkk]);
           meResyMatchedTID[ilay]->Fill(rechitmatchedresy[kkk]);
         }
       for(int l = 0; l<Tidnumrechitmatched2; l++){
         meChi2MatchedTID[ilay]->Fill(chi2matched[l]);
       }
       
     }
       }
       if (detid.subdetId() == int(StripSubdetector::TEC)){
     TECDetId tecid(myid);
     int Tecnumrechitrphi    = numrechitrphi;
     int Tecnumrechitsas     = numrechitsas;
     int Tecnumrechitmatched = numrechitmatched;
     totTecnumrechitrphi +=numrechitrphi;
     totTecnumrechitsas +=numrechitsas;
     totTecnumrechitmatched +=numrechitmatched;
 
     int Tecnumrechitrphi2    = i2;
     int Tecnumrechitsas2     = j2;
     int Tecnumrechitmatched2 = k2;
 
     int ilay = tecid.ring() - 1; //for histogram filling
     
     if(tecid.stereo()==0){
       for(int k = 0; k<Tecnumrechitrphi; k++){
         meNstpRphiTEC[ilay]->Fill(clusizrphi[k]);
         meAdcRphiTEC[ilay]->Fill(cluchgrphi[k]);
         mePosxRphiTEC[ilay]->Fill(rechitrphix[k]);
         meErrxRphiTEC[ilay]->Fill(sqrt(rechitrphierrx[k]));
         meResRphiTEC[ilay]->Fill(rechitrphires[k]);
         mePullLFRphiTEC[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
         mePullMFRphiTEC[ilay]->Fill(rechitrphipullMF[k]);
       }
       for(int l = 0; l<Tecnumrechitrphi2; l++){
         meChi2RphiTEC[ilay]->Fill(chi2rphi[l]);
       }
 
     } else  if(tecid.stereo()==1){
       for(int kk = 0; kk < Tecnumrechitsas; kk++)       
         {
           meNstpSasTEC[ilay]->Fill(clusizsas[kk]);
           meAdcSasTEC[ilay]->Fill(cluchgsas[kk]);
           mePosxSasTEC[ilay]->Fill(rechitsasx[kk]);
           meErrxSasTEC[ilay]->Fill(sqrt(rechitsaserrx[kk]));
           meResSasTEC[ilay]->Fill(rechitsasres[kk]);
           mePullLFSasTEC[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
           mePullMFSasTEC[ilay]->Fill(rechitsaspullMF[kk]);
         }     
       for(int l = 0; l<Tecnumrechitsas2; l++){
         meChi2SasTEC[ilay]->Fill(chi2sas[l]);
       }
 
     }
     if(Tecnumrechitmatched>0){
       for(int kkk = 0; kkk<Tecnumrechitmatched; kkk++)
         {
           mePosxMatchedTEC[ilay]->Fill(rechitmatchedx[kkk]);
           mePosyMatchedTEC[ilay]->Fill(rechitmatchedy[kkk]);
           meErrxMatchedTEC[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
           meErryMatchedTEC[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
           meResxMatchedTEC[ilay]->Fill(rechitmatchedresx[kkk]);
           meResyMatchedTEC[ilay]->Fill(rechitmatchedresy[kkk]);
         }     
       for(int l = 0; l<Tecnumrechitmatched2; l++){
         meChi2MatchedTEC[ilay]->Fill(chi2matched[l]);
       }
 
     }
       }
 
     }
   }
 
   //now fill the cumulative histograms of the hits
 
   meNumRphiTIB->Fill(totTibnumrechitrphi);
   meNumSasTIB->Fill(totTibnumrechitsas);
   meNumMatchedTIB->Fill(totTibnumrechitmatched);
   meNumRphiTOB->Fill(totTobnumrechitrphi);
   meNumSasTOB->Fill(totTobnumrechitsas);
   meNumMatchedTOB->Fill(totTobnumrechitmatched);
   meNumRphiTID->Fill(totTidnumrechitrphi);
   meNumSasTID->Fill(totTidnumrechitsas);
   meNumMatchedTID->Fill(totTidnumrechitmatched);
   meNumRphiTEC->Fill(totTecnumrechitrphi);
   meNumSasTEC->Fill(totTecnumrechitsas);
   meNumMatchedTEC->Fill(totTecnumrechitmatched);
 
   meNumTotRphi->Fill(totrechitrphi);
   meNumTotSas->Fill(totrechitsas);
   meNumTotMatched->Fill(totrechitmatched);
 
  
 }

void SiStripRecHitsValid::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 323 of file SiStripRecHitsValid.cc.

                                   {
 
 }

void SiStripRecHitsValid::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 327 of file SiStripRecHitsValid.cc.

References dbe_, outputFile_, and DQMStore::save().

                                  {
  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
 }

std::pair< LocalPoint, LocalVector > SiStripRecHitsValid::projectHit	(	const PSimHit &	hit,
		const StripGeomDetUnit *	stripDet,
		const BoundPlane &	plane
	)

private

Definition at line 946 of file SiStripRecHitsValid.cc.

References funct::cos(), dir, PSimHit::localDirection(), PSimHit::localPosition(), funct::sin(), StripGeomDetUnit::specificTopology(), StripTopology::strip(), StripTopology::stripAngle(), GeomDet::surface(), Surface::toGlobal(), GloballyPositioned< T >::toLocal(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

 {
   //  const StripGeomDetUnit* stripDet = dynamic_cast<const StripGeomDetUnit*>(hit.det());
   //if (stripDet == 0) throw MeasurementDetException("HitMatcher hit is not on StripGeomDetUnit");
   
   const StripTopology& topol = stripDet->specificTopology();
   GlobalPoint globalpos= stripDet->surface().toGlobal(hit.localPosition());
   LocalPoint localHit = plane.toLocal(globalpos);
   //track direction
   LocalVector locdir=hit.localDirection();
   //rotate track in new frame
   
   GlobalVector globaldir= stripDet->surface().toGlobal(locdir);
   LocalVector dir=plane.toLocal(globaldir);
   float scale = -localHit.z() / dir.z();
   
   LocalPoint projectedPos = localHit + scale*dir;
   
   //  std::cout << "projectedPos " << projectedPos << std::endl;
   
   float selfAngle = topol.stripAngle( topol.strip( hit.localPosition()));
   
   LocalVector stripDir( sin(selfAngle), cos(selfAngle), 0); // vector along strip in hit frame
   
   LocalVector localStripDir( plane.toLocal(stripDet->surface().toGlobal( stripDir)));
   
   return std::pair<LocalPoint,LocalVector>( projectedPos, localStripDir);
 }