#include <L2TauAnalyzer.h>

Inheritance diagram for L2TauAnalyzer:

Public Member Functions
	L2TauAnalyzer (const edm::ParameterSet &)

	~L2TauAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

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

virtual void	beginJob ()

virtual void	endJob ()

MatchElementL2	match (const reco::Jet &, const LVColl &)

MatchElementL2	match (const reco::Jet &, const l1extra::L1JetParticleCollection &)

Private Attributes
float	cl_drRMS

float	cl_etaRMS

int	cl_Nclusters

float	cl_phiRMS

float	ecalIsol_Et

bool	IsSignal_

float	jetEMF

float	JetEt

float	JetEta

float	L1et

float	L1eta

edm::InputTag	l1Jets_

edm::InputTag	l1Taus_

TFile *	l2file

edm::InputTag	l2TauInfoAssoc_

TTree *	l2tree

double	matchDR_

edm::InputTag	mcColl_

float	MCet

float	MCeta

std::string	rootFile_

float	seedTowerEt

float	towerIsol_Et

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)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 39 of file L2TauAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 7 of file L2TauAnalyzer.cc.

References cl_drRMS, cl_etaRMS, cl_Nclusters, cl_phiRMS, ecalIsol_Et, jetEMF, JetEt, JetEta, L1et, L1eta, l2file, l2tree, MCet, MCeta, rootFile_, seedTowerEt, and towerIsol_Et.

                                                           :
  l2TauInfoAssoc_(iConfig.getParameter<edm::InputTag>("L2InfoAssociationInput")),
  l1Taus_(iConfig.getParameter<edm::InputTag>("L1TauCollection")),
  l1Jets_(iConfig.getParameter<edm::InputTag>("L1JetCollection")),
  rootFile_(iConfig.getParameter<std::string>("outputFileName")),
  IsSignal_(iConfig.getParameter<bool>("IsSignal")),
  mcColl_(iConfig.getParameter<edm::InputTag>("MatchedCollection"))
 {
   //File Setup
   l2file = new TFile(rootFile_.c_str(),"recreate");
   //Tree Setup
   l2tree = new TTree("l2tree","Level 2 Tau Tree");
 
 
   //Initialize the vars
   ecalIsol_Et=0.;
   towerIsol_Et=0.;
   cl_etaRMS=0.;
   cl_phiRMS=0.;
   cl_drRMS=0.;
   MCeta=0.;
   MCet=0.;
   cl_Nclusters=0;
   seedTowerEt = 0.;
   JetEt=0.;
   JetEta=0.;
   L1et=0.;
   L1eta=0.;
   jetEMF = 0.;
 
   //Setup Branches
   l2tree->Branch("ecalIsolEt",&ecalIsol_Et,"ecalIsolEt/F");
   l2tree->Branch("jetEMF",&jetEMF,"jetEMF/F");
   l2tree->Branch("towerIsolEt",&towerIsol_Et,"towerIsolEt/F");
   l2tree->Branch("clEtaRMS",&cl_etaRMS,"clEtaRMS/F");
   l2tree->Branch("clPhiRMS",&cl_phiRMS,"clPhiRMS/F");
   l2tree->Branch("clDrRMS",&cl_drRMS,"clDrRMS/F");
   l2tree->Branch("mcEta",&MCeta,"mcEta/F");
   l2tree->Branch("mcEt",&MCet,"mcEt/F");
   l2tree->Branch("clNclusters",&cl_Nclusters,"clNclusters/I");
   l2tree->Branch("seedTowerEt",&seedTowerEt,"seedTowerEt/F");
   l2tree->Branch("jetEt",&JetEt,"jetEt/F");
   l2tree->Branch("jetEta",&JetEta,"jetEta/F");
   l2tree->Branch("L1Et",&L1et,"L1Et/F");
   l2tree->Branch("L1Eta",&L1eta,"L1Eta/F");
  
 }

L2TauAnalyzer::~L2TauAnalyzer ( )

Definition at line 56 of file L2TauAnalyzer.cc.

57 {

58 }

Member Function Documentation

void L2TauAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 63 of file L2TauAnalyzer.cc.

