AnalyzerMinbias Class Reference

Inheritance diagram for AnalyzerMinbias:

Classes
struct	myInfo

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	AnalyzerMinbias (const edm::ParameterSet &)
void	beginJob () override
void	endJob () override
	~AnalyzerMinbias () override=default
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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
const std::string	fOutputFileName
TH1D *	h_Noise [4]
TH1D *	h_Signal [4]
std::string	hcalfile_
TFile *	hOutputFile_
int	ieta
const 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
const bool	runNZS_
const bool	theRecalib_
const edm::EDGetTokenT< HBHERecHitCollection >	tok_hbheNormal_
const edm::EDGetTokenT< HBHERecHitCollection >	tok_hbherecoMB_
const edm::EDGetTokenT< HBHERecHitCollection >	tok_hbherecoNoise_
const edm::EDGetTokenT< HFRecHitCollection >	tok_hfrecoMB_
const edm::EDGetTokenT< HFRecHitCollection >	tok_hfrecoNoise_
const edm::EDGetTokenT< L1GlobalTriggerObjectMapRecord >	tok_hltL1GtMap_
const edm::EDGetTokenT< HORecHitCollection >	tok_horecoMB_
const edm::EDGetTokenT< HORecHitCollection >	tok_horecoNoise_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_respCorr_
int	trigbit

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 62 of file AnalyzerMinbias.cc.

Constructor & Destructor Documentation

AnalyzerMinbias()

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

explicit

Definition at line 118 of file AnalyzerMinbias.cc.

References TFileService::kSharedResource.

    : fOutputFileName(iConfig.getUntrackedParameter<std::string>("HistOutFile")),
      runNZS_(iConfig.getUntrackedParameter<bool>("RunNZS", true)),
      theRecalib_(iConfig.getParameter<bool>("Recalib")),
      ignoreL1_(iConfig.getUntrackedParameter<bool>("IgnoreL1", true)),
      tok_hbherecoMB_(consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"))),
      tok_hbherecoNoise_(consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputNoise"))),
      tok_horecoMB_(consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputMB"))),
      tok_horecoNoise_(consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputNoise"))),
      tok_hfrecoMB_(consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"))),
      tok_hfrecoNoise_(consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputNoise"))),
      tok_hbheNormal_(consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"))),
      tok_hltL1GtMap_(consumes<L1GlobalTriggerObjectMapRecord>(edm::InputTag("hltL1GtObjectMap"))),
      tok_respCorr_(esConsumes<HcalRespCorrs, HcalRespCorrsRcd>()) {
  usesResource(TFileService::kSharedResource);
  // get name of output file with histogramms
  // get token names of modules, producing object collections
}

~AnalyzerMinbias()

AnalyzerMinbias::~AnalyzerMinbias ( )

overridedefault

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 241 of file AnalyzerMinbias.cc.

References analyzeHcal(), ntuplemaker::fill, nano_mu_digi_cff::float, edm::EventSetup::getData(), L1GlobalTriggerObjectMapRecord::gtObjectMap(), iEvent, ignoreL1_, cuy::ii, edm::HandleBase::isValid(), convertSQLiteXML::ok, edm::Handle< T >::product(), rnnum_, runNZS_, edm::SortedCollection< T, SORT >::size(), AlCaHLTBitMon_QueryRunRegistry::string, theRecalib_, tok_hbheNormal_, tok_hbherecoMB_, tok_hbherecoNoise_, tok_hfrecoMB_, tok_hfrecoNoise_, tok_hltL1GtMap_, and tok_respCorr_.

                                                                                 {
  rnnum_ = (float)iEvent.run();
  const HcalRespCorrs* myRecalib = nullptr;
  if (theRecalib_) {
    myRecalib = &iSetup.getData(tok_respCorr_);
  }  // theRecalib

  const edm::Handle<HBHERecHitCollection> hbheNormal = iEvent.getHandle(tok_hbheNormal_);
  if (!hbheNormal.isValid()) {
    edm::LogVerbatim("AnalyzerMinimumBias") << " hbheNormal failed";
  } else {
    edm::LogVerbatim("AnalyzerMinimumBias") << " The size of the normal collection " << hbheNormal->size();
  }

  const edm::Handle<HBHERecHitCollection> hbheNS = iEvent.getHandle(tok_hbherecoNoise_);
  if (!hbheNS.isValid()) {
    edm::LogWarning("AnalyzerMinimumBias") << "HcalCalibAlgos: Error! can't get hbheNoise product!";
    return;
  }
  const HBHERecHitCollection HithbheNS = *(hbheNS.product());
  edm::LogVerbatim("AnalyzerMinimumBias") << "HBHE NS size of collection " << HithbheNS.size();
  if (runNZS_ && HithbheNS.size() != 5184) {
    edm::LogWarning("AnalyzerMinimumBias") << "HBHE NS problem " << rnnum_ << " size " << HithbheNS.size();
    return;
  }

  const edm::Handle<HBHERecHitCollection> hbheMB = iEvent.getHandle(tok_hbherecoMB_);
  if (!hbheMB.isValid()) {
    edm::LogWarning("AnalyzerMinimumBias") << "HcalCalibAlgos: Error! can't get hbhe product!";
    return;
  }
  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
  edm::LogVerbatim("AnalyzerMinimumBias") << "HBHE MB size of collection " << HithbheMB.size();
  if (runNZS_ && HithbheMB.size() != 5184) {
    edm::LogWarning("AnalyzerMinimumBias") << "HBHE problem " << rnnum_ << " size " << HithbheMB.size();
    return;
  }

  const edm::Handle<HFRecHitCollection> hfNS = iEvent.getHandle(tok_hfrecoNoise_);
  if (!hfNS.isValid()) {
    edm::LogWarning("AnalyzerMinimumBias") << "HcalCalibAlgos: Error! can't get hfNoise product!";
    return;
  }
  const HFRecHitCollection HithfNS = *(hfNS.product());
  edm::LogVerbatim("AnalyzerMinimumBias") << "HF NS size of collection " << HithfNS.size();
  if (runNZS_ && HithfNS.size() != 1728) {
    edm::LogWarning("AnalyzerMinimumBias") << "HF NS problem " << rnnum_ << " size " << HithfNS.size();
    return;
  }

  const edm::Handle<HFRecHitCollection> hfMB = iEvent.getHandle(tok_hfrecoMB_);
  if (!hfMB.isValid()) {
    edm::LogWarning("AnalyzerMinimumBias") << "HcalCalibAlgos: Error! can't get hf product!";
    return;
  }
  const HFRecHitCollection HithfMB = *(hfMB.product());
  edm::LogVerbatim("AnalyzerMinimumBias") << "HF MB size of collection " << HithfMB.size();
  if (runNZS_ && HithfMB.size() != 1728) {
    edm::LogWarning("AnalyzerMinimumBias") << "HF problem " << rnnum_ << " size " << HithfMB.size();
    return;
  }

  if (ignoreL1_) {
    analyzeHcal(myRecalib, HithbheNS, HithbheMB, HithfNS, HithfMB, 1, true);
  } else {
    const edm::Handle<L1GlobalTriggerObjectMapRecord> gtObjectMapRecord = iEvent.getHandle(tok_hltL1GtMap_);
    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();
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("AnalyzerMinimumBias")
              << "Trigger[" << ii << "] " << algoNameStr << " bit " << algoBit << " entered";
#endif
        }
      }
      if (!ok)
        edm::LogVerbatim("AnalyzerMinimumBias") << "No passed L1 Triggers";
    }
  }
}

analyzeHcal()

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 332 of file AnalyzerMinbias.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), ntuplemaker::fill, HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), h_Noise, h_Signal, info(), myMap_, DetId::rawId(), rnnum_, HcalDetId::subdet(), 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);
    }
  }
}

