cms::Analyzer_minbias Class Reference

#include <Analyzer_minbias.h>

Inheritance diagram for cms::Analyzer_minbias:

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
	Analyzer_minbias (const edm::ParameterSet &)
void	beginJob () override
void	beginRun (edm::Run const &, edm::EventSetup const &) override
void	endJob () override
void	endRun (edm::Run const &, edm::EventSetup const &) override
	~Analyzer_minbias () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const 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
bool	wantsStreamRuns () const
	~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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
int	depth
float	eta
std::string	fOutputFileName
edm::Service< TFileService >	fs
TH1F *	hbheNoiseE
TH1F *	hbheSignalE
std::string	hcalfile_
TH1F *	hCalo1 [73][43]
TH1F *	hCalo1mom2 [73][43]
TH1F *	hCalo2 [73][43]
TH1F *	hCalo2mom2 [73][43]
TH1F *	hfNoiseE
TH1F *	hfSignalE
TH2F *	hHBHEsize_vs_run
TH2F *	hHFsize_vs_run
TFile *	hOutputFile
int	ieta
int	iphi
double	meannoise_min [73][43]
double	meannoise_pl [73][43]
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::ofstream *	myout_hcal
int	mysubd
TTree *	myTree
double	nevent
int	nevent_run
double	noise_min [73][43]
double	noise_pl [73][43]
float	occup
float	phi
double	theDFFillDetMapMin0 [5][5][73][43]
double	theDFFillDetMapMin1 [5][5][73][43]
double	theDFFillDetMapMin2 [5][5][73][43]
double	theDFFillDetMapPl0 [5][5][73][43]
double	theDFFillDetMapPl1 [5][5][73][43]
double	theDFFillDetMapPl2 [5][5][73][43]
double	theMBFillDetMapMin0 [5][5][73][43]
double	theMBFillDetMapMin1 [5][5][73][43]
double	theMBFillDetMapMin2 [5][5][73][43]
double	theMBFillDetMapMin4 [5][5][73][43]
double	theMBFillDetMapPl0 [5][5][73][43]
double	theMBFillDetMapPl1 [5][5][73][43]
double	theMBFillDetMapPl2 [5][5][73][43]
double	theMBFillDetMapPl4 [5][5][73][43]
double	theNSFillDetMapMin0 [5][5][73][43]
double	theNSFillDetMapMin1 [5][5][73][43]
double	theNSFillDetMapMin2 [5][5][73][43]
double	theNSFillDetMapMin4 [5][5][73][43]
double	theNSFillDetMapPl0 [5][5][73][43]
double	theNSFillDetMapPl1 [5][5][73][43]
double	theNSFillDetMapPl2 [5][5][73][43]
double	theNSFillDetMapPl4 [5][5][73][43]
bool	theRecalib
edm::EDGetTokenT< FEDRawDataCollection >	tok_data_
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	tok_gtRec_
edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection >	tok_hbheNoise_
edm::EDGetTokenT< HBHERecHitCollection >	tok_hbheNorm_
edm::EDGetTokenT< HFRecHitCollection >	tok_hf_
edm::EDGetTokenT< HFRecHitCollection >	tok_hfNoise_
edm::EDGetTokenT< HORecHitCollection >	tok_ho_
edm::EDGetTokenT< HORecHitCollection >	tok_hoNoise_
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_respCorr_

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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 61 of file Analyzer_minbias.h.

Constructor & Destructor Documentation

Analyzer_minbias()

cms::Analyzer_minbias::Analyzer_minbias ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 43 of file Analyzer_minbias.cc.

                                                                   {
    usesResource(TFileService::kSharedResource);
    // get name of output file with histogramms
    fOutputFileName = iConfig.getUntrackedParameter<std::string>("HistOutFile");
    // get names of modules, producing object collections
 
    tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"));
    tok_ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputMB"));
    tok_hf_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"));
    tok_data_ = consumes<FEDRawDataCollection>(edm::InputTag(iConfig.getParameter<std::string>("InputLabel")));
 
    tok_hbheNoise_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputNoise"));
    tok_hoNoise_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputNoise"));
    tok_hfNoise_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputNoise"));
 
    // this was hardcodded..
    tok_gtRec_ = consumes<L1GlobalTriggerReadoutRecord>(edm::InputTag("gtDigisAlCaMB"));
    tok_hbheNorm_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
 
    tok_respCorr_ = esConsumes<HcalRespCorrs, HcalRespCorrsRcd>();
 
    theRecalib = iConfig.getParameter<bool>("Recalib");
 
    //
    //
    for (int i = 0; i < 73; i++) {
      for (int j = 0; j < 43; j++) {
        noise_min[i][j] = 0.;
        noise_pl[i][j] = 0.;
      }
    }
    //
    //
  }

References fOutputFileName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, dqmiolumiharvest::j, TFileService::kSharedResource, noise_min, noise_pl, AlCaHLTBitMon_QueryRunRegistry::string, theRecalib, tok_data_, tok_gtRec_, tok_hbhe_, tok_hbheNoise_, tok_hbheNorm_, tok_hf_, tok_hfNoise_, tok_ho_, tok_hoNoise_, and tok_respCorr_.

