SiStripElectronAnalyzer Class Reference

Inheritance diagram for SiStripElectronAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	endJob (void) override
virtual void	initNtuple (void)
	SiStripElectronAnalyzer (const edm::ParameterSet &)
	~SiStripElectronAnalyzer () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
double	unwrapPhi (double phi) const

Private Attributes
std::string	basicClusterCollection_
std::string	basicClusterProducer_
std::string	eBRecHitCollection_
std::string	eBRecHitProducer_
std::string	electronCollection_
std::string	electronProducer_
TH1F *	emaxSuperClusters_
TH1F *	emaxSuperClustersEl_
TH1F *	energySuperClusters_
TH1F *	energySuperClustersEl_
TH1F *	energySuperClustersFailed
TH1F *	energySuperClustersPassed
TH1F *	eOverPFailed
TH1F *	eOverPPassed
float	EShower_ [1000]
TFile *	file_
std::string	fileName_
int	MatchedDetector_ [1000]
float	MatchedHitCorr_ [1000]
float	MatchedHitNoise_ [1000]
float	MatchedHitPhi_ [1000]
float	MatchedHitR_ [1000]
float	MatchedHitSignal_ [1000]
float	MatchedHitSigX_ [1000]
float	MatchedHitSigY_ [1000]
float	MatchedHitTheta_ [1000]
int	MatchedHitWidth_ [1000]
float	MatchedHitX_ [1000]
float	MatchedHitY_ [1000]
float	MatchedHitZ_ [1000]
int	MatchedLayer_ [1000]
std::string	mctruthCollection_
std::string	mctruthProducer_
int	MonoDetector_ [1000]
float	MonoHitCorr_ [1000]
float	MonoHitNoise_ [1000]
float	MonoHitPhi_ [1000]
float	MonoHitR_ [1000]
float	MonoHitSignal_ [1000]
float	MonoHitSigX_ [1000]
float	MonoHitSigY_ [1000]
float	MonoHitTheta_ [1000]
int	MonoHitWidth_ [1000]
float	MonoHitX_ [1000]
float	MonoHitY_ [1000]
float	MonoHitZ_ [1000]
int	MonoLayer_ [1000]
TTree *	myTree_
int	NMatchedHits_
int	NMonoHits_
int	NShowers_
int	NStereoHits_
TH1F *	numCand_
TH1F *	numElectrons_
TH1F *	numSiMatchedHits_
TH1F *	numSiMonoHits_
TH1F *	numSiStereoHits_
TH1F *	numSuperClusters_
TH1F *	pDiff
TH1F *	pElectronFailed
TH1F *	pElectronPassed
TH1F *	phiWidthSuperClusters_
TH1F *	phiWidthSuperClustersEl_
TH1F *	ptDiff
TH1F *	ptElectronFailed
TH1F *	ptElectronPassed
std::string	siElectronCollection_
std::string	siElectronProducer_
std::string	siHitProducer_
std::string	siMatchedHitCollection_
std::string	siRphiHitCollection_
std::string	siStereoHitCollection_
TH1F *	sizeSuperClusters_
TH1F *	sizeSuperClustersEl_
TH1F *	sizeSuperClustersFailed
TH1F *	sizeSuperClustersPassed
int	StereoDetector_ [1000]
float	StereoHitCorr_ [1000]
float	StereoHitNoise_ [1000]
float	StereoHitPhi_ [1000]
float	StereoHitR_ [1000]
float	StereoHitSignal_ [1000]
float	StereoHitSigX_ [1000]
float	StereoHitSigY_ [1000]
float	StereoHitTheta_ [1000]
int	StereoHitWidth_ [1000]
float	StereoHitX_ [1000]
float	StereoHitY_ [1000]
float	StereoHitZ_ [1000]
int	StereoLayer_ [1000]
std::string	superClusterCollection_
std::string	superClusterProducer_
const edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	topoToken_
const edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerToken_
float	XShower_ [1000]
float	YShower_ [1000]
float	ZShower_ [1000]

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 57 of file SiStripElectronAnalyzer.cc.

Constructor & Destructor Documentation

SiStripElectronAnalyzer::SiStripElectronAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 209 of file SiStripElectronAnalyzer.cc.

References eBRecHitCollection_, eBRecHitProducer_, electronCollection_, electronProducer_, emaxSuperClusters_, emaxSuperClustersEl_, energySuperClusters_, energySuperClustersEl_, energySuperClustersFailed, energySuperClustersPassed, eOverPFailed, eOverPPassed, file_, fileName_, edm::ParameterSet::getParameter(), mctruthCollection_, mctruthProducer_, numCand_, numElectrons_, numSiMatchedHits_, numSiMonoHits_, numSiStereoHits_, numSuperClusters_, pDiff, pElectronFailed, pElectronPassed, phiWidthSuperClusters_, phiWidthSuperClustersEl_, ptDiff, ptElectronFailed, ptElectronPassed, siElectronCollection_, siElectronProducer_, siHitProducer_, siMatchedHitCollection_, siRphiHitCollection_, siStereoHitCollection_, sizeSuperClusters_, sizeSuperClustersEl_, sizeSuperClustersFailed, sizeSuperClustersPassed, AlCaHLTBitMon_QueryRunRegistry::string, superClusterCollection_, and superClusterProducer_.

    : topoToken_(esConsumes()), trackerToken_(esConsumes()) {
  //now do what ever initialization is needed
  fileName_ = iConfig.getParameter<std::string>("fileName");

  file_ = new TFile(fileName_.c_str(), "RECREATE");
  numCand_ = new TH1F("numCandidates", "Number of candidates found", 10, -0.5, 9.5);
  numElectrons_ = new TH1F("numElectrons", "Number of Electrons found", 10, -0.5, 9.5);
  numSuperClusters_ = new TH1F("numSuperClusters", "Number of Ecal SuperClusters", 50, 0, 50);

  energySuperClusters_ = new TH1F("energySuperClusters", "Super Cluster Energy - all ", 200, 0, 2000.);
  energySuperClustersEl_ = new TH1F("energySuperClustersEl", "Super Cluster Energy - Electron Cands ", 200, 0., 2000.);

  sizeSuperClusters_ = new TH1F("sizeSuperClusters", "Super Cluster Size - all ", 20, 0, 19);
  sizeSuperClustersEl_ = new TH1F("sizeSuperClustersEl", "Super Cluster Size - Electron Cands ", 20, 0, 19);

  emaxSuperClusters_ = new TH1F("emaxSuperClusters", "Super Cluster Emax - all ", 200, 0, 2000.);
  emaxSuperClustersEl_ = new TH1F("emaxSuperClustersEl", "Super Cluster Emax - Electron Cands ", 200, 0, 2000.);

  phiWidthSuperClusters_ = new TH1F("phiWidthSuperClusters", "Super Cluster Width - all ", 20, 0., 0.05);
  phiWidthSuperClustersEl_ = new TH1F("phiWidthSuperClustersEl", "Super Cluster Width - Electron Cands ", 20, 0., 0.05);

  ptDiff = new TH1F("ptDiff", " ptDiff ", 20, -10., 10.);
  pDiff = new TH1F("pDiff", " pDiff ", 100, -50., 50.);

  pElectronFailed = new TH1F("pElectronFailed", " pElectronFailed ", 55, 0., 110.);
  ptElectronFailed = new TH1F("ptElectronFailed", " ptElectronFailed ", 55, 0., 110.);

  pElectronPassed = new TH1F("pElectronPassed", " pElectronPassed ", 55, 0., 110.);
  ptElectronPassed = new TH1F("ptElectronPassed", " ptElectronPassed ", 55, 0., 110.);

  sizeSuperClustersFailed = new TH1F("sizeSuperClustersFailed", "Super Cluster Size - Failed ", 20, 0, 19);
  sizeSuperClustersPassed = new TH1F("sizeSuperClustersPassed", "Super Cluster Size - Passed ", 20, 0, 19);

  energySuperClustersPassed = new TH1F("energySuperClustersPassed", "Super Cluster Energy - Passed ", 125, 0, 250.);
  energySuperClustersFailed = new TH1F("energySuperClustersFailed", "Super Cluster Energy - Failed ", 125, 0, 250.);

  eOverPFailed = new TH1F("eOverPFailed", " E over P - Failed ", 50, 0, 10.);
  eOverPPassed = new TH1F("eOverPPassed", " E over P - Passed ", 50, 0, 10.);

  numSiStereoHits_ = new TH1F("numSiStereoHits", "Number of Si StereoHits", 100, 0, 1000);
  numSiMonoHits_ = new TH1F("numSiMonoHits", "Number of Si MonoHits", 100, 0, 1000);
  numSiMatchedHits_ = new TH1F("numSiMatchedHits", "Number of Si MatchedHits", 100, 0, 1000);


  mctruthProducer_ = iConfig.getParameter<std::string>("mctruthProducer");
  mctruthCollection_ = iConfig.getParameter<std::string>("mctruthCollection");

  superClusterProducer_ = iConfig.getParameter<std::string>("superClusterProducer");
  superClusterCollection_ = iConfig.getParameter<std::string>("superClusterCollection");

  eBRecHitProducer_ = iConfig.getParameter<std::string>("recHitProducer");
  eBRecHitCollection_ = iConfig.getParameter<std::string>("recHitCollection");

  siElectronProducer_ = iConfig.getParameter<std::string>("siElectronProducer");
  siElectronCollection_ = iConfig.getParameter<std::string>("siElectronCollection");

  electronProducer_ = iConfig.getParameter<std::string>("electronProducer");
  electronCollection_ = iConfig.getParameter<std::string>("electronCollection");

  siHitProducer_ = iConfig.getParameter<std::string>("siHitProducer");
  siRphiHitCollection_ = iConfig.getParameter<std::string>("siRphiHitCollection");
  siStereoHitCollection_ = iConfig.getParameter<std::string>("siStereoHitCollection");
  siMatchedHitCollection_ = iConfig.getParameter<std::string>("siMatchedHitCollection");
}

