Mods/interface/HKFactorProducer.h

//--------------------------------------------------------------------------------------------------
// $Id: HKFactorProducer.h,v 1.6 2011/06/10 10:42:36 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; }        
      void               SetDoHiggsPtReweighting(Bool_t b) { fDoHiggsPtReweighting = b; }        
      void               SetOutputName(const char *f)      { fOutputName           = f; }

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

      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)
      Bool_t             fDoHiggsPtReweighting; //=true then it does reweighting
      HWWKfactorList    *fPt_histo;             //!histogram with weights read from input file
      const MCEventInfo *fMCEventInfo;          //!event info branch pointer
      TH1D              *hDHKFactor[10];        //!output histograms
      TFile             *fOutputFile;           //output file handle
      TString            fOutputName;           //output file name

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

    ClassDef(HKFactorProducer, 1) // Module to produce k factors
  };
}
#endif
Revision:	1.7
Committed:	Tue Jul 26 15:01:20 2011 UTC (13 years, 9 months ago) by sixie
Content type:	text/plain
Branch:	MAIN
CVS Tags:	Mit_025c, Mit_025b, Mit_025a, Mit_025, Mit_025pre2, Mit_024b, Mit_025pre1, Mit_024a, Mit_024
Changes since 1.6:	+23 -21 lines
Log Message:	add bool option to turn off loading kfactor
#	User	Rev	Content
1	loizides	1.2	//--------------------------------------------------------------------------------------------------
2	sixie	1.7	// $Id: HKFactorProducer.h,v 1.6 2011/06/10 10:42:36 ceballos Exp $
3	ceballos	1.1	//
4			// HKFactorProducer
5			//
6	loizides	1.2	// Produces the k factors from LO to NNLO.
7	ceballos	1.1	//
8			// Authors: G. Gomez-Ceballos
9	loizides	1.2	//--------------------------------------------------------------------------------------------------
10	ceballos	1.1
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	sixie	1.7	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			void SetDoHiggsPtReweighting(Bool_t b) { fDoHiggsPtReweighting = b; }
37			void SetOutputName(const char *f) { fOutputName = f; }
38	ceballos	1.1
39			protected:
40	loizides	1.2	void Process();
41			void SlaveBegin();
42	ceballos	1.6	void SlaveTerminate();
43	ceballos	1.1
44	sixie	1.7	Int_t fProcessID; //process id (from pythia)
45			TString fInputFileName; //input file name
46			TString fMCBosonsName; //boson collection input name
47			TString fMCEvInfoName; //event info branch name
48			Bool_t fIsData; //=true then it does nothing (def=0)
49			Bool_t fMakePDFNtuple; //=true then it does nothing (def=0)
50			Bool_t fDoHiggsPtReweighting; //=true then it does reweighting
51			HWWKfactorList *fPt_histo; //!histogram with weights read from input file
52			const MCEventInfo *fMCEventInfo; //!event info branch pointer
53			TH1D *hDHKFactor[10]; //!output histograms
54			TFile *fOutputFile; //output file handle
55			TString fOutputName; //output file name
56	ceballos	1.1
57	sixie	1.7	float fTreeVariables[8]; //Ntuple variables
58			TTree *fTree; //ntuple tree
59	ceballos	1.5
60	loizides	1.2	ClassDef(HKFactorProducer, 1) // Module to produce k factors
61	ceballos	1.1	};
62			}
63			#endif