~Analyzer_minbias()

cms::Analyzer_minbias::~Analyzer_minbias ( )

override

Definition at line 78 of file Analyzer_minbias.cc.

                                      {
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
  }

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 321 of file Analyzer_minbias.cc.

                                                                                    {
    edm::LogInfo("AnalyzerMB") << " Start Analyzer_minbias::analyze " << nevent;
    nevent++;
    nevent_run++;
    using namespace edm;
 
    float rnnum = (float)iEvent.run();
 
    std::vector<StableProvenance const*> theProvenance;
    iEvent.getAllStableProvenance(theProvenance);
 
    for (auto const& provenance : theProvenance) {
      edm::LogInfo("AnalyzerMB") << " Print all process/modulelabel/product names " << provenance->processName()
                                 << " , " << provenance->moduleLabel() << " , " << provenance->productInstanceName();
    }
    /*
      edm::Handle<FEDRawDataCollection> rawdata;  
      iEvent.getByToken(tok_data_,rawdata);
 
      if (!rawdata.isValid()) {
      edm::LogInfo("AnalyzerMB")<<" No valid collection ";
      } else {
      edm::LogInfo("AnalyzerMB")<<" Valid collection ";
      int calibType = -1 ; int numEmptyFEDs = 0 ; 
      std::vector<int> calibTypeCounter(8,0) ; 
      for (int i=FEDNumbering::MINHCALFEDID;
        i<=FEDNumbering::MAXHCALFEDID; i++) {
    const FEDRawData& fedData = rawdata->FEDData(i) ; 
    if ( fedData.size() < 24 ) numEmptyFEDs++ ; 
    if ( fedData.size() < 24 ) continue ; 
    //      int value = ((const HcalDCCHeader*)(fedData.data()))->getCalibType() ; 
    //      calibTypeCounter.at(value)++ ; // increment the counter for this calib type
    }
    edm::LogInfo("AnalyzerMB")<<" NumFed "<<numEmptyFEDs<<" "<<calibType;
    }
    */
    /*
      std::vector<edm::Handle<FEDRawDataCollection> > rawdata1;
      iEvent.getManyByType(rawdata1);
      
      for(std::vector<edm::Handle<FEDRawDataCollection> >::const_iterator it = rawdata1.begin();it != rawdata1.end(); it++) {
 
      edm::LogInfo("AnalyzerMB")<<" Many by Type product name "<< (*it).provenance()->processName()<<
      " "<<(*it).provenance()->moduleLabel();
 
      if((*it).provenance()->moduleLabel() == "hltHcalCalibrationRaw") {
      int calibType = -1 ; int numEmptyFEDs = 0 ; 
  
      for (int i=FEDNumbering::MINHCALFEDID;
        i<=FEDNumbering::MAXHCALFEDID; i++) {
    const FEDRawData& fedData = (*it)->FEDData(i) ; 
    edm::LogInfo("AnalyzerMB")<<" FED size "<<fedData.size();
    if ( fedData.size() < 24 ) numEmptyFEDs++ ; 
    if ( fedData.size() < 24 ) continue ; 
    int value = ((const HcalDCCHeader*)(fedData.data()))->getCalibType() ; 
    edm::LogInfo("AnalyzerMB")<<" Value "<<value;
      }
      edm::LogInfo("AnalyzerMB")<<" Many by Type NumFed "<<numEmptyFEDs<<" "<<calibType;
     }
   }
 
   
    */
 
    // Geometry
    //    edm::ESHandle<CaloGeometry> pG;
    //    iSetup.get<CaloGeometryRecord>().get(pG);
    // ======
 
    /*
      edm::ESHandle<L1GtTriggerMenu> menuRcd;
      iSetup.get<L1GtTriggerMenuRcd>().get(menuRcd) ;
      const L1GtTriggerMenu* menu = menuRcd.product();
      const AlgorithmMap& bitMap = menu->gtAlgorithmMap();
 
      edm::Handle<L1GlobalTriggerReadoutRecord> gtRecord;
      iEvent.getByToken(tok_gtRec_, gtRecord);
 
      if (!gtRecord.isValid()) {
 
      //     LogDebug("L1GlobalTriggerRecordProducer")
      //       << "\n\n Error: no L1GlobalTriggerReadoutRecord found with input tag "
      //       << m_l1GtReadoutRecord
      //       << "\n Returning empty L1GlobalTriggerRecord.\n\n";
      edm::LogInfo("AnalyzerMB")<<" No L1 trigger record ";
      } else {
 
      const DecisionWord dWord = gtRecord->decisionWord();
 
      for (CItAlgo itAlgo = bitMap.begin(); itAlgo != bitMap.end(); itAlgo++) {
      bool decision=menu->gtAlgorithmResult(itAlgo->first,dWord);
      if(decision == 1) edm::LogInfo("AnalyzerMB")<<" Trigger "<<itAlgo->first<<" "<<decision;
      }
      
    }
    */
 
    const HcalRespCorrs* myRecalib = nullptr;
    if (theRecalib) {
      myRecalib = &iSetup.getData(tok_respCorr_);
    }  // theRecalib
 
    // Noise part for HB HE
 
    double tmpNSFillDetMapPl1[5][5][73][43];
    double tmpNSFillDetMapMin1[5][5][73][43];
 
    for (int i = 0; i < 5; i++) {
      for (int j = 0; j < 5; j++) {
        for (int k = 0; k < 73; k++) {
          for (int l = 0; l < 43; l++) {
            tmpNSFillDetMapPl1[i][j][k][l] = 0.;
            tmpNSFillDetMapMin1[i][j][k][l] = 0.;
          }
        }
      }
    }
 
    edm::Handle<HBHERecHitCollection> hbheNormal;
    iEvent.getByToken(tok_hbheNorm_, hbheNormal);
    if (!hbheNormal.isValid()) {
      edm::LogWarning("AnalyzerMB") << " hbheNormal failed ";
    } else {
      edm::LogInfo("AnalyzerMB") << " The size of the normal collection " << hbheNormal->size();
    }
 
    edm::Handle<HBHERecHitCollection> hbheNS;
    iEvent.getByToken(tok_hbheNoise_, hbheNS);
 
    if (!hbheNS.isValid()) {
      edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe"
                                    << " product! No HBHE MS ";
      return;
    }
 
    const HBHERecHitCollection HithbheNS = *(hbheNS.product());
    edm::LogInfo("AnalyzerMB") << " HBHE NS size of collection " << HithbheNS.size();
    hHBHEsize_vs_run->Fill(rnnum, (float)HithbheNS.size());
 
    if (HithbheNS.size() != 5184) {
      edm::LogWarning("AnalyzerMB") << " HBHE problem " << rnnum << " " << HithbheNS.size();
      //          return;
    }
    edm::Handle<HBHERecHitCollection> hbheMB;
    iEvent.getByToken(tok_hbhe_, hbheMB);
 
    if (!hbheMB.isValid()) {
      edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe"
                                    << " product! No HBHE MB";
      //   return ;
    }
 
    const HBHERecHitCollection HithbheMB = *(hbheMB.product());
    edm::LogInfo("AnalyzerMB") << " HBHE MB size of collection " << HithbheMB.size();
    if (HithbheMB.size() != 5184) {
      edm::LogWarning("AnalyzerMB") << " HBHE problem " << rnnum << " " << HithbheMB.size();
      //   return;
    }
 
    edm::Handle<HFRecHitCollection> hfNS;
    iEvent.getByToken(tok_hfNoise_, hfNS);
 
    if (!hfNS.isValid()) {
      edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf"
                                    << " product! No HF NS ";
      // return ;
    }
 
    const HFRecHitCollection HithfNS = *(hfNS.product());
    edm::LogInfo("AnalyzerMB") << " HFE NS size of collection " << HithfNS.size();
    hHFsize_vs_run->Fill(rnnum, (float)HithfNS.size());
    if (HithfNS.size() != 1728) {
      edm::LogWarning("AnalyzerMB") << " HF problem " << rnnum << " " << HithfNS.size();
      //    return;
    }
 
    edm::Handle<HFRecHitCollection> hfMB;
    iEvent.getByToken(tok_hf_, hfMB);
 
    if (!hfMB.isValid()) {
      edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf"
                                    << " product! No HF MB";
      // return ;
    }
 
    const HFRecHitCollection HithfMB = *(hfMB.product());
    edm::LogInfo("AnalyzerMB") << " HF MB size of collection " << HithfMB.size();
    if (HithfMB.size() != 1728) {
      edm::LogWarning("AnalyzerMB") << " HF problem " << rnnum << " " << HithfMB.size();
      //      return;
    }
 
    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);
 
      int mysu = ((hid).rawId() >> 25) & 0x7;
      if (hid.ieta() > 0) {
        theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
        theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
        theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
        theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
 
        tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
 
      } else {
        theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
        theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
        theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
        theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
 
        tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
      }
 
      if (hid.depth() == 1) {
        hbheNoiseE->Fill(energyhit);
 
        if (energyhit < -2.)
          edm::LogInfo("AnalyzerMB") << " Run " << rnnum << " ieta,iphi " << hid.ieta() << " " << hid.iphi()
                                     << energyhit;
 
        // if( hid.ieta() > 0 ) {
        //  hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit-noise_pl[hid.iphi()][hid.ieta()]);
        //  hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
        // } else {
        //  hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit-noise_min[hid.iphi()][abs(hid.ieta())]);
        //  hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
        // } // eta><0
 
      }  // depth=1
 
    }  // 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);
 
      int mysu = ((hid).rawId() >> 25) & 0x7;
      if (hid.ieta() > 0) {
        theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] += 1.;
        theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] += energyhit;
        theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] += pow(energyhit, 2);
        theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] += pow(energyhit, 4);
        float mydiff = energyhit - tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
 
        theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + mydiff;
        theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(mydiff, 2);
      } else {
        theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
        theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
        theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
        theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
 
        float mydiff = energyhit - tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
        theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + mydiff;
        theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(mydiff, 2);
      }
 
      if (hid.depth() == 1) {
        hbheSignalE->Fill(energyhit);
 
        if (hid.ieta() > 0) {
          hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
          hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit, 2));
        } else {
          hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
          hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit, 2));
        }  // eta><0
 
      }  // depth=1
 
    }  // 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
      //
      DetId id = (*hbheItr).detid();
      HcalDetId hid = HcalDetId(id);
 
      if (fabs(energyhit) > 40.)
        continue;
 
      int mysu = hid.subdetId();
      if (hid.ieta() > 0) {
        theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
        theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
        theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
        theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
 
        tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
 
      } else {
        theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
        theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
        theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
        theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
 
        tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
      }
 
      if (hid.depth() == 1) {
        hfNoiseE->Fill(energyhit);
 
        //if( hid.ieta() > 0 ) {
        // hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit-noise_pl[hid.iphi()][hid.ieta()]);
        // hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
        //} else {
        // hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit-noise_min[hid.iphi()][abs(hid.ieta())]);
        // hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
        //} // eta><0
 
      }  // depth=1
 
    }  // HBHE_NS
 
    // Signal part for HB HE
 
    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);
 
      int mysu = ((hid).rawId() >> 25) & 0x7;
      if (hid.ieta() > 0) {
        theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
        theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
        theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
        theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
 
        theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit -
            tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
        theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] +
            pow((energyhit - tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]), 2);
      } else {
        theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
        theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
        theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
        theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
            theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
 
        theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit -
            tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
        theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
            theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] +
            pow((energyhit - tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()]), 2);
      }
 
      if (hid.depth() == 1) {
        hfSignalE->Fill(energyhit);
 
        if (hid.ieta() > 0) {
          hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
          hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit, 2));
        } else {
          hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
          hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit, 2));
        }  // eta><0
 
      }  // depth=1
 
    }  // HF_MB
 
    edm::LogInfo("AnalyzerMB") << " Event is finished ";
  }

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), HcalDetId::depth(), edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), HcalObjRepresent::Fill(), dqmMemoryStats::float, edm::EventSetup::getData(), HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, mps_fire::i, HcalDetId::ieta(), iEvent, HcalDetId::iphi(), edm::HandleBase::isValid(), dqmiolumiharvest::j, dqmdumpme::k, cmsLHEtoEOSManager::l, nevent, nevent_run, funct::pow(), edm::Handle< T >::product(), DetId::rawId(), edm::SortedCollection< T, SORT >::size(), DetId::subdetId(), theDFFillDetMapMin1, theDFFillDetMapMin2, theDFFillDetMapPl1, theDFFillDetMapPl2, theMBFillDetMapMin0, theMBFillDetMapMin1, theMBFillDetMapMin2, theMBFillDetMapMin4, theMBFillDetMapPl0, theMBFillDetMapPl1, theMBFillDetMapPl2, theMBFillDetMapPl4, theNSFillDetMapMin0, theNSFillDetMapMin1, theNSFillDetMapMin2, theNSFillDetMapMin4, theNSFillDetMapPl0, theNSFillDetMapPl1, theNSFillDetMapPl2, theNSFillDetMapPl4, theRecalib, tok_hbhe_, tok_hbheNoise_, tok_hbheNorm_, tok_hf_, tok_hfNoise_, and tok_respCorr_.