SiStripElectronAnalyzer::~SiStripElectronAnalyzer ( )

override

Definition at line 281 of file SiStripElectronAnalyzer.cc.

References file_.

                                                  {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)

  file_->Write();
  file_->Close();
}

Member Function Documentation

void SiStripElectronAnalyzer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overridevirtual

For events w/ more than 1 electron candidate, try to plot m(e,e)

//////////////////////////////////////////////////////////////////////////////// ////////////////// Now for tracker hits: ///////////////////////////////////////

loop again to get all info into myTree

Implements edm::one::EDAnalyzerBase.

Definition at line 439 of file SiStripElectronAnalyzer.cc.

References funct::abs(), SiStripRecHit2D::cluster(), counter, alignCSCRings::e, eBRecHitCollection_, eBRecHitProducer_, electronCollection_, electronProducer_, HI_PhotonSkim_cff::electrons, emaxSuperClusters_, emaxSuperClustersEl_, relval_parameters_module::energy, energySuperClusters_, energySuperClustersEl_, EShower_, edm::Event::getByLabel(), edm::EventSetup::getHandle(), mps_fire::i, gpuClustering::id, initNtuple(), edm::Ref< C, T, F >::isNonnull(), BaseTrackerRecHit::localPositionError(), LogDebug, mag(), MatchedDetector_, MatchedHitCorr_, MatchedHitNoise_, MatchedHitPhi_, MatchedHitR_, MatchedHitSignal_, MatchedHitSigX_, MatchedHitSigY_, MatchedHitTheta_, MatchedHitWidth_, MatchedHitX_, MatchedHitY_, MatchedHitZ_, MatchedLayer_, MonoDetector_, MonoHitCorr_, MonoHitNoise_, MonoHitPhi_, MonoHitR_, MonoHitSignal_, MonoHitSigX_, MonoHitSigY_, MonoHitTheta_, MonoHitWidth_, MonoHitX_, MonoHitY_, MonoHitZ_, MonoLayer_, myMaxHits, myTree_, NMatchedHits_, NMonoHits_, NShowers_, NStereoHits_, numCand_, numElectrons_, numSiMatchedHits_, numSiMonoHits_, numSiStereoHits_, numSuperClusters_, pDiff, pElectronFailed, pElectronPassed, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), phiWidthSuperClusters_, phiWidthSuperClustersEl_, position, funct::pow(), ptDiff, ptElectronFailed, ptElectronPassed, siElectronCollection_, siElectronProducer_, siHitProducer_, siMatchedHitCollection_, siRphiHitCollection_, siStereoHitCollection_, sizeSuperClusters_, sizeSuperClustersEl_, sizeSuperClustersFailed, sizeSuperClustersPassed, edmScanValgrind::skipThis, mathSSE::sqrt(), StereoDetector_, StereoHitCorr_, StereoHitNoise_, StereoHitPhi_, StereoHitR_, StereoHitSignal_, StereoHitSigX_, StereoHitSigY_, StereoHitTheta_, StereoHitWidth_, StereoHitX_, StereoHitY_, StereoHitZ_, StereoLayer_, str, superClusterCollection_, superClusterProducer_, PV3DBase< T, PVType, FrameType >::theta(), StripSubdetector::TIB, StripSubdetector::TOB, topoToken_, trackerToken_, PV3DBase< T, PVType, FrameType >::x(), XShower_, LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), YShower_, LocalError::yy(), PV3DBase< T, PVType, FrameType >::z(), and ZShower_.

                                                                                         {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopo = iSetup.getHandle(topoToken_);

  using namespace std;  // so you can say "cout" and "endl"

  initNtuple();

  // http://cmsdoc.cern.ch/swdev/lxr/CMSSW/source/clhep/CLHEP/HepMC/GenParticle.h
  // http://cmsdoc.cern.ch/swdev/lxr/CMSSW/source/clhep/CLHEP/HepMC/GenVertex.h
  // removed by JED - causes trouble in release post 0_9_0
  //   edm::Handle<edm::HepMCProduct> mctruthHandle;
  //   iEvent.getByLabel(mctruthProducer_, mctruthCollection_, mctruthHandle);
  //   HepMC::GenEvent mctruth = mctruthHandle->getHepMCData();

  //   for (HepMC::GenEvent::particle_const_iterator partIter = mctruth.particles_begin();
  //    partIter != mctruth.particles_end();
  //    ++partIter) {
  // //    for (HepMC::GenEvent::vertex_const_iterator vertIter = mctruth.vertices_begin();
  // //     vertIter != mctruth.vertices_end();
  // //     ++vertIter) {
  //    CLHEP::HepLorentzVector creation = (*partIter)->CreationVertex();
  //    CLHEP::HepLorentzVector momentum = (*partIter)->Momentum();
  //   HepPDT::ParticleID id = (*partIter)->particleID();  // electrons and positrons are 11 and -11
  //     edm::LogInfo("") << "MC particle id " << id.pid() << ", creationVertex " << creation << " cm, initialMomentum " << momentum << " GeV/c" << endl;
  //  }

  // load the rechits for the Ecal
  edm::Handle<EcalRecHitCollection> pRecHits;
  iEvent.getByLabel(eBRecHitProducer_, eBRecHitCollection_, pRecHits);
  // Create a pointer to the RecHits  - unused for now
  //  const EcalRecHitCollection *hitCollection = pRecHits.product();

  // http://cmsdoc.cern.ch/swdev/lxr/CMSSW/source/self/DataFormats/EgammaReco/interface/SuperCluster.h
  edm::Handle<reco::SuperClusterCollection> clusterHandle;
  iEvent.getByLabel(superClusterProducer_, superClusterCollection_, clusterHandle);


  LogDebug("") << " Start loop over " << clusterHandle->end() - clusterHandle->begin() << "  superClusters ";

  for (reco::SuperClusterCollection::const_iterator clusterIter = clusterHandle->begin();
       clusterIter != clusterHandle->end();
       ++clusterIter) {
    double energy = clusterIter->energy();
    math::XYZPoint position = clusterIter->position();
    std::ostringstream str;

    str << " SuperCluster " << energy << " GeV, position " << position << " cm"
        << "\n";

    energySuperClusters_->Fill(energy);
    sizeSuperClusters_->Fill(clusterIter->clustersSize());
    // this only makes sense for hybrid superclusters

    // try to point to the constituent clusters for this SuperCluster

    str << "About to loop over basicClusters"
        << "\n";

    double emaxSuperCluster = 0.;
    double phibar = 0.;
    double phi2bar = 0.;
    double eTotSuperCluster = 0.;

    for (reco::CaloCluster_iterator basicClusterIter = clusterIter->clustersBegin();
         basicClusterIter != clusterIter->clustersEnd();
         ++basicClusterIter) {
      //std::vector<DetId> theIds= (*basicClusterIter)->getHitsByDetId();

      str << " basicCluster Energy " << (*basicClusterIter)->energy() << " Position " << (*basicClusterIter)->position()
          << " \n"
          << "        Position phi " << (*basicClusterIter)->position().phi() << " recHits "
          << (*basicClusterIter)->size() << " \n";

      double eCluster = (*basicClusterIter)->energy();
      if (eCluster > emaxSuperCluster) {
        emaxSuperCluster = eCluster;
      }
      eTotSuperCluster += eCluster;
      double phiCluster = (*basicClusterIter)->position().phi();
      phibar += eCluster * phiCluster;
      phi2bar += eCluster * phiCluster * phiCluster;

    }  // end of basicClusterIter loop

    phibar = phibar / eTotSuperCluster;
    phi2bar = phi2bar / eTotSuperCluster;
    double phiWidth = phi2bar - phibar * phibar;
    if (phiWidth > 0.) {
      phiWidth = std::pow(phiWidth, 0.5);
    } else {
      phiWidth = 0.;
    }
    str << " SuperCluster stats "
        << "\n";
    str << "phibar " << phibar << " phi2bar " << phi2bar << " eTotSuperCluster " << eTotSuperCluster << " phiWidth "
        << phiWidth << std::endl;

    phiWidthSuperClusters_->Fill(phiWidth);

    emaxSuperClusters_->Fill(emaxSuperCluster);

    str << " Done with this SuperCluster " << std::endl;

    LogDebug("") << str.str();

  }  // end of loop over superClusters

  LogDebug("") << " End loop over superClusters ";


  //
  // loop over all EcalRecHits and print out their x,y,z,E
  //  edm::LogInfo("") << " Dumping all recHits in this event " << endl ;
  // for(EcalRecHitCollection::const_iterator _blah = hitCollection->begin();
  //     _blah != hitCollection->end() ; ++_blah ) {
  //   edm::LogInfo("") << "Ecal RecHit Energy: " << _blah->energy() << endl ;
  //   //      " Position " << _blah.position() << endl ;
  //  }
  //
  //  edm::LogInfo("") << "Dump finished " << endl ;
  //

  // DataFormats/EgammaCandidates/src/SiStripElectron.cc
  edm::Handle<reco::SiStripElectronCollection> siStripElectronHandle;
  iEvent.getByLabel(siElectronProducer_, siElectronCollection_, siStripElectronHandle);


  LogDebug("") << " Dumping Algo's guess of SiStripElectron Candidate Info ";
  int numberOfElectrons = 0;
  // need to check if fit succeeded
  LogDebug("") << " Number of SiStripElectrons  " << siStripElectronHandle->size();

  for (reco::SiStripElectronCollection::const_iterator electronIter = siStripElectronHandle->begin();
       electronIter != siStripElectronHandle->end();
       ++electronIter) {
    LogDebug("") << "about to get stuff from electroncandidate " << numberOfElectrons << "\n"
                 << "supercluster energy = " << electronIter->superCluster()->energy() << "\n"
                 << "fit results are phi(r) = " << electronIter->phiAtOrigin() << " + " << electronIter->phiVsRSlope()
                 << "*r"
                 << "\n"
                 << " chi2 " << electronIter->chi2() << " ndof " << electronIter->ndof() << "\n"
                 << " Pt " << electronIter->pt() << "\n"
                 << "P, Px, Py, Pz  " << electronIter->p() << " " << electronIter->px() << " " << electronIter->py()
                 << " " << electronIter->pz() << "\n"
                 << "you get the idea...";

    // make plots for supercluster that an electron has been associ w/. here
    energySuperClustersEl_->Fill(electronIter->superCluster()->energy());
    sizeSuperClustersEl_->Fill(electronIter->superCluster()->clustersSize());

    // loop over basicClusters to get energy
    double emaxSuperCluster = 0.;
    double phibar = 0.;
    double phi2bar = 0.;
    double eTotSuperCluster = 0.;

    for (reco::CaloCluster_iterator basicClusterIter = electronIter->superCluster()->clustersBegin();
         basicClusterIter != electronIter->superCluster()->clustersEnd();
         ++basicClusterIter) {
      //std::vector<DetId> theIds= (*basicClusterIter)->getHitsByDetId();

      double eCluster = (*basicClusterIter)->energy();
      if (eCluster > emaxSuperCluster) {
        emaxSuperCluster = eCluster;
      }
      eTotSuperCluster += eCluster;
      double phiCluster = (*basicClusterIter)->position().phi();
      phibar += eCluster * phiCluster;
      phi2bar += eCluster * phiCluster * phiCluster;
    }

    phibar = phibar / eTotSuperCluster;
    phi2bar = phi2bar / eTotSuperCluster;
    double phiWidth = phi2bar - phibar * phibar;
    if (phiWidth > 0.) {
      phiWidth = std::pow(phiWidth, 0.5);
    } else {
      phiWidth = 0.;
    }

    phiWidthSuperClustersEl_->Fill(phiWidth);

    emaxSuperClustersEl_->Fill(emaxSuperCluster);

    numberOfElectrons++;
  }

  numCand_->Fill(siStripElectronHandle->size());


  // Now loop over the electrons (ie the fitted things.)

  LogDebug("") << " About to check Electrons";

  edm::Handle<reco::ElectronCollection> electrons;
  iEvent.getByLabel(electronProducer_, electronCollection_, electrons);

  numElectrons_->Fill(electrons->end() - electrons->begin());

  // set up vector of bool for SiStrips having or not having Electrons
  // this causes a warning because of variable array size at compilation time ;
  // BAD  bool hasElectron_[siStripElectronHandle->end()- siStripElectronHandle->begin()] ;
  bool* hasElectron_ = new bool[siStripElectronHandle->end() - siStripElectronHandle->begin()];
  for (int icount = 0; icount < siStripElectronHandle->end() - siStripElectronHandle->begin(); ++icount) {
    hasElectron_[icount] = false;
  }

  // also set up a counter to associate the ith electron to the jth strippy
  // Electron_to_strippy[i] = j: i-th Electron is j-th strippy
  // BAD  unsigned int Electron_to_strippy[electrons->end()- electrons->begin()];
  unsigned int* Electron_to_strippy = new unsigned int[electrons->end() - electrons->begin()];
  for (int icount = 0; icount < electrons->end() - electrons->begin(); ++icount) {
    Electron_to_strippy[icount] = 0;
  }

  unsigned int ecount = 0;
  for (reco::ElectronCollection::const_iterator electronIter = electrons->begin(); electronIter != electrons->end();
       ++electronIter) {
    LogDebug("") << " Associating Electrons to Strippies ";
    LogDebug("") << " PT is " << electronIter->track()->pt();

    reco::TrackRef tr = (*electronIter).track();
    uint32_t id = (*electronIter->track()->recHitsBegin())->geographicalId().rawId();
    LocalPoint pos = (*electronIter->track()->recHitsBegin())->localPosition();

    unsigned int icount = 0;
    LogDebug("") << " About to loop over Strippies "
                 << " \n "
                 << " icount " << icount << " max " << siStripElectronHandle->end() - siStripElectronHandle->begin();

    for (reco::SiStripElectronCollection::const_iterator strippyiter = siStripElectronHandle->begin();
         strippyiter != siStripElectronHandle->end();
         ++strippyiter) {
      bool hitInCommon = false;
      // loop over rphi hits
      for (std::vector<SiStripRecHit2D>::const_iterator hiter = strippyiter->rphiRecHits().begin();
           hiter != strippyiter->rphiRecHits().end();
           ++hiter) {
        if (hiter->geographicalId().rawId() == id && (hiter->localPosition() - pos).mag() < 1e-10) {
          hitInCommon = true;
          break;
        }
      }

      for (std::vector<SiStripRecHit2D>::const_iterator hiter = strippyiter->stereoRecHits().begin();
           hiter != strippyiter->stereoRecHits().end();
           ++hiter) {
        if (hiter->geographicalId().rawId() == id && (hiter->localPosition() - pos).mag() < 1e-10) {
          hitInCommon = true;
          break;
        }
      }
      if (hitInCommon) {  //this Electron belongs to this SiStripElectron.
        hasElectron_[icount] = true;
        Electron_to_strippy[ecount] = icount;
        ptDiff->Fill(std::abs(electronIter->track()->pt()) - std::abs(strippyiter->pt()));
        pDiff->Fill(std::abs(electronIter->track()->p()) - std::abs(strippyiter->p()));
      }
      icount++;
    }  // Sistrip loop
    ecount++;
  }  // Electrons

  LogDebug("") << " Done looping over Electrons ";

  unsigned int counter = 0;
  for (reco::SiStripElectronCollection::const_iterator strippyIter = siStripElectronHandle->begin();
       strippyIter != siStripElectronHandle->end();
       ++strippyIter) {
    bool skipThis = !hasElectron_[counter];
    if (skipThis) {
      // plot stuff for SIStripElectrons that don't have fits associated

      LogDebug("") << " SiStrip Failed Electron "
                   << " \n "
                   << " p " << strippyIter->p() << " \n "
                   << " pt " << strippyIter->pt() << " \n "
                   << " SuperClust size " << strippyIter->superCluster()->clustersSize();

      pElectronFailed->Fill(std::abs(strippyIter->p()));
      ptElectronFailed->Fill(std::abs(strippyIter->pt()));
      sizeSuperClustersFailed->Fill(strippyIter->superCluster()->clustersSize());
      LogDebug("") << " done filling Failed histos ";
      //      energySuperClustersFailed->Fill(strippyIter->superCluster()->energy());
      //       if(strippyIter->p()>0.) {
      //    eOverPFailed->Fill(strippyIter->superCluster()->energy()/strippyIter->p());
      //       }else {
      //    eOverPFailed->Fill(-1.0);
      //       }

    } else {
      LogDebug("") << " SiStrip Passed Electron "
                   << " \n "
                   << " p " << strippyIter->p() << " \n "
                   << " pt " << strippyIter->pt() << " \n "
                   << " SuperClust size " << strippyIter->superCluster()->clustersSize();
      pElectronPassed->Fill(std::abs(strippyIter->p()));
      ptElectronPassed->Fill(std::abs(strippyIter->pt()));
      sizeSuperClustersPassed->Fill(strippyIter->superCluster()->clustersSize());
      LogDebug("") << " done filling passed histos ";
      //      energySuperClustersPassed->Fill(strippyIter->superCluster()->energy());
      //       if(strippyIter->p()>0.) {
      //    eOverPPassed->Fill(strippyIter->superCluster()->energy()/strippyIter->p());
      //       }else {
      //    eOverPPassed->Fill(-1.0);
      //       }

    }  // skipThis
    counter++;
  }

  LogDebug("") << "Dump info for all electrons ";

  for (reco::ElectronCollection::const_iterator electronIter1 = electrons->begin(); electronIter1 != electrons->end();
       ++electronIter1) {
    reco::TrackRef tr1 = (*electronIter1).track();
    // let's find its associated SiStripElectron and SuperCluster
    unsigned int ecount1 = electronIter1 - electrons->begin();
    unsigned int stripCount1 = 0;
    reco::SiStripElectronCollection::const_iterator strippyIter1;
    for (reco::SiStripElectronCollection::const_iterator strippyIter = siStripElectronHandle->begin();
         strippyIter != siStripElectronHandle->end();
         ++strippyIter) {
      if (Electron_to_strippy[ecount1] == stripCount1) {
        strippyIter1 = strippyIter;
        break;
      }
      stripCount1++;
    }  // strippy loop
    ecount1++;

    std::ostringstream str;

    str << " SiStripElect p , px, py, pz " << strippyIter1->p() << "  " << strippyIter1->px() << "  "
        << strippyIter1->py() << "  " << strippyIter1->pz() << "\n " << std::endl;

    str << " Electron p px, py, pz,  = " << tr1->p() << "  " << tr1->px() << "  " << tr1->py() << "  " << tr1->pz()
        << "\n"
        << std::endl;

    double EClust1 = strippyIter1->superCluster()->energy();
    double XClust1 = strippyIter1->superCluster()->x();
    double YClust1 = strippyIter1->superCluster()->y();
    double ZClust1 = strippyIter1->superCluster()->z();

    double rho1 = sqrt(XClust1 * XClust1 + YClust1 * YClust1 + ZClust1 * ZClust1);
    double costheta1 = ZClust1 / rho1;
    double sintheta1 = sqrt(1 - costheta1 * costheta1);
    if (ZClust1 < 0) {
      sintheta1 = -sintheta1;
    }
    double cosphi1 = XClust1 / sqrt(XClust1 * XClust1 + YClust1 * YClust1);
    double sinphi1 = YClust1 / sqrt(XClust1 * XClust1 + YClust1 * YClust1);

    str << " Ecal for electron E, px, py, pz " << EClust1 << " " << EClust1 * sintheta1 * cosphi1 << " "
        << EClust1 * sintheta1 * sinphi1 << " " << EClust1 * costheta1 << "\n"
        << std::endl;

    LogDebug("") << str.str();

  }  // loop over electrons
  LogDebug("") << "Done Dumping info for all electrons ";

  // LogDebug("")<< " Checking Electrons" ;
  //  LogDebug("")<< " PT is " << electronIter->track()->pt() ;
  //    reco::TrackRef tr =(*electronIter).track();
  if (electrons->end() - electrons->begin() > 1) {
    edm::LogInfo("") << " Two electrons in this event " << std::endl;
    for (reco::ElectronCollection::const_iterator electronIter1 = electrons->begin();
         electronIter1 != electrons->end() - 1;
         ++electronIter1) {
      reco::TrackRef tr1 = (*electronIter1).track();

      // let's find its associated SiStripElectron and SuperCluster
      // use the Electron_to_strippy[] array
      unsigned int ecount1 = electronIter1 - electrons->begin();
      // loop over strippies to find the corresponding one
      unsigned int stripCount1 = 0;
      reco::SiStripElectronCollection::const_iterator strippyIter1;
      for (reco::SiStripElectronCollection::const_iterator strippyIter = siStripElectronHandle->begin();
           strippyIter != siStripElectronHandle->end();
           ++strippyIter) {
        if (Electron_to_strippy[ecount1] == stripCount1) {
          strippyIter1 = strippyIter;
          break;
        }
        stripCount1++;
      }  // strippy loop

      double EClust1 = strippyIter1->superCluster()->energy();
      double XClust1 = strippyIter1->superCluster()->x();
      double YClust1 = strippyIter1->superCluster()->y();
      double ZClust1 = strippyIter1->superCluster()->z();

      for (reco::ElectronCollection::const_iterator electronIter2 = electronIter1 + 1;
           electronIter2 != electrons->end();
           ++electronIter2) {
        reco::TrackRef tr2 = (*electronIter2).track();

        unsigned int ecount2 = electronIter2 - electrons->begin();
        unsigned int stripCount2 = 0;
        reco::SiStripElectronCollection::const_iterator strippyIter2;
        for (reco::SiStripElectronCollection::const_iterator strippyIter = siStripElectronHandle->begin();
             strippyIter != siStripElectronHandle->end();
             ++strippyIter) {
          if (Electron_to_strippy[ecount2] == stripCount2) {
            strippyIter2 = strippyIter;
            break;
          }
          stripCount2++;
        }  // strippy loop

        double EClust2 = strippyIter2->superCluster()->energy();
        double XClust2 = strippyIter2->superCluster()->x();
        double YClust2 = strippyIter2->superCluster()->y();
        double ZClust2 = strippyIter2->superCluster()->z();

        // now get supercluster from this:

        edm::LogInfo("") << " Electron p1 = " << tr1->p() << " p1x " << tr1->px() << " p1y " << tr1->py() << " p1z "
                         << tr1->pz() << std::endl;

        edm::LogInfo("") << " Electron p2 = " << tr2->p() << " p2x " << tr2->px() << " p2y " << tr2->py() << " p2z "
                         << tr2->pz() << std::endl;

        // combine the two in an (e,e) pair
        double Zpx = tr1->px() + tr2->px();
        double Zpy = tr1->py() + tr2->py();
        double Zpz = tr1->pz() + tr2->pz();
        double Ze = std::abs(tr1->p()) + std::abs(tr2->p());
        edm::LogInfo("") << " Z mass " << sqrt(Ze * Ze - Zpx * Zpx - Zpy * Zpy - Zpz * Zpz) << std::endl;

        // combine the SuperClusts into a Z
        double rho1 = sqrt(XClust1 * XClust1 + YClust1 * YClust1 + ZClust1 * ZClust1);
        double costheta1 = ZClust1 / rho1;
        double sintheta1 = sqrt(1 - costheta1 * costheta1);
        if (ZClust1 < 0) {
          sintheta1 = -sintheta1;
        }
        double cosphi1 = XClust1 / sqrt(XClust1 * XClust1 + YClust1 * YClust1);
        double sinphi1 = YClust1 / sqrt(XClust1 * XClust1 + YClust1 * YClust1);

        double rho2 = sqrt(XClust2 * XClust2 + YClust2 * YClust2 + ZClust2 * ZClust2);
        double costheta2 = ZClust2 / rho2;
        double sintheta2 = sqrt(1 - costheta2 * costheta2);
        if (ZClust2 < 0) {
          sintheta2 = -sintheta2;
        }
        double cosphi2 = XClust2 / sqrt(XClust2 * XClust2 + YClust2 * YClust2);
        double sinphi2 = YClust2 / sqrt(XClust2 * XClust2 + YClust2 * YClust2);

        edm::LogInfo("") << "Energy of supercluster for 1st electron " << EClust1 << " "
                         << EClust1 * sintheta1 * cosphi1 << " " << EClust1 * sintheta1 * sinphi1 << " "
                         << EClust1 * costheta1 << " " << std::endl;

        edm::LogInfo("") << "Energy of supercluster for 2nd electron " << EClust2 << " "
                         << EClust2 * sintheta2 * cosphi2 << " " << EClust2 * sintheta2 * sinphi2 << " "
                         << EClust2 * costheta2 << " " << std::endl;

        // get the supercluster pair
        double Zgpx = EClust1 * sintheta1 * cosphi1 + EClust2 * sintheta2 * cosphi2;
        double Zgpy = EClust1 * sintheta1 * sinphi1 + EClust2 * sintheta2 * sinphi2;
        double Zgpz = EClust1 * costheta1 + EClust2 * costheta2;
        double ZgE = EClust1 + EClust2;

        edm::LogInfo("") << " Z mass from ECAL " << sqrt(ZgE * ZgE - Zgpx * Zgpx - Zgpy * Zgpy - Zgpz * Zgpz)
                         << std::endl;

      }  //inner loop
    }    // outer loop
  }      // m(ee) loop

  delete[] hasElectron_;
  delete[] Electron_to_strippy;

  LogDebug("") << " About to dump tracker info ";

  edm::ESHandle<TrackerGeometry> trackerHandle = iSetup.getHandle(trackerToken_);

  edm::Handle<SiStripRecHit2DCollection> rphiHitsHandle;
  iEvent.getByLabel(siHitProducer_, siRphiHitCollection_, rphiHitsHandle);

  edm::Handle<SiStripRecHit2DCollection> stereoHitsHandle;
  iEvent.getByLabel(siHitProducer_, siStereoHitCollection_, stereoHitsHandle);

  edm::Handle<SiStripMatchedRecHit2DCollection> matchedHitsHandle;
  iEvent.getByLabel(siHitProducer_, siMatchedHitCollection_, matchedHitsHandle);


  NShowers_ = 0;
  for (reco::SuperClusterCollection::const_iterator clusterIter = clusterHandle->begin();
       clusterIter != clusterHandle->end();
       ++clusterIter) {
    double energy = clusterIter->energy();
    math::XYZPoint position = clusterIter->position();
    if (NShowers_ < myMaxHits) {
      EShower_[NShowers_] = energy;
      XShower_[NShowers_] = position.x();
      YShower_[NShowers_] = position.y();
      ZShower_[NShowers_] = position.z();
      ++NShowers_;
    }
    // Loop over all crystals in this supercluster - see
    // RecoEcal/EgamaClusterProducers/src/EgammaSimpleAnalyzer.cc
    // Look also at DataFormats/EgammaReco/interface/SuperCluster.h
  }
  numSuperClusters_->Fill(NShowers_);

  LogDebug("") << " Looping over stereo hits ";

  int myHits = 0;
  for (SiStripRecHit2DCollection::DataContainer::const_iterator hit = stereoHitsHandle->data().begin(),
                                                                hitend = stereoHitsHandle->data().end();
       hit != hitend;
       ++hit) {
    DetId id(hit->geographicalId());
    if ((hit->geographicalId()).subdetId() == StripSubdetector::TIB ||
        (hit->geographicalId()).subdetId() == StripSubdetector::TOB) {
      GlobalPoint position = trackerHandle->idToDet(hit->geographicalId())->surface().toGlobal(hit->localPosition());
      //from RecoLocalTracker/SiStripClusterizer/test/TestCluster.cc
      // cf also TrackHitAssociator.cc SiStripRecHitMatcher.cc SiStrip1DMeasurementTransformator.cc (KalmanUpdators)
      SiStripRecHit2D const rechit = *hit;
      //    LocalPoint myposition = rechit.localPosition() ;
      LocalError myerror = rechit.localPositionError();

      // Get layer and subdetector ID here for this hit
      // see SiStripRecHitConverter/test/ValHit.cc
      Int_t siLayerNum = 0;
      Int_t siDetNum = 0;
      string siDetName = "";
      if ((hit->geographicalId()).subdetId() == StripSubdetector::TIB) {
        //     siLayerNum = tTopo->tibLayer(rechit->geographicalID());
        siLayerNum = tTopo->tibLayer(id);
        siDetNum = 1;
        siDetName = "TIB";
      } else if ((hit->geographicalId()).subdetId() == StripSubdetector::TOB) {
        siLayerNum = tTopo->tobLayer(id);
        siDetNum = 2;
        siDetName = "TOB";
        //      } else if ( (hit->geographicalId()).subdetId() == StripSubdetector::TID ){
        //    // should we use side/wheel/ring/module/stereo() ?
        //    siLayerNum = tTopo->tidWheel(id);
        //    siDetNum = 3 ;
        //    siDetName = "TID" ;
        //  }else if ( (hit->geographicalId()).subdetId() == StripSubdetector::TEC ){
        //    //choices are side/petal/wheel/ring/module/glued/stereo
        //    siLayerNum = tTopo->tecWheel(id);
        //    siDetNum = 4 ;
        //    siDetName = "TEC" ;
      } else {
        siLayerNum = -999;
        siDetNum = -999;
        siDetName = "NULL";
      }
      //    LogDebug("") << siDetName << " " << siLayerNum ;

      const SiStripRecHit2D::ClusterRef& clust = rechit.cluster();
      double Signal = 0;
      double Noise2 = 0;
      int StripCount = 0;
      if (clust.isNonnull()) {
        //    LogDebug("") << " barycenter " << clust->barycenter() ;
        //    const std::vector<uint16_t> amplitudes=clust->amplitudes();
        const auto& amplitudes = clust->amplitudes();
        for (size_t i = 0; i < amplitudes.size(); i++) {
          Signal += amplitudes[i];
          //ignore for now       Noise2 +=SiStripNoiseService_.getNoise(detid,clust->firstStrip()+i)*SiStripNoiseService_.getNoise(detid,clust->firstStrip()+i);
          StripCount++;
        }
      } else {
        LogDebug("") << " null cluster ";
      }
      //    LogDebug("") << "Signal " << Signal << " Noise2 " << Noise2 << " StripCount " << StripCount ;
      // Dump position
      //    LogDebug("") << " Stereo "
      //             << "local position: "<<myposition.x()<<" "
      //             << myposition.y()<<" "<<myposition.z()<<"\n"
      //             << "local error: "<<myerror.xx()<<" "
      //             << myerror.xy()<<" "<<myerror.yy() << "\n"
      //             << "global position: " << position.x() << " "
      //             <<  position.y()<<" "<< position.z()<<"\n"
      //             << " siDetNum " << siDetNum
      //             << " siLayerNum " << siLayerNum ;

      if (myHits < myMaxHits) {
        StereoHitX_[myHits] = position.x();
        StereoHitY_[myHits] = position.y();
        StereoHitZ_[myHits] = position.z();

        StereoHitR_[myHits] = position.perp();
        StereoHitPhi_[myHits] = position.phi();
        StereoHitTheta_[myHits] = position.theta();

        StereoHitSigX_[myHits] = sqrt(myerror.xx());
        StereoHitSigY_[myHits] = sqrt(myerror.yy());
        StereoHitCorr_[myHits] = myerror.xy() / sqrt(myerror.xx() * myerror.yy());

        StereoHitSignal_[myHits] = Signal;
        StereoHitNoise_[myHits] = Noise2;
        StereoHitWidth_[myHits] = StripCount;

        StereoDetector_[myHits] = siDetNum;
        StereoLayer_[myHits] = siLayerNum;

        ++myHits;
      }
    }  // end if this is the right subdetector
  }    // end loop over hits
  NStereoHits_ = myHits;

  numSiStereoHits_->Fill(NStereoHits_);


  LogDebug("") << " Looping over Mono Hits ";
  myHits = 0;
  for (SiStripRecHit2DCollection::DataContainer::const_iterator hit = rphiHitsHandle->data().begin(),
                                                                hitend = rphiHitsHandle->data().end();
       hit != hitend;
       ++hit) {
    DetId id(hit->geographicalId());

    if ((hit->geographicalId()).subdetId() == StripSubdetector::TIB ||
        (hit->geographicalId()).subdetId() == StripSubdetector::TOB) {
      GlobalPoint position = trackerHandle->idToDet(hit->geographicalId())->surface().toGlobal(hit->localPosition());
      //from RecoLocalTracker/SiStripClusterizer/test/TestCluster.cc
      // cf also TrackHitAssociator.cc SiStripRecHitMatcher.cc SiStrip1DMeasurementTransformator.cc (KalmanUpdators)
      SiStripRecHit2D const rechit = *hit;
      //    LocalPoint myposition = rechit.localPosition() ;
      LocalError myerror = rechit.localPositionError();

      // Get layer and subdetector ID here for this hit
      // see SiStripRecHitConverter/test/ValHit.cc
      Int_t siLayerNum = 0;
      Int_t siDetNum = 0;
      string siDetName = "";
      if ((hit->geographicalId()).subdetId() == StripSubdetector::TIB) {
        //     siLayerNum = tTopo->tibLayer(rechit->geographicalID());
        siLayerNum = tTopo->tibLayer(id);
        siDetNum = 1;
        siDetName = "TIB";
      } else if ((hit->geographicalId()).subdetId() == StripSubdetector::TOB) {
        siLayerNum = tTopo->tobLayer(id);
        siDetNum = 2;
        siDetName = "TOB";
        //  } else if ( (hit->geographicalId()).subdetId() == StripSubdetector::TID ){
        //    // should we use side/wheel/ring/module/stereo() ?
        //    siLayerNum = tTopo->tidWheel(id);
        //    siDetNum = 3 ;
        //    siDetName = "TID" ;
        //  }else if ( (hit->geographicalId()).subdetId() == StripSubdetector::TEC ){
        //    //choices are side/petal/wheel/ring/module/glued/stereo
        //    siLayerNum = tTopo->tecWheel(id);
        //    siDetNum = 4 ;
        //    siDetName = "TEC"
        ;
      } else {
        siLayerNum = -999;
        siDetNum = -999;
        siDetName = "NULL";
      }
      //    LogDebug("") << siDetName << " " << siLayerNum ;
      const SiStripRecHit2D::ClusterRef& clust = rechit.cluster();
      double Signal = 0;
      double Noise2 = 0;
      int StripCount = 0;
      if (clust.isNonnull()) {
        //    LogDebug("") << " barycenter " << clust->barycenter() ;
        //    const std::vector<uint16_t> amplitudes=clust->amplitudes();
        const auto& amplitudes = clust->amplitudes();
        for (size_t i = 0; i < amplitudes.size(); i++) {
          Signal += amplitudes[i];
          //ignore for now       Noise2 +=SiStripNoiseService_.getNoise(detid,clust->firstStrip()+i)*SiStripNoiseService_.getNoise(detid,clust->firstStrip()+i);
          StripCount++;
        }
      } else {
        LogDebug("") << " null cluster ";
      }
      //    LogDebug("") << "Signal " << Signal << " Noise2 " << Noise2 << " StripCount " << StripCount ;

      // Dump position info
      //    LogDebug("") << " Mono "
      //             << "local position: "<<myposition.x()<<" "
      //             << myposition.y()<<" "<<myposition.z()<<"\n"
      //             <<"local error: "<<myerror.xx()<<" "
      //             << myerror.xy()<<" "<<myerror.yy() << "\n"
      //             << "global position: " << position.x() << " "
      //             <<  position.y()<<" "<< position.z()<<"\n"
      //             << " siDetNum " << siDetNum
      //             << " siLayerNum " << siLayerNum ;

      if (myHits < myMaxHits) {
        MonoHitX_[myHits] = position.x();
        MonoHitY_[myHits] = position.y();
        MonoHitZ_[myHits] = position.z();

        MonoHitR_[myHits] = position.perp();
        MonoHitPhi_[myHits] = position.phi();
        MonoHitTheta_[myHits] = position.theta();

        MonoHitSigX_[myHits] = sqrt(myerror.xx());
        MonoHitSigY_[myHits] = sqrt(myerror.yy());
        MonoHitCorr_[myHits] = myerror.xy() / sqrt(myerror.xx() * myerror.yy());

        MonoHitSignal_[myHits] = Signal;
        MonoHitNoise_[myHits] = Noise2;
        MonoHitWidth_[myHits] = StripCount;

        MonoDetector_[myHits] = siDetNum;
        MonoLayer_[myHits] = siLayerNum;

        ++myHits;
      }  // of  if(myHits < myMaxHits)
         // LogDebug("")<< "end of myHits < myMaxHits " ;
    }    // end if this is the right subdetector
         //      LogDebug("")<< "end of TIB/TOB check " ;
  }      // end loop over hits
  //    LogDebug("")<< " end of loop over hits  " ;
  NMonoHits_ = myHits;

  numSiMonoHits_->Fill(NMonoHits_);


  LogDebug("") << "  Loop over Matched Hits ";

  myHits = 0;
  for (SiStripMatchedRecHit2DCollection::DataContainer::const_iterator hit = matchedHitsHandle->data().begin(),
                                                                       hitend = matchedHitsHandle->data().end();
       hit != hitend;
       ++hit) {
    DetId id(hit->geographicalId());
    if ((hit->geographicalId()).subdetId() == StripSubdetector::TIB ||
        (hit->geographicalId()).subdetId() == StripSubdetector::TOB) {
      GlobalPoint position = trackerHandle->idToDet(hit->geographicalId())->surface().toGlobal(hit->localPosition());
      SiStripMatchedRecHit2D const rechit = *hit;
      //    LocalPoint myposition = rechit.localPosition() ;
      LocalError myerror = rechit.localPositionError();

      // Get layer and subdetector ID here for this hit
      // see SiStripRecHitConverter/test/ValHit.cc
      Int_t siLayerNum = 0;
      Int_t siDetNum = 0;
      string siDetName = "";
      if ((hit->geographicalId()).subdetId() == StripSubdetector::TIB) {
        siLayerNum = tTopo->tibLayer(id);
        siDetNum = 1;
        siDetName = "TIB";
      } else if ((hit->geographicalId()).subdetId() == StripSubdetector::TOB) {
        siLayerNum = tTopo->tobLayer(id);
        siDetNum = 2;
        siDetName = "TOB";
        //      } else if ( (hit->geographicalId()).subdetId() == StripSubdetector::TID ){
        //        // should we use side/wheel/ring/module/stereo() ?
        //        siLayerNum = tTopo->tidWheel(id);
        //        siDetNum = 3 ;
        //        siDetName = "TID" ;
        //      }else if ( (hit->geographicalId()).subdetId() == StripSubdetector::TEC ){
        //        //choices are side/petal/wheel/ring/module/glued/stereo
        //        siLayerNum = tTopo->tecWheel(id);
        //        siDetNum = 4 ;
        //        siDetName = "TEC" ;
      } else {
        siLayerNum = -999;
        siDetNum = -999;
        siDetName = "NULL";
      }
      //    const edm::Ref<edm::DetSetVector<SiStripCluster>, SiStripCluster, edm::refhelper::FindForDetSetVector<SiStripCluster> > clust=rechit.cluster();
      double Signal = 0;
      double Noise2 = 0;
      int StripCount = 0;
      //    if(clust.isNonnull()) {
      //      LogDebug("") << " barycenter " << clust->barycenter() ;
      //      const std::vector<uint16_t> amplitudes=clust->amplitudes();
      //      for(size_t i = 0 ; i<amplitudes.size(); i++ ){
      //        Signal +=amplitudes[i] ;
      //        //ignore for now         Noise2 +=SiStripNoiseService_.getNoise(detid,clust->firstStrip()+i)*SiStripNoiseService_.getNoise(detid,clust->firstStrip()+i);
      //        StripCount++;
      //      }
      //    } else {
      //      LogDebug("") << " null cluster " ;
      //    }
      //    LogDebug("") << "Signal " << Signal << " Noise2 " << Noise2 << " StripCount " << StripCount ;

      // Dump position info
      //    LogDebug("") << " Matched "
      //             << "local position: "<<myposition.x()<<" "
      //             << myposition.y()<<" "<<myposition.z()<<"\n"
      //             << "local error: "<<myerror.xx()<<" "
      //             << myerror.xy()<<" "<<myerror.yy() << "\n"
      //             << "global position: " << position.x() << " "
      //             <<  position.y()<<" "<< position.z()<<"\n"
      //             << " siDetNum " << siDetNum
      //             << " siLayerNum " << siLayerNum ;

      if (myHits < myMaxHits) {
        MatchedHitX_[myHits] = position.x();
        MatchedHitY_[myHits] = position.y();
        MatchedHitZ_[myHits] = position.z();

        MatchedHitR_[myHits] = position.perp();
        MatchedHitPhi_[myHits] = position.phi();
        MatchedHitTheta_[myHits] = position.theta();

        MatchedHitSigX_[myHits] = sqrt(myerror.xx());
        MatchedHitSigY_[myHits] = sqrt(myerror.yy());
        MatchedHitCorr_[myHits] = myerror.xy() / sqrt(myerror.xx() * myerror.yy());

        MatchedHitSignal_[myHits] = Signal;
        MatchedHitNoise_[myHits] = Noise2;
        MatchedHitWidth_[myHits] = StripCount;

        MatchedDetector_[myHits] = siDetNum;
        MatchedLayer_[myHits] = siLayerNum;

        ++myHits;
      }
    }  // end if this is the right subdetector (TIB/TOB)
  }    // end loop over hits
  NMatchedHits_ = myHits;

  numSiMatchedHits_->Fill(NMatchedHits_);


  LogDebug("") << "Writing to myTree with " << NShowers_ << " Showers " << NStereoHits_ << " Si StereoHits "
               << NMonoHits_ << " Si MonoHits " << NMatchedHits_ << " Si MatchedHits ";

  myTree_->Fill();

}  // end of Analyzer

void SiStripElectronAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 305 of file SiStripElectronAnalyzer.cc.

References EShower_, MatchedDetector_, MatchedHitCorr_, MatchedHitNoise_, MatchedHitPhi_, MatchedHitR_, MatchedHitSignal_, MatchedHitSigX_, MatchedHitSigY_, MatchedHitTheta_, MatchedHitWidth_, MatchedHitX_, MatchedHitY_, MatchedHitZ_, MatchedLayer_, MonoDetector_, MonoHitCorr_, MonoHitNoise_, MonoHitPhi_, MonoHitR_, MonoHitSignal_, MonoHitSigX_, MonoHitSigY_, MonoHitTheta_, MonoHitWidth_, MonoHitX_, MonoHitY_, MonoHitZ_, MonoLayer_, myTree_, NMatchedHits_, NMonoHits_, NShowers_, NStereoHits_, StereoDetector_, StereoHitCorr_, StereoHitNoise_, StereoHitPhi_, StereoHitR_, StereoHitSignal_, StereoHitSigX_, StereoHitSigY_, StereoHitTheta_, StereoHitWidth_, StereoHitX_, StereoHitY_, StereoHitZ_, StereoLayer_, XShower_, YShower_, and ZShower_.

                                       {
  myTree_ = new TTree("myTree", "my first Tree example");

  myTree_->Branch("NShowers", &NShowers_, "NShowers/I");

  // first specify the ECAL clusters
  // need to explicitly include array length.
  myTree_->Branch("EShower", &EShower_, "EShower[1000]/F");
  myTree_->Branch("XShower", &XShower_, "XShower[1000]/F");
  myTree_->Branch("YShower", &YShower_, "YShower[1000]/F");
  myTree_->Branch("ZShower", &ZShower_, "ZShower[1000]/F");

  // second specify the Si Stereo Hits
  myTree_->Branch("NStereoHits", &NStereoHits_, "NStereoHits/I");
  myTree_->Branch("StereoHitX", &StereoHitX_, "StereoHitX[1000]/F");
  myTree_->Branch("StereoHitY", &StereoHitY_, "StereoHitY[1000]/F");
  myTree_->Branch("StereoHitZ", &StereoHitZ_, "StereoHitZ[1000]/F");

  myTree_->Branch("StereoHitR", &StereoHitR_, "StereoHitR[1000]/F");
  myTree_->Branch("StereoHitPhi", &StereoHitPhi_, "StereoHitPhi[1000]/F");
  myTree_->Branch("StereoHitTheta", &StereoHitTheta_, "StereoHitTheta[1000]/F");

  myTree_->Branch("StereoHitSigX", &StereoHitSigX_, "StereoHitSigX[1000]/F");
  myTree_->Branch("StereoHitSigY", &StereoHitSigY_, "StereoHitSigY[1000]/F");
  myTree_->Branch("StereoHitCorr", &StereoHitCorr_, "StereoHitCorr[1000]/F");

  myTree_->Branch("StereoHitSignal", &StereoHitSignal_, "StereoHitSignal[1000]/F");
  myTree_->Branch("StereoHitNoise", &StereoHitNoise_, "StereoHitNoise[1000]/F");
  myTree_->Branch("StereoHitWidth", &StereoHitWidth_, "StereoHitWidth[1000]/I");

  myTree_->Branch("StereoDetector", &StereoDetector_, "StereoDetector[1000]/I");
  myTree_->Branch("StereoLayer", &StereoLayer_, "StereoLayer[1000]/I");

  // specify the Si mono (rphi) hits
  myTree_->Branch("NMonoHits", &NMonoHits_, "NMonoHits/I");
  myTree_->Branch("MonoHitX", &MonoHitX_, "MonoHitX[1000]/F");
  myTree_->Branch("MonoHitY", &MonoHitY_, "MonoHitY[1000]/F");
  myTree_->Branch("MonoHitZ", &MonoHitZ_, "MonoHitZ[1000]/F");

  myTree_->Branch("MonoHitR", &MonoHitR_, "MonoHitR[1000]/F");
  myTree_->Branch("MonoHitPhi", &MonoHitPhi_, "MonoHitPhi[1000]/F");
  myTree_->Branch("MonoHitTheta", &MonoHitTheta_, "MonoHitTheta[1000]/F");

  myTree_->Branch("MonoHitSigX", &MonoHitSigX_, "MonoHitSigX[1000]/F");
  myTree_->Branch("MonoHitSigY", &MonoHitSigY_, "MonoHitSigY[1000]/F");
  myTree_->Branch("MonoHitCorr", &MonoHitCorr_, "MonoHitCorr[1000]/F");

  myTree_->Branch("MonoHitSignal", &MonoHitSignal_, "MonoHitSignal[1000]/F");
  myTree_->Branch("MonoHitNoise", &MonoHitNoise_, "MonoHitNoise[1000]/F");
  myTree_->Branch("MonoHitWidth", &MonoHitWidth_, "MonoHitWidth[1000]/I");

  myTree_->Branch("MonoDetector", &MonoDetector_, "MonoDetector[1000]/I");
  myTree_->Branch("MonoLayer", &MonoLayer_, "MonoLayer[1000]/I");

  // specify the Si matched (rphi) hits
  myTree_->Branch("NMatchedHits", &NMatchedHits_, "NMatchedHits/I");
  myTree_->Branch("MatchedHitX", &MatchedHitX_, "MatchedHitX[1000]/F");
  myTree_->Branch("MatchedHitY", &MatchedHitY_, "MatchedHitY[1000]/F");
  myTree_->Branch("MatchedHitZ", &MatchedHitZ_, "MatchedHitZ[1000]/F");

  myTree_->Branch("MatchedHitR", &MatchedHitR_, "MatchedHitR[1000]/F");
  myTree_->Branch("MatchedHitPhi", &MatchedHitPhi_, "MatchedHitPhi[1000]/F");
  myTree_->Branch("MatchedHitTheta", &MatchedHitTheta_, "MatchedHitTheta[1000]/F");

  myTree_->Branch("MatchedHitSigX", &MatchedHitSigX_, "MatchedHitSigX[1000]/F");
  myTree_->Branch("MatchedHitSigY", &MatchedHitSigY_, "MatchedHitSigY[1000]/F");
  myTree_->Branch("MatchedHitCorr", &MatchedHitCorr_, "MatchedHitCorr[1000]/F");

  myTree_->Branch("MatchedHitSignal", &MatchedHitSignal_, "MatchedHitSignal[1000]/F");
  myTree_->Branch("MatchedHitNoise", &MatchedHitNoise_, "MatchedHitNoise[1000]/F");
  myTree_->Branch("MatchedHitWidth", &MatchedHitWidth_, "MatchedHitWidth[1000]/I");

  myTree_->Branch("MatchedDetector", &MatchedDetector_, "MatchedDetector[1000]/I");
  myTree_->Branch("MatchedLayer", &MatchedLayer_, "MatchedLayer[1000]/I");
}

void SiStripElectronAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 1295 of file SiStripElectronAnalyzer.cc.

References emaxSuperClusters_, emaxSuperClustersEl_, energySuperClusters_, energySuperClustersEl_, file_, LogDebug, myTree_, numCand_, numElectrons_, numSiMatchedHits_, numSiMonoHits_, numSiStereoHits_, numSuperClusters_, pDiff, pElectronFailed, pElectronPassed, phiWidthSuperClusters_, phiWidthSuperClustersEl_, ptDiff, ptElectronFailed, ptElectronPassed, sizeSuperClusters_, sizeSuperClustersEl_, sizeSuperClustersFailed, and sizeSuperClustersPassed.

                                     {
  LogDebug("") << "Entering endJob ";
  file_->cd();
  numCand_->Write();
  numElectrons_->Write();
  numSuperClusters_->Write();

  energySuperClusters_->Write();
  sizeSuperClusters_->Write();
  emaxSuperClusters_->Write();
  phiWidthSuperClusters_->Write();

  energySuperClustersEl_->Write();
  sizeSuperClustersEl_->Write();
  emaxSuperClustersEl_->Write();
  phiWidthSuperClustersEl_->Write();

  ptDiff->Write();
  pDiff->Write();
  pElectronFailed->Write();
  ptElectronFailed->Write();
  pElectronPassed->Write();
  ptElectronPassed->Write();
  sizeSuperClustersPassed->Write();
  sizeSuperClustersFailed->Write();
  //   energySuperClustersPassed->Write();
  //   energySuperClustersFailed->Write();
  //   eOverPPassed->Write();
  //   eOverPFailed->Write();

  numSiStereoHits_->Write();
  numSiMonoHits_->Write();
  numSiMatchedHits_->Write();

  // disable for large dataset
  LogDebug("") << " Writing out ntuple is disabled for now ";
  myTree_->Write();

  file_->Close();
}

