CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
myFastSimVal Class Reference

#include <myFastSimVal.h>

Inheritance diagram for myFastSimVal:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 myFastSimVal (const edm::ParameterSet &)
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void beginJob () override
 
void endJob () override
 

Private Attributes

std::string CaloJetAlgorithm1
 
std::string CaloJetAlgorithm2
 
std::string CaloJetAlgorithm3
 
std::string CaloJetAlgorithm4
 
TH1F dEta1
 
TH1F dEta12
 
TH1F dEta2
 
TH1F dEta3
 
TH1F dEta4
 
TH1F dEtaPar1
 
TH1F dEtaPar2
 
TH1F dEtaPar3
 
TH1F dEtaPar4
 
TH1F dijetMass1
 
TH1F dijetMass101
 
TH1F dijetMass102
 
TH1F dijetMass103
 
TH1F dijetMass12
 
TH1F dijetMass13
 
TH1F dijetMass131
 
TH1F dijetMass132
 
TH1F dijetMass133
 
TH1F dijetMass2
 
TH1F dijetMass22
 
TH1F dijetMass23
 
TH1F dijetMass3
 
TH1F dijetMass32
 
TH1F dijetMass33
 
TH1F dijetMass4
 
TH1F dijetMass401
 
TH1F dijetMass402
 
TH1F dijetMass403
 
TH1F dijetMass42
 
TH1F dijetMass43
 
TH1F dijetMass_2000_101
 
TH1F dijetMass_2000_131
 
TH1F dijetMass_2000_401
 
TH1F dijetMass_5000_101
 
TH1F dijetMass_5000_131
 
TH1F dijetMass_5000_401
 
TH1F dijetMass_700_101
 
TH1F dijetMass_700_131
 
TH1F dijetMass_700_401
 
TH1F dijetMassCor1
 
TH1F dijetMassCor101
 
TH1F dijetMassCor131
 
TH1F dijetMassCor401
 
TH1F dijetMassCor_2000_1
 
TH1F dijetMassCor_2000_101
 
TH1F dijetMassCor_2000_131
 
TH1F dijetMassCor_2000_401
 
TH1F dijetMassCor_5000_1
 
TH1F dijetMassCor_5000_101
 
TH1F dijetMassCor_5000_131
 
TH1F dijetMassCor_5000_401
 
TH1F dijetMassCor_700_1
 
TH1F dijetMassCor_700_101
 
TH1F dijetMassCor_700_131
 
TH1F dijetMassCor_700_401
 
TH1F dijetMassP1
 
TH1F dijetMassP101
 
TH1F dijetMassP131
 
TH1F dijetMassP2
 
TH1F dijetMassP3
 
TH1F dijetMassP401
 
TH1F dijetMassP_2000_101
 
TH1F dijetMassP_2000_131
 
TH1F dijetMassP_2000_401
 
TH1F dijetMassP_5000_101
 
TH1F dijetMassP_5000_131
 
TH1F dijetMassP_5000_401
 
TH1F dijetMassP_700_101
 
TH1F dijetMassP_700_131
 
TH1F dijetMassP_700_401
 
TH1F dPhi1
 
TH1F dPhi12
 
TH1F dPhi2
 
TH1F dPhi3
 
TH1F dPhi4
 
TH1F dPhiPar1
 
TH1F dPhiPar2
 
TH1F dPhiPar3
 
TH1F dPhiPar4
 
TH1F dPt1
 
TH1F dPt100Frac1
 
TH1F dPt100Frac2
 
TH1F dPt100Frac3
 
TH1F dPt100Frac4
 
TH1F dPt12
 
TH1F dPt2
 
TH1F dPt20Frac1
 
TH1F dPt20Frac2
 
TH1F dPt20Frac3
 
TH1F dPt20Frac4
 
TH1F dPt3
 
TH1F dPt4
 
TH1F dPt40Frac1
 
TH1F dPt40Frac2
 
TH1F dPt40Frac3
 
TH1F dPt40Frac4
 
TH1F dPt80Frac1
 
TH1F dPt80Frac2
 
TH1F dPt80Frac3
 
TH1F dPt80Frac4
 
TH1F dPtFrac1
 
TH1F dPtFrac2
 
TH1F dPtFrac3
 
TH1F dPtFrac4
 
TH1F dPtPar1
 
TH1F dPtPar2
 
TH1F dPtPar3
 
TH1F dPtPar4
 
TH1F dR1
 
TH1F dR12
 
TH1F dR2
 
TH1F dR3
 
TH1F dR4
 
TH1F dRPar1
 
TH1F dRPar2
 
TH1F dRPar3
 
TH1F dRPar4
 
TH1F dRParton
 
TH1F dRPartonMin
 
TH1F emEneLeadJet1
 
TH1F emEneLeadJet12
 
TH1F emEneLeadJet13
 
TH1F emEneLeadJet2
 
TH1F emEneLeadJet22
 
TH1F emEneLeadJet23
 
TH1F emEneLeadJet3
 
TH1F emEneLeadJet32
 
TH1F emEneLeadJet33
 
TH1F emEneLeadJetEta1_1
 
TH1F emEneLeadJetEta1_2
 
TH1F emEneLeadJetEta2_1
 
TH1F emEneLeadJetEta2_2
 
TH1F emEneLeadJetEta3_1
 
TH1F emEneLeadJetEta3_2
 
std::string GenJetAlgorithm1
 
std::string GenJetAlgorithm2
 
std::string GenJetAlgorithm3
 
std::string GenJetAlgorithm4
 
TH1F h_etaCal1
 
TH1F h_etaCal2
 
TH1F h_etaCal3
 
TH1F h_etaCal4
 
TH1F h_etaCalL1
 
TH1F h_etaCalL2
 
TH1F h_etaCalL3
 
TH1F h_etaCalL4
 
TH1F h_etaGen1
 
TH1F h_etaGen2
 
TH1F h_etaGen3
 
TH1F h_etaGen4
 
TH1F h_etaGenL1
 
TH1F h_etaGenL2
 
TH1F h_etaGenL3
 
TH1F h_jet1Pt1
 
TH1F h_jet1Pt2
 
TH1F h_jet1Pt3
 
TH1F h_jet1Pt4
 
TH1F h_jet2Pt1
 
TH1F h_jet2Pt2
 
TH1F h_jet2Pt3
 
TH1F h_jet2Pt4
 
TH1F h_jet3Pt1
 
TH1F h_jet3Pt2
 
TH1F h_jet3Pt3
 
TH1F h_jet3Pt4
 
TH1F h_jet4Pt1
 
TH1F h_jet4Pt2
 
TH1F h_jet4Pt3
 
TH1F h_jet4Pt4
 
TH1F h_jet5Pt1
 
TH1F h_jet5Pt2
 
TH1F h_jet5Pt3
 
TH1F h_jet5Pt4
 
TH1F h_jet6Pt1
 
TH1F h_jet6Pt2
 
TH1F h_jet6Pt3
 
TH1F h_jet6Pt4
 
TH1F h_jet7Pt1
 
TH1F h_jet7Pt2
 
TH1F h_jet7Pt3
 
TH1F h_jet7Pt4
 
TH1F h_jetEt1
 
TH1F h_jetEt2
 
TH1F h_jetEt3
 
TH1F h_lowPtCal11
 
TH1F h_lowPtCal12
 
TH1F h_lowPtCal13
 
TH1F h_lowPtCal1c11
 
TH1F h_lowPtCal1c12
 
TH1F h_lowPtCal1c13
 
TH1F h_lowPtCal21
 
TH1F h_lowPtCal22
 
TH1F h_lowPtCal23
 
TH1F h_lowPtCal2c11
 
TH1F h_lowPtCal2c12
 
TH1F h_lowPtCal2c13
 
TH1F h_lowPtCal31
 
TH1F h_lowPtCal32
 
TH1F h_lowPtCal33
 
TH1F h_lowPtCal3c11
 
TH1F h_lowPtCal3c12
 
TH1F h_lowPtCal3c13
 
TH1F h_lowPtCal41
 
TH1F h_lowPtCal42
 
TH1F h_lowPtCal43
 
TH1F h_lowPtCal4c11
 
TH1F h_lowPtCal4c12
 
TH1F h_lowPtCal4c13
 
TH1F h_missEt1
 
TH1F h_missEt1s
 
TH1F h_missEt2
 
TH1F h_missEt2s
 
TH1F h_missEt3
 
TH1F h_missEt3s
 
TH1F h_nCalJets1
 
TH1F h_nCalJets2
 
TH1F h_nCalJets3
 
TH1F h_nCalJets4
 
TH1F h_nGenJets1
 
TH1F h_nGenJets2
 
TH1F h_nGenJets3
 
TH1F h_nGenJets4
 
TH1F h_phiCal1
 
TH1F h_phiCal2
 
TH1F h_phiCal3
 
TH1F h_phiCal4
 
TH1F h_phiCalL1
 
TH1F h_phiCalL2
 
TH1F h_phiCalL3
 
TH1F h_phiCalL4
 
TH1F h_phiGen1
 
TH1F h_phiGen2
 
TH1F h_phiGen3
 
TH1F h_phiGen4
 
TH1F h_phiGenL1
 
TH1F h_phiGenL2
 
TH1F h_phiGenL3
 
TH1F h_ptCal1
 
TH1F h_ptCal12
 
TH1F h_ptCal13
 
TH1F h_ptCal2
 
TH1F h_ptCal22
 
TH1F h_ptCal23
 
TH1F h_ptCal3
 
TH1F h_ptCal32
 
TH1F h_ptCal33
 
TH1F h_ptCal4
 
TH1F h_ptCal42
 
TH1F h_ptCal43
 
TH1F h_ptCalL1
 
TH1F h_ptCalL12
 
TH1F h_ptCalL13
 
TH1F h_ptCalL2
 
TH1F h_ptCalL22
 
TH1F h_ptCalL23
 
TH1F h_ptCalL3
 
TH1F h_ptCalL32
 
TH1F h_ptCalL33
 
TH1F h_ptCalL4
 
TH1F h_ptCalL42
 
TH1F h_ptCalL43
 
TH1F h_ptGen1
 
TH1F h_ptGen12
 
TH1F h_ptGen13
 
TH1F h_ptGen2
 
TH1F h_ptGen22
 
TH1F h_ptGen23
 
TH1F h_ptGen3
 
TH1F h_ptGen32
 
TH1F h_ptGen33
 
TH1F h_ptGen4
 
TH1F h_ptGen42
 
TH1F h_ptGen43
 
TH1F h_ptGenL1
 
TH1F h_ptGenL12
 
TH1F h_ptGenL13
 
TH1F h_ptGenL2
 
TH1F h_ptGenL22
 
TH1F h_ptGenL23
 
TH1F h_ptGenL3
 
TH1F h_ptGenL32
 
TH1F h_ptGenL33
 
TH1F h_totMissEt1
 
TH1F h_totMissEt2
 
TH1F h_totMissEt3
 
TH1F hadEneLeadJet1
 
TH1F hadEneLeadJet12
 
TH1F hadEneLeadJet13
 
TH1F hadEneLeadJet2
 
TH1F hadEneLeadJet22
 
TH1F hadEneLeadJet23
 
TH1F hadEneLeadJet3
 
TH1F hadEneLeadJet32
 
TH1F hadEneLeadJet33
 
TH1F hadEneLeadJetEta1_1
 
TH1F hadEneLeadJetEta1_2
 
TH1F hadEneLeadJetEta2_1
 
TH1F hadEneLeadJetEta2_2
 
TH1F hadEneLeadJetEta3_1
 
TH1F hadEneLeadJetEta3_2
 
TH1F hadFracEta11
 
TH1F hadFracEta12
 
TH1F hadFracEta13
 
TH1F hadFracEta21
 
TH1F hadFracEta22
 
TH1F hadFracEta23
 
TH1F hadFracEta31
 
TH1F hadFracEta32
 
TH1F hadFracEta33
 
TH1F hadFracLeadJet1
 
TH1F hadFracLeadJet2
 
TH1F hadFracLeadJet3
 
TProfile hf_nJet1
 
TProfile hf_nJet11
 
TProfile hf_nJet1s
 
TProfile hf_nJet2
 
TProfile hf_nJet21
 
TProfile hf_nJet2s
 
TProfile hf_nJet3
 
TProfile hf_nJet31
 
TProfile hf_nJet3s
 
TProfile hf_nJet4
 
TProfile hf_nJet41
 
TProfile hf_nJet4s
 
TProfile hf_PtResponse1
 
TProfile hf_PtResponse2
 
TProfile hf_PtResponse3
 
TProfile hf_PtResponse4
 
TH1F hf_sumTowerAllEx
 
TH1F hf_sumTowerAllEy
 
TProfile hf_TowerDelR1
 
TProfile hf_TowerDelR12
 
TProfile hf_TowerDelR2
 
TProfile hf_TowerDelR22
 
TProfile hf_TowerDelR3
 
TProfile hf_TowerDelR32
 
TH1F hf_TowerJetEt1
 
std::string JetCorrectionService
 
TFile * m_file
 
TH1F matchedAllPt11
 
TH1F matchedAllPt12
 
TH1F matchedAllPt13
 
TH1F matchedAllPt21
 
TH1F matchedAllPt22
 
TH1F matchedAllPt23
 
TH1F matchedAllPt31
 
TH1F matchedAllPt32
 
TH1F matchedAllPt33
 
TH1F matchedAllPt41
 
TH1F matchedAllPt42
 
TH1F matchedAllPt43
 
TH1F matchedPt11
 
TH1F matchedPt12
 
TH1F matchedPt13
 
TH1F matchedPt21
 
TH1F matchedPt22
 
TH1F matchedPt23
 
TH1F matchedPt31
 
TH1F matchedPt32
 
TH1F matchedPt33
 
TH1F matchedPt41
 
TH1F matchedPt42
 
TH1F matchedPt43
 
TH1F MET1
 
TH1F MET12
 
TH1F MET13
 
TH1F nTowers1
 
TH1F nTowers2
 
TH1F nTowers3
 
TH1F nTowers4
 
TH1F nTowersLeadJet1
 
TH1F nTowersLeadJet2
 
TH1F nTowersLeadJet3
 
TH1F nTowersLeadJetPt1
 
TH1F nTowersLeadJetPt2
 
TH1F nTowersLeadJetPt3
 
TH1F nTowersLeadJetPt4
 
TH1F nTowersSecondJet1
 
TH1F nTowersSecondJet2
 
TH1F nTowersSecondJet3
 
TH1F ParMatch1
 
TH1F ParMatch2
 
TH1F ParMatch3
 
TH1F SumEt1
 
TH1F SumEt12
 
TH1F SumEt13
 
TH1F tbarMass
 
TH1F tbarMassGen
 
TH1F tMass
 
TH1F tMassGen
 
TH1F topMass1
 
TH1F topMass2
 
TH1F topMass3
 
TH1F topMassParton
 
TH1F totEneLeadJetEta1_1
 
TH1F totEneLeadJetEta1_2
 
TH1F totEneLeadJetEta2_1
 
TH1F totEneLeadJetEta2_2
 
TH1F totEneLeadJetEta3_1
 
TH1F totEneLeadJetEta3_2
 
TH1F TowerEtLeadJet1
 
TH1F TowerEtLeadJet12
 
TH1F TowerEtLeadJet13
 
TH1F TowerEtLeadJet2
 
TH1F TowerEtLeadJet22
 
TH1F TowerEtLeadJet23
 
TH1F TowerEtLeadJet3
 
TH1F TowerEtLeadJet32
 
TH1F TowerEtLeadJet33
 
TH1F ZpMass
 
TH1F ZpMass_2000_10
 
TH1F ZpMass_2000_13
 
TH1F ZpMass_2000_40
 
TH1F ZpMass_5000_10
 
TH1F ZpMass_5000_13
 
TH1F ZpMass_5000_40
 
TH1F ZpMass_700_10
 
TH1F ZpMass_700_13
 
TH1F ZpMass_700_40
 
TH1F ZpMassGen
 
TH1F ZpMassGen10
 
TH1F ZpMassGen13
 
TH1F ZpMassGen40
 
TH1F ZpMassGen_2000_10
 
TH1F ZpMassGen_2000_13
 
TH1F ZpMassGen_2000_40
 
TH1F ZpMassGen_5000_10
 
TH1F ZpMassGen_5000_13
 
TH1F ZpMassGen_5000_40
 
TH1F ZpMassGen_700_10
 
TH1F ZpMassGen_700_13
 
TH1F ZpMassGen_700_40
 
TH1F ZpMassMatched1
 
TH1F ZpMassMatched2
 
TH1F ZpMassMatched3
 
TH1F ZpMassRes101
 
TH1F ZpMassRes102
 
TH1F ZpMassRes103
 
TH1F ZpMassRes131
 
TH1F ZpMassRes132
 
TH1F ZpMassRes133
 
TH1F ZpMassRes401
 
TH1F ZpMassRes402
 
TH1F ZpMassRes403
 
TH1F ZpMassResL101
 
TH1F ZpMassResL102
 
TH1F ZpMassResL103
 
TH1F ZpMassResL131
 
TH1F ZpMassResL132
 
TH1F ZpMassResL133
 
TH1F ZpMassResL401
 
TH1F ZpMassResL402
 
TH1F ZpMassResL403
 
TH1F ZpMassResPRL101
 
TH1F ZpMassResPRL102
 
TH1F ZpMassResPRL103
 
TH1F ZpMassResPRL131
 
TH1F ZpMassResPRL132
 
TH1F ZpMassResPRL133
 
TH1F ZpMassResPRL401
 
TH1F ZpMassResPRL402
 
TH1F ZpMassResPRL403
 
TH1F ZpMassResRL101
 
TH1F ZpMassResRL102
 
TH1F ZpMassResRL103
 
TH1F ZpMassResRL131
 
TH1F ZpMassResRL132
 
TH1F ZpMassResRL133
 
TH1F ZpMassResRL401
 
TH1F ZpMassResRL402
 
TH1F ZpMassResRL403
 
TH1F ZpMassResRLoP101
 
TH1F ZpMassResRLoP102
 
TH1F ZpMassResRLoP103
 
TH1F ZpMassResRLoP131
 
TH1F ZpMassResRLoP132
 
TH1F ZpMassResRLoP133
 
TH1F ZpMassResRLoP401
 
TH1F ZpMassResRLoP402
 
TH1F ZpMassResRLoP403
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 16 of file myFastSimVal.h.

Constructor & Destructor Documentation

myFastSimVal::myFastSimVal ( const edm::ParameterSet &  cfg)

Definition at line 34 of file myFastSimVal.cc.

34  :
35  CaloJetAlgorithm1( cfg.getParameter<string>( "CaloJetAlgorithm1" ) ),
36  CaloJetAlgorithm2( cfg.getParameter<string>( "CaloJetAlgorithm2" ) ),
37  CaloJetAlgorithm3( cfg.getParameter<string>( "CaloJetAlgorithm3" ) ),
38  CaloJetAlgorithm4( cfg.getParameter<string>( "CaloJetAlgorithm4" ) ),
39  GenJetAlgorithm1( cfg.getParameter<string>( "GenJetAlgorithm1" ) ),
40  GenJetAlgorithm2( cfg.getParameter<string>( "GenJetAlgorithm2" ) ),
41  GenJetAlgorithm3( cfg.getParameter<string>( "GenJetAlgorithm3" ) ),
42  GenJetAlgorithm4( cfg.getParameter<string>( "GenJetAlgorithm4" ) ),
43  JetCorrectionService( cfg.getParameter<string>( "JetCorrectionService" ) )
44 {
45 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
myFastSimVal::GenJetAlgorithm3
std::string GenJetAlgorithm3
Definition: myFastSimVal.h:26
myFastSimVal::CaloJetAlgorithm3
std::string CaloJetAlgorithm3
Definition: myFastSimVal.h:24
myFastSimVal::GenJetAlgorithm1
std::string GenJetAlgorithm1
Definition: myFastSimVal.h:26
myFastSimVal::CaloJetAlgorithm4
std::string CaloJetAlgorithm4
Definition: myFastSimVal.h:24
myFastSimVal::GenJetAlgorithm4
std::string GenJetAlgorithm4
Definition: myFastSimVal.h:26
myFastSimVal::CaloJetAlgorithm1
std::string CaloJetAlgorithm1
Definition: myFastSimVal.h:24
myFastSimVal::GenJetAlgorithm2
std::string GenJetAlgorithm2
Definition: myFastSimVal.h:26
myFastSimVal::CaloJetAlgorithm2
std::string CaloJetAlgorithm2
Definition: myFastSimVal.h:24
myFastSimVal::JetCorrectionService
std::string JetCorrectionService
Definition: myFastSimVal.h:28

Member Function Documentation

void myFastSimVal::analyze ( const edm::Event &  evt,
const edm::EventSetup &  es 
)
overrideprivate

Definition at line 654 of file myFastSimVal.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), CaloJetAlgorithm1, CaloJetAlgorithm2, CaloJetAlgorithm3, CaloJetAlgorithm4, isolatedTracks_cfi::caloJets, eleHcalExtractorBlocks_cff::caloTowers, funct::cos(), gather_cfg::cout, reco::Candidate::daughter(), hiPixelPairStep_cff::deltaPhi, boostedElectronIsolation_cff::deltaR, dEta1, dEta12, dEta2, dEta3, dEta4, dEtaPar1, dEtaPar2, dEtaPar3, dijetMass1, dijetMass101, dijetMass102, dijetMass103, dijetMass12, dijetMass13, dijetMass131, dijetMass132, dijetMass133, dijetMass2, dijetMass22, dijetMass23, dijetMass3, dijetMass32, dijetMass33, dijetMass4, dijetMass401, dijetMass402, dijetMass403, dijetMass42, dijetMass43, dijetMass_2000_101, dijetMass_2000_131, dijetMass_2000_401, dijetMass_5000_101, dijetMass_5000_131, dijetMass_5000_401, dijetMass_700_101, dijetMass_700_131, dijetMass_700_401, dijetMassCor1, dijetMassCor101, dijetMassCor131, dijetMassCor401, dijetMassCor_2000_1, dijetMassCor_2000_101, dijetMassCor_2000_131, dijetMassCor_2000_401, dijetMassCor_5000_1, dijetMassCor_5000_101, dijetMassCor_5000_131, dijetMassCor_5000_401, dijetMassCor_700_1, dijetMassCor_700_101, dijetMassCor_700_131, dijetMassCor_700_401, dijetMassP1, dijetMassP101, dijetMassP131, dijetMassP2, dijetMassP3, dijetMassP401, dijetMassP_2000_101, dijetMassP_2000_131, dijetMassP_5000_101, dijetMassP_5000_131, dijetMassP_5000_401, dijetMassP_700_101, dijetMassP_700_131, dPhi1, dPhi12, dPhi2, dPhi3, dPhi4, dPhiPar1, dPhiPar2, dPhiPar3, dPt1, dPt100Frac1, dPt100Frac2, dPt100Frac3, dPt100Frac4, dPt12, dPt2, dPt20Frac1, dPt20Frac2, dPt20Frac3, dPt20Frac4, dPt3, dPt4, dPt40Frac1, dPt40Frac2, dPt40Frac3, dPt40Frac4, dPt80Frac1, dPt80Frac2, dPt80Frac3, dPt80Frac4, dPtFrac1, dPtFrac2, dPtFrac3, dPtFrac4, dPtPar1, dPtPar2, dPtPar3, dR1, dR12, dR2, dR3, dR4, dRPar1, dRPar2, dRPar3, dRParton, dRPartonMin, emEneLeadJet1, emEneLeadJet12, emEneLeadJet13, emEneLeadJetEta1_1, emEneLeadJetEta1_2, emEneLeadJetEta2_1, emEneLeadJetEta2_2, emEneLeadJetEta3_1, emEneLeadJetEta3_2, edm::SortedCollection< T, SORT >::end(), stringResolutionProvider_cfi::et, PVValHelper::eta, relval_steps::gen(), GenJetAlgorithm1, GenJetAlgorithm2, GenJetAlgorithm3, GenJetAlgorithm4, ttbarCategorization_cff::genJets, GenHFHadronMatcher_cfi::genParticles, edm::Event::getByLabel(), h_etaCal1, h_etaCal2, h_etaCal3, h_etaCal4, h_etaCalL1, h_etaCalL2, h_etaCalL3, h_etaCalL4, h_etaGen1, h_etaGen2, h_etaGen3, h_etaGenL1, h_etaGenL2, h_etaGenL3, h_jet1Pt1, h_jet1Pt2, h_jet1Pt3, h_jet1Pt4, h_jet2Pt1, h_jet2Pt2, h_jet2Pt3, h_jet2Pt4, h_jet3Pt1, h_jet3Pt2, h_jet3Pt3, h_jet3Pt4, h_jet4Pt1, h_jet4Pt2, h_jet4Pt3, h_jet4Pt4, h_jet5Pt1, h_jet5Pt2, h_jet5Pt3, h_jet5Pt4, h_jet6Pt1, h_jet6Pt2, h_jet6Pt3, h_jet6Pt4, h_jet7Pt1, h_jet7Pt2, h_jet7Pt3, h_jet7Pt4, h_jetEt1, h_jetEt2, h_jetEt3, h_lowPtCal11, h_lowPtCal12, h_lowPtCal13, h_lowPtCal1c11, h_lowPtCal1c12, h_lowPtCal1c13, h_lowPtCal21, h_lowPtCal22, h_lowPtCal23, h_lowPtCal2c11, h_lowPtCal2c12, h_lowPtCal2c13, h_lowPtCal31, h_lowPtCal32, h_lowPtCal33, h_lowPtCal3c11, h_lowPtCal3c12, h_lowPtCal3c13, h_lowPtCal41, h_lowPtCal42, h_lowPtCal43, h_lowPtCal4c11, h_lowPtCal4c12, h_lowPtCal4c13, h_missEt1, h_missEt1s, h_missEt2, h_missEt2s, h_missEt3, h_missEt3s, h_nCalJets1, h_nCalJets2, h_nCalJets3, h_nCalJets4, h_nGenJets1, h_nGenJets2, h_nGenJets3, h_phiCal1, h_phiCal2, h_phiCal3, h_phiCal4, h_phiCalL1, h_phiCalL2, h_phiCalL3, h_phiCalL4, h_phiGen1, h_phiGen2, h_phiGen3, h_phiGenL1, h_phiGenL2, h_phiGenL3, h_ptCal1, h_ptCal12, h_ptCal13, h_ptCal2, h_ptCal22, h_ptCal23, h_ptCal3, h_ptCal32, h_ptCal33, h_ptCal4, h_ptCal42, h_ptCal43, h_ptCalL1, h_ptCalL12, h_ptCalL13, h_ptCalL2, h_ptCalL22, h_ptCalL23, h_ptCalL3, h_ptCalL32, h_ptCalL33, h_ptCalL4, h_ptCalL42, h_ptCalL43, h_ptGen1, h_ptGen12, h_ptGen13, h_ptGen2, h_ptGen22, h_ptGen23, h_ptGen3, h_ptGen32, h_ptGen33, h_ptGenL1, h_ptGenL12, h_ptGenL13, h_ptGenL2, h_ptGenL22, h_ptGenL23, h_ptGenL3, h_ptGenL32, h_ptGenL33, h_totMissEt1, h_totMissEt2, h_totMissEt3, hadEneLeadJet1, hadEneLeadJet12, hadEneLeadJet13, hadEneLeadJetEta1_1, hadEneLeadJetEta1_2, hadEneLeadJetEta2_1, hadEneLeadJetEta2_2, hadEneLeadJetEta3_1, hadEneLeadJetEta3_2, hadFracEta11, hadFracEta21, hadFracEta31, hadFracLeadJet1, hf_nJet1, hf_nJet11, hf_nJet1s, hf_nJet2, hf_nJet21, hf_nJet2s, hf_nJet3, hf_nJet31, hf_nJet3s, hf_nJet4, hf_nJet41, hf_nJet4s, hf_PtResponse1, hf_PtResponse2, hf_PtResponse3, hf_PtResponse4, hf_sumTowerAllEx, hf_sumTowerAllEy, hf_TowerDelR1, hf_TowerDelR12, hf_TowerJetEt1, mps_fire::i, globals_cff::id1, globals_cff::id2, cuy::ii, reco::btau::jetPhi, reco::btau::jetPt, fwrapper::jets, findQualityFiles::jj, ResonanceBuilder::mass, reco::Candidate::mass(), matchedAllPt11, matchedAllPt12, matchedAllPt13, matchedAllPt21, matchedAllPt22, matchedAllPt23, matchedAllPt31, matchedAllPt32, matchedAllPt33, matchedAllPt41, matchedAllPt42, matchedAllPt43, matchedPt11, matchedPt12, matchedPt13, matchedPt21, matchedPt22, matchedPt23, matchedPt31, matchedPt32, matchedPt33, matchedPt41, matchedPt42, matchedPt43, MAXJETS, MET1, MET12, MET13, ALCARECOTkAlWMuNu_cff::minJetPt, reco::Candidate::mother(), custom_jme_cff::nConstituents, nEvent, njet, nPart(), nTowers1, nTowers2, nTowers3, nTowers4, nTowersLeadJet1, nTowersLeadJetPt1, nTowersLeadJetPt2, nTowersLeadJetPt3, nTowersLeadJetPt4, nTowersSecondJet1, reco::Candidate::numberOfDaughters(), AlCaHLTBitMon_ParallelJobs::p, reco::Candidate::p4(), ParMatch1, ParMatch2, ParMatch3, common_cff::pdgId, reco::Candidate::pdgId(), phi, EnergyCorrector::pt, Scenarios_cff::scale, funct::sin(), edm::OwnVector< T, P >::size(), mathSSE::sqrt(), SumEt1, SumEt12, SumEt13, protons_cff::t, tbarMass, tbarMassGen, tMass, tMassGen, topMass1, topMass2, topMass3, topMassParton, totEneLeadJetEta1_1, totEneLeadJetEta1_2, totEneLeadJetEta2_1, totEneLeadJetEta2_2, totEneLeadJetEta3_1, totEneLeadJetEta3_2, TowerEtLeadJet1, TowerEtLeadJet12, TowerEtLeadJet13, reco::Candidate::vx(), reco::Candidate::vy(), reco::Candidate::vz(), ZpMass, ZpMass_2000_10, ZpMass_2000_13, ZpMass_2000_40, ZpMass_5000_10, ZpMass_5000_13, ZpMass_5000_40, ZpMass_700_10, ZpMass_700_13, ZpMass_700_40, ZpMassGen, ZpMassGen10, ZpMassGen13, ZpMassGen40, ZpMassGen_2000_10, ZpMassGen_2000_13, ZpMassGen_2000_40, ZpMassGen_5000_10, ZpMassGen_5000_13, ZpMassGen_5000_40, ZpMassGen_700_10, ZpMassGen_700_13, ZpMassGen_700_40, ZpMassMatched1, ZpMassMatched2, ZpMassMatched3, ZpMassRes101, ZpMassRes131, ZpMassRes401, ZpMassRes402, ZpMassRes403, ZpMassResL101, ZpMassResL102, ZpMassResL103, ZpMassResL131, ZpMassResL132, ZpMassResL133, ZpMassResL401, ZpMassResL402, ZpMassResL403, ZpMassResPRL101, ZpMassResPRL102, ZpMassResPRL103, ZpMassResPRL131, ZpMassResPRL132, ZpMassResPRL133, ZpMassResPRL401, ZpMassResPRL402, ZpMassResPRL403, ZpMassResRL101, ZpMassResRL102, ZpMassResRL103, ZpMassResRL131, ZpMassResRL132, ZpMassResRL133, ZpMassResRL401, ZpMassResRL402, ZpMassResRL403, ZpMassResRLoP101, ZpMassResRLoP102, ZpMassResRLoP103, ZpMassResRLoP131, ZpMassResRLoP132, ZpMassResRLoP133, ZpMassResRLoP401, ZpMassResRLoP402, and ZpMassResRLoP403.