beginJob()

void cms::Analyzer_minbias::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 88 of file Analyzer_minbias.cc.

                                  {
    //  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("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("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("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");
 
    edm::LogInfo("AnalyzerMB") << " Before ordering Histos ";
 
    char str0[32];
    char str1[32];
 
    char str10[32];
    char str11[32];
 
    int k = 0;
    nevent = 0;
    // Size of collections
 
    hHBHEsize_vs_run =
        fs->make<TH2F>("hHBHEsize_vs_run", "hHBHEsize_vs_run", 500, 111500., 112000., 6101, -100.5, 6000.5);
    hHFsize_vs_run = fs->make<TH2F>("hHFsize_vs_run", "hHFsize_vs_run", 500, 111500., 112000., 6101, -100.5, 6000.5);
 
    for (int i = 1; i < 73; i++) {
      for (int j = 1; j < 43; j++) {
        meannoise_pl[i][j] = 0.;
        meannoise_min[i][j] = 0.;
 
        //     for(int l=1;l<5;l++){
        k = i * 1000 + j;
        sprintf(str0, "mpl%d", k);
        sprintf(str1, "mmin%d", k);
 
        sprintf(str10, "vpl%d", k);
        sprintf(str11, "vmin%d", k);
        //      edm::LogInfo("AnalyzerMB")<<" "<<i<<" "<<j;
        if (j < 30) {
          // first order moment
          hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
          hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
 
          // second order moment
          hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 20.);
          hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 20.);
        } else {
          // HF
          // first order moment
          //   edm::LogInfo("AnalyzerMB")<<" "<<i<<" "<<j<<" "<<k;
          if (j < 40) {
            hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
            hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
            //
            // second order moment
            hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 40.);
            hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 40.);
          } else {
            hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
            hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
 
            // second order moment
            hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 120.);
            hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 120.);
          }
        }  // HE/HF boundary
           //     } // l
      }    // j
    }      // i
 
    hbheNoiseE = fs->make<TH1F>("hbheNoiseE", "hbheNoiseE", 320, -10., 10.);
    hfNoiseE = fs->make<TH1F>("hfNoiseE", "hfNoiseE", 320, -10., 10.);
    hbheSignalE = fs->make<TH1F>("hbheSignalE", "hbheSignalE", 320, -10., 10.);
    hfSignalE = fs->make<TH1F>("hfSignalE", "hfSignalE", 320, -10., 10.);
 
    edm::LogInfo("AnalyzerMB") << " After ordering Histos ";
 
    std::string ccc = "noise_0.dat";
 
    myout_hcal = new std::ofstream(ccc.c_str());
    if (!myout_hcal)
      edm::LogInfo("AnalyzerMB") << " Output file not open!!! ";
 
    //
    for (int i = 0; i < 5; i++) {
      for (int j = 0; j < 5; j++) {
        for (int k = 0; k < 73; k++) {
          for (int l = 0; l < 43; l++) {
            theMBFillDetMapPl0[i][j][k][l] = 0.;
            theMBFillDetMapPl1[i][j][k][l] = 0.;
            theMBFillDetMapPl2[i][j][k][l] = 0.;
            theMBFillDetMapPl4[i][j][k][l] = 0.;
 
            theMBFillDetMapMin0[i][j][k][l] = 0.;
            theMBFillDetMapMin1[i][j][k][l] = 0.;
            theMBFillDetMapMin2[i][j][k][l] = 0.;
            theMBFillDetMapMin4[i][j][k][l] = 0.;
 
            theNSFillDetMapPl0[i][j][k][l] = 0.;
            theNSFillDetMapPl1[i][j][k][l] = 0.;
            theNSFillDetMapPl2[i][j][k][l] = 0.;
            theNSFillDetMapPl4[i][j][k][l] = 0.;
 
            theNSFillDetMapMin0[i][j][k][l] = 0.;
            theNSFillDetMapMin1[i][j][k][l] = 0.;
            theNSFillDetMapMin2[i][j][k][l] = 0.;
            theNSFillDetMapMin4[i][j][k][l] = 0.;
 
            theDFFillDetMapPl0[i][j][k][l] = 0.;
            theDFFillDetMapPl1[i][j][k][l] = 0.;
            theDFFillDetMapPl2[i][j][k][l] = 0.;
            theDFFillDetMapMin0[i][j][k][l] = 0.;
            theDFFillDetMapMin1[i][j][k][l] = 0.;
            theDFFillDetMapMin2[i][j][k][l] = 0.;
          }
        }
      }
    }
 
    return;
  }