void SiStripElectronAnalyzer::initNtuple ( void  )

virtual

Definition at line 381 of file SiStripElectronAnalyzer.cc.

References EShower_, init, LogDebug, MatchedHitNoise_, MatchedHitPhi_, MatchedHitR_, MatchedHitSignal_, MatchedHitTheta_, MatchedHitWidth_, MatchedHitX_, MatchedHitY_, MatchedHitZ_, MonoHitNoise_, MonoHitPhi_, MonoHitR_, MonoHitSignal_, MonoHitTheta_, MonoHitWidth_, MonoHitX_, MonoHitY_, MonoHitZ_, myMaxHits, NMatchedHits_, NMonoHits_, NShowers_, NStereoHits_, StereoHitNoise_, StereoHitPhi_, StereoHitR_, StereoHitSignal_, StereoHitTheta_, StereoHitWidth_, StereoHitX_, StereoHitY_, StereoHitZ_, XShower_, YShower_, and ZShower_.

Referenced by analyze().

                                         {
  LogDebug("") << " In initNtuple ";

  NShowers_ = -999;
  for (int init = 0; init < myMaxHits; ++init) {
    EShower_[init] = -999.;
    XShower_[init] = -999.;
    YShower_[init] = -999.;
    ZShower_[init] = -999.;
  }
  NStereoHits_ = -999;

  for (int init = 0; init < myMaxHits; ++init) {
    StereoHitX_[init] = -999.;
    StereoHitY_[init] = -999.;
    StereoHitZ_[init] = -999.;
    StereoHitR_[init] = -999.;
    StereoHitPhi_[init] = -999.;
    StereoHitTheta_[init] = -999.;

    StereoHitSignal_[init] = -999.;
    StereoHitNoise_[init] = -999.;
    StereoHitWidth_[init] = -999;
    ;
  }

  NMonoHits_ = -999;
  for (int init = 0; init < myMaxHits; ++init) {
    MonoHitX_[init] = -999.;
    MonoHitY_[init] = -999.;
    MonoHitZ_[init] = -999.;
    MonoHitR_[init] = -999.;
    MonoHitPhi_[init] = -999.;
    MonoHitTheta_[init] = -999.;

    MonoHitSignal_[init] = -999.;
    MonoHitNoise_[init] = -999.;
    MonoHitWidth_[init] = -999;
    ;
  }

  NMatchedHits_ = -999;
  for (int init = 0; init < myMaxHits; ++init) {
    MatchedHitX_[init] = -999.;
    MatchedHitY_[init] = -999.;
    MatchedHitZ_[init] = -999.;
    MatchedHitR_[init] = -999.;
    MatchedHitPhi_[init] = -999.;
    MatchedHitTheta_[init] = -999.;

    MatchedHitSignal_[init] = -999.;
    MatchedHitNoise_[init] = -999.;
    MatchedHitWidth_[init] = -999;
    ;
  }
}

