SiStripElectronAnalyzer Class Reference

#include <RecoEgamma/Examples/interface/SiStripElectronAnalyzer.h>

Inheritance diagram for SiStripElectronAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob (void)
virtual void	initNtuple (void)
	SiStripElectronAnalyzer (const edm::ParameterSet &)
	~SiStripElectronAnalyzer ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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_
float	XShower_ [1000]
float	YShower_ [1000]
float	ZShower_ [1000]

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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 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)

Detailed Description

Description: <one line="" class="" summary>="">

Usage: <usage>

Definition at line 46 of file SiStripElectronAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 67 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_.

{
  //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

(

)

Definition at line 158 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  )

virtual

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::EDAnalyzer.

Definition at line 325 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::false, edm::EventSetup::get(), edm::Event::getByLabel(), i, 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_, Signal, siHitProducer_, siMatchedHitCollection_, siRphiHitCollection_, siStereoHitCollection_, sizeSuperClusters_, sizeSuperClustersEl_, sizeSuperClustersFailed, sizeSuperClustersPassed, mathSSE::sqrt(), StereoDetector_, StereoHitCorr_, StereoHitNoise_, StereoHitPhi_, StereoHitR_, StereoHitSignal_, StereoHitSigX_, StereoHitSigY_, StereoHitTheta_, StereoHitWidth_, StereoHitX_, StereoHitY_, StereoHitZ_, StereoLayer_, superClusterCollection_, superClusterProducer_, PV3DBase< T, PVType, FrameType >::theta(), StripSubdetector::TIB, StripSubdetector::TOB, funct::true, 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.get<TrackerTopologyRcd>().get(tTopo);


  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.get<TrackerDigiGeometryRecord>().get(trackerHandle);

  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  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 184 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  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1255 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 265 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 58 of file SiStripElectronAnalyzer.h.

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 180 of file SiStripElectronAnalyzer.h.

std::string SiStripElectronAnalyzer::basicClusterProducer_

private

Definition at line 179 of file SiStripElectronAnalyzer.h.

std::string SiStripElectronAnalyzer::eBRecHitCollection_

private

Definition at line 182 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::eBRecHitProducer_

private

Definition at line 181 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::electronCollection_

private

Definition at line 186 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::electronProducer_

private

Definition at line 185 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::emaxSuperClusters_

private

Definition at line 75 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::emaxSuperClustersEl_

private

Definition at line 80 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::energySuperClusters_

private

Definition at line 73 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::energySuperClustersEl_

private

Definition at line 78 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::energySuperClustersFailed

private

Definition at line 92 of file SiStripElectronAnalyzer.h.

Referenced by SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::energySuperClustersPassed

private

Definition at line 91 of file SiStripElectronAnalyzer.h.

Referenced by SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::eOverPFailed

private

Definition at line 93 of file SiStripElectronAnalyzer.h.

Referenced by SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::eOverPPassed

private

Definition at line 94 of file SiStripElectronAnalyzer.h.

Referenced by SiStripElectronAnalyzer().

float SiStripElectronAnalyzer::EShower_[1000]

private

Definition at line 104 of file SiStripElectronAnalyzer.h.

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

TFile* SiStripElectronAnalyzer::file_

private

Definition at line 69 of file SiStripElectronAnalyzer.h.

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

std::string SiStripElectronAnalyzer::fileName_

private

Definition at line 66 of file SiStripElectronAnalyzer.h.

Referenced by SiStripElectronAnalyzer().

int SiStripElectronAnalyzer::MatchedDetector_[1000]

private

Definition at line 171 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitCorr_[1000]

private

Definition at line 165 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitNoise_[1000]

private

Definition at line 168 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MatchedHitPhi_[1000]

private

Definition at line 159 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MatchedHitR_[1000]

private

Definition at line 158 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MatchedHitSignal_[1000]

private

Definition at line 167 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MatchedHitSigX_[1000]

private

Definition at line 163 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitSigY_[1000]

private

Definition at line 164 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MatchedHitTheta_[1000]

private

Definition at line 160 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::MatchedHitWidth_[1000]

private

Definition at line 169 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MatchedHitX_[1000]

private

Definition at line 154 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MatchedHitY_[1000]

private

Definition at line 155 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MatchedHitZ_[1000]

private

Definition at line 156 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::MatchedLayer_[1000]

private

Definition at line 172 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

std::string SiStripElectronAnalyzer::mctruthCollection_

private

Definition at line 176 of file SiStripElectronAnalyzer.h.

