Mods/interface/HKFactorProducer.h

//--------------------------------------------------------------------------------------------------
// $Id: HKFactorProducer.h,v 1.4 2009/12/07 11:39:50 ceballos Exp $
//
// HKFactorProducer
//
// Produces the k factors from LO to NNLO.
//
// Authors: G. Gomez-Ceballos
//--------------------------------------------------------------------------------------------------

#ifndef MITPHYSICS_MODS_HKFACTORPRODUCER_H
#define MITPHYSICS_MODS_HKFACTORPRODUCER_H

#include "MitAna/TreeMod/interface/BaseMod.h" 
#include "MitAna/DataTree/interface/MCEventInfo.h"
#include "MitPhysics/Mods/interface/HWWKFactorList.h"

class TH1D;
class TH2D;

namespace mithep 
{
  class HKFactorProducer : public BaseMod
  {
    public:
      HKFactorProducer(const char *name="HKFactorProducer", 
                   const char *title="KFactor information module");
      ~HKFactorProducer();

      void               SetProcessID(Int_t d)             { fProcessID     = d; }
      void               SetInputFilename(const char *s)   { fInputFileName = s; }
      void               SetMCBosonsName(const char *s)    { fMCBosonsName  = s; }
      void               SetMCEventInfoName(const char *s) { fMCEvInfoName  = s; }
      void               SetIsData(Bool_t b)               { fIsData        = b; }       
      void               SetMakePDFNtuple(Bool_t b)        { fMakePDFNtuple = b; }       

    protected:
      void               Process();
      void               SlaveBegin();

      Int_t              fProcessID;        //process id (from pythia)
      TString            fInputFileName;    //input file name
      TString            fMCBosonsName;     //boson collection input name
      TString            fMCEvInfoName;     //event info branch name
      Bool_t             fIsData;           //=true then it does nothing (def=0)
      Bool_t             fMakePDFNtuple;    //=true then it does nothing (def=0)
      HWWKfactorList    *fPt_histo;         //!histogram with weights read from input file
      const MCEventInfo *fMCEventInfo;      //!event info branch pointer
      TH1D              *hDHKFactor[10];    //!output histograms

      float             fTreeVariables[8];            //Ntuple variables
      TTree            *fTree;                        //ntuple tree

    ClassDef(HKFactorProducer, 1) // Module to produce k factors
  };
}
#endif
Revision:	1.5
Committed:	Tue Oct 19 22:25:25 2010 UTC (14 years, 6 months ago) by ceballos
Content type:	text/plain
Branch:	MAIN
CVS Tags:	Mit_020d, TMit_020d, Mit_020c, Mit_021, Mit_021pre2, Mit_021pre1, Mit_020b, Mit_020a, Mit_020, Mit_020pre1, Mit_018, Mit_017, Mit_017pre3, Mit_017pre2, Mit_017pre1, Mit_016, Mit_015b, Mit_015a
Changes since 1.4:	+6 -1 lines
Log Message:	new stuff for PDFs
#	Content
1	//--------------------------------------------------------------------------------------------------
2	// $Id: HKFactorProducer.h,v 1.4 2009/12/07 11:39:50 ceballos Exp $
3	//
4	// HKFactorProducer
5	//
6	// Produces the k factors from LO to NNLO.
7	//
8	// Authors: G. Gomez-Ceballos
9	//--------------------------------------------------------------------------------------------------
10
11	#ifndef MITPHYSICS_MODS_HKFACTORPRODUCER_H
12	#define MITPHYSICS_MODS_HKFACTORPRODUCER_H
13
14	#include "MitAna/TreeMod/interface/BaseMod.h"
15	#include "MitAna/DataTree/interface/MCEventInfo.h"
16	#include "MitPhysics/Mods/interface/HWWKFactorList.h"
17
18	class TH1D;
19	class TH2D;
20
21	namespace mithep
22	{
23	class HKFactorProducer : public BaseMod
24	{
25	public:
26	HKFactorProducer(const char *name="HKFactorProducer",
27	const char *title="KFactor information module");
28	~HKFactorProducer();
29
30	void SetProcessID(Int_t d) { fProcessID = d; }
31	void SetInputFilename(const char *s) { fInputFileName = s; }
32	void SetMCBosonsName(const char *s) { fMCBosonsName = s; }
33	void SetMCEventInfoName(const char *s) { fMCEvInfoName = s; }
34	void SetIsData(Bool_t b) { fIsData = b; }
35	void SetMakePDFNtuple(Bool_t b) { fMakePDFNtuple = b; }
36
37	protected:
38	void Process();
39	void SlaveBegin();
40
41	Int_t fProcessID; //process id (from pythia)
42	TString fInputFileName; //input file name
43	TString fMCBosonsName; //boson collection input name
44	TString fMCEvInfoName; //event info branch name
45	Bool_t fIsData; //=true then it does nothing (def=0)
46	Bool_t fMakePDFNtuple; //=true then it does nothing (def=0)
47	HWWKfactorList *fPt_histo; //!histogram with weights read from input file
48	const MCEventInfo *fMCEventInfo; //!event info branch pointer
49	TH1D *hDHKFactor[10]; //!output histograms
50
51	float fTreeVariables[8]; //Ntuple variables
52	TTree *fTree; //ntuple tree
53
54	ClassDef(HKFactorProducer, 1) // Module to produce k factors
55	};
56	}
57	#endif