double SiStripElectronAnalyzer::unwrapPhi ( double  phi ) const

inlineprivate

Definition at line 68 of file SiStripElectronAnalyzer.cc.

References M_PI, and phi.

                                     {
    while (phi > M_PI) {
      phi -= 2. * M_PI;
    }
    while (phi < -M_PI) {
      phi += 2. * M_PI;
    }
    return phi;
  }

Member Data Documentation

std::string SiStripElectronAnalyzer::basicClusterCollection_

private

Definition at line 193 of file SiStripElectronAnalyzer.cc.

std::string SiStripElectronAnalyzer::basicClusterProducer_

private

Definition at line 192 of file SiStripElectronAnalyzer.cc.

std::string SiStripElectronAnalyzer::eBRecHitCollection_

private

Definition at line 195 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::eBRecHitProducer_

private

Definition at line 194 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::electronCollection_

private

Definition at line 199 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::electronProducer_

private

Definition at line 198 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::emaxSuperClusters_

private

Definition at line 90 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::emaxSuperClustersEl_

private

Definition at line 95 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::energySuperClusters_

private

Definition at line 88 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::energySuperClustersEl_

private

Definition at line 93 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::energySuperClustersFailed

private

Definition at line 107 of file SiStripElectronAnalyzer.cc.

Referenced by SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::energySuperClustersPassed

