AnalyzerMinbias Class Reference

Inheritance diagram for AnalyzerMinbias:

Classes
struct	myInfo

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override
	AnalyzerMinbias (const edm::ParameterSet &)
virtual void	beginJob ()
virtual void	endJob ()
	~AnalyzerMinbias ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
void	analyzeHcal (const HcalRespCorrs *myRecalib, const HBHERecHitCollection &HithbheNS, const HBHERecHitCollection &HithbheMB, const HFRecHitCollection &HithfNS, const HFRecHitCollection &HithfMB, int algoBit, bool fill)

Private Attributes
int	cells
int	depth
float	eta
std::string	fOutputFileName
TH1D *	h_Noise [4]
TH1D *	h_Signal [4]
std::string	hcalfile_
TFile *	hOutputFile
int	ieta
bool	ignoreL1_
int	iphi
float	mom0_Diff
float	mom0_MB
float	mom0_Noise
float	mom1_Diff
float	mom1_MB
float	mom1_Noise
float	mom2_Diff
float	mom2_MB
float	mom2_Noise
float	mom3_Diff
float	mom3_MB
float	mom3_Noise
float	mom4_Diff
float	mom4_MB
float	mom4_Noise
int	mydet
std::map< std::pair< int, HcalDetId >, myInfo >	myMap_
std::ofstream *	myout_hcal
int	mysubd
TTree *	myTree
float	occup
float	phi
double	rnnum
double	rnnumber
bool	runNZS_
bool	theRecalib_
edm::EDGetTokenT < HBHERecHitCollection >	tok_hbheNormal_
edm::EDGetTokenT < HBHERecHitCollection >	tok_hbherecoMB_
edm::EDGetTokenT < HBHERecHitCollection >	tok_hbherecoNoise_
edm::EDGetTokenT < HFRecHitCollection >	tok_hfrecoMB_
edm::EDGetTokenT < HFRecHitCollection >	tok_hfrecoNoise_
edm::EDGetTokenT < L1GlobalTriggerObjectMapRecord >	tok_hltL1GtMap_
edm::EDGetTokenT < HORecHitCollection >	tok_horecoMB_
edm::EDGetTokenT < HORecHitCollection >	tok_horecoNoise_
int	trigbit

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 66 of file AnalyzerMinbias.cc.

Constructor & Destructor Documentation

AnalyzerMinbias::AnalyzerMinbias ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 120 of file AnalyzerMinbias.cc.

References fOutputFileName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ignoreL1_, HLT_FULL_cff::InputTag, runNZS_, AlCaHLTBitMon_QueryRunRegistry::string, theRecalib_, tok_hbheNormal_, tok_hbherecoMB_, tok_hbherecoNoise_, tok_hfrecoMB_, tok_hfrecoNoise_, tok_hltL1GtMap_, tok_horecoMB_, and tok_horecoNoise_.

                                                               {

  // get name of output file with histogramms
  fOutputFileName = iConfig.getUntrackedParameter<std::string>("HistOutFile");
  
  // get token names of modules, producing object collections
  tok_hbherecoMB_   = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"));
  tok_horecoMB_     = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputMB"));
  tok_hfrecoMB_     = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"));

  tok_hbherecoNoise_= consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputNoise"));
  tok_horecoNoise_  = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputNoise"));
  tok_hfrecoNoise_  = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputNoise"));

  theRecalib_       = iConfig.getParameter<bool>("Recalib"); 
  ignoreL1_         = iConfig.getUntrackedParameter<bool>("IgnoreL1", true);
  runNZS_           = iConfig.getUntrackedParameter<bool>("RunNZS", true);


  tok_hbheNormal_   = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
  tok_hltL1GtMap_   = consumes<L1GlobalTriggerObjectMapRecord>(edm::InputTag("hltL1GtObjectMap"));
}

AnalyzerMinbias::~AnalyzerMinbias

(

)

Definition at line 143 of file AnalyzerMinbias.cc.

{ }

Member Function Documentation

void AnalyzerMinbias::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 250 of file AnalyzerMinbias.cc.