References depth, eta, fs, hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, mps_fire::i, ieta, iphi, dqmiolumiharvest::j, dqmdumpme::k, cmsLHEtoEOSManager::l, TFileService::make(), meannoise_min, meannoise_pl, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom4_MB, mom4_Noise, mydet, myout_hcal, mysubd, myTree, nevent, occup, phi, AlCaHLTBitMon_QueryRunRegistry::string, theDFFillDetMapMin0, theDFFillDetMapMin1, theDFFillDetMapMin2, theDFFillDetMapPl0, theDFFillDetMapPl1, theDFFillDetMapPl2, theMBFillDetMapMin0, theMBFillDetMapMin1, theMBFillDetMapMin2, theMBFillDetMapMin4, theMBFillDetMapPl0, theMBFillDetMapPl1, theMBFillDetMapPl2, theMBFillDetMapPl4, theNSFillDetMapMin0, theNSFillDetMapMin1, theNSFillDetMapMin2, theNSFillDetMapMin4, theNSFillDetMapPl0, theNSFillDetMapPl1, theNSFillDetMapPl2, and theNSFillDetMapPl4.

beginRun()

void cms::Analyzer_minbias::beginRun ( edm::Run const &  , edm::EventSetup const &    )

override

Definition at line 83 of file Analyzer_minbias.cc.

