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
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)
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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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 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

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 332 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(), mathSSE::return(), 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 92 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, mathSSE::return(), 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 85 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 238 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, mathSSE::return(), 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 88 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 150 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 147 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 138 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 142 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 148 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 140 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 144 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 139 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 143 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().