#include <ZeePlots.h>

Inheritance diagram for ZeePlots:

Public Member Functions
void	bookEleClassesPlots ()
void	bookEleHistograms ()
void	bookEleMCHistograms ()
void	bookHLTHistograms ()
void	bookZHistograms ()
void	bookZMCHistograms ()
void	fillEleClassesPlots (calib::CalibElectron *)
void	fillEleInfo (const reco::GsfElectronCollection *)
void	fillEleMCInfo (const HepMC::GenEvent *)
void	fillHLTInfo (edm::Handle< edm::TriggerResults >)
void	fillZInfo (std::pair< calib::CalibElectron , calib::CalibElectron > myZeeCandidate)
void	fillZMCInfo (const HepMC::GenEvent *)
void	openFile ()
void	writeEleClassesPlots ()
void	writeEleHistograms ()
void	writeHLTHistograms ()
void	writeMCEleHistograms ()
void	writeMCZHistograms ()
void	writeZHistograms ()
	ZeePlots (const char *)
	ZeePlots (const edm::ParameterSet &)
	~ZeePlots ()
	~ZeePlots ()
Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
bool	CheckCut (const pat::Electron *wenu, int i)
bool	CheckCut2 (const pat::Electron *wenu, int i)
bool	CheckCut2Inv (const pat::Electron *wenu, int i)
bool	CheckCutInv (const pat::Electron *wenu, int i)
bool	CheckCuts (const pat::Electron *ele)
bool	CheckCuts2 (const pat::Electron *ele)
bool	CheckCuts2Inverse (const pat::Electron *ele)
bool	CheckCuts2NminusOne (const pat::Electron *ele, int jj)
bool	CheckCutsInverse (const pat::Electron *ele)
bool	CheckCutsNminusOne (const pat::Electron *ele, int jj)
virtual void	endJob ()
double	ReturnCandVar (const pat::Electron *ele, int i)
Private Attributes
std::vector< Double_t >	CutVars2_
std::vector< Double_t >	CutVars_
Double_t	deta2_EB_
Double_t	deta2_EE_
Double_t	deta_EB_
bool	deta_EB_inv
Double_t	deta_EE_
bool	deta_EE_inv
Double_t	dphi2_EB_
Double_t	dphi2_EE_
Double_t	dphi_EB_
bool	dphi_EB_inv
Double_t	dphi_EE_
bool	dphi_EE_inv
Double_t	ecalIso2_EB_
Double_t	ecalIso2_EE_
Double_t	ecalIso_EB_
bool	ecalIso_EB_inv
Double_t	ecalIso_EE_
bool	ecalIso_EE_inv
TFile *	file_
const char *	fileName_
TH1F *	h1_FiredTriggers_
TH1F *	h1_gen_ZEta_
TH1F *	h1_gen_ZMass_
TH1F *	h1_gen_ZPhi_
TH1F *	h1_gen_ZPt_
TH1F *	h1_gen_ZRapidity_
TH1F *	h1_HLT1Electron_FiredEvents_
TH1F *	h1_HLT1Electron_HLT2Electron_FiredEvents_
TH1F *	h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_
TH1F *	h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents_
TH1F *	h1_HLT2Electron_FiredEvents_
TH1F *	h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_
TH1F *	h1_HLT2ElectronRelaxed_FiredEvents_
TH1F *	h1_HLTVisitedEvents_
TH1F *	h1_mcEle_Energy_
TH1F *	h1_mcEleEta_
TH1F *	h1_mcElePhi_
TH1F *	h1_mcElePt_
TH1F *	h1_nEleReco_
TH1F *	h1_occupancyVsEtaCrack_
TH1F *	h1_occupancyVsEtaGold_
TH1F *	h1_occupancyVsEtaShower_
TH1F *	h1_occupancyVsEtaSilver_
TH1F *	h1_reco_ZEta_
TH1F *	h1_reco_ZPhi_
TH1F *	h1_reco_ZPt_
TH1F *	h1_reco_ZRapidity_
TH1F *	h1_reco_ZTheta_
TH1F *	h1_recoEleEnergy_
TH1F *	h1_recoEleEta_
TH1F *	h1_recoElePhi_
TH1F *	h1_recoElePt_
TH1F *	h_e_ETA
TH1F *	h_e_PHI
TH1F *	h_e_PT
TH1F *	h_EB_deta
TH1F *	h_EB_dphi
TH1F *	h_EB_ecaliso
TH1F *	h_EB_hcaliso
TH1F *	h_EB_HoE
TH1F *	h_EB_sIetaIeta
TH1F *	h_EB_trkiso
TH1F *	h_EE_deta
TH1F *	h_EE_dphi
TH1F *	h_EE_ecaliso
TH1F *	h_EE_hcaliso
TH1F *	h_EE_HoE
TH1F *	h_EE_sIetaIeta
TH1F *	h_EE_trkiso
TH1F *	h_mee
TH1F *	h_mee_EBEB
TH1F *	h_mee_EBEE
TH1F *	h_mee_EEEE
TH1F *	h_Zcand_PT
TH1F *	h_Zcand_Y
Double_t	hcalIso2_EB_
Double_t	hcalIso2_EE_
Double_t	hcalIso_EB_
bool	hcalIso_EB_inv
Double_t	hcalIso_EE_
bool	hcalIso_EE_inv
TFile *	histofile
Double_t	hoe2_EB_
Double_t	hoe2_EE_
Double_t	hoe_EB_
bool	hoe_EB_inv
Double_t	hoe_EE_
bool	hoe_EE_inv
std::vector< Double_t >	InvVars_
int	nBarrelVars_
std::string	outputFile_
Double_t	sihih2_EB_
Double_t	sihih2_EE_
Double_t	sihih_EB_
bool	sihih_EB_inv
Double_t	sihih_EE_
bool	sihih_EE_inv
Double_t	trackIso2_EB_
Double_t	trackIso2_EE_
Double_t	trackIso_EB_
bool	trackIso_EB_inv
Double_t	trackIso_EE_
bool	trackIso_EE_inv
bool	useDifferentSecondLegSelection_
Double_t	userIso2_EB_
Double_t	userIso2_EE_
Double_t	userIso_EB_
bool	userIso_EB_inv
Double_t	userIso_EE_
bool	userIso_EE_inv
edm::InputTag	zeeCollectionTag_

Detailed Description

Definition at line 53 of file ZeePlots.h.

Constructor & Destructor Documentation

ZeePlots::ZeePlots ( const char * fileName )

Definition at line 35 of file ZeePlots.cc.

References file_, convertXMLtoSQLite_cfg::fileName, and fileName_.

{

  fileName_ = fileName;
  file_ = new TFile(fileName_, "RECREATE");
}

ZeePlots::~ZeePlots ( )

Definition at line 43 of file ZeePlots.cc.

References file_.

{
  
  file_->Close();
  
  delete file_;

}

ZeePlots::ZeePlots ( const edm::ParameterSet & iConfig ) [explicit]

Definition at line 21 of file ZeePlots.cc.