private

Definition at line 106 of file SiStripElectronAnalyzer.cc.

Referenced by SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::eOverPFailed

private

Definition at line 108 of file SiStripElectronAnalyzer.cc.

Referenced by SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::eOverPPassed

private

Definition at line 109 of file SiStripElectronAnalyzer.cc.

Referenced by SiStripElectronAnalyzer().

float SiStripElectronAnalyzer::EShower_[1000]

private

Definition at line 118 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

TFile* SiStripElectronAnalyzer::file_

private

Definition at line 84 of file SiStripElectronAnalyzer.cc.

Referenced by endJob(), SiStripElectronAnalyzer(), and ~SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::fileName_

private

Definition at line 82 of file SiStripElectronAnalyzer.cc.

Referenced by SiStripElectronAnalyzer().

int SiStripElectronAnalyzer::MatchedDetector_[1000]

private

Definition at line 185 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitCorr_[1000]

private

Definition at line 179 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitNoise_[1000]

private

Definition at line 182 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MatchedHitPhi_[1000]

private

Definition at line 173 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MatchedHitR_[1000]

private

Definition at line 172 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MatchedHitSignal_[1000]

private

Definition at line 181 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MatchedHitSigX_[1000]

private

Definition at line 177 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitSigY_[1000]

private

Definition at line 178 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitTheta_[1000]

