ZMuMuAnalyzer Class Reference

Inheritance diagram for ZMuMuAnalyzer:

Public Member Functions
	ZMuMuAnalyzer (const edm::ParameterSet &pset)
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup)
virtual void	endJob ()

Private Attributes
double	etacut_
TH1D *	h_zMuMu_mass_
TH1D *	h_zMuMuMatched_mass_
TH1D *	h_zMuSingleStandAlone_mass_
TH1D *	h_zMuSingleStandAloneMatched_mass_
TH1D *	h_zMuSingleStandAloneOverlap_mass_
TH1D *	h_zMuSingleStandAloneOverlapMatched_mass_
TH1D *	h_zMuSingleTrack_mass_
TH1D *	h_zMuSingleTrackMatched_mass_
double	isocut_
double	maxZmass_
double	minZmass_
InputTag	muIso_
OverlapChecker	overlap_
double	ptcut_
double	ptSTAcut_
InputTag	standAloneIso_
InputTag	trackIso_
InputTag	zMuMu_
InputTag	zMuMuMap_
InputTag	zMuStandAlone_
InputTag	zMuStandAloneMap_
InputTag	zMuTrack_
InputTag	zMuTrackMap_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 41 of file ZMuMuAnalyzer.cc.

Constructor & Destructor Documentation

ZMuMuAnalyzer::ZMuMuAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 62 of file ZMuMuAnalyzer.cc.

References h_zMuMu_mass_, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleStandAloneOverlap_mass_, h_zMuSingleStandAloneOverlapMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, and TFileDirectory::make().

                                                        : 
  zMuMu_( pset.getParameter<InputTag>( "zMuMu" ) ),
  zMuTrack_( pset.getParameter<InputTag>( "zMuTrack" ) ),
  zMuStandAlone_( pset.getParameter<InputTag>( "zMuStandAlone" ) ),
  muIso_( pset.getParameter<InputTag>( "muIso" ) ),
  trackIso_( pset.getParameter<InputTag>( "trackIso" ) ),
  standAloneIso_( pset.getParameter<InputTag>( "standAloneIso" ) ),
  zMuMuMap_( pset.getParameter<InputTag>( "zMuMuMap" ) ),
  zMuTrackMap_( pset.getParameter<InputTag>( "zMuTrackMap" ) ),
  zMuStandAloneMap_( pset.getParameter<InputTag>( "zMuStandAloneMap" ) ),
  isocut_( pset.getParameter<double>( "isocut" ) ),
  etacut_( pset.getParameter<double>( "etacut" ) ),
  ptcut_( pset.getParameter<double>( "ptcut" ) ),
  ptSTAcut_( pset.getParameter<double>( "ptSTAcut" ) ),
  
  minZmass_( pset.getParameter<double>( "minZmass" )),
  maxZmass_( pset.getParameter<double>( "maxZmass" )) {
  
  Service<TFileService> fs;
  h_zMuMu_mass_ = fs->make<TH1D>( "ZMuMumass", "ZMuMu mass(GeV)", 200,  0., 200. );
  h_zMuSingleTrack_mass_ = fs->make<TH1D>( "ZMuSingleTrackmass", "ZMuSingleTrack mass(GeV)", 100,  0., 200. );
  h_zMuSingleStandAlone_mass_ = fs->make<TH1D>( "ZMuSingleStandAlonemass", "ZMuSingleStandAlone mass(GeV)", 50,  0., 200. );
  h_zMuSingleStandAloneOverlap_mass_ = fs->make<TH1D>( "ZMuSingleStandAloneOverlapmass", "ZMuSingleStandAloneOverlap  mass(GeV)", 50,  0., 200. );
  
  
  h_zMuMuMatched_mass_ = fs->make<TH1D>( "ZMuMuMatchedmass", "ZMuMu Matched  mass(GeV)", 200,  0., 200. );
  h_zMuSingleTrackMatched_mass_ = fs->make<TH1D>( "ZMuSingleTrackmassMatched", "ZMuSingleTrackMatched mass(GeV)", 100,  0., 200. );
  h_zMuSingleStandAloneMatched_mass_ = fs->make<TH1D>( "ZMuSingleStandAlonemassMatched", "ZMuSingleStandAloneMatched mass(GeV)", 50,  0., 200. );
  h_zMuSingleStandAloneOverlapMatched_mass_ = fs->make<TH1D>( "ZMuSingleStandAloneOverlapmassMatched", "ZMuSingleStandAloneMatched Overlap  mass(GeV)", 50,  0., 200. );
}