References gather_cfg::cout, deta2_EB_, deta2_EE_, deta_EB_, deta_EB_inv, deta_EE_, deta_EE_inv, dphi2_EB_, dphi2_EE_, dphi_EB_, dphi_EB_inv, dphi_EE_, dphi_EE_inv, ecalIso2_EB_, ecalIso2_EE_, ecalIso_EB_, ecalIso_EB_inv, ecalIso_EE_, ecalIso_EE_inv, edm::ParameterSet::getUntrackedParameter(), hcalIso2_EB_, hcalIso2_EE_, hcalIso_EB_, hcalIso_EB_inv, hcalIso_EE_, hcalIso_EE_inv, hoe2_EB_, hoe2_EE_, hoe_EB_, hoe_EB_inv, hoe_EE_, hoe_EE_inv, outputFile_, sihih2_EB_, sihih2_EE_, sihih_EB_, sihih_EB_inv, sihih_EE_, sihih_EE_inv, trackIso2_EB_, trackIso2_EE_, trackIso_EB_, trackIso_EB_inv, trackIso_EE_, trackIso_EE_inv, useDifferentSecondLegSelection_, userIso2_EB_, userIso2_EE_, userIso_EB_, userIso_EB_inv, userIso_EE_, userIso_EE_inv, and zeeCollectionTag_.

{
//                   I N P U T      P A R A M E T E R S
//
//  ZEE COLLECTION   //////////////////////////////////////////////////////
//
  
  zeeCollectionTag_ = iConfig.getUntrackedParameter<edm::InputTag>
    ("zeeCollectionTag");
  //
  // code parameters
  //
  std::string outputFile_D = "histos.root";
  outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile",
                                                           outputFile_D);
  //
  //
  // the selection cuts:
  trackIso_EB_ = iConfig.getUntrackedParameter<Double_t>("trackIso_EB");
  ecalIso_EB_ = iConfig.getUntrackedParameter<Double_t>("ecalIso_EB");
  hcalIso_EB_ = iConfig.getUntrackedParameter<Double_t>("hcalIso_EB");
  //
  trackIso_EE_ = iConfig.getUntrackedParameter<Double_t>("trackIso_EE");
  ecalIso_EE_ = iConfig.getUntrackedParameter<Double_t>("ecalIso_EE");
  hcalIso_EE_ = iConfig.getUntrackedParameter<Double_t>("hcalIso_EE");
  //
  sihih_EB_ = iConfig.getUntrackedParameter<Double_t>("sihih_EB");
  dphi_EB_ = iConfig.getUntrackedParameter<Double_t>("dphi_EB");
  deta_EB_ = iConfig.getUntrackedParameter<Double_t>("deta_EB");
  hoe_EB_ = iConfig.getUntrackedParameter<Double_t>("hoe_EB");
  userIso_EB_ = iConfig.getUntrackedParameter<Double_t>("userIso_EB",1000.);
  //
  sihih_EE_ = iConfig.getUntrackedParameter<Double_t>("sihih_EE");
  dphi_EE_ = iConfig.getUntrackedParameter<Double_t>("dphi_EE");
  deta_EE_ = iConfig.getUntrackedParameter<Double_t>("deta_EE");
  hoe_EE_ = iConfig.getUntrackedParameter<Double_t>("hoe_EE");
  userIso_EE_ = iConfig.getUntrackedParameter<Double_t>("userIso_EE",1000.);
  //
  trackIso_EB_inv = iConfig.getUntrackedParameter<Bool_t>("trackIso_EB_inv", 
                                                            false);
  ecalIso_EB_inv = iConfig.getUntrackedParameter<Bool_t>("ecalIso_EB_inv",
                                                           false);
  hcalIso_EB_inv = iConfig.getUntrackedParameter<Bool_t>("hcalIso_EB_inv",
                                                           false);
  //
  trackIso_EE_inv = iConfig.getUntrackedParameter<Bool_t>("trackIso_EE_inv",
                                                            false);
  ecalIso_EE_inv = iConfig.getUntrackedParameter<Bool_t>("ecalIso_EE_inv",
                                                           false);
  hcalIso_EE_inv = iConfig.getUntrackedParameter<Bool_t>("hcalIso_EE_inv",
                                                           false);
  //
  sihih_EB_inv = iConfig.getUntrackedParameter<Bool_t>("sihih_EB_inv", false);
  dphi_EB_inv = iConfig.getUntrackedParameter<Bool_t>("dphi_EB_inv",false);
  deta_EB_inv = iConfig.getUntrackedParameter<Bool_t>("deta_EB_inv",false);
  hoe_EB_inv = iConfig.getUntrackedParameter<Bool_t>("hoe_EB_inv",false);
  userIso_EB_inv=iConfig.getUntrackedParameter<Bool_t>("userIso_EB_inv",false);
  //
  sihih_EE_inv = iConfig.getUntrackedParameter<Bool_t>("sihih_EE_inv", false);
  dphi_EE_inv = iConfig.getUntrackedParameter<Bool_t>("dphi_EE_inv", false);
  deta_EE_inv = iConfig.getUntrackedParameter<Bool_t>("deta_EE_inv",false);
  hoe_EE_inv = iConfig.getUntrackedParameter<Bool_t>("hoe_EE_inv",false);
  userIso_EE_inv=iConfig.getUntrackedParameter<Bool_t>("userIso_EE_inv",false);
  useDifferentSecondLegSelection_ = iConfig.getUntrackedParameter<Bool_t>
    ("useDifferentSecondLegSelection",false);
  if (useDifferentSecondLegSelection_) {
    std::cout << "ZeePlots: WARNING: you have chosen to use a different "
              << " selection for the 2nd leg of the Z" << std::endl;
    trackIso2_EB_ = iConfig.getUntrackedParameter<Double_t>("trackIso2_EB");
    ecalIso2_EB_ = iConfig.getUntrackedParameter<Double_t>("ecalIso2_EB");
    hcalIso2_EB_ = iConfig.getUntrackedParameter<Double_t>("hcalIso2_EB");
    //
    trackIso2_EE_ = iConfig.getUntrackedParameter<Double_t>("trackIso2_EE");
    ecalIso2_EE_ = iConfig.getUntrackedParameter<Double_t>("ecalIso2_EE");
    hcalIso2_EE_ = iConfig.getUntrackedParameter<Double_t>("hcalIso2_EE");
    //
    sihih2_EB_ = iConfig.getUntrackedParameter<Double_t>("sihih2_EB");
    dphi2_EB_ = iConfig.getUntrackedParameter<Double_t>("dphi2_EB");
    deta2_EB_ = iConfig.getUntrackedParameter<Double_t>("deta2_EB");
    hoe2_EB_ = iConfig.getUntrackedParameter<Double_t>("hoe2_EB");
    userIso2_EB_=iConfig.getUntrackedParameter<Double_t>("userIso2_EB", 1000.);
    //
    sihih2_EE_ = iConfig.getUntrackedParameter<Double_t>("sihih2_EE");
    dphi2_EE_ = iConfig.getUntrackedParameter<Double_t>("dphi2_EE");
    deta2_EE_ = iConfig.getUntrackedParameter<Double_t>("deta2_EE");
    hoe2_EE_ = iConfig.getUntrackedParameter<Double_t>("hoe2_EE");
    userIso2_EE_=iConfig.getUntrackedParameter<Double_t>("userIso2_EE", 1000.);
  }
  else {
    trackIso2_EB_ = trackIso_EB_;
    ecalIso2_EB_ = ecalIso_EB_;
    hcalIso2_EB_ = hcalIso_EB_;
    //
    trackIso2_EE_ = trackIso_EE_;
    ecalIso2_EE_ = ecalIso_EE_;
    hcalIso2_EE_ = hcalIso_EE_;
    //
    sihih2_EB_ = sihih_EB_;
    dphi2_EB_ = dphi_EB_;
    deta2_EB_ = deta_EB_;
    hoe2_EB_ = hoe_EB_;
    userIso2_EB_ = userIso_EB_;
    //
    sihih2_EE_ = sihih_EE_;
    dphi2_EE_ = dphi_EE_;
    deta2_EE_ = deta_EE_;
    hoe2_EE_ = hoe_EE_;
    userIso2_EE_ = userIso_EE_;
  }

}

