cms::DiJetAnalyzer Class Reference

#include <DiJetAnalyzer.h>

Inheritance diagram for cms::DiJetAnalyzer:

Public Member Functions
	DiJetAnalyzer (const edm::ParameterSet &)
	~DiJetAnalyzer ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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

Private Attributes
bool	allowMissingInputs_
TClonesArray *	cells
Float_t	emEnergy
Float_t	etVetoJet
UInt_t	eventNumber
std::string	fOutputFileName
TFile *	hOutputFile
Int_t	iEtaHit
UInt_t	iPhiHit
HcalRespCorrs *	oldRespCorrs
Float_t	probeJetEmFrac
TLorentzVector *	probeJetP4
Float_t	PxTrkHcal
Float_t	PyTrkHcal
Float_t	PzTrkHcal
UInt_t	runNumber
Float_t	tagJetEmFrac
TLorentzVector *	tagJetP4
Float_t	targetE
edm::EDGetTokenT < EcalRecHitCollection >	tok_ec_
edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_
edm::EDGetTokenT < HFRecHitCollection >	tok_hf_
edm::EDGetTokenT < HORecHitCollection >	tok_ho_
edm::EDGetTokenT < reco::CaloJetCollection >	tok_jets_
TTree *	tree
Float_t	xTrkHcal
Float_t	yTrkHcal
Float_t	zTrkHcal

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 51 of file DiJetAnalyzer.h.

Constructor & Destructor Documentation

cms::DiJetAnalyzer::DiJetAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 33 of file DiJetAnalyzer.cc.

References allowMissingInputs_, fOutputFileName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), AlCaHLTBitMon_QueryRunRegistry::string, tok_ec_, tok_hbhe_, tok_hf_, tok_ho_, and tok_jets_.

{
  tok_jets_ = consumes<reco::CaloJetCollection>(iConfig.getParameter<edm::InputTag>("jetsInput"));
  tok_ec_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("ecInput"));
  tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
  tok_ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"));
  tok_hf_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInput"));

  // get name of output file with histogramms
  fOutputFileName = iConfig.getUntrackedParameter<std::string>("HistOutFile");

   allowMissingInputs_ = true;
}

cms::DiJetAnalyzer::~DiJetAnalyzer

(

)

Definition at line 48 of file DiJetAnalyzer.cc.

{

}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 60 of file DiJetAnalyzer.cc.