References analyzeHcal(), lumiContext::fill, edm::EventSetup::get(), edm::Event::getByToken(), ignoreL1_, cuy::ii, edm::HandleBase::isValid(), LogDebug, convertSQLiteXML::ok, edm::Handle< T >::product(), edm::ESHandle< class >::product(), reco::return(), rnnum, edm::Event::run(), runNZS_, edm::SortedCollection< T, SORT >::size(), AlCaHLTBitMon_QueryRunRegistry::string, theRecalib_, tok_hbheNormal_, tok_hbherecoMB_, tok_hbherecoNoise_, tok_hfrecoMB_, tok_hfrecoNoise_, and tok_hltL1GtMap_.

                                                                                 {
    
  rnnum = (float)iEvent.run(); 
  const HcalRespCorrs* myRecalib=0;
  if (theRecalib_) {
    edm::ESHandle <HcalRespCorrs> recalibCorrs;
    iSetup.get<HcalRespCorrsRcd>().get("recalibrate",recalibCorrs);
    myRecalib = recalibCorrs.product();
  } // theRecalib
    
  edm::Handle<HBHERecHitCollection> hbheNormal;
  iEvent.getByToken(tok_hbheNormal_, hbheNormal);
  if (!hbheNormal.isValid()) {  
    edm::LogInfo("AnalyzerMB") << " hbheNormal failed";
  } else {
    edm::LogInfo("AnalyzerMB") << " The size of the normal collection "
                   << hbheNormal->size();
  }

  edm::Handle<HBHERecHitCollection> hbheNS;
  iEvent.getByToken(tok_hbherecoNoise_, hbheNS);
  if (!hbheNS.isValid()) {
    edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbheNoise product!";
    return ;
  }
  const HBHERecHitCollection HithbheNS = *(hbheNS.product());
  edm::LogInfo("AnalyzerMB") << "HBHE NS size of collection " << HithbheNS.size();
  if (runNZS_  && HithbheNS.size() != 5184) {
    edm::LogWarning("AnalyzerMB") << "HBHE NS problem " << rnnum << " size "
                  << HithbheNS.size();
    return;
  }
  
  edm::Handle<HBHERecHitCollection> hbheMB;
  iEvent.getByToken(tok_hbherecoMB_, hbheMB);
  if (!hbheMB.isValid()) {
    edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe product!";
    return ;
  }
  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
  edm::LogInfo("AnalyzerMB") << "HBHE MB size of collection " << HithbheMB.size();
  if (runNZS_  && HithbheMB.size() != 5184) {
    edm::LogWarning("AnalyzerMB") << "HBHE problem " << rnnum << " size "
                  << HithbheMB.size();
    return;
  }
    
  edm::Handle<HFRecHitCollection> hfNS;
  iEvent.getByToken(tok_hfrecoNoise_, hfNS);
  if (!hfNS.isValid()) {
    edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hfNoise product!";
    return ;
  }
  const HFRecHitCollection HithfNS = *(hfNS.product());
  edm::LogInfo("AnalyzerMB") << "HF NS size of collection "<< HithfNS.size();
  if (runNZS_  && HithfNS.size() != 1728) {
    edm::LogWarning("AnalyzerMB") << "HF NS problem " << rnnum << " size "
                  << HithfNS.size();
    return;
  }
  
  edm::Handle<HFRecHitCollection> hfMB;
  iEvent.getByToken(tok_hfrecoMB_, hfMB);
  if (!hfMB.isValid()) {
    edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf product!";
    return ;
  }
  const HFRecHitCollection HithfMB = *(hfMB.product());
  edm::LogInfo("AnalyzerMB") << "HF MB size of collection " << HithfMB.size();
  if(runNZS_  && HithfMB.size() != 1728) {
    edm::LogWarning("AnalyzerMB") << "HF problem " << rnnum << " size "
                  << HithfMB.size();
    return;
  }
    
  if (ignoreL1_) {
    analyzeHcal(myRecalib,HithbheNS,HithbheMB,HithfNS,HithfMB,1,true);
  } else {
    edm::Handle<L1GlobalTriggerObjectMapRecord> gtObjectMapRecord;
    iEvent.getByToken(tok_hltL1GtMap_, gtObjectMapRecord);
    if (gtObjectMapRecord.isValid()) {
      const std::vector<L1GlobalTriggerObjectMap>& objMapVec = gtObjectMapRecord->gtObjectMap();
      int  ii(0);
      bool ok(false), fill(true);
      for (std::vector<L1GlobalTriggerObjectMap>::const_iterator itMap = objMapVec.begin();
       itMap != objMapVec.end(); ++itMap, ++ii) {
    bool resultGt = (*itMap).algoGtlResult();
    if (resultGt == 1) {
      ok = true;
      int algoBit = (*itMap).algoBitNumber();
      analyzeHcal(myRecalib,HithbheNS,HithbheMB,HithfNS,HithfMB,algoBit,fill);
      fill = false;
      std::string algoNameStr = (*itMap).algoName();
      LogDebug("AnalyzerMB") << "Trigger[" << ii << "] " << algoNameStr
                 << " bit " << algoBit << " entered";
    }
      }
      if (!ok) edm::LogInfo("AnalyzerMB") << "No passed L1 Triggers";
    }
  }
}