ZeePlots::~ZeePlots ( )

Member Function Documentation

void ZeePlots::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	es
	)		`[private, virtual]`

Implements edm::EDAnalyzer.

Definition at line 154 of file ZeePlots.cc.

References CheckCuts(), CheckCuts2(), gather_cfg::cout, pat::EcalIso, edm::Event::getByLabel(), h_e_ETA, h_e_PHI, h_e_PT, h_EB_deta, h_EB_dphi, h_EB_ecaliso, h_EB_hcaliso, h_EB_HoE, h_EB_sIetaIeta, h_EB_trkiso, h_EE_deta, h_EE_dphi, h_EE_ecaliso, h_EE_hcaliso, h_EE_HoE, h_EE_sIetaIeta, h_EE_trkiso, h_mee, h_mee_EBEB, h_mee_EBEE, h_mee_EEEE, h_Zcand_PT, h_Zcand_Y, pat::HcalIso, pat::TrackIso, Gflash::Z, and zeeCollectionTag_.

{
  using namespace std;
  //  cout << "In analyzer now..." << endl;
  //
  //  Get the collections here
  //
  edm::Handle<pat::CompositeCandidateCollection> ZeeCands;
  iEvent.getByLabel(zeeCollectionTag_, ZeeCands);

  if (not ZeeCands.isValid()) {
    cout << "Warning: no zee candidates in this event..." << endl;
    return;
  }
  //
  //
  const pat::CompositeCandidateCollection *zcands = ZeeCands.product();
  const pat::CompositeCandidateCollection::const_iterator 
    zeeIter = zcands->begin();
  const pat::CompositeCandidate zee = *zeeIter;
  //
  // get the parts of the composite candidate:
  const pat::Electron * myElec1=
    dynamic_cast<const pat::Electron*> (zee.daughter("electron1"));
  const pat::Electron * myElec2=
    dynamic_cast<const pat::Electron*> (zee.daughter("electron2"));
  // you can access MET like that if needed in some application
  //  const pat::MET * myMet=
  //  dynamic_cast<const pat::MET*> (zee.daughter("met"));

  //double met = myMet->et();


  TLorentzVector e1;
  TLorentzVector e2;

//  math::XYZVector p1 =   myElec1->trackMomentumAtVtx();
//  math::XYZVector p2 =   myElec2->trackMomentumAtVtx();
//  e1.SetPxPyPzE(p1.X(), p1.Y(), p1.Z(), myElec1->caloEnergy());
//  e2.SetPxPyPzE(p2.X(), p2.Y(), p2.Z(), myElec2->caloEnergy());

  // Use directly the et,eta,phi from pat::Electron; assume e mass = 0.0
  e1.SetPtEtaPhiM(myElec1->et(),myElec1->eta(),myElec1->phi(),0.0);
  e2.SetPtEtaPhiM(myElec2->et(),myElec2->eta(),myElec2->phi(),0.0);

 
  TLorentzVector Z = e1+e2;
  double mee = Z.M();
  // the selection plots:
  bool pass = CheckCuts(myElec1) && CheckCuts2(myElec2);
  //cout << "This event passes? " << pass << ", mee is: " << mee
  //   << " and the histo is filled." << endl;
  if (not pass) return;

  h_mee->Fill(mee);
  if(fabs(e1.Eta())<1.479 && fabs(e2.Eta())<1.479)h_mee_EBEB->Fill(mee);
  if(fabs(e1.Eta())<1.479 && fabs(e2.Eta())>1.479)h_mee_EBEE->Fill(mee);
  if(fabs(e1.Eta())>1.479 && fabs(e2.Eta())<1.479)h_mee_EBEE->Fill(mee);
  if(fabs(e1.Eta())>1.479 && fabs(e2.Eta())>1.479)h_mee_EEEE->Fill(mee);
  
  h_Zcand_PT->Fill(Z.Pt());
  h_Zcand_Y->Fill(Z.Rapidity());

  h_e_PT->Fill(e1.Pt()); h_e_PT->Fill(e2.Pt()); 
  h_e_ETA->Fill(e1.Eta()); h_e_ETA->Fill(e2.Eta()); 
  h_e_PHI->Fill(e1.Phi()); h_e_PHI->Fill(e2.Phi()); 

   if(fabs(myElec1->eta())<1.479){
      h_EB_trkiso->Fill( myElec1->userIsolation(pat::TrackIso) );
      h_EB_ecaliso->Fill( myElec1->userIsolation(pat::EcalIso) );
      h_EB_hcaliso->Fill( myElec1->userIsolation(pat::HcalIso) );
      h_EB_sIetaIeta->Fill( myElec1->scSigmaIEtaIEta() );
      h_EB_dphi->Fill( myElec1->deltaPhiSuperClusterTrackAtVtx() );
      h_EB_deta->Fill( myElec1->deltaEtaSuperClusterTrackAtVtx() );
      h_EB_HoE->Fill( myElec1->hadronicOverEm() );
    }
    else{
      h_EE_trkiso->Fill( myElec1->userIsolation(pat::TrackIso) );
      h_EE_ecaliso->Fill( myElec1->userIsolation(pat::EcalIso) );
      h_EE_hcaliso->Fill( myElec1->userIsolation(pat::HcalIso) );
      h_EE_sIetaIeta->Fill( myElec1->scSigmaIEtaIEta() );
      h_EE_dphi->Fill( myElec1->deltaPhiSuperClusterTrackAtVtx() );
      h_EE_deta->Fill( myElec1->deltaEtaSuperClusterTrackAtVtx() );
      h_EE_HoE->Fill( myElec1->hadronicOverEm() );
    }


 if(fabs(myElec2->eta())<1.479){
      h_EB_trkiso->Fill( myElec2->userIsolation(pat::TrackIso) );
      h_EB_ecaliso->Fill( myElec2->userIsolation(pat::EcalIso) );
      h_EB_hcaliso->Fill( myElec2->userIsolation(pat::HcalIso) );
      h_EB_sIetaIeta->Fill( myElec2->scSigmaIEtaIEta() );
      h_EB_dphi->Fill( myElec2->deltaPhiSuperClusterTrackAtVtx() );
      h_EB_deta->Fill( myElec2->deltaEtaSuperClusterTrackAtVtx() );
      h_EB_HoE->Fill( myElec2->hadronicOverEm() );
    }
    else{
      h_EE_trkiso->Fill( myElec2->userIsolation(pat::TrackIso) );
      h_EE_ecaliso->Fill( myElec2->userIsolation(pat::EcalIso) );
      h_EE_hcaliso->Fill( myElec2->userIsolation(pat::HcalIso) );
      h_EE_sIetaIeta->Fill( myElec2->scSigmaIEtaIEta() );
      h_EE_dphi->Fill( myElec2->deltaPhiSuperClusterTrackAtVtx() );
      h_EE_deta->Fill( myElec2->deltaEtaSuperClusterTrackAtVtx() );
      h_EE_HoE->Fill( myElec2->hadronicOverEm() );
    }




  //double scEta = myElec->superCluster()->eta();
  //double scPhi = myElec->superCluster()->phi();
  //double scEt = myElec->superCluster()->energy()/cosh(scEta);

}