Member Function Documentation

void ZMuMuAnalyzer::analyze ( const edm::Event &  event, const edm::EventSetup &  setup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 93 of file ZMuMuAnalyzer.cc.

References edm::RefToBase< T >::castTo(), reco::Candidate::daughter(), etacut_, h_zMuMu_mass_, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleStandAloneOverlap_mass_, h_zMuSingleStandAloneOverlapMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, i, edm::ProductID::id(), edm::HandleBase::id(), edm::Ref< C, T, F >::id(), isocut_, j, edm::Ref< C, T, F >::key(), m, reco::Candidate::mass(), reco::Candidate::masterClone(), maxZmass_, minZmass_, muIso_, overlap_, ptcut_, ptSTAcut_, standAloneIso_, trackIso_, ZMuMuAnalysisNtupler_cff::zMuMu, zMuMu_, zMuMuMap_, zMuStandAlone_, zMuStandAloneMap_, zMuTrack_, and zMuTrackMap_.

                                                                             {
  Handle<CandidateCollection> zMuMu;
  event.getByLabel(zMuMu_, zMuMu);
  Handle<CandidateCollection> zMuTrack;
  event.getByLabel( zMuTrack_, zMuTrack );
  Handle<CandidateCollection> zMuStandAlone;
  event.getByLabel( zMuStandAlone_, zMuStandAlone );  

  unsigned int nZMuMu = zMuMu->size();
  unsigned int nZTrackMu = zMuTrack->size();
  unsigned int nZStandAloneMu = zMuStandAlone->size();
  static const double zMass = 91.1876; // PDG Z mass

  //  cout << "nZMuMu = " << nZMuMu << endl;
  //  cout << "nZTrackMu = " << nZTrackMu << endl;
  //  cout << "nZStandAloneMu = " << nZStandAloneMu << endl;
 
  Handle<CandMatchMap> zMuMuMap;
  if( nZMuMu > 0 ) {
    event.getByLabel(zMuMuMap_, zMuMuMap);
  }

  Handle<CandMatchMap> zMuTrackMap;
  if( nZTrackMu > 0 ) {
    event.getByLabel( zMuTrackMap_, zMuTrackMap );
  }

  Handle<CandMatchMap> zMuStandAloneMap;
  if( nZStandAloneMu > 0 ) {
    event.getByLabel( zMuStandAloneMap_, zMuStandAloneMap );  
  }    

  Handle<IsolationCollection> muIso;
  event.getByLabel(muIso_, muIso);
  ProductID muIsoId = muIso->keyProduct().id();
  Handle<IsolationCollection> trackIso;
  event.getByLabel(trackIso_, trackIso);
  ProductID trackIsoId = trackIso->keyProduct().id();
  
  Handle<IsolationCollection> standAloneIso;
  event.getByLabel(standAloneIso_, standAloneIso);
  ProductID standAloneIsoId = standAloneIso->keyProduct().id();
  
  if (nZMuMu > 0) {
    double mass = 1000000.;
    for( unsigned int i = 0; i < nZMuMu; i++ ) {
      const Candidate & zmmCand = (*zMuMu)[ i ];
      CandidateRef CandRef(zMuMu,i);
      CandidateRef lep1 = zmmCand.daughter( 0 )->masterClone().castTo<CandidateRef>();
      CandidateRef lep2 = zmmCand.daughter( 1 )->masterClone().castTo<CandidateRef>();
      
      const  double iso1 = muIso->value( lep1.key() );  
      const  double iso2 = muIso->value( lep2.key() );  
      
      double m = zmmCand.mass();
      if (lep1->pt()>ptcut_ && lep2->pt()>ptcut_ &&    
      fabs(lep1->eta())<etacut_ && fabs(lep2->eta())<etacut_ &&
      m>minZmass_ && m<maxZmass_ && iso1 < isocut_ && iso2 <isocut_) {
    if ( fabs( mass - zMass ) > fabs( m - zMass ) ) {
      mass = m;
    }
    
    h_zMuMu_mass_->Fill( mass );      
    CandMatchMap::const_iterator m0 = zMuMuMap->find(CandRef);
    if( m0 != zMuMuMap->end()) {
        h_zMuMuMatched_mass_->Fill( mass );
    }
      }
    }
  }
  
  //ZmuSingleTRack
  if (nZMuMu ==0 && nZTrackMu>0) {
    for( unsigned int j = 0; j < nZTrackMu; j++ ) {
      const Candidate & ztmCand = (*zMuTrack)[ j ];
      CandidateRef CandRef(zMuTrack,j);
      CandidateRef lep1 = ztmCand.daughter( 0 )->masterClone().castTo<CandidateRef>();
      CandidateRef lep2 = ztmCand.daughter( 1 )->masterClone().castTo<CandidateRef>();
      
      ProductID id1 = lep1.id();
      ProductID id2 = lep2.id();
      double iso1 = -1;
      double iso2 = -1;
      
      if( id1 == muIsoId ) 
    iso1 = muIso->value( lep1.key() );  
      else if ( id1 == trackIsoId )
    iso1 = trackIso->value( lep1.key() );   
      
      if( id2 == muIsoId ) 
          iso2 = muIso->value( lep2.key() );    
      else if ( id2 == trackIsoId )
    iso2 = trackIso->value( lep2.key() );   
 
      double mt = ztmCand.mass();
      if (lep1->pt()>ptcut_ && lep2->pt()>ptcut_ &&    
      fabs(lep1->eta())<etacut_ && fabs(lep2->eta())<etacut_ &&
      mt>minZmass_ && mt<maxZmass_ && iso1<isocut_ && iso2 <isocut_) {
    h_zMuSingleTrack_mass_->Fill( mt );
    CandMatchMap::const_iterator m0 = zMuTrackMap->find(CandRef);
    if( m0 != zMuTrackMap->end()) {
      h_zMuSingleTrackMatched_mass_->Fill( mt );
    }
      }
    }
  }   
  
  //ZmuSingleStandAlone
  if (nZMuMu ==0 && nZStandAloneMu>0) {
    //      unsigned int index = 1000;
    for( unsigned int j = 0; j < nZStandAloneMu; j++ ) {
      const Candidate & zsmCand = (*zMuStandAlone)[ j ];
      CandidateRef CandRef(zMuStandAlone,j);
      CandidateRef lep1 = zsmCand.daughter( 0 )->masterClone().castTo<CandidateRef>();
      CandidateRef lep2 = zsmCand.daughter( 1 )->masterClone().castTo<CandidateRef>();
      
      ProductID id1 = lep1.id();
      ProductID id2 = lep2.id();
      double iso1 = -1;
      double iso2 = -1;
      
      if( id1 == muIsoId ) 
    iso1 = muIso->value( lep1.key() );  
      else if ( id1 == standAloneIsoId )
    iso1 = standAloneIso->value( lep1.key() );  
      
      if( id2 == muIsoId ) 
    iso2 = muIso->value( lep2.key() );  
      else if ( id2 == standAloneIsoId )
    iso2 = standAloneIso->value( lep2.key() );  
      
      double ms = zsmCand.mass();
      if (lep1->pt()>ptSTAcut_ && lep2->pt()>ptSTAcut_ &&    
      fabs(lep1->eta())<etacut_ && fabs(lep2->eta())<etacut_ &&
      ms>minZmass_ && ms<maxZmass_ && iso1<isocut_ && iso2 <isocut_) {
    h_zMuSingleStandAlone_mass_->Fill( ms );
    CandMatchMap::const_iterator m0 = zMuStandAloneMap->find(CandRef);
    if( m0 != zMuStandAloneMap->end()) {
      h_zMuSingleStandAloneMatched_mass_->Fill( ms );
    }
    
    bool noOverlap = true;
    for( unsigned int j = 0; j < zMuTrack->size(); j++ ) {
      const Candidate & ztmCand = (*zMuTrack)[ j ];
      CandidateRef CandReft(zMuTrack,j);
      
      CandidateRef lep1 = ztmCand.daughter( 0 )->masterClone().castTo<CandidateRef>();
      CandidateRef lep2 = ztmCand.daughter( 1 )->masterClone().castTo<CandidateRef>();
      
      ProductID id1 = lep1.id();
      ProductID id2 = lep2.id();
      double iso1 = -1;
      double iso2 = -1;
      
      if( id1 == muIsoId ) 
        iso1 = muIso->value( lep1.key() );  
      else if ( id1 == trackIsoId )
        iso1 = trackIso->value( lep1.key() );   
      
      if( id2 == muIsoId ) 
        iso2 = muIso->value( lep2.key() );  
      else if ( id2 == trackIsoId )
        iso2 = trackIso->value( lep2.key() );   
      
      double mt = ztmCand.mass();
      if (lep1->pt()>ptcut_ && lep2->pt()>ptcut_ &&    
          fabs(lep1->eta())<etacut_ && fabs(lep2->eta())<etacut_ &&
          mt>minZmass_ && mt<maxZmass_ && iso1<isocut_ && iso2 <isocut_) {
        
        if ( overlap_( ztmCand, zsmCand ) ) {
          noOverlap = false; 
          break;
        }
        if (!noOverlap ) {
          h_zMuSingleStandAloneOverlap_mass_->Fill( ms ); 
          CandMatchMap::const_iterator m1 = zMuTrackMap->find(CandReft);
          CandMatchMap::const_iterator m2 = zMuStandAloneMap->find(CandRef);
          
          if( m1 != zMuTrackMap->end()  && m2 != zMuStandAloneMap->end()  ) {
        h_zMuSingleStandAloneOverlapMatched_mass_->Fill( ms );
          }
        }
      }         
    }
      }
    }
  }
}  

void ZMuMuAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 282 of file ZMuMuAnalyzer.cc.

References gather_cfg::cout, h_zMuMu_mass_, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleStandAloneOverlap_mass_, h_zMuSingleStandAloneOverlapMatched_mass_, h_zMuSingleTrack_mass_, and h_zMuSingleTrackMatched_mass_.

                           {
  double Nzmm = h_zMuMu_mass_->GetEntries() ;
  double Nzsm = h_zMuSingleStandAlone_mass_->GetEntries()  ;
  double Nzsnom = h_zMuSingleStandAloneOverlap_mass_->GetEntries()  ;
  double Nztm = h_zMuSingleTrack_mass_->GetEntries();
  
  double NzmmMatch = h_zMuMuMatched_mass_->GetEntries() ;
  double NzsmMatch = h_zMuSingleStandAloneMatched_mass_->GetEntries()  ;
  double NzsnomMatch = h_zMuSingleStandAloneOverlapMatched_mass_->GetEntries()  ;
  double NztmMatch = h_zMuSingleTrackMatched_mass_->GetEntries();

  cout<<"-- N SingleTrackMu = "<<Nztm<<endl;
  cout<<"-----N SinglStandAloneMu = "<<Nzsm<<endl;
  cout<<"-----N SingleStandAloneOverlapMu = "<<Nzsnom<<endl;
  cout<<"------- N MuMu = "<<Nzmm<<endl;
  
  cout<<"-- N SingleTrackMuMatched = "<<NztmMatch<<endl;
  cout<<"-----N SinglStandAloneMuMatched = "<<NzsmMatch<<endl;
  cout<<"-----N SingleStandAloneOverlapMuMatched  = "<<NzsnomMatch<<endl;
  cout<<"------- N MuMu Matched  = "<<NzmmMatch<<endl;
}

Member Data Documentation

double ZMuMuAnalyzer::etacut_

private

Definition at line 54 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

TH1D* ZMuMuAnalyzer::h_zMuMu_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

TH1D * ZMuMuAnalyzer::h_zMuMuMatched_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

TH1D * ZMuMuAnalyzer::h_zMuSingleStandAlone_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

TH1D * ZMuMuAnalyzer::h_zMuSingleStandAloneMatched_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

TH1D * ZMuMuAnalyzer::h_zMuSingleStandAloneOverlap_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

TH1D * ZMuMuAnalyzer::h_zMuSingleStandAloneOverlapMatched_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

TH1D * ZMuMuAnalyzer::h_zMuSingleTrack_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

TH1D * ZMuMuAnalyzer::h_zMuSingleTrackMatched_mass_

private

Definition at line 55 of file ZMuMuAnalyzer.cc.

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

double ZMuMuAnalyzer::isocut_

private

Definition at line 54 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

double ZMuMuAnalyzer::maxZmass_

private

Definition at line 54 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

double ZMuMuAnalyzer::minZmass_

private

Definition at line 54 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::muIso_

private

Definition at line 52 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

OverlapChecker ZMuMuAnalyzer::overlap_

private

Definition at line 48 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

double ZMuMuAnalyzer::ptcut_

private

Definition at line 54 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

double ZMuMuAnalyzer::ptSTAcut_

private

Definition at line 54 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::standAloneIso_

private

Definition at line 52 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::trackIso_

private

Definition at line 52 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::zMuMu_

private

Definition at line 49 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::zMuMuMap_

private

Definition at line 53 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::zMuStandAlone_

private

Definition at line 51 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::zMuStandAloneMap_

private

Definition at line 53 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::zMuTrack_

private

Definition at line 50 of file ZMuMuAnalyzer.cc.

Referenced by analyze().

InputTag ZMuMuAnalyzer::zMuTrackMap_

private

Definition at line 53 of file ZMuMuAnalyzer.cc.

Referenced by analyze().