void AnalyzerMinbias::analyzeHcal ( const HcalRespCorrs *  myRecalib, const HBHERecHitCollection &  HithbheNS, const HBHERecHitCollection &  HithbheMB, const HFRecHitCollection &  HithfNS, const HFRecHitCollection &  HithfMB, int  algoBit, bool  fill  )

private

Definition at line 352 of file AnalyzerMinbias.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), h_Noise, h_Signal, info(), myMap_, DetId::rawId(), rnnum, HcalDetId::subdet(), AnalyzerMinbias::myInfo::theNS0, and theRecalib_.

Referenced by analyze().

                                          {

  // Noise part for HB HE
  std::map<std::pair<int,HcalDetId>,myInfo> tmpMap;
  tmpMap.clear();
  
  for (HBHERecHitCollection::const_iterator hbheItr=HithbheNS.begin(); 
       hbheItr!=HithbheNS.end(); hbheItr++) {
    
    // Recalibration of energy
    float icalconst=1.;  
    DetId mydetid = hbheItr->id().rawId();
    if (theRecalib_) icalconst=myRecalib->getValues(mydetid)->getValue();
    
    HBHERecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
    double energyhit = aHit.energy();
    
    DetId id = (*hbheItr).detid(); 
    HcalDetId hid=HcalDetId(id);
    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    if (itr1 == myMap_.end()) {
      myInfo info;
      myMap_[std::pair<int,HcalDetId>(algoBit,hid)] = info;
      itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    }
    itr1->second.theNS0++;
    itr1->second.theNS1 += energyhit;
    itr1->second.theNS2 += (energyhit*energyhit);
    itr1->second.theNS3 += (energyhit*energyhit*energyhit);
    itr1->second.theNS4 += (energyhit*energyhit*energyhit*energyhit);
    itr1->second.runcheck = rnnum;
    if (fill) h_Noise[hid.subdet()-1]->Fill(energyhit);

    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr2 = tmpMap.find(std::pair<int,HcalDetId>(algoBit,hid));
    if (itr2 == tmpMap.end()) {
      myInfo info;
      tmpMap[std::pair<int,HcalDetId>(algoBit,hid)] = info;
      itr2 = tmpMap.find(std::pair<int,HcalDetId>(algoBit,hid));
    }
    itr2->second.theNS0++;
    itr2->second.theNS1 += energyhit;
    itr2->second.theNS2 += (energyhit*energyhit);
    itr2->second.theNS3 += (energyhit*energyhit*energyhit);
    itr2->second.theNS4 += (energyhit*energyhit*energyhit*energyhit);
    itr2->second.runcheck = rnnum;
      
  } // HBHE_NS
  
    // Signal part for HB HE
  
  for (HBHERecHitCollection::const_iterator hbheItr=HithbheMB.begin(); 
       hbheItr!=HithbheMB.end(); hbheItr++) {
    // Recalibration of energy
    float icalconst=1.;  
    DetId mydetid = hbheItr->id().rawId();
    if (theRecalib_) icalconst=myRecalib->getValues(mydetid)->getValue();
    
    HBHERecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
    double energyhit = aHit.energy();
    
    DetId id = (*hbheItr).detid(); 
    HcalDetId hid=HcalDetId(id);
          
    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr2 = tmpMap.find(std::pair<int,HcalDetId>(algoBit,hid));
          
    if (itr1 == myMap_.end()) {
      myInfo info;
      myMap_[std::pair<int,HcalDetId>(algoBit,hid)] = info;
      itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    } 
    itr1->second.theMB0++;
    itr1->second.theDif0 = 0;
    itr1->second.theMB1 += energyhit;
    itr1->second.theMB2 += (energyhit*energyhit);
    itr1->second.theMB3 += (energyhit*energyhit*energyhit);
    itr1->second.theMB4 += (energyhit*energyhit*energyhit*energyhit);
    itr1->second.runcheck = rnnum;
    float mydiff = 0.0;
    if (itr2 !=tmpMap.end()) {
      mydiff = energyhit - (itr2->second.theNS1);
      itr1->second.theDif0++;
      itr1->second.theDif1 += mydiff;
      itr1->second.theDif2 += (mydiff*mydiff);
      if (fill) h_Signal[hid.subdet()-1]->Fill(mydiff);
    }
  } // HBHE_MB
      
    // HF
  
  for (HFRecHitCollection::const_iterator hbheItr=HithfNS.begin(); 
       hbheItr!=HithfNS.end(); hbheItr++) {
    // Recalibration of energy
    float icalconst=1.;  
    DetId mydetid = hbheItr->id().rawId();
    if (theRecalib_) icalconst=myRecalib->getValues(mydetid)->getValue();
    
    HFRecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
    double energyhit = aHit.energy();
    // Remove PMT hits
    if(fabs(energyhit) > 40. ) continue;
    DetId id = (*hbheItr).detid(); 
    HcalDetId hid=HcalDetId(id);
    
    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    
    if (itr1 == myMap_.end()) {
      myInfo info;
      myMap_[std::pair<int,HcalDetId>(algoBit,hid)] = info;
      itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    } 
    itr1->second.theNS0++;
    itr1->second.theNS1 += energyhit;
    itr1->second.theNS2 += (energyhit*energyhit);
    itr1->second.theNS3 += (energyhit*energyhit*energyhit);
    itr1->second.theNS4 += (energyhit*energyhit*energyhit*energyhit);
    itr1->second.runcheck = rnnum;
    if (fill) h_Noise[hid.subdet()-1]->Fill(energyhit);
 
    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr2 = tmpMap.find(std::pair<int,HcalDetId>(algoBit,hid));
    if (itr2 == tmpMap.end()) {
      myInfo info;
      tmpMap[std::pair<int,HcalDetId>(algoBit,hid)] = info;
      itr2 = tmpMap.find(std::pair<int,HcalDetId>(algoBit,hid));
    }
    itr2->second.theNS0++;
    itr2->second.theNS1 += energyhit;
    itr2->second.theNS2 += (energyhit*energyhit);
    itr2->second.theNS3 += (energyhit*energyhit*energyhit);
    itr2->second.theNS4 += (energyhit*energyhit*energyhit*energyhit);
    itr2->second.runcheck = rnnum;
    
  } // HF_NS
      
      
    // Signal part for HF
      
  for (HFRecHitCollection::const_iterator hbheItr=HithfMB.begin(); 
       hbheItr!=HithfMB.end(); hbheItr++) {
    // Recalibration of energy
    float icalconst=1.;  
    DetId mydetid = hbheItr->id().rawId();
    if (theRecalib_) icalconst=myRecalib->getValues(mydetid)->getValue();
    HFRecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
    
    double energyhit = aHit.energy();
    // Remove PMT hits
    if(fabs(energyhit) > 40. ) continue;
          
    DetId id = (*hbheItr).detid(); 
    HcalDetId hid=HcalDetId(id);
          
    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    std::map<std::pair<int,HcalDetId>,myInfo>::iterator itr2 = tmpMap.find(std::pair<int,HcalDetId>(algoBit,hid));
    
    if (itr1 == myMap_.end()) {
      myInfo info;
      myMap_[std::pair<int,HcalDetId>(algoBit,hid)] = info;
      itr1 = myMap_.find(std::pair<int,HcalDetId>(algoBit,hid));
    }
    itr1->second.theMB0++;
    itr1->second.theDif0 = 0;
    itr1->second.theMB1 += energyhit;
    itr1->second.theMB2 += (energyhit*energyhit);
    itr1->second.theMB3 += (energyhit*energyhit*energyhit);
    itr1->second.theMB4 += (energyhit*energyhit*energyhit*energyhit);
    itr1->second.runcheck = rnnum;
    float mydiff = 0.0;
    if (itr2 !=tmpMap.end()) {
      mydiff = energyhit - (itr2->second.theNS1);
      itr1->second.theDif0++;
      itr1->second.theDif1 += mydiff;
      itr1->second.theDif2 += (mydiff*mydiff);
      if (fill) h_Signal[hid.subdet()-1]->Fill(mydiff);
    }
  }
}