void ZeePlots::beginJob ( void ) [private, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 400 of file ZeePlots.cc.

References CutVars2_, CutVars_, deta2_EB_, deta2_EE_, deta_EB_, deta_EB_inv, deta_EE_, deta_EE_inv, dphi2_EB_, dphi2_EE_, dphi_EB_, dphi_EB_inv, dphi_EE_, dphi_EE_inv, ecalIso2_EB_, ecalIso2_EE_, ecalIso_EB_, ecalIso_EB_inv, ecalIso_EE_, ecalIso_EE_inv, h_e_ETA, h_e_PHI, h_e_PT, h_EB_deta, h_EB_dphi, h_EB_ecaliso, h_EB_hcaliso, h_EB_HoE, h_EB_sIetaIeta, h_EB_trkiso, h_EE_deta, h_EE_dphi, h_EE_ecaliso, h_EE_hcaliso, h_EE_HoE, h_EE_sIetaIeta, h_EE_trkiso, h_mee, h_mee_EBEB, h_mee_EBEE, h_mee_EEEE, h_Zcand_PT, h_Zcand_Y, hcalIso2_EB_, hcalIso2_EE_, hcalIso_EB_, hcalIso_EB_inv, hcalIso_EE_, hcalIso_EE_inv, hoe2_EB_, hoe2_EE_, hoe_EB_, hoe_EB_inv, hoe_EE_, hoe_EE_inv, InvVars_, nBarrelVars_, sihih2_EB_, sihih2_EE_, sihih_EB_, sihih_EB_inv, sihih_EE_, sihih_EE_inv, trackIso2_EB_, trackIso2_EE_, trackIso_EB_, trackIso_EB_inv, trackIso_EE_, trackIso_EE_inv, userIso2_EB_, userIso2_EE_, userIso_EB_, userIso_EB_inv, userIso_EE_, and userIso_EE_inv.

{
  //std::cout << "In beginJob()" << std::endl;
  h_mee      = new TH1F("h_mee",      "h_mee",       200, 0, 200);
  h_mee_EBEB = new TH1F("h_mee_EBEB", "h_mee_EBEB", 200, 0, 200);
  h_mee_EBEE = new TH1F("h_mee_EBEE", "h_mee_EBEE", 200, 0, 200);
  h_mee_EEEE = new TH1F("h_mee_EEEE", "h_mee_EEEE", 200, 0, 200);

  h_Zcand_PT = new TH1F("h_Zcand_PT", "h_Zcand_PT", 200,  0, 100);
  h_Zcand_Y  = new TH1F("h_Zcand_Y",  "h_Zcand_Y" , 200, -5, 5);

  h_e_PT  = new TH1F("h_e_PT", "h_e_PT", 200,  0, 100);
  h_e_ETA = new TH1F("h_e_ETA","h_e_ETA",200, -3, 3);
  h_e_PHI = new TH1F("h_e_PHI","h_e_PHI",200, -4, 4);


  //VALIDATION PLOTS
  //EB
  h_EB_trkiso = new TH1F("h_EB_trkiso","h_EB_trkiso",200 , 0.0, 9.0);
  h_EB_ecaliso = new TH1F("h_EB_ecaliso","h_EB_ecaliso",200, 0.0 , 9.0);
  h_EB_hcaliso = new TH1F("h_EB_hcaliso","h_EB_hcaliso",200, 0.0 , 9.0);
  h_EB_sIetaIeta = new TH1F("h_EB_sIetaIeta","h_EB_sIetaIeta",200, 0.0 , 0.02 );
  h_EB_dphi = new TH1F("h_EB_dphi","h_EB_dphi",200, -0.03 , 0.03 );
  h_EB_deta = new TH1F("h_EB_deta","h_EB_deta",200, -0.01 , 0.01) ;
  h_EB_HoE = new TH1F("h_EB_HoE","h_EB_HoE",200, 0.0 , 0.2 );
  //EE
  h_EE_trkiso = new TH1F("h_EE_trkiso","h_EE_trkiso",200 , 0.0, 9.0);
  h_EE_ecaliso = new TH1F("h_EE_ecaliso","h_EE_ecaliso",200, 0.0 , 9.0);
  h_EE_hcaliso = new TH1F("h_EE_hcaliso","h_EE_hcaliso",200, 0.0 , 9.0);
  h_EE_sIetaIeta = new TH1F("h_EE_sIetaIeta","h_EE_sIetaIeta",200, 0.0 , 0.1 );
  h_EE_dphi = new TH1F("h_EE_dphi","h_EE_dphi",200, -0.03 , 0.03 );
  h_EE_deta = new TH1F("h_EE_deta","h_EE_deta",200, -0.01 , 0.01) ;
  h_EE_HoE = new TH1F("h_EE_HoE","h_EE_HoE",200, 0.0 , 0.2 );





  //
  // if you add some new variable change the nBarrelVars_ accordingly
  nBarrelVars_ = 8;
  //
  // Put EB variables together and EE variables together
  // number of barrel variables = number of endcap variable
  // if you don't want to use some variable put a very high cut
  CutVars_.push_back( trackIso_EB_ );//0
  CutVars_.push_back( ecalIso_EB_ ); //1
  CutVars_.push_back( hcalIso_EB_ ); //2
  CutVars_.push_back( sihih_EB_ );   //3
  CutVars_.push_back( dphi_EB_ );    //4
  CutVars_.push_back( deta_EB_ );    //5
  CutVars_.push_back( hoe_EB_ );     //6
  CutVars_.push_back( userIso_EB_ ); //7

  CutVars_.push_back( trackIso_EE_);//0
  CutVars_.push_back( ecalIso_EE_); //1
  CutVars_.push_back( hcalIso_EE_); //2
  CutVars_.push_back( sihih_EE_);   //3
  CutVars_.push_back( dphi_EE_);    //4
  CutVars_.push_back( deta_EE_);    //5
  CutVars_.push_back( hoe_EE_ );    //6 
  CutVars_.push_back( userIso_EE_ );//7 
  //
  // 2nd leg variables
  CutVars2_.push_back( trackIso2_EB_ );//0
  CutVars2_.push_back( ecalIso2_EB_ ); //1
  CutVars2_.push_back( hcalIso2_EB_ ); //2
  CutVars2_.push_back( sihih2_EB_ );   //3
  CutVars2_.push_back( dphi2_EB_ );    //4
  CutVars2_.push_back( deta2_EB_ );    //5
  CutVars2_.push_back( hoe2_EB_ );     //6
  CutVars2_.push_back( userIso2_EB_ ); //7

  CutVars2_.push_back( trackIso2_EE_);//0
  CutVars2_.push_back( ecalIso2_EE_); //1
  CutVars2_.push_back( hcalIso2_EE_); //2
  CutVars2_.push_back( sihih2_EE_);   //3
  CutVars2_.push_back( dphi2_EE_);    //4
  CutVars2_.push_back( deta2_EE_);    //5
  CutVars2_.push_back( hoe2_EE_ );    //6 
  CutVars2_.push_back( userIso2_EE_ );//7 
  //...........................................
  InvVars_.push_back( trackIso_EB_inv);//0
  InvVars_.push_back( ecalIso_EB_inv); //1
  InvVars_.push_back( hcalIso_EB_inv); //2
  InvVars_.push_back( sihih_EB_inv);   //3
  InvVars_.push_back( dphi_EB_inv);    //4
  InvVars_.push_back( deta_EB_inv);    //5
  InvVars_.push_back( hoe_EB_inv);     //6
  InvVars_.push_back( userIso_EB_inv); //7
  //
  InvVars_.push_back( trackIso_EE_inv);//0
  InvVars_.push_back( ecalIso_EE_inv); //1
  InvVars_.push_back( hcalIso_EE_inv); //2
  InvVars_.push_back( sihih_EE_inv);   //3
  InvVars_.push_back( dphi_EE_inv);    //4
  InvVars_.push_back( deta_EE_inv);    //5
  InvVars_.push_back( hoe_EE_inv);     //6
  InvVars_.push_back( userIso_EE_inv); //7
  //



}

void ZeePlots::bookEleClassesPlots ( )

Definition at line 418 of file ZeePlots.cc.

References file_, h1_occupancyVsEtaCrack_, h1_occupancyVsEtaGold_, h1_occupancyVsEtaShower_, and h1_occupancyVsEtaSilver_.

                                  {

  file_->cd();

  h1_occupancyVsEtaGold_ = new TH1F("occupancyVsEtaGold","occupancyVsEtaGold", 200, -4.,4.);
  h1_occupancyVsEtaGold_->SetYTitle("Electron statistics");
  h1_occupancyVsEtaGold_->SetXTitle("Eta channel");

  h1_occupancyVsEtaSilver_ = new TH1F("occupancyVsEtaSilver","occupancyVsEtaSilver", 200, -4.,4.);
  h1_occupancyVsEtaSilver_->SetYTitle("Electron statistics");
  h1_occupancyVsEtaSilver_->SetXTitle("Eta channel");

  h1_occupancyVsEtaShower_ = new TH1F("occupancyVsEtaShower","occupancyVsEtaShower", 200, -4.,4.);
  h1_occupancyVsEtaShower_->SetYTitle("Electron statistics");
  h1_occupancyVsEtaShower_->SetXTitle("Eta channel");

  h1_occupancyVsEtaCrack_ = new TH1F("occupancyVsEtaCrack","occupancyVsEtaCrack", 200, -4.,4.);
  h1_occupancyVsEtaCrack_->SetYTitle("Electron statistics");
  h1_occupancyVsEtaCrack_->SetXTitle("Eta channel");

}

void ZeePlots::bookEleHistograms ( )

Definition at line 231 of file ZeePlots.cc.

References file_, h1_nEleReco_, h1_recoEleEnergy_, h1_recoEleEta_, h1_recoElePhi_, and h1_recoElePt_.

Referenced by ZeeCalibration::bookHistograms().

                                {

  file_->cd();

  h1_nEleReco_ = new TH1F("h1_nEleReco", "h1_nEleReco", 20, 0, 20 );
  h1_nEleReco_->SetXTitle("Num. of reco electrons");
  
  h1_recoEleEnergy_ = new TH1F("recoEleEnergy","EleEnergy from SC",300,0.,300.);
  h1_recoEleEnergy_->SetXTitle("eleSCEnergy(GeV)");
  h1_recoEleEnergy_->SetYTitle("events");
  
  h1_recoElePt_ = new TH1F("recoElePt","p_{T} of reco electrons",300,0.,300.);
  h1_recoElePt_->SetXTitle("p_{T}(GeV/c)");
  h1_recoElePt_->SetYTitle("events");

  h1_recoEleEta_ = new TH1F("recoEleEta","Eta of reco electrons",100,-4.,4.);
  h1_recoEleEta_->SetXTitle("#eta");
  h1_recoEleEta_->SetYTitle("events");

   
  h1_recoElePhi_ = new TH1F("recoElePhi","Phi of reco electrons",100,-4.,4.);
  h1_recoElePhi_->SetXTitle("#phi");
  h1_recoElePhi_->SetYTitle("events");



}

void ZeePlots::bookEleMCHistograms ( )

Definition at line 186 of file ZeePlots.cc.

References file_, h1_mcEle_Energy_, h1_mcEleEta_, h1_mcElePhi_, and h1_mcElePt_.

Referenced by ZeeCalibration::bookHistograms().

                                  {

  file_->cd();

  h1_mcEle_Energy_ = new TH1F("mcEleEnergy","mc EleEnergy",300,0.,300.);
  h1_mcEle_Energy_->SetXTitle("E (GeV)");
  h1_mcEle_Energy_->SetYTitle("events");

  h1_mcElePt_ = new TH1F("mcElePt","p_{T} of MC electrons",300,0.,300.);
  h1_mcElePt_->SetXTitle("p_{T}(GeV/c)");
  h1_mcElePt_->SetYTitle("events");
  
  h1_mcEleEta_ = new TH1F("mcEleEta","Eta of MC electrons",100,-4.,4.);
  h1_mcEleEta_->SetXTitle("#eta");
  h1_mcEleEta_->SetYTitle("events");

  h1_mcElePhi_ = new TH1F("mcElePhi","Phi of MC electrons",100,-4.,4.);
  h1_mcElePhi_->SetXTitle("#phi");
  h1_mcElePhi_->SetYTitle("events");

}

void ZeePlots::bookHLTHistograms ( )

Definition at line 317 of file ZeePlots.cc.

References file_, h1_FiredTriggers_, h1_HLT1Electron_FiredEvents_, h1_HLT1Electron_HLT2Electron_FiredEvents_, h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_, h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents_, h1_HLT2Electron_FiredEvents_, h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_, h1_HLT2ElectronRelaxed_FiredEvents_, and h1_HLTVisitedEvents_.

                                {

  file_->cd();

  h1_FiredTriggers_= new TH1F("h1_FiredTriggers", "h1_FiredTriggers", 5,0,5);
 
  h1_HLTVisitedEvents_ = new TH1F("h1_HLTVisitedEvents", "h1_HLTVisitedEvents", 5,0,5);

  h1_HLT1Electron_FiredEvents_ =  new TH1F("h1_HLT1Electron_FiredEvents", "h1_HLT1Electron_FiredEvents", 5,0,5);
  h1_HLT2Electron_FiredEvents_ =  new TH1F("h1_HLT2Electron_FiredEvents", "h1_HLT2Electron_FiredEvents", 5,0,5);
  h1_HLT2ElectronRelaxed_FiredEvents_ =  new TH1F("h1_HLT2ElectronRelaxed_FiredEvents", "h1_HLT2ElectronRelaxed_FiredEvents", 5,0,5);

  h1_HLT1Electron_HLT2Electron_FiredEvents_ =  new TH1F("h1_HLT1Electron_HLT2Electron_FiredEvents", "h1_HLT1Electron_HLT2Electron_FiredEvents", 5,0,5);
  h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents_ =  new TH1F("h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents", "h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents", 5,0,5);
  h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_ =  new TH1F("h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents", "h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents", 5,0,5);

  h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_ =  new TH1F("h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents", "h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents", 5,0,5);


}

void ZeePlots::bookZHistograms ( )

Definition at line 89 of file ZeePlots.cc.

References file_, h1_reco_ZEta_, h1_reco_ZPhi_, h1_reco_ZPt_, h1_reco_ZRapidity_, and h1_reco_ZTheta_.

Referenced by ZeeCalibration::bookHistograms().

                              {

  file_ -> cd();

  h1_reco_ZEta_ = new TH1F("reco_ZEta","Eta of reco Z",200,-6.,6.);
  h1_reco_ZEta_->SetXTitle("#eta");
  h1_reco_ZEta_->SetYTitle("events");
  
  h1_reco_ZTheta_ = new TH1F("reco_ZTheta","Theta of reco Z",200, 0., 4.);
  h1_reco_ZTheta_->SetXTitle("#theta");
  h1_reco_ZTheta_->SetYTitle("events");
  
  h1_reco_ZRapidity_ = new TH1F("reco_ZRapidity","Rapidity of reco Z",200,-6.,6.);
  h1_reco_ZRapidity_->SetXTitle("Y");
  h1_reco_ZRapidity_->SetYTitle("events");
  
  h1_reco_ZPhi_ = new TH1F("reco_ZPhi","Phi of reco Z",100,-4.,4.);
  h1_reco_ZPhi_->SetXTitle("#phi");
  h1_reco_ZPhi_->SetYTitle("events");
  
  h1_reco_ZPt_ = new TH1F("reco_ZPt","Pt of reco Z",200,0.,100.);
  h1_reco_ZPt_->SetXTitle("p_{T} (GeV/c)");
  h1_reco_ZPt_->SetYTitle("events");
  
  
}

void ZeePlots::bookZMCHistograms ( )

Definition at line 62 of file ZeePlots.cc.

References file_, h1_gen_ZEta_, h1_gen_ZMass_, h1_gen_ZPhi_, h1_gen_ZPt_, and h1_gen_ZRapidity_.

Referenced by ZeeCalibration::bookHistograms().

                                {

  file_ -> cd();

  h1_gen_ZMass_ = new TH1F("gen_ZMass","Generated Z mass",200,0.,150.);
  h1_gen_ZMass_->SetXTitle("gen_ZMass (GeV)");
  h1_gen_ZMass_->SetYTitle("events");

  h1_gen_ZEta_ = new TH1F("gen_ZEta","Eta of gen Z",200,-6.,6.);
  h1_gen_ZEta_->SetXTitle("#eta");
  h1_gen_ZEta_->SetYTitle("events");

  h1_gen_ZPhi_ = new TH1F("gen_ZPhi","Phi of gen Z",200,-4.,4.);
  h1_gen_ZPhi_->SetXTitle("#phi");
  h1_gen_ZPhi_->SetYTitle("events");

  h1_gen_ZRapidity_ = new TH1F("gen_ZRapidity","Rapidity of gen Z",200,-6.,6.);
  h1_gen_ZRapidity_->SetXTitle("Y");
  h1_gen_ZRapidity_->SetYTitle("events");

  h1_gen_ZPt_ = new TH1F("gen_ZPt","Pt of gen Z",200, 0.,100.);
  h1_gen_ZPt_->SetXTitle("p_{T} (GeV/c)");
  h1_gen_ZPt_->SetYTitle("events");


}

bool ZeePlots::CheckCut	(	const pat::Electron *	wenu,
		int	i
	)		`[private]`

Definition at line 312 of file ZeePlots.cc.

References CutVars_, i, nBarrelVars_, ReturnCandVar(), and pat::Electron::superCluster().

Referenced by CheckCuts(), and CheckCutsNminusOne().

                                                     {
  double fabseta = fabs(ele->superCluster()->eta());
  if ( fabseta<1.479) {
    return fabs(ReturnCandVar(ele, i)) < CutVars_[i];
  }
  return fabs(ReturnCandVar(ele, i)) < CutVars_[i+nBarrelVars_];
}

bool ZeePlots::CheckCut2	(	const pat::Electron *	wenu,
		int	i
	)		`[private]`

Definition at line 376 of file ZeePlots.cc.

References CutVars2_, i, nBarrelVars_, ReturnCandVar(), and pat::Electron::superCluster().

Referenced by CheckCuts2(), and CheckCuts2NminusOne().

                                                      {
  double fabseta = fabs(ele->superCluster()->eta());
  if ( fabseta<1.479) {
    return fabs(ReturnCandVar(ele, i)) < CutVars2_[i];
  }
  return fabs(ReturnCandVar(ele, i)) < CutVars2_[i+nBarrelVars_];
}

bool ZeePlots::CheckCut2Inv	(	const pat::Electron *	wenu,
		int	i
	)		`[private]`

Definition at line 384 of file ZeePlots.cc.

References CutVars2_, i, InvVars_, nBarrelVars_, ReturnCandVar(), and pat::Electron::superCluster().

Referenced by CheckCuts2Inverse().

                                                         {
  double fabseta = fabs(ele->superCluster()->eta());
  if ( fabseta<1.479) {
    if (InvVars_[i]) return fabs(ReturnCandVar(ele, i))>CutVars2_[i];
    return fabs(ReturnCandVar(ele, i)) < CutVars2_[i];
  }
  if (InvVars_[i+nBarrelVars_]) {
    if (InvVars_[i])
      return fabs(ReturnCandVar(ele, i))>CutVars2_[i+nBarrelVars_];
  }
  return fabs(ReturnCandVar(ele, i)) < CutVars2_[i+nBarrelVars_];
}

bool ZeePlots::CheckCutInv	(	const pat::Electron *	wenu,
		int	i
	)		`[private]`

Definition at line 320 of file ZeePlots.cc.

References CutVars_, i, InvVars_, nBarrelVars_, ReturnCandVar(), and pat::Electron::superCluster().

Referenced by CheckCutsInverse().

                                                        {
  double fabseta = fabs(ele->superCluster()->eta());
  if ( fabseta<1.479) {
    if (InvVars_[i]) return fabs(ReturnCandVar(ele, i))>CutVars_[i];
    return fabs(ReturnCandVar(ele, i)) < CutVars_[i];
  }
  if (InvVars_[i+nBarrelVars_]) {
    if (InvVars_[i])
      return fabs(ReturnCandVar(ele, i))>CutVars_[i+nBarrelVars_];
  }
  return fabs(ReturnCandVar(ele, i)) < CutVars_[i+nBarrelVars_];
}

bool ZeePlots::CheckCuts ( const pat::Electron * ele ) [private]

Definition at line 285 of file ZeePlots.cc.

References CheckCut(), i, and nBarrelVars_.

Referenced by analyze().

{
  for (int i=0; i<nBarrelVars_; ++i) {
    if (not CheckCut(ele, i)) return false;
  }
  return true;
}

bool ZeePlots::CheckCuts2 ( const pat::Electron * ele ) [private]

Definition at line 349 of file ZeePlots.cc.

References CheckCut2(), i, and nBarrelVars_.

Referenced by analyze().

{
  for (int i=0; i<nBarrelVars_; ++i) {
    if (not CheckCut2(ele, i)) return false;
  }
  return true;
}

bool ZeePlots::CheckCuts2Inverse ( const pat::Electron * ele ) [private]

Definition at line 358 of file ZeePlots.cc.

References CheckCut2Inv(), i, and nBarrelVars_.

{
  for (int i=0; i<nBarrelVars_; ++i){
    if ( CheckCut2Inv(ele, i) == false) return false;
  }
  return true;

}

bool ZeePlots::CheckCuts2NminusOne	(	const pat::Electron *	ele,
		int	jj
	)		`[private]`

Definition at line 367 of file ZeePlots.cc.

References CheckCut2(), i, and nBarrelVars_.

{
  for (int i=0; i<nBarrelVars_; ++i){
    if (i==jj) continue;
    if ( CheckCut2(ele, i) == false) return false;
  }
  return true;
}

bool ZeePlots::CheckCutsInverse ( const pat::Electron * ele ) [private]

Definition at line 294 of file ZeePlots.cc.

References CheckCutInv(), i, and nBarrelVars_.

{
  for (int i=0; i<nBarrelVars_; ++i){
    if ( CheckCutInv(ele, i) == false) return false;
  }
  return true;

}

bool ZeePlots::CheckCutsNminusOne	(	const pat::Electron *	ele,
		int	jj
	)		`[private]`

Definition at line 303 of file ZeePlots.cc.

References CheckCut(), i, and nBarrelVars_.

{
  for (int i=0; i<nBarrelVars_; ++i){
    if (i==jj) continue;
    if ( CheckCut(ele, i) == false) return false;
  }
  return true;
}

void ZeePlots::endJob ( void ) [private, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 507 of file ZeePlots.cc.

References h_e_ETA, h_e_PHI, h_e_PT, h_EB_deta, h_EB_dphi, h_EB_ecaliso, h_EB_hcaliso, h_EB_HoE, h_EB_sIetaIeta, h_EB_trkiso, h_EE_deta, h_EE_dphi, h_EE_ecaliso, h_EE_hcaliso, h_EE_HoE, h_EE_sIetaIeta, h_EE_trkiso, h_mee, h_mee_EBEB, h_mee_EBEE, h_mee_EEEE, h_Zcand_PT, h_Zcand_Y, runonSM::newfile, and outputFile_.

                 {
  TFile * newfile = new TFile(TString(outputFile_),"RECREATE");
  //
  h_mee->Write();
  h_mee_EBEB->Write();
  h_mee_EBEE->Write();
  h_mee_EEEE->Write();
  h_Zcand_PT->Write();
  h_Zcand_Y->Write();

  h_e_PT->Write();
  h_e_ETA->Write();
  h_e_PHI->Write();


  h_EB_trkiso->Write();
  h_EB_ecaliso->Write();
  h_EB_hcaliso->Write();
  h_EB_sIetaIeta->Write();
  h_EB_dphi->Write();
  h_EB_deta->Write();
  h_EB_HoE->Write();

  h_EE_trkiso->Write();
  h_EE_ecaliso->Write();
  h_EE_hcaliso->Write();
  h_EE_sIetaIeta->Write();
  h_EE_dphi->Write();
  h_EE_deta->Write();
  h_EE_HoE->Write();

  //
  newfile->Close();

}

void ZeePlots::fillEleClassesPlots ( calib::CalibElectron * myEle )

Definition at line 389 of file ZeePlots.cc.

References reco::GsfElectron::classification(), gather_cfg::cout, reco::LeafCandidate::eta(), calib::CalibElectron::getRecoElectron(), h1_occupancyVsEtaCrack_, h1_occupancyVsEtaGold_, h1_occupancyVsEtaShower_, and h1_occupancyVsEtaSilver_.

                                                             {

  int myClass = myEle->getRecoElectron()->classification();

  float myEta = myEle->getRecoElectron()->eta();
  
  if(myClass==0 || myClass==100)
    h1_occupancyVsEtaGold_->Fill(myEta);
  
  std::cout<< "[ZeePlots::fillEleClassesPlots]Done gold"<< std::endl;
  
  if(myClass==40 || myClass==140)
    h1_occupancyVsEtaCrack_->Fill(myEta);
  
  std::cout<< "[ZeePlots::fillEleClassesPlots]Done crack"<< std::endl;
  
  if( (myClass>=30 && myClass<=34) || (myClass>=130 && myClass<=134) )
    h1_occupancyVsEtaShower_->Fill(myEta);
  
  std::cout<< "[ZeePlots::fillEleClassesPlots]Done shower"<< std::endl;
  
  if( myClass==10 || myClass==20 || myClass==110 || myClass ==120)
    h1_occupancyVsEtaSilver_->Fill(myEta);
  
  std::cout<< "[ZeePlots::fillEleClassesPlots]Done"<< std::endl;
  
}

void ZeePlots::fillEleInfo ( const reco::GsfElectronCollection * electronCollection )

Definition at line 261 of file ZeePlots.cc.

References file_, h1_nEleReco_, h1_recoEleEnergy_, h1_recoEleEta_, h1_recoElePhi_, and h1_recoElePt_.

Referenced by ZeeCalibration::duringLoop().

                                                                              {

  file_->cd();

  h1_nEleReco_->Fill(electronCollection->size());
  
  for(reco::GsfElectronCollection::const_iterator eleIt = electronCollection->begin();   eleIt != electronCollection->end(); eleIt++)
    {
      
  file_->cd();

      h1_recoEleEnergy_->Fill( eleIt->superCluster()->energy() );
      h1_recoElePt_->Fill( eleIt->pt() );
      h1_recoEleEta_->Fill( eleIt->eta() );
      h1_recoElePhi_->Fill( eleIt->phi() );
      
    }//end loop on electrons   

}

void ZeePlots::fillEleMCInfo ( const HepMC::GenEvent * myGenEvent )

Definition at line 210 of file ZeePlots.cc.

References abs, file_, h1_mcEle_Energy_, h1_mcEleEta_, h1_mcElePhi_, h1_mcElePt_, and L1TEmulatorMonitor_cff::p.

Referenced by ZeeCalibration::duringLoop().

                                                              {

  file_->cd();
  
  for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
        p != myGenEvent->particles_end(); ++p ) {
    
    if (  abs( (*p)->pdg_id() ) == 11 )
      {
        h1_mcEle_Energy_->Fill( (*p)->momentum().e() );
        h1_mcElePt_->Fill( (*p)->momentum().perp() );
        h1_mcEleEta_->Fill( (*p)->momentum().eta() );
        h1_mcElePhi_->Fill( (*p)->momentum().phi() );
        
      }//matches if (  abs( (*p)->pdg_id() ) == 11 )

  }//end loop over MC particles
  
}

void ZeePlots::fillHLTInfo ( edm::Handle< edm::TriggerResults > hltTriggerResultHandle )

Definition at line 341 of file ZeePlots.cc.

References file_, h1_FiredTriggers_, h1_HLT1Electron_FiredEvents_, h1_HLT1Electron_HLT2Electron_FiredEvents_, h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_, h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents_, h1_HLT2Electron_FiredEvents_, h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_, h1_HLT2ElectronRelaxed_FiredEvents_, h1_HLTVisitedEvents_, and i.

Referenced by ZeeCalibration::duringLoop().

                                                                               {

  file_->cd();
  
  int hltCount = hltTriggerResultHandle->size();

  bool aHLTResults[200] = { false };
    
  for(int i = 0 ; i < hltCount ; i++) {

    aHLTResults[i] = hltTriggerResultHandle->accept(i);
    if(aHLTResults[i])
      h1_FiredTriggers_->Fill(i);

    //HLT bit 32 = HLT1Electron
    //HLT bit 34 = HLT2Electron
    //HLT bit 35 = HLT2ElectronRelaxed

  }

  h1_HLTVisitedEvents_->Fill(1);
    
  if(aHLTResults[32] && !aHLTResults[34] && !aHLTResults[35])
    h1_HLT1Electron_FiredEvents_->Fill(1);

  if(aHLTResults[34] && !aHLTResults[32] && !aHLTResults[35])
    h1_HLT2Electron_FiredEvents_->Fill(1);

  if(aHLTResults[35] && !aHLTResults[32] && !aHLTResults[34])
    h1_HLT2ElectronRelaxed_FiredEvents_->Fill(1);

  if(aHLTResults[32] && aHLTResults[34] && !aHLTResults[35])
    h1_HLT1Electron_HLT2Electron_FiredEvents_->Fill(1);

  if(aHLTResults[32] && aHLTResults[35] && !aHLTResults[34])
    h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents_->Fill(1);

  if(aHLTResults[34] && aHLTResults[35] && !aHLTResults[32])
    h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_->Fill(1);

  if(aHLTResults[32] && aHLTResults[34] && aHLTResults[35])
    h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents_->Fill(1);



}

void ZeePlots::fillZInfo ( std::pair< calib::CalibElectron *, calib::CalibElectron * > myZeeCandidate )

Definition at line 118 of file ZeePlots.cc.

References ZeeKinematicTools::calculateZEta(), ZeeKinematicTools::calculateZPhi(), ZeeKinematicTools::calculateZPt(), ZeeKinematicTools::calculateZRapidity(), ZeeKinematicTools::calculateZTheta(), h1_reco_ZEta_, h1_reco_ZPhi_, h1_reco_ZPt_, h1_reco_ZRapidity_, and h1_reco_ZTheta_.

Referenced by ZeeCalibration::duringLoop().

                                                                                          {

  h1_reco_ZEta_->Fill( ZeeKinematicTools::calculateZEta(myZeeCandidate) );
  h1_reco_ZTheta_->Fill( ZeeKinematicTools::calculateZTheta(myZeeCandidate) );
  h1_reco_ZRapidity_->Fill( ZeeKinematicTools::calculateZRapidity(myZeeCandidate) );
  h1_reco_ZPhi_->Fill( ZeeKinematicTools::calculateZPhi(myZeeCandidate) );
  h1_reco_ZPt_->Fill( ZeeKinematicTools::calculateZPt(myZeeCandidate) );

}

void ZeePlots::fillZMCInfo ( const HepMC::GenEvent * myGenEvent )

Definition at line 158 of file ZeePlots.cc.

References file_, h1_gen_ZEta_, h1_gen_ZMass_, h1_gen_ZPhi_, h1_gen_ZPt_, h1_gen_ZRapidity_, funct::log(), and L1TEmulatorMonitor_cff::p.

Referenced by ZeeCalibration::duringLoop().

                                                            {

  file_->cd();

  for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
        p != myGenEvent->particles_end(); ++p ) {//loop over MC particles
 
    if ( (*p)->pdg_id() == 23 && (*p)->status() == 2 ){
      
      h1_gen_ZMass_->Fill( (*p)->momentum().m() );
      h1_gen_ZEta_->Fill( (*p)->momentum().eta() );
      
      float genZ_Y = 0.5 * log ( ( (*p)->momentum().e() + (*p)->momentum().pz() ) /  ( (*p)->momentum().e() - (*p)->momentum().pz() ) )   ;
      
      h1_gen_ZRapidity_->Fill( genZ_Y );
      h1_gen_ZPt_->Fill((*p)->momentum().perp());
      h1_gen_ZPhi_->Fill((*p)->momentum().phi());

      

    }
  }//end loop over MC particles

  return;  
}