654  {
655 
656  int EtaOk10, EtaOk13, EtaOk40;
657 
658  double ZpM, ZpMG, ZpMM;
659  double LeadMass1, LeadMass2, LeadMass3, LeadMass4;
660  double LeadMassP1, LeadMassP2, LeadMassP3;
661 
662 
663  float minJetPt = 30.;
664  float minJetPt10 = 10.;
665  int jetInd, allJetInd;
666  int usedInd = -1;
667  // double matchedDelR = 0.1;
668  double matchedDelR = 0.3;
669 
670  ZpMG = 0;
671  LeadMass1 = -1;
672  LeadMass2 = -1;
673  LeadMass3 = -1;
674 
675  math::XYZTLorentzVector p4tmp[2], p4cortmp[2];
676  nEvent++;
677 
678  // ********************************
679  // **** Get the CaloJet1 collection
680  // ********************************
681 
682  Handle<CaloJetCollection> caloJets1;
683  evt.getByLabel( CaloJetAlgorithm1, caloJets1 );
684 
685  // Count Jets above Pt cut
686  for (int istep = 0; istep < 100; ++istep) {
687  int njet = 0;
688  float ptStep = (istep * (5000./100.));
689  for ( CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++ cal ) {
690  if ( cal->pt() > ptStep ) njet++;
691  }
692 
693  hf_nJet1.Fill( ptStep, njet );
694  }
695 
696  // Count Jets above Pt cut
697  for (int istep = 0; istep < 100; ++istep) {
698  int njet = 0;
699  float ptStep = (istep * (200./100.));
700  for ( CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++ cal ) {
701  if ( cal->pt() > ptStep ) njet++;
702  }
703 
704  hf_nJet1s.Fill( ptStep, njet );
705  }
706 
707  // Count Jets above Pt cut
708  for (int istep = 0; istep < 100; ++istep) {
709  int njet = 0;
710  float ptStep = (istep * (3000./100.));
711  for ( CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++ cal ) {
712  if ( cal->pt() > ptStep ) njet++;
713  }
714 
715  hf_nJet11.Fill( ptStep, njet );
716  }
717 
718 
719  //Loop over the two leading CaloJets and fill some histograms
720  jetInd = 0;
721  allJetInd = 0;
722  EtaOk10 = 0;
723  EtaOk13 = 0;
724  EtaOk40 = 0;
725 
726  // const JetCorrector* corrector =
727  // JetCorrector::getJetCorrector (JetCorrectionService, es);
728 
729  double highestPt;
730  double nextPt;
731 
732  highestPt = 0.0;
733  nextPt = 0.0;
734 
735 
736  for( CaloJetCollection::const_iterator cal = caloJets1->begin(); cal != caloJets1->end(); ++ cal ) {
737 
738  // double scale = corrector->correction (*cal);
739  double scale = 1.0;
740  double corPt = scale*cal->pt();
741  // double corPt = cal->pt();
742 
743 
744  if (corPt>highestPt) {
745  nextPt = highestPt;
746  p4cortmp[1] = p4cortmp[0];
747  highestPt = corPt;
748  p4cortmp[0] = scale*cal->p4();
749  } else if (corPt>nextPt) {
750  nextPt = corPt;
751  p4cortmp[1] = scale*cal->p4();
752  }
753 
754  /***
755  std::cout << ">>> Corr Jet: corPt = "
756  << corPt << ", scale = " << scale
757  << " pt = " << cal->pt()
758  << " highestPt = " << highestPt
759  << " nextPt = " << nextPt
760  << std::endl;
761  ****/
762 
763  allJetInd++;
764  if (allJetInd == 1) {
765  h_jet1Pt1.Fill( cal->pt() );
766  p4tmp[0] = cal->p4();
767  if ( fabs(cal->eta()) < 1.0) EtaOk10++;
768  if ( fabs(cal->eta()) < 1.3) EtaOk13++;
769  if ( fabs(cal->eta()) < 4.0) EtaOk40++;
770  }
771  if (allJetInd == 2) {
772  h_jet2Pt1.Fill( cal->pt() );
773  p4tmp[1] = cal->p4();
774  if ( fabs(cal->eta()) < 1.0) EtaOk10++;
775  if ( fabs(cal->eta()) < 1.3) EtaOk13++;
776  if ( fabs(cal->eta()) < 4.0) EtaOk40++;
777  }
778  if ( (allJetInd == 1) || (allJetInd == 2) ) {
779 
780  h_ptCalL1.Fill( cal->pt() );
781  h_ptCalL12.Fill( cal->pt() );
782  h_ptCalL13.Fill( cal->pt() );
783 
784  h_etaCalL1.Fill( cal->eta() );
785  h_phiCalL1.Fill( cal->phi() );
786  }
787 
788  if (allJetInd == 3) h_jet3Pt1.Fill( cal->pt() );
789  if (allJetInd == 4) h_jet4Pt1.Fill( cal->pt() );
790  if (allJetInd == 5) h_jet5Pt1.Fill( cal->pt() );
791  if (allJetInd == 6) h_jet6Pt1.Fill( cal->pt() );
792  if (allJetInd == 7) h_jet7Pt1.Fill( cal->pt() );
793 
794  if ( fabs(cal->eta()) < 1.3) {
795  h_lowPtCal11.Fill( cal->pt() );
796  if ( cal->pt() > 10.) {
797  h_lowPtCal1c11.Fill( cal->pt() );
798  }
799  }
800  if ( (fabs(cal->eta())> 1.3) && ( fabs(cal->eta()) < 3.) ) {
801  h_lowPtCal12.Fill( cal->pt() );
802  if ( cal->pt() > 10.) {
803  h_lowPtCal1c12.Fill( cal->pt() );
804  }
805  }
806  if ( fabs(cal->eta()) > 3.0) {
807  h_lowPtCal13.Fill( cal->pt() );
808  if ( cal->pt() > 10.) {
809  h_lowPtCal1c13.Fill( cal->pt() );
810  }
811  }
812 
813 
814 
815  if ( cal->pt() > minJetPt) {
816  // std::cout << "CALO JET1 #" << jetInd << std::endl << cal->print() << std::endl;
817  h_ptCal1.Fill( cal->pt() );
818  h_ptCal12.Fill( cal->pt() );
819  h_ptCal13.Fill( cal->pt() );
820 
821  h_etaCal1.Fill( cal->eta() );
822  h_phiCal1.Fill( cal->phi() );
823  jetInd++;
824  }
825  }
826 
827  // h_nCalJets1.Fill( caloJets1->size() );
828  h_nCalJets1.Fill( jetInd );
829  if (jetInd > 1) {
830  LeadMass1 = (p4tmp[0]+p4tmp[1]).mass();
831  dijetMass1.Fill( LeadMass1 );
832  dijetMass12.Fill( LeadMass1 );
833  dijetMass13.Fill( LeadMass1 );
834  if (EtaOk10 == 2) {
835  dijetMass101.Fill( LeadMass1 );
836  dijetMass102.Fill( LeadMass1 );
837  dijetMass103.Fill( LeadMass1 );
838  dijetMass_700_101.Fill( LeadMass1 );
839  dijetMass_2000_101.Fill( LeadMass1 );
840  dijetMass_5000_101.Fill( LeadMass1 );
841  }
842  if (EtaOk13 == 2) {
843  dijetMass131.Fill( LeadMass1 );
844  dijetMass132.Fill( LeadMass1 );
845  dijetMass133.Fill( LeadMass1 );
846  dijetMass_700_131.Fill( LeadMass1 );
847  dijetMass_2000_131.Fill( LeadMass1 );
848  dijetMass_5000_131.Fill( LeadMass1 );
849  }
850  if (EtaOk40 == 2) {
851  dijetMass401.Fill( LeadMass1 );
852  dijetMass402.Fill( LeadMass1 );
853  dijetMass403.Fill( LeadMass1 );
854  dijetMass_700_401.Fill( LeadMass1 );
855  dijetMass_2000_401.Fill( LeadMass1 );
856  dijetMass_5000_401.Fill( LeadMass1 );
857  }
858 
859  LeadMass1 = (p4cortmp[0]+p4cortmp[1]).mass();
860 
861  /****
862  if (LeadMass1 < 30.) {
863  std::cout << " XXX Low Mass "
864  << (p4tmp[0]+p4tmp[1]).mass()
865  << " / "
866  << (p4cortmp[0]+p4cortmp[1]).mass()
867  << std::endl;
868 
869  std::cout << " p4 1 = " << p4tmp[0]
870  << " p4 2 = " << p4tmp[1]
871  << " p4 cor 1 = " << p4cortmp[0]
872  << " p4 cor 2 = " << p4cortmp[0]
873  << endl;
874 
875  }
876  ****/
877 
878  /****
879  dijetMassCor1.Fill( LeadMass1 );
880  dijetMassCor_700_1.Fill( LeadMass1 );
881  dijetMassCor_2000_1.Fill( LeadMass1 );
882  dijetMassCor_5000_1.Fill( LeadMass1 );
883 
884  if (EtaOk10 == 2) {
885  dijetMassCor101.Fill( LeadMass1 );
886  dijetMassCor_700_101.Fill( LeadMass1 );
887  dijetMassCor_2000_101.Fill( LeadMass1 );
888  dijetMassCor_5000_101.Fill( LeadMass1 );
889  }
890  if (EtaOk13 == 2) {
891  dijetMassCor131.Fill( LeadMass1 );
892  dijetMassCor_700_131.Fill( LeadMass1 );
893  dijetMassCor_2000_131.Fill( LeadMass1 );
894  dijetMassCor_5000_131.Fill( LeadMass1 );
895  }
896  if (EtaOk40 == 2) {
897  dijetMassCor401.Fill( LeadMass1 );
898  dijetMassCor_700_401.Fill( LeadMass1 );
899  dijetMassCor_2000_401.Fill( LeadMass1 );
900  dijetMassCor_5000_401.Fill( LeadMass1 );
901  }
902  ****/
903 
904  }
905 
906 
907  // ********************************
908  // **** Get the CaloJet2 collection
909  // ********************************
910  Handle<CaloJetCollection> caloJets2;
911  evt.getByLabel( CaloJetAlgorithm2, caloJets2 );
912 
913  // Count Jets above Pt cut
914  for (int istep = 0; istep < 100; ++istep) {
915  int njet = 0;
916  float ptStep = (istep * (5000./100.));
917 
918  for ( CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++ cal )
919  if ( cal->pt() > ptStep ) njet++;
920 
921  hf_nJet2.Fill( ptStep, njet );
922  }
923 
924  for (int istep = 0; istep < 100; ++istep) {
925  int njet = 0;
926  float ptStep = (istep * (200./100.));
927 
928  for ( CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++ cal )
929  if ( cal->pt() > ptStep ) njet++;
930 
931  hf_nJet2s.Fill( ptStep, njet );
932  }
933 
934 
935  for (int istep = 0; istep < 100; ++istep) {
936  int njet = 0;
937  float ptStep = (istep * (3000./100.));
938 
939  for ( CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++ cal )
940  if ( cal->pt() > ptStep ) njet++;
941 
942  hf_nJet21.Fill( ptStep, njet );
943  }
944 
945 
946 
947 
948 
949  //Loop over the two leading CaloJets and fill some histograms
950  jetInd = 0;
951  allJetInd = 0;
952  for( CaloJetCollection::const_iterator cal = caloJets2->begin(); cal != caloJets2->end(); ++cal ) {
953 
954  allJetInd++;
955  if (allJetInd == 1) {
956  h_jet1Pt2.Fill( cal->pt() );
957  p4tmp[0] = cal->p4();
958  }
959  if (allJetInd == 2) {
960  h_jet2Pt2.Fill( cal->pt() );
961  p4tmp[1] = cal->p4();
962  }
963  if ( (allJetInd == 1) || (allJetInd == 2) ) {
964  h_ptCalL2.Fill( cal->pt() );
965  h_ptCalL22.Fill( cal->pt() );
966  h_ptCalL23.Fill( cal->pt() );
967 
968  h_etaCalL2.Fill( cal->eta() );
969  h_phiCalL2.Fill( cal->phi() );
970  }
971  if (allJetInd == 3) h_jet3Pt2.Fill( cal->pt() );
972  if (allJetInd == 4) h_jet4Pt2.Fill( cal->pt() );
973  if (allJetInd == 5) h_jet5Pt2.Fill( cal->pt() );
974  if (allJetInd == 6) h_jet6Pt2.Fill( cal->pt() );
975  if (allJetInd == 7) h_jet7Pt2.Fill( cal->pt() );
976 
977  if ( fabs(cal->eta()) < 1.3) {
978  h_lowPtCal21.Fill( cal->pt() );
979  if ( cal->pt() > 10.) {
980  h_lowPtCal2c11.Fill( cal->pt() );
981  }
982  }
983  if ( (fabs(cal->eta())> 1.3) && ( fabs(cal->eta()) < 3.) ) {
984  h_lowPtCal22.Fill( cal->pt() );
985  if ( cal->pt() > 10.) {
986  h_lowPtCal2c12.Fill( cal->pt() );
987  }
988  }
989  if ( fabs(cal->eta()) > 3.0) {
990  h_lowPtCal23.Fill( cal->pt() );
991  if ( cal->pt() > 10.) {
992  h_lowPtCal2c13.Fill( cal->pt() );
993  }
994  }
995 
996 
997  if ( cal->pt() > minJetPt) {
998  h_ptCal2.Fill( cal->pt() );
999  h_ptCal22.Fill( cal->pt() );
1000  h_ptCal23.Fill( cal->pt() );
1001 
1002  h_etaCal2.Fill( cal->eta() );
1003  h_phiCal2.Fill( cal->phi() );
1004  jetInd++;
1005  }
1006  }
1007  // h_nCalJets2.Fill( caloJets2->size() );
1008  h_nCalJets2.Fill( jetInd );
1009  if (jetInd > 1) {
1010  LeadMass2 = (p4tmp[0]+p4tmp[1]).mass();
1011  dijetMass2.Fill( LeadMass2 );
1012  dijetMass22.Fill( LeadMass2 );
1013  dijetMass23.Fill( LeadMass2 );
1014 
1015 
1016  dijetMassCor1.Fill( LeadMass2 );
1017  dijetMassCor_700_1.Fill( LeadMass2 );
1018  dijetMassCor_2000_1.Fill( LeadMass2 );
1019  dijetMassCor_5000_1.Fill( LeadMass2 );
1020 
1021  if (EtaOk10 == 2) {
1022  dijetMassCor101.Fill( LeadMass2 );
1023  dijetMassCor_700_101.Fill( LeadMass2 );
1024  dijetMassCor_2000_101.Fill( LeadMass2 );
1025  dijetMassCor_5000_101.Fill( LeadMass2 );
1026  }
1027  if (EtaOk13 == 2) {
1028  dijetMassCor131.Fill( LeadMass2 );
1029  dijetMassCor_700_131.Fill( LeadMass2 );
1030  dijetMassCor_2000_131.Fill( LeadMass2 );
1031  dijetMassCor_5000_131.Fill( LeadMass2 );
1032  }
1033  if (EtaOk40 == 2) {
1034  dijetMassCor401.Fill( LeadMass2 );
1035  dijetMassCor_700_401.Fill( LeadMass2 );
1036  dijetMassCor_2000_401.Fill( LeadMass2 );
1037  dijetMassCor_5000_401.Fill( LeadMass2 );
1038  }
1039 
1040 
1041  }
1042 
1043 
1044 
1045  // ********************************
1046  // **** Get the CaloJet3 collection
1047  // ********************************
1048  Handle<CaloJetCollection> caloJets3;
1049  evt.getByLabel( CaloJetAlgorithm3, caloJets3 );
1050 
1051  //Loop over the two leading CaloJets and fill some histograms
1052  jetInd = 0;
1053  allJetInd = 0;
1054 
1055  // Count Jets above Pt cut
1056  for (int istep = 0; istep < 100; ++istep) {
1057  int njet = 0;
1058  float ptStep = (istep * (5000./100.));
1059 
1060  for ( CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++ cal )
1061  if ( cal->pt() > ptStep ) njet++;
1062 
1063 
1064  hf_nJet3.Fill( ptStep, njet );
1065  }
1066 
1067  for (int istep = 0; istep < 100; ++istep) {
1068  int njet = 0;
1069  float ptStep = (istep * (200./100.));
1070 
1071  for ( CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++ cal )
1072  if ( cal->pt() > ptStep ) njet++;
1073 
1074  hf_nJet3s.Fill( ptStep, njet );
1075  }
1076 
1077  for (int istep = 0; istep < 100; ++istep) {
1078  int njet = 0;
1079  float ptStep = (istep * (3000./100.));
1080 
1081  for ( CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++ cal )
1082  if ( cal->pt() > ptStep ) njet++;
1083 
1084  hf_nJet31.Fill( ptStep, njet );
1085  }
1086 
1087 
1088  for( CaloJetCollection::const_iterator cal = caloJets3->begin(); cal != caloJets3->end(); ++ cal ) {
1089 
1090  allJetInd++;
1091  if (allJetInd == 1) {
1092  h_jet1Pt3.Fill( cal->pt() );
1093  p4tmp[0] = cal->p4();
1094  }
1095  if (allJetInd == 2) {
1096  h_jet2Pt3.Fill( cal->pt() );
1097  p4tmp[1] = cal->p4();
1098  }
1099  if ( (allJetInd == 1) || (allJetInd == 2) ) {
1100  h_ptCalL3.Fill( cal->pt() );
1101  h_ptCalL32.Fill( cal->pt() );
1102  h_ptCalL33.Fill( cal->pt() );
1103 
1104  h_etaCalL3.Fill( cal->eta() );
1105  h_phiCalL3.Fill( cal->phi() );
1106  }
1107  if (allJetInd == 3) h_jet3Pt3.Fill( cal->pt() );
1108  if (allJetInd == 4) h_jet4Pt3.Fill( cal->pt() );
1109  if (allJetInd == 5) h_jet5Pt3.Fill( cal->pt() );
1110  if (allJetInd == 6) h_jet6Pt3.Fill( cal->pt() );
1111  if (allJetInd == 7) h_jet7Pt3.Fill( cal->pt() );
1112 
1113 
1114  if ( fabs(cal->eta()) < 1.3) {
1115  h_lowPtCal31.Fill( cal->pt() );
1116  if ( cal->pt() > 10.) {
1117  h_lowPtCal3c11.Fill( cal->pt() );
1118  }
1119  }
1120  if ( (fabs(cal->eta())> 1.3) && ( fabs(cal->eta()) < 3.) ) {
1121  h_lowPtCal32.Fill( cal->pt() );
1122  if ( cal->pt() > 10.) {
1123  h_lowPtCal3c12.Fill( cal->pt() );
1124  }
1125  }
1126  if ( fabs(cal->eta()) > 3.0) {
1127  h_lowPtCal33.Fill( cal->pt() );
1128  if ( cal->pt() > 10.) {
1129  h_lowPtCal3c13.Fill( cal->pt() );
1130  }
1131  }
1132 
1133 
1134 
1135  if ( cal->pt() > minJetPt) {
1136  // std::cout << "CALO JET3 #" << jetInd << std::endl << cal->print() << std::endl;
1137  h_ptCal3.Fill( cal->pt() );
1138  h_ptCal32.Fill( cal->pt() );
1139  h_ptCal33.Fill( cal->pt() );
1140 
1141  h_etaCal3.Fill( cal->eta() );
1142  h_phiCal3.Fill( cal->phi() );
1143  jetInd++;
1144  }
1145  }
1146  // h_nCalJets3.Fill( caloJets3->size() );
1147  h_nCalJets3.Fill( jetInd );
1148  if (jetInd > 1) {
1149  LeadMass3 = (p4tmp[0]+p4tmp[1]).mass();
1150  dijetMass3.Fill( LeadMass3 );
1151  dijetMass32.Fill( LeadMass3 );
1152  dijetMass33.Fill( LeadMass3 );
1153  }
1154 
1155 
1156 
1157  // ********************************
1158  // **** Get the CaloJet4 collection
1159  // ********************************
1160 
1161  Handle<CaloJetCollection> caloJets4;
1162  evt.getByLabel( CaloJetAlgorithm4, caloJets4 );
1163 
1164  //Loop over the two leading CaloJets and fill some histograms
1165  jetInd = 0;
1166  allJetInd = 0;
1167 
1168  // Count Jets above Pt cut
1169  for (int istep = 0; istep < 100; ++istep) {
1170  int njet = 0;
1171  float ptStep = (istep * (5000./100.));
1172 
1173  for ( CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++ cal )
1174  if ( cal->pt() > ptStep ) njet++;
1175 
1176 
1177  hf_nJet4.Fill( ptStep, njet );
1178  }
1179 
1180  for (int istep = 0; istep < 100; ++istep) {
1181  int njet = 0;
1182  float ptStep = (istep * (200./100.));
1183 
1184  for ( CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++ cal )
1185  if ( cal->pt() > ptStep ) njet++;
1186 
1187  hf_nJet4s.Fill( ptStep, njet );
1188  }
1189 
1190  for (int istep = 0; istep < 100; ++istep) {
1191  int njet = 0;
1192  float ptStep = (istep * (3000./100.));
1193 
1194  for ( CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++ cal )
1195  if ( cal->pt() > ptStep ) njet++;
1196 
1197  hf_nJet41.Fill( ptStep, njet );
1198  }
1199 
1200 
1201  for( CaloJetCollection::const_iterator cal = caloJets4->begin(); cal != caloJets4->end(); ++ cal ) {
1202 
1203  allJetInd++;
1204  if (allJetInd == 1) {
1205  h_jet1Pt4.Fill( cal->pt() );
1206  p4tmp[0] = cal->p4();
1207  }
1208  if (allJetInd == 2) {
1209  h_jet2Pt4.Fill( cal->pt() );
1210  p4tmp[1] = cal->p4();
1211  }
1212  if ( (allJetInd == 1) || (allJetInd == 2) ) {
1213  h_ptCalL4.Fill( cal->pt() );
1214  h_ptCalL42.Fill( cal->pt() );
1215  h_ptCalL43.Fill( cal->pt() );
1216 
1217  h_etaCalL4.Fill( cal->eta() );
1218  h_phiCalL4.Fill( cal->phi() );
1219  }
1220  if (allJetInd == 3) h_jet3Pt4.Fill( cal->pt() );
1221  if (allJetInd == 4) h_jet4Pt4.Fill( cal->pt() );
1222  if (allJetInd == 5) h_jet5Pt4.Fill( cal->pt() );
1223  if (allJetInd == 6) h_jet6Pt4.Fill( cal->pt() );
1224  if (allJetInd == 7) h_jet7Pt4.Fill( cal->pt() );
1225 
1226 
1227  if ( fabs(cal->eta()) < 1.3) {
1228  h_lowPtCal41.Fill( cal->pt() );
1229  if ( cal->pt() > 10.) {
1230  h_lowPtCal4c11.Fill( cal->pt() );
1231  }
1232  }
1233  if ( (fabs(cal->eta())> 1.3) && ( fabs(cal->eta()) < 3.) ) {
1234  h_lowPtCal42.Fill( cal->pt() );
1235  if ( cal->pt() > 10.) {
1236  h_lowPtCal4c12.Fill( cal->pt() );
1237  }
1238  }
1239  if ( fabs(cal->eta()) > 3.0) {
1240  h_lowPtCal43.Fill( cal->pt() );
1241  if ( cal->pt() > 10.) {
1242  h_lowPtCal4c13.Fill( cal->pt() );
1243  }
1244  }
1245 
1246 
1247  if ( cal->pt() > minJetPt) {
1248  // std::cout << "CALO JET4 #" << jetInd << std::endl << cal->print() << std::endl;
1249  h_ptCal4.Fill( cal->pt() );
1250  h_ptCal42.Fill( cal->pt() );
1251  h_ptCal43.Fill( cal->pt() );
1252 
1253  h_etaCal4.Fill( cal->eta() );
1254  h_phiCal4.Fill( cal->phi() );
1255  jetInd++;
1256  }
1257  }
1258  // h_nCalJets4.Fill( caloJets4->size() );
1259  h_nCalJets4.Fill( jetInd );
1260  if (jetInd > 1) {
1261  LeadMass4 = (p4tmp[0]+p4tmp[1]).mass();
1262  dijetMass4.Fill( LeadMass4 );
1263  dijetMass42.Fill( LeadMass4 );
1264  dijetMass43.Fill( LeadMass4 );
1265  }
1266 
1267 
1268 
1269  // *********************************************
1270  // *********************************************
1271 
1272 
1273  //**** Get the GenJet1 collection
1274  Handle<GenJetCollection> genJets1;
1275  evt.getByLabel( GenJetAlgorithm1, genJets1 );
1276 
1277 
1278  //Loop over the two leading GenJets and fill some histograms
1279  jetInd = 0;
1280  allJetInd = 0;
1281  for( GenJetCollection::const_iterator gen = genJets1->begin(); gen != genJets1->end(); ++ gen ) {
1282  allJetInd++;
1283  if (allJetInd == 1) {
1284  p4tmp[0] = gen->p4();
1285  }
1286  if (allJetInd == 2) {
1287  p4tmp[1] = gen->p4();
1288  }
1289 
1290  if ( (allJetInd == 1) || (allJetInd == 2) ) {
1291  h_ptGenL1.Fill( gen->pt() );
1292  h_ptGenL12.Fill( gen->pt() );
1293  h_ptGenL13.Fill( gen->pt() );
1294 
1295  h_etaGenL1.Fill( gen->eta() );
1296  h_phiGenL1.Fill( gen->phi() );
1297  }
1298 
1299  if ( gen->pt() > minJetPt) {
1300  // std::cout << "GEN JET1 #" << jetInd << std::endl << gen->print() << std::endl;
1301  h_ptGen1.Fill( gen->pt() );
1302  h_ptGen12.Fill( gen->pt() );
1303  h_ptGen13.Fill( gen->pt() );
1304 
1305  h_etaGen1.Fill( gen->eta() );
1306  h_phiGen1.Fill( gen->phi() );
1307  jetInd++;
1308  }
1309  }
1310 
1311  LeadMassP1 = (p4tmp[0]+p4tmp[1]).mass();
1312  dijetMassP1.Fill( LeadMassP1 );
1313  if (EtaOk10 == 2) {
1314  dijetMassP101.Fill( LeadMassP1 );
1315  dijetMassP_700_101.Fill( LeadMassP1 );
1316  dijetMassP_2000_101.Fill( LeadMassP1 );
1317  dijetMassP_5000_101.Fill( LeadMassP1 );
1318  }
1319  if (EtaOk13 == 2) {
1320  dijetMassP131.Fill( LeadMassP1 );
1321  dijetMassP_700_131.Fill( LeadMassP1 );
1322  dijetMassP_2000_131.Fill( LeadMassP1 );
1323  dijetMassP_5000_131.Fill( LeadMassP1 );
1324 
1325  }
1326  if (EtaOk40 == 2) {
1327  dijetMassP401.Fill( LeadMassP1 );
1328  dijetMassP_5000_401.Fill( LeadMassP1 );
1329  dijetMassP_5000_401.Fill( LeadMassP1 );
1330  dijetMassP_5000_401.Fill( LeadMassP1 );
1331  }
1332 
1333  // h_nGenJets1.Fill( genJets1->size() );
1334  h_nGenJets1.Fill( jetInd );
1335 
1336  //**** Get the GenJet2 collection
1337  Handle<GenJetCollection> genJets2;
1338  evt.getByLabel( GenJetAlgorithm2, genJets2 );
1339 
1340  //Loop over the two leading GenJets and fill some histograms
1341  jetInd = 0;
1342  allJetInd = 0;
1343  for( GenJetCollection::const_iterator gen = genJets2->begin(); gen != genJets2->end(); ++ gen ) {
1344  allJetInd++;
1345  if (allJetInd == 1) {
1346  p4tmp[0] = gen->p4();
1347  }
1348  if (allJetInd == 2) {
1349  p4tmp[1] = gen->p4();
1350  }
1351  if ( (allJetInd == 1) || (allJetInd == 2) ) {
1352  h_ptGenL2.Fill( gen->pt() );
1353  h_ptGenL22.Fill( gen->pt() );
1354  h_ptGenL23.Fill( gen->pt() );
1355 
1356  h_etaGenL2.Fill( gen->eta() );
1357  h_phiGenL2.Fill( gen->phi() );
1358  }
1359 
1360  if ( gen->pt() > minJetPt) {
1361  // std::cout << "GEN JET2 #" << jetInd << std::endl << gen->print() << std::endl;
1362  h_ptGen2.Fill( gen->pt() );
1363  h_ptGen22.Fill( gen->pt() );
1364  h_ptGen23.Fill( gen->pt() );
1365 
1366  h_etaGen2.Fill( gen->eta() );
1367  h_phiGen2.Fill( gen->phi() );
1368  jetInd++;
1369  }
1370  }
1371 
1372  LeadMassP2 = (p4tmp[0]+p4tmp[1]).mass();
1373  dijetMassP2.Fill( LeadMassP2 );
1374 
1375 
1376  // h_nGenJets2.Fill( genJets2->size() );
1377  h_nGenJets2.Fill( jetInd );
1378 
1379  //**** Get the GenJet3 collection
1380  Handle<GenJetCollection> genJets3;
1381  evt.getByLabel( GenJetAlgorithm3, genJets3 );
1382 
1383  //Loop over the two leading GenJets and fill some histograms
1384  jetInd = 0;
1385  allJetInd = 0;
1386  for( GenJetCollection::const_iterator gen = genJets3->begin(); gen != genJets3->end(); ++ gen ) {
1387  allJetInd++;
1388  if (allJetInd == 1) {
1389  p4tmp[0] = gen->p4();
1390  }
1391  if (allJetInd == 2) {
1392  p4tmp[1] = gen->p4();
1393  }
1394  if ( (allJetInd == 1) || (allJetInd == 2) ) {
1395  h_ptGenL3.Fill( gen->pt() );
1396  h_ptGenL32.Fill( gen->pt() );
1397  h_ptGenL33.Fill( gen->pt() );
1398 
1399  h_etaGenL3.Fill( gen->eta() );
1400  h_phiGenL3.Fill( gen->phi() );
1401  }
1402 
1403  if ( gen->pt() > minJetPt) {
1404  // std::cout << "GEN JET3 #" << jetInd << std::endl << gen->print() << std::endl;
1405  h_ptGen3.Fill( gen->pt() );
1406  h_ptGen32.Fill( gen->pt() );
1407  h_ptGen33.Fill( gen->pt() );
1408 
1409  h_etaGen3.Fill( gen->eta() );
1410  h_phiGen3.Fill( gen->phi() );
1411  jetInd++;
1412  }
1413  }
1414 
1415  LeadMassP3 = (p4tmp[0]+p4tmp[1]).mass();
1416  dijetMassP3.Fill( LeadMassP3 );
1417 
1418 
1419  // h_nGenJets3.Fill( genJets3->size() );
1420  h_nGenJets3.Fill( jetInd );
1421 
1422 
1423  // *********************
1424  // MidPoint Jet Matching
1425 
1426  Handle<GenJetCollection> genJets;
1427  Handle<CaloJetCollection> caloJets;
1428 
1429  // evt.getByLabel( "midPointCone5GenJets", genJets );
1430  // evt.getByLabel( "midPointCone5CaloJets", caloJets );
1431  evt.getByLabel( GenJetAlgorithm1, genJets );
1432  evt.getByLabel( CaloJetAlgorithm1, caloJets );
1433 
1434 
1435  int maxJets = MAXJETS;
1436 
1437  jetInd = 0;
1438  double dRmin[MAXJETS];
1439  math::XYZTLorentzVector p4gen[MAXJETS], p4cal[MAXJETS],
1440  p4par[MAXJETS], p4Zp[MAXJETS], p4part[MAXJETS];
1441 
1442  int used[MAXJETS];
1443  int nj;
1444 
1445  for( int i=0; i<maxJets; ++i ) used[i] = 0;
1446 
1447  // cout << ">>>>>>>>> " << endl;
1448 
1449 
1450 
1451  for( GenJetCollection::const_iterator gen = genJets->begin();
1452  gen != genJets->end() && jetInd<maxJets; ++ gen ) {
1453 
1454  p4gen[jetInd] = gen->p4(); //Gen 4-vector
1455  dRmin[jetInd] = 1000.0;
1456 
1457  nj = 0;
1458  usedInd = -1;
1459 
1460  for( CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++ cal ) {
1461 
1462  double delR = deltaR( cal->eta(), cal->phi(), gen->eta(), gen->phi() );
1463 
1464  if ( (delR<dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
1465  dRmin[jetInd] = delR; // delta R of match
1466  p4cal[jetInd] = cal->p4(); // Matched Cal 4-vector
1467  usedInd = nj;
1468  }
1469 
1470  nj++;
1471  }
1472 
1473  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1474  matchedAllPt11.Fill(p4gen[jetInd].pt());
1475  }
1476  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1477  matchedAllPt12.Fill(p4gen[jetInd].pt());
1478  }
1479  if (fabs(p4gen[jetInd].eta()) > 3.) {
1480  matchedAllPt13.Fill(p4gen[jetInd].pt());
1481  }
1482 
1483  if (usedInd != -1) {
1484 
1485  used[usedInd] = 1;
1486 
1487  if (p4cal[jetInd].pt() > minJetPt10)
1488  hf_PtResponse1.Fill(p4cal[jetInd].eta(), p4cal[jetInd].pt()/p4gen[jetInd].pt());
1489 
1490  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1491  matchedPt11.Fill(p4gen[jetInd].pt());
1492  }
1493  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1494  matchedPt12.Fill(p4gen[jetInd].pt());
1495  }
1496  if (fabs(p4gen[jetInd].eta()) > 3.) {
1497  matchedPt13.Fill(p4gen[jetInd].pt());
1498  }
1499  /***
1500  cout << " >>>SISCone "
1501  << jetInd << " "
1502  << p4gen[jetInd].pt() << " "
1503  << p4gen[jetInd].eta() << " "
1504  << p4gen[jetInd].phi() << " "
1505  << p4cal[jetInd].pt() << " "
1506  << p4cal[jetInd].eta() << " "
1507  << p4cal[jetInd].phi()
1508  << endl;
1509  ***/
1510 
1511  double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1512  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1513  if ( (p4gen[jetInd].pt() > 20) && (p4gen[jetInd].pt() < 40) ) {
1514  dPt20Frac1.Fill(dpt/p4gen[jetInd].pt());
1515  }
1516  if ( (p4gen[jetInd].pt() > 40) && (p4gen[jetInd].pt() < 60) ) {
1517  dPt40Frac1.Fill(dpt/p4gen[jetInd].pt());
1518  }
1519  if ( (p4gen[jetInd].pt() > 80) && (p4gen[jetInd].pt() < 100) ) {
1520  dPt80Frac1.Fill(dpt/p4gen[jetInd].pt());
1521  }
1522  if ( (p4gen[jetInd].pt() > 100) && (p4gen[jetInd].pt() < 120) ) {
1523  dPt100Frac1.Fill(dpt/p4gen[jetInd].pt());
1524  }
1525  }
1526  if ( (p4gen[jetInd].pt() > minJetPt10) && (fabs(p4gen[jetInd].eta()) < 1.3) ) {
1527 
1528  dR1.Fill(dRmin[jetInd]);
1529  double dphi = deltaPhi(p4cal[jetInd].phi(), p4gen[jetInd].phi());
1530  dPhi1.Fill(dphi);
1531  double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1532  dEta1.Fill(deta);
1533  dPt1.Fill(dpt);
1534  dPtFrac1.Fill(dpt/p4gen[jetInd].pt());
1535 
1536  if ( ( (dpt/p4gen[jetInd].pt()) < -0.5 ) && ( fabs(dpt) > 90. ) ) {
1537 
1538  cout << " deltaR min = " << dRmin[jetInd]
1539  << " Ind = " << jetInd << " / " << usedInd << " / " << used[nj]
1540  << " Del pt / frac = " << dpt
1541  << " / " << dpt/p4gen[jetInd].pt()
1542  << " cal/gen pt = " << p4cal[jetInd].pt()
1543  << " / " << p4gen[jetInd].pt()
1544  << " cal/gen eta = " << p4cal[jetInd].eta()
1545  << " / " << p4gen[jetInd].eta()
1546  << " cal/gen phi = " << p4cal[jetInd].phi()
1547  << " / " << p4gen[jetInd].phi()
1548  << endl;
1549  }
1550 
1551  }
1552 
1553  jetInd++;
1554 
1555  }
1556 
1557  }
1558 
1559 
1560 
1561  // *********************
1562  // Seedless Jet Matching
1563 
1564  // Handle<GenJetCollection> genJets;
1565  // Handle<CaloJetCollection> caloJets;
1566 
1567  // evt.getByLabel( "sisCone5GenJets", genJets );
1568  // evt.getByLabel( "sisCone5CaloJets", caloJets );
1569  evt.getByLabel( GenJetAlgorithm2, genJets );
1570  evt.getByLabel( CaloJetAlgorithm2, caloJets );
1571 
1572  // int maxJets = 20;
1573  jetInd = 0;
1574  // double dRmin[20];
1575  // math::XYZTLorentzVector p4jet[20], p4gen[20], p4cal[20], p4cor[20];
1576 
1577  for( int i=0; i<maxJets; ++i ) used[i] = 0;
1578  for( GenJetCollection::const_iterator gen = genJets->begin();
1579  gen != genJets->end() && jetInd<maxJets; ++ gen ) {
1580  p4gen[jetInd] = gen->p4(); //Gen 4-vector
1581  dRmin[jetInd] = 1000.0;
1582 
1583  nj = 0;
1584  usedInd = -1;
1585 
1586  for( CaloJetCollection::const_iterator cal = caloJets->begin();
1587  cal != caloJets->end(); ++ cal ) {
1588  double delR = deltaR( cal->eta(), cal->phi(), gen->eta(), gen->phi() );
1589 
1590  if ( (delR<dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
1591  dRmin[jetInd] = delR; // delta R of match
1592  p4cal[jetInd] = cal->p4(); // Matched Cal 4-vector
1593  usedInd = nj;
1594  }
1595  nj++;
1596  }
1597 
1598  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1599  matchedAllPt21.Fill(p4gen[jetInd].pt());
1600  }
1601  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1602  matchedAllPt22.Fill(p4gen[jetInd].pt());
1603  }
1604  if (fabs(p4gen[jetInd].eta()) > 3.) {
1605  matchedAllPt23.Fill(p4gen[jetInd].pt());
1606  }
1607 
1608 
1609  if (usedInd != -1) {
1610 
1611  used[usedInd] = 1;
1612 
1613  if (p4cal[jetInd].pt() > minJetPt10)
1614  hf_PtResponse2.Fill(p4cal[jetInd].eta(), p4cal[jetInd].pt()/p4gen[jetInd].pt());
1615 
1616  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1617  matchedPt21.Fill(p4gen[jetInd].pt());
1618  }
1619  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1620  matchedPt22.Fill(p4gen[jetInd].pt());
1621  }
1622  if (fabs(p4gen[jetInd].eta()) > 3.) {
1623  matchedPt23.Fill(p4gen[jetInd].pt());
1624  }
1625  /***
1626  cout << " >>>IterCone "
1627  << jetInd << " "
1628  << p4gen[jetInd].pt() << " "
1629  << p4gen[jetInd].eta() << " "
1630  << p4gen[jetInd].phi() << " "
1631  << p4cal[jetInd].pt() << " "
1632  << p4cal[jetInd].eta() << " "
1633  << p4cal[jetInd].phi()
1634  << endl;
1635  ***/
1636 
1637 
1638  double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1639  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1640  if ( (p4gen[jetInd].pt() > 20) && (p4gen[jetInd].pt() < 40) ) {
1641  dPt20Frac2.Fill(dpt/p4gen[jetInd].pt());
1642  }
1643  if ( (p4gen[jetInd].pt() > 40) && (p4gen[jetInd].pt() < 60) ) {
1644  dPt40Frac2.Fill(dpt/p4gen[jetInd].pt());
1645  }
1646  if ( (p4gen[jetInd].pt() > 80) && (p4gen[jetInd].pt() < 100) ) {
1647  dPt80Frac2.Fill(dpt/p4gen[jetInd].pt());
1648  }
1649  if ( (p4gen[jetInd].pt() > 100) && (p4gen[jetInd].pt() < 120) ) {
1650  dPt100Frac2.Fill(dpt/p4gen[jetInd].pt());
1651  }
1652  }
1653  if ( (p4gen[jetInd].pt() > minJetPt10) && (fabs(p4gen[jetInd].eta()) < 1.3) ) {
1654 
1655  dR2.Fill(dRmin[jetInd]);
1656  double dphi = deltaPhi(p4cal[jetInd].phi(), p4gen[jetInd].phi());
1657  dPhi2.Fill(dphi);
1658  double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1659  dEta2.Fill(deta);
1660  dPt2.Fill(dpt);
1661  dPtFrac2.Fill(dpt/p4gen[jetInd].pt());
1662 
1663  }
1664 
1665  jetInd++;
1666  }
1667 
1668  }
1669 
1670  // *********************
1671  // Kt Jet Matching
1672 
1673  // Handle<GenJetCollection> genJets;
1674  // Handle<CaloJetCollection> caloJets;
1675 
1676  // evt.getByLabel( "sisCone5GenJets", genJets );
1677  // evt.getByLabel( "sisCone5CaloJets", caloJets );
1678  evt.getByLabel( GenJetAlgorithm3, genJets );
1679  evt.getByLabel( CaloJetAlgorithm3, caloJets );
1680 
1681  // int maxJets = 20;
1682  jetInd = 0;
1683  // double dRmin[20];
1684  // math::XYZTLorentzVector p4jet[20], p4gen[20], p4cal[20], p4cor[20];
1685 
1686  for( int i=0; i<maxJets; ++i ) used[i] = 0;
1687  for( GenJetCollection::const_iterator gen = genJets->begin();
1688  gen != genJets->end() && jetInd<maxJets; ++ gen ) {
1689  p4gen[jetInd] = gen->p4(); //Gen 4-vector
1690  dRmin[jetInd] = 1000.0;
1691 
1692  nj = 0;
1693  usedInd = -1;
1694 
1695  for( CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++ cal ) {
1696  double delR = deltaR( cal->eta(), cal->phi(), gen->eta(), gen->phi() );
1697 
1698  if ( (delR<dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
1699  dRmin[jetInd] = delR; // delta R of match
1700  p4cal[jetInd] = cal->p4(); // Matched Cal 4-vector
1701  usedInd = nj;
1702  }
1703  nj++;
1704  }
1705 
1706  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1707  matchedAllPt31.Fill(p4gen[jetInd].pt());
1708  }
1709  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1710  matchedAllPt32.Fill(p4gen[jetInd].pt());
1711  }
1712  if (fabs(p4gen[jetInd].eta()) > 3.) {
1713  matchedAllPt33.Fill(p4gen[jetInd].pt());
1714  }
1715 
1716 
1717  if (usedInd != -1) {
1718  used[usedInd] = 1;
1719 
1720  if (p4cal[jetInd].pt() > minJetPt10)
1721  hf_PtResponse3.Fill(p4cal[jetInd].eta(), p4cal[jetInd].pt()/p4gen[jetInd].pt());
1722 
1723  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1724  matchedPt31.Fill(p4gen[jetInd].pt());
1725  }
1726  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1727  matchedPt32.Fill(p4gen[jetInd].pt());
1728  }
1729  if (fabs(p4gen[jetInd].eta()) > 3.) {
1730  matchedPt33.Fill(p4gen[jetInd].pt());
1731  }
1732  /***
1733  cout << " >>>MidPoint "
1734  << jetInd << " "
1735  << p4gen[jetInd].pt() << " "
1736  << p4gen[jetInd].eta() << " "
1737  << p4gen[jetInd].phi() << " "
1738  << p4cal[jetInd].pt() << " "
1739  << p4cal[jetInd].eta() << " "
1740  << p4cal[jetInd].phi()
1741  << endl;
1742  ***/
1743 
1744  double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1745  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1746  if ( (p4gen[jetInd].pt() > 20) && (p4gen[jetInd].pt() < 40) ) {
1747  dPt20Frac3.Fill(dpt/p4gen[jetInd].pt());
1748  }
1749  if ( (p4gen[jetInd].pt() > 40) && (p4gen[jetInd].pt() < 60) ) {
1750  dPt40Frac3.Fill(dpt/p4gen[jetInd].pt());
1751  }
1752  if ( (p4gen[jetInd].pt() > 80) && (p4gen[jetInd].pt() < 100) ) {
1753  dPt80Frac3.Fill(dpt/p4gen[jetInd].pt());
1754  }
1755  if ( (p4gen[jetInd].pt() > 100) && (p4gen[jetInd].pt() < 120) ) {
1756  dPt100Frac3.Fill(dpt/p4gen[jetInd].pt());
1757  }
1758  }
1759 
1760  if ( (p4gen[jetInd].pt() > minJetPt10) && (fabs(p4gen[jetInd].eta()) < 1.3) ) {
1761 
1762  dR3.Fill(dRmin[jetInd]);
1763  double dphi = deltaPhi(p4cal[jetInd].phi(), p4gen[jetInd].phi());
1764  dPhi3.Fill(dphi);
1765  double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1766  dEta3.Fill(deta);
1767  dPt3.Fill(dpt);
1768  dPtFrac3.Fill(dpt/p4gen[jetInd].pt());
1769 
1770 
1771  }
1772 
1773  jetInd++;
1774  }
1775  }
1776 
1777  // *********************
1778  // Jet Algorithm 4
1779 
1780  evt.getByLabel( GenJetAlgorithm4, genJets );
1781  evt.getByLabel( CaloJetAlgorithm4, caloJets );
1782 
1783  // int maxJets = 20;
1784  jetInd = 0;
1785  // double dRmin[20];
1786  // math::XYZTLorentzVector p4jet[20], p4gen[20], p4cal[20], p4cor[20];
1787 
1788  for( int i=0; i<maxJets; ++i ) used[i] = 0;
1789  for( GenJetCollection::const_iterator gen = genJets->begin();
1790  gen != genJets->end() && jetInd<maxJets; ++ gen ) {
1791  p4gen[jetInd] = gen->p4(); //Gen 4-vector
1792  dRmin[jetInd] = 1000.0;
1793 
1794  nj = 0;
1795  usedInd = -1;
1796 
1797  for( CaloJetCollection::const_iterator cal = caloJets->begin();
1798  cal != caloJets->end(); ++ cal ) {
1799  double delR = deltaR( cal->eta(), cal->phi(), gen->eta(), gen->phi() );
1800 
1801  if ( (delR<dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
1802  dRmin[jetInd] = delR; // delta R of match
1803  p4cal[jetInd] = cal->p4(); // Matched Cal 4-vector
1804  usedInd = nj;
1805  }
1806  nj++;
1807  }
1808 
1809  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1810  matchedAllPt41.Fill(p4gen[jetInd].pt());
1811  }
1812  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1813  matchedAllPt42.Fill(p4gen[jetInd].pt());
1814  }
1815  if (fabs(p4gen[jetInd].eta()) > 3.) {
1816  matchedAllPt43.Fill(p4gen[jetInd].pt());
1817  }
1818 
1819 
1820  if (usedInd != -1) {
1821 
1822  used[usedInd] = 1;
1823 
1824  if (p4cal[jetInd].pt() > minJetPt10)
1825  hf_PtResponse4.Fill(p4cal[jetInd].eta(), p4cal[jetInd].pt()/p4gen[jetInd].pt());
1826 
1827  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1828  matchedPt41.Fill(p4gen[jetInd].pt());
1829  }
1830  if ((fabs(p4gen[jetInd].eta()) > 1.3) && (fabs(p4gen[jetInd].eta()) < 3.) ) {
1831  matchedPt42.Fill(p4gen[jetInd].pt());
1832  }
1833  if (fabs(p4gen[jetInd].eta()) > 3.) {
1834  matchedPt43.Fill(p4gen[jetInd].pt());
1835  }
1836 
1837  /***
1838  cout << " >>>KtClus "
1839  << jetInd << " "
1840  << p4gen[jetInd].pt() << " "
1841  << p4gen[jetInd].eta() << " "
1842  << p4gen[jetInd].phi() << " "
1843  << p4cal[jetInd].pt() << " "
1844  << p4cal[jetInd].eta() << " "
1845  << p4cal[jetInd].phi()
1846  << endl;
1847  ***/
1848 
1849 
1850  double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1851  if (fabs(p4gen[jetInd].eta()) < 1.3) {
1852  if ( (p4gen[jetInd].pt() > 20) && (p4gen[jetInd].pt() < 40) ) {
1853  dPt20Frac4.Fill(dpt/p4gen[jetInd].pt());
1854  }
1855  if ( (p4gen[jetInd].pt() > 40) && (p4gen[jetInd].pt() < 60) ) {
1856  dPt40Frac4.Fill(dpt/p4gen[jetInd].pt());
1857  }
1858  if ( (p4gen[jetInd].pt() > 80) && (p4gen[jetInd].pt() < 100) ) {
1859  dPt80Frac4.Fill(dpt/p4gen[jetInd].pt());
1860  }
1861  if ( (p4gen[jetInd].pt() > 100) && (p4gen[jetInd].pt() < 120) ) {
1862  dPt100Frac4.Fill(dpt/p4gen[jetInd].pt());
1863  }
1864  }
1865 
1866  if ( (p4gen[jetInd].pt() > minJetPt10) && (fabs(p4gen[jetInd].eta()) < 1.3) ) {
1867 
1868  dR4.Fill(dRmin[jetInd]);
1869  double dphi = deltaPhi(p4cal[jetInd].phi(), p4gen[jetInd].phi());
1870  dPhi4.Fill(dphi);
1871  double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1872  dEta4.Fill(deta);
1873  dPt4.Fill(dpt);
1874  dPtFrac4.Fill(dpt/p4gen[jetInd].pt());
1875 
1876  }
1877 
1878  jetInd++;
1879  }
1880 
1881  }
1882 
1883 
1884 
1885  // *********************
1886  // MidPoint - Seedless Jet Matching
1887 
1888  Handle<CaloJetCollection> calo1Jets;
1889  Handle<CaloJetCollection> calo2Jets;
1890  Handle<CaloJetCollection> calo3Jets;
1891 
1892  evt.getByLabel( CaloJetAlgorithm1, calo1Jets );
1893  evt.getByLabel( CaloJetAlgorithm2, calo2Jets );
1894  evt.getByLabel( CaloJetAlgorithm3, calo3Jets );
1895 
1896  jetInd = 0;
1897 
1898  for( int i=0; i<maxJets; ++i ) used[i] = 0;
1899  for( CaloJetCollection::const_iterator cal1 = calo1Jets->begin();
1900  cal1 != calo1Jets->end() && jetInd<maxJets; ++cal1 ) {
1901 
1902  p4gen[jetInd] = cal1->p4(); //Gen 4-vector
1903  dRmin[jetInd] = 1000.0;
1904 
1905  nj = 0;
1906  for( CaloJetCollection::const_iterator cal2 = calo2Jets->begin(); cal2 != calo2Jets->end(); ++cal2 ) {
1907 
1908  double delR = deltaR( cal1->eta(), cal1->phi(), cal2->eta(), cal2->phi() );
1909  if ( (delR<dRmin[jetInd]) && (used[nj] == 0) ) {
1910  dRmin[jetInd] = delR; // delta R of match
1911  p4cal[jetInd] = cal2->p4(); // Matched Cal 4-vector
1912  usedInd = nj;
1913  }
1914  nj++;
1915  }
1916  used[usedInd] = 1;
1917 
1918  if (p4gen[jetInd].pt() > minJetPt) {
1919  dR12.Fill(dRmin[jetInd]);
1920  double dphi = deltaPhi(p4cal[jetInd].phi(), p4gen[jetInd].phi());
1921  dPhi12.Fill(dphi);
1922  double deta = p4cal[jetInd].eta() - p4gen[jetInd].eta();
1923  dEta12.Fill(deta);
1924  double dpt = p4cal[jetInd].pt() - p4gen[jetInd].pt();
1925  dPt12.Fill(dpt);
1926  }
1927 
1928  jetInd++;
1929  }
1930 
1931  // ******************************************
1932  // ******************************************
1933 
1934 
1935  Handle<CandidateCollection> genParticles;
1936  evt.getByLabel("genParticleCandidates",genParticles);
1937 
1938 
1939  // *********************
1940  // Partons (Z')
1941 
1942  int nPart = 0;
1943  for (size_t i =0;i< genParticles->size(); i++) {
1944 
1945  const Candidate &p = (*genParticles)[i];
1946  // int Status = p.status();
1947  // bool ParticleIsStable = Status==1;
1948  int id = p.pdgId();
1949 
1950  if (id == 32) {
1951 
1952  if (p.numberOfDaughters() != 0) {
1953  /***
1954  cout << "Z': part = " << i << " id = " << id
1955  << " daughters = " << p.numberOfDaughters()
1956  << " mass = " << p.mass()
1957  << endl;
1958  ***/
1959  ZpMG = p.mass();
1960  ZpMassGen.Fill( ZpMG );
1961  if (EtaOk10 == 2) {
1962  ZpMassGen10.Fill( ZpMG );
1963  ZpMassGen_700_10.Fill( ZpMG );
1964  ZpMassGen_2000_10.Fill( ZpMG );
1965  ZpMassGen_5000_10.Fill( ZpMG );
1966  }
1967  if (EtaOk13 == 2) {
1968  ZpMassGen13.Fill( ZpMG );
1969  ZpMassGen_700_13.Fill( ZpMG );
1970  ZpMassGen_2000_13.Fill( ZpMG );
1971  ZpMassGen_5000_13.Fill( ZpMG );
1972  }
1973  if (EtaOk40 == 2) {
1974  ZpMassGen40.Fill( ZpMG );
1975  ZpMassGen_700_40.Fill( ZpMG );
1976  ZpMassGen_2000_40.Fill( ZpMG );
1977  ZpMassGen_5000_40.Fill( ZpMG );
1978  }
1979  }
1980 
1981  for( int id1=0, nd1=p.numberOfDaughters(); id1 < nd1; ++id1 ) {
1982 
1983  const Candidate * d1 = p.daughter(id1);
1984 
1985  if ( abs(d1->pdgId()) != 32) {
1986  math::XYZTLorentzVector momentum=d1->p4();
1987  p4Zp[nPart] = momentum=d1->p4();
1988  nPart++;
1989  }
1990 
1991  }
1992  }
1993 
1994  }
1995 
1996  // *********************
1997  // Match jets to Zp
1998  int genInd = 0;
1999 
2000  if (nPart == 2) {
2001 
2002  ZpM = (p4Zp[0]+p4Zp[1]).mass();
2003  ZpMass.Fill( ZpM );
2004 
2005  if (EtaOk10 == 2) {
2006  ZpMass_700_10.Fill( ZpM );
2007  ZpMass_2000_10.Fill( ZpM );
2008  ZpMass_5000_10.Fill( ZpM );
2009  }
2010  if (EtaOk13 == 2) {
2011  ZpMass_700_13.Fill( ZpM );
2012  ZpMass_2000_13.Fill( ZpM );
2013  ZpMass_5000_13.Fill( ZpM );
2014  }
2015  if (EtaOk40 == 2) {
2016  ZpMass_700_40.Fill( ZpM );
2017  ZpMass_2000_40.Fill( ZpM );
2018  ZpMass_5000_40.Fill( ZpM );
2019  }
2020 
2021  int usedInd;
2022 
2023  // ***********
2024  // **** Calor1
2025  usedInd = -1;
2026  jetInd = 0;
2027 
2028  for( int i=0; i<maxJets; ++i ) used[i] = 0;
2029  for( int i=0; i<2; ++i ) {
2030 
2031  dRmin[jetInd] = 1000.0;
2032 
2033  int nj = 0;
2034  for( CaloJetCollection::const_iterator cal1 = calo1Jets->begin();
2035  cal1 != calo1Jets->end() && jetInd<maxJets; ++cal1 ) {
2036 
2037  double delR = deltaR( cal1->eta(), cal1->phi(), p4Zp[i].eta(), p4Zp[i].phi() );
2038 
2039  // if ( (delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
2040  if ( (delR < dRmin[jetInd]) && (used[nj] == 0) ) {
2041  dRmin[jetInd] = delR; // delta R of match
2042  p4cal[jetInd] = cal1->p4(); // Matched Cal 4-vector
2043  usedInd = nj;
2044  genInd = i;
2045  }
2046 
2047  /****
2048  cout << "Delta R = " << delR
2049  << " deltaR min = " << dRmin[jetInd]
2050  << " Ind = " << jetInd << " / " << nj << " / " << used[nj]
2051  << " cal1 eta = " << cal1->eta()
2052  << " p4par eta = " << p4Zp[i].eta()
2053  << " cal1 phi = " << cal1->phi()
2054  << " p4par phi = " << p4Zp[i].phi()
2055  << endl;
2056  cout << " "
2057  << " p4par = " << p4Zp[i]
2058  << " p4cal = " << cal1->p4()
2059  << endl;
2060  ***/
2061 
2062  nj++;
2063  }
2064 
2065  // Found matched jet
2066  if (usedInd != -1) {
2067  used[usedInd] = 1;
2068  jetInd++;
2069  }
2070 
2071  }
2072 
2073  ZpMM = (p4cal[0]+p4cal[1]).mass();
2074  ZpMassMatched1.Fill( ZpMM );
2075 
2076  if ((ZpMG != 0) && (EtaOk40 == 2)) {
2077  ZpMassRes401.Fill( (ZpMM - ZpMG) / ZpMG );
2078 
2079  ZpMassResL401.Fill( (LeadMass1 - ZpMG) / ZpMG );
2080  ZpMassResL402.Fill( (LeadMass2 - ZpMG) / ZpMG );
2081  ZpMassResL403.Fill( (LeadMass3 - ZpMG) / ZpMG );
2082 
2083  ZpMassResRL401.Fill( LeadMass1 / ZpMG );
2084  ZpMassResRL402.Fill( LeadMass2 / ZpMG );
2085  ZpMassResRL403.Fill( LeadMass3 / ZpMG );
2086 
2087  ZpMassResRLoP401.Fill( LeadMass1 / LeadMassP1 );
2088  ZpMassResRLoP402.Fill( LeadMass2 / LeadMassP2 );
2089  ZpMassResRLoP403.Fill( LeadMass3 / LeadMassP2 );
2090 
2091  ZpMassResPRL401.Fill( LeadMassP1 / ZpMG );
2092  ZpMassResPRL402.Fill( LeadMassP2 / ZpMG );
2093  ZpMassResPRL403.Fill( LeadMassP3 / ZpMG );
2094 
2095  }
2096 
2097  if ((ZpMG != 0) && (EtaOk10 == 2)) {
2098  ZpMassRes101.Fill( (ZpMM - ZpMG) / ZpMG );
2099 
2100  ZpMassResL101.Fill( (LeadMass1 - ZpMG) / ZpMG );
2101  ZpMassResL102.Fill( (LeadMass2 - ZpMG) / ZpMG );
2102  ZpMassResL103.Fill( (LeadMass3 - ZpMG) / ZpMG );
2103 
2104  ZpMassResRL101.Fill( LeadMass1 / ZpMG );
2105  ZpMassResRL102.Fill( LeadMass2 / ZpMG );
2106  ZpMassResRL103.Fill( LeadMass3 / ZpMG );
2107 
2108  ZpMassResRLoP101.Fill( LeadMass1 / LeadMassP1 );
2109  ZpMassResRLoP102.Fill( LeadMass2 / LeadMassP2 );
2110  ZpMassResRLoP103.Fill( LeadMass3 / LeadMassP2 );
2111 
2112  ZpMassResPRL101.Fill( LeadMassP1 / ZpMG );
2113  ZpMassResPRL102.Fill( LeadMassP2 / ZpMG );
2114  ZpMassResPRL103.Fill( LeadMassP3 / ZpMG );
2115 
2116  }
2117 
2118  if ((ZpMG != 0) && (EtaOk13 == 2)) {
2119  ZpMassRes131.Fill( (ZpMM - ZpMG) / ZpMG );
2120 
2121  ZpMassResL131.Fill( (LeadMass1 - ZpMG) / ZpMG );
2122  ZpMassResL132.Fill( (LeadMass2 - ZpMG) / ZpMG );
2123  ZpMassResL133.Fill( (LeadMass3 - ZpMG) / ZpMG );
2124 
2125  ZpMassResRL131.Fill( LeadMass1 / ZpMG );
2126  ZpMassResRL132.Fill( LeadMass2 / ZpMG );
2127  ZpMassResRL133.Fill( LeadMass3 / ZpMG );
2128 
2129  ZpMassResRLoP131.Fill( LeadMass1 / LeadMassP1 );
2130  ZpMassResRLoP132.Fill( LeadMass2 / LeadMassP2 );
2131  ZpMassResRLoP133.Fill( LeadMass3 / LeadMassP2 );
2132 
2133  ZpMassResPRL131.Fill( LeadMassP1 / ZpMG );
2134  ZpMassResPRL132.Fill( LeadMassP2 / ZpMG );
2135  ZpMassResPRL133.Fill( LeadMassP3 / ZpMG );
2136 
2137  }
2138 
2139 
2140 
2141  // ***********
2142  // **** Calor2
2143  usedInd = -1;
2144  jetInd = 0;
2145 
2146  for( int i=0; i<maxJets; ++i ) used[i] = 0;
2147  for( int i=0; i<2; ++i ) {
2148 
2149  dRmin[jetInd] = 1000.0;
2150 
2151  int nj = 0;
2152  for( CaloJetCollection::const_iterator cal2 = calo2Jets->begin();
2153  cal2 != calo2Jets->end() && jetInd<maxJets; ++cal2 ) {
2154 
2155  double delR = deltaR( cal2->eta(), cal2->phi(), p4Zp[i].eta(), p4Zp[i].phi() );
2156 
2157  if ( (delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
2158  dRmin[jetInd] = delR; // delta R of match
2159  p4cal[jetInd] = cal2->p4(); // Matched Cal 4-vector
2160  usedInd = nj;
2161  }
2162 
2163  /****
2164  cout << "Delta R = " << delR
2165  << " deltaR min = " << dRmin[jetInd]
2166  << " Ind = " << jetInd << " / " << nj << " / " << used[nj]
2167  << " p4par = " << p4par[i]
2168  << " p4cal = " << cal1->p4()
2169  << " cal1 eta = " << cal1->eta()
2170  << " p4par eta = " << p4par[i].eta()
2171  << endl;
2172  ****/
2173 
2174  nj++;
2175  }
2176 
2177  // Found matched jet
2178  if (usedInd != -1) {
2179  used[usedInd] = 1;
2180  jetInd++;
2181  }
2182 
2183  }
2184 
2185  ZpMM = (p4cal[0]+p4cal[1]).mass();
2186  ZpMassMatched2.Fill( ZpMM );
2187  ZpMassRes402.Fill( (ZpMM - ZpM) / ZpM );
2188 
2189 
2190  // ***********
2191  // **** Calor3
2192  usedInd = -1;
2193  jetInd = 0;
2194 
2195  for( int i=0; i<maxJets; ++i ) used[i] = 0;
2196  for( int i=0; i<2; ++i ) {
2197 
2198  dRmin[jetInd] = 1000.0;
2199 
2200  int nj = 0;
2201  for( CaloJetCollection::const_iterator cal3 = calo3Jets->begin();
2202  cal3 != calo3Jets->end() && jetInd<maxJets; ++cal3 ) {
2203 
2204  double delR = deltaR( cal3->eta(), cal3->phi(), p4Zp[i].eta(), p4Zp[i].phi() );
2205 
2206  if ( (delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
2207  dRmin[jetInd] = delR; // delta R of match
2208  p4cal[jetInd] = cal3->p4(); // Matched Cal 4-vector
2209  usedInd = nj;
2210  }
2211 
2212  /****
2213  cout << "Delta R = " << delR
2214  << " deltaR min = " << dRmin[jetInd]
2215  << " Ind = " << jetInd << " / " << nj << " / " << used[nj]
2216  << " p4par = " << p4par[i]
2217  << " p4cal = " << cal1->p4()
2218  << " cal1 eta = " << cal1->eta()
2219  << " p4par eta = " << p4par[i].eta()
2220  << endl;
2221  ****/
2222 
2223  nj++;
2224  }
2225 
2226 
2227  // Found matched jet
2228  if (usedInd != -1) {
2229  used[usedInd] = 1;
2230  jetInd++;
2231  }
2232 
2233  }
2234 
2235  ZpMM = (p4cal[0]+p4cal[1]).mass();
2236  ZpMassMatched3.Fill( ZpMM );
2237  ZpMassRes403.Fill( (ZpMM - ZpM) / ZpM );
2238 
2239  } else {
2240  cout << "Z' (3): nPart = " << nPart << endl;
2241  }
2242 
2243 
2244 
2245  // *********************
2246  // Partons (ttbar) Jet Matching
2247 
2248  // cout << ">>> Begin MC list" << endl;
2249  int nJet = 0;
2250 
2251 
2252  int ii = 1;
2253  int jj = 4;
2254 
2255  for (size_t i =0;i< genParticles->size(); i++) {
2256 
2257  const Candidate &p = (*genParticles)[i];
2258  // int Status = p.status();
2259  // bool ParticleIsStable = Status==1;
2260  int id = p.pdgId();
2261 
2262  // Top Quarks
2263  if (abs(id) == 6) {
2264  cout << "TOP: id = " << id << " mass = " << p.mass() << endl;
2265 
2266  topMassParton.Fill(p.mass());
2267 
2268  if (id == 6) tMassGen.Fill(p.mass());
2269  if (id == -6) tbarMassGen.Fill(p.mass());
2270 
2271  for( size_t id1=0, nd1=p.numberOfDaughters(); id1 < nd1; ++id1 ) {
2272 
2273  const Candidate * d1 = p.daughter(id1);
2274 
2275  // b - quark
2276  if ( abs(d1->pdgId()) == 5) {
2277 
2278  math::XYZTLorentzVector momentum=d1->p4();
2279  p4par[nJet++] = momentum=d1->p4();
2280 
2281  cout << "Daughter1: id = " << d1->pdgId()
2282  << " daughters = " << d1->numberOfDaughters()
2283  << " mother 1 = " << (d1->mother())->pdgId()
2284  << " Momentum " << momentum << " GeV/c"
2285  << endl;
2286 
2287 
2288  if ( (d1->mother())->pdgId() == 6 ) {
2289  p4part[0] = momentum=d1->p4();
2290  cout << ">>> part0 = " << p4part[0] << endl;
2291  }
2292  if ( (d1->mother())->pdgId() == -6) {
2293  p4part[3] = momentum=d1->p4();
2294  cout << ">>> part3 = " << p4part[3] << endl;
2295  }
2296 
2297  }
2298 
2299 
2300  // W
2301  // Check for fully hadronic decay
2302 
2303  if ( abs(d1->pdgId()) == 24) {
2304 
2305  for( size_t id2=0, nd2=d1->numberOfDaughters(); id2 < nd2; ++id2 ) {
2306 
2307  const Candidate * d2 = d1->daughter(id2);
2308 
2309  if (abs(d2->pdgId()) < 9) {
2310 
2311  math::XYZVector vertex(d2->vx(),d2->vy(),d2->vz());
2312  math::XYZTLorentzVector momentum=d2->p4();
2313  p4par[nJet++] = momentum=d2->p4();
2314 
2315  if ( (d1->mother())->pdgId() == 6 ) {
2316  p4part[ii] = momentum=d2->p4();
2317  cout << ">>> part" << ii << " = " << p4part[ii] << endl;
2318  ii++;
2319  }
2320  if ( (d1->mother())->pdgId() == -6 ) {
2321  p4part[jj] = momentum=d2->p4();
2322  cout << ">>> part" << jj << " = " << p4part[jj] << endl;
2323  jj++;
2324  }
2325 
2326  cout << "Daughter2: id = " << d2->pdgId()
2327  << " daughters = " << d2->numberOfDaughters()
2328  << " mother 2 = " << (d2->mother())->pdgId()
2329  << " Momentum " << momentum << " GeV/c"
2330  << endl;
2331  }
2332 
2333  }
2334  }
2335 
2336  // if ( pdgId == d->pdgId() && d->status() == 1 ) {
2337  // }
2338  }
2339  }
2340 
2341  }
2342  // cout << ">>> N Jets = " << nJet << endl;
2343 
2344  if (nJet == 6) {
2345 
2346  double tmass = (p4part[0]+p4part[1]+p4part[2]).mass();
2347  double tbarmass = (p4part[3]+p4part[4]+p4part[5]).mass();
2348 
2349  tMass.Fill(tmass);
2350  tbarMass.Fill(tbarmass);
2351 
2352  cout << ">>> T Mass = " << tmass << " / " << tbarmass << endl;
2353 
2354  double mindR = 1000.;
2355  for( size_t i=0; i<6; ++i ) {
2356  for( size_t j=0; j<6; ++j ) {
2357  if (j > i) {
2358  double delR = deltaR( p4par[i].eta(), p4par[i].phi(), p4par[j].eta(), p4par[j].phi() );
2359  if (delR < mindR) mindR = delR;
2360  dRParton.Fill(delR);
2361  }
2362  }
2363  }
2364  dRPartonMin.Fill(mindR);
2365 
2366  int usedInd;
2367  usedInd = -1;
2368  jetInd = 0;
2369 
2370  for( int i=0; i<maxJets; ++i ) used[i] = 0;
2371  for( int i=0; i<6; ++i ) {
2372 
2373  dRmin[jetInd] = 1000.0;
2374 
2375  int nj = 0;
2376  for( CaloJetCollection::const_iterator cal1 = calo1Jets->begin();
2377  cal1 != calo1Jets->end() && jetInd<maxJets; ++cal1 ) {
2378 
2379  double delR = deltaR( cal1->eta(), cal1->phi(), p4par[i].eta(), p4par[i].phi() );
2380 
2381  if ( (delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
2382  dRmin[jetInd] = delR; // delta R of match
2383  p4cal[jetInd] = cal1->p4(); // Matched Cal 4-vector
2384  usedInd = nj;
2385  genInd = i;
2386  }
2387 
2388  /****
2389  cout << "Delta R = " << delR
2390  << " deltaR min = " << dRmin[jetInd]
2391  << " Ind = " << jetInd << " / " << nj << " / " << used[nj]
2392  << " p4par = " << p4par[i]
2393  << " p4cal = " << cal1->p4()
2394  << " cal1 eta = " << cal1->eta()
2395  << " p4par eta = " << p4par[i].eta()
2396  << endl;
2397  ****/
2398 
2399 
2400  nj++;
2401  }
2402 
2403 
2404  // Found matched jet
2405  if (usedInd != -1) {
2406  used[usedInd] = 1;
2407 
2408  dRPar1.Fill(dRmin[jetInd]);
2409  double dphi = deltaPhi(p4cal[jetInd].phi(), p4par[genInd].phi());
2410  dPhiPar1.Fill(dphi);
2411  double deta = p4cal[jetInd].eta() - p4par[genInd].eta();
2412  dEtaPar1.Fill(deta);
2413  double dpt = p4cal[jetInd].pt() - p4par[genInd].pt();
2414  dPtPar1.Fill(dpt);
2415  jetInd++;
2416  }
2417 
2418  }
2419  ParMatch1.Fill(jetInd);
2420  if (jetInd == 6) {
2421  topMass1.Fill( (p4cal[0]+p4cal[1]+p4cal[2]).mass() );
2422  topMass1.Fill( (p4cal[3]+p4cal[4]+p4cal[5]).mass() );
2423  }
2424 
2425  /***
2426  cout << "Collection Size = " << calo1Jets->size()
2427  << " / " << jetInd
2428  << endl;
2429  ***/
2430 
2431  // ***********************
2432  jetInd = 0;
2433  usedInd = -1;
2434 
2435  for( int i=0; i<maxJets; ++i ) used[i] = 0;
2436  for( int i=0; i<6; ++i ) {
2437 
2438  dRmin[jetInd] = 1000.0;
2439 
2440  int nj = 0;
2441  for( CaloJetCollection::const_iterator cal2 = calo2Jets->begin();
2442  cal2 != calo2Jets->end() && jetInd<maxJets; ++cal2 ) {
2443 
2444  double delR = deltaR( cal2->eta(), cal2->phi(), p4par[i].eta(), p4par[i].phi() );
2445 
2446  if ( (delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
2447  dRmin[jetInd] = delR; // delta R of match
2448  p4cal[jetInd] = cal2->p4(); // Matched Cal 4-vector
2449  usedInd = nj;
2450  genInd = i;
2451  }
2452 
2453  /****
2454  cout << "Delta R = " << delR
2455  << " deltaR min = " << dRmin[jetInd]
2456  << " Ind = " << jetInd << " / " << nj << " / " << used[nj]
2457  << " cal2 eta = " << cal2->eta()
2458  << " p4par eta = " << p4par[i].eta()
2459  << endl;
2460  ****/
2461 
2462  nj++;
2463  }
2464  if (usedInd != -1) {
2465  used[usedInd] = 1;
2466 
2467  dRPar2.Fill(dRmin[jetInd]);
2468  double dphi = deltaPhi(p4cal[jetInd].phi(), p4par[genInd].phi());
2469  dPhiPar2.Fill(dphi);
2470  double deta = p4cal[jetInd].eta() - p4par[genInd].eta();
2471  dEtaPar2.Fill(deta);
2472  double dpt = p4cal[jetInd].pt() - p4par[genInd].pt();
2473  dPtPar2.Fill(dpt);
2474 
2475  jetInd++;
2476  }
2477 
2478  }
2479  ParMatch2.Fill(jetInd);
2480  if (jetInd == 6) {
2481  topMass2.Fill( (p4cal[0]+p4cal[1]+p4cal[2]).mass() );
2482  topMass2.Fill( (p4cal[3]+p4cal[4]+p4cal[5]).mass() );
2483  }
2484 
2485 
2486  /***
2487  cout << "Collection Size = " << calo2Jets->size()
2488  << " / " << jetInd
2489  << endl;
2490  ***/
2491 
2492 
2493  // ***********************
2494  jetInd = 0;
2495  usedInd = -1;
2496 
2497  for( int i=0; i<maxJets; ++i ) used[i] = 0;
2498  for( int i=0; i<6; ++i ) {
2499 
2500  dRmin[jetInd] = 1000.0;
2501 
2502  int nj = 0;
2503  for( CaloJetCollection::const_iterator cal3 = calo3Jets->begin();
2504  cal3 != calo3Jets->end() && jetInd<maxJets; ++cal3 ) {
2505 
2506  double delR = deltaR( cal3->eta(), cal3->phi(), p4par[i].eta(), p4par[i].phi() );
2507 
2508  if ( (delR < dRmin[jetInd]) && (delR < matchedDelR) && (used[nj] == 0) ) {
2509  dRmin[jetInd] = delR; // delta R of match
2510  p4cal[jetInd] = cal3->p4(); // Matched Cal 4-vector
2511  usedInd = nj;
2512  genInd = i;
2513  }
2514  /****
2515  cout << "Delta R = " << delR
2516  << " deltaR min = " << dRmin[jetInd]
2517  << " Ind = " << jetInd << " / " << nj << " / " << used[nj]
2518  << " cal3 eta = " << cal3->eta()
2519  << " p4par eta = " << p4par[i].eta()
2520  << endl;
2521  ****/
2522 
2523  nj++;
2524  }
2525  if (usedInd != -1) {
2526  used[usedInd] = 1;
2527 
2528  dRPar3.Fill(dRmin[jetInd]);
2529  double dphi = deltaPhi(p4cal[jetInd].phi(), p4par[genInd].phi());
2530  dPhiPar3.Fill(dphi);
2531  double deta = p4cal[jetInd].eta() - p4par[genInd].eta();
2532  dEtaPar3.Fill(deta);
2533  double dpt = p4cal[jetInd].pt() - p4par[genInd].pt();
2534  dPtPar3.Fill(dpt);
2535 
2536  jetInd++;
2537  }
2538  }
2539  ParMatch3.Fill(jetInd);
2540  if (jetInd == 6) {
2541  topMass3.Fill( (p4cal[0]+p4cal[1]+p4cal[2]).mass() );
2542  topMass3.Fill( (p4cal[3]+p4cal[4]+p4cal[5]).mass() );
2543  }
2544 
2545  /***
2546  cout << "Collection Size = " << calo3Jets->size()
2547  << " / " << jetInd
2548  << endl;
2549  ***/
2550 
2551  }
2552 
2553  int nTow1, nTow2, nTow3, nTow4;
2554 
2555  Handle<CaloJetCollection> jets;
2556 
2557  // *********************
2558  // Jet Properties
2559  // *********************
2560 
2561 
2562 
2563  // --- Loop over jets and make a list of all the used towers
2564  evt.getByLabel( CaloJetAlgorithm1, jets );
2565  int jjet = 0;
2566  for ( CaloJetCollection::const_iterator ijet=jets->begin(); ijet!=jets->end(); ijet++) {
2567  jjet++;
2568 
2569  float hadEne = ijet->hadEnergyInHB() + ijet->hadEnergyInHO() +
2570  ijet->hadEnergyInHE() + ijet->hadEnergyInHF();
2571  float emEne = ijet->emEnergyInEB() + ijet->emEnergyInEE() + ijet->emEnergyInHF();
2572  float had = ijet->energyFractionHadronic();
2573 
2574  float j_et = ijet->et();
2575 
2576  if (fabs(ijet->eta()) < 1.3) {
2577  totEneLeadJetEta1_1.Fill(hadEne+emEne);
2578  hadEneLeadJetEta1_1.Fill(hadEne);
2579  emEneLeadJetEta1_1.Fill(emEne);
2580 
2581  totEneLeadJetEta1_2.Fill(hadEne+emEne);
2582  hadEneLeadJetEta1_2.Fill(hadEne);
2583  emEneLeadJetEta1_2.Fill(emEne);
2584 
2585  if (ijet->pt() > minJetPt10)
2586  hadFracEta11.Fill(had);
2587  }
2588  if ((fabs(ijet->eta()) > 1.3) && (fabs(ijet->eta()) < 3.) ) {
2589 
2590  totEneLeadJetEta2_1.Fill(hadEne+emEne);
2591  hadEneLeadJetEta2_1.Fill(hadEne);
2592  emEneLeadJetEta2_1.Fill(emEne);
2593 
2594  totEneLeadJetEta2_2.Fill(hadEne+emEne);
2595  hadEneLeadJetEta2_2.Fill(hadEne);
2596  emEneLeadJetEta2_2.Fill(emEne);
2597 
2598  if (ijet->pt() > minJetPt10)
2599  hadFracEta21.Fill(had);
2600  }
2601  if (fabs(ijet->eta()) > 3.) {
2602 
2603  totEneLeadJetEta3_1.Fill(hadEne+emEne);
2604  hadEneLeadJetEta3_1.Fill(hadEne);
2605  emEneLeadJetEta3_1.Fill(emEne);
2606 
2607  totEneLeadJetEta3_2.Fill(hadEne+emEne);
2608  hadEneLeadJetEta3_2.Fill(hadEne);
2609  emEneLeadJetEta3_2.Fill(emEne);
2610 
2611  if (ijet->pt() > minJetPt10)
2612  hadFracEta31.Fill(had);
2613  }
2614 
2615 
2616  if (jjet == 1) {
2617  hadFracLeadJet1.Fill(had);
2618  hadEneLeadJet1.Fill(hadEne);
2619  hadEneLeadJet12.Fill(hadEne);
2620  hadEneLeadJet13.Fill(hadEne);
2621  emEneLeadJet1.Fill(emEne);
2622  emEneLeadJet12.Fill(emEne);
2623  emEneLeadJet13.Fill(emEne);
2624  }
2625 
2626  const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2627  int nConstituents = jetCaloRefs.size();
2628 
2629  if (jjet == 1) {
2630 
2631  nTow1 = nTow2 = nTow3 = nTow4 = 0;
2632  for (int i = 0; i <nConstituents ; i++){
2633 
2634  float et = jetCaloRefs[i]->et();
2635 
2636  if (et > 0.5) nTow1++;
2637  if (et > 1.0) nTow2++;
2638  if (et > 1.5) nTow3++;
2639  if (et > 2.0) nTow4++;
2640 
2641  hf_TowerJetEt1.Fill(et/j_et);
2642 
2643  }
2644 
2645  nTowersLeadJetPt1.Fill(nTow1);
2646  nTowersLeadJetPt2.Fill(nTow2);
2647  nTowersLeadJetPt3.Fill(nTow3);
2648  nTowersLeadJetPt4.Fill(nTow4);
2649 
2650  }
2651 
2652  if ( (jjet == 1) && (fabs(ijet->eta()) < 1.3) ) {
2653  nTowersLeadJet1.Fill( nConstituents );
2654 
2655  for (int i = 0; i <nConstituents ; i++){
2656  float t_et = jetCaloRefs[i]->et();
2657  double delR = deltaR( ijet->eta(), ijet->phi(), jetCaloRefs[i]->eta(), jetCaloRefs[i]->phi() );
2658  hf_TowerDelR1.Fill( delR, t_et/j_et);
2659  hf_TowerDelR12.Fill( delR, t_et/j_et);
2660  TowerEtLeadJet1.Fill( t_et );
2661  TowerEtLeadJet12.Fill( t_et );
2662  TowerEtLeadJet13.Fill( t_et );
2663  }
2664  }
2665 
2666 
2667  if ( (jjet == 2) && (fabs(ijet->eta()) < 1.3) ) {
2668  nTowersSecondJet1.Fill( nConstituents );
2669  }
2670 
2671 
2672 
2673  }
2674 
2675  // *********************
2676  // Unclustered Energy
2677  // *********************
2678 
2679  double SumPtJet(0);
2680 
2681  double SumEtNotJets(0);
2682  double SumEtJets(0);
2683  double SumEtTowers(0);
2684 
2685  double sumJetPx(0);
2686  double sumJetPy(0);
2687 
2688  double sumTowerAllPx(0);
2689  double sumTowerAllPy(0);
2690 
2691  double sumTowerAllEx(0);
2692  double sumTowerAllEy(0);
2693 
2694  std::vector<CaloTowerPtr> UsedTowerList;
2695  std::vector<CaloTower> TowerUsedInJets;
2696  std::vector<CaloTower> TowerNotUsedInJets;
2697 
2698 
2699  // *********************
2700  // Towers
2701 
2702  Handle<CaloTowerCollection> caloTowers;
2703  evt.getByLabel( "towerMaker", caloTowers );
2704 
2705 
2706 
2707  nTow1 = nTow2 = nTow3 = nTow4 = 0;
2708 
2709  double sum_et = 0.0;
2710  double sum_ex = 0.0;
2711  double sum_ey = 0.0;
2712  // double sum_ez = 0.0;
2713 
2714  // --- Loop over towers and make a lists of used and unused towers
2715  for (CaloTowerCollection::const_iterator tower = caloTowers->begin();
2716  tower != caloTowers->end(); tower++) {
2717 
2718  Double_t et = tower->et();
2719 
2720  if (et > 0.5) nTow1++;
2721  if (et > 1.0) nTow2++;
2722  if (et > 1.5) nTow3++;
2723  if (et > 2.0) nTow4++;
2724 
2725  if(et>0.5) {
2726 
2727  // ********
2728  double phix = tower->phi();
2729  // double theta = tower->theta();
2730  // double e = tower->energy();
2731  // double et = e*sin(theta);
2732  // double et = tower->emEt() + tower->hadEt();
2733  double et = tower->et();
2734 
2735  // sum_ez += e*cos(theta);
2736  sum_et += et;
2737  sum_ex += et*cos(phix);
2738  sum_ey += et*sin(phix);
2739  // ********
2740 
2741  Double_t phi = tower->phi();
2742  SumEtTowers += tower->et();
2743 
2744  sumTowerAllEx += et*cos(phi);
2745  sumTowerAllEy += et*sin(phi);
2746 
2747  }
2748 
2749  }
2750 
2751  SumEt1.Fill(sum_et);
2752  SumEt12.Fill(sum_et);
2753  SumEt13.Fill(sum_et);
2754 
2755  MET1.Fill(sqrt( sum_ex*sum_ex + sum_ey*sum_ey));
2756  MET12.Fill(sqrt( sum_ex*sum_ex + sum_ey*sum_ey));
2757  MET13.Fill(sqrt( sum_ex*sum_ex + sum_ey*sum_ey));
2758 
2759  // met->mex = -sum_ex;
2760  // met->mey = -sum_ey;
2761  // met->mez = -sum_ez;
2762  // met->met = sqrt( sum_ex*sum_ex + sum_ey*sum_ey );
2763  // cout << "MET = " << met->met << endl;
2764  // met->sumet = sum_et;
2765  // met->phi = atan2( -sum_ey, -sum_ex );
2766 
2767 
2768 
2769  hf_sumTowerAllEx.Fill(sumTowerAllEx);
2770  hf_sumTowerAllEy.Fill(sumTowerAllEy);
2771 
2772  nTowers1.Fill(nTow1);
2773  nTowers2.Fill(nTow2);
2774  nTowers3.Fill(nTow3);
2775  nTowers4.Fill(nTow4);
2776 
2777  // *********************
2778  // MidPoint
2779  //
2780  UsedTowerList.clear();
2781  TowerUsedInJets.clear();
2782  TowerNotUsedInJets.clear();
2783 
2784  // --- Loop over jets and make a list of all the used towers
2785  evt.getByLabel( CaloJetAlgorithm1, jets );
2786  for ( CaloJetCollection::const_iterator ijet=jets->begin(); ijet!=jets->end(); ijet++) {
2787 
2788  Double_t jetPt = ijet->pt();
2789  Double_t jetPhi = ijet->phi();
2790 
2791  // if (jetPt>5.0) {
2792 
2793  Double_t jetPx = jetPt*cos(jetPhi);
2794  Double_t jetPy = jetPt*sin(jetPhi);
2795 
2796  sumJetPx +=jetPx;
2797  sumJetPy +=jetPy;
2798 
2799  const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2800  int nConstituents = jetCaloRefs.size();
2801  for (int i = 0; i <nConstituents ; i++){
2802  UsedTowerList.push_back(jetCaloRefs[i]);
2803  }
2804 
2805  SumPtJet +=jetPt;
2806  // }
2807 
2808  }
2809 
2810 
2811  int NTowersUsed = UsedTowerList.size();
2812 
2813  // --- Loop over towers and make a lists of used and unused towers
2814  for (CaloTowerCollection::const_iterator tower = caloTowers->begin();
2815  tower != caloTowers->end(); tower++) {
2816 
2817  CaloTower t = *tower;
2818  Double_t et = tower->et();
2819 
2820  if(et>0) {
2821 
2822  Double_t phi = tower->phi();
2823  SumEtTowers += tower->et();
2824 
2825  sumTowerAllPx += et*cos(phi);
2826  sumTowerAllPy += et*sin(phi);
2827 
2828  bool used = false;
2829 
2830  for(int i=0; i<NTowersUsed; i++){
2831  if(tower->id() == UsedTowerList[i]->id()){
2832  used=true;
2833  break;
2834  }
2835  }
2836 
2837  if (used) {
2838  TowerUsedInJets.push_back(t);
2839  } else {
2840  TowerNotUsedInJets.push_back(t);
2841  }
2842 
2843  }
2844 
2845  }
2846 
2847  int nUsed = TowerUsedInJets.size();
2848  int nNotUsed = TowerNotUsedInJets.size();
2849 
2850  SumEtJets = 0;
2851  SumEtNotJets = 0;
2852 
2853  for(int i=0;i<nUsed;i++){
2854  SumEtJets += TowerUsedInJets[i].et();
2855  }
2856  h_jetEt1.Fill(SumEtJets);
2857 
2858  for(int i=0;i<nNotUsed;i++){
2859  if (TowerNotUsedInJets[i].et() > 0.5)
2860  SumEtNotJets += TowerNotUsedInJets[i].et();
2861  h_missEt1.Fill(TowerNotUsedInJets[i].et());
2862  h_missEt1s.Fill(TowerNotUsedInJets[i].et());
2863  }
2864  h_totMissEt1.Fill(SumEtNotJets);
2865 
2866 
2867 
2868  // *********************
2869  // SISCone
2870  //
2871  UsedTowerList.clear();
2872  TowerUsedInJets.clear();
2873  TowerNotUsedInJets.clear();
2874 
2875  // --- Loop over jets and make a list of all the used towers
2876  evt.getByLabel( CaloJetAlgorithm2, jets );
2877  for ( CaloJetCollection::const_iterator ijet=jets->begin(); ijet!=jets->end(); ijet++) {
2878 
2879  Double_t jetPt = ijet->pt();
2880  Double_t jetPhi = ijet->phi();
2881 
2882  // if (jetPt>5.0) {
2883 
2884  Double_t jetPx = jetPt*cos(jetPhi);
2885  Double_t jetPy = jetPt*sin(jetPhi);
2886 
2887  sumJetPx +=jetPx;
2888  sumJetPy +=jetPy;
2889 
2890  const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2891  int nConstituents = jetCaloRefs.size();
2892  for (int i = 0; i <nConstituents ; i++){
2893  UsedTowerList.push_back(jetCaloRefs[i]);
2894  }
2895 
2896  SumPtJet +=jetPt;
2897  // }
2898 
2899  }
2900 
2901 
2902  // Handle<CaloTowerCollection> caloTowers;
2903  // evt.getByLabel( "towerMaker", caloTowers );
2904 
2905  NTowersUsed = UsedTowerList.size();
2906 
2907  // --- Loop over towers and make a lists of used and unused towers
2908  for (CaloTowerCollection::const_iterator tower = caloTowers->begin();
2909  tower != caloTowers->end(); tower++) {
2910 
2911  CaloTower t = *tower;
2912  Double_t et = tower->et();
2913 
2914  if(et>0) {
2915 
2916  Double_t phi = tower->phi();
2917 
2918  SumEtTowers += tower->et();
2919 
2920  sumTowerAllPx += et*cos(phi);
2921  sumTowerAllPy += et*sin(phi);
2922 
2923  bool used = false;
2924 
2925  for(int i=0; i<NTowersUsed; i++){
2926  if(tower->id() == UsedTowerList[i]->id()){
2927  used=true;
2928  break;
2929  }
2930  }
2931 
2932  if (used) {
2933  TowerUsedInJets.push_back(t);
2934  } else {
2935  TowerNotUsedInJets.push_back(t);
2936  }
2937 
2938  }
2939 
2940  }
2941 
2942  nUsed = TowerUsedInJets.size();
2943  nNotUsed = TowerNotUsedInJets.size();
2944 
2945  SumEtJets = 0;
2946  SumEtNotJets = 0;
2947 
2948  for(int i=0;i<nUsed;i++){
2949  SumEtJets += TowerUsedInJets[i].et();
2950  }
2951  h_jetEt2.Fill(SumEtJets);
2952 
2953  for(int i=0;i<nNotUsed;i++){
2954  if (TowerNotUsedInJets[i].et() > 0.5)
2955  SumEtNotJets += TowerNotUsedInJets[i].et();
2956  h_missEt2.Fill(TowerNotUsedInJets[i].et());
2957  h_missEt2s.Fill(TowerNotUsedInJets[i].et());
2958  }
2959  h_totMissEt2.Fill(SumEtNotJets);
2960 
2961 
2962  // *********************
2963  // KtClus
2964  //
2965  UsedTowerList.clear();
2966  TowerUsedInJets.clear();
2967  TowerNotUsedInJets.clear();
2968 
2969  // --- Loop over jets and make a list of all the used towers
2970  evt.getByLabel( CaloJetAlgorithm3, jets );
2971  for ( CaloJetCollection::const_iterator ijet=jets->begin(); ijet!=jets->end(); ijet++) {
2972 
2973  Double_t jetPt = ijet->pt();
2974  Double_t jetPhi = ijet->phi();
2975 
2976  // if (jetPt>5.0) {
2977 
2978  Double_t jetPx = jetPt*cos(jetPhi);
2979  Double_t jetPy = jetPt*sin(jetPhi);
2980 
2981  sumJetPx +=jetPx;
2982  sumJetPy +=jetPy;
2983 
2984  const std::vector<CaloTowerPtr> jetCaloRefs = ijet->getCaloConstituents();
2985  int nConstituents = jetCaloRefs.size();
2986  for (int i = 0; i <nConstituents ; i++){
2987  UsedTowerList.push_back(jetCaloRefs[i]);
2988  }
2989 
2990  SumPtJet +=jetPt;
2991  // }
2992 
2993  }
2994 
2995 
2996  // Handle<CaloTowerCollection> caloTowers;
2997  // evt.getByLabel( "towerMaker", caloTowers );
2998 
2999  NTowersUsed = UsedTowerList.size();
3000 
3001  // --- Loop over towers and make a lists of used and unused towers
3002  for (CaloTowerCollection::const_iterator tower = caloTowers->begin();
3003  tower != caloTowers->end(); tower++) {
3004 
3005  CaloTower t = *tower;
3006  Double_t et = tower->et();
3007 
3008  if(et>0) {
3009 
3010  // Double_t phi = tower->phi();
3011 
3012  // SumEtTowers += tower->et();
3013  // sumTowerAllPx += et*cos(phi);
3014  // sumTowerAllPy += et*sin(phi);
3015 
3016  bool used = false;
3017 
3018  for(int i=0; i<NTowersUsed; i++){
3019  if(tower->id() == UsedTowerList[i]->id()){
3020  used=true;
3021  break;
3022  }
3023  }
3024 
3025  if (used) {
3026  TowerUsedInJets.push_back(t);
3027  } else {
3028  TowerNotUsedInJets.push_back(t);
3029  }
3030 
3031  }
3032 
3033  }
3034 
3035  nUsed = TowerUsedInJets.size();
3036  nNotUsed = TowerNotUsedInJets.size();
3037 
3038  SumEtJets = 0;
3039  SumEtNotJets = 0;
3040 
3041  for(int i=0;i<nUsed;i++){
3042  SumEtJets += TowerUsedInJets[i].et();
3043  }
3044  h_jetEt3.Fill(SumEtJets);
3045 
3046  for(int i=0;i<nNotUsed;i++){
3047  if (TowerNotUsedInJets[i].et() > 0.5)
3048  SumEtNotJets += TowerNotUsedInJets[i].et();
3049  h_missEt3.Fill(TowerNotUsedInJets[i].et());
3050  h_missEt3s.Fill(TowerNotUsedInJets[i].et());
3051  }
3052  h_totMissEt3.Fill(SumEtNotJets);
3053 
3054 }
myFastSimVal::ZpMassResRL132
TH1F ZpMassResRL132
Definition: myFastSimVal.h:91
myFastSimVal::totEneLeadJetEta1_1
TH1F totEneLeadJetEta1_1
Definition: myFastSimVal.h:45
myFastSimVal::ZpMassResL403
TH1F ZpMassResL403
Definition: myFastSimVal.h:89
myFastSimVal::ZpMassResRLoP102
TH1F ZpMassResRLoP102
Definition: myFastSimVal.h:93
myFastSimVal::ZpMassGen_2000_40
TH1F ZpMassGen_2000_40
Definition: myFastSimVal.h:80
custom_jme_cff.nConstituents
nConstituents
Definition: custom_jme_cff.py:137
myFastSimVal::ZpMass_5000_40
TH1F ZpMass_5000_40
Definition: myFastSimVal.h:81
myFastSimVal::ZpMassGen_5000_13
TH1F ZpMassGen_5000_13
Definition: myFastSimVal.h:82
myFastSimVal::hf_nJet2
TProfile hf_nJet2
Definition: myFastSimVal.h:32
myFastSimVal::hadEneLeadJetEta1_2
TH1F hadEneLeadJetEta1_2
Definition: myFastSimVal.h:49
myFastSimVal::ZpMassResL133
TH1F ZpMassResL133
Definition: myFastSimVal.h:88
myFastSimVal::hadEneLeadJetEta3_2
TH1F hadEneLeadJetEta3_2
Definition: myFastSimVal.h:49
myFastSimVal::dijetMass_700_101
TH1F dijetMass_700_101
Definition: myFastSimVal.h:114
myFastSimVal::matchedAllPt31
TH1F matchedAllPt31
Definition: myFastSimVal.h:54
myFastSimVal::hadFracEta31
TH1F hadFracEta31
Definition: myFastSimVal.h:67
myFastSimVal::GenJetAlgorithm3
std::string GenJetAlgorithm3
Definition: myFastSimVal.h:26
myFastSimVal::hadFracEta11
TH1F hadFracEta11
Definition: myFastSimVal.h:67
Scenarios_cff.scale
scale
Definition: Scenarios_cff.py:2186
mps_fire.i
i
Definition: mps_fire.py:338
myFastSimVal::dijetMassCor_700_101
TH1F dijetMassCor_700_101
Definition: myFastSimVal.h:103
myFastSimVal::ZpMassResPRL402
TH1F ZpMassResPRL402
Definition: myFastSimVal.h:98
myFastSimVal::dijetMass_5000_131
TH1F dijetMass_5000_131
Definition: myFastSimVal.h:116
globals_cff.id1
id1
Definition: globals_cff.py:33
myFastSimVal::ZpMassResPRL131
TH1F ZpMassResPRL131
Definition: myFastSimVal.h:97
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
myFastSimVal::matchedAllPt12
TH1F matchedAllPt12
Definition: myFastSimVal.h:52
reco::Candidate::vx
virtual double vx() const =0
x coordinate of vertex position
myFastSimVal::h_lowPtCal1c12
TH1F h_lowPtCal1c12
Definition: myFastSimVal.h:178
myFastSimVal::ZpMassResRLoP131
TH1F ZpMassResRLoP131
Definition: myFastSimVal.h:94
myFastSimVal::ZpMassResPRL103
TH1F ZpMassResPRL103
Definition: myFastSimVal.h:96
myFastSimVal::CaloJetAlgorithm3
std::string CaloJetAlgorithm3
Definition: myFastSimVal.h:24
myFastSimVal::ZpMass_700_13
TH1F ZpMass_700_13
Definition: myFastSimVal.h:77
myFastSimVal::GenJetAlgorithm1
std::string GenJetAlgorithm1
Definition: myFastSimVal.h:26
myFastSimVal::hf_TowerDelR1
TProfile hf_TowerDelR1
Definition: myFastSimVal.h:30
myFastSimVal::tbarMassGen
TH1F tbarMassGen
Definition: myFastSimVal.h:86
myFastSimVal::h_lowPtCal2c12
TH1F h_lowPtCal2c12
Definition: myFastSimVal.h:178
myFastSimVal::ZpMassResL101
TH1F ZpMassResL101
Definition: myFastSimVal.h:87
myFastSimVal::totEneLeadJetEta3_2
TH1F totEneLeadJetEta3_2
Definition: myFastSimVal.h:48
myFastSimVal::ZpMassResRLoP402
TH1F ZpMassResRLoP402
Definition: myFastSimVal.h:95
myFastSimVal::matchedAllPt43
TH1F matchedAllPt43
Definition: myFastSimVal.h:55
myFastSimVal::CaloJetAlgorithm4
std::string CaloJetAlgorithm4
Definition: myFastSimVal.h:24
reco::Candidate::daughter
virtual const Candidate * daughter(size_type i) const =0
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
myFastSimVal::dijetMassP_5000_101
TH1F dijetMassP_5000_101
Definition: myFastSimVal.h:122
isolatedTracks_cfi.caloJets
caloJets
Definition: isolatedTracks_cfi.py:33
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
myFastSimVal::hadEneLeadJetEta3_1
TH1F hadEneLeadJetEta3_1
Definition: myFastSimVal.h:46
edm::OwnVector::size
size_type size() const
Definition: OwnVector.h:264
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:82
myFastSimVal::hf_nJet31
TProfile hf_nJet31
Definition: myFastSimVal.h:34
myFastSimVal::hadEneLeadJet12
TH1F hadEneLeadJet12
Definition: myFastSimVal.h:62
myFastSimVal::ZpMassGen10
TH1F ZpMassGen10
Definition: myFastSimVal.h:76
myFastSimVal::dijetMassCor_5000_101
TH1F dijetMassCor_5000_101
Definition: myFastSimVal.h:105
myFastSimVal::h_lowPtCal1c11
TH1F h_lowPtCal1c11
Definition: myFastSimVal.h:177
reco::Candidate::mother
virtual const Candidate * mother(size_type i=0) const =0
return pointer to mother
myFastSimVal::h_lowPtCal3c12
TH1F h_lowPtCal3c12
Definition: myFastSimVal.h:178
eleHcalExtractorBlocks_cff.caloTowers
caloTowers
Definition: eleHcalExtractorBlocks_cff.py:15
myFastSimVal::dijetMassP101
TH1F dijetMassP101
Definition: myFastSimVal.h:118
myFastSimVal::dijetMassCor_2000_131
TH1F dijetMassCor_2000_131
Definition: myFastSimVal.h:104
myFastSimVal::ZpMassResL401
TH1F ZpMassResL401
Definition: myFastSimVal.h:89
myFastSimVal::dijetMass_700_131
TH1F dijetMass_700_131
Definition: myFastSimVal.h:114
reco::Candidate::vy
virtual double vy() const =0
y coordinate of vertex position
myFastSimVal::dijetMassCor_5000_131
TH1F dijetMassCor_5000_131
Definition: myFastSimVal.h:105
myFastSimVal::emEneLeadJetEta3_1
TH1F emEneLeadJetEta3_1
Definition: myFastSimVal.h:47
myFastSimVal::ZpMass_700_10
TH1F ZpMass_700_10
Definition: myFastSimVal.h:77
myFastSimVal::h_lowPtCal4c13
TH1F h_lowPtCal4c13
Definition: myFastSimVal.h:179
myFastSimVal::nTowersSecondJet1
TH1F nTowersSecondJet1
Definition: myFastSimVal.h:70
myFastSimVal::ZpMassResRL103
TH1F ZpMassResRL103
Definition: myFastSimVal.h:90
njet
int njet
Definition: HydjetWrapper.h:95
myFastSimVal::ZpMassGen_2000_10
TH1F ZpMassGen_2000_10
Definition: myFastSimVal.h:80
myFastSimVal::dijetMass_2000_131
TH1F dijetMass_2000_131
Definition: myFastSimVal.h:115
myFastSimVal::hadEneLeadJet13
TH1F hadEneLeadJet13
Definition: myFastSimVal.h:63
myFastSimVal::dijetMassP_2000_101
TH1F dijetMassP_2000_101
Definition: myFastSimVal.h:121
myFastSimVal::ZpMassMatched1
TH1F ZpMassMatched1
Definition: myFastSimVal.h:107
nEvent
int nEvent
Definition: myFastSimVal.cc:49
myFastSimVal::hf_nJet3
TProfile hf_nJet3
Definition: myFastSimVal.h:32
MAXJETS
#define MAXJETS
Definition: myFastSimVal.cc:30
myFastSimVal::dijetMass_700_401
TH1F dijetMass_700_401
Definition: myFastSimVal.h:114
myFastSimVal::nTowersLeadJet1
TH1F nTowersLeadJet1
Definition: myFastSimVal.h:70
myFastSimVal::matchedAllPt32
TH1F matchedAllPt32
Definition: myFastSimVal.h:54
myFastSimVal::ZpMassResRLoP403
TH1F ZpMassResRLoP403
Definition: myFastSimVal.h:95
myFastSimVal::ZpMassGen_5000_10
TH1F ZpMassGen_5000_10
Definition: myFastSimVal.h:82
myFastSimVal::matchedAllPt23
TH1F matchedAllPt23
Definition: myFastSimVal.h:53
myFastSimVal::ZpMassResPRL101
TH1F ZpMassResPRL101
Definition: myFastSimVal.h:96
myFastSimVal::topMassParton
TH1F topMassParton
Definition: myFastSimVal.h:83
myFastSimVal::emEneLeadJetEta2_1
TH1F emEneLeadJetEta2_1
Definition: myFastSimVal.h:47
myFastSimVal::dijetMassCor_2000_101
TH1F dijetMassCor_2000_101
Definition: myFastSimVal.h:104
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
myFastSimVal::ZpMassResPRL401
TH1F ZpMassResPRL401
Definition: myFastSimVal.h:98
myFastSimVal::matchedPt23
TH1F matchedPt23
Definition: myFastSimVal.h:57
myFastSimVal::matchedPt33
TH1F matchedPt33
Definition: myFastSimVal.h:58
myFastSimVal::emEneLeadJet13
TH1F emEneLeadJet13
Definition: myFastSimVal.h:66
myFastSimVal::hf_PtResponse1
TProfile hf_PtResponse1
Definition: myFastSimVal.h:35
myFastSimVal::matchedPt32
TH1F matchedPt32
Definition: myFastSimVal.h:58
myFastSimVal::ZpMassResL103
TH1F ZpMassResL103
Definition: myFastSimVal.h:87
myFastSimVal::matchedAllPt21
TH1F matchedAllPt21
Definition: myFastSimVal.h:53
myFastSimVal::dijetMassCor_2000_401
TH1F dijetMassCor_2000_401
Definition: myFastSimVal.h:104
myFastSimVal::dijetMassP_5000_401
TH1F dijetMassP_5000_401
Definition: myFastSimVal.h:122
myFastSimVal::hf_nJet4s
TProfile hf_nJet4s
Definition: myFastSimVal.h:33
myFastSimVal::nTowersLeadJetPt1
TH1F nTowersLeadJetPt1
Definition: myFastSimVal.h:44
myFastSimVal::ZpMassResRL401
TH1F ZpMassResRL401
Definition: myFastSimVal.h:92
myFastSimVal::ZpMass_5000_13
TH1F ZpMass_5000_13
Definition: myFastSimVal.h:81
myFastSimVal::ZpMassResL402
TH1F ZpMassResL402
Definition: myFastSimVal.h:89
myFastSimVal::nTowersLeadJetPt3
TH1F nTowersLeadJetPt3
Definition: myFastSimVal.h:44
myFastSimVal::hf_nJet1
TProfile hf_nJet1
Definition: myFastSimVal.h:32
myFastSimVal::ZpMassMatched3
TH1F ZpMassMatched3
Definition: myFastSimVal.h:107
myFastSimVal::ZpMassResRL131
TH1F ZpMassResRL131
Definition: myFastSimVal.h:91
myFastSimVal::dijetMassP_5000_131
TH1F dijetMassP_5000_131
Definition: myFastSimVal.h:122
myFastSimVal::matchedPt21
TH1F matchedPt21
Definition: myFastSimVal.h:57
myFastSimVal::h_lowPtCal4c12
TH1F h_lowPtCal4c12
Definition: myFastSimVal.h:178
myFastSimVal::hf_nJet1s
TProfile hf_nJet1s
Definition: myFastSimVal.h:33
myFastSimVal::hadFracEta21
TH1F hadFracEta21
Definition: myFastSimVal.h:67
myFastSimVal::hf_nJet3s
TProfile hf_nJet3s
Definition: myFastSimVal.h:33
myFastSimVal::ZpMassResRLoP133
TH1F ZpMassResRLoP133
Definition: myFastSimVal.h:94
myFastSimVal::matchedPt41
TH1F matchedPt41
Definition: myFastSimVal.h:59
boostedElectronIsolation_cff.deltaR
deltaR
Definition: boostedElectronIsolation_cff.py:77
reco::Candidate::p4
virtual const LorentzVector & p4() const =0
four-momentum Lorentz vector
myFastSimVal::nTowersLeadJetPt4
TH1F nTowersLeadJetPt4
Definition: myFastSimVal.h:44
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
protons_cff.t
t
Definition: protons_cff.py:34
myFastSimVal::ZpMassResRL102
TH1F ZpMassResRL102
Definition: myFastSimVal.h:90
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:42
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
myFastSimVal::ZpMassResPRL102
TH1F ZpMassResPRL102
Definition: myFastSimVal.h:96
myFastSimVal::dijetMassP_700_101
TH1F dijetMassP_700_101
Definition: myFastSimVal.h:120
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
myFastSimVal::hf_PtResponse2
TProfile hf_PtResponse2
Definition: myFastSimVal.h:35
myFastSimVal::dijetMassCor_5000_1
TH1F dijetMassCor_5000_1
Definition: myFastSimVal.h:105
myFastSimVal::dijetMass_5000_101
TH1F dijetMass_5000_101
Definition: myFastSimVal.h:116
myFastSimVal::dijetMassCor_700_131
TH1F dijetMassCor_700_131
Definition: myFastSimVal.h:103
myFastSimVal::hf_PtResponse3
TProfile hf_PtResponse3
Definition: myFastSimVal.h:35
myFastSimVal::totEneLeadJetEta3_1
TH1F totEneLeadJetEta3_1
Definition: myFastSimVal.h:45
myFastSimVal::hf_TowerDelR12
TProfile hf_TowerDelR12
Definition: myFastSimVal.h:31
myFastSimVal::matchedAllPt41
TH1F matchedAllPt41
Definition: myFastSimVal.h:55
myFastSimVal::ZpMassResPRL133
TH1F ZpMassResPRL133
Definition: myFastSimVal.h:97
myFastSimVal::ZpMassGen_2000_13
TH1F ZpMassGen_2000_13
Definition: myFastSimVal.h:80
GenHFHadronMatcher_cfi.genParticles
genParticles
Definition: GenHFHadronMatcher_cfi.py:4
myFastSimVal::ZpMass_2000_40
TH1F ZpMass_2000_40
Definition: myFastSimVal.h:79
myFastSimVal::ZpMassResL131
TH1F ZpMassResL131
Definition: myFastSimVal.h:88
relval_steps.gen
def gen(fragment, howMuch)
Production test section ####.
Definition: relval_steps.py:490
myFastSimVal::ZpMassResRLoP401
TH1F ZpMassResRLoP401
Definition: myFastSimVal.h:95
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:480
myFastSimVal::ZpMassRes101
TH1F ZpMassRes101
Definition: myFastSimVal.h:99
myFastSimVal::ZpMass_700_40
TH1F ZpMass_700_40
Definition: myFastSimVal.h:77
myFastSimVal::hadEneLeadJetEta2_2
TH1F hadEneLeadJetEta2_2
Definition: myFastSimVal.h:49
myFastSimVal::ZpMassResL132
TH1F ZpMassResL132
Definition: myFastSimVal.h:88
myFastSimVal::ZpMassGen_700_10
TH1F ZpMassGen_700_10
Definition: myFastSimVal.h:78
myFastSimVal::ZpMassResRLoP103
TH1F ZpMassResRLoP103
Definition: myFastSimVal.h:93
cuy.ii
ii
Definition: cuy.py:590
myFastSimVal::ZpMassResL102
TH1F ZpMassResL102
Definition: myFastSimVal.h:87
myFastSimVal::ZpMassGen_5000_40
TH1F ZpMassGen_5000_40
Definition: myFastSimVal.h:82
myFastSimVal::TowerEtLeadJet12
TH1F TowerEtLeadJet12
Definition: myFastSimVal.h:74
myFastSimVal::matchedAllPt13
TH1F matchedAllPt13
Definition: myFastSimVal.h:52
myFastSimVal::h_lowPtCal3c11
TH1F h_lowPtCal3c11
Definition: myFastSimVal.h:177
myFastSimVal::hadEneLeadJetEta1_1
TH1F hadEneLeadJetEta1_1
Definition: myFastSimVal.h:46
myFastSimVal::ZpMassGen_700_13
TH1F ZpMassGen_700_13
Definition: myFastSimVal.h:78
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
myFastSimVal::ZpMass_2000_13
TH1F ZpMass_2000_13
Definition: myFastSimVal.h:79
myFastSimVal::matchedAllPt22
TH1F matchedAllPt22
Definition: myFastSimVal.h:53
myFastSimVal::dijetMassCor_700_1
TH1F dijetMassCor_700_1
Definition: myFastSimVal.h:103
myFastSimVal::dijetMassCor_2000_1
TH1F dijetMassCor_2000_1
Definition: myFastSimVal.h:104
myFastSimVal::hf_nJet41
TProfile hf_nJet41
Definition: myFastSimVal.h:34
myFastSimVal::totEneLeadJetEta2_1
TH1F totEneLeadJetEta2_1
Definition: myFastSimVal.h:45
myFastSimVal::matchedPt13
TH1F matchedPt13
Definition: myFastSimVal.h:56
myFastSimVal::dijetMassCor_5000_401
TH1F dijetMassCor_5000_401
Definition: myFastSimVal.h:105
myFastSimVal::GenJetAlgorithm4
std::string GenJetAlgorithm4
Definition: myFastSimVal.h:26
nPart
TString nPart(Int_t part, TString string, TString delimit=";", Bool_t removerest=true)
myFastSimVal::dijetMassCor_700_401
TH1F dijetMassCor_700_401
Definition: myFastSimVal.h:103
myFastSimVal::ZpMassMatched2
TH1F ZpMassMatched2
Definition: myFastSimVal.h:107
myFastSimVal::matchedPt12
TH1F matchedPt12
Definition: myFastSimVal.h:56
myFastSimVal::ZpMass_2000_10
TH1F ZpMass_2000_10
Definition: myFastSimVal.h:79
myFastSimVal::matchedPt22
TH1F matchedPt22
Definition: myFastSimVal.h:57
myFastSimVal::dijetMassP401
TH1F dijetMassP401
Definition: myFastSimVal.h:118
myFastSimVal::hadFracLeadJet1
TH1F hadFracLeadJet1
Definition: myFastSimVal.h:70
common_cff.pdgId
pdgId
Definition: common_cff.py:54
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
myFastSimVal::ZpMassGen13
TH1F ZpMassGen13
Definition: myFastSimVal.h:76
reco::Candidate::mass
virtual double mass() const =0
mass
myFastSimVal::dijetMassCor101
TH1F dijetMassCor101
Definition: myFastSimVal.h:124
myFastSimVal::h_lowPtCal4c11
TH1F h_lowPtCal4c11
Definition: myFastSimVal.h:177
myFastSimVal::dijetMassCor131
TH1F dijetMassCor131
Definition: myFastSimVal.h:124
myFastSimVal::hf_sumTowerAllEy
TH1F hf_sumTowerAllEy
Definition: myFastSimVal.h:37
myFastSimVal::emEneLeadJetEta2_2
TH1F emEneLeadJetEta2_2
Definition: myFastSimVal.h:50
myFastSimVal::emEneLeadJetEta1_1
TH1F emEneLeadJetEta1_1
Definition: myFastSimVal.h:47
myFastSimVal::dijetMassP_700_131
TH1F dijetMassP_700_131
Definition: myFastSimVal.h:120
globals_cff.id2
id2
Definition: globals_cff.py:34
myFastSimVal::emEneLeadJet12
TH1F emEneLeadJet12
Definition: myFastSimVal.h:65
myFastSimVal::h_lowPtCal2c13
TH1F h_lowPtCal2c13
Definition: myFastSimVal.h:179
stringResolutionProvider_cfi.et
et
define resolution functions of each parameter
Definition: stringResolutionProvider_cfi.py:13
myFastSimVal::hf_nJet21
TProfile hf_nJet21
Definition: myFastSimVal.h:34
myFastSimVal::dijetMass_2000_401
TH1F dijetMass_2000_401
Definition: myFastSimVal.h:115
myFastSimVal::emEneLeadJetEta3_2
TH1F emEneLeadJetEta3_2
Definition: myFastSimVal.h:50
myFastSimVal::emEneLeadJetEta1_2
TH1F emEneLeadJetEta1_2
Definition: myFastSimVal.h:50
myFastSimVal::dijetMass_2000_101
TH1F dijetMass_2000_101
Definition: myFastSimVal.h:115
myFastSimVal::ZpMassResRL133
TH1F ZpMassResRL133
Definition: myFastSimVal.h:91
myFastSimVal::ZpMassResPRL132
TH1F ZpMassResPRL132
Definition: myFastSimVal.h:97
myFastSimVal::CaloJetAlgorithm1
std::string CaloJetAlgorithm1
Definition: myFastSimVal.h:24
myFastSimVal::totEneLeadJetEta2_2
TH1F totEneLeadJetEta2_2
Definition: myFastSimVal.h:48
myFastSimVal::ZpMassGen40
TH1F ZpMassGen40
Definition: myFastSimVal.h:76
myFastSimVal::emEneLeadJet1
TH1F emEneLeadJet1
Definition: myFastSimVal.h:64
myFastSimVal::hf_sumTowerAllEx
TH1F hf_sumTowerAllEx
Definition: myFastSimVal.h:37
myFastSimVal::hf_TowerJetEt1
TH1F hf_TowerJetEt1
Definition: myFastSimVal.h:41
myFastSimVal::hf_nJet2s
TProfile hf_nJet2s
Definition: myFastSimVal.h:33
myFastSimVal::hadEneLeadJetEta2_1
TH1F hadEneLeadJetEta2_1
Definition: myFastSimVal.h:46
myFastSimVal::ZpMassResRL402
TH1F ZpMassResRL402
Definition: myFastSimVal.h:92
myFastSimVal::matchedAllPt11
TH1F matchedAllPt11
Definition: myFastSimVal.h:52
myFastSimVal::ZpMass_5000_10
TH1F ZpMass_5000_10
Definition: myFastSimVal.h:81
hiPixelPairStep_cff.deltaPhi
deltaPhi
Definition: hiPixelPairStep_cff.py:87
myFastSimVal::matchedAllPt33
TH1F matchedAllPt33
Definition: myFastSimVal.h:54
myFastSimVal::matchedPt42
TH1F matchedPt42
Definition: myFastSimVal.h:59
myFastSimVal::ZpMassGen_700_40
TH1F ZpMassGen_700_40
Definition: myFastSimVal.h:78
CaloTower::et
double et(double vtxZ) const
Definition: CaloTower.h:132
myFastSimVal::hf_nJet11
TProfile hf_nJet11
Definition: myFastSimVal.h:34
gather_cfg.cout
cout
Definition: gather_cfg.py:144
myFastSimVal::nTowersLeadJetPt2
TH1F nTowersLeadJetPt2
Definition: myFastSimVal.h:44
myFastSimVal::TowerEtLeadJet1
TH1F TowerEtLeadJet1
Definition: myFastSimVal.h:73
reco::Candidate::vz
virtual double vz() const =0
z coordinate of vertex position
myFastSimVal::h_lowPtCal1c13
TH1F h_lowPtCal1c13
Definition: myFastSimVal.h:179
myFastSimVal::ZpMassResRL101
TH1F ZpMassResRL101
Definition: myFastSimVal.h:90
myFastSimVal::ZpMassResRL403
TH1F ZpMassResRL403
Definition: myFastSimVal.h:92
myFastSimVal::GenJetAlgorithm2
std::string GenJetAlgorithm2
Definition: myFastSimVal.h:26
myFastSimVal::hf_PtResponse4
TProfile hf_PtResponse4
Definition: myFastSimVal.h:35
myFastSimVal::dijetMassCor1
TH1F dijetMassCor1
Definition: myFastSimVal.h:124
myFastSimVal::dijetMassCor401
TH1F dijetMassCor401
Definition: myFastSimVal.h:124
myFastSimVal::hf_nJet4
TProfile hf_nJet4
Definition: myFastSimVal.h:32
myFastSimVal::CaloJetAlgorithm2
std::string CaloJetAlgorithm2
Definition: myFastSimVal.h:24
myFastSimVal::dijetMassP131
TH1F dijetMassP131
Definition: myFastSimVal.h:118
myFastSimVal::h_lowPtCal2c11
TH1F h_lowPtCal2c11
Definition: myFastSimVal.h:177
reco::Candidate::numberOfDaughters
virtual size_type numberOfDaughters() const =0
number of daughters
myFastSimVal::totEneLeadJetEta1_2
TH1F totEneLeadJetEta1_2
Definition: myFastSimVal.h:48
myFastSimVal::ZpMassResPRL403
TH1F ZpMassResPRL403
Definition: myFastSimVal.h:98
myFastSimVal::TowerEtLeadJet13
TH1F TowerEtLeadJet13
Definition: myFastSimVal.h:75
myFastSimVal::matchedPt43
TH1F matchedPt43
Definition: myFastSimVal.h:59
myFastSimVal::dijetMassP_2000_131
TH1F dijetMassP_2000_131
Definition: myFastSimVal.h:121
myFastSimVal::matchedPt31
TH1F matchedPt31
Definition: myFastSimVal.h:58
myFastSimVal::matchedAllPt42
TH1F matchedAllPt42
Definition: myFastSimVal.h:55
ALCARECOTkAlWMuNu_cff.minJetPt
minJetPt
Definition: ALCARECOTkAlWMuNu_cff.py:47
myFastSimVal::hadEneLeadJet1
TH1F hadEneLeadJet1
Definition: myFastSimVal.h:61
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
myFastSimVal::ZpMassResRLoP132
TH1F ZpMassResRLoP132
Definition: myFastSimVal.h:94
myFastSimVal::matchedPt11
TH1F matchedPt11
Definition: myFastSimVal.h:56
myFastSimVal::dijetMass_5000_401
TH1F dijetMass_5000_401
Definition: myFastSimVal.h:116
myFastSimVal::h_lowPtCal3c13
TH1F h_lowPtCal3c13
Definition: myFastSimVal.h:179
myFastSimVal::ZpMassResRLoP101
TH1F ZpMassResRLoP101
Definition: myFastSimVal.h:93
ttbarCategorization_cff.genJets
genJets
Definition: ttbarCategorization_cff.py:29
void myFastSimVal::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file myFastSimVal.cc.

References dEta1, dEta12, dEta2, dEta3, dEta4, dEtaPar1, dEtaPar2, dEtaPar3, dEtaPar4, dijetMass1, dijetMass101, dijetMass102, dijetMass103, dijetMass12, dijetMass13, dijetMass131, dijetMass132, dijetMass133, dijetMass2, dijetMass22, dijetMass23, dijetMass3, dijetMass32, dijetMass33, dijetMass4, dijetMass401, dijetMass402, dijetMass403, dijetMass42, dijetMass43, dijetMass_2000_101, dijetMass_2000_131, dijetMass_2000_401, dijetMass_5000_101, dijetMass_5000_131, dijetMass_5000_401, dijetMass_700_101, dijetMass_700_131, dijetMass_700_401, dijetMassCor1, dijetMassCor101, dijetMassCor131, dijetMassCor401, dijetMassCor_2000_1, dijetMassCor_2000_101, dijetMassCor_2000_131, dijetMassCor_2000_401, dijetMassCor_5000_1, dijetMassCor_5000_101, dijetMassCor_5000_131, dijetMassCor_5000_401, dijetMassCor_700_1, dijetMassCor_700_101, dijetMassCor_700_131, dijetMassCor_700_401, dijetMassP1, dijetMassP101, dijetMassP131, dijetMassP2, dijetMassP3, dijetMassP401, dijetMassP_2000_101, dijetMassP_2000_131, dijetMassP_2000_401, dijetMassP_5000_101, dijetMassP_5000_131, dijetMassP_5000_401, dijetMassP_700_101, dijetMassP_700_131, dijetMassP_700_401, dPhi1, dPhi12, dPhi2, dPhi3, dPhi4, dPhiPar1, dPhiPar2, dPhiPar3, dPhiPar4, dPt1, dPt100Frac1, dPt100Frac2, dPt100Frac3, dPt100Frac4, dPt12, dPt2, dPt20Frac1, dPt20Frac2, dPt20Frac3, dPt20Frac4, dPt3, dPt4, dPt40Frac1, dPt40Frac2, dPt40Frac3, dPt40Frac4, dPt80Frac1, dPt80Frac2, dPt80Frac3, dPt80Frac4, dPtFrac1, dPtFrac2, dPtFrac3, dPtFrac4, dPtPar1, dPtPar2, dPtPar3, dPtPar4, dR1, dR12, dR2, dR3, dR4, dRPar1, dRPar2, dRPar3, dRPar4, dRParton, dRPartonMin, emEneLeadJet1, emEneLeadJet12, emEneLeadJet13, emEneLeadJet2, emEneLeadJet22, emEneLeadJet23, emEneLeadJet3, emEneLeadJet32, emEneLeadJet33, emEneLeadJetEta1_1, emEneLeadJetEta1_2, emEneLeadJetEta2_1, emEneLeadJetEta2_2, emEneLeadJetEta3_1, emEneLeadJetEta3_2, h_etaCal1, h_etaCal2, h_etaCal3, h_etaCal4, h_etaCalL1, h_etaCalL2, h_etaCalL3, h_etaCalL4, h_etaGen1, h_etaGen2, h_etaGen3, h_etaGenL1, h_etaGenL2, h_etaGenL3, h_jet1Pt1, h_jet1Pt2, h_jet1Pt3, h_jet1Pt4, h_jet2Pt1, h_jet2Pt2, h_jet2Pt3, h_jet2Pt4, h_jet3Pt1, h_jet3Pt2, h_jet3Pt3, h_jet3Pt4, h_jet4Pt1, h_jet4Pt2, h_jet4Pt3, h_jet4Pt4, h_jet5Pt1, h_jet5Pt2, h_jet5Pt3, h_jet5Pt4, h_jet6Pt1, h_jet6Pt2, h_jet6Pt3, h_jet6Pt4, h_jet7Pt1, h_jet7Pt2, h_jet7Pt3, h_jet7Pt4, h_jetEt1, h_jetEt2, h_jetEt3, h_lowPtCal11, h_lowPtCal12, h_lowPtCal13, h_lowPtCal1c11, h_lowPtCal1c12, h_lowPtCal1c13, h_lowPtCal21, h_lowPtCal22, h_lowPtCal23, h_lowPtCal2c11, h_lowPtCal2c12, h_lowPtCal2c13, h_lowPtCal31, h_lowPtCal32, h_lowPtCal33, h_lowPtCal3c11, h_lowPtCal3c12, h_lowPtCal3c13, h_lowPtCal41, h_lowPtCal42, h_lowPtCal43, h_lowPtCal4c11, h_lowPtCal4c12, h_lowPtCal4c13, h_missEt1, h_missEt1s, h_missEt2, h_missEt2s, h_missEt3, h_missEt3s, h_nCalJets1, h_nCalJets2, h_nCalJets3, h_nCalJets4, h_nGenJets1, h_nGenJets2, h_nGenJets3, h_nGenJets4, h_phiCal1, h_phiCal2, h_phiCal3, h_phiCal4, h_phiCalL1, h_phiCalL2, h_phiCalL3, h_phiCalL4, h_phiGen1, h_phiGen2, h_phiGen3, h_phiGenL1, h_phiGenL2, h_phiGenL3, h_ptCal1, h_ptCal12, h_ptCal13, h_ptCal2, h_ptCal22, h_ptCal23, h_ptCal3, h_ptCal32, h_ptCal33, h_ptCal4, h_ptCal42, h_ptCal43, h_ptCalL1, h_ptCalL12, h_ptCalL13, h_ptCalL2, h_ptCalL22, h_ptCalL23, h_ptCalL3, h_ptCalL32, h_ptCalL33, h_ptCalL4, h_ptCalL42, h_ptCalL43, h_ptGen1, h_ptGen12, h_ptGen13, h_ptGen2, h_ptGen22, h_ptGen23, h_ptGen3, h_ptGen32, h_ptGen33, h_ptGen4, h_ptGen42, h_ptGen43, h_ptGenL1, h_ptGenL12, h_ptGenL13, h_ptGenL2, h_ptGenL22, h_ptGenL23, h_ptGenL3, h_ptGenL32, h_ptGenL33, h_totMissEt1, h_totMissEt2, h_totMissEt3, hadEneLeadJet1, hadEneLeadJet12, hadEneLeadJet13, hadEneLeadJet2, hadEneLeadJet22, hadEneLeadJet23, hadEneLeadJet3, hadEneLeadJet32, hadEneLeadJet33, hadEneLeadJetEta1_1, hadEneLeadJetEta1_2, hadEneLeadJetEta2_1, hadEneLeadJetEta2_2, hadEneLeadJetEta3_1, hadEneLeadJetEta3_2, hadFracEta11, hadFracEta12, hadFracEta13, hadFracEta21, hadFracEta22, hadFracEta23, hadFracEta31, hadFracEta32, hadFracEta33, hadFracLeadJet1, hadFracLeadJet2, hadFracLeadJet3, hf_nJet1, hf_nJet11, hf_nJet1s, hf_nJet2, hf_nJet21, hf_nJet2s, hf_nJet3, hf_nJet31, hf_nJet3s, hf_nJet4, hf_nJet41, hf_nJet4s, hf_PtResponse1, hf_PtResponse2, hf_PtResponse3, hf_PtResponse4, hf_sumTowerAllEx, hf_sumTowerAllEy, hf_TowerDelR1, hf_TowerDelR12, hf_TowerDelR2, hf_TowerDelR22, hf_TowerDelR3, hf_TowerDelR32, hf_TowerJetEt1, m_file, M_PI, matchedAllPt11, matchedAllPt12, matchedAllPt13, matchedAllPt21, matchedAllPt22, matchedAllPt23, matchedAllPt31, matchedAllPt32, matchedAllPt33, matchedAllPt41, matchedAllPt42, matchedAllPt43, matchedPt11, matchedPt12, matchedPt13, matchedPt21, matchedPt22, matchedPt23, matchedPt31, matchedPt32, matchedPt33, matchedPt41, matchedPt42, matchedPt43, MET1, MET12, MET13, nTowers1, nTowers2, nTowers3, nTowers4, nTowersLeadJet1, nTowersLeadJet2, nTowersLeadJet3, nTowersLeadJetPt1, nTowersLeadJetPt2, nTowersLeadJetPt3, nTowersLeadJetPt4, nTowersSecondJet1, nTowersSecondJet2, nTowersSecondJet3, ParMatch1, ParMatch2, ParMatch3, SumEt1, SumEt12, SumEt13, tbarMass, tbarMassGen, tMass, tMassGen, topMass1, topMass2, topMass3, topMassParton, totEneLeadJetEta1_1, totEneLeadJetEta1_2, totEneLeadJetEta2_1, totEneLeadJetEta2_2, totEneLeadJetEta3_1, totEneLeadJetEta3_2, TowerEtLeadJet1, TowerEtLeadJet12, TowerEtLeadJet13, TowerEtLeadJet2, TowerEtLeadJet22, TowerEtLeadJet23, TowerEtLeadJet3, TowerEtLeadJet32, TowerEtLeadJet33, ZpMass, ZpMass_2000_10, ZpMass_2000_13, ZpMass_2000_40, ZpMass_5000_10, ZpMass_5000_13, ZpMass_5000_40, ZpMass_700_10, ZpMass_700_13, ZpMass_700_40, ZpMassGen, ZpMassGen10, ZpMassGen13, ZpMassGen40, ZpMassGen_2000_10, ZpMassGen_2000_13, ZpMassGen_2000_40, ZpMassGen_5000_10, ZpMassGen_5000_13, ZpMassGen_5000_40, ZpMassGen_700_10, ZpMassGen_700_13, ZpMassGen_700_40, ZpMassMatched1, ZpMassMatched2, ZpMassMatched3, ZpMassRes101, ZpMassRes102, ZpMassRes103, ZpMassRes131, ZpMassRes132, ZpMassRes133, ZpMassRes401, ZpMassRes402, ZpMassRes403, ZpMassResL101, ZpMassResL102, ZpMassResL103, ZpMassResL131, ZpMassResL132, ZpMassResL133, ZpMassResL401, ZpMassResL402, ZpMassResL403, ZpMassResPRL101, ZpMassResPRL102, ZpMassResPRL103, ZpMassResPRL131, ZpMassResPRL132, ZpMassResPRL133, ZpMassResPRL401, ZpMassResPRL402, ZpMassResPRL403, ZpMassResRL101, ZpMassResRL102, ZpMassResRL103, ZpMassResRL131, ZpMassResRL132, ZpMassResRL133, ZpMassResRL401, ZpMassResRL402, ZpMassResRL403, ZpMassResRLoP101, ZpMassResRLoP102, ZpMassResRLoP103, ZpMassResRLoP131, ZpMassResRLoP132, ZpMassResRLoP133, ZpMassResRLoP401, ZpMassResRLoP402, and ZpMassResRLoP403.

51  {
52 
53  // Open the histogram file and book some associated histograms
54  m_file=new TFile("histo.root","RECREATE");
55 
56  tMassGen = TH1F("tMassGen","T Mass Gen",100,0,200);
57  tbarMassGen = TH1F("tbarMassGen","Tbar Mass Gen",100,0,200);
58 
59  tMass = TH1F("tMass","T Mass",100,0,200);
60  tbarMass = TH1F("tbarMass","Tbar Mass",100,0,200);
61 
62  topMassParton = TH1F("topMassParton","Top Mass Parton",100,0,200);
63  topMass1 = TH1F("topMass1","Top Mass 1",100,0,200);
64  topMass2 = TH1F("topMass2","Top Mass 2",100,0,200);
65  topMass3 = TH1F("topMass3","Top Mass 3",100,0,200);
66 
67  ZpMass = TH1F("ZpMass","Generated Zp Mass",160,0,8000);
68  ZpMassGen = TH1F("ZpMassGen","Gen Zp Mass",160,0,8000);
69  ZpMassMatched1 = TH1F("ZpMassMatched1","Calor Zp Mass 1",160,0,8000);
70  ZpMassMatched2 = TH1F("ZpMassMatched2","Calor Zp Mass 2",160,0,8000);
71  ZpMassMatched3 = TH1F("ZpMassMatched3","Calor Zp Mass 3",160,0,8000);
72 
73  ZpMassGen10 = TH1F("ZpMassGen10","Gen Zp Mass",160,0,8000);
74  ZpMassGen13 = TH1F("ZpMassGen13","Gen Zp Mass",160,0,8000);
75  ZpMassGen40 = TH1F("ZpMassGen40","Gen Zp Mass",160,0,8000);
76 
77  ZpMass_700_10 = TH1F("ZpMass_700_10","Parton Zp Mass",100,0,1000);
78  ZpMass_700_13 = TH1F("ZpMass_700_13","Parton Zp Mass",100,0,1000);
79  ZpMass_700_40 = TH1F("ZpMass_700_40","Parton Zp Mass",100,0,1000);
80 
81  ZpMassGen_700_10 = TH1F("ZpMassGen_700_10","Gen Zp Mass",100,0,1000);
82  ZpMassGen_700_13 = TH1F("ZpMassGen_700_13","Gen Zp Mass",100,0,1000);
83  ZpMassGen_700_40 = TH1F("ZpMassGen_700_40","Gen Zp Mass",100,0,1000);
84 
85  ZpMassGen_2000_10 = TH1F("ZpMassGen_2000_10","Gen Zp Mass",100,1500,2500);
86  ZpMassGen_2000_13 = TH1F("ZpMassGen_2000_13","Gen Zp Mass",100,1500,2500);
87  ZpMassGen_2000_40 = TH1F("ZpMassGen_2000_40","Gen Zp Mass",100,1500,2500);
88 
89  ZpMass_2000_10 = TH1F("ZpMass_2000_10","Parton Zp Mass",100,1500,2500);
90  ZpMass_2000_13 = TH1F("ZpMass_2000_13","Parton Zp Mass",100,1500,2500);
91  ZpMass_2000_40 = TH1F("ZpMass_2000_40","Parton Zp Mass",100,1500,2500);
92 
93  ZpMassGen_5000_10 = TH1F("ZpMassGen_5000_10","Gen Zp Mass",150,4000,5500);
94  ZpMassGen_5000_13 = TH1F("ZpMassGen_5000_13","Gen Zp Mass",150,4000,5500);
95  ZpMassGen_5000_40 = TH1F("ZpMassGen_5000_40","Gen Zp Mass",150,4000,5500);
96 
97  ZpMass_5000_10 = TH1F("ZpMass_5000_10","Parton Zp Mass",150,4000,5500);
98  ZpMass_5000_13 = TH1F("ZpMass_5000_13","Parton Zp Mass",150,4000,5500);
99  ZpMass_5000_40 = TH1F("ZpMass_5000_40","Parton Zp Mass",150,4000,5500);
100 
101  ZpMassRes101 = TH1F("ZpMassRes101","Zp Mass Resolution 1",100,-2,2);
102  ZpMassRes102 = TH1F("ZpMassRes102","Zp Mass Resolution 2",100,-2,2);
103  ZpMassRes103 = TH1F("ZpMassRes103","Zp Mass Resolution 3",100,-2,2);
104 
105  ZpMassRes131 = TH1F("ZpMassRes131","Zp Mass Resolution 1",100,-2,2);
106  ZpMassRes132 = TH1F("ZpMassRes132","Zp Mass Resolution 2",100,-2,2);
107  ZpMassRes133 = TH1F("ZpMassRes133","Zp Mass Resolution 3",100,-2,2);
108 
109  ZpMassRes401 = TH1F("ZpMassRes401","Zp Mass Resolution 1",100,-2,2);
110  ZpMassRes402 = TH1F("ZpMassRes402","Zp Mass Resolution 2",100,-2,2);
111  ZpMassRes403 = TH1F("ZpMassRes403","Zp Mass Resolution 3",100,-2,2);
112 
113  ZpMassResL101 = TH1F("ZpMassResL101","Zp Mass Resolution Leading Jets 1",100,0,2);
114  ZpMassResL102 = TH1F("ZpMassResL102","Zp Mass Resolution Leading Jets 2",100,0,2);
115  ZpMassResL103 = TH1F("ZpMassResL103","Zp Mass Resolution Leading Jets 3",100,0,2);
116 
117  ZpMassResRL101 = TH1F("ZpMassResRL101","Zp Mass Res. Ratio Leading Jets 1",100,0,2);
118  ZpMassResRL102 = TH1F("ZpMassResRL102","Zp Mass Res. Ratio Leading Jets 2",100,0,2);
119  ZpMassResRL103 = TH1F("ZpMassResRL103","Zp Mass Res. Ratio Leading Jets 3",100,0,2);
120 
121  ZpMassResRLoP101 = TH1F("ZpMassResRLoP101","Zp Mass RLoP Ratio Leading Jets 1",100,0,2);
122  ZpMassResRLoP102 = TH1F("ZpMassResRLoP102","Zp Mass RLoP Ratio Leading Jets 2",100,0,2);
123  ZpMassResRLoP103 = TH1F("ZpMassResRLoP103","Zp Mass RLoP Ratio Leading Jets 3",100,0,2);
124 
125  ZpMassResPRL101 = TH1F("ZpMassResPRL101","Zp Mass Res. P Ratio Leading Jets 1",100,0,2);
126  ZpMassResPRL102 = TH1F("ZpMassResPRL102","Zp Mass Res. P Ratio Leading Jets 2",100,0,2);
127  ZpMassResPRL103 = TH1F("ZpMassResPRL103","Zp Mass Res. P Ratio Leading Jets 3",100,0,2);
128 
129 
130  ZpMassResL131 = TH1F("ZpMassResL131","Zp Mass Resolution Leading Jets 1",100,0,2);
131  ZpMassResL132 = TH1F("ZpMassResL132","Zp Mass Resolution Leading Jets 2",100,0,2);
132  ZpMassResL133 = TH1F("ZpMassResL133","Zp Mass Resolution Leading Jets 3",100,0,2);
133 
134  ZpMassResRL131 = TH1F("ZpMassResRL131","Zp Mass Res. Ratio Leading Jets 1",100,0,2);
135  ZpMassResRL132 = TH1F("ZpMassResRL132","Zp Mass Res. Ratio Leading Jets 2",100,0,2);
136  ZpMassResRL133 = TH1F("ZpMassResRL133","Zp Mass Res. Ratio Leading Jets 3",100,0,2);
137 
138  ZpMassResRLoP131 = TH1F("ZpMassResRLoP131","Zp Mass RLoP Ratio Leading Jets 1",100,0,2);
139  ZpMassResRLoP132 = TH1F("ZpMassResRLoP132","Zp Mass RLoP Ratio Leading Jets 2",100,0,2);
140  ZpMassResRLoP133 = TH1F("ZpMassResRLoP133","Zp Mass RLoP Ratio Leading Jets 3",100,0,2);
141 
142  ZpMassResPRL131 = TH1F("ZpMassResPRL131","Zp Mass Res. P Ratio Leading Jets 1",100,0,2);
143  ZpMassResPRL132 = TH1F("ZpMassResPRL132","Zp Mass Res. P Ratio Leading Jets 2",100,0,2);
144  ZpMassResPRL133 = TH1F("ZpMassResPRL133","Zp Mass Res. P Ratio Leading Jets 3",100,0,2);
145 
146 
147  ZpMassResL401 = TH1F("ZpMassResL401","Zp Mass Resolution Leading Jets 1",100,0,2);
148  ZpMassResL402 = TH1F("ZpMassResL402","Zp Mass Resolution Leading Jets 2",100,0,2);
149  ZpMassResL403 = TH1F("ZpMassResL403","Zp Mass Resolution Leading Jets 3",100,0,2);
150 
151  ZpMassResRL401 = TH1F("ZpMassResRL401","Zp Mass Res. Ratio Leading Jets 1",100,0,2);
152  ZpMassResRL402 = TH1F("ZpMassResRL402","Zp Mass Res. Ratio Leading Jets 2",100,0,2);
153  ZpMassResRL403 = TH1F("ZpMassResRL403","Zp Mass Res. Ratio Leading Jets 3",100,0,2);
154 
155  ZpMassResRLoP401 = TH1F("ZpMassResRLoP401","Zp Mass RLoP Ratio Leading Jets 1",100,0,2);
156  ZpMassResRLoP402 = TH1F("ZpMassResRLoP402","Zp Mass RLoP Ratio Leading Jets 2",100,0,2);
157  ZpMassResRLoP403 = TH1F("ZpMassResRLoP403","Zp Mass RLoP Ratio Leading Jets 3",100,0,2);
158 
159  ZpMassResPRL401 = TH1F("ZpMassResPRL401","Zp Mass Res. P Ratio Leading Jets 1",100,0,2);
160  ZpMassResPRL402 = TH1F("ZpMassResPRL402","Zp Mass Res. P Ratio Leading Jets 2",100,0,2);
161  ZpMassResPRL403 = TH1F("ZpMassResPRL403","Zp Mass Res. P Ratio Leading Jets 3",100,0,2);
162 
163  dijetMass1 = TH1F("dijetMass1","DiJet Mass 1",100,0,4000);
164  dijetMass12 = TH1F("dijetMass12","DiJet Mass 1 2",100,0,6000);
165  dijetMass13 = TH1F("dijetMass13","DiJet Mass 1 3",100,0,12000);
166  dijetMass2 = TH1F("dijetMass2","DiJet Mass 2",100,0,4000);
167  dijetMass22 = TH1F("dijetMass22","DiJet Mass 2 2",100,0,6000);
168  dijetMass23 = TH1F("dijetMass23","DiJet Mass 2 3",100,0,12000);
169  dijetMass3 = TH1F("dijetMass3","DiJet Mass 3",100,0,4000);
170  dijetMass32 = TH1F("dijetMass32","DiJet Mass 3 2",100,0,6000);
171  dijetMass33 = TH1F("dijetMass33","DiJet Mass 3 3",100,0,12000);
172  dijetMass4 = TH1F("dijetMass4","DiJet Mass 4",100,0,4000);
173  dijetMass42 = TH1F("dijetMass42","DiJet Mass 4 2",100,0,6000);
174  dijetMass43 = TH1F("dijetMass43","DiJet Mass 4 3",100,0,12000);
175 
176  dijetMass101 = TH1F("dijetMass101","DiJet Mass 1",100,0,6000);
177  dijetMass131 = TH1F("dijetMass131","DiJet Mass 1",100,0,6000);
178  dijetMass401 = TH1F("dijetMass401","DiJet Mass 1",100,0,6000);
179 
180  dijetMass102 = TH1F("dijetMass102","DiJet Mass 2",100,0,6000);
181  dijetMass132 = TH1F("dijetMass132","DiJet Mass 2",100,0,6000);
182  dijetMass402 = TH1F("dijetMass402","DiJet Mass 2",100,0,6000);
183 
184  dijetMass103 = TH1F("dijetMass103","DiJet Mass 3",100,0,10000);
185  dijetMass133 = TH1F("dijetMass133","DiJet Mass 3",100,0,10000);
186  dijetMass403 = TH1F("dijetMass403","DiJet Mass 3",100,0,10000);
187 
188  dijetMass_700_101 = TH1F("dijetMass_700_101","DiJet Mass 1",100,0,1000);
189  dijetMass_700_131 = TH1F("dijetMass_700_131","DiJet Mass 1",100,0,1000);
190  dijetMass_700_401 = TH1F("dijetMass_700_401","DiJet Mass 1",100,0,1000);
191 
192  dijetMass_2000_101 = TH1F("dijetMass_2000_101","DiJet Mass 1",100,1500,2500);
193  dijetMass_2000_131 = TH1F("dijetMass_2000_131","DiJet Mass 1",100,1500,2500);
194  dijetMass_2000_401 = TH1F("dijetMass_2000_401","DiJet Mass 1",100,1500,2500);
195 
196  dijetMass_5000_101 = TH1F("dijetMass_5000_101","DiJet Mass 1",150,4000,5500);
197  dijetMass_5000_131 = TH1F("dijetMass_5000_131","DiJet Mass 1",150,4000,5500);
198  dijetMass_5000_401 = TH1F("dijetMass_5000_401","DiJet Mass 1",150,4000,5500);
199 
200 
201  dijetMassCor1 = TH1F("dijetMassCor1","DiJet Mass 1",160,0,8000);
202  dijetMassCor101 = TH1F("dijetMassCor101","DiJet Mass Cor 101",160,0,8000);
203  dijetMassCor131 = TH1F("dijetMassCor131","DiJet Mass Cor 131",160,0,8000);
204  dijetMassCor401 = TH1F("dijetMassCor401","DiJet Mass Cor 401",160,0,8000);
205 
206  dijetMassCor_700_1 = TH1F("dijetMassCor_700_1","DiJet Mass 1",100,0,1000);
207  dijetMassCor_700_101 = TH1F("dijetMassCor_700_101","DiJet Mass Cor 101",100,0,1000);
208  dijetMassCor_700_131 = TH1F("dijetMassCor_700_131","DiJet Mass Cor 131",100,0,1000);
209  dijetMassCor_700_401 = TH1F("dijetMassCor_700_401","DiJet Mass Cor 401",100,0,1000);
210 
211  dijetMassCor_2000_1 = TH1F("dijetMassCor_2000_1","DiJet Mass 1",100,1500,2500);
212  dijetMassCor_2000_101 = TH1F("dijetMassCor_2000_101","DiJet Mass Cor 101",100,1500,2500);
213  dijetMassCor_2000_131 = TH1F("dijetMassCor_2000_131","DiJet Mass Cor 131",100,1500,2500);
214  dijetMassCor_2000_401 = TH1F("dijetMassCor_2000_401","DiJet Mass Cor 401",100,1500,2500);
215 
216  dijetMassCor_5000_1 = TH1F("dijetMassCor_5000_1","DiJet Mass 1",150,4000,5500);
217  dijetMassCor_5000_101 = TH1F("dijetMassCor_5000_101","DiJet Mass Cor 101",150,4000,5500);
218  dijetMassCor_5000_131 = TH1F("dijetMassCor_5000_131","DiJet Mass Cor 131",150,4000,5500);
219  dijetMassCor_5000_401 = TH1F("dijetMassCor_5000_401","DiJet Mass Cor 401",150,4000,5500);
220 
221  dijetMassP1 = TH1F("dijetMassP1","DiJet Mass P 1",160,0,8000);
222  dijetMassP2 = TH1F("dijetMassP2","DiJet Mass P 2",160,0,8000);
223  dijetMassP3 = TH1F("dijetMassP3","DiJet Mass P 3",160,0,8000);
224 
225 
226  dijetMassP101 = TH1F("dijetMassP101","DiJet Mass P 1",160,0,8000);
227  dijetMassP131 = TH1F("dijetMassP131","DiJet Mass P 1",160,0,8000);
228  dijetMassP401 = TH1F("dijetMassP401","DiJet Mass P 1",160,0,8000);
229 
230  dijetMassP_700_101 = TH1F("dijetMassP_700_101","DiJet Mass P 1",100,0,1000);
231  dijetMassP_700_131 = TH1F("dijetMassP_700_131","DiJet Mass P 1",100,0,1000);
232  dijetMassP_700_401 = TH1F("dijetMassP_700_401","DiJet Mass P 1",100,0,1000);
233 
234  dijetMassP_2000_101 = TH1F("dijetMassP_2000_101","DiJet Mass P 1",100,1500,2500);
235  dijetMassP_2000_131 = TH1F("dijetMassP_2000_131","DiJet Mass P 1",100,1500,2500);
236  dijetMassP_2000_401 = TH1F("dijetMassP_2000_401","DiJet Mass P 1",100,1500,2500);
237 
238  dijetMassP_5000_101 = TH1F("dijetMassP_5000_101","DiJet Mass P 1",150,4000,5500);
239  dijetMassP_5000_131 = TH1F("dijetMassP_5000_131","DiJet Mass P 1",150,4000,5500);
240  dijetMassP_5000_401 = TH1F("dijetMassP_5000_401","DiJet Mass P 1",150,4000,5500);
241 
242  totEneLeadJetEta1_1 = TH1F("totEneLeadJetEta1_1","Total Energy Lead Jet Eta1 1",100,0,1500);
243  totEneLeadJetEta2_1 = TH1F("totEneLeadJetEta2_1","Total Energy Lead Jet Eta2 1",100,0,1500);
244  totEneLeadJetEta3_1 = TH1F("totEneLeadJetEta3_1","Total Energy Lead Jet Eta3 1",100,0,1500);
245  hadEneLeadJetEta1_1 = TH1F("hadEneLeadJetEta1_1","Hadronic Energy Lead Jet Eta1 1",100,0,1500);
246  hadEneLeadJetEta2_1 = TH1F("hadEneLeadJetEta2_1","Hadronic Energy Lead Jet Eta2 1",100,0,1500);
247  hadEneLeadJetEta3_1 = TH1F("hadEneLeadJetEta3_1","Hadronic Energy Lead Jet Eta3 1",100,0,1500);
248  emEneLeadJetEta1_1 = TH1F("emEneLeadJetEta1_1","EM Energy Lead Jet Eta1 1",100,0,1500);
249  emEneLeadJetEta2_1 = TH1F("emEneLeadJetEta2_1","EM Energy Lead Jet Eta2 1",100,0,1500);
250  emEneLeadJetEta3_1 = TH1F("emEneLeadJetEta3_1","EM Energy Lead Jet Eta3 1",100,0,1500);
251 
252  totEneLeadJetEta1_2 = TH1F("totEneLeadJetEta1_2","Total Energy Lead Jet Eta1 2",100,0,6000);
253  totEneLeadJetEta2_2 = TH1F("totEneLeadJetEta2_2","Total Energy Lead Jet Eta2 2",100,0,6000);
254  totEneLeadJetEta3_2 = TH1F("totEneLeadJetEta3_2","Total Energy Lead Jet Eta3 2",100,0,6000);
255  hadEneLeadJetEta1_2 = TH1F("hadEneLeadJetEta1_2","Hadronic Energy Lead Jet Eta1 2",100,0,6000);
256  hadEneLeadJetEta2_2 = TH1F("hadEneLeadJetEta2_2","Hadronic Energy Lead Jet Eta2 2",100,0,6000);
257  hadEneLeadJetEta3_2 = TH1F("hadEneLeadJetEta3_2","Hadronic Energy Lead Jet Eta3 2",100,0,6000);
258  emEneLeadJetEta1_2 = TH1F("emEneLeadJetEta1_2","EM Energy Lead Jet Eta1 2",100,0,5000);
259  emEneLeadJetEta2_2 = TH1F("emEneLeadJetEta2_2","EM Energy Lead Jet Eta2 2",100,0,5000);
260  emEneLeadJetEta3_2 = TH1F("emEneLeadJetEta3_2","EM Energy Lead Jet Eta3 2",100,0,5000);
261 
262  hadEneLeadJet1 = TH1F("hadEneLeadJet1","Hadronic Energy Lead Jet 1",100,0,3000);
263  hadEneLeadJet12 = TH1F("hadEneLeadJet12","Hadronic Energy Lead Jet 1 2",100,0,4000);
264  hadEneLeadJet13 = TH1F("hadEneLeadJet13","Hadronic Energy Lead Jet 1 3",100,0,6000);
265  hadEneLeadJet2 = TH1F("hadEneLeadJet2","Hadronic Energy Lead Jet 2",100,0,3000);
266  hadEneLeadJet22 = TH1F("hadEneLeadJet22","Hadronic Energy Lead Jet 2 2",100,0,4000);
267  hadEneLeadJet23 = TH1F("hadEneLeadJet23","Hadronic Energy Lead Jet 2 3",100,0,6000);
268  hadEneLeadJet3 = TH1F("hadEneLeadJet3","Hadronic Energy Lead Jet 3",100,0,3000);
269  hadEneLeadJet32 = TH1F("hadEneLeadJet32","Hadronic Energy Lead Jet 3 2",100,0,4000);
270  hadEneLeadJet33 = TH1F("hadEneLeadJet33","Hadronic Energy Lead Jet 3 3",100,0,6000);
271 
272  emEneLeadJet1 = TH1F("emEneLeadJet1","EM Energy Lead Jet 1",100,0,1500);
273  emEneLeadJet12 = TH1F("emEneLeadJet12","EM Energy Lead Jet 1 2",100,0,3000);
274  emEneLeadJet13 = TH1F("emEneLeadJet13","EM Energy Lead Jet 1 3",100,0,5000);
275  emEneLeadJet2 = TH1F("emEneLeadJet2","EM Energy Lead Jet 2",100,0,1500);
276  emEneLeadJet22 = TH1F("emEneLeadJet22","EM Energy Lead Jet 2 2",100,0,3000);
277  emEneLeadJet23 = TH1F("emEneLeadJet23","EM Energy Lead Jet 2 3",100,0,5000);
278  emEneLeadJet3 = TH1F("emEneLeadJet3","EM Energy Lead Jet 3",100,0,1500);
279  emEneLeadJet32 = TH1F("emEneLeadJet32","EM Energy Lead Jet 3 2",100,0,3000);
280  emEneLeadJet33 = TH1F("emEneLeadJet33","EM Energy Lead Jet 3 3",100,0,5000);
281 
282  hadFracEta11 = TH1F("hadFracEta11","Hadronic Fraction Eta1 Jet 1",100,0,1);
283  hadFracEta21 = TH1F("hadFracEta21","Hadronic Fraction Eta2 Jet 1",100,0,1);
284  hadFracEta31 = TH1F("hadFracEta31","Hadronic Fraction Eta3 Jet 1",100,0,1);
285 
286  hadFracEta12 = TH1F("hadFracEta12","Hadronic Fraction Eta1 Jet 2",100,0,1);
287  hadFracEta22 = TH1F("hadFracEta22","Hadronic Fraction Eta2 Jet 2",100,0,1);
288  hadFracEta32 = TH1F("hadFracEta32","Hadronic Fraction Eta3 Jet 2",100,0,1);
289 
290  hadFracEta13 = TH1F("hadFracEta13","Hadronic Fraction Eta1 Jet 3",100,0,1);
291  hadFracEta23 = TH1F("hadFracEta23","Hadronic Fraction Eta2 Jet 3",100,0,1);
292  hadFracEta33 = TH1F("hadFracEta33","Hadronic Fraction Eta3 Jet 3",100,0,1);
293 
294  hadFracLeadJet1 = TH1F("hadFracLeadJet1","Hadronic Fraction Lead Jet 1",100,0,1);
295  hadFracLeadJet2 = TH1F("hadFracLeadJet2","Hadronic Fraction Lead Jet 2",100,0,1);
296  hadFracLeadJet3 = TH1F("hadFracLeadJet3","Hadronic Fraction Lead Jet 3",100,0,1);
297 
298  SumEt1 = TH1F("SumEt1","SumEt 1",100,0,1000);
299  SumEt12 = TH1F("SumEt12","SumEt 1 2",100,0,4000);
300  SumEt13 = TH1F("SumEt13","SumEt 1 3",100,0,15000);
301 
302  MET1 = TH1F("MET1", "MET 1",100,0,200);
303  MET12 = TH1F("MET12", "MET 1 2",100,0,1000);
304  MET13 = TH1F("MET13", "MET 1 3",100,0,3000);
305 
306  nTowersLeadJet1 = TH1F("nTowersLeadJet1","Number of Towers Lead Jet 1",100,0,100);
307  nTowersLeadJet2 = TH1F("nTowersLeadJet2","Number of Towers Lead Jet 2",100,0,100);
308  nTowersLeadJet3 = TH1F("nTowersLeadJet3","Number of Towers Lead Jet 3",100,0,100);
309 
310  nTowersSecondJet1 = TH1F("nTowersSecondJet1","Number of Towers Second Jet 1",100,0,100);
311  nTowersSecondJet2 = TH1F("nTowersSecondJet2","Number of Towers Second Jet 2",100,0,100);
312  nTowersSecondJet3 = TH1F("nTowersSecondJet3","Number of Towers Second Jet 3",100,0,100);
313 
314  hf_PtResponse1 = TProfile("PtResponse1","Pt Response 1", 100, -5, 5, 0, 10);
315  hf_PtResponse2 = TProfile("PtResponse2","Pt Response 2", 100, -5, 5, 0, 10);
316  hf_PtResponse3 = TProfile("PtResponse3","Pt Response 3", 100, -5, 5, 0, 10);
317  hf_PtResponse4 = TProfile("PtResponse4","Pt Response 4", 100, -5, 5, 0, 10);
318 
319  hf_TowerDelR1 = TProfile("hf_TowerDelR1","Tower Del R 1", 100, 0, 2, 0, 10);
320  hf_TowerDelR12 = TProfile("hf_TowerDelR12","Tower Del R 1", 80, 0, 0.8, 0, 10);
321  hf_TowerDelR2 = TProfile("hf_TowerDelR2","Tower Del R 2", 100, 0, 2, 0, 10);
322  hf_TowerDelR22 = TProfile("hf_TowerDelR22","Tower Del R 2", 80, 0, 0.8, 0, 10);
323  hf_TowerDelR3 = TProfile("hf_TowerDelR3","Tower Del R 3", 100, 0, 2, 0, 10);
324  hf_TowerDelR32 = TProfile("hf_TowerDelR32","Tower Del R 3", 80, 0, 0.8, 0, 10);
325 
326  hf_sumTowerAllEx = TH1F("sumTowerAllEx","Tower Ex",100,-1000,1000);
327  hf_sumTowerAllEy = TH1F("sumTowerAllEy","Tower Ey",100,-1000,1000);
328 
329  hf_TowerJetEt1 = TH1F("TowerJetEt1","Tower/Jet Et 1",50,0,1);
330 
331  nTowers1 = TH1F("nTowers1","Number of Towers pt 0.5",100,0,500);
332  nTowers2 = TH1F("nTowers2","Number of Towers pt 1.0",100,0,500);
333  nTowers3 = TH1F("nTowers3","Number of Towers pt 1.5",100,0,500);
334  nTowers4 = TH1F("nTowers4","Number of Towers pt 2.0",100,0,500);
335 
336  nTowersLeadJetPt1 = TH1F("nTowersLeadJetPt1","Number of Towers in Lead Jet pt 0.5",100,0,200);
337  nTowersLeadJetPt2 = TH1F("nTowersLeadJetPt2","Number of Towers in Lead Jet pt 1.0",100,0,200);
338  nTowersLeadJetPt3 = TH1F("nTowersLeadJetPt3","Number of Towers in Lead Jet pt 1.5",100,0,200);
339  nTowersLeadJetPt4 = TH1F("nTowersLeadJetPt4","Number of Towers in Lead Jet pt 2.0",100,0,200);
340 
341  TowerEtLeadJet1 = TH1F("TowerEtLeadJet1","Towers Et Lead Jet 1",100,0,2000);
342  TowerEtLeadJet12 = TH1F("TowerEtLeadJet12","Towers Et Lead Jet 1 2",100,0,6000);
343  TowerEtLeadJet13 = TH1F("TowerEtLeadJet13","Towers Et Lead Jet 1 3",100,0,300);
344  TowerEtLeadJet2 = TH1F("TowerEtLeadJet2","Towers Et Lead Jet 2",100,0,2000);
345  TowerEtLeadJet22 = TH1F("TowerEtLeadJet22","Towers Et Lead Jet 2 2",100,0,6000);
346  TowerEtLeadJet23 = TH1F("TowerEtLeadJet23","Towers Et Lead Jet 2 3",100,0,300);
347  TowerEtLeadJet3 = TH1F("TowerEtLeadJet3","Towers Et Lead Jet 3",100,0,2000);
348  TowerEtLeadJet32 = TH1F("TowerEtLeadJet32","Towers Et Lead Jet 3 2",100,0,6000);
349  TowerEtLeadJet33 = TH1F("TowerEtLeadJet33","Towers Et Lead Jet 3 3",100,0,300);
350 
351  hf_nJet1 = TProfile("hf_nJet1", "Num Jets 1", 100, 0, 5000, 0, 50);
352  hf_nJet2 = TProfile("hf_nJet2", "Num Jets 2", 100, 0, 5000, 0, 50);
353  hf_nJet3 = TProfile("hf_nJet3", "Num Jets 3", 100, 0, 5000, 0, 50);
354  hf_nJet4 = TProfile("hf_nJet4", "Num Jets 4", 100, 0, 5000, 0, 50);
355 
356  hf_nJet1s = TProfile("hf_nJet1s", "Num Jets 1", 100, 0, 200, 0, 50);
357  hf_nJet2s = TProfile("hf_nJet2s", "Num Jets 2", 100, 0, 200, 0, 50);
358  hf_nJet3s = TProfile("hf_nJet3s", "Num Jets 3", 100, 0, 200, 0, 50);
359  hf_nJet4s = TProfile("hf_nJet4s", "Num Jets 4", 100, 0, 200, 0, 50);
360 
361  hf_nJet11 = TProfile("hf_nJet11", "Num Jets 1 1", 100, 0, 3000, 0, 50);
362  hf_nJet21 = TProfile("hf_nJet21", "Num Jets 2 1", 100, 0, 3000, 0, 50);
363  hf_nJet31 = TProfile("hf_nJet31", "Num Jets 3 1", 100, 0, 3000, 0, 50);
364  hf_nJet41 = TProfile("hf_nJet41", "Num Jets 4 1", 100, 0, 3000, 0, 50);
365 
366  dRPar1 = TH1F("dRPar1","Parton dR with matched CaloJet1",100,0,0.5);
367  dPhiPar1 = TH1F("dPhiPar1","Parton dPhi with matched CaloJet1",200,-0.5,0.5);
368  dEtaPar1 = TH1F("dEtaPar1","Parton dEta with matched CaloJet1",200,-0.5,0.5);
369  dPtPar1 = TH1F("dPtPar1","Parton dPt with matched CaloJet1",200,-50,50);
370 
371  dRPar2 = TH1F("dRPar2","Parton dR with matched CaloJet2",100,0,0.5);
372  dPhiPar2 = TH1F("dPhiPar2","Parton dPhi with matched CaloJet2",200,-0.5,0.5);
373  dEtaPar2 = TH1F("dEtaPar2","Parton dEta with matched CaloJet2",200,-0.5,0.5);
374  dPtPar2 = TH1F("dPtPar2","Parton dPt with matched CaloJet2",200,-50,50);
375 
376  dRPar3 = TH1F("dRPar3","Parton dR with matched CaloJet3",100,0,0.5);
377  dPhiPar3 = TH1F("dPhiPar3","Parton dPhi with matched CaloJet3",200,-0.5,0.5);
378  dEtaPar3 = TH1F("dEtaPar3","Parton dEta with matched CaloJet3",200,-0.5,0.5);
379  dPtPar3 = TH1F("dPtPar3","Parton dPt with matched CaloJet3",200,-50,50);
380 
381  dRPar4 = TH1F("dRPar4","Parton dR with matched CaloJet4",100,0,0.5);
382  dPhiPar4 = TH1F("dPhiPar4","Parton dPhi with matched CaloJet4",200,-0.5,0.5);
383  dEtaPar4 = TH1F("dEtaPar4","Parton dEta with matched CaloJet4",200,-0.5,0.5);
384  dPtPar4 = TH1F("dPtPar4","Parton dPt with matched CaloJet4",200,-50,50);
385 
386  dRParton = TH1F("dRParton","dR Parton",100,0,10.0);
387  dRPartonMin = TH1F("dRPartonMin","Min dR Parton",100,0,2.0);
388 
389  dR1 = TH1F("dR1","GenJets dR with matched CaloJet",100,0,0.5);
390  dPhi1 = TH1F("dPhi1","GenJets dPhi with matched CaloJet",200,-0.5,0.5);
391  dEta1 = TH1F("dEta1","GenJets dEta with matched CaloJet",200,-0.5,0.5);
392  dPt1 = TH1F("dPt1","GenJets dPt with matched CaloJet",200,-100,100);
393  dPtFrac1 = TH1F("dPtFrac1","GenJets dPt frac with matched CaloJet",100,-1,1);
394  dPt20Frac1 = TH1F("dPt20Frac1","GenJets dPt frac with matched CaloJet",100,-1,1);
395  dPt40Frac1 = TH1F("dPt40Frac1","GenJets dPt frac with matched CaloJet",100,-1,1);
396  dPt80Frac1 = TH1F("dPt80Frac1","GenJets dPt frac with matched CaloJet",100,-1,1);
397  dPt100Frac1 = TH1F("dPt100Frac1","GenJets dPt frac with matched CaloJet",100,-1,1);
398 
399  dR2 = TH1F("dR2","GenJets dR with matched CaloJet",100,0,0.5);
400  dPhi2 = TH1F("dPhi2","GenJets dPhi with matched CaloJet",200,-0.5,0.5);
401  dEta2 = TH1F("dEta2","GenJets dEta with matched CaloJet",200,-0.5,0.5);
402  dPt2 = TH1F("dPt2","GenJets dPt with matched CaloJet",200,-100,100);
403  dPtFrac2 = TH1F("dPtFrac2","GenJets dPt frac with matched CaloJet",100,-1,1);
404  dPt20Frac2 = TH1F("dPt20Frac2","GenJets dPt frac with matched CaloJet",100,-1,1);
405  dPt40Frac2 = TH1F("dPt40Frac2","GenJets dPt frac with matched CaloJet",100,-1,1);
406  dPt80Frac2 = TH1F("dPt80Frac2","GenJets dPt frac with matched CaloJet",100,-1,1);
407  dPt100Frac2 = TH1F("dPt100Frac2","GenJets dPt frac with matched CaloJet",100,-1,1);
408 
409  dR3 = TH1F("dR3","GenJets dR with matched CaloJet",100,0,0.5);
410  dPhi3 = TH1F("dPhi3","GenJets dPhi with matched CaloJet",200,-0.5,0.5);
411  dEta3 = TH1F("dEta3","GenJets dEta with matched CaloJet",200,-0.5,0.5);
412  dPt3 = TH1F("dPt3","GenJets dPt with matched CaloJet",200,-100,100);
413  dPtFrac3 = TH1F("dPtFrac3","GenJets dPt frac with matched CaloJet",100,-1,1);
414  dPt20Frac3 = TH1F("dPt20Frac3","GenJets dPt frac with matched CaloJet",100,-1,1);
415  dPt40Frac3 = TH1F("dPt40Frac3","GenJets dPt frac with matched CaloJet",100,-1,1);
416  dPt80Frac3 = TH1F("dPt80Frac3","GenJets dPt frac with matched CaloJet",100,-1,1);
417  dPt100Frac3 = TH1F("dPt100Frac3","GenJets dPt frac with matched CaloJet",100,-1,1);
418 
419  dR4 = TH1F("dR4","GenJets dR with matched CaloJet",100,0,0.5);
420  dPhi4 = TH1F("dPhi4","GenJets dPhi with matched CaloJet",200,-0.5,0.5);
421  dEta4 = TH1F("dEta4","GenJets dEta with matched CaloJet",200,-0.5,0.5);
422  dPt4 = TH1F("dPt4","GenJets dPt with matched CaloJet",200,-100,100);
423  dPtFrac4 = TH1F("dPtFrac4","GenJets dPt frac with matched CaloJet",100,-1,1);
424  dPt20Frac4 = TH1F("dPt20Frac4","GenJets dPt frac with matched CaloJet",100,-1,1);
425  dPt40Frac4 = TH1F("dPt40Frac4","GenJets dPt frac with matched CaloJet",100,-1,1);
426  dPt80Frac4 = TH1F("dPt80Frac4","GenJets dPt frac with matched CaloJet",100,-1,1);
427  dPt100Frac4 = TH1F("dPt100Frac4","GenJets dPt frac with matched CaloJet",100,-1,1);
428 
429  dR12 = TH1F("dR12","dR MidPoint - SISCone",100,0,0.5);
430  dPhi12 = TH1F("dPhi12","dPhi MidPoint - SISCone",200,-0.5,0.5);
431  dEta12 = TH1F("dEta12","dEta MidPoint - SISCone",200,-0.5,0.5);
432  dPt12 = TH1F("dPt12","dPt MidPoint - SISCone",200,-100,100);
433 
434 
435 
436  h_nCalJets1 = TH1F( "nCalJets1", "Number of CalJets1", 20, 0, 20 );
437  h_nCalJets2 = TH1F( "nCalJets2", "Number of CalJets2", 20, 0, 20 );
438  h_nCalJets3 = TH1F( "nCalJets3", "Number of CalJets3", 20, 0, 20 );
439  h_nCalJets4 = TH1F( "nCalJets4", "Number of CalJets4", 20, 0, 20 );
440 
441  h_lowPtCal11 = TH1F( "lowPtCal11", "Low p_{T} of CalJet1 Eta1", 100, 0, 100 );
442  h_lowPtCal21 = TH1F( "lowPtCal21", "Low p_{T} of CalJet2 Eta1", 100, 0, 100 );
443  h_lowPtCal31 = TH1F( "lowPtCal31", "Low p_{T} of CalJet3 Eta1", 100, 0, 100 );
444  h_lowPtCal41 = TH1F( "lowPtCal41", "Low p_{T} of CalJet4 Eta1", 100, 0, 100 );
445 
446  h_lowPtCal12 = TH1F( "lowPtCal12", "Low p_{T} of CalJet1 Eta2", 100, 0, 100 );
447  h_lowPtCal22 = TH1F( "lowPtCal22", "Low p_{T} of CalJet2 Eta2", 100, 0, 100 );
448  h_lowPtCal32 = TH1F( "lowPtCal32", "Low p_{T} of CalJet3 Eta2", 100, 0, 100 );
449  h_lowPtCal42 = TH1F( "lowPtCal42", "Low p_{T} of CalJet4 Eta2", 100, 0, 100 );
450 
451  h_lowPtCal13 = TH1F( "lowPtCal13", "Low p_{T} of CalJet1 Eta3", 100, 0, 100 );
452  h_lowPtCal23 = TH1F( "lowPtCal23", "Low p_{T} of CalJet2 Eta3", 100, 0, 100 );
453  h_lowPtCal33 = TH1F( "lowPtCal33", "Low p_{T} of CalJet3 Eta3", 100, 0, 100 );
454  h_lowPtCal43 = TH1F( "lowPtCal43", "Low p_{T} of CalJet4 Eta3", 100, 0, 100 );
455 
456  h_lowPtCal1c11 = TH1F( "lowPtCal1c11", "Low p_{T} of CalJet1 c1 Eta1", 100, 0, 100 );
457  h_lowPtCal2c11 = TH1F( "lowPtCal2c11", "Low p_{T} of CalJet2 c1 Eta1", 100, 0, 100 );
458  h_lowPtCal3c11 = TH1F( "lowPtCal3c11", "Low p_{T} of CalJet3 c1 Eta1", 100, 0, 100 );
459  h_lowPtCal4c11 = TH1F( "lowPtCal4c11", "Low p_{T} of CalJet4 c1 Eta1", 100, 0, 100 );
460 
461  h_lowPtCal1c12 =