/data/refman/pasoursint/CMSSW_5_3_9_patch3/src/Validation/EventGenerator/interface/TauValidation.h

Go to the documentation of this file.

#ifndef TauValidation_H
#define TauValidation_H

/*class TauValidation
 *  
 *  Class to fill Event Generator dqm monitor elements; works on HepMCProduct
 *
 *  $Date: 2012/02/10 14:51:48 $
 *  $Revision: 1.10 $
 *
 */

// framework & common header files
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Run.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"

//DQM services
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"

#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"
#include "TLorentzVector.h"

#include "Validation/EventGenerator/interface/WeightManager.h"

class TauValidation : public edm::EDAnalyzer
{
    public:
        // tau decays
        enum  {undetermined,
               electron,
               muon,
               pi,
               rho,
               a1,
               K,
               Kstar,
               pi1pi0,
               pinpi0,
               tripi,
//               tripinpi0,
               stable};
        // tau mother particles 
        enum  {other,
               B,
               D,
               gamma,
               Z,
               W,
               HSM,
               H0,
               A0,
               Hpm};

    public:
        explicit TauValidation(const edm::ParameterSet&);
        virtual ~TauValidation();
        virtual void beginJob();
        virtual void endJob();  
        virtual void analyze(const edm::Event&, const edm::EventSetup&);
        virtual void beginRun(const edm::Run&, const edm::EventSetup&);
        virtual void endRun(const edm::Run&, const edm::EventSetup&);

    private:
        int tauMother(const HepMC::GenParticle*, double weight);
        int tauProngs(const HepMC::GenParticle*, double weight);
        int tauDecayChannel(const HepMC::GenParticle*, double weight);
        int findMother(const HepMC::GenParticle*);
        int findTauDecayChannel(const HepMC::GenParticle*);
        void rtau(const HepMC::GenParticle*,int,int, double weight);
        void spinEffects(const HepMC::GenParticle*,int,int, double weight);
        void spinEffectsZ(const HepMC::GenParticle*, double weight);
        double leadingPionMomentum(const HepMC::GenParticle*, double weight);
        double visibleTauEnergy(const HepMC::GenParticle*);
        TLorentzVector leadingPionP4(const HepMC::GenParticle*);
        TLorentzVector motherP4(const HepMC::GenParticle*);
        void photons(const HepMC::GenParticle*, double weight);

        WeightManager _wmanager;

        edm::InputTag hepmcCollection_;

        double tauEtCut;

        edm::ESHandle<HepPDT::ParticleDataTable> fPDGTable ;
  
        DQMStore *dbe;

        MonitorElement *nEvt;
        MonitorElement *TauPt, *TauEta, *TauPhi, *TauProngs, *TauDecayChannels, *TauMothers, 
                       *TauRtauW, *TauRtauHpm,
                       *TauSpinEffectsW, *TauSpinEffectsHpm, *TauSpinEffectsZ,
                       *TauPhotonsN,*TauPhotonsPt;
        unsigned int NJAKID;
        MonitorElement *JAKID;
        std::vector<std::vector<MonitorElement *> > JAKInvMass;

};

#endif