void ZeePlots::openFile ( )

Definition at line 54 of file ZeePlots.cc.

References file_.

                       {


  file_ -> cd();

}

double ZeePlots::ReturnCandVar	(	const pat::Electron *	ele,
		int	i
	)		`[private]`

Definition at line 333 of file ZeePlots.cc.

References gather_cfg::cout, reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), pat::EcalIso, reco::GsfElectron::hadronicOverEm(), pat::HcalIso, reco::GsfElectron::scSigmaIEtaIEta(), pat::TrackIso, pat::User1Iso, and pat::Lepton< LeptonType >::userIsolation().

Referenced by CheckCut(), CheckCut2(), CheckCut2Inv(), and CheckCutInv().

                                                            {
  if (i==0) return ele->userIsolation(pat::TrackIso);
  else if (i==1) return ele->userIsolation(pat::EcalIso);
  else if (i==2) return ele->userIsolation(pat::HcalIso);
  else if (i==3) return ele->scSigmaIEtaIEta();
  else if (i==4) return ele->deltaPhiSuperClusterTrackAtVtx();
  else if (i==5) return ele->deltaEtaSuperClusterTrackAtVtx();
  else if (i==6) return ele->hadronicOverEm();
  else if (i==7) return ele->userIsolation(pat::User1Iso);
  std::cout << "Error in ZeePlots::ReturnCandVar" << std::endl;
  return -1.;

}

