---

SUSYControlHighPtPhotonSkim.cc

Go to the documentation of this file.

/* \class SUSYControlHighPtPhotonSkim
 *
 * High Energy Photon SUSY Skim (control sample)
 * one photon and one electron > xx GeV in barrel + isolation 
 *
 * $Date: 2007/12/09 10:54:57 $
 * $Revision: 1.6 $
 *
 * \author Daniele del Re - Univ. La Sapienza & INFN
 *
 */

#include <iostream>
#include <string>
#include <list>
#include <cmath>
#include <cstdio>
#include <vector>
#include <memory>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/Common/interface/Handle.h"    

#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"

#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Common/interface/AssociationVector.h"

#include "SUSYBSMAnalysis/CSA07Skims/interface/SUSYControlHighPtPhotonSkim.h"

using namespace edm;
using namespace std;
using namespace reco;

SUSYControlHighPtPhotonSkim::SUSYControlHighPtPhotonSkim( const edm::ParameterSet& iConfig ) :
  nEvents_(0), nAccepted_(0)
{
  Photonsrc_ = iConfig.getParameter<InputTag>( "Photonsrc" );
  Electronsrc_ = iConfig.getParameter<InputTag>( "Electronsrc" );
  PhotonPtmin_ = 
    iConfig.getParameter<double>( "PhotonPtmin");
  ElectronPtmin_ = 
    iConfig.getParameter<double>( "ElectronPtmin");
  IsIsolated_ = iConfig.getParameter<bool>( "IsIsolated");
  IsolationCut_ = iConfig.getParameter<double>( "IsolationCut");
}

/*------------------------------------------------------------------------*/

SUSYControlHighPtPhotonSkim::~SUSYControlHighPtPhotonSkim() 
{}

/*------------------------------------------------------------------------*/

bool SUSYControlHighPtPhotonSkim::filter( edm::Event& iEvent, 
                                       const edm::EventSetup& iSetup )
{
  nEvents_++;

  typedef AssociationVector<RefProd<CandidateCollection>, vector<double> > PhotonMapCollection;
  Handle<PhotonMapCollection>  PhotonHandle;

  iEvent.getByLabel( Photonsrc_, PhotonHandle );

  if ( PhotonHandle->empty() ) return false;

  Handle<GsfElectronCollection>  ElectronHandle;

  iEvent.getByLabel( Electronsrc_, ElectronHandle );

  if ( ElectronHandle->empty() ) return false;

  int nPhoton = 0;
  int nElectron = 0;

  for ( PhotonMapCollection::const_iterator it = PhotonHandle->begin(); 
        it != PhotonHandle->end(); it++ ) {

    bool iso = it->second < IsolationCut_;
    if(!IsIsolated_) iso = 1; 

    if (iso && fabs(it->first->eta()) < 1.479 && it->first->pt() > PhotonPtmin_) { 

      int overlap(0);
      for ( GsfElectronCollection::const_iterator itel = ElectronHandle->begin(); 
            itel != ElectronHandle->end(); itel++ ) {
        
        if( itel->pt() > ElectronPtmin_ ){
          double dr = sqrt( pow(it->first->eta() - itel->eta(),2) +
                            pow(it->first->phi() - itel->phi(),2) );
          if ( dr<0.1 ) overlap = 1;
        }
      }
      if (!overlap) nPhoton++;

    }

  }

  for ( GsfElectronCollection::const_iterator itel = ElectronHandle->begin(); 
        itel != ElectronHandle->end(); itel++ ) {

    if (fabs(itel->eta()) < 1.479 && itel->pt() > ElectronPtmin_) 
      nElectron++;
 }

  if (!nPhoton) return false;
  if (!nElectron) return false;

  nAccepted_++;

  return true;
}

/*------------------------------------------------------------------------*/

void SUSYControlHighPtPhotonSkim::endJob()
{
  edm::LogVerbatim( "SUSYControlHighPtPhotonSkim" ) 
    << "Events read " << nEvents_
    << " Events accepted " << nAccepted_
    << "\nEfficiency " << (double)(nAccepted_)/(double)(nEvents_) 
    << endl;
}

---