References cl_drRMS, cl_etaRMS, cl_Nclusters, cl_phiRMS, MatchElementL2::deltar, reco::L2TauIsolationInfo::ecalClusterShape(), ecalIsol_Et, reco::L2TauIsolationInfo::ecalIsolEt(), reco::CaloJet::emEnergyFraction(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), edm::Event::getByLabel(), reco::L2TauIsolationInfo::hcalIsolEt(), IsSignal_, metsig::jet, jetEMF, JetEt, JetEta, L1et, L1eta, l1Jets_, l1Taus_, l2TauInfoAssoc_, l2tree, match(), MatchElementL2::matched, mcColl_, MatchElementL2::mcEt, MCet, MatchElementL2::mcEta, MCeta, reco::L2TauIsolationInfo::nEcalHits(), AlCaHLTBitMon_ParallelJobs::p, dt_dqm_sourceclient_common_cff::reco, reco::L2TauIsolationInfo::seedHcalHitEt(), seedTowerEt, and towerIsol_Et.

 {
    using namespace edm;
    using namespace reco;
 
    Handle<L2TauInfoAssociation> l2TauInfoAssoc; //Handle to the input (L2 Tau Info Association)
    Handle<LVColl> McInfo; //Handle To The Truth!!!!
    Handle<l1extra::L1JetParticleCollection> L1Taus; //Handle To The L1 Taus
    Handle<l1extra::L1JetParticleCollection> L1Jets; //Handle To The L1 jets
 
         if(iEvent.getByLabel(l2TauInfoAssoc_,l2TauInfoAssoc))//get the handle
       {
         if(l2TauInfoAssoc->size()>0)
           for(L2TauInfoAssociation::const_iterator p = l2TauInfoAssoc->begin();p!=l2TauInfoAssoc->end();++p)
         {
           const L2TauIsolationInfo l2info = p->val;
               const CaloJet& jet =*(p->key);
        
           MatchElementL2 mcMatch;
           mcMatch.matched=false;
           mcMatch.mcEt=0;
           mcMatch.mcEta=0;
           mcMatch.deltar=0;
 
           if(IsSignal_) //Get Collection and match it
             {
                   if(iEvent.getByLabel(mcColl_,McInfo))
               mcMatch=match(jet,*McInfo);
             }
 
           if((mcMatch.matched&&IsSignal_)||(!IsSignal_))
             {
               //Fill variables
               jetEMF = jet.emEnergyFraction();
               ecalIsol_Et=l2info.ecalIsolEt();
               towerIsol_Et=l2info.hcalIsolEt();
               cl_Nclusters=l2info.nEcalHits();
               cl_etaRMS=l2info.ecalClusterShape()[0];
               cl_phiRMS=l2info.ecalClusterShape()[1];
               cl_drRMS=l2info.ecalClusterShape()[2];
               seedTowerEt = l2info.seedHcalHitEt();
               MCeta =mcMatch.mcEta;
               MCet=mcMatch.mcEt;
               JetEt = jet.et();
               JetEta = jet.eta();
 
               //Match with L1 and fill
               L1et=0;
               L1eta=0;
               if(iEvent.getByLabel(l1Taus_,L1Taus))
             {
                 MatchElementL2 l1Match;
                 l1Match.matched=false;
                 l1Match.mcEt=0;
                 l1Match.mcEta=0;
                 l1Match.deltar=0;
                 l1Match=match(jet,*L1Taus);
                 if(l1Match.matched)
                   {
                 L1et=l1Match.mcEt;
                 L1eta=l1Match.mcEta;
                   }
                 //If not matched look at the jet collection
                 else
                   {
                 if(iEvent.getByLabel(l1Jets_,L1Jets))
                   {
                     l1Match=match(jet,*L1Taus);
                     if(l1Match.matched)
                       {
                     L1et=l1Match.mcEt;
                     L1eta=l1Match.mcEta;
                       }
 
                   }
                   }
 
             }
               //Fill Tree
               l2tree->Fill();
             }
        
         }
       } 
 }

void L2TauAnalyzer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 152 of file L2TauAnalyzer.cc.

 {
 
 }

void L2TauAnalyzer::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 159 of file L2TauAnalyzer.cc.

References l2file.

                       {
   l2file->Write();
 
 }

MatchElementL2 L2TauAnalyzer::match	(	const reco::Jet &	jet,
		const LVColl &	McInfo
	)

private

Definition at line 165 of file L2TauAnalyzer.cc.

References delta, MatchElementL2::deltar, MatchElementL2::matched, MatchElementL2::mcEt, MatchElementL2::mcEta, and reco::LeafCandidate::p4().

Referenced by analyze(), and match().

 {
 
   //Loop On the Collection and see if your tau jet is matched to one there
  //Also find the nearest Matched MC Particle to your Jet (to be complete)
  
  bool matched=false;
  double delta_min=100.;
  double mceta=0;
  double mcet=0;
  
  double matchingDR=0.3;
 
 
 
 
  if(McInfo.size()>0)
   for(std::vector<LV>::const_iterator it = McInfo.begin();it!=McInfo.end();++it)
    {
           double delta = ROOT::Math::VectorUtil::DeltaR(jet.p4().Vect(),*it);
       if(delta<matchingDR)
         {
           matched=true;
           if(delta<delta_min)
         {
           delta_min=delta;
           mceta=it->eta();
           mcet=it->Et();
         }
         }
    }
 
   //Create Struct and send it out!
   MatchElementL2 match;
   match.matched=matched;
   match.deltar=delta_min;
   match.mcEta = mceta;
   match.mcEt = mcet;
 
 
  return match;
 }

MatchElementL2 L2TauAnalyzer::match	(	const reco::Jet &	jet,
		const l1extra::L1JetParticleCollection &	McInfo
	)

private

Definition at line 209 of file L2TauAnalyzer.cc.

References delta, MatchElementL2::deltar, match(), MatchElementL2::matched, MatchElementL2::mcEt, MatchElementL2::mcEta, and reco::LeafCandidate::p4().

 {
 
   //Loop On the Collection and see if your tau jet is matched to one there
  //Also find the nearest Matched MC Particle to your Jet (to be complete)
  
  bool matched=false;
  double delta_min=100.;
  double mceta=0;
  double mcet=0;
  
  double matchingDR=0.5;
 
 
 
 
  if(McInfo.size()>0)
   for(l1extra::L1JetParticleCollection::const_iterator it = McInfo.begin();it!=McInfo.end();++it)
    {
           double delta = ROOT::Math::VectorUtil::DeltaR(jet.p4().Vect(),it->p4().Vect());
       if(delta<matchingDR)
         {
           matched=true;
           if(delta<delta_min)
         {
           delta_min=delta;
           mceta=it->eta();
           mcet=it->et();
         }
         }
    }
 
   //Create Struct and send it out!
   MatchElementL2 match;
   match.matched=matched;
   match.deltar=delta_min;
   match.mcEta = mceta;
   match.mcEt = mcet;
 
 
  return match;
 }