void ZeePlots::writeEleClassesPlots ( )

Definition at line 440 of file ZeePlots.cc.

References file_, h1_occupancyVsEtaCrack_, h1_occupancyVsEtaGold_, h1_occupancyVsEtaShower_, and h1_occupancyVsEtaSilver_.

                                   {

  file_->cd();

  h1_occupancyVsEtaGold_->Write();
  h1_occupancyVsEtaSilver_->Write();
  h1_occupancyVsEtaShower_->Write();
  h1_occupancyVsEtaCrack_->Write();

}

void ZeePlots::writeEleHistograms ( )

Definition at line 283 of file ZeePlots.cc.

References gather_cfg::cout, file_, h1_recoEleEnergy_, h1_recoEleEta_, h1_recoElePhi_, and h1_recoElePt_.

Referenced by ZeeCalibration::endOfJob().

                                 {

  file_->cd();

  std::cout << "Start with ZeePlots::writeEleHistograms(), done file_->cd(); " << std::endl;
 
  h1_recoEleEnergy_->Write();
  h1_recoElePt_->Write();
  h1_recoEleEta_->Write();
  h1_recoElePhi_->Write();

  std::cout << "Done with ZeePlots::writeEleHistograms() " << std::endl;

}

void ZeePlots::writeHLTHistograms ( )