Referenced by SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::mctruthProducer_

private

Definition at line 175 of file SiStripElectronAnalyzer.h.

Referenced by SiStripElectronAnalyzer().

int SiStripElectronAnalyzer::MonoDetector_[1000]

private

Definition at line 149 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitCorr_[1000]

private

Definition at line 143 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitNoise_[1000]

private

Definition at line 146 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MonoHitPhi_[1000]

private

Definition at line 137 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MonoHitR_[1000]

private

Definition at line 136 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MonoHitSignal_[1000]

private

Definition at line 145 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MonoHitSigX_[1000]

private

Definition at line 141 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitSigY_[1000]

private

Definition at line 142 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::MonoHitTheta_[1000]

private

Definition at line 138 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::MonoHitWidth_[1000]

private

Definition at line 147 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MonoHitX_[1000]

private

Definition at line 132 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MonoHitY_[1000]

private

Definition at line 133 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::MonoHitZ_[1000]

private

Definition at line 134 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::MonoLayer_[1000]

private

Definition at line 150 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

TTree* SiStripElectronAnalyzer::myTree_

private

Definition at line 101 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::NMatchedHits_

private

Definition at line 153 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::NMonoHits_

private

Definition at line 131 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::NShowers_

private

Definition at line 103 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::NStereoHits_

private

Definition at line 109 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::numCand_

private

Definition at line 70 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::numElectrons_

private

Definition at line 71 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::numSiMatchedHits_

private

Definition at line 99 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::numSiMonoHits_

private

Definition at line 98 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::numSiStereoHits_

private

Definition at line 97 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::numSuperClusters_

private

Definition at line 72 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::pDiff

private

Definition at line 84 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::pElectronFailed

private

Definition at line 85 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::pElectronPassed

private

Definition at line 87 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::phiWidthSuperClusters_

private

Definition at line 76 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::phiWidthSuperClustersEl_

private

Definition at line 81 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::ptDiff

private

Definition at line 83 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::ptElectronFailed

private

Definition at line 86 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::ptElectronPassed

private

Definition at line 88 of file SiStripElectronAnalyzer.h.

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

std::string SiStripElectronAnalyzer::siElectronCollection_

private

Definition at line 184 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siElectronProducer_

private

Definition at line 183 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siHitProducer_

private

Definition at line 187 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siMatchedHitCollection_

private

Definition at line 190 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siRphiHitCollection_

private

Definition at line 188 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::siStereoHitCollection_

private

Definition at line 189 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

TH1F* SiStripElectronAnalyzer::sizeSuperClusters_

private

Definition at line 74 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::sizeSuperClustersEl_

private

Definition at line 79 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::sizeSuperClustersFailed

private

Definition at line 90 of file SiStripElectronAnalyzer.h.

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

TH1F* SiStripElectronAnalyzer::sizeSuperClustersPassed

private

Definition at line 89 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::StereoDetector_[1000]

private

Definition at line 127 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitCorr_[1000]

private

Definition at line 121 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitNoise_[1000]

private

Definition at line 124 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::StereoHitPhi_[1000]

private

Definition at line 115 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::StereoHitR_[1000]

private

Definition at line 114 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::StereoHitSignal_[1000]

private

Definition at line 123 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::StereoHitSigX_[1000]

private

Definition at line 119 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitSigY_[1000]

private

Definition at line 120 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

float SiStripElectronAnalyzer::StereoHitTheta_[1000]

private

Definition at line 116 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::StereoHitWidth_[1000]

private

Definition at line 125 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::StereoHitX_[1000]

private

Definition at line 110 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::StereoHitY_[1000]

private

Definition at line 111 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::StereoHitZ_[1000]

private

Definition at line 112 of file SiStripElectronAnalyzer.h.

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

int SiStripElectronAnalyzer::StereoLayer_[1000]

private

Definition at line 128 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

std::string SiStripElectronAnalyzer::superClusterCollection_

private

Definition at line 178 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

std::string SiStripElectronAnalyzer::superClusterProducer_

private

Definition at line 177 of file SiStripElectronAnalyzer.h.

Referenced by analyze(), and SiStripElectronAnalyzer().

float SiStripElectronAnalyzer::XShower_[1000]

private

Definition at line 105 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::YShower_[1000]

private

Definition at line 106 of file SiStripElectronAnalyzer.h.

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

float SiStripElectronAnalyzer::ZShower_[1000]

private

Definition at line 107 of file SiStripElectronAnalyzer.h.

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