| 54 |
|
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h" |
| 55 |
|
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h" |
| 56 |
|
|
| 57 |
< |
|
| 57 |
> |
#include "DataFormats/CaloRecHit/interface/CaloRecHit.h" |
| 58 |
> |
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h" |
| 59 |
> |
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h" |
| 60 |
|
|
| 61 |
|
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h" |
| 62 |
|
#include "SimDataFormats/TrackingHit/interface/PSimHit.h" |
| 249 |
|
OpenGLDisplayer::endJob() { |
| 250 |
|
|
| 251 |
|
MyEvents->Save((char*) OutputFile.c_str()); |
| 252 |
< |
MyEvents->Load((char*) OutputFile.c_str()); |
| 252 |
> |
// MyEvents->Load((char*) OutputFile.c_str()); |
| 253 |
|
} |
| 254 |
|
|
| 255 |
|
|
| 371 |
|
MyEvent->MyRecoTrackCollection.push_back(MyrecoTrack); |
| 372 |
|
} |
| 373 |
|
|
| 374 |
+ |
|
| 375 |
+ |
|
| 376 |
+ |
edm::Handle<EcalRecHitCollection > h_EcalEB_RecHits; |
| 377 |
+ |
iEvent.getByLabel("ecalRecHit","EcalRecHitsEB", h_EcalEB_RecHits); |
| 378 |
+ |
EcalRecHitCollection EcalEB_RecHits = *h_EcalEB_RecHits.product(); |
| 379 |
+ |
|
| 380 |
+ |
for(unsigned int eh=0;eh<EcalEB_RecHits.size();eh++){ |
| 381 |
+ |
MyCaloHit temp_EcalEBHit; |
| 382 |
+ |
temp_EcalEBHit.E = EcalEB_RecHits[eh].energy(); |
| 383 |
+ |
temp_EcalEBHit.t = EcalEB_RecHits[eh].time(); |
| 384 |
+ |
temp_EcalEBHit.DetId = (EcalEB_RecHits[eh].detid()).rawId(); |
| 385 |
+ |
|
| 386 |
+ |
MyEvent->MyEcalCaloHitCollection.push_back(temp_EcalEBHit); |
| 387 |
+ |
} |
| 388 |
+ |
|
| 389 |
+ |
|
| 390 |
+ |
edm::Handle<EcalRecHitCollection > h_EcalEE_RecHits; |
| 391 |
+ |
iEvent.getByLabel("ecalRecHit","EcalRecHitsEE", h_EcalEE_RecHits); |
| 392 |
+ |
EcalRecHitCollection EcalEE_RecHits = *h_EcalEE_RecHits.product(); |
| 393 |
+ |
|
| 394 |
+ |
for(unsigned int eh=0;eh<EcalEE_RecHits.size();eh++){ |
| 395 |
+ |
MyCaloHit temp_EcalEEHit; |
| 396 |
+ |
temp_EcalEEHit.E = EcalEE_RecHits[eh].energy(); |
| 397 |
+ |
temp_EcalEEHit.t = EcalEE_RecHits[eh].time(); |
| 398 |
+ |
temp_EcalEEHit.DetId = (EcalEE_RecHits[eh].detid()).rawId(); |
| 399 |
+ |
|
| 400 |
+ |
MyEvent->MyEcalCaloHitCollection.push_back(temp_EcalEEHit); |
| 401 |
+ |
} |
| 402 |
+ |
|
| 403 |
+ |
|
| 404 |
+ |
edm::Handle<EcalRecHitCollection > h_EcalES_RecHits; |
| 405 |
+ |
iEvent.getByLabel("ecalPreshowerRecHit","EcalRecHitsES", h_EcalES_RecHits); |
| 406 |
+ |
EcalRecHitCollection EcalES_RecHits = *h_EcalES_RecHits.product(); |
| 407 |
+ |
|
| 408 |
+ |
for(unsigned int eh=0;eh<EcalES_RecHits.size();eh++){ |
| 409 |
+ |
MyCaloHit temp_EcalESHit; |
| 410 |
+ |
temp_EcalESHit.E = EcalES_RecHits[eh].energy(); |
| 411 |
+ |
temp_EcalESHit.t = EcalES_RecHits[eh].time(); |
| 412 |
+ |
temp_EcalESHit.DetId = (EcalES_RecHits[eh].detid()).rawId(); |
| 413 |
+ |
|
| 414 |
+ |
MyEvent->MyEcalCaloHitCollection.push_back(temp_EcalESHit); |
| 415 |
+ |
} |
| 416 |
+ |
|
| 417 |
+ |
|
| 418 |
+ |
|
| 419 |
+ |
edm::Handle<HBHERecHitCollection > h_HcalHBHE_RecHits; |
| 420 |
+ |
iEvent.getByLabel("hbhereco", h_HcalHBHE_RecHits); |
| 421 |
+ |
HBHERecHitCollection HcalHBHE_RecHits = *h_HcalHBHE_RecHits.product(); |
| 422 |
+ |
|
| 423 |
+ |
for(unsigned int hh=0;hh<HcalHBHE_RecHits.size();hh++){ |
| 424 |
+ |
MyCaloHit temp_HcalHBHEHit; |
| 425 |
+ |
temp_HcalHBHEHit.E = HcalHBHE_RecHits[hh].energy(); |
| 426 |
+ |
temp_HcalHBHEHit.t = HcalHBHE_RecHits[hh].time(); |
| 427 |
+ |
temp_HcalHBHEHit.DetId = (HcalHBHE_RecHits[hh].detid()).rawId(); |
| 428 |
+ |
|
| 429 |
+ |
MyEvent->MyHcalCaloHitCollection.push_back(temp_HcalHBHEHit); |
| 430 |
+ |
} |
| 431 |
+ |
|
| 432 |
+ |
|
| 433 |
+ |
edm::Handle<HORecHitCollection > h_HcalHO_RecHits; |
| 434 |
+ |
iEvent.getByLabel("horeco", h_HcalHO_RecHits); |
| 435 |
+ |
HORecHitCollection HcalHO_RecHits = *h_HcalHO_RecHits.product(); |
| 436 |
+ |
|
| 437 |
+ |
for(unsigned int hh=0;hh<HcalHO_RecHits.size();hh++){ |
| 438 |
+ |
MyCaloHit temp_HcalHOHit; |
| 439 |
+ |
temp_HcalHOHit.E = HcalHO_RecHits[hh].energy(); |
| 440 |
+ |
temp_HcalHOHit.t = HcalHO_RecHits[hh].time(); |
| 441 |
+ |
temp_HcalHOHit.DetId = (HcalHO_RecHits[hh].detid()).rawId(); |
| 442 |
+ |
|
| 443 |
+ |
MyEvent->MyHcalCaloHitCollection.push_back(temp_HcalHOHit); |
| 444 |
+ |
} |
| 445 |
+ |
|
| 446 |
+ |
edm::Handle<HFRecHitCollection > h_HcalHF_RecHits; |
| 447 |
+ |
iEvent.getByLabel("hfreco", h_HcalHF_RecHits); |
| 448 |
+ |
HFRecHitCollection HcalHF_RecHits = *h_HcalHF_RecHits.product(); |
| 449 |
+ |
|
| 450 |
+ |
for(unsigned int hh=0;hh<HcalHF_RecHits.size();hh++){ |
| 451 |
+ |
MyCaloHit temp_HcalHFHit; |
| 452 |
+ |
temp_HcalHFHit.E = HcalHF_RecHits[hh].energy(); |
| 453 |
+ |
temp_HcalHFHit.t = HcalHF_RecHits[hh].time(); |
| 454 |
+ |
temp_HcalHFHit.DetId = (HcalHF_RecHits[hh].detid()).rawId(); |
| 455 |
+ |
|
| 456 |
+ |
MyEvent->MyHcalCaloHitCollection.push_back(temp_HcalHFHit); |
| 457 |
+ |
} |
| 458 |
+ |
|
| 459 |
+ |
|
| 460 |
+ |
|
| 461 |
+ |
|
| 462 |
+ |
|
| 463 |
+ |
|
| 464 |
+ |
|
| 465 |
+ |
|
| 466 |
+ |
|
| 467 |
+ |
|
| 468 |
+ |
|
| 469 |
|
MyEvents->Events.push_back(MyEvent); |
| 470 |
|
|
| 471 |
|
} |
