/data/refman/pasoursint/CMSSW_5_3_10_patch1/src/RecoJets/JetAnalyzers/interface/JetAnaPythia.h

Go to the documentation of this file.

// Template class: JetAnaPythia
// Description:  Example of simple analyzer for jets produced by Pythia
// Author: R. Harris
// Date:  28 - October - 2008
#ifndef JetAnaPythia_h
#define JetAnaPythia_h
#include <TTree.h>
#include <TH1.h>
#include <TFile.h>
#include "TNamed.h"
#include <vector>
#include <map>
#include "FWCore/Framework/interface/EDAnalyzer.h"

template<class Jet>
class JetAnaPythia : public edm::EDAnalyzer 
   {
     public:
       JetAnaPythia(edm::ParameterSet const& cfg);
     private:
       typedef std::vector<Jet> JetCollection;
       void FillHist1D(const TString& histName, const Double_t& x, const Double_t& wt);
       void beginJob();
       void analyze(edm::Event const& e, edm::EventSetup const& iSetup);
       void endJob();
       std::map<TString, TH1*> m_HistNames1D;  
       TTree* mcTruthTree_;
       float xsec;
       float weight;
       float pt_hat;
       int   nJets;
       float etaJet1, etaJet2;
       float ptJet1,  ptJet2;
       float etaPart1, etaPart2;
       float ptPart1,  ptPart2;
       float diJetMass;
       float diPartMass; 
       TFile* m_file;
       std::string JetAlgorithm;
       std::string HistoFileName;
       int NJets;    
       bool debug;
       int eventsGen;
       std::string anaLevel; 
        std::vector<double> xsecGen;
        std::vector<double> ptHatEdges;

   };
#endif