void ZeePlots::writeMCEleHistograms ( )

Definition at line 300 of file ZeePlots.cc.

References gather_cfg::cout, file_, h1_mcEle_Energy_, h1_mcEleEta_, h1_mcElePhi_, and h1_mcElePt_.

Referenced by ZeeCalibration::endOfJob().

                                   {

  file_->cd();

  std::cout << "Start with ZeePlots::writeMCEleHistograms(), done file_->cd(); " << std::endl;

  h1_mcEle_Energy_->Write();
  h1_mcElePt_->Write();
  h1_mcEleEta_->Write();
  h1_mcElePhi_->Write();

  std::cout << "Done with ZeePlots::writeMCEleHistograms() " << std::endl;

}

void ZeePlots::writeMCZHistograms ( )

Definition at line 146 of file ZeePlots.cc.

References file_, h1_gen_ZPhi_, h1_gen_ZPt_, and h1_gen_ZRapidity_.

Referenced by ZeeCalibration::endOfJob().

                                  {

  file_->cd();
  
  h1_gen_ZRapidity_->Write();
  h1_gen_ZPt_->Write();
  h1_gen_ZPhi_->Write();

}

void ZeePlots::writeZHistograms ( )

Definition at line 132 of file ZeePlots.cc.

References file_, h1_reco_ZEta_, h1_reco_ZPhi_, h1_reco_ZPt_, h1_reco_ZRapidity_, and h1_reco_ZTheta_.

Referenced by ZeeCalibration::endOfJob().

                                {

  file_->cd();

  h1_reco_ZEta_->Write();
  h1_reco_ZTheta_->Write();
  h1_reco_ZRapidity_->Write();
  h1_reco_ZPhi_->Write();
  h1_reco_ZPt_->Write();
  
}