private

Definition at line 174 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::MatchedHitWidth_[1000]

private

Definition at line 183 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MatchedHitX_[1000]

private

Definition at line 168 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MatchedHitY_[1000]

private

Definition at line 169 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MatchedHitZ_[1000]

private

Definition at line 170 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::MatchedLayer_[1000]

private

Definition at line 186 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::string SiStripElectronAnalyzer::mctruthCollection_

private

Definition at line 189 of file SiStripElectronAnalyzer.cc.

Referenced by SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::mctruthProducer_

private

Definition at line 188 of file SiStripElectronAnalyzer.cc.

Referenced by SiStripElectronAnalyzer().

int SiStripElectronAnalyzer::MonoDetector_[1000]

private

Definition at line 163 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitCorr_[1000]

private

Definition at line 157 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitNoise_[1000]

private

Definition at line 160 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MonoHitPhi_[1000]

private

Definition at line 151 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MonoHitR_[1000]

private

Definition at line 150 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MonoHitSignal_[1000]

private

Definition at line 159 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MonoHitSigX_[1000]

private

Definition at line 155 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitSigY_[1000]

private

Definition at line 156 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitTheta_[1000]

private

Definition at line 152 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::MonoHitWidth_[1000]

private

Definition at line 161 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MonoHitX_[1000]

private

Definition at line 146 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MonoHitY_[1000]

private

Definition at line 147 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::MonoHitZ_[1000]

private

Definition at line 148 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::MonoLayer_[1000]

private

Definition at line 164 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

TTree* SiStripElectronAnalyzer::myTree_

private

Definition at line 115 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and endJob().

int SiStripElectronAnalyzer::NMatchedHits_

private

Definition at line 167 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::NMonoHits_

private

Definition at line 145 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::NShowers_

private

Definition at line 117 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::NStereoHits_

private

Definition at line 123 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

TH1F* SiStripElectronAnalyzer::numCand_

private

Definition at line 85 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::numElectrons_

private

Definition at line 86 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::numSiMatchedHits_

private

Definition at line 113 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::numSiMonoHits_

private

Definition at line 112 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::numSiStereoHits_

private

Definition at line 111 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::numSuperClusters_

private

Definition at line 87 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::pDiff

private

Definition at line 99 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::pElectronFailed

private

Definition at line 100 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::pElectronPassed

private

Definition at line 102 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::phiWidthSuperClusters_

private

Definition at line 91 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::phiWidthSuperClustersEl_

private

Definition at line 96 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::ptDiff

private

Definition at line 98 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::ptElectronFailed

private

Definition at line 101 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::ptElectronPassed

private

Definition at line 103 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siElectronCollection_

private

Definition at line 197 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siElectronProducer_

private

Definition at line 196 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siHitProducer_

private

Definition at line 200 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siMatchedHitCollection_

private

Definition at line 203 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siRphiHitCollection_

private

Definition at line 201 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siStereoHitCollection_

private

Definition at line 202 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::sizeSuperClusters_

private

Definition at line 89 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::sizeSuperClustersEl_

private

Definition at line 94 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::sizeSuperClustersFailed

private

Definition at line 105 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::sizeSuperClustersPassed

private

Definition at line 104 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), endJob(), and SiStripElectronAnalyzer().

int SiStripElectronAnalyzer::StereoDetector_[1000]

private

Definition at line 141 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitCorr_[1000]

private

Definition at line 135 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitNoise_[1000]

private

Definition at line 138 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::StereoHitPhi_[1000]

private

Definition at line 129 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::StereoHitR_[1000]

private

Definition at line 128 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::StereoHitSignal_[1000]

private

Definition at line 137 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::StereoHitSigX_[1000]

private

Definition at line 133 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitSigY_[1000]

private

Definition at line 134 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitTheta_[1000]

private

Definition at line 130 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::StereoHitWidth_[1000]

private

Definition at line 139 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::StereoHitX_[1000]

private

Definition at line 124 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::StereoHitY_[1000]

private

Definition at line 125 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::StereoHitZ_[1000]

private

Definition at line 126 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

int SiStripElectronAnalyzer::StereoLayer_[1000]

private

Definition at line 142 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::string SiStripElectronAnalyzer::superClusterCollection_

private

Definition at line 191 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::superClusterProducer_

private

Definition at line 190 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), and SiStripElectronAnalyzer().

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiStripElectronAnalyzer::topoToken_

private

Definition at line 79 of file SiStripElectronAnalyzer.cc.

Referenced by analyze().

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiStripElectronAnalyzer::trackerToken_

private

Definition at line 80 of file SiStripElectronAnalyzer.cc.

Referenced by analyze().

float SiStripElectronAnalyzer::XShower_[1000]

private

Definition at line 119 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::YShower_[1000]

private

Definition at line 120 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().

float SiStripElectronAnalyzer::ZShower_[1000]

private

Definition at line 121 of file SiStripElectronAnalyzer.cc.

Referenced by analyze(), beginJob(), and initNtuple().