void AnalyzerMinbias::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 145 of file AnalyzerMinbias.cc.

References cells, depth, eta, fOutputFileName, h_Noise, h_Signal, hOutputFile, ieta, iphi, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom3_MB, mom4_MB, mom4_Noise, mydet, myMap_, mysubd, myTree, mergeVDriftHistosByStation::name, occup, phi, rnnumber, AlCaHLTBitMon_QueryRunRegistry::string, indexGen::title, and trigbit.

                               {

  std::string det[4] = {"HB", "HE", "HO", "HF"};
  char        name[80], title[80];
  for (int subd=0; subd<4; ++subd) {
    sprintf(name,"Noise_%s",det[subd].c_str());
    sprintf(title,"Energy Distribution for Noise in %s",det[subd].c_str());
    h_Noise[subd] = new TH1D(name,title,100,-10.,10.);
    sprintf(name,"Signal_%s",det[subd].c_str());
    sprintf(title,"Energy Distribution for Signal in %s",det[subd].c_str());
    h_Signal[subd] = new TH1D(name,title,100,-10.,10.);
  }

  hOutputFile   = new TFile(fOutputFileName.c_str(), "RECREATE") ;
  myTree = new TTree("RecJet","RecJet Tree");
  myTree->Branch("mydet",       &mydet,      "mydet/I");
  myTree->Branch("mysubd",      &mysubd,     "mysubd/I");
  myTree->Branch("cells",       &cells,      "cells");
  myTree->Branch("depth",       &depth,      "depth/I");
  myTree->Branch("ieta",        &ieta,       "ieta/I");
  myTree->Branch("iphi",        &iphi,       "iphi/I");
  myTree->Branch("eta",         &eta,        "eta/F");
  myTree->Branch("phi",         &phi,        "phi/F");
  myTree->Branch("mom0_MB",     &mom0_MB,    "mom0_MB/F");
  myTree->Branch("mom1_MB",     &mom1_MB,    "mom1_MB/F");
  myTree->Branch("mom2_MB",     &mom2_MB,    "mom2_MB/F");
  myTree->Branch("mom3_MB",     &mom3_MB,    "mom3_MB/F");
  myTree->Branch("mom4_MB",     &mom4_MB,    "mom4_MB/F");
  myTree->Branch("mom0_Noise",  &mom0_Noise, "mom0_Noise/F");
  myTree->Branch("mom1_Noise",  &mom1_Noise, "mom1_Noise/F");
  myTree->Branch("mom2_Noise",  &mom2_Noise, "mom2_Noise/F");
  myTree->Branch("mom3_Noise",  &mom2_Noise, "mom3_Noise/F");
  myTree->Branch("mom4_Noise",  &mom4_Noise, "mom4_Noise/F");
  myTree->Branch("mom0_Diff",   &mom0_Diff,  "mom0_Diff/F");
  myTree->Branch("mom1_Diff",   &mom1_Diff,  "mom1_Diff/F");
  myTree->Branch("mom2_Diff",   &mom2_Diff,  "mom2_Diff/F");
  myTree->Branch("occup",       &occup,      "occup/F");
  myTree->Branch("trigbit",     &trigbit,    "trigbit/I");
  myTree->Branch("rnnumber",    &rnnumber,   "rnnumber/D");

  myMap_.clear();
}

