d1/dc2/SiStripElectronAnalyzer_8h_source.html

#ifndef RecoEgamma_Examples_SiStripElectronAnalyzer_h

#define RecoEgamma_Examples_SiStripElectronAnalyzer_h

// -*- C++ -*-

//

// Package:     RecoEgamma/Examples

// Class  :     SiStripElectronAnalyzer

//

//

// Original Author:

//         Created:  Fri May 26 16:52:45 EDT 2006

//


// system include files

#include <memory>

#include <map>

#include <cmath>


// user include files


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

#include "FWCore/Framework/interface/one/EDAnalyzer.h"


#include "TFile.h"

#include "TH1F.h"


#include "TNtuple.h"

#include "TTree.h"

#include "TBranch.h"


// forward declarations


#define myMaxHits 1000


//

// class decleration

//


class SiStripElectronAnalyzer : public edm::one::EDAnalyzer<> {

public:

  explicit SiStripElectronAnalyzer(const edm::ParameterSet&);

  ~SiStripElectronAnalyzer() override;


  void analyze(const edm::Event&, const edm::EventSetup&) override;

  void beginJob() override;

  virtual void initNtuple(void);

  void endJob(void) override;


private:

  double unwrapPhi(double phi) const {

    while (phi > M_PI) {

      phi -= 2. * M_PI;

    }

    while (phi < -M_PI) {

      phi += 2. * M_PI;

    }

    return phi;

  }


  // ----------member data ---------------------------

  std::string fileName_;


  TFile* file_;

  TH1F* numCand_;

  TH1F* numElectrons_;

  TH1F* numSuperClusters_;

  TH1F* energySuperClusters_;

  TH1F* sizeSuperClusters_;

  TH1F* emaxSuperClusters_;

  TH1F* phiWidthSuperClusters_;


  TH1F* energySuperClustersEl_;

  TH1F* sizeSuperClustersEl_;

  TH1F* emaxSuperClustersEl_;

  TH1F* phiWidthSuperClustersEl_;


  TH1F* ptDiff;

  TH1F* pDiff;

  TH1F* pElectronFailed;

  TH1F* ptElectronFailed;

  TH1F* pElectronPassed;

  TH1F* ptElectronPassed;

  TH1F* sizeSuperClustersPassed;

  TH1F* sizeSuperClustersFailed;

  TH1F* energySuperClustersPassed;

  TH1F* energySuperClustersFailed;

  TH1F* eOverPFailed;

  TH1F* eOverPPassed;


  TH1F* numSiStereoHits_;

  TH1F* numSiMonoHits_;

  TH1F* numSiMatchedHits_;


  TTree* myTree_;


  int NShowers_;

  float EShower_[myMaxHits];

  float XShower_[myMaxHits];

  float YShower_[myMaxHits];

  float ZShower_[myMaxHits];


  int NStereoHits_;

  float StereoHitX_[myMaxHits];

  float StereoHitY_[myMaxHits];

  float StereoHitZ_[myMaxHits];


  float StereoHitR_[myMaxHits];

  float StereoHitPhi_[myMaxHits];

  float StereoHitTheta_[myMaxHits];


  // errors in local coords

  float StereoHitSigX_[myMaxHits];

  float StereoHitSigY_[myMaxHits];

  float StereoHitCorr_[myMaxHits];


  float StereoHitSignal_[myMaxHits];

  float StereoHitNoise_[myMaxHits];

  int StereoHitWidth_[myMaxHits];


  int StereoDetector_[myMaxHits];

  int StereoLayer_[myMaxHits];


  // mono corresponds to "rphi" only hits

  int NMonoHits_;

  float MonoHitX_[myMaxHits];

  float MonoHitY_[myMaxHits];

  float MonoHitZ_[myMaxHits];


  float MonoHitR_[myMaxHits];

  float MonoHitPhi_[myMaxHits];

  float MonoHitTheta_[myMaxHits];


  // errors in local coords

  float MonoHitSigX_[myMaxHits];

  float MonoHitSigY_[myMaxHits];

  float MonoHitCorr_[myMaxHits];


  float MonoHitSignal_[myMaxHits];

  float MonoHitNoise_[myMaxHits];

  int MonoHitWidth_[myMaxHits];


  int MonoDetector_[myMaxHits];

  int MonoLayer_[myMaxHits];


  // matched hits

  int NMatchedHits_;

  float MatchedHitX_[myMaxHits];

  float MatchedHitY_[myMaxHits];

  float MatchedHitZ_[myMaxHits];


  float MatchedHitR_[myMaxHits];

  float MatchedHitPhi_[myMaxHits];

  float MatchedHitTheta_[myMaxHits];


  // errors in local coords

  float MatchedHitSigX_[myMaxHits];

  float MatchedHitSigY_[myMaxHits];

  float MatchedHitCorr_[myMaxHits];


  float MatchedHitSignal_[myMaxHits];

  float MatchedHitNoise_[myMaxHits];

  int MatchedHitWidth_[myMaxHits];


  int MatchedDetector_[myMaxHits];

  int MatchedLayer_[myMaxHits];


  std::string mctruthProducer_;

  std::string mctruthCollection_;

  std::string superClusterProducer_;

  std::string superClusterCollection_;

  std::string basicClusterProducer_;

  std::string basicClusterCollection_;

  std::string eBRecHitProducer_;

  std::string eBRecHitCollection_;

  std::string siElectronProducer_;

  std::string siElectronCollection_;

  std::string electronProducer_;

  std::string electronCollection_;

  std::string siHitProducer_;

  std::string siRphiHitCollection_;

  std::string siStereoHitCollection_;

  std::string siMatchedHitCollection_;

};


#endif