{ nevent_run = 0; }

References nevent_run.

endJob()

void cms::Analyzer_minbias::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 230 of file Analyzer_minbias.cc.

                                {
    int ii = 0;
 
    for (int i = 1; i < 5; i++) {
      for (int j = 1; j < 5; j++) {
        for (int k = 1; k < 73; k++) {
          for (int l = 1; l < 43; l++) {
            if (theMBFillDetMapPl0[i][j][k][l] > 0) {
              mom0_MB = theMBFillDetMapPl0[i][j][k][l];
              mom1_MB = theMBFillDetMapPl1[i][j][k][l];
              mom2_MB = theMBFillDetMapPl2[i][j][k][l];
              mom4_MB = theMBFillDetMapPl4[i][j][k][l];
              mom0_Noise = theNSFillDetMapPl0[i][j][k][l];
              mom1_Noise = theNSFillDetMapPl1[i][j][k][l];
              mom2_Noise = theNSFillDetMapPl2[i][j][k][l];
              mom4_Noise = theNSFillDetMapPl4[i][j][k][l];
              mom0_Diff = theDFFillDetMapPl0[i][j][k][l];
              mom1_Diff = theDFFillDetMapPl1[i][j][k][l];
              mom2_Diff = theDFFillDetMapPl2[i][j][k][l];
 
              mysubd = i;
              depth = j;
              ieta = l;
              iphi = k;
              edm::LogInfo("AnalyzerMB") << " Result Plus= " << mysubd << " " << ieta << " " << iphi << " mom0  "
                                         << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB;
              myTree->Fill();
              ii++;
            }  // Pl > 0
 
            if (theMBFillDetMapMin0[i][j][k][l] > 0) {
              mom0_MB = theMBFillDetMapMin0[i][j][k][l];
              mom1_MB = theMBFillDetMapMin1[i][j][k][l];
              mom2_MB = theMBFillDetMapMin2[i][j][k][l];
              mom4_MB = theMBFillDetMapMin4[i][j][k][l];
              mom0_Noise = theNSFillDetMapMin0[i][j][k][l];
              mom1_Noise = theNSFillDetMapMin1[i][j][k][l];
              mom2_Noise = theNSFillDetMapMin2[i][j][k][l];
              mom4_Noise = theNSFillDetMapMin4[i][j][k][l];
              mom0_Diff = theDFFillDetMapMin0[i][j][k][l];
              mom1_Diff = theDFFillDetMapMin1[i][j][k][l];
              mom2_Diff = theDFFillDetMapMin2[i][j][k][l];
 
              mysubd = i;
              depth = j;
              ieta = -1 * l;
              iphi = k;
              edm::LogInfo("AnalyzerMB") << " Result Minus= " << mysubd << " " << ieta << " " << iphi << " mom0  "
                                         << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB;
              myTree->Fill();
              ii++;
 
            }  // Min>0
          }    // ieta
        }      // iphi
      }        // depth
    }          //subd
 
    edm::LogInfo("AnalyzerMB") << " Number of cells " << ii;
 
    //  hOutputFile->Write();
    //  hOutputFile->cd();
    //  myTree->Write();
    //  hHBHEsize_vs_run->Write() ;
    //  hHFsize_vs_run->Write() ;
 
    for (int i = 1; i < 73; i++) {
      for (int j = 1; j < 43; j++) {
        hCalo1[i][j]->Write();
        hCalo2[i][j]->Write();
        hCalo1mom2[i][j]->Write();
        hCalo2mom2[i][j]->Write();
      }
    }
 
    //  hbheNoiseE->Write() ;
    //  hfNoiseE->Write() ;
    //  hbheSignalE->Write() ;
    //  hfSignalE->Write() ;
    //  hOutputFile->Close() ;
 
    edm::LogInfo("AnalyzerMB") << " File is closed ";
 
    return;
  }