void AnalyzerMinbias::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 190 of file AnalyzerMinbias.cc.

References cells, depth, h_Noise, h_Signal, hOutputFile, i, ieta, cuy::ii, info(), iphi, LogDebug, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom3_MB, mom3_Noise, mom4_MB, mom4_Noise, myMap_, mysubd, myTree, rnnumber, AnalyzerMinbias::myInfo::runcheck, AnalyzerMinbias::myInfo::theDif0, AnalyzerMinbias::myInfo::theDif1, AnalyzerMinbias::myInfo::theDif2, AnalyzerMinbias::myInfo::theMB0, AnalyzerMinbias::myInfo::theMB1, AnalyzerMinbias::myInfo::theMB2, AnalyzerMinbias::myInfo::theMB3, AnalyzerMinbias::myInfo::theMB4, AnalyzerMinbias::myInfo::theNS0, AnalyzerMinbias::myInfo::theNS1, AnalyzerMinbias::myInfo::theNS2, AnalyzerMinbias::myInfo::theNS3, AnalyzerMinbias::myInfo::theNS4, and trigbit.

                             {
   
  int ii=0;
  for (std::map<std::pair<int,HcalDetId>,myInfo>::const_iterator itr=myMap_.begin(); itr != myMap_.end(); ++itr) {
    LogDebug("AnalyzerMB") << "Fired trigger bit number " << itr->first.first; 
    myInfo info = itr->second;
    if (info.theMB0 > 0) { 
      mom0_MB    = info.theMB0;
      mom1_MB    = info.theMB1;
      mom2_MB    = info.theMB2;
      mom3_MB    = info.theMB3;
      mom4_MB    = info.theMB4;
      mom0_Noise = info.theNS0;
      mom1_Noise = info.theNS1;
      mom2_Noise = info.theNS2;
      mom3_Noise = info.theNS3;
      mom4_Noise = info.theNS4;
      mom0_Diff  = info.theDif0;
      mom1_Diff  = info.theDif1;
      mom2_Diff  = info.theDif2;
      rnnumber   = info.runcheck;
      trigbit    = itr->first.first; 
      mysubd     = itr->first.second.subdet();
      depth      = itr->first.second.depth();
      ieta       = itr->first.second.ieta();
      iphi       = itr->first.second.iphi();
      
      LogDebug("AnalyzerMB") << " Result=  " << trigbit << " " << mysubd
                             << " " << ieta << " " << iphi << " mom0  "
                             << mom0_MB << " mom1 " << mom1_MB << " mom2 "
                             << mom2_MB << " mom3 " << mom3_MB << " mom4 " 
                 << mom4_MB << " mom0_Noise " << mom0_Noise 
                             << " mom1_Noise " << mom1_Noise << " mom2_Noise "
                             << mom2_Noise << " mom3_Noise " << mom3_Noise 
                             << " mom4_Noise " << mom4_Noise << " mom0_Diff "
                             << mom0_Diff << " mom1_Diff " << mom1_Diff
                             << " mom2_Diff " << mom2_Diff;
      myTree->Fill();
      ii++;
    }
  }
  cells = ii; 
  LogDebug("AnalyzerMB") << "cells" << " " << cells;    
  hOutputFile->Write();   
  hOutputFile->cd();
  myTree->Write();
  for(int i=0; i<4; i++){
    h_Noise[i]->Write();
    h_Signal[i]->Write();
  }
  hOutputFile->Close() ;
}