beginJob()

void AnalyzerMinbias::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 137 of file AnalyzerMinbias.cc.

References cells, depth, eta, fOutputFileName, compareTotals::fs, 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, mom3_Noise, mom4_MB, mom4_Noise, mydet, myMap_, mysubd, myTree_, Skims_PA_cff::name, occup, phi, rnnumber, AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::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.);
  }

  edm::Service<TFileService> fs;
  hOutputFile_ = new TFile(fOutputFileName.c_str(), "RECREATE");
  myTree_ = fs->make<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", &mom3_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();
}

endJob()

void AnalyzerMinbias::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 182 of file AnalyzerMinbias.cc.

References cells, depth, h_Noise, h_Signal, hOutputFile_, mps_fire::i, ieta, cuy::ii, info(), iphi, genParticles_cff::map, 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, and trigbit.

Referenced by o2olib.O2ORunMgr::executeJob().

                             {
  int ii = 0;
  for (std::map<std::pair<int, HcalDetId>, myInfo>::const_iterator itr = myMap_.begin(); itr != myMap_.end(); ++itr) {
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("AnalyzerMinimumBias") << "Fired trigger bit number " << itr->first.first;
#endif
    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();
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("AnalyzerMinimumBias")
          << " 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;
#endif
      myTree_->Fill();
      ii++;
    }
  }
  cells = ii;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("AnalyzerMinimumBias") << "cells " << cells;
#endif
  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

cells

int AnalyzerMinbias::cells

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

depth

int AnalyzerMinbias::depth

private

Definition at line 103 of file AnalyzerMinbias.cc.

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

eta

float AnalyzerMinbias::eta

private

Definition at line 104 of file AnalyzerMinbias.cc.