References depth, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, mps_fire::i, ieta, cuy::ii, iphi, dqmiolumiharvest::j, dqmdumpme::k, cmsLHEtoEOSManager::l, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom4_MB, mom4_Noise, mysubd, myTree, theDFFillDetMapMin0, theDFFillDetMapMin1, theDFFillDetMapMin2, theDFFillDetMapPl0, theDFFillDetMapPl1, theDFFillDetMapPl2, theMBFillDetMapMin0, theMBFillDetMapMin1, theMBFillDetMapMin2, theMBFillDetMapMin4, theMBFillDetMapPl0, theMBFillDetMapPl1, theMBFillDetMapPl2, theMBFillDetMapPl4, theNSFillDetMapMin0, theNSFillDetMapMin1, theNSFillDetMapMin2, theNSFillDetMapMin4, theNSFillDetMapPl0, theNSFillDetMapPl1, theNSFillDetMapPl2, and theNSFillDetMapPl4.

Referenced by o2olib.O2ORunMgr::executeJob().

endRun()

void cms::Analyzer_minbias::endRun ( edm::Run const &  r, edm::EventSetup const &    )

override

Definition at line 84 of file Analyzer_minbias.cc.

                                                                     {
    edm::LogInfo("AnalyzerMB") << " Runnumber " << r.run() << " Nevents  " << nevent_run;
  }

References nevent_run, and alignCSCRings::r.

Member Data Documentation

depth

int cms::Analyzer_minbias::depth

private

Definition at line 98 of file Analyzer_minbias.h.

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

eta

float cms::Analyzer_minbias::eta

private

Definition at line 99 of file Analyzer_minbias.h.

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

fOutputFileName

std::string cms::Analyzer_minbias::fOutputFileName

private

Definition at line 74 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

fs

edm::Service<TFileService> cms::Analyzer_minbias::fs

private

Definition at line 81 of file Analyzer_minbias.h.

Referenced by beginJob().

hbheNoiseE

TH1F* cms::Analyzer_minbias::hbheNoiseE

private

Definition at line 89 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

hbheSignalE

TH1F* cms::Analyzer_minbias::hbheSignalE

private

Definition at line 90 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

hcalfile_

std::string cms::Analyzer_minbias::hcalfile_

private

