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
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 const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (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 const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
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_

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)
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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 59 of file Analyzer_minbias.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 43 of file Analyzer_minbias.cc.

References fOutputFileName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, 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_, and tok_hoNoise_.

                                                                   {
    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"));

    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.;
      }
    }
    //
    //
  }

cms::Analyzer_minbias::~Analyzer_minbias ( )

override

Definition at line 76 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

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 319 of file Analyzer_minbias.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), DEFINE_FWK_MODULE, HcalDetId::depth(), edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), HcalObjRepresent::Fill(), objects.autophobj::float, edm::EventSetup::get(), edm::Event::getAllStableProvenance(), edm::Event::getByToken(), HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, mps_fire::i, HcalDetId::ieta(), HcalDetId::iphi(), edm::HandleBase::isValid(), gen::k, checklumidiff::l, nevent, nevent_run, funct::pow(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), DetId::rawId(), edm::Event::run(), 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_, and tok_hfNoise_.

                                                                                    {
    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) {
      // Radek:
      edm::ESHandle<HcalRespCorrs> recalibCorrs;
      iSetup.get<HcalRespCorrsRcd>().get("recalibrate", recalibCorrs);
      myRecalib = recalibCorrs.product();
      // end
    }  // 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 ";
  }

void cms::Analyzer_minbias::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 86 of file Analyzer_minbias.cc.

References depth, eta, fs, hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, mps_fire::i, ieta, iphi, gen::k, checklumidiff::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.

                                  {
    //  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[15];
    char str1[15];

    char str10[15];
    char str11[15];

    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;
  }

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

override

Definition at line 81 of file Analyzer_minbias.cc.

References nevent_run.

{ nevent_run = 0; }

void cms::Analyzer_minbias::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 228 of file Analyzer_minbias.cc.

References depth, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, mps_fire::i, ieta, cuy::ii, iphi, gen::k, checklumidiff::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().

                                {
    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;
  }

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

override

Definition at line 82 of file Analyzer_minbias.cc.

References nevent_run, and edm::RunBase::run().

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

Member Data Documentation

int cms::Analyzer_minbias::depth

private

Definition at line 96 of file Analyzer_minbias.h.

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

float cms::Analyzer_minbias::eta

private

Definition at line 97 of file Analyzer_minbias.h.

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

std::string cms::Analyzer_minbias::fOutputFileName

private

Definition at line 72 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 79 of file Analyzer_minbias.h.

Referenced by beginJob().

TH1F* cms::Analyzer_minbias::hbheNoiseE

private

Definition at line 87 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

TH1F* cms::Analyzer_minbias::hbheSignalE

private

Definition at line 88 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

std::string cms::Analyzer_minbias::hcalfile_

private

Definition at line 73 of file Analyzer_minbias.h.

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

private

Definition at line 83 of file Analyzer_minbias.h.

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

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

private

Definition at line 85 of file Analyzer_minbias.h.

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

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

private

Definition at line 84 of file Analyzer_minbias.h.

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

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

private

Definition at line 86 of file Analyzer_minbias.h.

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

TH1F* cms::Analyzer_minbias::hfNoiseE

private

Definition at line 89 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

TH1F* cms::Analyzer_minbias::hfSignalE

private

Definition at line 90 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

TH2F* cms::Analyzer_minbias::hHBHEsize_vs_run

private

Definition at line 92 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

TH2F* cms::Analyzer_minbias::hHFsize_vs_run

private

Definition at line 93 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

TFile* cms::Analyzer_minbias::hOutputFile

private

Definition at line 81 of file Analyzer_minbias.h.

int cms::Analyzer_minbias::ieta

private

Definition at line 96 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

int cms::Analyzer_minbias::iphi

private

Definition at line 96 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

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

private

Definition at line 104 of file Analyzer_minbias.h.

Referenced by beginJob().

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

private

Definition at line 104 of file Analyzer_minbias.h.

Referenced by beginJob().

float cms::Analyzer_minbias::mom0_Diff

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom0_MB

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom0_Noise

private

Definition at line 99 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom1_Diff

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom1_MB

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom1_Noise

private

Definition at line 99 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom2_Diff

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom2_MB

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom2_Noise

private

Definition at line 99 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom3_Diff

private

Definition at line 100 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom3_MB

private

Definition at line 98 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom3_Noise

private

Definition at line 99 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom4_Diff

private

Definition at line 100 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom4_MB

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom4_Noise

private

Definition at line 99 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

int cms::Analyzer_minbias::mydet

private

Definition at line 96 of file Analyzer_minbias.h.

Referenced by beginJob().

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

private

Definition at line 74 of file Analyzer_minbias.h.

Referenced by beginJob().

int cms::Analyzer_minbias::mysubd

private

Definition at line 96 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

TTree* cms::Analyzer_minbias::myTree

private

Definition at line 82 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

double cms::Analyzer_minbias::nevent

private

Definition at line 109 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

int cms::Analyzer_minbias::nevent_run

private

Definition at line 95 of file Analyzer_minbias.h.

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

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

private

Definition at line 105 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 105 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

float cms::Analyzer_minbias::occup

private

Definition at line 98 of file Analyzer_minbias.h.

Referenced by beginJob().

float cms::Analyzer_minbias::phi

private

Definition at line 97 of file Analyzer_minbias.h.

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

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

private

Definition at line 133 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

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

private

Definition at line 134 of file Analyzer_minbias.h.

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

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

private

Definition at line 135 of file Analyzer_minbias.h.

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

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

private

Definition at line 130 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

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

private

Definition at line 131 of file Analyzer_minbias.h.

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

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

private

Definition at line 132 of file Analyzer_minbias.h.

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

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

private

Definition at line 115 of file Analyzer_minbias.h.

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

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

private

Definition at line 116 of file Analyzer_minbias.h.

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

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

private

Definition at line 117 of file Analyzer_minbias.h.

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

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

private

Definition at line 118 of file Analyzer_minbias.h.

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

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

private

Definition at line 110 of file Analyzer_minbias.h.

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

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

private

Definition at line 111 of file Analyzer_minbias.h.

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

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

private

Definition at line 112 of file Analyzer_minbias.h.

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

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

private

Definition at line 113 of file Analyzer_minbias.h.

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

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

private

Definition at line 125 of file Analyzer_minbias.h.

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

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

private

Definition at line 126 of file Analyzer_minbias.h.

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

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

private

Definition at line 127 of file Analyzer_minbias.h.

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

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

private

Definition at line 128 of file Analyzer_minbias.h.

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

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

private

Definition at line 120 of file Analyzer_minbias.h.

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

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

private

Definition at line 121 of file Analyzer_minbias.h.

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

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

private

Definition at line 122 of file Analyzer_minbias.h.

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

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

private

Definition at line 123 of file Analyzer_minbias.h.

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

bool cms::Analyzer_minbias::theRecalib

private

Definition at line 149 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 76 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 146 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 137 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 141 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 147 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 139 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 143 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 138 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 142 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().