DPNtupleMakerAndAnalyzer/GenStudy.h at master · TENorbert/DPNtupleMakerAndAnalyzer · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#ifndef GenStudy_H
#define GenStudy_H

// -*- C++ -*-
//
// Package:    GenStudy
// Class:      GenStudy
//
/**\class GenStudy GenStudy.cc EXO/GenStudy/src/GenStudy.cc

 Description: [one line class summary]

 Implementation:
     [Notes on implementation]
*/
//
// Original Author:  Shih-Chuan Kao
//         Created:  Thu Mar 15 12:51:07 CDT 2012
// $Id: GenStudy.h,v 1.15 2012/12/31 22:45:14 sckao Exp $
//
//


// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/Candidate/interface/Particle.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "HepMC/GenParticle.h"
#include "HepMC/SimpleVector.h"

#include <TMath.h>
#include <Math/VectorUtil.h>
#include <TRandom.h>
#include "TFile.h"
#include "TTree.h"
#include "TLorentzVector.h"

#include "Ntuple.h"
//
// class declaration
//
using namespace std ;
typedef math::XYZTLorentzVector LorentzVector;

//class GenStudy : public edm::EDAnalyzer {
class GenStudy {
   public:
      explicit GenStudy(const edm::ParameterSet& );
      ~GenStudy();

      void GetGen(const edm::Event&, Ntuple& leaves );
      void GetGenEvent(const edm::Event&, Ntuple& leaves, bool debug = false );
      void PrintGenEvent(const edm::Event& ) ;

   private:

      bool Propagator( LorentzVector v, double& x, double& y, double& z, double& t , double ctaugamma = 999999. );

      TRandom* expPDF ;

      // ----------member data ---------------------------
      double tau ;
      edm::InputTag genSrc ;

};

#endif