Definition at line 75 of file Analyzer_minbias.h.

hCalo1

TH1F* cms::Analyzer_minbias::hCalo1[73][43]

private

Definition at line 85 of file Analyzer_minbias.h.

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

hCalo1mom2

TH1F* cms::Analyzer_minbias::hCalo1mom2[73][43]

private

Definition at line 87 of file Analyzer_minbias.h.

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

hCalo2

TH1F* cms::Analyzer_minbias::hCalo2[73][43]

private

Definition at line 86 of file Analyzer_minbias.h.

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

hCalo2mom2

TH1F* cms::Analyzer_minbias::hCalo2mom2[73][43]

private

Definition at line 88 of file Analyzer_minbias.h.

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

hfNoiseE

TH1F* cms::Analyzer_minbias::hfNoiseE

private

Definition at line 91 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

hfSignalE

TH1F* cms::Analyzer_minbias::hfSignalE

private

Definition at line 92 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

hHBHEsize_vs_run

TH2F* cms::Analyzer_minbias::hHBHEsize_vs_run

private

Definition at line 94 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

hHFsize_vs_run

TH2F* cms::Analyzer_minbias::hHFsize_vs_run

private

Definition at line 95 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

hOutputFile

TFile* cms::Analyzer_minbias::hOutputFile

private

Definition at line 83 of file Analyzer_minbias.h.

ieta

int cms::Analyzer_minbias::ieta

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

iphi

int cms::Analyzer_minbias::iphi

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

meannoise_min

double cms::Analyzer_minbias::meannoise_min[73][43]

private

Definition at line 106 of file Analyzer_minbias.h.

Referenced by beginJob().

meannoise_pl

double cms::Analyzer_minbias::meannoise_pl[73][43]

private

Definition at line 106 of file Analyzer_minbias.h.

Referenced by beginJob().

mom0_Diff

float cms::Analyzer_minbias::mom0_Diff

private

Definition at line 102 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom0_MB

float cms::Analyzer_minbias::mom0_MB

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom0_Noise

float cms::Analyzer_minbias::mom0_Noise

private

Definition at line 101 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom1_Diff

float cms::Analyzer_minbias::mom1_Diff

private

Definition at line 102 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom1_MB

float cms::Analyzer_minbias::mom1_MB

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom1_Noise

float cms::Analyzer_minbias::mom1_Noise

private

Definition at line 101 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom2_Diff

float cms::Analyzer_minbias::mom2_Diff

private

Definition at line 102 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom2_MB

float cms::Analyzer_minbias::mom2_MB

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom2_Noise

float cms::Analyzer_minbias::mom2_Noise

private

Definition at line 101 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom3_Diff

float cms::Analyzer_minbias::mom3_Diff

private

Definition at line 102 of file Analyzer_minbias.h.

mom3_MB

float cms::Analyzer_minbias::mom3_MB

private

Definition at line 100 of file Analyzer_minbias.h.

mom3_Noise

float cms::Analyzer_minbias::mom3_Noise

private

Definition at line 101 of file Analyzer_minbias.h.

mom4_Diff

float cms::Analyzer_minbias::mom4_Diff

private

Definition at line 102 of file Analyzer_minbias.h.

mom4_MB

float cms::Analyzer_minbias::mom4_MB

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mom4_Noise

float cms::Analyzer_minbias::mom4_Noise

private

Definition at line 101 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

mydet

int cms::Analyzer_minbias::mydet

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob().

myout_hcal

std::ofstream* cms::Analyzer_minbias::myout_hcal

private

Definition at line 76 of file Analyzer_minbias.h.

Referenced by beginJob().

mysubd

int cms::Analyzer_minbias::mysubd

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

myTree

TTree* cms::Analyzer_minbias::myTree

private

Definition at line 84 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

nevent

double cms::Analyzer_minbias::nevent

private

Definition at line 111 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

nevent_run

int cms::Analyzer_minbias::nevent_run

private

Definition at line 97 of file Analyzer_minbias.h.

Referenced by analyze(), beginRun(), and endRun().

noise_min

double cms::Analyzer_minbias::noise_min[73][43]

private

Definition at line 107 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

noise_pl

double cms::Analyzer_minbias::noise_pl[73][43]

private

Definition at line 107 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

occup

float cms::Analyzer_minbias::occup

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob().

phi

float cms::Analyzer_minbias::phi

private

Definition at line 99 of file Analyzer_minbias.h.

Referenced by Particle.Particle::__str__(), beginJob(), and ntupleDataFormat.Track::phiPull().

theDFFillDetMapMin0

double cms::Analyzer_minbias::theDFFillDetMapMin0[5][5][73][43]

private

Definition at line 135 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

theDFFillDetMapMin1

double cms::Analyzer_minbias::theDFFillDetMapMin1[5][5][73][43]

private

Definition at line 136 of file Analyzer_minbias.h.

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

theDFFillDetMapMin2

double cms::Analyzer_minbias::theDFFillDetMapMin2[5][5][73][43]

private

Definition at line 137 of file Analyzer_minbias.h.

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

theDFFillDetMapPl0