References allowMissingInputs_, HcalDetId::depth(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, alignCSCRings::e, emEnergy, PV3DBase< T, PVType, FrameType >::eta(), etVetoJet, edm::EventID::event(), eventNumber, edm::EventSetup::get(), edm::Event::getByToken(), CaloGeometry::getValidDetIds(), DetId::Hcal, i, edm::EventBase::id(), HcalDetId::ieta(), iEtaHit, HcalDetId::iphi(), iPhiHit, metsig::jet, fwrapper::jets, PV3DBase< T, PVType, FrameType >::phi(), pi, probeJetEmFrac, probeJetP4, PxTrkHcal, PyTrkHcal, PzTrkHcal, dt_dqm_sourceclient_common_cff::reco, edm::EventID::run(), runNumber, funct::sin(), mathSSE::sqrt(), tagJetP4, targetE, tok_ec_, tok_hbhe_, tok_hf_, tok_ho_, tok_jets_, and tree.

{

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

   const double pi = 4.*atan(1.);

  // event number & run number 
   eventNumber = iEvent.id().event(); 
   runNumber = iEvent.id().run();
 

  // read jets
   CaloJet jet1, jet2, jet3; 
   try {
   edm::Handle<CaloJetCollection> jets;
   iEvent.getByToken(tok_jets_,jets);
   if(jets->size()>1){ 
    jet1 = (*jets)[0]; 
    jet2 = (*jets)[1];
     if(fabs(jet1.eta())>fabs(jet2.eta())){
       CaloJet jet = jet1; 
       jet1 = jet2; 
       jet2 = jet; 
     } 
     //     if(fabs(jet1.eta())>eta_1 || (fabs(jet2.eta())-jet_R) < eta_2){ return;}
   } else {return;}
   tagJetP4->SetPxPyPzE(jet1.px(), jet1.py(), jet1.pz(), jet1.energy());
   probeJetP4->SetPxPyPzE(jet2.px(), jet2.py(), jet2.pz(), jet2.energy());
   if(jets->size()>2){
     jet3 = (*jets)[2];
     etVetoJet = jet3.et();
   } else { etVetoJet = 0.;}
   }catch (cms::Exception& e) { // can't find it!
     if (!allowMissingInputs_) { throw e; }  
   }


   double dR = 1000.; 
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   const CaloGeometry* geo = pG.product();
   vector<DetId> vid = geo->getValidDetIds();
   for(vector<DetId>::const_iterator idItr = vid.begin(); idItr != vid.end(); idItr++)
     {
      if( (*idItr).det() == DetId::Hcal ) {
        GlobalPoint pos = geo->getPosition(*idItr); 
        double deta = fabs(jet2.eta() - pos.eta());
        double dphi = fabs(jet2.phi() - pos.phi()); 
        if(dphi>pi){dphi=2*pi-dphi;}
        double dR_candidate = sqrt(deta*deta + dphi*dphi); 
        int ieta = HcalDetId(*idItr).ieta();
        int iphi = HcalDetId(*idItr).iphi();
        int idepth = HcalDetId(*idItr).depth();
        if(dR_candidate < dR && idepth == 1){
      dR = dR_candidate; 
          iEtaHit = ieta; 
          iPhiHit = iphi; 
    }
      }
     }

   targetE = jet1.et()/sin(jet2.theta()); 

   std::map<DetId,int> mapId; 
   vector<CaloTowerPtr> towers_fw = jet2.getCaloConstituents();
   vector<CaloTowerPtr>::const_iterator towersItr;
   for(towersItr = towers_fw.begin(); towersItr != towers_fw.end(); towersItr++){
      size_t tower_size = (*towersItr)->constituentsSize();
      for(size_t i=0; i<tower_size; i++){
        DetId id = (*towersItr)->constituent(i);
        mapId[id] = 1;
      }
   }


   //   probeJetEmFrac = 0.;
   double emEnergy = 0.; 
   try {
      Handle<EcalRecHitCollection> ec;
      iEvent.getByToken(tok_ec_,ec);
      for(EcalRecHitCollection::const_iterator ecItr = (*ec).begin();
                                                ecItr != (*ec).end(); ++ecItr)
      {
        DetId id = ecItr->detid();
        if(mapId[id]==1){
      emEnergy += ecItr->energy(); 
    }
      }
   } catch (cms::Exception& e) { // can't find it!
      if (!allowMissingInputs_) throw e;
   }

   targetE += -emEnergy;
   probeJetEmFrac = emEnergy/jet2.energy(); 

   int nHits = 0; 
   try {
      Handle<HBHERecHitCollection> hbhe;
      iEvent.getByToken(tok_hbhe_, hbhe);
      for(HBHERecHitCollection::const_iterator hbheItr=hbhe->begin(); 
                                                 hbheItr!=hbhe->end(); hbheItr++)
      {
         DetId id = hbheItr->detid();     
         if(mapId[id]==1){
          TCell* cell = new TCell(id.rawId(),hbheItr->energy()); 
          (*cells)[nHits] = cell;
          nHits++;  
         }       
      }
   } catch (cms::Exception& e) { // can't find it!
      if (!allowMissingInputs_) throw e;
   }
 

   try {
      Handle<HORecHitCollection> ho;
      iEvent.getByToken(tok_ho_, ho);
      for(HORecHitCollection::const_iterator hoItr=ho->begin(); 
                                               hoItr!=ho->end(); hoItr++)
      {
         DetId id = hoItr->detid();
         if(mapId[id]==1){
          TCell* cell = new TCell(id.rawId(),hoItr->energy()); 
          (*cells)[nHits] = cell;
          nHits++;  
         }
      }
   } catch (cms::Exception& e) { // can't find it!
      if (!allowMissingInputs_) throw e;
   }

   try {
      Handle<HFRecHitCollection> hf;
      iEvent.getByToken(tok_hf_, hf);
      for(HFRecHitCollection::const_iterator hfItr=hf->begin(); 
                                               hfItr!=hf->end(); hfItr++)
      {
         DetId id = hfItr->detid();
         if(mapId[id]==1){
          TCell* cell = new TCell(id.rawId(),hfItr->energy()); 
          (*cells)[nHits] = cell; 
          nHits++; 
         }
      }
   } catch (cms::Exception& e) { // can't find it!
      if (!allowMissingInputs_) throw e;
   }
 
   probeJetEmFrac = emEnergy/jet2.energy();
   
PxTrkHcal = 0;
PyTrkHcal = 0;
PzTrkHcal = 0;
  
 tree->Fill(); 
 
}

void cms::DiJetAnalyzer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 224 of file DiJetAnalyzer.cc.

References cells, emEnergy, etVetoJet, eventNumber, fOutputFileName, hOutputFile, iEtaHit, iPhiHit, probeJetEmFrac, probeJetP4, PxTrkHcal, PyTrkHcal, PzTrkHcal, runNumber, tagJetEmFrac, tagJetP4, targetE, tree, xTrkHcal, yTrkHcal, and zTrkHcal.

{

  hOutputFile   = new TFile( fOutputFileName.c_str(), "RECREATE" ) ;

    tree = new TTree("hcalCalibTree", "Tree for IsoTrack Calibration");

    cells = new TClonesArray("TCell");
    tagJetP4 = new TLorentzVector();
    probeJetP4 = new TLorentzVector();
 

    tree->Branch("eventNumber", &eventNumber, "eventNumber/i");
    tree->Branch("runNumber", &runNumber, "runNumber/i");  
    tree->Branch("iEtaHit", &iEtaHit, "iEtaHit/I");
    tree->Branch("iPhiHit", &iPhiHit, "iPhiHit/i");    

    
    tree->Branch("xTrkHcal", &xTrkHcal, "xTrkHcal/F");
    tree->Branch("yTrkHcal", &yTrkHcal, "yTrkHcal/F");
    tree->Branch("zTrkHcal", &zTrkHcal, "zTrkHcal/F");

    tree->Branch("PxTrkHcal", &PxTrkHcal, "PxTrkHcal/F");
    tree->Branch("PyTrkHcal", &PyTrkHcal, "PyTrkHcal/F");
    tree->Branch("PzTrkHcal", &PzTrkHcal, "PzTrkHcal/F");

    tree->Branch("cells", &cells, 64000); 
    tree->Branch("emEnergy", &emEnergy, "emEnergy/F"); 
    tree->Branch("targetE", &targetE, "targetE/F");
    tree->Branch("etVetoJet", &etVetoJet, "etVetoJet/F");
    tree->Branch("tagJetP4", "TLorentzVector", &tagJetP4);
    tree->Branch("probeJetP4", "TLorentzVector", &probeJetP4);
    tree->Branch("tagJetEmFrac", &tagJetEmFrac,"tagJetEmFrac/F"); 
    tree->Branch("probeJetEmFrac", &probeJetEmFrac,"probeJetEmFrac/F");    
   
}

void cms::DiJetAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 263 of file DiJetAnalyzer.cc.

References hOutputFile, and tree.

                      {

   hOutputFile->SetCompressionLevel(2);
   hOutputFile->cd();
   tree->Write(); 
   hOutputFile->Close() ;

}

Member Data Documentation

bool cms::DiJetAnalyzer::allowMissingInputs_

private

Definition at line 106 of file DiJetAnalyzer.h.

Referenced by analyze(), and DiJetAnalyzer().

TClonesArray* cms::DiJetAnalyzer::cells

private

Definition at line 94 of file DiJetAnalyzer.h.

Referenced by beginJob().

Float_t cms::DiJetAnalyzer::emEnergy

private

Definition at line 96 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

Float_t cms::DiJetAnalyzer::etVetoJet

private

Definition at line 99 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

UInt_t cms::DiJetAnalyzer::eventNumber

private

Definition at line 81 of file DiJetAnalyzer.h.

Referenced by analyze(), beginJob(), and Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::setFilterBranches().

std::string cms::DiJetAnalyzer::fOutputFileName

private

Definition at line 74 of file DiJetAnalyzer.h.

Referenced by beginJob(), and DiJetAnalyzer().

TFile* cms::DiJetAnalyzer::hOutputFile

private

Definition at line 77 of file DiJetAnalyzer.h.

Referenced by beginJob(), and endJob().

Int_t cms::DiJetAnalyzer::iEtaHit

private

Definition at line 83 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

UInt_t cms::DiJetAnalyzer::iPhiHit

private

Definition at line 84 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

HcalRespCorrs* cms::DiJetAnalyzer::oldRespCorrs

private

Definition at line 108 of file DiJetAnalyzer.h.

Float_t cms::DiJetAnalyzer::probeJetEmFrac

private

Definition at line 103 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

TLorentzVector* cms::DiJetAnalyzer::probeJetP4

private

Definition at line 101 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

Float_t cms::DiJetAnalyzer::PxTrkHcal

private

Definition at line 90 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

Float_t cms::DiJetAnalyzer::PyTrkHcal

private

Definition at line 91 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

Float_t cms::DiJetAnalyzer::PzTrkHcal

private

Definition at line 92 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

UInt_t cms::DiJetAnalyzer::runNumber

private

Definition at line 82 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

Float_t cms::DiJetAnalyzer::tagJetEmFrac

private

Definition at line 102 of file DiJetAnalyzer.h.

Referenced by beginJob().

TLorentzVector* cms::DiJetAnalyzer::tagJetP4

private

Definition at line 100 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

Float_t cms::DiJetAnalyzer::targetE

private

Definition at line 97 of file DiJetAnalyzer.h.

Referenced by analyze(), and beginJob().

edm::EDGetTokenT<EcalRecHitCollection> cms::DiJetAnalyzer::tok_ec_

private

Definition at line 66 of file DiJetAnalyzer.h.

Referenced by analyze(), and DiJetAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> cms::DiJetAnalyzer::tok_hbhe_

private

Definition at line 67 of file DiJetAnalyzer.h.

Referenced by analyze(), and DiJetAnalyzer().

edm::EDGetTokenT<HFRecHitCollection> cms::DiJetAnalyzer::tok_hf_

private

Definition at line 69 of file DiJetAnalyzer.h.

Referenced by analyze(), and DiJetAnalyzer().

edm::EDGetTokenT<HORecHitCollection> cms::DiJetAnalyzer::tok_ho_

private

Definition at line 68 of file DiJetAnalyzer.h.

Referenced by analyze(), and DiJetAnalyzer().

edm::EDGetTokenT<reco::CaloJetCollection> cms::DiJetAnalyzer::tok_jets_

private

Definition at line 65 of file DiJetAnalyzer.h.

Referenced by analyze(), and DiJetAnalyzer().

TTree* cms::DiJetAnalyzer::tree

private

Definition at line 79 of file DiJetAnalyzer.h.

Referenced by svgfig.SVG::__str__(), analyze(), beginJob(), endJob(), and python.cmstools.EventTree::SetAlias().

Float_t cms::DiJetAnalyzer::xTrkHcal

private

Definition at line 86 of file DiJetAnalyzer.h.

Referenced by beginJob().

Float_t cms::DiJetAnalyzer::yTrkHcal

private

Definition at line 87 of file DiJetAnalyzer.h.

Referenced by beginJob().

Float_t cms::DiJetAnalyzer::zTrkHcal

private

Definition at line 88 of file DiJetAnalyzer.h.

Referenced by beginJob().