362 float rnnum = (float)iEvent.
run();
383 if (!gtRecord.isValid()) {
390 cout<<
" No L1 trigger record "<<endl;
395 for (
CItAlgo itAlgo = bitMap.begin(); itAlgo != bitMap.end(); itAlgo++)
397 bool decision=
menu->gtAlgorithmResult(itAlgo->first,dWord);
398 if(decision == 1)
std::cout<<
" Trigger "<<itAlgo->first<<
" "<<decision<<std::endl;
408 myRecalib = recalibCorrs.
product();
414 double tmpNSFillDetMapPl1[5][5][73][43];
415 double tmpNSFillDetMapMin1[5][5][73][43];
417 for (
int i=0;
i<5;
i++)
419 for (
int j=0;
j<5;
j++)
421 for (
int k=0;
k<73;
k++)
423 for (
int l=0;
l<43;
l++)
425 tmpNSFillDetMapPl1[
i][
j][
k][
l] = 0.;
426 tmpNSFillDetMapMin1[
i][
j][
k][
l] = 0.;
433 if(!hbheNormal.isValid()){
434 cout<<
" hbheNormal failed "<<endl;
436 cout<<
" The size of the normal collection "<<hbheNormal->size()<<endl;
444 if(!hbheNS.isValid()){
445 LogDebug(
"") <<
"HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
446 cout<<
" No HBHE MS "<<endl;
452 cout<<
" HBHE NS size of collection "<<HithbheNS.size()<<endl;
455 if(HithbheNS.size()!= 5184) {
456 cout<<
" HBHE problem "<<rnnum<<
" "<<HithbheNS.size()<<endl;
462 if(!hbheMB.isValid()){
463 LogDebug(
"") <<
"HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
464 cout<<
" No HBHE MB"<<endl;
469 cout<<
" HBHE MB size of collection "<<HithbheMB.size()<<endl;
470 if(HithbheMB.size()!= 5184) {
471 cout<<
" HBHE problem "<<rnnum<<
" "<<HithbheMB.size()<<endl;
479 LogDebug(
"") <<
"HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
480 cout<<
" No HF NS "<<endl;
485 cout<<
" HFE NS size of collection "<<HithfNS.size()<<endl;
487 if(HithfNS.size()!= 1728) {
488 cout<<
" HF problem "<<rnnum<<
" "<<HithfNS.size()<<endl;
496 LogDebug(
"") <<
"HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
497 cout<<
" No HBHE MB"<<endl;
502 cout<<
" HF MB size of collection "<<HithfMB.size()<<endl;
503 if(HithfMB.size()!= 1728) {
504 cout<<
" HF problem "<<rnnum<<
" "<<HithfMB.size()<<endl;
515 if(
theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
517 HBHERecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
519 double energyhit = aHit.
energy();
521 DetId id = (*hbheItr).detid();
526 int mysu = ((hid).rawId()>>25)&0x7;
527 if( hid.
ieta() > 0 ) {
533 tmpNSFillDetMapPl1[mysu][hid.
depth()][hid.
iphi()][hid.
ieta()] = energyhit;
543 tmpNSFillDetMapMin1[mysu][hid.
depth()][hid.
iphi()][hid.
ieta()] = energyhit;
547 if(hid.
depth() == 1) {
550 if(energyhit<-2.)
std::cout<<
" Run "<<rnnum<<
" ieta,iphi "<<hid.
ieta()<<
" "<<hid.
iphi()<<energyhit<<std::endl;
574 if(
theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
576 HBHERecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
578 double energyhit = aHit.
energy();
580 DetId id = (*hbheItr).detid();
583 int mysu = ((hid).rawId()>>25)&0x7;
584 if( hid.
ieta() > 0 ) {
589 float mydiff = energyhit - tmpNSFillDetMapPl1[mysu][hid.
depth()][hid.
iphi()][hid.
ieta()];
601 float mydiff = energyhit - tmpNSFillDetMapMin1[mysu][hid.
depth()][hid.
iphi()][hid.
ieta()];
607 if(hid.
depth() == 1) {
611 if( hid.
ieta() > 0 ) {
631 if(
theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
633 HFRecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
635 double energyhit = aHit.
energy();
639 if(fabs(energyhit) > 40. )
continue;
641 DetId id = (*hbheItr).detid();
644 int mysu = ((hid).rawId()>>25)&0x7;
645 if( hid.
ieta() > 0 ) {
651 tmpNSFillDetMapPl1[mysu][hid.
depth()][hid.
iphi()][hid.
ieta()] = energyhit;
661 tmpNSFillDetMapMin1[mysu][hid.
depth()][hid.
iphi()][hid.
ieta()] = energyhit;
665 if(hid.
depth() == 1) {
688 if(
theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
690 HFRecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
692 double energyhit = aHit.
energy();
696 if(fabs(energyhit) > 40. )
continue;
698 DetId id = (*hbheItr).detid();
701 int mysu = ((hid).rawId()>>25)&0x7;
702 if( hid.
ieta() > 0 ) {
724 if(hid.
depth() == 1) {
741 std::cout<<
" Event is finished "<<std::endl;
double theMBFillDetMapMin0[5][5][73][43]
double theMBFillDetMapMin2[5][5][73][43]
double theDFFillDetMapPl2[5][5][73][43]
double theMBFillDetMapMin1[5][5][73][43]
double theNSFillDetMapPl4[5][5][73][43]
double theNSFillDetMapPl0[5][5][73][43]
double theNSFillDetMapMin0[5][5][73][43]
double theMBFillDetMapPl0[5][5][73][43]
bool getByToken(EDGetToken token, Handle< PROD > &result) const
TH1F * hCalo2mom2[73][43]
std::vector< HBHERecHit >::const_iterator const_iterator
double theDFFillDetMapPl1[5][5][73][43]
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > tok_gtRec_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNoise_
double theNSFillDetMapPl2[5][5][73][43]
double theDFFillDetMapMin2[5][5][73][43]
TH1F * hCalo1mom2[73][43]
edm::EDGetTokenT< HFRecHitCollection > tok_hfNoise_
double theNSFillDetMapPl1[5][5][73][43]
uint32_t rawId() const
get the raw id
double theMBFillDetMapPl1[5][5][73][43]
int depth() const
get the tower depth
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
std::vector< bool > DecisionWord
typedefs
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
int ieta() const
get the cell ieta
Abs< T >::type abs(const T &t)
double theMBFillDetMapPl4[5][5][73][43]
double theNSFillDetMapMin1[5][5][73][43]
double theNSFillDetMapMin4[5][5][73][43]
double theMBFillDetMapPl2[5][5][73][43]
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNorm_
int iphi() const
get the cell iphi
T const * product() const
double theMBFillDetMapMin4[5][5][73][43]
return(e1-e2)*(e1-e2)+dp *dp
double theNSFillDetMapMin2[5][5][73][43]
double theDFFillDetMapMin1[5][5][73][43]
Power< A, B >::type pow(const A &a, const B &b)