double cms::Analyzer_minbias::theDFFillDetMapPl0[5][5][73][43]

private

Definition at line 132 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

theDFFillDetMapPl1

double cms::Analyzer_minbias::theDFFillDetMapPl1[5][5][73][43]

private

Definition at line 133 of file Analyzer_minbias.h.

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

theDFFillDetMapPl2

double cms::Analyzer_minbias::theDFFillDetMapPl2[5][5][73][43]

private

Definition at line 134 of file Analyzer_minbias.h.

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

theMBFillDetMapMin0

double cms::Analyzer_minbias::theMBFillDetMapMin0[5][5][73][43]

private

Definition at line 117 of file Analyzer_minbias.h.

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

theMBFillDetMapMin1

double cms::Analyzer_minbias::theMBFillDetMapMin1[5][5][73][43]

private

Definition at line 118 of file Analyzer_minbias.h.

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

theMBFillDetMapMin2

double cms::Analyzer_minbias::theMBFillDetMapMin2[5][5][73][43]

private

Definition at line 119 of file Analyzer_minbias.h.

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

theMBFillDetMapMin4

double cms::Analyzer_minbias::theMBFillDetMapMin4[5][5][73][43]

private

Definition at line 120 of file Analyzer_minbias.h.

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

theMBFillDetMapPl0

double cms::Analyzer_minbias::theMBFillDetMapPl0[5][5][73][43]

private

Definition at line 112 of file Analyzer_minbias.h.

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

theMBFillDetMapPl1

double cms::Analyzer_minbias::theMBFillDetMapPl1[5][5][73][43]

private

Definition at line 113 of file Analyzer_minbias.h.

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

theMBFillDetMapPl2

double cms::Analyzer_minbias::theMBFillDetMapPl2[5][5][73][43]

private

Definition at line 114 of file Analyzer_minbias.h.

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

theMBFillDetMapPl4

double cms::Analyzer_minbias::theMBFillDetMapPl4[5][5][73][43]

private

Definition at line 115 of file Analyzer_minbias.h.

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

theNSFillDetMapMin0

double cms::Analyzer_minbias::theNSFillDetMapMin0[5][5][73][43]

private

Definition at line 127 of file Analyzer_minbias.h.

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

theNSFillDetMapMin1

double cms::Analyzer_minbias::theNSFillDetMapMin1[5][5][73][43]

private

Definition at line 128 of file Analyzer_minbias.h.

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

theNSFillDetMapMin2

double cms::Analyzer_minbias::theNSFillDetMapMin2[5][5][73][43]

private

Definition at line 129 of file Analyzer_minbias.h.

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

theNSFillDetMapMin4

double cms::Analyzer_minbias::theNSFillDetMapMin4[5][5][73][43]

private

Definition at line 130 of file Analyzer_minbias.h.

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

theNSFillDetMapPl0

double cms::Analyzer_minbias::theNSFillDetMapPl0[5][5][73][43]

private

Definition at line 122 of file Analyzer_minbias.h.

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

theNSFillDetMapPl1

double cms::Analyzer_minbias::theNSFillDetMapPl1[5][5][73][43]

private

Definition at line 123 of file Analyzer_minbias.h.

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

theNSFillDetMapPl2

double cms::Analyzer_minbias::theNSFillDetMapPl2[5][5][73][43]

private

Definition at line 124 of file Analyzer_minbias.h.

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

theNSFillDetMapPl4

double cms::Analyzer_minbias::theNSFillDetMapPl4[5][5][73][43]

private

Definition at line 125 of file Analyzer_minbias.h.

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

theRecalib

bool cms::Analyzer_minbias::theRecalib

private

Definition at line 153 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

tok_data_

edm::EDGetTokenT<FEDRawDataCollection> cms::Analyzer_minbias::tok_data_

private

Definition at line 78 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

tok_gtRec_

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> cms::Analyzer_minbias::tok_gtRec_

private

Definition at line 148 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> cms::Analyzer_minbias::tok_hbhe_

private

Definition at line 139 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

tok_hbheNoise_

edm::EDGetTokenT<HBHERecHitCollection> cms::Analyzer_minbias::tok_hbheNoise_

private

Definition at line 143 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

tok_hbheNorm_

edm::EDGetTokenT<HBHERecHitCollection> cms::Analyzer_minbias::tok_hbheNorm_

private

Definition at line 149 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

tok_hf_

edm::EDGetTokenT<HFRecHitCollection> cms::Analyzer_minbias::tok_hf_

private

Definition at line 141 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

tok_hfNoise_

edm::EDGetTokenT<HFRecHitCollection> cms::Analyzer_minbias::tok_hfNoise_

private

Definition at line 145 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

tok_ho_

edm::EDGetTokenT<HORecHitCollection> cms::Analyzer_minbias::tok_ho_

private

Definition at line 140 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

tok_hoNoise_

edm::EDGetTokenT<HORecHitCollection> cms::Analyzer_minbias::tok_hoNoise_

private

Definition at line 144 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

tok_respCorr_

edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> cms::Analyzer_minbias::tok_respCorr_

private

Definition at line 151 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().