| 1 |
loizides |
1.1 |
//--------------------------------------------------------------------------------------------------
|
| 2 |
loizides |
1.20 |
// $Id: GeneratorMod.h,v 1.19 2009/02/13 14:33:31 loizides Exp $
|
| 3 |
loizides |
1.1 |
//
|
| 4 |
|
|
// GeneratorMod
|
| 5 |
|
|
//
|
| 6 |
loizides |
1.4 |
// This module collects interesting generator information and publishes collections
|
| 7 |
|
|
// for subsequent modules.
|
| 8 |
loizides |
1.1 |
//
|
| 9 |
loizides |
1.4 |
// Authors: G.Ceballos
|
| 10 |
loizides |
1.1 |
//--------------------------------------------------------------------------------------------------
|
| 11 |
|
|
|
| 12 |
loizides |
1.4 |
#ifndef MITPHYSICS_MODS_GENERATORMOD_H
|
| 13 |
|
|
#define MITPHYSICS_MODS_GENERATORMOD_H
|
| 14 |
loizides |
1.1 |
|
| 15 |
|
|
#include "MitAna/TreeMod/interface/BaseMod.h"
|
| 16 |
|
|
#include "MitAna/DataTree/interface/Collections.h"
|
| 17 |
|
|
|
| 18 |
|
|
class TH1D;
|
| 19 |
|
|
class TH2D;
|
| 20 |
|
|
|
| 21 |
|
|
namespace mithep
|
| 22 |
|
|
{
|
| 23 |
|
|
class GeneratorMod : public BaseMod
|
| 24 |
|
|
{
|
| 25 |
|
|
public:
|
| 26 |
|
|
GeneratorMod(const char *name="GeneratorMod",
|
| 27 |
loizides |
1.4 |
const char *title="Generator information module");
|
| 28 |
|
|
|
| 29 |
loizides |
1.12 |
const char *GetMCPartName() const { return fMCPartName; }
|
| 30 |
loizides |
1.19 |
const char *GetMCLeptonsName() const { return fMCLeptonsName; }
|
| 31 |
loizides |
1.12 |
const char *GetMCAllLeptonsName() const { return fMCAllLeptonsName; }
|
| 32 |
|
|
const char *GetMCTausName() const { return fMCTausName; }
|
| 33 |
|
|
const char *GetMCNeutrinosName() const { return fMCNeutrinosName; }
|
| 34 |
|
|
const char *GetMCQuarksName() const { return fMCQuarksName; }
|
| 35 |
|
|
const char *GetMCqqHsName() const { return fMCqqHsName; }
|
| 36 |
|
|
const char *GetMCBosonsName() const { return fMCBosonsName; }
|
| 37 |
|
|
const char *GetMCPhotonsName() const { return fMCPhotonsName; }
|
| 38 |
|
|
void SetMCPartName(const char *s) { fMCPartName = s; }
|
| 39 |
|
|
void SetMCLeptonsName(const char * s) { fMCLeptonsName = s; }
|
| 40 |
|
|
void SetMCAllLeptonsName(const char * s) { fMCAllLeptonsName = s; }
|
| 41 |
|
|
void SetMCTausName(const char *s) { fMCTausName = s; }
|
| 42 |
|
|
void SetMCNeutrinosName(const char *s) { fMCNeutrinosName = s; }
|
| 43 |
|
|
void SetMCQuarksName(const char *s) { fMCQuarksName = s; }
|
| 44 |
|
|
void SetMCqqHsName(const char *s) { fMCqqHsName = s; }
|
| 45 |
|
|
void SetMCBosonsName(const char *s) { fMCBosonsName = s; }
|
| 46 |
|
|
void SetMCPhotonsName(const char *s) { fMCPhotonsName = s; }
|
| 47 |
|
|
void SetPtLeptonMin(Double_t x) { fPtLeptonMin = x; }
|
| 48 |
|
|
void SetEtaLeptonMax(Double_t x) { fEtaLeptonMax = x; }
|
| 49 |
|
|
void SetPtPhotonMin(Double_t x) { fPtPhotonMin = x; }
|
| 50 |
|
|
void SetEtaPhotonMax(Double_t x) { fEtaPhotonMax = x; }
|
| 51 |
loizides |
1.19 |
void SetPdgIdCut(UInt_t d) { fPdgIdCut = d; }
|
| 52 |
ceballos |
1.18 |
void SetMassMinCut(Double_t x) { fMassMinCut = x; }
|
| 53 |
|
|
void SetMassMaxCut(Double_t x) { fMassMaxCut = x; }
|
| 54 |
loizides |
1.4 |
|
| 55 |
loizides |
1.1 |
protected:
|
| 56 |
loizides |
1.14 |
void Process();
|
| 57 |
|
|
void SlaveBegin();
|
| 58 |
|
|
|
| 59 |
loizides |
1.12 |
Bool_t fFillHist; //=true then fill histos (def=0)
|
| 60 |
|
|
TString fMCPartName; //name of MCParticle branch
|
| 61 |
|
|
TString fMCLeptonsName; //name of lepton coll (from W)
|
| 62 |
|
|
TString fMCAllLeptonsName; //name of lepton coll (all)
|
| 63 |
|
|
TString fMCTausName; //name of tau coll (hadronic decays)
|
| 64 |
|
|
TString fMCNeutrinosName; //name of neutrinos coll
|
| 65 |
|
|
TString fMCQuarksName; //name of quarks coll
|
| 66 |
|
|
TString fMCqqHsName; //name of qqH coll
|
| 67 |
|
|
TString fMCBosonsName; //name of bosons coll
|
| 68 |
|
|
TString fMCPhotonsName; //name of photons coll
|
| 69 |
|
|
Double_t fPtLeptonMin; //pt min for leptons
|
| 70 |
|
|
Double_t fEtaLeptonMax; //eta max for leptons
|
| 71 |
|
|
Double_t fPtPhotonMin; //pt min for photons
|
| 72 |
|
|
Double_t fEtaPhotonMax; //eta max for photons
|
| 73 |
loizides |
1.19 |
UInt_t fPdgIdCut; //pdg id for particle used to select on mass (0=off)
|
| 74 |
|
|
Double_t fMassMinCut; //mass min for given PdgId particle
|
| 75 |
|
|
Double_t fMassMaxCut; //mass max for given PdgId particle
|
| 76 |
loizides |
1.12 |
const MCParticleCol *fParticles; //!MCParticle branch
|
| 77 |
ceballos |
1.18 |
|
| 78 |
ceballos |
1.17 |
TH1D *hDGenLeptons[40]; //!histos for W leptons
|
| 79 |
loizides |
1.12 |
TH1D *hDGenAllLeptons[20]; //!histos for all leptons
|
| 80 |
|
|
TH1D *hDGenTaus[20]; //!histos for taus
|
| 81 |
|
|
TH1D *hDGenNeutrinos[20]; //!histos for neutrinos
|
| 82 |
|
|
TH1D *hDGenQuarks[20]; //!histos for quarks
|
| 83 |
|
|
TH1D *hDGenWBF[20]; //!histos for WBF
|
| 84 |
|
|
TH1D *hDGenBosons[20]; //!histos for bosons
|
| 85 |
|
|
TH1D *hDGenPhotons[20]; //!histos for photons
|
| 86 |
loizides |
1.13 |
TH1D *hDVMass[20]; //!histos for auxiliar work
|
| 87 |
ceballos |
1.17 |
|
| 88 |
loizides |
1.16 |
ClassDef(GeneratorMod, 1) // Module to gather generator information
|
| 89 |
loizides |
1.1 |
};
|
| 90 |
|
|
}
|
| 91 |
|
|
#endif
|
