|
|
|
#include <Calibration/IsolatedParticles/plugins/IsolatedTracksNxN.cc>
Public Member Functions | |
| IsolatedTracksNxN (const edm::ParameterSet &) | |
| ~IsolatedTracksNxN () | |
Private Member Functions | |
| void | analyze (const edm::Event &, const edm::EventSetup &) |
| void | beginJob () |
| void | BookHistograms () |
| void | clearTreeVectors () |
| double | DeltaPhi (double v1, double v2) |
| double | DeltaR (double eta1, double phi1, double eta2, double phi2) |
| void | endJob () |
| void | printTrack (const reco::Track *pTrack) |
Private Attributes | |
| std::string | algoBitToName [nL1BitsMax] |
| const MagneticField * | bField |
| bool | debugL1Info_ |
| int | debugTrks_ |
| bool | doMC |
| edm::Service< TFileService > | fs |
| double | genPartEtaBins [4] |
| double | genPartPBins [16] |
| TH1I * | h_L1AlgoNames |
| TH1F * | h_maxNearP15x15 [NPBins][NEtaBins] |
| TH1F * | h_maxNearP21x21 [NPBins][NEtaBins] |
| TH1F * | h_maxNearP25x25 [NPBins][NEtaBins] |
| TH1F * | h_maxNearP31x31 [NPBins][NEtaBins] |
| TH1F * | h_nTracks |
| TH1F * | h_PVTracksWt |
| TH1F * | h_recEta_0 |
| TH1F * | h_recEta_1 |
| TH1F * | h_recEta_2 |
| TH2F * | h_recEtaP_0 |
| TH2F * | h_recEtaP_1 |
| TH2F * | h_recEtaP_2 |
| TH2F * | h_recEtaPt_0 |
| TH2F * | h_recEtaPt_1 |
| TH2F * | h_recEtaPt_2 |
| TH1F * | h_recP_0 |
| TH1F * | h_recP_1 |
| TH1F * | h_recP_2 |
| TH1F * | h_recPhi_0 |
| TH1F * | h_recPhi_1 |
| TH1F * | h_recPhi_2 |
| TH1F * | h_recPt_0 |
| TH1F * | h_recPt_1 |
| TH1F * | h_recPt_2 |
| bool | initL1 |
| edm::InputTag | JetExtender_ |
| edm::InputTag | JetSrc_ |
| edm::InputTag | L1extraCenJetSource_ |
| edm::InputTag | L1extraFwdJetSource_ |
| edm::InputTag | L1extraIsoEmSource_ |
| edm::InputTag | L1extraMuonSource_ |
| edm::InputTag | L1extraNonIsoEmSource_ |
| edm::InputTag | L1extraTauJetSource_ |
| edm::InputTag | L1GTObjectMapRcdSource_ |
| edm::InputTag | L1GTReadoutRcdSource_ |
| std::map< std::string, bool > | l1TriggerMap |
| double | maxTrackEta_ |
| double | minTrackP_ |
| int | myverbose_ |
| int | nEventProc |
| bool | printTrkHitPattern_ |
| double | pvTracksPtMin_ |
| int | t_Bunch |
| std::vector< double > * | t_e11x11 |
| std::vector< double > * | t_e11x11_10Sig |
| std::vector< double > * | t_e11x11_15Sig |
| std::vector< double > * | t_e11x11_20Sig |
| std::vector< double > * | t_e11x11_25Sig |
| std::vector< double > * | t_e11x11_30Sig |
| std::vector< double > * | t_e13x13 |
| std::vector< double > * | t_e15x15 |
| std::vector< double > * | t_e15x15_10Sig |
| std::vector< double > * | t_e15x15_15Sig |
| std::vector< double > * | t_e15x15_20Sig |
| std::vector< double > * | t_e15x15_25Sig |
| std::vector< double > * | t_e15x15_30Sig |
| std::vector< double > * | t_e21x21 |
| std::vector< double > * | t_e25x25 |
| std::vector< double > * | t_e31x31 |
| std::vector< double > * | t_e3x3 |
| std::vector< double > * | t_e5x5 |
| std::vector< double > * | t_e7x7 |
| std::vector< double > * | t_e7x7_10Sig |
| std::vector< double > * | t_e7x7_15Sig |
| std::vector< double > * | t_e7x7_20Sig |
| std::vector< double > * | t_e7x7_25Sig |
| std::vector< double > * | t_e7x7_30Sig |
| std::vector< double > * | t_e9x9 |
| std::vector< double > * | t_e9x9_10Sig |
| std::vector< double > * | t_e9x9_15Sig |
| std::vector< double > * | t_e9x9_20Sig |
| std::vector< double > * | t_e9x9_25Sig |
| std::vector< double > * | t_e9x9_30Sig |
| std::vector< int > * | t_ecalSpike11x11 |
| std::vector< double > * | t_esim11x11 |
| std::vector< double > * | t_esim11x11CharHad |
| std::vector< double > * | t_esim11x11Matched |
| std::vector< double > * | t_esim11x11NeutHad |
| std::vector< double > * | t_esim11x11Photon |
| std::vector< double > * | t_esim11x11Rest |
| std::vector< double > * | t_esim13x13 |
| std::vector< double > * | t_esim13x13CharHad |
| std::vector< double > * | t_esim13x13Matched |
| std::vector< double > * | t_esim13x13NeutHad |
| std::vector< double > * | t_esim13x13Photon |
| std::vector< double > * | t_esim13x13Rest |
| std::vector< double > * | t_esim15x15 |
| std::vector< double > * | t_esim15x15CharHad |
| std::vector< double > * | t_esim15x15Matched |
| std::vector< double > * | t_esim15x15NeutHad |
| std::vector< double > * | t_esim15x15Photon |
| std::vector< double > * | t_esim15x15Rest |
| std::vector< double > * | t_esim21x21 |
| std::vector< double > * | t_esim21x21CharHad |
| std::vector< double > * | t_esim21x21Matched |
| std::vector< double > * | t_esim21x21NeutHad |
| std::vector< double > * | t_esim21x21Photon |
| std::vector< double > * | t_esim21x21Rest |
| std::vector< double > * | t_esim25x25 |
| std::vector< double > * | t_esim25x25CharHad |
| std::vector< double > * | t_esim25x25Matched |
| std::vector< double > * | t_esim25x25NeutHad |
| std::vector< double > * | t_esim25x25Photon |
| std::vector< double > * | t_esim25x25Rest |
| std::vector< double > * | t_esim31x31 |
| std::vector< double > * | t_esim31x31CharHad |
| std::vector< double > * | t_esim31x31Matched |
| std::vector< double > * | t_esim31x31NeutHad |
| std::vector< double > * | t_esim31x31Photon |
| std::vector< double > * | t_esim31x31Rest |
| std::vector< double > * | t_esim3x3 |
| std::vector< double > * | t_esim3x3CharHad |
| std::vector< double > * | t_esim3x3Matched |
| std::vector< double > * | t_esim3x3NeutHad |
| std::vector< double > * | t_esim3x3Photon |
| std::vector< double > * | t_esim3x3Rest |
| std::vector< double > * | t_esim5x5 |
| std::vector< double > * | t_esim5x5CharHad |
| std::vector< double > * | t_esim5x5Matched |
| std::vector< double > * | t_esim5x5NeutHad |
| std::vector< double > * | t_esim5x5Photon |
| std::vector< double > * | t_esim5x5Rest |
| std::vector< double > * | t_esim7x7 |
| std::vector< double > * | t_esim7x7CharHad |
| std::vector< double > * | t_esim7x7Matched |
| std::vector< double > * | t_esim7x7NeutHad |
| std::vector< double > * | t_esim7x7Photon |
| std::vector< double > * | t_esim7x7Rest |
| std::vector< double > * | t_esim9x9 |
| std::vector< double > * | t_esim9x9CharHad |
| std::vector< double > * | t_esim9x9Matched |
| std::vector< double > * | t_esim9x9NeutHad |
| std::vector< double > * | t_esim9x9Photon |
| std::vector< double > * | t_esim9x9Rest |
| std::vector< double > * | t_esimPdgId |
| int | t_EvtNo |
| std::vector< double > * | t_h3x3 |
| std::vector< double > * | t_h3x3Sig |
| std::vector< double > * | t_h5x5 |
| std::vector< double > * | t_h5x5Sig |
| std::vector< double > * | t_h7x7 |
| std::vector< double > * | t_h7x7Sig |
| std::vector< double > * | t_hsim3x3 |
| std::vector< double > * | t_hsim3x3CharHad |
| std::vector< double > * | t_hsim3x3Matched |
| std::vector< double > * | t_hsim3x3NeutHad |
| std::vector< double > * | t_hsim3x3Photon |
| std::vector< double > * | t_hsim3x3Rest |
| std::vector< double > * | t_hsim5x5 |
| std::vector< double > * | t_hsim5x5CharHad |
| std::vector< double > * | t_hsim5x5Matched |
| std::vector< double > * | t_hsim5x5NeutHad |
| std::vector< double > * | t_hsim5x5Photon |
| std::vector< double > * | t_hsim5x5Rest |
| std::vector< double > * | t_hsim7x7 |
| std::vector< double > * | t_hsim7x7CharHad |
| std::vector< double > * | t_hsim7x7Matched |
| std::vector< double > * | t_hsim7x7NeutHad |
| std::vector< double > * | t_hsim7x7Photon |
| std::vector< double > * | t_hsim7x7Rest |
| std::vector< int > * | t_infoHcal |
| std::vector< double > * | t_jetEta |
| std::vector< double > * | t_jetPhi |
| std::vector< double > * | t_jetPt |
| std::vector< double > * | t_L1CenJetEta |
| std::vector< double > * | t_L1CenJetPhi |
| std::vector< double > * | t_L1CenJetPt |
| std::vector< int > * | t_L1Decision |
| std::vector< double > * | t_L1FwdJetEta |
| std::vector< double > * | t_L1FwdJetPhi |
| std::vector< double > * | t_L1FwdJetPt |
| std::vector< double > * | t_L1IsoEMEta |
| std::vector< double > * | t_L1IsoEMPhi |
| std::vector< double > * | t_L1IsoEMPt |
| std::vector< double > * | t_L1METEta |
| std::vector< double > * | t_L1METPhi |
| std::vector< double > * | t_L1METPt |
| std::vector< double > * | t_L1MuonEta |
| std::vector< double > * | t_L1MuonPhi |
| std::vector< double > * | t_L1MuonPt |
| std::vector< double > * | t_L1NonIsoEMEta |
| std::vector< double > * | t_L1NonIsoEMPhi |
| std::vector< double > * | t_L1NonIsoEMPt |
| std::vector< double > * | t_L1TauJetEta |
| std::vector< double > * | t_L1TauJetPhi |
| std::vector< double > * | t_L1TauJetPt |
| int | t_Lumi |
| std::vector< double > * | t_maxNearHcalP3x3 |
| std::vector< double > * | t_maxNearHcalP5x5 |
| std::vector< double > * | t_maxNearHcalP7x7 |
| std::vector< double > * | t_maxNearP11x11 |
| std::vector< double > * | t_maxNearP13x13 |
| std::vector< double > * | t_maxNearP15x15 |
| std::vector< double > * | t_maxNearP21x21 |
| std::vector< double > * | t_maxNearP25x25 |
| std::vector< double > * | t_maxNearP31x31 |
| std::vector< double > * | t_maxNearP7x7 |
| std::vector< double > * | t_maxNearP9x9 |
| std::vector< int > * | t_NLayersCrossed |
| int | t_nTracks |
| std::vector< double > * | t_nTrksJetCalo |
| std::vector< double > * | t_nTrksJetVtx |
| std::vector< double > * | t_PVisValid |
| std::vector< double > * | t_PVndof |
| std::vector< double > * | t_PVNTracks |
| std::vector< double > * | t_PVNTracksHP |
| std::vector< double > * | t_PVNTracksHPWt |
| std::vector< double > * | t_PVNTracksWt |
| std::vector< double > * | t_PVTracksSumPt |
| std::vector< double > * | t_PVTracksSumPtHP |
| std::vector< double > * | t_PVTracksSumPtHPWt |
| std::vector< double > * | t_PVTracksSumPtWt |
| std::vector< double > * | t_PVx |
| std::vector< double > * | t_PVy |
| std::vector< double > * | t_PVz |
| int | t_RunNo |
| std::vector< double > * | t_simTrackP |
| std::vector< double > * | t_trackChiSq |
| std::vector< double > * | t_trackChiSqAll |
| std::vector< double > * | t_trackDxy |
| std::vector< double > * | t_trackDxyAll |
| std::vector< double > * | t_trackDxyBS |
| std::vector< double > * | t_trackDxyPV |
| std::vector< double > * | t_trackDxyPVAll |
| std::vector< double > * | t_trackDz |
| std::vector< double > * | t_trackDzAll |
| std::vector< double > * | t_trackDzBS |
| std::vector< double > * | t_trackDzPV |
| std::vector< double > * | t_trackDzPVAll |
| std::vector< double > * | t_trackEcalEta |
| std::vector< double > * | t_trackEcalPhi |
| std::vector< double > * | t_trackEta |
| std::vector< double > * | t_trackEtaAll |
| std::vector< double > * | t_trackHcalEta |
| std::vector< double > * | t_trackHcalPhi |
| std::vector< int > * | t_trackHitInMeasTEC |
| std::vector< int > * | t_trackHitInMeasTIB |
| std::vector< int > * | t_trackHitInMeasTID |
| std::vector< int > * | t_trackHitInMeasTOB |
| std::vector< int > * | t_trackHitInMissTEC |
| std::vector< int > * | t_trackHitInMissTIB |
| std::vector< int > * | t_trackHitInMissTIBTID |
| std::vector< int > * | t_trackHitInMissTID |
| std::vector< int > * | t_trackHitInMissTOB |
| std::vector< int > * | t_trackHitOutMeasTEC |
| std::vector< int > * | t_trackHitOutMeasTIB |
| std::vector< int > * | t_trackHitOutMeasTID |
| std::vector< int > * | t_trackHitOutMeasTOB |
| std::vector< int > * | t_trackHitOutMissTEC |
| std::vector< int > * | t_trackHitOutMissTIB |
| std::vector< int > * | t_trackHitOutMissTID |
| std::vector< int > * | t_trackHitOutMissTOB |
| std::vector< int > * | t_trackHitOutMissTOBTEC |
| std::vector< int > * | t_trackHitsTEC |
| std::vector< int > * | t_trackHitsTOB |
| std::vector< double > * | t_trackL |
| std::vector< int > * | t_trackNOuterHits |
| std::vector< double > * | t_trackOutPosOutHitDr |
| std::vector< double > * | t_trackP |
| std::vector< double > * | t_trackPAll |
| std::vector< double > * | t_trackPdgIdAll |
| std::vector< double > * | t_trackPhi |
| std::vector< double > * | t_trackPhiAll |
| std::vector< double > * | t_trackPt |
| std::vector< double > * | t_trackPtAll |
| std::vector< int > * | t_trackPVIdx |
| std::vector< double > * | t_trkEcalEne |
| std::vector< double > * | t_trkHcalEne |
| double | tMaxE_ |
| double | tMaxH_ |
| double | tMinE_ |
| double | tMinH_ |
| TTree * | tree |
| std::map< std::string, bool > ::iterator | trig_iter |
Static Private Attributes | |
| static const size_t | NEtaBins = 3 |
| static const size_t | nL1BitsMax = 128 |
| static const size_t | NPBins = 15 |
Description: <one line="" class="" summary>="">
Implementation: <Notes on="" implementation>="">
Definition at line 111 of file IsolatedTracksNxN.h.
| IsolatedTracksNxN::IsolatedTracksNxN | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 30 of file IsolatedTracksNxN.cc.
References gather_cfg::cout, debugL1Info_, debugTrks_, doMC, funct::false, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), initL1, JetExtender_, JetSrc_, L1extraCenJetSource_, L1extraFwdJetSource_, L1extraIsoEmSource_, L1extraMuonSource_, L1extraNonIsoEmSource_, L1extraTauJetSource_, L1GTObjectMapRcdSource_, L1GTReadoutRcdSource_, maxTrackEta_, minTrackP_, myverbose_, printTrkHitPattern_, pvTracksPtMin_, tMaxE_, tMaxH_, tMinE_, and tMinH_.
{
//now do what ever initialization is needed
doMC = iConfig.getUntrackedParameter<bool> ("DoMC", false);
myverbose_ = iConfig.getUntrackedParameter<int> ("Verbosity", 5 );
pvTracksPtMin_ = iConfig.getUntrackedParameter<double>("PVTracksPtMin", 1.0 );
debugTrks_ = iConfig.getUntrackedParameter<int> ("DebugTracks" );
printTrkHitPattern_ = iConfig.getUntrackedParameter<bool> ("PrintTrkHitPattern" );
minTrackP_ = iConfig.getUntrackedParameter<double>("minTrackP", 1.0 );
maxTrackEta_ = iConfig.getUntrackedParameter<double>("maxTrackEta", 5.0 );
debugL1Info_ = iConfig.getUntrackedParameter<bool> ("DebugL1Info", false );
L1extraTauJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraTauJetSource" );
L1extraCenJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraCenJetSource" );
L1extraFwdJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraFwdJetSource" );
L1extraMuonSource_ = iConfig.getParameter<edm::InputTag> ("L1extraMuonSource" );
L1extraIsoEmSource_ = iConfig.getParameter<edm::InputTag> ("L1extraIsoEmSource" );
L1extraNonIsoEmSource_ = iConfig.getParameter<edm::InputTag> ("L1extraNonIsoEmSource" );
L1GTReadoutRcdSource_ = iConfig.getParameter<edm::InputTag> ("L1GTReadoutRcdSource" );
L1GTObjectMapRcdSource_= iConfig.getParameter<edm::InputTag> ("L1GTObjectMapRcdSource");
JetSrc_ = iConfig.getParameter<edm::InputTag> ("JetSource");
JetExtender_ = iConfig.getParameter<edm::InputTag> ("JetExtender");
tMinE_ = iConfig.getUntrackedParameter<double>("TimeMinCutECAL", -500.);
tMaxE_ = iConfig.getUntrackedParameter<double>("TimeMaxCutECAL", 500.);
tMinH_ = iConfig.getUntrackedParameter<double>("TimeMinCutHCAL", -500.);
tMaxH_ = iConfig.getUntrackedParameter<double>("TimeMaxCutHCAL", 500.);
//edm::ParameterSet parameters = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
//parameters_.loadParameters( parameters );
//trackAssociator_ = new TrackDetectorAssociator();
//trackAssociator_->useDefaultPropagator();
if(myverbose_>=0) {
std::cout <<"Parameters read from config file \n"
<<" doMC " << doMC
<<"\t myverbose_ " << myverbose_
<<"\t minTrackP_ " << minTrackP_
<< "\t maxTrackEta_ " << maxTrackEta_
<< "\t tMinE_ " << tMinE_
<< "\t tMaxE_ " << tMaxE_
<< "\t tMinH_ " << tMinH_
<< "\t tMaxH_ " << tMaxH_ << "\n"
<< std::endl;
}
initL1 = false;
}
| IsolatedTracksNxN::~IsolatedTracksNxN | ( | ) |
Definition at line 78 of file IsolatedTracksNxN.cc.
{
//delete trackAssociator_;
}
| void IsolatedTracksNxN::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 82 of file IsolatedTracksNxN.cc.
References abs, accept(), ExpressReco_HICollisions_FallBack::algo, algoBitToName, bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationEcal(), spr::chargeIsolationHcal(), clearTreeVectors(), gather_cfg::cout, debugL1Info_, debugTrks_, DeltaR(), diffTreeTool::diff, doMC, reco::TrackBase::dxy(), reco::TrackBase::dz(), EcalBarrel, spr::eCaloSimInfo(), spr::caloSimInfo::eChargedHad, spr::eECALmatrix(), spr::eECALSimInfo(), spr::caloSimInfo::eGamma, spr::eHCALmatrix(), spr::eHCALmatrixCell(), spr::eHCALSimInfo(), spr::caloSimInfo::eMatched, spr::caloSimInfo::eNeutralHad, spr::caloSimInfo::eRest, eta(), reco::TrackBase::eta(), spr::caloSimInfo::eTotal, edm::EventID::event(), genPartEtaBins, genPartPBins, edm::EventSetup::get(), edm::Event::getByLabel(), h_L1AlgoNames, h_maxNearP15x15, h_maxNearP21x21, h_maxNearP25x25, h_maxNearP31x31, h_nTracks, h_PVTracksWt, h_recEta_0, h_recEta_1, h_recEta_2, h_recEtaP_0, h_recEtaP_1, h_recEtaP_2, h_recEtaPt_0, h_recEtaPt_1, h_recEtaPt_2, h_recP_0, h_recP_1, h_recP_2, h_recPhi_0, h_recPhi_1, h_recPhi_2, h_recPt_0, h_recPt_1, h_recPt_2, reco::TrackBase::hitPattern(), edm::EventBase::id(), initL1, edm::HandleBase::isValid(), JetExtender_, analyzePatCleaning_cfg::jets, JetSrc_, L1extraCenJetSource_, L1extraFwdJetSource_, L1extraIsoEmSource_, L1extraMuonSource_, L1extraNonIsoEmSource_, L1extraTauJetSource_, L1GTReadoutRcdSource_, l1TriggerMap, edm::EventBase::luminosityBlock(), spr::matchedSimTrack(), maxTrackEta_, minTrackP_, reco::TrackBase::momentum(), myverbose_, NEtaBins, nEventProc, reco::TrackBase::normalizedChi2(), NPBins, reco::Track::outerPosition(), reco::TrackBase::p(), p1, spr::caloSimInfo::pdgMatched, phi, edm::ESHandle< T >::product(), edm::Handle< T >::product(), spr::propagateCALO(), spr::propagateTrackerEnd(), reco::TrackBase::pt(), ExpressReco_HICollisions_FallBack::pt, pvTracksPtMin_, reco::TrackBase::quality(), reco::TrackBase::qualityByName(), edm::EventID::run(), reco::HitPattern::stripTECLayersWithMeasurement(), reco::HitPattern::stripTECLayersWithoutMeasurement(), reco::HitPattern::stripTIBLayersWithMeasurement(), reco::HitPattern::stripTIBLayersWithoutMeasurement(), reco::HitPattern::stripTIDLayersWithMeasurement(), reco::HitPattern::stripTIDLayersWithoutMeasurement(), reco::HitPattern::stripTOBLayersWithMeasurement(), reco::HitPattern::stripTOBLayersWithoutMeasurement(), t_Bunch, t_e11x11, t_e11x11_10Sig, t_e11x11_15Sig, t_e11x11_20Sig, t_e11x11_25Sig, t_e11x11_30Sig, t_e15x15, t_e15x15_10Sig, t_e15x15_15Sig, t_e15x15_20Sig, t_e15x15_25Sig, t_e15x15_30Sig, t_e7x7, t_e7x7_10Sig, t_e7x7_15Sig, t_e7x7_20Sig, t_e7x7_25Sig, t_e7x7_30Sig, t_e9x9, t_e9x9_10Sig, t_e9x9_15Sig, t_e9x9_20Sig, t_e9x9_25Sig, t_e9x9_30Sig, t_ecalSpike11x11, t_esim11x11, t_esim11x11CharHad, t_esim11x11Matched, t_esim11x11NeutHad, t_esim11x11Photon, t_esim11x11Rest, t_esim15x15, t_esim15x15CharHad, t_esim15x15Matched, t_esim15x15NeutHad, t_esim15x15Photon, t_esim15x15Rest, t_esim7x7, t_esim7x7CharHad, t_esim7x7Matched, t_esim7x7NeutHad, t_esim7x7Photon, t_esim7x7Rest, t_esim9x9, t_esim9x9CharHad, t_esim9x9Matched, t_esim9x9NeutHad, t_esim9x9Photon, t_esim9x9Rest, t_esimPdgId, t_EvtNo, t_h3x3, t_h3x3Sig, t_h5x5, t_h5x5Sig, t_h7x7, t_h7x7Sig, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NeutHad, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NeutHad, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim7x7, t_hsim7x7CharHad, t_hsim7x7Matched, t_hsim7x7NeutHad, t_hsim7x7Photon, t_hsim7x7Rest, t_infoHcal, t_jetEta, t_jetPhi, t_jetPt, t_L1CenJetEta, t_L1CenJetPhi, t_L1CenJetPt, t_L1Decision, t_L1FwdJetEta, t_L1FwdJetPhi, t_L1FwdJetPt, t_L1IsoEMEta, t_L1IsoEMPhi, t_L1IsoEMPt, t_L1MuonEta, t_L1MuonPhi, t_L1MuonPt, t_L1NonIsoEMEta, t_L1NonIsoEMPhi, t_L1NonIsoEMPt, t_L1TauJetEta, t_L1TauJetPhi, t_L1TauJetPt, t_Lumi, t_maxNearHcalP3x3, t_maxNearHcalP5x5, t_maxNearHcalP7x7, t_maxNearP21x21, t_maxNearP31x31, t_NLayersCrossed, t_nTracks, t_nTrksJetCalo, t_nTrksJetVtx, t_PVisValid, t_PVndof, t_PVNTracks, t_PVNTracksHP, t_PVNTracksHPWt, t_PVNTracksWt, t_PVTracksSumPt, t_PVTracksSumPtHP, t_PVTracksSumPtHPWt, t_PVTracksSumPtWt, t_PVx, t_PVy, t_PVz, t_RunNo, t_simTrackP, t_trackChiSq, t_trackChiSqAll, t_trackDxy, t_trackDxyAll, t_trackDxyBS, t_trackDxyPV, t_trackDxyPVAll, t_trackDz, t_trackDzAll, t_trackDzBS, t_trackDzPV, t_trackDzPVAll, t_trackEcalEta, t_trackEcalPhi, t_trackEta, t_trackEtaAll, t_trackHcalEta, t_trackHcalPhi, t_trackHitInMeasTEC, t_trackHitInMeasTIB, t_trackHitInMeasTID, t_trackHitInMeasTOB, t_trackHitInMissTEC, t_trackHitInMissTIB, t_trackHitInMissTIBTID, t_trackHitInMissTID, t_trackHitInMissTOB, t_trackHitOutMeasTEC, t_trackHitOutMeasTIB, t_trackHitOutMeasTID, t_trackHitOutMeasTOB, t_trackHitOutMissTEC, t_trackHitOutMissTIB, t_trackHitOutMissTID, t_trackHitOutMissTOB, t_trackHitOutMissTOBTEC, t_trackHitsTEC, t_trackHitsTOB, t_trackL, t_trackNOuterHits, t_trackOutPosOutHitDr, t_trackP, t_trackPAll, t_trackPdgIdAll, t_trackPhi, t_trackPhiAll, t_trackPt, t_trackPtAll, t_trackPVIdx, t_trkEcalEne, t_trkHcalEne, tMaxE_, tMaxH_, tMinE_, tMinH_, reco::TrackBase::trackerExpectedHitsInner(), reco::TrackBase::trackerExpectedHitsOuter(), reco::HitPattern::trackerLayersWithMeasurement(), tree, trig_iter, CommonMethods::weight(), ExpressReco_HICollisions_FallBack::x, ExpressReco_HICollisions_FallBack::y, and z.
{
// std::cout<< tMinE_ << " " << tMaxE_ << " HCAL " << tMinH_ << " " << tMaxH_ << std::endl;
edm::ESHandle<MagneticField> bFieldH;
iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
bField = bFieldH.product();
clearTreeVectors();
t_RunNo = iEvent.id().run();
t_EvtNo = iEvent.id().event();
t_Lumi = iEvent.luminosityBlock();
t_Bunch = iEvent.bunchCrossing();
nEventProc++;
edm::Handle<reco::TrackCollection> trkCollection;
iEvent.getByLabel("generalTracks", trkCollection);
reco::TrackCollection::const_iterator trkItr;
if(debugTrks_>1){
std::cout << "Track Collection: " << std::endl;
std::cout << "Number of Tracks " << trkCollection->size() << std::endl;
}
std::string theTrackQuality = "highPurity";
reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality);
//===================== save L1 Trigger information =======================
// get L1TriggerReadout records
edm::Handle<L1GlobalTriggerReadoutRecord> gtRecord;
iEvent.getByLabel(L1GTReadoutRcdSource_, gtRecord);
// sanity check on L1 Trigger Records
if (!gtRecord.isValid()) {
std::cout << "\nL1GlobalTriggerReadoutRecord with \n \nnot found"
"\n --> returning false by default!\n" << std::endl;
}
edm::ESHandle<L1GtTriggerMenu> gtOMRec;
iSetup.get<L1GtTriggerMenuRcd>().get(gtOMRec) ;
// L1 decision word
const DecisionWord dWord = gtRecord->decisionWord();
unsigned int numberTriggerBits= dWord.size();
// just print the L1Bit number and AlgoName in first event
/*
if ( !initL1){
initL1=true;
std::cout << "\n Number of Trigger bits " << numberTriggerBits << "\n\n";
std::cout << "\tBit \t L1 Algorithm " << std::endl;
// get ObjectMaps from ObjectMapRecord
const std::vector<L1GlobalTriggerObjectMap>& objMapVec = gtOMRec->gtObjectMap();
for (std::vector<L1GlobalTriggerObjectMap>::const_iterator itMap = objMapVec.begin();
itMap != objMapVec.end(); ++itMap) {
// Get trigger bits
int itrig = (*itMap).algoBitNumber();
// Get trigger names
algoBitToName[itrig] = (*itMap).algoName();
std::cout << "\t" << itrig << "\t" << algoBitToName[itrig] << std::endl;
// store the algoNames as bin labels of a histogram
h_L1AlgoNames->GetXaxis()->SetBinLabel(itrig+1, algoBitToName[itrig].c_str());
} // end of for loop
} // end of initL1
*/
if ( !initL1){
initL1=true;
// std::cout << "\n Number of Trigger bits " << numberTriggerBits << "\n\n";
//std::cout << "\tBit \t L1 Algorithm " << std::endl;
int itrig=0;
for (CItAlgo algo = (*gtOMRec.product()).gtAlgorithmMap().begin(); algo!=(*gtOMRec.product()).gtAlgorithmMap().end(); ++algo, itrig++) {
// std::cout << "Name: " << (algo->second).algoName() << " Alias: " << (algo->second).algoAlias() << std::endl;
algoBitToName[itrig] = (algo->second).algoName() ;
}
} // end of initL1
for (unsigned int iBit = 0; iBit < numberTriggerBits; ++iBit) {
bool accept = dWord[iBit];
//if (accept) h_L1Results2->Fill(float(iBit));
// fill the trigger map
typedef std::map<std::string,bool>::value_type valType;
trig_iter=l1TriggerMap.find(algoBitToName[iBit]);
if (trig_iter==l1TriggerMap.end())
l1TriggerMap.insert(valType(algoBitToName[iBit],accept));
else
trig_iter->second=accept;
}
// save L1 decision for each event
for (unsigned int iBit = 0; iBit < numberTriggerBits; ++iBit) {
bool accept = dWord[iBit];
t_L1Decision->push_back(accept);
// fill the trigger map
if(debugL1Info_) std::cout << "Bit " << iBit << " " << algoBitToName[iBit] << " " << accept << std::endl;
if(accept) h_L1AlgoNames->Fill(iBit);
}
//===================
edm::Handle<reco::VertexCollection> recVtxs;
iEvent.getByLabel("offlinePrimaryVertices",recVtxs);
std::vector<reco::Track> svTracks;
math::XYZPoint leadPV(0,0,0);
double sumPtMax = -1.0;
for(unsigned int ind=0; ind < recVtxs->size(); ind++) {
if (!((*recVtxs)[ind].isFake())) {
double vtxTrkSumPt=0.0, vtxTrkSumPtWt=0.0;
int vtxTrkNWt =0;
double vtxTrkSumPtHP=0.0, vtxTrkSumPtHPWt=0.0;
int vtxTrkNHP =0, vtxTrkNHPWt =0;
reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
for (vtxTrack = (*recVtxs)[ind].tracks_begin(); vtxTrack!=(*recVtxs)[ind].tracks_end(); vtxTrack++) {
bool trkQuality = (*vtxTrack)->quality(trackQuality_);
if((*vtxTrack)->pt()<pvTracksPtMin_) continue;
vtxTrkSumPt += (*vtxTrack)->pt();
if( trkQuality ) {
vtxTrkSumPtHP += (*vtxTrack)->pt();
vtxTrkNHP++;
}
double weight = (*recVtxs)[ind].trackWeight(*vtxTrack);
h_PVTracksWt ->Fill(weight);
if( weight>0.5) {
vtxTrkSumPtWt += (*vtxTrack)->pt();
vtxTrkNWt++;
if( trkQuality ) {
vtxTrkSumPtHPWt += (*vtxTrack)->pt();
vtxTrkNHPWt++;
}
}
}
if(vtxTrkSumPt>sumPtMax) {
sumPtMax = vtxTrkSumPt;
leadPV = math::XYZPoint( (*recVtxs)[ind].x(),(*recVtxs)[ind].y(), (*recVtxs)[ind].z() );
}
t_PVx ->push_back( (*recVtxs)[ind].x() );
t_PVy ->push_back( (*recVtxs)[ind].y() );
t_PVz ->push_back( (*recVtxs)[ind].z() );
t_PVisValid ->push_back( (*recVtxs)[ind].isValid() );
t_PVNTracks ->push_back( (*recVtxs)[ind].tracksSize() );
t_PVndof ->push_back( (*recVtxs)[ind].ndof() );
t_PVNTracksWt ->push_back( vtxTrkNWt );
t_PVTracksSumPt ->push_back( vtxTrkSumPt );
t_PVTracksSumPtWt->push_back( vtxTrkSumPtWt );
t_PVNTracksHP ->push_back( vtxTrkNHP );
t_PVNTracksHPWt ->push_back( vtxTrkNHPWt );
t_PVTracksSumPtHP ->push_back( vtxTrkSumPtHP );
t_PVTracksSumPtHPWt->push_back( vtxTrkSumPtHPWt );
if(myverbose_==4) {
std::cout<<"PV "<<ind<<" isValid "<<(*recVtxs)[ind].isValid()<<" isFake "<<(*recVtxs)[ind].isFake()
<<" hasRefittedTracks() "<<ind<<" "<<(*recVtxs)[ind].hasRefittedTracks()
<<" refittedTrksSize "<<(*recVtxs)[ind].refittedTracks().size()
<<" tracksSize() "<<(*recVtxs)[ind].tracksSize()<<" sumPt "<<vtxTrkSumPt
<<std::endl;
}
}
}
// Get the beamspot
edm::Handle<reco::BeamSpot> beamSpotH;
iEvent.getByLabel("offlineBeamSpot", beamSpotH);
math::XYZPoint bspot;
bspot = ( beamSpotH.isValid() ) ? beamSpotH->position() : math::XYZPoint(0, 0, 0);
//===================
// L1Taus
edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
iEvent.getByLabel(L1extraTauJetSource_,l1TauHandle);
l1extra::L1JetParticleCollection::const_iterator itr;
for(itr = l1TauHandle->begin(); itr != l1TauHandle->end(); ++itr ) {
t_L1TauJetPt ->push_back( itr->pt() );
t_L1TauJetEta ->push_back( itr->eta() );
t_L1TauJetPhi ->push_back( itr->phi() );
if(debugL1Info_) {
std::cout << "tauJ p/pt " << itr->momentum() << " " << itr->pt()
<< " eta/phi " << itr->eta() << " " << itr->phi()
<< std::endl;
}
}
// L1 Central Jets
edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
iEvent.getByLabel(L1extraCenJetSource_,l1CenJetHandle);
for( itr = l1CenJetHandle->begin(); itr != l1CenJetHandle->end(); ++itr ) {
t_L1CenJetPt ->push_back( itr->pt() );
t_L1CenJetEta ->push_back( itr->eta() );
t_L1CenJetPhi ->push_back( itr->phi() );
if(debugL1Info_) {
std::cout << "cenJ p/pt " << itr->momentum() << " " << itr->pt()
<< " eta/phi " << itr->eta() << " " << itr->phi()
<< std::endl;
}
}
// L1 Forward Jets
edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
iEvent.getByLabel(L1extraFwdJetSource_,l1FwdJetHandle);
for( itr = l1FwdJetHandle->begin(); itr != l1FwdJetHandle->end(); ++itr ) {
t_L1FwdJetPt ->push_back( itr->pt() );
t_L1FwdJetEta ->push_back( itr->eta() );
t_L1FwdJetPhi ->push_back( itr->phi() );
if(debugL1Info_) {
std::cout << "fwdJ p/pt " << itr->momentum() << " " << itr->pt()
<< " eta/phi " << itr->eta() << " " << itr->phi()
<< std::endl;
}
}
// L1 Isolated EM onjects
l1extra::L1EmParticleCollection::const_iterator itrEm;
edm::Handle<l1extra::L1EmParticleCollection> l1IsoEmHandle ;
iEvent.getByLabel(L1extraIsoEmSource_, l1IsoEmHandle);
for( itrEm = l1IsoEmHandle->begin(); itrEm != l1IsoEmHandle->end(); ++itrEm ) {
t_L1IsoEMPt ->push_back( itrEm->pt() );
t_L1IsoEMEta ->push_back( itrEm->eta() );
t_L1IsoEMPhi ->push_back( itrEm->phi() );
if(debugL1Info_) {
std::cout << "isoEm p/pt " << itrEm->momentum() << " " << itrEm->pt()
<< " eta/phi " << itrEm->eta() << " " << itrEm->phi()
<< std::endl;
}
}
// L1 Non-Isolated EM onjects
edm::Handle<l1extra::L1EmParticleCollection> l1NonIsoEmHandle ;
iEvent.getByLabel(L1extraNonIsoEmSource_, l1NonIsoEmHandle);
for( itrEm = l1NonIsoEmHandle->begin(); itrEm != l1NonIsoEmHandle->end(); ++itrEm ) {
t_L1NonIsoEMPt ->push_back( itrEm->pt() );
t_L1NonIsoEMEta ->push_back( itrEm->eta() );
t_L1NonIsoEMPhi ->push_back( itrEm->phi() );
if(debugL1Info_) {
std::cout << "nonIsoEm p/pt " << itrEm->momentum() << " " << itrEm->pt()
<< " eta/phi " << itrEm->eta() << " " << itrEm->phi()
<< std::endl;
}
}
// L1 Muons
l1extra::L1MuonParticleCollection::const_iterator itrMu;
edm::Handle<l1extra::L1MuonParticleCollection> l1MuHandle ;
iEvent.getByLabel(L1extraMuonSource_, l1MuHandle);
for( itrMu = l1MuHandle->begin(); itrMu != l1MuHandle->end(); ++itrMu ) {
t_L1MuonPt ->push_back( itrMu->pt() );
t_L1MuonEta ->push_back( itrMu->eta() );
t_L1MuonPhi ->push_back( itrMu->phi() );
if(debugL1Info_) {
std::cout << "l1muon p/pt " << itrMu->momentum() << " " << itrMu->pt()
<< " eta/phi " << itrMu->eta() << " " << itrMu->phi()
<< std::endl;
}
}
//=====================================================================
edm::Handle<reco::CaloJetCollection> jets;
iEvent.getByLabel(JetSrc_,jets);
edm::Handle<reco::JetExtendedAssociation::Container> jetExtender;
iEvent.getByLabel(JetExtender_,jetExtender);
for(unsigned int ijet=0;ijet<(*jets).size();ijet++) {
t_jetPt ->push_back( (*jets)[ijet].pt() );
t_jetEta ->push_back( (*jets)[ijet].eta() );
t_jetPhi ->push_back( (*jets)[ijet].phi() );
//t_nTrksJetVtx ->push_back( reco::JetExtendedAssociation::tracksAtVertexNumber(*jetExtender,(*jets)[ijet]) );
//t_nTrksJetCalo ->push_back( reco::JetExtendedAssociation::tracksAtCaloNumber (*jetExtender,(*jets)[ijet]) );
t_nTrksJetVtx ->push_back( -1.0);
t_nTrksJetCalo ->push_back( -1.0);
}
//=====================================================================
// get handles to calogeometry and calotopology
edm::ESHandle<CaloGeometry> pG;
iSetup.get<CaloGeometryRecord>().get(pG);
const CaloGeometry* geo = pG.product();
edm::ESHandle<CaloTopology> theCaloTopology;
iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
const CaloTopology *caloTopology = theCaloTopology.product();
edm::ESHandle<HcalTopology> htopo;
iSetup.get<IdealGeometryRecord>().get(htopo);
const HcalTopology* theHBHETopology = htopo.product();
edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
iEvent.getByLabel("ecalRecHit","EcalRecHitsEB",barrelRecHitsHandle);
iEvent.getByLabel("ecalRecHit","EcalRecHitsEE",endcapRecHitsHandle);
// Retrieve the good/bad ECAL channels from the DB
edm::ESHandle<EcalChannelStatus> ecalChStatus;
iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();
// Retrieve trigger tower map
//const edm::ESHandle<EcalTrigTowerConstituentsMap> hTtmap;
edm::ESHandle<EcalTrigTowerConstituentsMap> hTtmap;
iSetup.get<IdealGeometryRecord>().get(hTtmap);
const EcalTrigTowerConstituentsMap& ttMap = *hTtmap;
edm::Handle<HBHERecHitCollection> hbhe;
iEvent.getByLabel("hbhereco",hbhe);
const HBHERecHitCollection Hithbhe = *(hbhe.product());
//get Handles to SimTracks and SimHits
edm::Handle<edm::SimTrackContainer> SimTk;
if (doMC) iEvent.getByLabel("g4SimHits",SimTk);
edm::SimTrackContainer::const_iterator simTrkItr;
edm::Handle<edm::SimVertexContainer> SimVtx;
if (doMC) iEvent.getByLabel("g4SimHits",SimVtx);
edm::SimVertexContainer::const_iterator vtxItr = SimVtx->begin();
//get Handles to PCaloHitContainers of eb/ee/hbhe
edm::Handle<edm::PCaloHitContainer> pcaloeb;
if (doMC) iEvent.getByLabel("g4SimHits", "EcalHitsEB", pcaloeb);
edm::Handle<edm::PCaloHitContainer> pcaloee;
if (doMC) iEvent.getByLabel("g4SimHits", "EcalHitsEE", pcaloee);
edm::Handle<edm::PCaloHitContainer> pcalohh;
if (doMC) iEvent.getByLabel("g4SimHits", "HcalHits", pcalohh);
//associates tracker rechits/simhits to a track
TrackerHitAssociator* associate=0;
if (doMC) associate = new TrackerHitAssociator(iEvent);
std::vector<int> ifGood(trkCollection->size(), 1);
h_nTracks->Fill(trkCollection->size());
int nTracks = 0;
t_nTracks = 0;
t_nTracks = trkCollection->size();
//decide the goodness of each track for this track collection
for( trkItr = trkCollection->begin(),nTracks=0; trkItr != trkCollection->end(); ++trkItr, nTracks++){
const reco::Track* pTrack = &(*trkItr);
bool trkQuality = pTrack->quality(trackQuality_);
if( !trkQuality ) ifGood[nTracks]=0;
}
// get the list of DetIds closest to the impact point of track on surface calorimeters
std::vector<spr::propagatedTrackID> trkCaloDets;
spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDets, false);
std::vector<spr::propagatedTrackID>::const_iterator trkDetItr;
if(myverbose_>2) {
for(trkDetItr = trkCaloDets.begin(); trkDetItr != trkCaloDets.end(); trkDetItr++){
std::cout<<trkDetItr->trkItr->p()<<" "<<trkDetItr->trkItr->eta()<<" "<<trkDetItr->okECAL<<" ";
if(trkDetItr->detIdECAL.subdetId() == EcalBarrel) std::cout << (EBDetId)trkDetItr->detIdECAL <<" ";
else std::cout << (EEDetId)trkDetItr->detIdECAL <<" ";
std::cout<<trkDetItr->okHCAL<<" ";
if(trkDetItr->okHCAL) std::cout<<(HcalDetId)trkDetItr->detIdHCAL;
std::cout << std::endl;
}
}
for(trkDetItr = trkCaloDets.begin(),nTracks=0; trkDetItr != trkCaloDets.end(); trkDetItr++,nTracks++){
const reco::Track* pTrack = &(*(trkDetItr->trkItr));
const reco::HitPattern& hitp = pTrack->hitPattern();
int nLayersCrossed = hitp.trackerLayersWithMeasurement() ;
int nOuterHits = hitp.stripTOBLayersWithMeasurement()+hitp.stripTECLayersWithMeasurement() ;
const reco::HitPattern& hitpIn = pTrack->trackerExpectedHitsInner();
const reco::HitPattern& hitpOut = pTrack->trackerExpectedHitsOuter();
double eta1 = pTrack->momentum().eta();
double phi1 = pTrack->momentum().phi();
double etaEcal1 = trkDetItr->etaECAL;
double phiEcal1 = trkDetItr->phiECAL;
double etaHcal1 = trkDetItr->etaHCAL;
double phiHcal1 = trkDetItr->phiHCAL;
double pt1 = pTrack->pt();
double p1 = pTrack->p();
double dxy1 = pTrack->dxy();
double dz1 = pTrack->dz();
double dxybs1 = beamSpotH.isValid() ? pTrack->dxy(bspot) : pTrack->dxy();
double dzbs1 = beamSpotH.isValid() ? pTrack->dz(bspot) : pTrack->dz();
double dxypv1 = (recVtxs->size()>0 && !((*recVtxs)[0].isFake())) ? pTrack->dxy(leadPV) : pTrack->dxy();
double dzpv1 = (recVtxs->size()>0 && !((*recVtxs)[0].isFake())) ? pTrack->dz(leadPV) : pTrack->dz();
double chisq1 = pTrack->normalizedChi2();
h_recEtaPt_0->Fill(eta1, pt1);
h_recEtaP_0 ->Fill(eta1, p1);
h_recPt_0 ->Fill(pt1);
h_recP_0 ->Fill(p1);
h_recEta_0 ->Fill(eta1);
h_recPhi_0 ->Fill(phi1);
if(ifGood[nTracks] && nLayersCrossed>7 ) {
h_recEtaPt_1->Fill(eta1, pt1);
h_recEtaP_1 ->Fill(eta1, p1);
h_recPt_1 ->Fill(pt1);
h_recP_1 ->Fill(p1);
h_recEta_1 ->Fill(eta1);
h_recPhi_1 ->Fill(phi1);
}
if( ! ifGood[nTracks] ) continue;
if( pt1>2.0 && nLayersCrossed>7) {
t_trackPAll ->push_back( p1 );
t_trackEtaAll ->push_back( eta1 );
t_trackPhiAll ->push_back( phi1 );
t_trackPtAll ->push_back( pt1 );
t_trackDxyAll ->push_back( dxy1 );
t_trackDzAll ->push_back( dz1 );
t_trackDxyPVAll ->push_back( dxypv1 );
t_trackDzPVAll ->push_back( dzpv1 );
t_trackChiSqAll ->push_back( chisq1 );
}
if (doMC) {
edm::SimTrackContainer::const_iterator matchedSimTrkAll = spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
if( matchedSimTrkAll != SimTk->end()) t_trackPdgIdAll->push_back( matchedSimTrkAll->type() );
}
if( pt1>minTrackP_ && std::abs(eta1)<maxTrackEta_ && trkDetItr->okECAL) {
double maxNearP31x31=999.0, maxNearP25x25=999.0, maxNearP21x21=999.0, maxNearP15x15=999.0;
//double maxNearP13x13=999.0, maxNearP11x11=999.0, maxNearP9x9 =999.0, maxNearP7x7 =999.0;
maxNearP31x31 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 15,15);
maxNearP25x25 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 12,12);
maxNearP21x21 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 10,10);
maxNearP15x15 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 7, 7);
//maxNearP13x13 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 6, 6);
//maxNearP11x11 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 5, 5);
//maxNearP9x9 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 4, 4);
//maxNearP7x7 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 3, 3);
int iTrkEtaBin=-1, iTrkMomBin=-1;
for(unsigned int ieta=0; ieta<NEtaBins; ieta++) {
if(std::abs(eta1)>genPartEtaBins[ieta] && std::abs(eta1)<genPartEtaBins[ieta+1] ) iTrkEtaBin = ieta;
}
for(unsigned int ipt=0; ipt<NPBins; ipt++) {
if( p1>genPartPBins[ipt] && p1<genPartPBins[ipt+1] ) iTrkMomBin = ipt;
}
if( iTrkMomBin>=0 && iTrkEtaBin>=0 ) {
h_maxNearP31x31[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP31x31 );
h_maxNearP25x25[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP25x25 );
h_maxNearP21x21[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP21x21 );
h_maxNearP15x15[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP15x15 );
}
if( maxNearP31x31<0.0 && nLayersCrossed>7 && nOuterHits>4) {
h_recEtaPt_2->Fill(eta1, pt1);
h_recEtaP_2 ->Fill(eta1, p1);
h_recPt_2 ->Fill(pt1);
h_recP_2 ->Fill(p1);
h_recEta_2 ->Fill(eta1);
h_recPhi_2 ->Fill(phi1);
}
// if isolated in 7x7 then store the further quantities
//if( maxNearP7x7<0.0) {
if( maxNearP31x31<0.0) {
// get the matching simTrack
double simTrackP = -1;
if (doMC) {
edm::SimTrackContainer::const_iterator matchedSimTrk = spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
if( matchedSimTrk != SimTk->end() )simTrackP = matchedSimTrk->momentum().P();
}
// get ECal Tranverse Profile
//double e3x3=0, e5x5=0, e7x7=0, e9x9=0, e11x11=0, e13x13=0, e15x15=0, e21x21=0, e25x25=0, e31x31=0;
//double e7x7_10Sig=0, e9x9_10Sig=0, e11x11_10Sig=0, e15x15_10Sig=0;
//double e7x7_15Sig=0, e9x9_15Sig=0, e11x11_15Sig=0, e15x15_15Sig=0;
//double e7x7_20Sig=0, e9x9_20Sig=0, e11x11_20Sig=0, e15x15_20Sig=0;
//double e7x7_25Sig=0, e9x9_25Sig=0, e11x11_25Sig=0, e15x15_25Sig=0;
//double e7x7_30Sig=0, e9x9_30Sig=0, e11x11_30Sig=0, e15x15_30Sig=0;
std::pair<double, bool> e3x3P, e5x5P, e7x7P, e9x9P, e11x11P, e13x13P, e15x15P, e21x21P, e25x25P, e31x31P;
std::pair<double, bool> e7x7_10SigP, e9x9_10SigP, e11x11_10SigP, e15x15_10SigP;
std::pair<double, bool> e7x7_15SigP, e9x9_15SigP, e11x11_15SigP, e15x15_15SigP;
std::pair<double, bool> e7x7_20SigP, e9x9_20SigP, e11x11_20SigP, e15x15_20SigP;
std::pair<double, bool> e7x7_25SigP, e9x9_25SigP, e11x11_25SigP, e15x15_25SigP;
std::pair<double, bool> e7x7_30SigP, e9x9_30SigP, e11x11_30SigP, e15x15_30SigP;
spr::caloSimInfo simInfo3x3, simInfo5x5, simInfo7x7, simInfo9x9;
spr::caloSimInfo simInfo11x11, simInfo13x13, simInfo15x15, simInfo21x21, simInfo25x25, simInfo31x31;
double trkEcalEne=0;
const DetId isoCell = trkDetItr->detIdECAL;
/*
e3x3 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,1,1, -100.0, -100.0);
e5x5 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,2,2, -100.0, -100.0);
e7x7 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,3,3, -100.0, -100.0);
e9x9 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,4,4, -100.0, -100.0);
e11x11 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,5,5, -100.0, -100.0);
e13x13 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,6,6, -100.0, -100.0);
e15x15 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,7,7, -100.0, -100.0);
e21x21 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,10,10, -100.0, -100.0);
e25x25 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,12,12, -100.0, -100.0);
e31x31 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,15,15, -100.0, -100.0);
e7x7_10Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,3,3, 0.030, 0.150);
e9x9_10Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,4,4, 0.030, 0.150);
e11x11_10Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,5,5, 0.030, 0.150);
e15x15_10Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,7,7, 0.030, 0.150);
e7x7_15Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,3,3, 0.045, 0.225);
e9x9_15Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,4,4, 0.045, 0.225);
e11x11_15Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,5,5, 0.045, 0.225);
e15x15_15Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,7,7, 0.045, 0.225);
e7x7_20Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,3,3, 0.060, 0.300);
e9x9_20Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,4,4, 0.060, 0.300);
e11x11_20Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,5,5, 0.060, 0.300);
e15x15_20Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,7,7, 0.060, 0.300);
e7x7_25Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,3,3, 0.075, 0.375);
e9x9_25Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,4,4, 0.075, 0.375);
e11x11_25Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,5,5, 0.075, 0.375);
e15x15_25Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,7,7, 0.075, 0.375);
e7x7_30Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,3,3, 0.090, 0.450);
e9x9_30Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,4,4, 0.090, 0.450);
e11x11_30Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,5,5, 0.090, 0.450);
e15x15_30Sig = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, geo, caloTopology,7,7, 0.090, 0.450);
std::cout << "default ecal rechit " << std::endl;
std::cout<<"e3x3 "<<e3x3<<" e9x9 "<<e9x9<<" e15x15 " << e15x15 << " e31x31 "<<e31x31<<std::endl;
std::cout<<"e7x7_10Sig "<<e7x7_10Sig<<" e11x11_10Sig "<<e11x11_10Sig<<" e15x15_10Sig "<<e15x15_10Sig<<std::endl;
std::pair<double, bool> eJunk1 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, 0.060, 0.300, false, true);
std::pair<double, bool> eJunk2 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, 0.20, 0.45, true, true);
*/
//e3x3P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,1,1, -100.0, -100.0);
//e5x5P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,2,2, -100.0, -100.0);
e7x7P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, -100.0, -100.0, tMinE_,tMaxE_);
e9x9P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,4,4, -100.0, -100.0, tMinE_,tMaxE_);
e11x11P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,5,5, -100.0, -100.0, tMinE_,tMaxE_);
//e13x13P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,6,6, -100.0, -100.0);
e15x15P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,7,7, -100.0, -100.0, tMinE_,tMaxE_);
//e21x21P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,10,10, -100.0, -100.0);
//e25x25P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,12,12, -100.0, -100.0);
//e31x31P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,15,15, -100.0, -100.0);
e7x7_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, 0.030, 0.150, tMinE_,tMaxE_);
e9x9_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,4,4, 0.030, 0.150, tMinE_,tMaxE_);
e11x11_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,5,5, 0.030, 0.150, tMinE_,tMaxE_);
e15x15_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,7,7, 0.030, 0.150, tMinE_,tMaxE_);
/*
e7x7_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, 0.045, 0.225);
e9x9_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,4,4, 0.045, 0.225);
e11x11_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,5,5, 0.045, 0.225);
e15x15_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,7,7, 0.045, 0.225);
*/
e7x7_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, ttMap, 3,3, 0.20,0.45, tMinE_,tMaxE_);
e9x9_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, ttMap, 4,4, 0.20,0.45, tMinE_,tMaxE_);
e11x11_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, ttMap, 5,5, 0.20,0.45, tMinE_,tMaxE_);
e15x15_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, ttMap, 7,7, 0.20,0.45, tMinE_,tMaxE_, false);
e7x7_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, 0.060, 0.300, tMinE_,tMaxE_);
e9x9_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,4,4, 0.060, 0.300, tMinE_,tMaxE_);
e11x11_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,5,5, 0.060, 0.300, tMinE_,tMaxE_);
e15x15_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,7,7, 0.060, 0.300, tMinE_,tMaxE_);
e7x7_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, 0.075, 0.375, tMinE_,tMaxE_);
e9x9_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,4,4, 0.075, 0.375, tMinE_,tMaxE_);
e11x11_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,5,5, 0.075, 0.375, tMinE_,tMaxE_);
e15x15_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,7,7, 0.075, 0.375, tMinE_,tMaxE_);
e7x7_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,3,3, 0.090, 0.450, tMinE_,tMaxE_);
e9x9_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,4,4, 0.090, 0.450, tMinE_,tMaxE_);
e11x11_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,5,5, 0.090, 0.450, tMinE_,tMaxE_);
e15x15_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,7,7, 0.090, 0.450, tMinE_,tMaxE_);
if(myverbose_ == 2) {
std::cout << "clean ecal rechit " << std::endl;
std::cout<<"e3x3 "<<e3x3P.first<<" e9x9 "<<e9x9P.first<<" e15x15 " << e15x15P.first << " e31x31 "<<e31x31P.first<<std::endl;
std::cout<<"e7x7_10Sig "<<e7x7_10SigP.first<<" e11x11_10Sig "<<e11x11_10SigP.first<<" e15x15_10Sig "<<e15x15_10SigP.first<<std::endl;
}
if (doMC) {
// check the energy from SimHits
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 1,1, simInfo3x3);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 2,2, simInfo5x5);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 3,3, simInfo7x7);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 4,4, simInfo9x9);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 5,5, simInfo11x11);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 6,6, simInfo13x13);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 7,7, simInfo15x15, 150.0,false);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 10,10, simInfo21x21);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 12,12, simInfo25x25);
spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 15,15, simInfo31x31);
trkEcalEne = spr::eCaloSimInfo(iEvent, geo, pcaloeb,pcaloee, SimTk, SimVtx, pTrack, *associate, 150.0, false);
if(myverbose_ == 1) {
std::cout << "Track momentum " << pt1 << std::endl;
std::cout << "ecal siminfo " << std::endl;
std::cout << "simInfo3x3: " << "eTotal " << simInfo3x3.eTotal << " eMatched " << simInfo3x3.eMatched << " eRest " << simInfo3x3.eRest << " eGamma "<<simInfo3x3.eGamma << " eNeutralHad " << simInfo3x3.eNeutralHad << " eChargedHad " << simInfo3x3.eChargedHad << std::endl;
std::cout << "simInfo5x5: " << "eTotal " << simInfo5x5.eTotal << " eMatched " << simInfo5x5.eMatched << " eRest " << simInfo5x5.eRest << " eGamma "<<simInfo5x5.eGamma << " eNeutralHad " << simInfo5x5.eNeutralHad << " eChargedHad " << simInfo5x5.eChargedHad << std::endl;
std::cout << "simInfo7x7: " << "eTotal " << simInfo7x7.eTotal << " eMatched " << simInfo7x7.eMatched << " eRest " << simInfo7x7.eRest << " eGamma "<<simInfo7x7.eGamma << " eNeutralHad " << simInfo7x7.eNeutralHad << " eChargedHad " << simInfo7x7.eChargedHad << std::endl;
std::cout << "simInfo9x9: " << "eTotal " << simInfo9x9.eTotal << " eMatched " << simInfo9x9.eMatched << " eRest " << simInfo9x9.eRest << " eGamma "<<simInfo9x9.eGamma << " eNeutralHad " << simInfo9x9.eNeutralHad << " eChargedHad " << simInfo9x9.eChargedHad << std::endl;
std::cout << "simInfo11x11: " << "eTotal " << simInfo11x11.eTotal << " eMatched " << simInfo11x11.eMatched << " eRest " << simInfo11x11.eRest << " eGamma "<<simInfo11x11.eGamma << " eNeutralHad " << simInfo11x11.eNeutralHad << " eChargedHad " << simInfo11x11.eChargedHad << std::endl;
std::cout << "simInfo15x15: " << "eTotal " << simInfo15x15.eTotal << " eMatched " << simInfo15x15.eMatched << " eRest " << simInfo15x15.eRest << " eGamma "<<simInfo15x15.eGamma << " eNeutralHad " << simInfo15x15.eNeutralHad << " eChargedHad " << simInfo15x15.eChargedHad << std::endl;
std::cout << "simInfo31x31: " << "eTotal " << simInfo31x31.eTotal << " eMatched " << simInfo31x31.eMatched << " eRest " << simInfo31x31.eRest << " eGamma "<<simInfo31x31.eGamma << " eNeutralHad " << simInfo31x31.eNeutralHad << " eChargedHad " << simInfo31x31.eChargedHad << std::endl;
std::cout << "trkEcalEne" << trkEcalEne << std::endl;
}
}
// ======= Get HCAL information
double hcalScale=1.0;
if( std::abs(pTrack->eta())<1.4 ) {
hcalScale=120.0;
} else {
hcalScale=135.0;
}
double maxNearHcalP3x3=-1, maxNearHcalP5x5=-1, maxNearHcalP7x7=-1;
maxNearHcalP3x3 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 1,1);
maxNearHcalP5x5 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 2,2);
maxNearHcalP7x7 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 3,3);
double h3x3=0, h5x5=0, h7x7=0;
double h3x3Sig=0, h5x5Sig=0, h7x7Sig=0;
double trkHcalEne = 0;
spr::caloSimInfo hsimInfo3x3, hsimInfo5x5, hsimInfo7x7;
if(trkDetItr->okHCAL) {
const DetId ClosestCell(trkDetItr->detIdHCAL);
// bool includeHO=false, bool algoNew=true, bool debug=false
h3x3 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,1,1, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_,tMaxH_);
h5x5 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,2,2, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_,tMaxH_);
h7x7 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,3,3, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_,tMaxH_);
h3x3Sig = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,1,1, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_);
h5x5Sig = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,2,2, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_);
h7x7Sig = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,3,3, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_);
if(myverbose_==2) {
std::cout << "HCAL 3x3 " << h3x3 << " " << h3x3Sig << " 5x5 " << h5x5 << " " << h5x5Sig << " 7x7 " << h7x7 << " " << h7x7Sig << std::endl;
}
if (doMC) {
spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, geo,pcalohh, SimTk, SimVtx, pTrack, *associate, 1,1, hsimInfo3x3);
spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, geo,pcalohh, SimTk, SimVtx, pTrack, *associate, 2,2, hsimInfo5x5);
spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, geo,pcalohh, SimTk, SimVtx, pTrack, *associate, 3,3, hsimInfo7x7, 150.0, false,false);
trkHcalEne = spr::eCaloSimInfo(iEvent, geo,pcalohh, SimTk, SimVtx, pTrack, *associate);
if(myverbose_ == 1) {
std::cout << "Hcal siminfo " << std::endl;
std::cout << "hsimInfo3x3: " << "eTotal " << hsimInfo3x3.eTotal << " eMatched " << hsimInfo3x3.eMatched << " eRest " << hsimInfo3x3.eRest << " eGamma "<<hsimInfo3x3.eGamma << " eNeutralHad " << hsimInfo3x3.eNeutralHad << " eChargedHad " << hsimInfo3x3.eChargedHad << std::endl;
std::cout << "hsimInfo5x5: " << "eTotal " << hsimInfo5x5.eTotal << " eMatched " << hsimInfo5x5.eMatched << " eRest " << hsimInfo5x5.eRest << " eGamma "<<hsimInfo5x5.eGamma << " eNeutralHad " << hsimInfo5x5.eNeutralHad << " eChargedHad " << hsimInfo5x5.eChargedHad << std::endl;
std::cout << "hsimInfo7x7: " << "eTotal " << hsimInfo7x7.eTotal << " eMatched " << hsimInfo7x7.eMatched << " eRest " << hsimInfo7x7.eRest << " eGamma "<<hsimInfo7x7.eGamma << " eNeutralHad " << hsimInfo7x7.eNeutralHad << " eChargedHad " << hsimInfo7x7.eChargedHad << std::endl;
std::cout << "trkHcalEne " << trkHcalEne << std::endl;
}
}
// debug the ecal and hcal matrix
if(myverbose_==4) {
std::cout<<"Run "<<iEvent.id().run()<<" Event "<<iEvent.id().event()<<std::endl;
std::vector<std::pair<DetId,double> > v7x7 = spr::eHCALmatrixCell(theHBHETopology, ClosestCell, hbhe,3,3, false, false);
double sumv=0.0;
for(unsigned int iv=0; iv<v7x7.size(); iv++) {
sumv += v7x7[iv].second;
}
std::cout<<"h7x7 "<<h7x7<<" v7x7 "<<sumv << " in " << v7x7.size() <<std::endl;
for(unsigned int iv=0; iv<v7x7.size(); iv++) {
HcalDetId id = v7x7[iv].first;
std::cout << " Cell " << iv << " 0x" << std::hex << v7x7[iv].first() << std::dec << " " << id << " Energy " << v7x7[iv].second << std::endl;
}
}
}
// ====================================================================================================
// get diff between track outermost hit position and the propagation point at outermost surface of tracker
std::pair<math::XYZPoint,double> point2_TK0 = spr::propagateTrackerEnd( pTrack, bField, false);
math::XYZPoint diff(pTrack->outerPosition().X()-point2_TK0.first.X(),
pTrack->outerPosition().Y()-point2_TK0.first.Y(),
pTrack->outerPosition().Z()-point2_TK0.first.Z() );
double trackOutPosOutHitDr = diff.R();
double trackL = point2_TK0.second;
//std::cout<<" propagted "<<point2_TK0.first<<" "<< point2_TK0.first.eta()<<" "<<point2_TK0.first.phi()<<std::endl;
//std::cout<<" outerPosition() "<< pTrack->outerPosition() << " "<< pTrack->outerPosition().eta()<< " " << pTrack->outerPosition().phi()<< std::endl;
//std::cout<<"diff " << diff << " diffR " <<diff.R()<<" diffR/L "<<diff.R()/point2_TK0.second <<std::endl;
for(unsigned int ind=0;ind<recVtxs->size();ind++) {
if (!((*recVtxs)[ind].isFake())) {
reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
if( DeltaR(eta1,phi1, (*vtxTrack)->eta(),(*vtxTrack)->phi()) < 0.01 ) t_trackPVIdx ->push_back( ind );
else t_trackPVIdx ->push_back( -1 );
}
}
// Fill the tree Branches here
t_trackP ->push_back( p1 );
t_trackPt ->push_back( pt1 );
t_trackEta ->push_back( eta1 );
t_trackPhi ->push_back( phi1 );
t_trackEcalEta ->push_back( etaEcal1 );
t_trackEcalPhi ->push_back( phiEcal1 );
t_trackHcalEta ->push_back( etaHcal1 );
t_trackHcalPhi ->push_back( phiHcal1 );
t_trackDxy ->push_back( dxy1 );
t_trackDz ->push_back( dz1 );
t_trackDxyBS ->push_back( dxybs1 );
t_trackDzBS ->push_back( dzbs1 );
t_trackDxyPV ->push_back( dxypv1 );
t_trackDzPV ->push_back( dzpv1 );
t_trackChiSq ->push_back( chisq1 );
t_trackNOuterHits ->push_back( nOuterHits );
t_NLayersCrossed ->push_back( nLayersCrossed );
t_trackHitsTOB ->push_back( hitp.stripTOBLayersWithMeasurement() );
t_trackHitsTEC ->push_back( hitp.stripTECLayersWithMeasurement() );
t_trackHitInMissTOB ->push_back( hitpIn.stripTOBLayersWithoutMeasurement() );
t_trackHitInMissTEC ->push_back( hitpIn.stripTECLayersWithoutMeasurement() );
t_trackHitInMissTIB ->push_back( hitpIn.stripTIBLayersWithoutMeasurement() );
t_trackHitInMissTID ->push_back( hitpIn.stripTIDLayersWithoutMeasurement() );
t_trackHitInMissTIBTID ->push_back( hitpIn.stripTIBLayersWithoutMeasurement() + hitpIn.stripTIDLayersWithoutMeasurement() );
t_trackHitOutMissTOB ->push_back( hitpOut.stripTOBLayersWithoutMeasurement() );
t_trackHitOutMissTEC ->push_back( hitpOut.stripTECLayersWithoutMeasurement() );
t_trackHitOutMissTIB ->push_back( hitpOut.stripTIBLayersWithoutMeasurement() );
t_trackHitOutMissTID ->push_back( hitpOut.stripTIDLayersWithoutMeasurement() );
t_trackHitOutMissTOBTEC ->push_back( hitpOut.stripTOBLayersWithoutMeasurement() + hitpOut.stripTECLayersWithoutMeasurement() );
t_trackHitInMeasTOB ->push_back( hitpIn.stripTOBLayersWithMeasurement() );
t_trackHitInMeasTEC ->push_back( hitpIn.stripTECLayersWithMeasurement() );
t_trackHitInMeasTIB ->push_back( hitpIn.stripTIBLayersWithMeasurement() );
t_trackHitInMeasTID ->push_back( hitpIn.stripTIDLayersWithMeasurement() );
t_trackHitOutMeasTOB ->push_back( hitpOut.stripTOBLayersWithMeasurement() );
t_trackHitOutMeasTEC ->push_back( hitpOut.stripTECLayersWithMeasurement() );
t_trackHitOutMeasTIB ->push_back( hitpOut.stripTIBLayersWithMeasurement() );
t_trackHitOutMeasTID ->push_back( hitpOut.stripTIDLayersWithMeasurement() );
t_trackOutPosOutHitDr ->push_back( trackOutPosOutHitDr );
t_trackL ->push_back( trackL );
t_maxNearP31x31 ->push_back( maxNearP31x31 );
//t_maxNearP25x25 ->push_back( maxNearP25x25 );
t_maxNearP21x21 ->push_back( maxNearP21x21 );
//t_maxNearP15x15 ->push_back( maxNearP15x15 );
//t_maxNearP13x13 ->push_back( maxNearP13x13 );
//t_maxNearP11x11 ->push_back( maxNearP11x11 );
//t_maxNearP9x9 ->push_back( maxNearP9x9 );
//t_maxNearP7x7 ->push_back( maxNearP7x7 );
t_ecalSpike11x11 ->push_back( e11x11P.second );
//t_e3x3 ->push_back( e3x3P.first );
//t_e5x5 ->push_back( e5x5P.first );
t_e7x7 ->push_back( e7x7P.first );
t_e9x9 ->push_back( e9x9P.first );
t_e11x11 ->push_back( e11x11P.first );
//t_e13x13 ->push_back( e13x13P.first );
t_e15x15 ->push_back( e15x15P.first );
//t_e21x21 ->push_back( e21x21P.first );
//t_e25x25 ->push_back( e25x25P.first );
//t_e31x31 ->push_back( e31x31P.first );
t_e7x7_10Sig ->push_back( e7x7_10SigP.first );
t_e9x9_10Sig ->push_back( e9x9_10SigP.first );
t_e11x11_10Sig ->push_back( e11x11_10SigP.first );
t_e15x15_10Sig ->push_back( e15x15_10SigP.first );
t_e7x7_15Sig ->push_back( e7x7_15SigP.first );
t_e9x9_15Sig ->push_back( e9x9_15SigP.first );
t_e11x11_15Sig ->push_back( e11x11_15SigP.first );
t_e15x15_15Sig ->push_back( e15x15_15SigP.first );
t_e7x7_20Sig ->push_back( e7x7_20SigP.first );
t_e9x9_20Sig ->push_back( e9x9_20SigP.first );
t_e11x11_20Sig ->push_back( e11x11_20SigP.first );
t_e15x15_20Sig ->push_back( e15x15_20SigP.first );
t_e7x7_25Sig ->push_back( e7x7_25SigP.first );
t_e9x9_25Sig ->push_back( e9x9_25SigP.first );
t_e11x11_25Sig ->push_back( e11x11_25SigP.first );
t_e15x15_25Sig ->push_back( e15x15_25SigP.first );
t_e7x7_30Sig ->push_back( e7x7_30SigP.first );
t_e9x9_30Sig ->push_back( e9x9_30SigP.first );
t_e11x11_30Sig ->push_back( e11x11_30SigP.first );
t_e15x15_30Sig ->push_back( e15x15_30SigP.first );
if (doMC) {
//t_esim3x3 ->push_back( simInfo3x3.eTotal );
//t_esim5x5 ->push_back( simInfo5x5.eTotal );
t_esim7x7 ->push_back( simInfo7x7.eTotal );
t_esim9x9 ->push_back( simInfo9x9.eTotal );
t_esim11x11 ->push_back( simInfo11x11.eTotal );
//t_esim13x13 ->push_back( simInfo13x13.eTotal );
t_esim15x15 ->push_back( simInfo15x15.eTotal );
//t_esim21x21 ->push_back( simInfo21x21.eTotal );
//t_esim25x25 ->push_back( simInfo25x25.eTotal );
//t_esim31x31 ->push_back( simInfo31x31.eTotal );
//t_esim3x3Matched ->push_back( simInfo3x3.eMatched );
//t_esim5x5Matched ->push_back( simInfo5x5.eMatched );
t_esim7x7Matched ->push_back( simInfo7x7.eMatched );
t_esim9x9Matched ->push_back( simInfo9x9.eMatched );
t_esim11x11Matched ->push_back( simInfo11x11.eMatched );
//t_esim13x13Matched ->push_back( simInfo13x13.eMatched );
t_esim15x15Matched ->push_back( simInfo15x15.eMatched );
//t_esim21x21Matched ->push_back( simInfo21x21.eMatched );
//t_esim25x25Matched ->push_back( simInfo25x25.eMatched );
//t_esim31x31Matched ->push_back( simInfo31x31.eMatched );
//t_esim3x3Rest ->push_back( simInfo3x3.eRest );
//t_esim5x5Rest ->push_back( simInfo5x5.eRest );
t_esim7x7Rest ->push_back( simInfo7x7.eRest );
t_esim9x9Rest ->push_back( simInfo9x9.eRest );
t_esim11x11Rest ->push_back( simInfo11x11.eRest );
//t_esim13x13Rest ->push_back( simInfo13x13.eRest );
t_esim15x15Rest ->push_back( simInfo15x15.eRest );
//t_esim21x21Rest ->push_back( simInfo21x21.eRest );
//t_esim25x25Rest ->push_back( simInfo25x25.eRest );
//t_esim31x31Rest ->push_back( simInfo31x31.eRest );
//t_esim3x3Photon ->push_back( simInfo3x3.eGamma );
//t_esim5x5Photon ->push_back( simInfo5x5.eGamma );
t_esim7x7Photon ->push_back( simInfo7x7.eGamma );
t_esim9x9Photon ->push_back( simInfo9x9.eGamma );
t_esim11x11Photon ->push_back( simInfo11x11.eGamma );
//t_esim13x13Photon ->push_back( simInfo13x13.eGamma );
t_esim15x15Photon ->push_back( simInfo15x15.eGamma );
//t_esim21x21Photon ->push_back( simInfo21x21.eGamma );
//t_esim25x25Photon ->push_back( simInfo25x25.eGamma );
//t_esim31x31Photon ->push_back( simInfo31x31.eGamma );
//t_esim3x3NeutHad ->push_back( simInfo3x3.eNeutralHad );
//t_esim5x5NeutHad ->push_back( simInfo5x5.eNeutralHad );
t_esim7x7NeutHad ->push_back( simInfo7x7.eNeutralHad );
t_esim9x9NeutHad ->push_back( simInfo9x9.eNeutralHad );
t_esim11x11NeutHad ->push_back( simInfo11x11.eNeutralHad );
//t_esim13x13NeutHad ->push_back( simInfo13x13.eNeutralHad );
t_esim15x15NeutHad ->push_back( simInfo15x15.eNeutralHad );
//t_esim21x21NeutHad ->push_back( simInfo21x21.eNeutralHad );
//t_esim25x25NeutHad ->push_back( simInfo25x25.eNeutralHad );
//t_esim31x31NeutHad ->push_back( simInfo31x31.eNeutralHad );
//t_esim3x3CharHad ->push_back( simInfo3x3.eChargedHad );
//t_esim5x5CharHad ->push_back( simInfo5x5.eChargedHad );
t_esim7x7CharHad ->push_back( simInfo7x7.eChargedHad );
t_esim9x9CharHad ->push_back( simInfo9x9.eChargedHad );
t_esim11x11CharHad ->push_back( simInfo11x11.eChargedHad );
//t_esim13x13CharHad ->push_back( simInfo13x13.eChargedHad );
t_esim15x15CharHad ->push_back( simInfo15x15.eChargedHad );
//t_esim21x21CharHad ->push_back( simInfo21x21.eChargedHad );
//t_esim25x25CharHad ->push_back( simInfo25x25.eChargedHad );
//t_esim31x31CharHad ->push_back( simInfo31x31.eChargedHad );
t_trkEcalEne ->push_back( trkEcalEne );
t_simTrackP ->push_back( simTrackP );
t_esimPdgId ->push_back( simInfo11x11.pdgMatched );
}
t_maxNearHcalP3x3 ->push_back( maxNearHcalP3x3 );
t_maxNearHcalP5x5 ->push_back( maxNearHcalP5x5 );
t_maxNearHcalP7x7 ->push_back( maxNearHcalP7x7 );
t_h3x3 ->push_back( h3x3 );
t_h5x5 ->push_back( h5x5 );
t_h7x7 ->push_back( h7x7 );
t_h3x3Sig ->push_back( h3x3Sig );
t_h5x5Sig ->push_back( h5x5Sig );
t_h7x7Sig ->push_back( h7x7Sig );
t_infoHcal ->push_back( trkDetItr->okHCAL );
if (doMC) {
t_trkHcalEne ->push_back( hcalScale*trkHcalEne );
t_hsim3x3 ->push_back( hcalScale*hsimInfo3x3.eTotal );
t_hsim5x5 ->push_back( hcalScale*hsimInfo5x5.eTotal );
t_hsim7x7 ->push_back( hcalScale*hsimInfo7x7.eTotal );
t_hsim3x3Matched ->push_back( hcalScale*hsimInfo3x3.eMatched );
t_hsim5x5Matched ->push_back( hcalScale*hsimInfo5x5.eMatched );
t_hsim7x7Matched ->push_back( hcalScale*hsimInfo7x7.eMatched );
t_hsim3x3Rest ->push_back( hcalScale*hsimInfo3x3.eRest );
t_hsim5x5Rest ->push_back( hcalScale*hsimInfo5x5.eRest );
t_hsim7x7Rest ->push_back( hcalScale*hsimInfo7x7.eRest );
t_hsim3x3Photon ->push_back( hcalScale*hsimInfo3x3.eGamma );
t_hsim5x5Photon ->push_back( hcalScale*hsimInfo5x5.eGamma );
t_hsim7x7Photon ->push_back( hcalScale*hsimInfo7x7.eGamma );
t_hsim3x3NeutHad ->push_back( hcalScale*hsimInfo3x3.eNeutralHad );
t_hsim5x5NeutHad ->push_back( hcalScale*hsimInfo5x5.eNeutralHad );
t_hsim7x7NeutHad ->push_back( hcalScale*hsimInfo7x7.eNeutralHad );
t_hsim3x3CharHad ->push_back( hcalScale*hsimInfo3x3.eChargedHad );
t_hsim5x5CharHad ->push_back( hcalScale*hsimInfo5x5.eChargedHad );
t_hsim7x7CharHad ->push_back( hcalScale*hsimInfo7x7.eChargedHad );
}
/*
if(hcalScale*hsimInfo3x3.eTotal > 50.0) {
std::cout << "Loosely Iso Track : eta " << eta1 << " Rec Mom " << p1 << " SimMom " << simTrackP << " h3x3 " << h3x3 << std::endl;
std::cout <<"Closest cell Hcal (atHCAL) " << (HcalDetId)ClosestCell << std::endl;
std::cout <<"trkHcalEne, etotal, matched, rest " <<hcalScale*trkHcalEne<<std::setw(15)<<hcalScale*hsimInfo3x3.eTotal
<<std::setw(15)<<hcalScale*hsimInfo3x3.eMatched<<std::setw(15)<<hcalScale*hsimInfo3x3.eRest
<<std::endl;
unsigned int nn = t_trkHcalEne->size();
std::cout <<"in Tree " << (*t_trkHcalEne)[nn-1] <<std::setw(15)<< (*t_hsim3x3)[nn-1]
<<std::setw(15)<< (*t_hsim3x3Matched)[nn-1] <<std::setw(15)<< (*t_hsim3x3Rest)[nn-1]
<< std::endl;
std::cout << "debug output \n" << std::endl;
spr::caloSimInfo hsimInfo3x3_debug;
spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, geo,pcalohh, SimTk, SimVtx, pTrack, *associate, 1,1, hsimInfo3x3_debug, 150.0, true);
}
*/
} // if loosely isolated track
} // check p1/eta
} // loop over track collection
//if (doMC) delete associate;
tree->Fill();
}
| void IsolatedTracksNxN::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 1022 of file IsolatedTracksNxN.cc.
References BookHistograms(), genPartEtaBins, genPartPBins, i, and nEventProc.
{
nEventProc=0;
// double tempgen_TH[21] = { 1.0, 2.0, 3.0, 4.0, 5.0,
double tempgen_TH[16] = { 0.0, 1.0, 2.0, 3.0, 4.0,
5.0, 6.0, 7.0, 9.0, 11.0,
15.0, 20.0, 30.0, 50.0, 75.0, 100.0};
for(int i=0; i<16; i++) genPartPBins[i] = tempgen_TH[i];
double tempgen_Eta[4] = {0.0, 1.131, 1.653, 2.172};
for(int i=0; i<4; i++) genPartEtaBins[i] = tempgen_Eta[i];
BookHistograms();
}
| void IsolatedTracksNxN::BookHistograms | ( | ) | [private] |
Definition at line 1317 of file IsolatedTracksNxN.cc.
References dir, doMC, fs, genPartEtaBins, genPartPBins, h_L1AlgoNames, h_maxNearP15x15, h_maxNearP21x21, h_maxNearP25x25, h_maxNearP31x31, h_nTracks, h_PVTracksWt, h_recEta_0, h_recEta_1, h_recEta_2, h_recEtaP_0, h_recEtaP_1, h_recEtaP_2, h_recEtaPt_0, h_recEtaPt_1, h_recEtaPt_2, h_recP_0, h_recP_1, h_recP_2, h_recPhi_0, h_recPhi_1, h_recPhi_2, h_recPt_0, h_recPt_1, h_recPt_2, TFileDirectory::make(), NEtaBins, NPBins, t_Bunch, t_e11x11, t_e11x11_10Sig, t_e11x11_15Sig, t_e11x11_20Sig, t_e11x11_25Sig, t_e11x11_30Sig, t_e13x13, t_e15x15, t_e15x15_10Sig, t_e15x15_15Sig, t_e15x15_20Sig, t_e15x15_25Sig, t_e15x15_30Sig, t_e21x21, t_e25x25, t_e31x31, t_e3x3, t_e5x5, t_e7x7, t_e7x7_10Sig, t_e7x7_15Sig, t_e7x7_20Sig, t_e7x7_25Sig, t_e7x7_30Sig, t_e9x9, t_e9x9_10Sig, t_e9x9_15Sig, t_e9x9_20Sig, t_e9x9_25Sig, t_e9x9_30Sig, t_ecalSpike11x11, t_esim11x11, t_esim11x11CharHad, t_esim11x11Matched, t_esim11x11NeutHad, t_esim11x11Photon, t_esim11x11Rest, t_esim13x13, t_esim13x13CharHad, t_esim13x13Matched, t_esim13x13NeutHad, t_esim13x13Photon, t_esim13x13Rest, t_esim15x15, t_esim15x15CharHad, t_esim15x15Matched, t_esim15x15NeutHad, t_esim15x15Photon, t_esim15x15Rest, t_esim21x21, t_esim21x21CharHad, t_esim21x21Matched, t_esim21x21NeutHad, t_esim21x21Photon, t_esim21x21Rest, t_esim25x25, t_esim25x25CharHad, t_esim25x25Matched, t_esim25x25NeutHad, t_esim25x25Photon, t_esim25x25Rest, t_esim31x31, t_esim31x31CharHad, t_esim31x31Matched, t_esim31x31NeutHad, t_esim31x31Photon, t_esim31x31Rest, t_esim3x3, t_esim3x3CharHad, t_esim3x3Matched, t_esim3x3NeutHad, t_esim3x3Photon, t_esim3x3Rest, t_esim5x5, t_esim5x5CharHad, t_esim5x5Matched, t_esim5x5NeutHad, t_esim5x5Photon, t_esim5x5Rest, t_esim7x7, t_esim7x7CharHad, t_esim7x7Matched, t_esim7x7NeutHad, t_esim7x7Photon, t_esim7x7Rest, t_esim9x9, t_esim9x9CharHad, t_esim9x9Matched, t_esim9x9NeutHad, t_esim9x9Photon, t_esim9x9Rest, t_esimPdgId, t_EvtNo, t_h3x3, t_h3x3Sig, t_h5x5, t_h5x5Sig, t_h7x7, t_h7x7Sig, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NeutHad, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NeutHad, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim7x7, t_hsim7x7CharHad, t_hsim7x7Matched, t_hsim7x7NeutHad, t_hsim7x7Photon, t_hsim7x7Rest, t_infoHcal, t_jetEta, t_jetPhi, t_jetPt, t_L1CenJetEta, t_L1CenJetPhi, t_L1CenJetPt, t_L1Decision, t_L1FwdJetEta, t_L1FwdJetPhi, t_L1FwdJetPt, t_L1IsoEMEta, t_L1IsoEMPhi, t_L1IsoEMPt, t_L1METEta, t_L1METPhi, t_L1METPt, t_L1MuonEta, t_L1MuonPhi, t_L1MuonPt, t_L1NonIsoEMEta, t_L1NonIsoEMPhi, t_L1NonIsoEMPt, t_L1TauJetEta, t_L1TauJetPhi, t_L1TauJetPt, t_Lumi, t_maxNearHcalP3x3, t_maxNearHcalP5x5, t_maxNearHcalP7x7, t_maxNearP11x11, t_maxNearP13x13, t_maxNearP15x15, t_maxNearP21x21, t_maxNearP25x25, t_maxNearP31x31, t_maxNearP7x7, t_maxNearP9x9, t_NLayersCrossed, t_nTracks, t_nTrksJetCalo, t_nTrksJetVtx, t_PVisValid, t_PVndof, t_PVNTracks, t_PVNTracksHP, t_PVNTracksHPWt, t_PVNTracksWt, t_PVTracksSumPt, t_PVTracksSumPtHP, t_PVTracksSumPtHPWt, t_PVTracksSumPtWt, t_PVx, t_PVy, t_PVz, t_RunNo, t_simTrackP, t_trackChiSq, t_trackChiSqAll, t_trackDxy, t_trackDxyAll, t_trackDxyBS, t_trackDxyPV, t_trackDxyPVAll, t_trackDz, t_trackDzAll, t_trackDzBS, t_trackDzPV, t_trackDzPVAll, t_trackEcalEta, t_trackEcalPhi, t_trackEta, t_trackEtaAll, t_trackHcalEta, t_trackHcalPhi, t_trackHitInMeasTEC, t_trackHitInMeasTIB, t_trackHitInMeasTID, t_trackHitInMeasTOB, t_trackHitInMissTEC, t_trackHitInMissTIB, t_trackHitInMissTIBTID, t_trackHitInMissTID, t_trackHitInMissTOB, t_trackHitOutMeasTEC, t_trackHitOutMeasTIB, t_trackHitOutMeasTID, t_trackHitOutMeasTOB, t_trackHitOutMissTEC, t_trackHitOutMissTIB, t_trackHitOutMissTID, t_trackHitOutMissTOB, t_trackHitOutMissTOBTEC, t_trackHitsTEC, t_trackHitsTOB, t_trackL, t_trackNOuterHits, t_trackOutPosOutHitDr, t_trackP, t_trackPAll, t_trackPdgIdAll, t_trackPhi, t_trackPhiAll, t_trackPt, t_trackPtAll, t_trackPVIdx, t_trkEcalEne, t_trkHcalEne, and tree.
Referenced by beginJob().
{
char hname[100], htit[100];
TFileDirectory dir = fs->mkdir("nearMaxTrackP");
for(unsigned int ieta=0; ieta<NEtaBins; ieta++) {
double lowEta=-5.0, highEta= 5.0;
lowEta = genPartEtaBins[ieta];
highEta = genPartEtaBins[ieta+1];
for(unsigned int ipt=0; ipt<NPBins; ipt++) {
double lowP=0.0, highP=300.0;
lowP = genPartPBins[ipt];
highP = genPartPBins[ipt+1];
sprintf(hname, "h_maxNearP31x31_ptBin%i_etaBin%i",ipt, ieta);
sprintf(htit, "maxNearP in 31x31 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
h_maxNearP31x31[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
h_maxNearP31x31[ipt][ieta] ->Sumw2();
sprintf(hname, "h_maxNearP25x25_ptBin%i_etaBin%i",ipt, ieta);
sprintf(htit, "maxNearP in 25x25 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
h_maxNearP25x25[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
h_maxNearP25x25[ipt][ieta] ->Sumw2();
sprintf(hname, "h_maxNearP21x21_ptBin%i_etaBin%i",ipt, ieta);
sprintf(htit, "maxNearP in 21x21 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
h_maxNearP21x21[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
h_maxNearP21x21[ipt][ieta] ->Sumw2();
sprintf(hname, "h_maxNearP15x15_ptBin%i_etaBin%i",ipt, ieta);
sprintf(htit, "maxNearP in 15x15 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
h_maxNearP15x15[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
h_maxNearP15x15[ipt][ieta] ->Sumw2();
}
}
h_L1AlgoNames = fs->make<TH1I>("h_L1AlgoNames", "h_L1AlgoNames:Bin Labels", 128, -0.5, 127.5);
// Reconstructed Tracks
h_PVTracksWt = fs->make<TH1F>("h_PVTracksWt", "h_PVTracksWt", 600, -0.1, 1.1);
h_nTracks = fs->make<TH1F>("h_nTracks", "h_nTracks", 1000, -0.5, 999.5);
sprintf(hname, "h_recEtaPt_0");
sprintf(htit, "h_recEtaPt (all tracks Eta vs pT)");
h_recEtaPt_0 = fs->make<TH2F>(hname, htit, 30, -3.0,3.0, 15, genPartPBins);
sprintf(hname, "h_recEtaP_0");
sprintf(htit, "h_recEtaP (all tracks Eta vs pT)");
h_recEtaP_0 = fs->make<TH2F>(hname, htit, 30, -3.0,3.0, 15, genPartPBins);
h_recPt_0 = fs->make<TH1F>("h_recPt_0", "Pt (all tracks)", 15, genPartPBins);
h_recP_0 = fs->make<TH1F>("h_recP_0", "P (all tracks)", 15, genPartPBins);
h_recEta_0 = fs->make<TH1F>("h_recEta_0", "Eta (all tracks)", 60, -3.0, 3.0);
h_recPhi_0 = fs->make<TH1F>("h_recPhi_0", "Phi (all tracks)", 100, -3.2, 3.2);
//-------------------------
sprintf(hname, "h_recEtaPt_1");
sprintf(htit, "h_recEtaPt (all good tracks Eta vs pT)");
h_recEtaPt_1 = fs->make<TH2F>(hname, htit, 30, -3.0,3.0, 15, genPartPBins);
sprintf(hname, "h_recEtaP_1");
sprintf(htit, "h_recEtaP (all good tracks Eta vs pT)");
h_recEtaP_1 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
h_recPt_1 = fs->make<TH1F>("h_recPt_1", "Pt (all good tracks)", 15, genPartPBins);
h_recP_1 = fs->make<TH1F>("h_recP_1", "P (all good tracks)", 15, genPartPBins);
h_recEta_1 = fs->make<TH1F>("h_recEta_1", "Eta (all good tracks)", 60, -3.0, 3.0);
h_recPhi_1 = fs->make<TH1F>("h_recPhi_1", "Phi (all good tracks)", 100, -3.2, 3.2);
//-------------------------
sprintf(hname, "h_recEtaPt_2");
sprintf(htit, "h_recEtaPt (charge isolation Eta vs pT)");
h_recEtaPt_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
sprintf(hname, "h_recEtaP_2");
sprintf(htit, "h_recEtaP (charge isolation Eta vs pT)");
h_recEtaP_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
h_recPt_2 = fs->make<TH1F>("h_recPt_2", "Pt (charge isolation)", 15, genPartPBins);
h_recP_2 = fs->make<TH1F>("h_recP_2", "P (charge isolation)", 15, genPartPBins);
h_recEta_2 = fs->make<TH1F>("h_recEta_2", "Eta (charge isolation)", 60, -3.0, 3.0);
h_recPhi_2 = fs->make<TH1F>("h_recPhi_2", "Phi (charge isolation)", 100, -3.2, 3.2);
tree = fs->make<TTree>("tree", "tree");
tree->SetAutoSave(10000);
tree->Branch("t_EvtNo" ,&t_EvtNo ,"t_EvtNo/I");
tree->Branch("t_RunNo" ,&t_RunNo ,"t_RunNo/I");
tree->Branch("t_Lumi" ,&t_Lumi ,"t_Lumi/I");
tree->Branch("t_Bunch" ,&t_Bunch ,"t_Bunch/I");
t_PVx = new std::vector<double>();
t_PVy = new std::vector<double>();
t_PVz = new std::vector<double>();
t_PVisValid = new std::vector<double>();
t_PVndof = new std::vector<double>();
t_PVNTracks = new std::vector<double>();
t_PVNTracksWt = new std::vector<double>();
t_PVTracksSumPt = new std::vector<double>();
t_PVTracksSumPtWt = new std::vector<double>();
t_PVNTracksHP = new std::vector<double>();
t_PVNTracksHPWt = new std::vector<double>();
t_PVTracksSumPtHP = new std::vector<double>();
t_PVTracksSumPtHPWt = new std::vector<double>();
tree->Branch("PVx" ,"vector<double>" ,&t_PVx);
tree->Branch("PVy" ,"vector<double>" ,&t_PVy);
tree->Branch("PVz" ,"vector<double>" ,&t_PVz);
tree->Branch("PVisValid" ,"vector<double>" ,&t_PVisValid);
tree->Branch("PVndof" ,"vector<double>" ,&t_PVndof);
tree->Branch("PVNTracks" ,"vector<double>" ,&t_PVNTracks);
tree->Branch("PVNTracksWt" ,"vector<double>" ,&t_PVNTracksWt);
tree->Branch("t_PVTracksSumPt" ,"vector<double>" ,&t_PVTracksSumPt);
tree->Branch("t_PVTracksSumPtWt" ,"vector<double>" ,&t_PVTracksSumPtWt);
tree->Branch("PVNTracksHP" ,"vector<double>" ,&t_PVNTracksHP);
tree->Branch("PVNTracksHPWt" ,"vector<double>" ,&t_PVNTracksHPWt);
tree->Branch("t_PVTracksSumPtHP" ,"vector<double>" ,&t_PVTracksSumPtHP);
tree->Branch("t_PVTracksSumPtHPWt" ,"vector<double>" ,&t_PVTracksSumPtHPWt);
//----- L1Trigger information
t_L1Decision = new std::vector<int>();
t_L1CenJetPt = new std::vector<double>();
t_L1CenJetEta = new std::vector<double>();
t_L1CenJetPhi = new std::vector<double>();
t_L1FwdJetPt = new std::vector<double>();
t_L1FwdJetEta = new std::vector<double>();
t_L1FwdJetPhi = new std::vector<double>();
t_L1TauJetPt = new std::vector<double>();
t_L1TauJetEta = new std::vector<double>();
t_L1TauJetPhi = new std::vector<double>();
t_L1MuonPt = new std::vector<double>();
t_L1MuonEta = new std::vector<double>();
t_L1MuonPhi = new std::vector<double>();
t_L1IsoEMPt = new std::vector<double>();
t_L1IsoEMEta = new std::vector<double>();
t_L1IsoEMPhi = new std::vector<double>();
t_L1NonIsoEMPt = new std::vector<double>();
t_L1NonIsoEMEta = new std::vector<double>();
t_L1NonIsoEMPhi = new std::vector<double>();
t_L1METPt = new std::vector<double>();
t_L1METEta = new std::vector<double>();
t_L1METPhi = new std::vector<double>();
tree->Branch("t_L1Decision", "vector<int>", &t_L1Decision);
tree->Branch("t_L1CenJetPt", "vector<double>", &t_L1CenJetPt);
tree->Branch("t_L1CenJetEta", "vector<double>", &t_L1CenJetEta);
tree->Branch("t_L1CenJetPhi", "vector<double>", &t_L1CenJetPhi);
tree->Branch("t_L1FwdJetPt", "vector<double>", &t_L1FwdJetPt);
tree->Branch("t_L1FwdJetEta", "vector<double>", &t_L1FwdJetEta);
tree->Branch("t_L1FwdJetPhi", "vector<double>", &t_L1FwdJetPhi);
tree->Branch("t_L1TauJetPt", "vector<double>", &t_L1TauJetPt);
tree->Branch("t_L1TauJetEta", "vector<double>", &t_L1TauJetEta);
tree->Branch("t_L1TauJetPhi", "vector<double>", &t_L1TauJetPhi);
tree->Branch("t_L1MuonPt", "vector<double>", &t_L1MuonPt);
tree->Branch("t_L1MuonEta", "vector<double>", &t_L1MuonEta);
tree->Branch("t_L1MuonPhi", "vector<double>", &t_L1MuonPhi);
tree->Branch("t_L1IsoEMPt", "vector<double>", &t_L1IsoEMPt);
tree->Branch("t_L1IsoEMEta", "vector<double>", &t_L1IsoEMEta);
tree->Branch("t_L1IsoEMPhi", "vector<double>", &t_L1IsoEMPhi);
tree->Branch("t_L1NonIsoEMPt", "vector<double>", &t_L1NonIsoEMPt);
tree->Branch("t_L1NonIsoEMEta", "vector<double>", &t_L1NonIsoEMEta);
tree->Branch("t_L1NonIsoEMPhi", "vector<double>", &t_L1NonIsoEMPhi);
tree->Branch("t_L1METPt", "vector<double>", &t_L1METPt);
tree->Branch("t_L1METEta", "vector<double>", &t_L1METEta);
tree->Branch("t_L1METPhi", "vector<double>", &t_L1METPhi);
t_jetPt = new std::vector<double>();
t_jetEta = new std::vector<double>();
t_jetPhi = new std::vector<double>();
t_nTrksJetCalo = new std::vector<double>();
t_nTrksJetVtx = new std::vector<double>();
tree->Branch("t_jetPt", "vector<double>",&t_jetPt);
tree->Branch("t_jetEta", "vector<double>",&t_jetEta);
tree->Branch("t_jetPhi", "vector<double>",&t_jetPhi);
tree->Branch("t_nTrksJetCalo", "vector<double>",&t_nTrksJetCalo);
tree->Branch("t_nTrksJetVtx", "vector<double>",&t_nTrksJetVtx);
t_trackPAll = new std::vector<double>();
t_trackEtaAll = new std::vector<double>();
t_trackPhiAll = new std::vector<double>();
t_trackPdgIdAll = new std::vector<double>();
t_trackPtAll = new std::vector<double>();
t_trackDxyAll = new std::vector<double>();
t_trackDzAll = new std::vector<double>();
t_trackDxyPVAll = new std::vector<double>();
t_trackDzPVAll = new std::vector<double>();
t_trackChiSqAll = new std::vector<double>();
tree->Branch("t_trackPAll", "vector<double>", &t_trackPAll );
tree->Branch("t_trackPhiAll", "vector<double>", &t_trackPhiAll );
tree->Branch("t_trackEtaAll", "vector<double>", &t_trackEtaAll );
tree->Branch("t_trackPtAll", "vector<double>", &t_trackPtAll );
tree->Branch("t_trackDxyAll", "vector<double>", &t_trackDxyAll );
tree->Branch("t_trackDzAll", "vector<double>", &t_trackDzAll );
tree->Branch("t_trackDxyPVAll", "vector<double>", &t_trackDxyPVAll );
tree->Branch("t_trackDzPVAll", "vector<double>", &t_trackDzPVAll );
tree->Branch("t_trackChiSqAll", "vector<double>", &t_trackChiSqAll );
//tree->Branch("t_trackPdgIdAll", "vector<double>", &t_trackPdgIdAll);
t_trackP = new std::vector<double>();
t_trackPt = new std::vector<double>();
t_trackEta = new std::vector<double>();
t_trackPhi = new std::vector<double>();
t_trackEcalEta = new std::vector<double>();
t_trackEcalPhi = new std::vector<double>();
t_trackHcalEta = new std::vector<double>();
t_trackHcalPhi = new std::vector<double>();
t_trackNOuterHits = new std::vector<int>();
t_NLayersCrossed = new std::vector<int>();
t_trackDxy = new std::vector<double>();
t_trackDxyBS = new std::vector<double>();
t_trackDz = new std::vector<double>();
t_trackDzBS = new std::vector<double>();
t_trackDxyPV = new std::vector<double>();
t_trackDzPV = new std::vector<double>();
t_trackPVIdx = new std::vector<int>();
t_trackChiSq = new std::vector<double>();
t_trackHitsTOB = new std::vector<int>();
t_trackHitsTEC = new std::vector<int>();
t_trackHitInMissTOB = new std::vector<int>();
t_trackHitInMissTEC = new std::vector<int>();
t_trackHitInMissTIB = new std::vector<int>();
t_trackHitInMissTID = new std::vector<int>();
t_trackHitInMissTIBTID= new std::vector<int>();
t_trackHitOutMissTOB = new std::vector<int>();
t_trackHitOutMissTEC = new std::vector<int>();
t_trackHitOutMissTIB = new std::vector<int>();
t_trackHitOutMissTID = new std::vector<int>();
t_trackHitOutMissTOBTEC= new std::vector<int>();
t_trackHitInMeasTOB = new std::vector<int>();
t_trackHitInMeasTEC = new std::vector<int>();
t_trackHitInMeasTIB = new std::vector<int>();
t_trackHitInMeasTID = new std::vector<int>();
t_trackHitOutMeasTOB = new std::vector<int>();
t_trackHitOutMeasTEC = new std::vector<int>();
t_trackHitOutMeasTIB = new std::vector<int>();
t_trackHitOutMeasTID = new std::vector<int>();
t_trackOutPosOutHitDr =new std::vector<double>();
t_trackL =new std::vector<double>();
tree->Branch("t_trackP", "vector<double>", &t_trackP );
tree->Branch("t_trackPt", "vector<double>", &t_trackPt );
tree->Branch("t_trackEta", "vector<double>", &t_trackEta );
tree->Branch("t_trackPhi", "vector<double>", &t_trackPhi );
tree->Branch("t_trackEcalEta", "vector<double>", &t_trackEcalEta );
tree->Branch("t_trackEcalPhi", "vector<double>", &t_trackEcalPhi );
tree->Branch("t_trackHcalEta", "vector<double>", &t_trackHcalEta );
tree->Branch("t_trackHcalPhi", "vector<double>", &t_trackHcalPhi );
tree->Branch("t_trackNOuterHits", "vector<int>", &t_trackNOuterHits );
tree->Branch("t_NLayersCrossed", "vector<int>", &t_NLayersCrossed );
tree->Branch("t_trackHitsTOB", "vector<int>", &t_trackHitsTOB );
tree->Branch("t_trackHitsTEC", "vector<int>", &t_trackHitsTEC );
tree->Branch("t_trackHitInMissTOB", "vector<int>", &t_trackHitInMissTOB );
tree->Branch("t_trackHitInMissTEC", "vector<int>", &t_trackHitInMissTEC );
tree->Branch("t_trackHitInMissTIB", "vector<int>", &t_trackHitInMissTIB );
tree->Branch("t_trackHitInMissTID", "vector<int>", &t_trackHitInMissTID );
tree->Branch("t_trackHitInMissTIBTID", "vector<int>", &t_trackHitInMissTIBTID );
tree->Branch("t_trackHitOutMissTOB", "vector<int>", &t_trackHitOutMissTOB);
tree->Branch("t_trackHitOutMissTEC", "vector<int>", &t_trackHitOutMissTEC);
tree->Branch("t_trackHitOutMissTIB", "vector<int>", &t_trackHitOutMissTIB);
tree->Branch("t_trackHitOutMissTID", "vector<int>", &t_trackHitOutMissTID);
tree->Branch("t_trackHitOutMissTOBTEC","vector<int>", &t_trackHitOutMissTOBTEC);
tree->Branch("t_trackHitInMeasTOB", "vector<int>", &t_trackHitInMeasTOB );
tree->Branch("t_trackHitInMeasTEC", "vector<int>", &t_trackHitInMeasTEC );
tree->Branch("t_trackHitInMeasTIB", "vector<int>", &t_trackHitInMeasTIB );
tree->Branch("t_trackHitInMeasTID", "vector<int>", &t_trackHitInMeasTID );
tree->Branch("t_trackHitOutMeasTOB", "vector<int>", &t_trackHitOutMeasTOB);
tree->Branch("t_trackHitOutMeasTEC", "vector<int>", &t_trackHitOutMeasTEC);
tree->Branch("t_trackHitOutMeasTIB", "vector<int>", &t_trackHitOutMeasTIB);
tree->Branch("t_trackHitOutMeasTID", "vector<int>", &t_trackHitOutMeasTID);
tree->Branch("t_trackOutPosOutHitDr", "vector<double>", &t_trackOutPosOutHitDr);
tree->Branch("t_trackL", "vector<double>", &t_trackL);
tree->Branch("t_trackDxy", "vector<double>", &t_trackDxy );
tree->Branch("t_trackDxyBS", "vector<double>", &t_trackDxyBS );
tree->Branch("t_trackDz", "vector<double>", &t_trackDz );
tree->Branch("t_trackDzBS", "vector<double>", &t_trackDzBS );
tree->Branch("t_trackDxyPV", "vector<double>", &t_trackDxyPV );
tree->Branch("t_trackDzPV", "vector<double>", &t_trackDzPV );
tree->Branch("t_trackChiSq", "vector<double>", &t_trackChiSq );
tree->Branch("t_trackPVIdx", "vector<int>", &t_trackPVIdx );
t_maxNearP31x31 = new std::vector<double>();
t_maxNearP25x25 = new std::vector<double>();
t_maxNearP21x21 = new std::vector<double>();
t_maxNearP15x15 = new std::vector<double>();
t_maxNearP13x13 = new std::vector<double>();
t_maxNearP11x11 = new std::vector<double>();
t_maxNearP9x9 = new std::vector<double>();
t_maxNearP7x7 = new std::vector<double>();
tree->Branch("t_maxNearP31x31", "vector<double>", &t_maxNearP31x31);
//tree->Branch("t_maxNearP25x25", "vector<double>", &t_maxNearP25x25);
tree->Branch("t_maxNearP21x21", "vector<double>", &t_maxNearP21x21);
//tree->Branch("t_maxNearP15x15", "vector<double>", &t_maxNearP15x15);
//tree->Branch("t_maxNearP13x13", "vector<double>", &t_maxNearP13x13);
//tree->Branch("t_maxNearP11x11", "vector<double>", &t_maxNearP11x11);
//tree->Branch("t_maxNearP9x9", "vector<double>", &t_maxNearP9x9);
//tree->Branch("t_maxNearP7x7", "vector<double>", &t_maxNearP7x7);
t_ecalSpike11x11 = new std::vector<int>();
t_e3x3 = new std::vector<double>();
t_e5x5 = new std::vector<double>();
t_e7x7 = new std::vector<double>();
t_e9x9 = new std::vector<double>();
t_e11x11 = new std::vector<double>();
t_e13x13 = new std::vector<double>();
t_e15x15 = new std::vector<double>();
t_e21x21 = new std::vector<double>();
t_e25x25 = new std::vector<double>();
t_e31x31 = new std::vector<double>();
//t_e11x11Xtals = new std::vector<std::vector<double> >();
tree->Branch("t_ecalSpike11x11", "vector<int>", &t_ecalSpike11x11);
//tree->Branch("t_e3x3", "vector<double>", &t_e3x3);
//tree->Branch("t_e5x5", "vector<double>", &t_e5x5);
tree->Branch("t_e7x7", "vector<double>", &t_e7x7);
tree->Branch("t_e9x9", "vector<double>", &t_e9x9);
tree->Branch("t_e11x11", "vector<double>", &t_e11x11);
//tree->Branch("t_e13x13", "vector<double>", &t_e13x13);
tree->Branch("t_e15x15", "vector<double>", &t_e15x15);
//tree->Branch("t_e21x21", "vector<double>", &t_e21x21);
//tree->Branch("t_e25x25", "vector<double>", &t_e25x25);
//tree->Branch("t_e31x31", "vector<double>", &t_e31x31);
//tree->Branch("t_e11x11Xtals", "vector<vector<double> >", &t_e11x11Xtals);
t_e7x7_10Sig = new std::vector<double>();
t_e9x9_10Sig = new std::vector<double>();
t_e11x11_10Sig = new std::vector<double>();
t_e15x15_10Sig = new std::vector<double>();
t_e7x7_15Sig = new std::vector<double>();
t_e9x9_15Sig = new std::vector<double>();
t_e11x11_15Sig = new std::vector<double>();
t_e15x15_15Sig = new std::vector<double>();
t_e7x7_20Sig = new std::vector<double>();
t_e9x9_20Sig = new std::vector<double>();
t_e11x11_20Sig = new std::vector<double>();
t_e15x15_20Sig = new std::vector<double>();
t_e7x7_25Sig = new std::vector<double>();
t_e9x9_25Sig = new std::vector<double>();
t_e11x11_25Sig = new std::vector<double>();
t_e15x15_25Sig = new std::vector<double>();
t_e7x7_30Sig = new std::vector<double>();
t_e9x9_30Sig = new std::vector<double>();
t_e11x11_30Sig = new std::vector<double>();
t_e15x15_30Sig = new std::vector<double>();
tree->Branch("t_e7x7_10Sig" ,"vector<double>", &t_e7x7_10Sig);
tree->Branch("t_e9x9_10Sig" ,"vector<double>", &t_e9x9_10Sig);
tree->Branch("t_e11x11_10Sig" ,"vector<double>", &t_e11x11_10Sig);
tree->Branch("t_e15x15_10Sig" ,"vector<double>", &t_e15x15_10Sig);
tree->Branch("t_e7x7_15Sig" ,"vector<double>", &t_e7x7_15Sig);
tree->Branch("t_e9x9_15Sig" ,"vector<double>", &t_e9x9_15Sig);
tree->Branch("t_e11x11_15Sig" ,"vector<double>", &t_e11x11_15Sig);
tree->Branch("t_e15x15_15Sig" ,"vector<double>", &t_e15x15_15Sig);
tree->Branch("t_e7x7_20Sig" ,"vector<double>", &t_e7x7_20Sig);
tree->Branch("t_e9x9_20Sig" ,"vector<double>", &t_e9x9_20Sig);
tree->Branch("t_e11x11_20Sig" ,"vector<double>", &t_e11x11_20Sig);
tree->Branch("t_e15x15_20Sig" ,"vector<double>", &t_e15x15_20Sig);
tree->Branch("t_e7x7_25Sig" ,"vector<double>", &t_e7x7_25Sig);
tree->Branch("t_e9x9_25Sig" ,"vector<double>", &t_e9x9_25Sig);
tree->Branch("t_e11x11_25Sig" ,"vector<double>", &t_e11x11_25Sig);
tree->Branch("t_e15x15_25Sig" ,"vector<double>", &t_e15x15_25Sig);
tree->Branch("t_e7x7_30Sig" ,"vector<double>", &t_e7x7_30Sig);
tree->Branch("t_e9x9_30Sig" ,"vector<double>", &t_e9x9_30Sig);
tree->Branch("t_e11x11_30Sig" ,"vector<double>", &t_e11x11_30Sig);
tree->Branch("t_e15x15_30Sig" ,"vector<double>", &t_e15x15_30Sig);
if (doMC) {
t_esim3x3 = new std::vector<double>();
t_esim5x5 = new std::vector<double>();
t_esim7x7 = new std::vector<double>();
t_esim9x9 = new std::vector<double>();
t_esim11x11 = new std::vector<double>();
t_esim13x13 = new std::vector<double>();
t_esim15x15 = new std::vector<double>();
t_esim21x21 = new std::vector<double>();
t_esim25x25 = new std::vector<double>();
t_esim31x31 = new std::vector<double>();
t_esim3x3Matched = new std::vector<double>();
t_esim5x5Matched = new std::vector<double>();
t_esim7x7Matched = new std::vector<double>();
t_esim9x9Matched = new std::vector<double>();
t_esim11x11Matched = new std::vector<double>();
t_esim13x13Matched = new std::vector<double>();
t_esim15x15Matched = new std::vector<double>();
t_esim21x21Matched = new std::vector<double>();
t_esim25x25Matched = new std::vector<double>();
t_esim31x31Matched = new std::vector<double>();
t_esim3x3Rest = new std::vector<double>();
t_esim5x5Rest = new std::vector<double>();
t_esim7x7Rest = new std::vector<double>();
t_esim9x9Rest = new std::vector<double>();
t_esim11x11Rest = new std::vector<double>();
t_esim13x13Rest = new std::vector<double>();
t_esim15x15Rest = new std::vector<double>();
t_esim21x21Rest = new std::vector<double>();
t_esim25x25Rest = new std::vector<double>();
t_esim31x31Rest = new std::vector<double>();
t_esim3x3Photon = new std::vector<double>();
t_esim5x5Photon = new std::vector<double>();
t_esim7x7Photon = new std::vector<double>();
t_esim9x9Photon = new std::vector<double>();
t_esim11x11Photon = new std::vector<double>();
t_esim13x13Photon = new std::vector<double>();
t_esim15x15Photon = new std::vector<double>();
t_esim21x21Photon = new std::vector<double>();
t_esim25x25Photon = new std::vector<double>();
t_esim31x31Photon = new std::vector<double>();
t_esim3x3NeutHad = new std::vector<double>();
t_esim5x5NeutHad = new std::vector<double>();
t_esim7x7NeutHad = new std::vector<double>();
t_esim9x9NeutHad = new std::vector<double>();
t_esim11x11NeutHad = new std::vector<double>();
t_esim13x13NeutHad = new std::vector<double>();
t_esim15x15NeutHad = new std::vector<double>();
t_esim21x21NeutHad = new std::vector<double>();
t_esim25x25NeutHad = new std::vector<double>();
t_esim31x31NeutHad = new std::vector<double>();
t_esim3x3CharHad = new std::vector<double>();
t_esim5x5CharHad = new std::vector<double>();
t_esim7x7CharHad = new std::vector<double>();
t_esim9x9CharHad = new std::vector<double>();
t_esim11x11CharHad = new std::vector<double>();
t_esim13x13CharHad = new std::vector<double>();
t_esim15x15CharHad = new std::vector<double>();
t_esim21x21CharHad = new std::vector<double>();
t_esim25x25CharHad = new std::vector<double>();
t_esim31x31CharHad = new std::vector<double>();
t_trkEcalEne = new std::vector<double>();
t_simTrackP = new std::vector<double>();
t_esimPdgId = new std::vector<double>();
//tree->Branch("t_esim3x3", "vector<double>", &t_esim3x3);
//tree->Branch("t_esim5x5", "vector<double>", &t_esim5x5);
tree->Branch("t_esim7x7", "vector<double>", &t_esim7x7);
tree->Branch("t_esim9x9", "vector<double>", &t_esim9x9);
tree->Branch("t_esim11x11", "vector<double>", &t_esim11x11);
//tree->Branch("t_esim13x13", "vector<double>", &t_esim13x13);
tree->Branch("t_esim15x15", "vector<double>", &t_esim15x15);
//tree->Branch("t_esim21x21", "vector<double>", &t_esim21x21);
//tree->Branch("t_esim25x25", "vector<double>", &t_esim25x25);
//tree->Branch("t_esim31x31", "vector<double>", &t_esim31x31);
//tree->Branch("t_esim3x3Matched", "vector<double>", &t_esim3x3Matched);
//tree->Branch("t_esim5x5Matched", "vector<double>", &t_esim5x5Matched);
tree->Branch("t_esim7x7Matched", "vector<double>", &t_esim7x7Matched);
tree->Branch("t_esim9x9Matched", "vector<double>", &t_esim9x9Matched);
tree->Branch("t_esim11x11Matched", "vector<double>", &t_esim11x11Matched);
//tree->Branch("t_esim13x13Matched", "vector<double>", &t_esim13x13Matched);
tree->Branch("t_esim15x15Matched", "vector<double>", &t_esim15x15Matched);
//tree->Branch("t_esim21x21Matched", "vector<double>", &t_esim21x21Matched);
//tree->Branch("t_esim25x25Matched", "vector<double>", &t_esim25x25Matched);
//tree->Branch("t_esim31x31Matched", "vector<double>", &t_esim31x31Matched);
//tree->Branch("t_esim3x3Rest", "vector<double>", &t_esim3x3Rest);
//tree->Branch("t_esim5x5Rest", "vector<double>", &t_esim5x5Rest);
tree->Branch("t_esim7x7Rest", "vector<double>", &t_esim7x7Rest);
tree->Branch("t_esim9x9Rest", "vector<double>", &t_esim9x9Rest);
tree->Branch("t_esim11x11Rest", "vector<double>", &t_esim11x11Rest);
//tree->Branch("t_esim13x13Rest", "vector<double>", &t_esim13x13Rest);
tree->Branch("t_esim15x15Rest", "vector<double>", &t_esim15x15Rest);
//tree->Branch("t_esim21x21Rest", "vector<double>", &t_esim21x21Rest);
//tree->Branch("t_esim25x25Rest", "vector<double>", &t_esim25x25Rest);
//tree->Branch("t_esim31x31Rest", "vector<double>", &t_esim31x31Rest);
//tree->Branch("t_esim3x3Photon", "vector<double>", &t_esim3x3Photon);
//tree->Branch("t_esim5x5Photon", "vector<double>", &t_esim5x5Photon);
tree->Branch("t_esim7x7Photon", "vector<double>", &t_esim7x7Photon);
tree->Branch("t_esim9x9Photon", "vector<double>", &t_esim9x9Photon);
tree->Branch("t_esim11x11Photon", "vector<double>", &t_esim11x11Photon);
//tree->Branch("t_esim13x13Photon", "vector<double>", &t_esim13x13Photon);
tree->Branch("t_esim15x15Photon", "vector<double>", &t_esim15x15Photon);
//tree->Branch("t_esim21x21Photon", "vector<double>", &t_esim21x21Photon);
//tree->Branch("t_esim25x25Photon", "vector<double>", &t_esim25x25Photon);
//tree->Branch("t_esim31x31Photon", "vector<double>", &t_esim31x31Photon);
//tree->Branch("t_esim3x3NeutHad", "vector<double>", &t_esim3x3NeutHad);
//tree->Branch("t_esim5x5NeutHad", "vector<double>", &t_esim5x5NeutHad);
tree->Branch("t_esim7x7NeutHad", "vector<double>", &t_esim7x7NeutHad);
tree->Branch("t_esim9x9NeutHad", "vector<double>", &t_esim9x9NeutHad);
tree->Branch("t_esim11x11NeutHad", "vector<double>", &t_esim11x11NeutHad);
//tree->Branch("t_esim13x13NeutHad", "vector<double>", &t_esim13x13NeutHad);
tree->Branch("t_esim15x15NeutHad", "vector<double>", &t_esim15x15NeutHad);
//tree->Branch("t_esim21x21NeutHad", "vector<double>", &t_esim21x21NeutHad);
//tree->Branch("t_esim25x25NeutHad", "vector<double>", &t_esim25x25NeutHad);
//tree->Branch("t_esim31x31NeutHad", "vector<double>", &t_esim31x31NeutHad);
//tree->Branch("t_esim3x3CharHad", "vector<double>", &t_esim3x3CharHad);
//tree->Branch("t_esim5x5CharHad", "vector<double>", &t_esim5x5CharHad);
tree->Branch("t_esim7x7CharHad", "vector<double>", &t_esim7x7CharHad);
tree->Branch("t_esim9x9CharHad", "vector<double>", &t_esim9x9CharHad);
tree->Branch("t_esim11x11CharHad", "vector<double>", &t_esim11x11CharHad);
//tree->Branch("t_esim13x13CharHad", "vector<double>", &t_esim13x13CharHad);
tree->Branch("t_esim15x15CharHad", "vector<double>", &t_esim15x15CharHad);
//tree->Branch("t_esim21x21CharHad", "vector<double>", &t_esim21x21CharHad);
//tree->Branch("t_esim25x25CharHad", "vector<double>", &t_esim25x25CharHad);
//tree->Branch("t_esim31x31CharHad", "vector<double>", &t_esim31x31CharHad);
tree->Branch("t_trkEcalEne", "vector<double>", &t_trkEcalEne);
tree->Branch("t_simTrackP", "vector<double>", &t_simTrackP);
tree->Branch("t_esimPdgId", "vector<double>", &t_esimPdgId);
}
t_maxNearHcalP3x3 = new std::vector<double>();
t_maxNearHcalP5x5 = new std::vector<double>();
t_maxNearHcalP7x7 = new std::vector<double>();
t_h3x3 = new std::vector<double>();
t_h5x5 = new std::vector<double>();
t_h7x7 = new std::vector<double>();
t_h3x3Sig = new std::vector<double>();
t_h5x5Sig = new std::vector<double>();
t_h7x7Sig = new std::vector<double>();
t_infoHcal = new std::vector<int>();
if (doMC) {
t_trkHcalEne = new std::vector<double>();
t_hsim3x3 = new std::vector<double>();
t_hsim5x5 = new std::vector<double>();
t_hsim7x7 = new std::vector<double>();
t_hsim3x3Matched = new std::vector<double>();
t_hsim5x5Matched = new std::vector<double>();
t_hsim7x7Matched = new std::vector<double>();
t_hsim3x3Rest = new std::vector<double>();
t_hsim5x5Rest = new std::vector<double>();
t_hsim7x7Rest = new std::vector<double>();
t_hsim3x3Photon = new std::vector<double>();
t_hsim5x5Photon = new std::vector<double>();
t_hsim7x7Photon = new std::vector<double>();
t_hsim3x3NeutHad = new std::vector<double>();
t_hsim5x5NeutHad = new std::vector<double>();
t_hsim7x7NeutHad = new std::vector<double>();
t_hsim3x3CharHad = new std::vector<double>();
t_hsim5x5CharHad = new std::vector<double>();
t_hsim7x7CharHad = new std::vector<double>();
}
tree->Branch("t_maxNearHcalP3x3", "vector<double>", &t_maxNearHcalP3x3);
tree->Branch("t_maxNearHcalP5x5", "vector<double>", &t_maxNearHcalP5x5);
tree->Branch("t_maxNearHcalP7x7", "vector<double>", &t_maxNearHcalP7x7);
tree->Branch("t_h3x3", "vector<double>", &t_h3x3);
tree->Branch("t_h5x5", "vector<double>", &t_h5x5);
tree->Branch("t_h7x7", "vector<double>", &t_h7x7);
tree->Branch("t_h3x3Sig", "vector<double>", &t_h3x3Sig);
tree->Branch("t_h5x5Sig", "vector<double>", &t_h5x5Sig);
tree->Branch("t_h7x7Sig", "vector<double>", &t_h7x7Sig);
tree->Branch("t_infoHcal", "vector<int>", &t_infoHcal);
if (doMC) {
tree->Branch("t_trkHcalEne", "vector<double>", &t_trkHcalEne);
tree->Branch("t_hsim3x3", "vector<double>", &t_hsim3x3);
tree->Branch("t_hsim5x5", "vector<double>", &t_hsim5x5);
tree->Branch("t_hsim7x7", "vector<double>", &t_hsim7x7);
tree->Branch("t_hsim3x3Matched", "vector<double>", &t_hsim3x3Matched);
tree->Branch("t_hsim5x5Matched", "vector<double>", &t_hsim5x5Matched);
tree->Branch("t_hsim7x7Matched", "vector<double>", &t_hsim7x7Matched);
tree->Branch("t_hsim3x3Rest", "vector<double>", &t_hsim3x3Rest);
tree->Branch("t_hsim5x5Rest", "vector<double>", &t_hsim5x5Rest);
tree->Branch("t_hsim7x7Rest", "vector<double>", &t_hsim7x7Rest);
tree->Branch("t_hsim3x3Photon", "vector<double>", &t_hsim3x3Photon);
tree->Branch("t_hsim5x5Photon", "vector<double>", &t_hsim5x5Photon);
tree->Branch("t_hsim7x7Photon", "vector<double>", &t_hsim7x7Photon);
tree->Branch("t_hsim3x3NeutHad", "vector<double>", &t_hsim3x3NeutHad);
tree->Branch("t_hsim5x5NeutHad", "vector<double>", &t_hsim5x5NeutHad);
tree->Branch("t_hsim7x7NeutHad", "vector<double>", &t_hsim7x7NeutHad);
tree->Branch("t_hsim3x3CharHad", "vector<double>", &t_hsim3x3CharHad);
tree->Branch("t_hsim5x5CharHad", "vector<double>", &t_hsim5x5CharHad);
tree->Branch("t_hsim7x7CharHad", "vector<double>", &t_hsim7x7CharHad);
}
tree->Branch("t_nTracks", &t_nTracks, "t_nTracks/I");
}
| void IsolatedTracksNxN::clearTreeVectors | ( | ) | [private] |
Definition at line 1067 of file IsolatedTracksNxN.cc.
References doMC, t_e11x11, t_e11x11_10Sig, t_e11x11_15Sig, t_e11x11_20Sig, t_e11x11_25Sig, t_e11x11_30Sig, t_e13x13, t_e15x15, t_e15x15_10Sig, t_e15x15_15Sig, t_e15x15_20Sig, t_e15x15_25Sig, t_e15x15_30Sig, t_e21x21, t_e25x25, t_e31x31, t_e3x3, t_e5x5, t_e7x7, t_e7x7_10Sig, t_e7x7_15Sig, t_e7x7_20Sig, t_e7x7_25Sig, t_e7x7_30Sig, t_e9x9, t_e9x9_10Sig, t_e9x9_15Sig, t_e9x9_20Sig, t_e9x9_25Sig, t_e9x9_30Sig, t_ecalSpike11x11, t_esim11x11, t_esim11x11CharHad, t_esim11x11Matched, t_esim11x11NeutHad, t_esim11x11Photon, t_esim11x11Rest, t_esim13x13, t_esim13x13CharHad, t_esim13x13Matched, t_esim13x13NeutHad, t_esim13x13Photon, t_esim13x13Rest, t_esim15x15, t_esim15x15CharHad, t_esim15x15Matched, t_esim15x15NeutHad, t_esim15x15Photon, t_esim15x15Rest, t_esim21x21, t_esim21x21CharHad, t_esim21x21Matched, t_esim21x21NeutHad, t_esim21x21Photon, t_esim21x21Rest, t_esim25x25, t_esim25x25CharHad, t_esim25x25Matched, t_esim25x25NeutHad, t_esim25x25Photon, t_esim25x25Rest, t_esim31x31, t_esim31x31CharHad, t_esim31x31Matched, t_esim31x31NeutHad, t_esim31x31Photon, t_esim31x31Rest, t_esim3x3, t_esim3x3CharHad, t_esim3x3Matched, t_esim3x3NeutHad, t_esim3x3Photon, t_esim3x3Rest, t_esim5x5, t_esim5x5CharHad, t_esim5x5Matched, t_esim5x5NeutHad, t_esim5x5Photon, t_esim5x5Rest, t_esim7x7, t_esim7x7CharHad, t_esim7x7Matched, t_esim7x7NeutHad, t_esim7x7Photon, t_esim7x7Rest, t_esim9x9, t_esim9x9CharHad, t_esim9x9Matched, t_esim9x9NeutHad, t_esim9x9Photon, t_esim9x9Rest, t_esimPdgId, t_h3x3, t_h3x3Sig, t_h5x5, t_h5x5Sig, t_h7x7, t_h7x7Sig, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NeutHad, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NeutHad, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim7x7, t_hsim7x7CharHad, t_hsim7x7Matched, t_hsim7x7NeutHad, t_hsim7x7Photon, t_hsim7x7Rest, t_infoHcal, t_jetEta, t_jetPhi, t_jetPt, t_L1CenJetEta, t_L1CenJetPhi, t_L1CenJetPt, t_L1Decision, t_L1FwdJetEta, t_L1FwdJetPhi, t_L1FwdJetPt, t_L1IsoEMEta, t_L1IsoEMPhi, t_L1IsoEMPt, t_L1METEta, t_L1METPhi, t_L1METPt, t_L1MuonEta, t_L1MuonPhi, t_L1MuonPt, t_L1NonIsoEMEta, t_L1NonIsoEMPhi, t_L1NonIsoEMPt, t_L1TauJetEta, t_L1TauJetPhi, t_L1TauJetPt, t_maxNearHcalP3x3, t_maxNearHcalP5x5, t_maxNearHcalP7x7, t_maxNearP11x11, t_maxNearP13x13, t_maxNearP15x15, t_maxNearP21x21, t_maxNearP25x25, t_maxNearP31x31, t_maxNearP7x7, t_maxNearP9x9, t_NLayersCrossed, t_nTrksJetCalo, t_nTrksJetVtx, t_PVisValid, t_PVndof, t_PVNTracks, t_PVNTracksHP, t_PVNTracksHPWt, t_PVNTracksWt, t_PVTracksSumPt, t_PVTracksSumPtHP, t_PVTracksSumPtHPWt, t_PVTracksSumPtWt, t_PVx, t_PVy, t_PVz, t_simTrackP, t_trackChiSq, t_trackChiSqAll, t_trackDxy, t_trackDxyAll, t_trackDxyBS, t_trackDxyPV, t_trackDxyPVAll, t_trackDz, t_trackDzAll, t_trackDzBS, t_trackDzPV, t_trackDzPVAll, t_trackEcalEta, t_trackEcalPhi, t_trackEta, t_trackEtaAll, t_trackHcalEta, t_trackHcalPhi, t_trackHitInMeasTEC, t_trackHitInMeasTIB, t_trackHitInMeasTID, t_trackHitInMeasTOB, t_trackHitInMissTEC, t_trackHitInMissTIB, t_trackHitInMissTIBTID, t_trackHitInMissTID, t_trackHitInMissTOB, t_trackHitOutMeasTEC, t_trackHitOutMeasTIB, t_trackHitOutMeasTID, t_trackHitOutMeasTOB, t_trackHitOutMissTEC, t_trackHitOutMissTIB, t_trackHitOutMissTID, t_trackHitOutMissTOB, t_trackHitOutMissTOBTEC, t_trackHitsTEC, t_trackHitsTOB, t_trackL, t_trackNOuterHits, t_trackOutPosOutHitDr, t_trackP, t_trackPAll, t_trackPdgIdAll, t_trackPhi, t_trackPhiAll, t_trackPt, t_trackPtAll, t_trackPVIdx, t_trkEcalEne, and t_trkHcalEne.
Referenced by analyze().
{
t_PVx ->clear();
t_PVy ->clear();
t_PVz ->clear();
t_PVisValid ->clear();
t_PVndof ->clear();
t_PVNTracks ->clear();
t_PVNTracksWt ->clear();
t_PVTracksSumPt ->clear();
t_PVTracksSumPtWt ->clear();
t_PVNTracksHP ->clear();
t_PVNTracksHPWt ->clear();
t_PVTracksSumPtHP ->clear();
t_PVTracksSumPtHPWt ->clear();
t_L1Decision ->clear();
t_L1CenJetPt ->clear();
t_L1CenJetEta ->clear();
t_L1CenJetPhi ->clear();
t_L1FwdJetPt ->clear();
t_L1FwdJetEta ->clear();
t_L1FwdJetPhi ->clear();
t_L1TauJetPt ->clear();
t_L1TauJetEta ->clear();
t_L1TauJetPhi ->clear();
t_L1MuonPt ->clear();
t_L1MuonEta ->clear();
t_L1MuonPhi ->clear();
t_L1IsoEMPt ->clear();
t_L1IsoEMEta ->clear();
t_L1IsoEMPhi ->clear();
t_L1NonIsoEMPt ->clear();
t_L1NonIsoEMEta ->clear();
t_L1NonIsoEMPhi ->clear();
t_L1METPt ->clear();
t_L1METEta ->clear();
t_L1METPhi ->clear();
t_jetPt ->clear();
t_jetEta ->clear();
t_jetPhi ->clear();
t_nTrksJetCalo ->clear();
t_nTrksJetVtx ->clear();
t_trackPAll ->clear();
t_trackEtaAll ->clear();
t_trackPhiAll ->clear();
t_trackPdgIdAll ->clear();
t_trackPtAll ->clear();
t_trackDxyAll ->clear();
t_trackDzAll ->clear();
t_trackDxyPVAll ->clear();
t_trackDzPVAll ->clear();
t_trackChiSqAll ->clear();
t_trackP ->clear();
t_trackPt ->clear();
t_trackEta ->clear();
t_trackPhi ->clear();
t_trackEcalEta ->clear();
t_trackEcalPhi ->clear();
t_trackHcalEta ->clear();
t_trackHcalPhi ->clear();
t_NLayersCrossed ->clear();
t_trackNOuterHits ->clear();
t_trackDxy ->clear();
t_trackDxyBS ->clear();
t_trackDz ->clear();
t_trackDzBS ->clear();
t_trackDxyPV ->clear();
t_trackDzPV ->clear();
t_trackChiSq ->clear();
t_trackPVIdx ->clear();
t_trackHitsTOB ->clear();
t_trackHitsTEC ->clear();
t_trackHitInMissTOB ->clear();
t_trackHitInMissTEC ->clear();
t_trackHitInMissTIB ->clear();
t_trackHitInMissTID ->clear();
t_trackHitInMissTIBTID ->clear();
t_trackHitOutMissTOB ->clear();
t_trackHitOutMissTEC ->clear();
t_trackHitOutMissTIB ->clear();
t_trackHitOutMissTID ->clear();
t_trackHitOutMissTOBTEC ->clear();
t_trackHitInMeasTOB ->clear();
t_trackHitInMeasTEC ->clear();
t_trackHitInMeasTIB ->clear();
t_trackHitInMeasTID ->clear();
t_trackHitOutMeasTOB ->clear();
t_trackHitOutMeasTEC ->clear();
t_trackHitOutMeasTIB ->clear();
t_trackHitOutMeasTID ->clear();
t_trackOutPosOutHitDr ->clear();
t_trackL ->clear();
t_maxNearP31x31 ->clear();
t_maxNearP25x25 ->clear();
t_maxNearP21x21 ->clear();
t_maxNearP15x15 ->clear();
t_maxNearP13x13 ->clear();
t_maxNearP11x11 ->clear();
t_maxNearP9x9 ->clear();
t_maxNearP7x7 ->clear();
t_ecalSpike11x11 ->clear();
t_e3x3 ->clear();
t_e5x5 ->clear();
t_e7x7 ->clear();
t_e9x9 ->clear();
t_e11x11 ->clear();
t_e13x13 ->clear();
t_e15x15 ->clear();
t_e21x21 ->clear();
t_e25x25 ->clear();
t_e31x31 ->clear();
// t_e11x11Xtals ->clear();
t_e7x7_10Sig ->clear();
t_e9x9_10Sig ->clear();
t_e11x11_10Sig ->clear();
t_e15x15_10Sig ->clear();
t_e7x7_15Sig ->clear();
t_e9x9_15Sig ->clear();
t_e11x11_15Sig ->clear();
t_e15x15_15Sig ->clear();
t_e7x7_20Sig ->clear();
t_e9x9_20Sig ->clear();
t_e11x11_20Sig ->clear();
t_e15x15_20Sig ->clear();
t_e7x7_25Sig ->clear();
t_e9x9_25Sig ->clear();
t_e11x11_25Sig ->clear();
t_e15x15_25Sig ->clear();
t_e7x7_30Sig ->clear();
t_e9x9_30Sig ->clear();
t_e11x11_30Sig ->clear();
t_e15x15_30Sig ->clear();
if (doMC) {
t_simTrackP ->clear();
t_esimPdgId ->clear();
t_trkEcalEne ->clear();
t_esim3x3 ->clear();
t_esim5x5 ->clear();
t_esim7x7 ->clear();
t_esim9x9 ->clear();
t_esim11x11 ->clear();
t_esim13x13 ->clear();
t_esim15x15 ->clear();
t_esim21x21 ->clear();
t_esim25x25 ->clear();
t_esim31x31 ->clear();
t_esim3x3Matched ->clear();
t_esim5x5Matched ->clear();
t_esim7x7Matched ->clear();
t_esim9x9Matched ->clear();
t_esim11x11Matched ->clear();
t_esim13x13Matched ->clear();
t_esim15x15Matched ->clear();
t_esim21x21Matched ->clear();
t_esim25x25Matched ->clear();
t_esim31x31Matched ->clear();
t_esim3x3Rest ->clear();
t_esim5x5Rest ->clear();
t_esim7x7Rest ->clear();
t_esim9x9Rest ->clear();
t_esim11x11Rest ->clear();
t_esim13x13Rest ->clear();
t_esim15x15Rest ->clear();
t_esim21x21Rest ->clear();
t_esim25x25Rest ->clear();
t_esim31x31Rest ->clear();
t_esim3x3Photon ->clear();
t_esim5x5Photon ->clear();
t_esim7x7Photon ->clear();
t_esim9x9Photon ->clear();
t_esim11x11Photon ->clear();
t_esim13x13Photon ->clear();
t_esim15x15Photon ->clear();
t_esim21x21Photon ->clear();
t_esim25x25Photon ->clear();
t_esim31x31Photon ->clear();
t_esim3x3NeutHad ->clear();
t_esim5x5NeutHad ->clear();
t_esim7x7NeutHad ->clear();
t_esim9x9NeutHad ->clear();
t_esim11x11NeutHad ->clear();
t_esim13x13NeutHad ->clear();
t_esim15x15NeutHad ->clear();
t_esim21x21NeutHad ->clear();
t_esim25x25NeutHad ->clear();
t_esim31x31NeutHad ->clear();
t_esim3x3CharHad ->clear();
t_esim5x5CharHad ->clear();
t_esim7x7CharHad ->clear();
t_esim9x9CharHad ->clear();
t_esim11x11CharHad ->clear();
t_esim13x13CharHad ->clear();
t_esim15x15CharHad ->clear();
t_esim21x21CharHad ->clear();
t_esim25x25CharHad ->clear();
t_esim31x31CharHad ->clear();
}
t_maxNearHcalP3x3 ->clear();
t_maxNearHcalP5x5 ->clear();
t_maxNearHcalP7x7 ->clear();
t_h3x3 ->clear();
t_h5x5 ->clear();
t_h7x7 ->clear();
t_h3x3Sig ->clear();
t_h5x5Sig ->clear();
t_h7x7Sig ->clear();
t_infoHcal ->clear();
if (doMC) {
t_trkHcalEne ->clear();
t_hsim3x3 ->clear();
t_hsim5x5 ->clear();
t_hsim7x7 ->clear();
t_hsim3x3Matched ->clear();
t_hsim5x5Matched ->clear();
t_hsim7x7Matched ->clear();
t_hsim3x3Rest ->clear();
t_hsim5x5Rest ->clear();
t_hsim7x7Rest ->clear();
t_hsim3x3Photon ->clear();
t_hsim5x5Photon ->clear();
t_hsim7x7Photon ->clear();
t_hsim3x3NeutHad ->clear();
t_hsim5x5NeutHad ->clear();
t_hsim7x7NeutHad ->clear();
t_hsim3x3CharHad ->clear();
t_hsim5x5CharHad ->clear();
t_hsim7x7CharHad ->clear();
}
}
| double IsolatedTracksNxN::DeltaPhi | ( | double | v1, |
| double | v2 | ||
| ) | [private] |
Definition at line 1899 of file IsolatedTracksNxN.cc.
References abs, corr, diffTreeTool::diff, and pi.
Referenced by DeltaR().
| double IsolatedTracksNxN::DeltaR | ( | double | eta1, |
| double | phi1, | ||
| double | eta2, | ||
| double | phi2 | ||
| ) | [private] |
Definition at line 1911 of file IsolatedTracksNxN.cc.
References DeltaPhi(), and mathSSE::sqrt().
Referenced by analyze().
| void IsolatedTracksNxN::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 1041 of file IsolatedTracksNxN.cc.
References algoBitToName, cont, gather_cfg::cout, h_L1AlgoNames, l1TriggerMap, ExpressReco_HICollisions_FallBack::nbins, nEventProc, and trig_iter.
{
std::cout << "Number of Events Processed " << nEventProc << std::endl;
if( h_L1AlgoNames ){
int nbins=h_L1AlgoNames->GetNbinsX();
for (int ibin=0; ibin<nbins; ++ibin){
double cont = (double)h_L1AlgoNames->GetBinContent(ibin+1);
if( cont>0 ) {
trig_iter=l1TriggerMap.find(algoBitToName[ibin]);
const char* trigName = trig_iter->first.c_str();
h_L1AlgoNames->GetXaxis()->SetBinLabel(ibin+1,trigName);
std::cout<<"===============> "<<ibin+1<<" "<<trigName<<" "<<(int)cont<< std::endl;
}
}
}
}
| void IsolatedTracksNxN::printTrack | ( | const reco::Track * | pTrack | ) | [private] |
Definition at line 1917 of file IsolatedTracksNxN.cc.
References reco::TrackBase::charge(), gather_cfg::cout, reco::TrackBase::d0(), reco::TrackBase::eta(), reco::TrackBase::hitPattern(), i, reco::TrackBase::momentum(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfHits(), reco::TrackBase::numberOfLostHits(), reco::TrackBase::numberOfValidHits(), reco::TrackBase::p(), L1TEmulatorMonitor_cff::p, p1, p2, reco::TrackBase::phi(), reco::HitPattern::pixelBarrelLayersWithMeasurement(), reco::HitPattern::pixelEndcapLayersWithMeasurement(), reco::HitPattern::pixelLayersWithMeasurement(), reco::HitPattern::printHitPattern(), printTrkHitPattern_, reco::TrackBase::pt(), reco::TrackBase::qoverp(), reco::TrackBase::qoverpError(), reco::TrackBase::quality(), reco::TrackBase::qualityByName(), reco::TrackBase::qualityName(), reco::TrackBase::referencePoint(), reco::HitPattern::stripLayersWithMeasurement(), reco::HitPattern::stripTECLayersWithMeasurement(), reco::HitPattern::stripTIBLayersWithMeasurement(), reco::HitPattern::stripTIDLayersWithMeasurement(), reco::HitPattern::stripTOBLayersWithMeasurement(), reco::TrackBase::trackerExpectedHitsInner(), reco::TrackBase::trackerExpectedHitsOuter(), and reco::HitPattern::trackerLayersWithMeasurement().
{
std::string theTrackQuality = "highPurity";
reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);
std::cout << " Reference Point " << pTrack->referencePoint() <<"\n"
<< " TrackMmentum " << pTrack->momentum()
<< " (pt,eta,phi)(" << pTrack->pt()<<","<<pTrack->eta()<<","<<pTrack->phi()<<")"
<< " p " << pTrack->p() << "\n"
<< " Normalized chi2 " << pTrack->normalizedChi2() <<" charge " << pTrack->charge()
<< " qoverp() " << pTrack->qoverp() <<"+-" << pTrack->qoverpError()
<< " d0 " << pTrack->d0() << "\n"
<< " NValidHits " << pTrack->numberOfValidHits() << " NLostHits " << pTrack->numberOfLostHits()
<< " TrackQuality " << pTrack->qualityName(trackQuality_) << " " << pTrack->quality(trackQuality_)
<< std::endl;
if( printTrkHitPattern_ ) {
const reco::HitPattern& p = pTrack->hitPattern();
const reco::HitPattern& p1 = pTrack->trackerExpectedHitsInner();
const reco::HitPattern& p2 = pTrack->trackerExpectedHitsOuter();
std::cout<<"default " << std::endl;
for (int i=0; i<p.numberOfHits(); i++) {
p.printHitPattern(i, std::cout);
}
std::cout<<"trackerExpectedHitsInner() " << std::endl;
for (int i=0; i<p1.numberOfHits(); i++) {
p1.printHitPattern(i, std::cout);
}
std::cout<<"trackerExpectedHitsOuter() " << std::endl;
for (int i=0; i<p2.numberOfHits(); i++) {
p2.printHitPattern(i, std::cout);
}
std::cout << "\n \t trackerLayersWithMeasurement() " << p.trackerLayersWithMeasurement()
<< "\n \t pixelLayersWithMeasurement() " << p.pixelLayersWithMeasurement()
<< "\n \t stripLayersWithMeasurement() " << p.stripLayersWithMeasurement()
<< "\n \t pixelBarrelLayersWithMeasurement() " << p.pixelBarrelLayersWithMeasurement()
<< "\n \t pixelEndcapLayersWithMeasurement() " << p.pixelEndcapLayersWithMeasurement()
<< "\n \t stripTIBLayersWithMeasurement() " << p.stripTIBLayersWithMeasurement()
<< "\n \t stripTIDLayersWithMeasurement() " << p.stripTIDLayersWithMeasurement()
<< "\n \t stripTOBLayersWithMeasurement() " << p.stripTOBLayersWithMeasurement()
<< "\n \t stripTECLayersWithMeasurement() " << p.stripTECLayersWithMeasurement()
<< std::endl;
}
}
std::string IsolatedTracksNxN::algoBitToName[nL1BitsMax] [private] |
Definition at line 137 of file IsolatedTracksNxN.h.
const MagneticField* IsolatedTracksNxN::bField [private] |
Definition at line 160 of file IsolatedTracksNxN.h.
Referenced by analyze().
bool IsolatedTracksNxN::debugL1Info_ [private] |
Definition at line 142 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
int IsolatedTracksNxN::debugTrks_ [private] |
Definition at line 143 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
bool IsolatedTracksNxN::doMC [private] |
Definition at line 135 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), clearTreeVectors(), and IsolatedTracksNxN().
edm::Service<TFileService> IsolatedTracksNxN::fs [private] |
Definition at line 283 of file IsolatedTracksNxN.h.
Referenced by BookHistograms().
double IsolatedTracksNxN::genPartEtaBins[4] [private] |
Definition at line 162 of file IsolatedTracksNxN.h.
Referenced by analyze(), beginJob(), and BookHistograms().
double IsolatedTracksNxN::genPartPBins[16] [private] |
Definition at line 162 of file IsolatedTracksNxN.h.
Referenced by analyze(), beginJob(), and BookHistograms().
TH1I* IsolatedTracksNxN::h_L1AlgoNames [private] |
Definition at line 172 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and endJob().
TH1F* IsolatedTracksNxN::h_maxNearP15x15[NPBins][NEtaBins] [private] |
Definition at line 167 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_maxNearP21x21[NPBins][NEtaBins] [private] |
Definition at line 167 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_maxNearP25x25[NPBins][NEtaBins] [private] |
Definition at line 167 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_maxNearP31x31[NPBins][NEtaBins] [private] |
Definition at line 167 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F* IsolatedTracksNxN::h_nTracks [private] |
Definition at line 175 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F* IsolatedTracksNxN::h_PVTracksWt [private] |
Definition at line 173 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recEta_0 [private] |
Definition at line 177 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recEta_1 [private] |
Definition at line 180 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recEta_2 [private] |
Definition at line 183 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH2F * IsolatedTracksNxN::h_recEtaP_0 [private] |
Definition at line 178 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH2F * IsolatedTracksNxN::h_recEtaP_1 [private] |
Definition at line 181 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH2F * IsolatedTracksNxN::h_recEtaP_2 [private] |
Definition at line 184 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH2F* IsolatedTracksNxN::h_recEtaPt_0 [private] |
Definition at line 178 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH2F* IsolatedTracksNxN::h_recEtaPt_1 [private] |
Definition at line 181 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH2F* IsolatedTracksNxN::h_recEtaPt_2 [private] |
Definition at line 184 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recP_0 [private] |
Definition at line 177 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recP_1 [private] |
Definition at line 180 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recP_2 [private] |
Definition at line 183 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recPhi_0 [private] |
Definition at line 177 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recPhi_1 [private] |
Definition at line 180 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F * IsolatedTracksNxN::h_recPhi_2 [private] |
Definition at line 183 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F* IsolatedTracksNxN::h_recPt_0 [private] |
Definition at line 177 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F* IsolatedTracksNxN::h_recPt_1 [private] |
Definition at line 180 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
TH1F* IsolatedTracksNxN::h_recPt_2 [private] |
Definition at line 183 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
bool IsolatedTracksNxN::initL1 [private] |
Definition at line 135 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
edm::InputTag IsolatedTracksNxN::JetExtender_ [private] |
Definition at line 149 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
edm::InputTag IsolatedTracksNxN::JetSrc_ [private] |
Definition at line 149 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
Definition at line 146 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
Definition at line 146 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
Definition at line 147 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
Definition at line 147 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
Definition at line 147 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
Definition at line 146 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
Definition at line 148 of file IsolatedTracksNxN.h.
Referenced by IsolatedTracksNxN().
Definition at line 148 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
std::map<std::string,bool> IsolatedTracksNxN::l1TriggerMap [private] |
Definition at line 138 of file IsolatedTracksNxN.h.
double IsolatedTracksNxN::maxTrackEta_ [private] |
Definition at line 151 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
double IsolatedTracksNxN::minTrackP_ [private] |
Definition at line 151 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
int IsolatedTracksNxN::myverbose_ [private] |
Definition at line 145 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
const size_t IsolatedTracksNxN::NEtaBins = 3 [static, private] |
Definition at line 165 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
int IsolatedTracksNxN::nEventProc [private] |
Definition at line 154 of file IsolatedTracksNxN.h.
Referenced by analyze(), beginJob(), and endJob().
const size_t IsolatedTracksNxN::nL1BitsMax = 128 [static, private] |
Definition at line 136 of file IsolatedTracksNxN.h.
const size_t IsolatedTracksNxN::NPBins = 15 [static, private] |
Definition at line 164 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
bool IsolatedTracksNxN::printTrkHitPattern_ [private] |
Definition at line 144 of file IsolatedTracksNxN.h.
Referenced by IsolatedTracksNxN(), and printTrack().
double IsolatedTracksNxN::pvTracksPtMin_ [private] |
Definition at line 141 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
int IsolatedTracksNxN::t_Bunch [private] |
Definition at line 192 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
std::vector<double> * IsolatedTracksNxN::t_e11x11 [private] |
Definition at line 237 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e11x11_10Sig [private] |
Definition at line 241 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e11x11_15Sig [private] |
Definition at line 242 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e11x11_20Sig [private] |
Definition at line 243 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e11x11_25Sig [private] |
Definition at line 244 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e11x11_30Sig [private] |
Definition at line 245 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_e13x13 [private] |
Definition at line 238 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e15x15 [private] |
Definition at line 238 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e15x15_10Sig [private] |
Definition at line 241 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e15x15_15Sig [private] |
Definition at line 242 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e15x15_20Sig [private] |
Definition at line 243 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e15x15_25Sig [private] |
Definition at line 244 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e15x15_30Sig [private] |
Definition at line 245 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e21x21 [private] |
Definition at line 238 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e25x25 [private] |
Definition at line 238 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e31x31 [private] |
Definition at line 238 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_e3x3 [private] |
Definition at line 237 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e5x5 [private] |
Definition at line 237 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e7x7 [private] |
Definition at line 237 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_e7x7_10Sig [private] |
Definition at line 241 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_e7x7_15Sig [private] |
Definition at line 242 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_e7x7_20Sig [private] |
Definition at line 243 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_e7x7_25Sig [private] |
Definition at line 244 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_e7x7_30Sig [private] |
Definition at line 245 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e9x9 [private] |
Definition at line 237 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e9x9_10Sig [private] |
Definition at line 241 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e9x9_15Sig [private] |
Definition at line 242 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e9x9_20Sig [private] |
Definition at line 243 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e9x9_25Sig [private] |
Definition at line 244 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_e9x9_30Sig [private] |
Definition at line 245 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_ecalSpike11x11 [private] |
Definition at line 236 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim11x11 [private] |
Definition at line 251 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim11x11CharHad [private] |
Definition at line 266 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim11x11Matched [private] |
Definition at line 254 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim11x11NeutHad [private] |
Definition at line 263 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim11x11Photon [private] |
Definition at line 260 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim11x11Rest [private] |
Definition at line 257 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim13x13 [private] |
Definition at line 252 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim13x13CharHad [private] |
Definition at line 267 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim13x13Matched [private] |
Definition at line 255 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim13x13NeutHad [private] |
Definition at line 264 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim13x13Photon [private] |
Definition at line 261 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim13x13Rest [private] |
Definition at line 258 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim15x15 [private] |
Definition at line 252 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim15x15CharHad [private] |
Definition at line 267 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim15x15Matched [private] |
Definition at line 255 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim15x15NeutHad [private] |
Definition at line 264 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim15x15Photon [private] |
Definition at line 261 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim15x15Rest [private] |
Definition at line 258 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim21x21 [private] |
Definition at line 252 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim21x21CharHad [private] |
Definition at line 267 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim21x21Matched [private] |
Definition at line 255 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim21x21NeutHad [private] |
Definition at line 264 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim21x21Photon [private] |
Definition at line 261 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim21x21Rest [private] |
Definition at line 258 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim25x25 [private] |
Definition at line 252 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim25x25CharHad [private] |
Definition at line 267 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim25x25Matched [private] |
Definition at line 255 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim25x25NeutHad [private] |
Definition at line 264 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim25x25Photon [private] |
Definition at line 261 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim25x25Rest [private] |
Definition at line 258 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim31x31 [private] |
Definition at line 252 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim31x31CharHad [private] |
Definition at line 267 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim31x31Matched [private] |
Definition at line 255 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim31x31NeutHad [private] |
Definition at line 264 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim31x31Photon [private] |
Definition at line 261 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim31x31Rest [private] |
Definition at line 258 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim3x3 [private] |
Definition at line 251 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim3x3CharHad [private] |
Definition at line 266 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim3x3Matched [private] |
Definition at line 254 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim3x3NeutHad [private] |
Definition at line 263 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim3x3Photon [private] |
Definition at line 260 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esim3x3Rest [private] |
Definition at line 257 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim5x5 [private] |
Definition at line 251 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim5x5CharHad [private] |
Definition at line 266 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim5x5Matched [private] |
Definition at line 254 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim5x5NeutHad [private] |
Definition at line 263 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim5x5Photon [private] |
Definition at line 260 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim5x5Rest [private] |
Definition at line 257 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim7x7 [private] |
Definition at line 251 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim7x7CharHad [private] |
Definition at line 266 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim7x7Matched [private] |
Definition at line 254 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim7x7NeutHad [private] |
Definition at line 263 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim7x7Photon [private] |
Definition at line 260 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim7x7Rest [private] |
Definition at line 257 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim9x9 [private] |
Definition at line 251 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim9x9CharHad [private] |
Definition at line 266 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim9x9Matched [private] |
Definition at line 254 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim9x9NeutHad [private] |
Definition at line 263 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim9x9Photon [private] |
Definition at line 260 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_esim9x9Rest [private] |
Definition at line 257 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_esimPdgId [private] |
Definition at line 247 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
int IsolatedTracksNxN::t_EvtNo [private] |
Definition at line 192 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
std::vector<double>* IsolatedTracksNxN::t_h3x3 [private] |
Definition at line 271 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_h3x3Sig [private] |
Definition at line 272 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_h5x5 [private] |
Definition at line 271 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_h5x5Sig [private] |
Definition at line 272 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_h7x7 [private] |
Definition at line 271 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_h7x7Sig [private] |
Definition at line 272 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_hsim3x3 [private] |
Definition at line 276 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_hsim3x3CharHad [private] |
Definition at line 281 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_hsim3x3Matched [private] |
Definition at line 277 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_hsim3x3NeutHad [private] |
Definition at line 280 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_hsim3x3Photon [private] |
Definition at line 279 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_hsim3x3Rest [private] |
Definition at line 278 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim5x5 [private] |
Definition at line 276 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim5x5CharHad [private] |
Definition at line 281 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim5x5Matched [private] |
Definition at line 277 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim5x5NeutHad [private] |
Definition at line 280 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim5x5Photon [private] |
Definition at line 279 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim5x5Rest [private] |
Definition at line 278 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim7x7 [private] |
Definition at line 276 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim7x7CharHad [private] |
Definition at line 281 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim7x7Matched [private] |
Definition at line 277 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim7x7NeutHad [private] |
Definition at line 280 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim7x7Photon [private] |
Definition at line 279 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_hsim7x7Rest [private] |
Definition at line 278 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_infoHcal [private] |
Definition at line 273 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_jetEta [private] |
Definition at line 205 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_jetPhi [private] |
Definition at line 205 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_jetPt [private] |
Definition at line 205 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1CenJetEta [private] |
Definition at line 197 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1CenJetPhi [private] |
Definition at line 197 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_L1CenJetPt [private] |
Definition at line 197 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_L1Decision [private] |
Definition at line 196 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1FwdJetEta [private] |
Definition at line 198 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1FwdJetPhi [private] |
Definition at line 198 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_L1FwdJetPt [private] |
Definition at line 198 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1IsoEMEta [private] |
Definition at line 201 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1IsoEMPhi [private] |
Definition at line 201 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_L1IsoEMPt [private] |
Definition at line 201 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1METEta [private] |
Definition at line 203 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1METPhi [private] |
Definition at line 203 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_L1METPt [private] |
Definition at line 203 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1MuonEta [private] |
Definition at line 200 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1MuonPhi [private] |
Definition at line 200 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_L1MuonPt [private] |
Definition at line 200 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1NonIsoEMEta [private] |
Definition at line 202 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1NonIsoEMPhi [private] |
Definition at line 202 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_L1NonIsoEMPt [private] |
Definition at line 202 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1TauJetEta [private] |
Definition at line 199 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_L1TauJetPhi [private] |
Definition at line 199 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_L1TauJetPt [private] |
Definition at line 199 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
int IsolatedTracksNxN::t_Lumi [private] |
Definition at line 192 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
std::vector<double>* IsolatedTracksNxN::t_maxNearHcalP3x3 [private] |
Definition at line 270 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_maxNearHcalP5x5 [private] |
Definition at line 270 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_maxNearHcalP7x7 [private] |
Definition at line 270 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP11x11 [private] |
Definition at line 232 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP13x13 [private] |
Definition at line 231 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP15x15 [private] |
Definition at line 230 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP21x21 [private] |
Definition at line 229 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP25x25 [private] |
Definition at line 228 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP31x31 [private] |
Definition at line 227 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP7x7 [private] |
Definition at line 234 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_maxNearP9x9 [private] |
Definition at line 233 of file IsolatedTracksNxN.h.
Referenced by BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_NLayersCrossed [private] |
Definition at line 219 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
int IsolatedTracksNxN::t_nTracks [private] |
Definition at line 190 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
std::vector<double>* IsolatedTracksNxN::t_nTrksJetCalo [private] |
Definition at line 206 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_nTrksJetVtx [private] |
Definition at line 206 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVisValid [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVndof [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVNTracks [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_PVNTracksHP [private] |
Definition at line 194 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVNTracksHPWt [private] |
Definition at line 194 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVNTracksWt [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVTracksSumPt [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVTracksSumPtHP [private] |
Definition at line 194 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVTracksSumPtHPWt [private] |
Definition at line 194 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVTracksSumPtWt [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_PVx [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVy [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_PVz [private] |
Definition at line 193 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
int IsolatedTracksNxN::t_RunNo [private] |
Definition at line 192 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
std::vector<double> * IsolatedTracksNxN::t_simTrackP [private] |
Definition at line 247 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackChiSq [private] |
Definition at line 216 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackChiSqAll [private] |
Definition at line 210 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackDxy [private] |
Definition at line 214 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackDxyAll [private] |
Definition at line 210 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackDxyBS [private] |
Definition at line 214 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackDxyPV [private] |
Definition at line 215 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackDxyPVAll [private] |
Definition at line 210 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackDz [private] |
Definition at line 214 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackDzAll [private] |
Definition at line 210 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackDzBS [private] |
Definition at line 214 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackDzPV [private] |
Definition at line 215 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackDzPVAll [private] |
Definition at line 210 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackEcalEta [private] |
Definition at line 213 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackEcalPhi [private] |
Definition at line 213 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackEta [private] |
Definition at line 212 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackEtaAll [private] |
Definition at line 208 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackHcalEta [private] |
Definition at line 213 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackHcalPhi [private] |
Definition at line 213 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitInMeasTEC [private] |
Definition at line 223 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitInMeasTIB [private] |
Definition at line 223 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitInMeasTID [private] |
Definition at line 223 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_trackHitInMeasTOB [private] |
Definition at line 223 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTEC [private] |
Definition at line 221 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTIB [private] |
Definition at line 221 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTIBTID [private] |
Definition at line 221 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTID [private] |
Definition at line 221 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_trackHitInMissTOB [private] |
Definition at line 221 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitOutMeasTEC [private] |
Definition at line 224 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitOutMeasTIB [private] |
Definition at line 224 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitOutMeasTID [private] |
Definition at line 224 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_trackHitOutMeasTOB [private] |
Definition at line 224 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTEC [private] |
Definition at line 222 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTIB [private] |
Definition at line 222 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTID [private] |
Definition at line 222 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_trackHitOutMissTOB [private] |
Definition at line 222 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTOBTEC [private] |
Definition at line 222 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackHitsTEC [private] |
Definition at line 220 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_trackHitsTOB [private] |
Definition at line 220 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackL [private] |
Definition at line 225 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int> * IsolatedTracksNxN::t_trackNOuterHits [private] |
Definition at line 219 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackOutPosOutHitDr [private] |
Definition at line 225 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackP [private] |
Definition at line 212 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackPAll [private] |
Definition at line 208 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackPdgIdAll [private] |
Definition at line 208 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackPhi [private] |
Definition at line 212 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackPhiAll [private] |
Definition at line 208 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double> * IsolatedTracksNxN::t_trackPt [private] |
Definition at line 212 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trackPtAll [private] |
Definition at line 209 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<int>* IsolatedTracksNxN::t_trackPVIdx [private] |
Definition at line 217 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trkEcalEne [private] |
Definition at line 249 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
std::vector<double>* IsolatedTracksNxN::t_trkHcalEne [private] |
Definition at line 275 of file IsolatedTracksNxN.h.
Referenced by analyze(), BookHistograms(), and clearTreeVectors().
double IsolatedTracksNxN::tMaxE_ [private] |
Definition at line 152 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
double IsolatedTracksNxN::tMaxH_ [private] |
Definition at line 152 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
double IsolatedTracksNxN::tMinE_ [private] |
Definition at line 152 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
double IsolatedTracksNxN::tMinH_ [private] |
Definition at line 152 of file IsolatedTracksNxN.h.
Referenced by analyze(), and IsolatedTracksNxN().
TTree* IsolatedTracksNxN::tree [private] |
Definition at line 188 of file IsolatedTracksNxN.h.
Referenced by analyze(), and BookHistograms().
std::map<std::string,bool>::iterator IsolatedTracksNxN::trig_iter [private] |
Definition at line 139 of file IsolatedTracksNxN.h.
1.7.3