Member Data Documentation

int AnalyzerMinbias::cells

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int AnalyzerMinbias::depth

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by dqm_interfaces.DirID::__eq__(), dqm_interfaces.DirID::__repr__(), beginJob(), endJob(), and dqm_interfaces.DirWalkerDB::run().

float AnalyzerMinbias::eta

private

Definition at line 107 of file AnalyzerMinbias.cc.

Referenced by Particle.Particle::__str__(), beginJob(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

std::string AnalyzerMinbias::fOutputFileName

private

Definition at line 96 of file AnalyzerMinbias.cc.

Referenced by AnalyzerMinbias(), and beginJob().

TH1D* AnalyzerMinbias::h_Noise[4]

private

Definition at line 100 of file AnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and endJob().

TH1D * AnalyzerMinbias::h_Signal[4]

private

Definition at line 100 of file AnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and endJob().

std::string AnalyzerMinbias::hcalfile_

private

Definition at line 96 of file AnalyzerMinbias.cc.

TFile* AnalyzerMinbias::hOutputFile

private

Definition at line 98 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int AnalyzerMinbias::ieta

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

bool AnalyzerMinbias::ignoreL1_

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

int AnalyzerMinbias::iphi

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom0_Diff

private

Definition at line 110 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom0_MB

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom0_Noise

private

Definition at line 109 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom1_Diff

private

Definition at line 110 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom1_MB

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom1_Noise

private

Definition at line 109 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom2_Diff

private

Definition at line 110 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom2_MB

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom2_Noise

private

Definition at line 109 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom3_Diff

private

Definition at line 110 of file AnalyzerMinbias.cc.

float AnalyzerMinbias::mom3_MB

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom3_Noise

private

Definition at line 109 of file AnalyzerMinbias.cc.

Referenced by endJob().

float AnalyzerMinbias::mom4_Diff

private

Definition at line 110 of file AnalyzerMinbias.cc.

float AnalyzerMinbias::mom4_MB

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::mom4_Noise

private

Definition at line 109 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int AnalyzerMinbias::mydet

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob().

std::map<std::pair<int,HcalDetId>,myInfo> AnalyzerMinbias::myMap_

private

Definition at line 112 of file AnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and endJob().

std::ofstream* AnalyzerMinbias::myout_hcal

private

Definition at line 97 of file AnalyzerMinbias.cc.

int AnalyzerMinbias::mysubd

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

TTree* AnalyzerMinbias::myTree

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float AnalyzerMinbias::occup

private

Definition at line 108 of file AnalyzerMinbias.cc.

Referenced by beginJob().

float AnalyzerMinbias::phi

private

Definition at line 107 of file AnalyzerMinbias.cc.

Referenced by Particle.Particle::__str__(), and beginJob().

double AnalyzerMinbias::rnnum

private

Definition at line 102 of file AnalyzerMinbias.cc.

Referenced by analyze(), and analyzeHcal().

double AnalyzerMinbias::rnnumber

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

bool AnalyzerMinbias::runNZS_

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

bool AnalyzerMinbias::theRecalib_

private

Definition at line 101 of file AnalyzerMinbias.cc.

Referenced by analyze(), analyzeHcal(), and AnalyzerMinbias().

edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbheNormal_

private

Definition at line 116 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbherecoMB_

private

Definition at line 113 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbherecoNoise_

private

Definition at line 113 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

edm::EDGetTokenT<HFRecHitCollection> AnalyzerMinbias::tok_hfrecoMB_

private

Definition at line 114 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

edm::EDGetTokenT<HFRecHitCollection> AnalyzerMinbias::tok_hfrecoNoise_

private

Definition at line 114 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> AnalyzerMinbias::tok_hltL1GtMap_

private

Definition at line 117 of file AnalyzerMinbias.cc.

Referenced by analyze(), and AnalyzerMinbias().

edm::EDGetTokenT<HORecHitCollection> AnalyzerMinbias::tok_horecoMB_

private

Definition at line 115 of file AnalyzerMinbias.cc.

Referenced by AnalyzerMinbias().

edm::EDGetTokenT<HORecHitCollection> AnalyzerMinbias::tok_horecoNoise_

private

Definition at line 115 of file AnalyzerMinbias.cc.

Referenced by AnalyzerMinbias().

int AnalyzerMinbias::trigbit

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().