| 1 |
//--------------------------------------------------------------------------------------------------
|
| 2 |
// $Id: GeneratorMod.h,v 1.35 2009/12/06 14:59:43 ceballos Exp $
|
| 3 |
//
|
| 4 |
// GeneratorMod
|
| 5 |
//
|
| 6 |
// This module collects interesting generator information and publishes collections
|
| 7 |
// for subsequent modules.
|
| 8 |
//
|
| 9 |
// Authors: G.Ceballos
|
| 10 |
//--------------------------------------------------------------------------------------------------
|
| 11 |
|
| 12 |
#ifndef MITPHYSICS_MODS_GENERATORMOD_H
|
| 13 |
#define MITPHYSICS_MODS_GENERATORMOD_H
|
| 14 |
|
| 15 |
#include "MitAna/TreeMod/interface/BaseMod.h"
|
| 16 |
#include "MitAna/DataTree/interface/MCParticleFwd.h"
|
| 17 |
#include "MitAna/DataTree/interface/MetFwd.h"
|
| 18 |
|
| 19 |
class TH1D;
|
| 20 |
class TH2D;
|
| 21 |
|
| 22 |
namespace mithep
|
| 23 |
{
|
| 24 |
class GeneratorMod : public BaseMod
|
| 25 |
{
|
| 26 |
public:
|
| 27 |
GeneratorMod(const char *name="GeneratorMod",
|
| 28 |
const char *title="Generator information module");
|
| 29 |
~GeneratorMod();
|
| 30 |
|
| 31 |
Bool_t GetIsData() const { return fIsData; }
|
| 32 |
Bool_t GetApplyISRFilter() const { return fApplyISRFilter; }
|
| 33 |
Bool_t GetCopyArrays() const { return fCopyArrays; }
|
| 34 |
const char *GetMCAllLeptonsName() const { return fMCAllLeptonsName; }
|
| 35 |
const char *GetMCBosonsName() const { return fMCBosonsName; }
|
| 36 |
const char *GetMCISRPhotonsName() const { return fMCISRPhotonsName; }
|
| 37 |
const char *GetMCLeptonsName() const { return fMCLeptonsName; }
|
| 38 |
const char *GetMCMETName() const { return fMCMETName; }
|
| 39 |
const char *GetMCNeutrinosName() const { return fMCNeutrinosName; }
|
| 40 |
const char *GetMCPartName() const { return fMCPartName; }
|
| 41 |
const char *GetMCPhotonsName() const { return fMCPhotonsName; }
|
| 42 |
const char *GetMCQuarksName() const { return fMCQuarksName; }
|
| 43 |
const char *GetMCRadPhotonsName() const { return fMCRadPhotonsName; }
|
| 44 |
const char *GetMCTausName() const { return fMCTausName; }
|
| 45 |
const char *GetMCqqHsName() const { return fMCqqHsName; }
|
| 46 |
Bool_t GetPrintDebug() const { return fPrintDebug; }
|
| 47 |
void SetIsData(Bool_t b) { fIsData = b; }
|
| 48 |
void SetApplyISRFilter(Bool_t b) { fApplyISRFilter = b; }
|
| 49 |
void SetCopyArrays(Bool_t b) { fCopyArrays = b; }
|
| 50 |
void SetEtaLeptonMax(Double_t x) { fEtaLeptonMax = x; }
|
| 51 |
void SetEtaPhotonMax(Double_t x) { fEtaPhotonMax = x; }
|
| 52 |
void SetEtaRadPhotonMax(Double_t x) { fEtaRadPhotonMax = x; }
|
| 53 |
void SetMCAllLeptonsName(const char * s) { fMCAllLeptonsName = s; }
|
| 54 |
void SetMCBosonsName(const char *s) { fMCBosonsName = s; }
|
| 55 |
void SetMCISRPhotonsName(const char *s) { fMCISRPhotonsName = s; }
|
| 56 |
void SetMCLeptonsName(const char * s) { fMCLeptonsName = s; }
|
| 57 |
void SetMCMETName(const char * s) { fMCMETName = s; }
|
| 58 |
void SetMCNeutrinosName(const char *s) { fMCNeutrinosName = s; }
|
| 59 |
void SetMCPartName(const char *s) { fMCPartName = s; }
|
| 60 |
void SetMCPhotonsName(const char *s) { fMCPhotonsName = s; }
|
| 61 |
void SetMCQuarksName(const char *s) { fMCQuarksName = s; }
|
| 62 |
void SetMCRadPhotonsName(const char *s) { fMCRadPhotonsName = s; }
|
| 63 |
void SetMCTausName(const char *s) { fMCTausName = s; }
|
| 64 |
void SetMCqqHsName(const char *s) { fMCqqHsName = s; }
|
| 65 |
void SetMassMaxCut(Double_t x) { fMassMaxCut = x; }
|
| 66 |
void SetMassMinCut(Double_t x) { fMassMinCut = x; }
|
| 67 |
void SetPdgIdCut(UInt_t d) { fPdgIdCut = d; }
|
| 68 |
void SetPrintDebug(bool b) { fPrintDebug = b; }
|
| 69 |
void SetPtLeptonMin(Double_t x) { fPtLeptonMin = x; }
|
| 70 |
void SetPtPhotonMin(Double_t x) { fPtPhotonMin = x; }
|
| 71 |
void SetPtRadPhotonMin(Double_t x) { fPtRadPhotonMin = x; }
|
| 72 |
|
| 73 |
protected:
|
| 74 |
void Process();
|
| 75 |
void SlaveBegin();
|
| 76 |
|
| 77 |
Bool_t fIsData; //=true then it does nothing (def=0)
|
| 78 |
Bool_t fPrintDebug; //=true then print debug info (def=0)
|
| 79 |
Bool_t fCopyArrays; //=true then copy array content for skimming (def=0)
|
| 80 |
TString fMCPartName; //name of MCParticle branch
|
| 81 |
TString fMCMETName; //name of met coll
|
| 82 |
TString fMCLeptonsName; //name of lepton coll (from W/Z/H)
|
| 83 |
TString fMCAllLeptonsName; //name of lepton coll (all)
|
| 84 |
TString fMCTausName; //name of tau coll (hadronic decays)
|
| 85 |
TString fMCNeutrinosName; //name of neutrinos coll
|
| 86 |
TString fMCQuarksName; //name of quarks coll
|
| 87 |
TString fMCqqHsName; //name of qqH coll
|
| 88 |
TString fMCBosonsName; //name of bosons coll
|
| 89 |
TString fMCPhotonsName; //name of photons coll
|
| 90 |
TString fMCRadPhotonsName; //name of rad photons coll
|
| 91 |
TString fMCISRPhotonsName; //name of ISR photons coll
|
| 92 |
Double_t fPtLeptonMin; //pt min for leptons
|
| 93 |
Double_t fEtaLeptonMax; //eta max for leptons
|
| 94 |
Double_t fPtPhotonMin; //pt min for photons
|
| 95 |
Double_t fEtaPhotonMax; //eta max for photons
|
| 96 |
Double_t fPtRadPhotonMin; //pt min for rad photons
|
| 97 |
Double_t fEtaRadPhotonMax; //eta max for rad photons
|
| 98 |
UInt_t fPdgIdCut; //pdg id for particle used to select on mass (def=0)
|
| 99 |
Double_t fMassMinCut; //mass min for given PdgId particle
|
| 100 |
Double_t fMassMaxCut; //mass max for given PdgId particle
|
| 101 |
Bool_t fApplyISRFilter; //=true then apply ISR filter (def=0)
|
| 102 |
const MCParticleCol *fParticles; //!MCParticle branch
|
| 103 |
TH1D *hDGenPtMin; //!histo for Pt min for leptons from W/Z
|
| 104 |
TH1D *hDGenMet[10]; //!histos for gen MET
|
| 105 |
TH1D *hDGenLeptons[40]; //!histos for W/Z/H leptons
|
| 106 |
TH1D *hDGenAllLeptons[20]; //!histos for all leptons
|
| 107 |
TH1D *hDGenTaus[20]; //!histos for taus
|
| 108 |
TH1D *hDGenNeutrinos[20]; //!histos for neutrinos
|
| 109 |
TH1D *hDGenQuarks[20]; //!histos for quarks
|
| 110 |
TH1D *hDGenWBF[20]; //!histos for WBF
|
| 111 |
TH1D *hDGenBosons[20]; //!histos for bosons
|
| 112 |
TH1D *hDGenPhotons[20]; //!histos for photons
|
| 113 |
TH1D *hDGenRadPhotons[20]; //!histos for rad photons
|
| 114 |
TH1D *hDGenISRPhotons[20]; //!histos for ISR photons
|
| 115 |
TH1D *hDVMass[20]; //!histos for auxiliar MG work
|
| 116 |
TH1D *hDVVMass[50]; //!histos for auxiliar VV work
|
| 117 |
MCParticleArr *fGenLeptons; //!copied owning array for skimming
|
| 118 |
MCParticleArr *fGenAllLeptons; //!copied owning array for skimming
|
| 119 |
MCParticleArr *fGenTaus; //!copied owning array for skimming
|
| 120 |
MCParticleArr *fGenNeutrinos; //!copied owning array for skimming
|
| 121 |
MCParticleArr *fGenQuarks; //!copied owning array for skimming
|
| 122 |
MCParticleArr *fGenqqHs; //!copied owning array for skimming
|
| 123 |
MCParticleArr *fGenBosons; //!copied owning array for skimming
|
| 124 |
MCParticleArr *fGenPhotons; //!copied owning array for skimming
|
| 125 |
MCParticleArr *fGenRadPhotons; //!copied owning array for skimming
|
| 126 |
MCParticleArr *fGenISRPhotons; //!copied owning array for skimming
|
| 127 |
|
| 128 |
ClassDef(GeneratorMod, 1) // Module to gather generator information
|
| 129 |
};
|
| 130 |
}
|
| 131 |
#endif
|