Referenced by Particle.Particle::__str__(), beginJob(), datamodel.Object::DeltaR(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

fOutputFileName

const std::string AnalyzerMinbias::fOutputFileName

private

Definition at line 92 of file AnalyzerMinbias.cc.

Referenced by beginJob().

h_Noise

TH1D* AnalyzerMinbias::h_Noise[4]

private

Definition at line 98 of file AnalyzerMinbias.cc.

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

h_Signal

TH1D * AnalyzerMinbias::h_Signal[4]

private

Definition at line 98 of file AnalyzerMinbias.cc.

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

hcalfile_

std::string AnalyzerMinbias::hcalfile_

private

Definition at line 94 of file AnalyzerMinbias.cc.

hOutputFile_

TFile* AnalyzerMinbias::hOutputFile_

private

Definition at line 96 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

ieta

int AnalyzerMinbias::ieta

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

ignoreL1_

const bool AnalyzerMinbias::ignoreL1_

private

Definition at line 93 of file AnalyzerMinbias.cc.

Referenced by analyze().

iphi

int AnalyzerMinbias::iphi

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom0_Diff

float AnalyzerMinbias::mom0_Diff

private

Definition at line 107 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom0_MB

float AnalyzerMinbias::mom0_MB

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom0_Noise

float AnalyzerMinbias::mom0_Noise

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom1_Diff

float AnalyzerMinbias::mom1_Diff

private

Definition at line 107 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom1_MB

float AnalyzerMinbias::mom1_MB

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom1_Noise

float AnalyzerMinbias::mom1_Noise

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom2_Diff

float AnalyzerMinbias::mom2_Diff

private

Definition at line 107 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom2_MB

float AnalyzerMinbias::mom2_MB

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom2_Noise

float AnalyzerMinbias::mom2_Noise

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom3_Diff

float AnalyzerMinbias::mom3_Diff

private

Definition at line 107 of file AnalyzerMinbias.cc.

mom3_MB

float AnalyzerMinbias::mom3_MB

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom3_Noise

float AnalyzerMinbias::mom3_Noise

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom4_Diff

float AnalyzerMinbias::mom4_Diff

private

Definition at line 107 of file AnalyzerMinbias.cc.

mom4_MB

float AnalyzerMinbias::mom4_MB

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mom4_Noise

float AnalyzerMinbias::mom4_Noise

private

Definition at line 106 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

mydet

int AnalyzerMinbias::mydet

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob().

myMap_

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

private

Definition at line 109 of file AnalyzerMinbias.cc.

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

myout_hcal

std::ofstream* AnalyzerMinbias::myout_hcal

private

Definition at line 95 of file AnalyzerMinbias.cc.

mysubd

int AnalyzerMinbias::mysubd

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

myTree_

TTree* AnalyzerMinbias::myTree_

private

Definition at line 97 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

occup

float AnalyzerMinbias::occup

private

Definition at line 105 of file AnalyzerMinbias.cc.

Referenced by beginJob().

phi

float AnalyzerMinbias::phi

private

Definition at line 104 of file AnalyzerMinbias.cc.

Referenced by Particle.Particle::__str__(), beginJob(), datamodel.Object::DeltaR(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

rnnum_

double AnalyzerMinbias::rnnum_

private

Definition at line 99 of file AnalyzerMinbias.cc.

Referenced by analyze(), and analyzeHcal().

rnnumber

double AnalyzerMinbias::rnnumber

private

Definition at line 102 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

runNZS_

const bool AnalyzerMinbias::runNZS_

private

Definition at line 93 of file AnalyzerMinbias.cc.

Referenced by analyze().

theRecalib_

const bool AnalyzerMinbias::theRecalib_

private

Definition at line 93 of file AnalyzerMinbias.cc.

Referenced by analyze(), and analyzeHcal().

tok_hbheNormal_

const edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbheNormal_

private

Definition at line 113 of file AnalyzerMinbias.cc.

Referenced by analyze().

tok_hbherecoMB_

const edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbherecoMB_

private

Definition at line 110 of file AnalyzerMinbias.cc.

Referenced by analyze().

tok_hbherecoNoise_

const edm::EDGetTokenT<HBHERecHitCollection> AnalyzerMinbias::tok_hbherecoNoise_

private

Definition at line 110 of file AnalyzerMinbias.cc.

Referenced by analyze().

tok_hfrecoMB_

const edm::EDGetTokenT<HFRecHitCollection> AnalyzerMinbias::tok_hfrecoMB_

private

Definition at line 112 of file AnalyzerMinbias.cc.

Referenced by analyze().

tok_hfrecoNoise_

const edm::EDGetTokenT<HFRecHitCollection> AnalyzerMinbias::tok_hfrecoNoise_

private

Definition at line 112 of file AnalyzerMinbias.cc.

Referenced by analyze().

tok_hltL1GtMap_

const edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> AnalyzerMinbias::tok_hltL1GtMap_

private

Definition at line 114 of file AnalyzerMinbias.cc.

Referenced by analyze().

tok_horecoMB_

const edm::EDGetTokenT<HORecHitCollection> AnalyzerMinbias::tok_horecoMB_

private

Definition at line 111 of file AnalyzerMinbias.cc.

tok_horecoNoise_

const edm::EDGetTokenT<HORecHitCollection> AnalyzerMinbias::tok_horecoNoise_

private

Definition at line 111 of file AnalyzerMinbias.cc.

tok_respCorr_

const edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> AnalyzerMinbias::tok_respCorr_

private

Definition at line 115 of file AnalyzerMinbias.cc.

Referenced by analyze().

trigbit

int AnalyzerMinbias::trigbit

private

Definition at line 103 of file AnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().