| 48 |
|
#include "FWCore/Utilities/interface/InputTag.h" |
| 49 |
|
#include "DataFormats/HeavyIonEvent/interface/CentralityBins.h" |
| 50 |
|
#include "DataFormats/CaloTowers/interface/CaloTower.h" |
| 51 |
< |
#include "DataFormats/HeavyIonEvent/interface/Centrality.h" |
| 51 |
> |
#include "DataFormats/HeavyIonEvent/interface/CentralityProvider.h" |
| 52 |
|
#include "SimDataFormats/HiGenData/interface/GenHIEvent.h" |
| 53 |
|
#include "DataFormats/PatCandidates/interface/Jet.h" |
| 54 |
|
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h" |
| 75 |
|
virtual void endJob() ; |
| 76 |
|
|
| 77 |
|
// ----------member data --------------------------- |
| 78 |
– |
|
| 79 |
– |
|
| 80 |
– |
edm::Handle<reco::Centrality> cent; |
| 78 |
|
edm::Handle<vector<double> > ktRhos; |
| 79 |
|
edm::Handle<vector<double> > akRhos; |
| 80 |
|
|
| 103 |
|
edm::InputTag EESrc1_; |
| 104 |
|
edm::InputTag EESrc2_; |
| 105 |
|
edm::InputTag signalTag_; |
| 109 |
– |
edm::InputTag centTag_; |
| 106 |
|
|
| 107 |
|
TNtuple* ntEB; |
| 108 |
|
TNtuple* ntEE; |
| 111 |
|
TNtuple* ntjet; |
| 112 |
|
|
| 113 |
|
double cone; |
| 114 |
+ |
bool jetsOnly_; |
| 115 |
|
|
| 116 |
|
edm::Service<TFileService> fs; |
| 117 |
< |
const CentralityBins * cbins_; |
| 117 |
> |
CentralityProvider* centrality_; |
| 118 |
|
const CaloGeometry *geo; |
| 119 |
|
}; |
| 120 |
|
|
| 130 |
|
// constructors and destructor |
| 131 |
|
// |
| 132 |
|
RecHitComparison::RecHitComparison(const edm::ParameterSet& iConfig) : |
| 133 |
< |
cbins_(0), |
| 133 |
> |
centrality_(0), |
| 134 |
|
geo(0), |
| 135 |
|
cone(0.5) |
| 136 |
|
{ |
| 137 |
|
//now do what ever initialization is needed |
| 138 |
< |
centTag_ = iConfig.getUntrackedParameter<edm::InputTag>("centrality",edm::InputTag("hiCentrality","","RECO")); |
| 142 |
< |
|
| 138 |
> |
jetsOnly_ = iConfig.getUntrackedParameter<bool>("jetsOnly",false); |
| 139 |
|
signalTag_ = iConfig.getUntrackedParameter<edm::InputTag>("signalJets",edm::InputTag("iterativeCone5CaloJets","","SIGNAL")); |
| 140 |
|
|
| 141 |
|
HcalRecHitHFSrc1_ = iConfig.getUntrackedParameter<edm::InputTag>("hcalHFRecHitSrc1",edm::InputTag("hfreco")); |
| 167 |
|
void |
| 168 |
|
RecHitComparison::analyze(const edm::Event& ev, const edm::EventSetup& iSetup) |
| 169 |
|
{ |
| 170 |
< |
|
| 175 |
< |
if(!cbins_) cbins_ = getCentralityBinsFromDB(iSetup); |
| 170 |
> |
if(!centrality_) centrality_ = new CentralityProvider(iSetup); |
| 171 |
|
if(!geo){ |
| 172 |
|
edm::ESHandle<CaloGeometry> pGeo; |
| 173 |
|
iSetup.get<CaloGeometryRecord>().get(pGeo); |
| 175 |
|
} |
| 176 |
|
|
| 177 |
|
using namespace edm; |
| 183 |
– |
|
| 184 |
– |
ev.getByLabel(centTag_,cent); |
| 178 |
|
ev.getByLabel(EBSrc1_,ebHits1); |
| 179 |
< |
ev.getByLabel(EBSrc1_,ebHits2); |
| 179 |
> |
ev.getByLabel(EBSrc2_,ebHits2); |
| 180 |
|
ev.getByLabel(signalTag_,signalJets); |
| 181 |
|
|
| 182 |
|
ev.getByLabel(HcalRecHitHFSrc1_,hfHits1); |
| 187 |
|
ev.getByLabel(EESrc1_,eeHits1); |
| 188 |
|
ev.getByLabel(EESrc2_,eeHits2); |
| 189 |
|
|
| 190 |
< |
double hf = cent->EtHFhitSum(); |
| 191 |
< |
double sumet = cent->EtMidRapiditySum(); |
| 192 |
< |
int run = ev.id().run(); |
| 200 |
< |
run = 1; |
| 201 |
< |
int bin = cbins_->getBin(hf); |
| 202 |
< |
int margin = 0; |
| 190 |
> |
centrality_->newEvent(ev,iSetup); |
| 191 |
> |
double hf = centrality_->centralityValue(ev); |
| 192 |
> |
int bin = centrality_->getBin(ev); |
| 193 |
|
|
| 194 |
|
vector<double> fFull; |
| 195 |
|
vector<double> f05; |
| 278 |
|
double eta2 = pos2.eta(); |
| 279 |
|
double phi2 = pos2.eta(); |
| 280 |
|
double et2 = e2*sin(pos2.theta()); |
| 281 |
< |
if(isjet) ntEB->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet); |
| 281 |
> |
if(!jetsOnly_ || isjet) ntEB->Fill(e1,et1,e2,et2,eta2,phi2,hf,bin,jetpt,drjet); |
| 282 |
|
} |
| 283 |
|
|
| 284 |
|
for(unsigned int i = 0; i < eeHits1->size(); ++i){ |
| 315 |
|
double eta2 = pos2.eta(); |
| 316 |
|
double phi2 = pos2.eta(); |
| 317 |
|
double et2 = e2*sin(pos2.theta()); |
| 318 |
< |
if(isjet) ntEE->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet); |
| 318 |
> |
if(!jetsOnly_ || isjet) ntEE->Fill(e1,et1,e2,et2,eta2,phi2,hf,bin,jetpt,drjet); |
| 319 |
|
} |
| 320 |
|
|
| 321 |
|
for(unsigned int i = 0; i < hbheHits1->size(); ++i){ |
| 352 |
|
double eta2 = pos2.eta(); |
| 353 |
|
double phi2 = pos2.eta(); |
| 354 |
|
double et2 = e2*sin(pos2.theta()); |
| 355 |
< |
if(isjet) ntHBHE->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet); |
| 355 |
> |
if(!jetsOnly_ || isjet) ntHBHE->Fill(e1,et1,e2,et2,eta2,phi2,hf,bin,jetpt,drjet); |
| 356 |
|
} |
| 357 |
|
|
| 358 |
|
for(unsigned int i = 0; i < hfHits1->size(); ++i){ |
| 389 |
|
double eta2 = pos2.eta(); |
| 390 |
|
double phi2 = pos2.eta(); |
| 391 |
|
double et2 = e2*sin(pos2.theta()); |
| 392 |
< |
if(isjet) ntHF->Fill(e1,et1,e2,et2,eta2,phi2,sumet,hf,bin,jetpt,drjet); |
| 392 |
> |
if(!jetsOnly_ || isjet) ntHF->Fill(e1,et1,e2,et2,eta2,phi2,hf,bin,jetpt,drjet); |
| 393 |
|
} |
| 394 |
|
|
| 395 |
|
for(unsigned int j1 = 0 ; j1 < signalJets->size(); ++j1){ |
| 408 |
|
void |
| 409 |
|
RecHitComparison::beginJob() |
| 410 |
|
{ |
| 411 |
< |
ntEB = fs->make<TNtuple>("ntEB","","e1:et1:e2:et2:eta:phi:sumet:hf:bin:ptjet:drjet"); |
| 412 |
< |
ntEE = fs->make<TNtuple>("ntEE","","e1:et1:e2:et2:eta:phi:sumet:hf:bin:ptjet:drjet"); |
| 413 |
< |
ntHBHE = fs->make<TNtuple>("ntHBHE","","e1:et1:e2:et2:eta:phi:sumet:hf:bin:ptjet:drjet"); |
| 414 |
< |
ntHF = fs->make<TNtuple>("ntHF","","e1:et1:e2:et2:eta:phi:sumet:hf:bin:ptjet:drjet"); |
| 411 |
> |
ntEB = fs->make<TNtuple>("ntEB","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet"); |
| 412 |
> |
ntEE = fs->make<TNtuple>("ntEE","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet"); |
| 413 |
> |
ntHBHE = fs->make<TNtuple>("ntHBHE","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet"); |
| 414 |
> |
ntHF = fs->make<TNtuple>("ntHF","","e1:et1:e2:et2:eta:phi:hf:bin:ptjet:drjet"); |
| 415 |
|
ntjet = fs->make<TNtuple>("ntjet","","bin:pt:eta:ethit:f05:f1:f15:f2:f25:f3:em:emf"); |
| 416 |
|
|
| 417 |
|
} |
