#include <Analyzer_minbias.h>

Inheritance diagram for cms::Analyzer_minbias:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

	Analyzer_minbias (const edm::ParameterSet &)

virtual void	beginJob ()

virtual void	beginRun (const edm::Run &r, const edm::EventSetup &iSetup)

virtual void	endJob ()

virtual void	endRun (const edm::Run &r, const edm::EventSetup &iSetup)

	~Analyzer_minbias ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
int	depth

float	eta

std::string	fOutputFileName

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 55 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, 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_.

                                                                    {
     // 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 ( )

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	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Definition at line 336 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(), reco::return(), edm::Event::run(), edm::SortedCollection< T, SORT >::size(), 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=0;
     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);
       // theGeometry->getGeometry(detId)->getPosition()).eta()
 
       //         edm::LogInfo("AnalyzerMB")<<hid.ieta()<<" "<<hid.iphi()<<" "<<hid.depth()<<" "<<(pG->getGeometry(id)->getPosition()).eta()<<" "<<(pG->getGeometry(id)->getPosition()).phi()
       //                      <<" "<<(pG->getGeometry(id)->getPosition()).perp()<<" "<<(pG->getGeometry(id)->getPosition()).z();
      
       if(fabs(energyhit) > 40. ) continue;
  
       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) {
     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 )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 90 of file Analyzer_minbias.cc.

References depth, eta, fOutputFileName, hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, hOutputFile, mps_fire::i, ieta, iphi, gen::k, checklumidiff::l, 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, reco::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 = new 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 = new TH2F("hHBHEsize_vs_run","hHBHEsize_vs_run",500,111500.,112000.,6101,-100.5,6000.5);
     hHFsize_vs_run = new 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] = new TH1F(str0, "h0", 320, -10., 10.);
       hCalo2[i][j] = new TH1F(str1, "h1", 320, -10., 10.);
       
       // second order moment
       hCalo1mom2[i][j] = new TH1F(str10, "h10", 320, 0., 20.);
       hCalo2mom2[i][j] = new TH1F(str11, "h11", 320, 0., 20.);
     } else {
       // HF
       // first order moment
       //   edm::LogInfo("AnalyzerMB")<<" "<<i<<" "<<j<<" "<<k;
       if(j < 40) {
         hCalo1[i][j] = new TH1F(str0, "h0", 320, -10., 10.);
         hCalo2[i][j] = new TH1F(str1, "h1", 320, -10., 10.);
         //
         // second order moment
         hCalo1mom2[i][j] = new TH1F(str10, "h10", 320, 0., 40.);
         hCalo2mom2[i][j] = new TH1F(str11, "h11", 320, 0., 40.);
       } else {
         hCalo1[i][j] = new TH1F(str0,"h0" , 320, -10., 10.);
         hCalo2[i][j] = new TH1F(str1, "h1", 320, -10., 10.);
 
         // second order moment
         hCalo1mom2[i][j] = new TH1F(str10, "h10", 320, 0., 120.);
         hCalo2mom2[i][j] = new TH1F(str11, "h11", 320, 0., 120.);
 
       }
     } // HE/HF boundary
     //     } // l
       } // j
     } // i
 
 
     hbheNoiseE = new TH1F("hbheNoiseE","hbheNoiseE", 320, -10., 10.);
     hfNoiseE = new TH1F("hfNoiseE","hfNoiseE", 320, -10., 10.);
     hbheSignalE = new TH1F("hbheSignalE","hbheSignalE", 320, -10., 10.);
     hfSignalE = new 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	(	const edm::Run &	r,
		const edm::EventSetup &	iSetup
	)

virtual

Definition at line 83 of file Analyzer_minbias.cc.

References nevent_run.

                                                                                {
     nevent_run = 0;
   }

void cms::Analyzer_minbias::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 237 of file Analyzer_minbias.cc.

References depth, hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, hOutputFile, 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, reco::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	(	const edm::Run &	r,
		const edm::EventSetup &	iSetup
	)

virtual

Definition at line 86 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 91 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 92 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 68 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias(), and beginJob().

TH1F* cms::Analyzer_minbias::hbheNoiseE

private

Definition at line 82 of file Analyzer_minbias.h.

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

TH1F* cms::Analyzer_minbias::hbheSignalE

private

Definition at line 83 of file Analyzer_minbias.h.

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

std::string cms::Analyzer_minbias::hcalfile_

private

Definition at line 69 of file Analyzer_minbias.h.

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

private

Definition at line 78 of file Analyzer_minbias.h.

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

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

private

Definition at line 80 of file Analyzer_minbias.h.

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

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

private

Definition at line 79 of file Analyzer_minbias.h.

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

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

private

Definition at line 81 of file Analyzer_minbias.h.

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

TH1F* cms::Analyzer_minbias::hfNoiseE

private

Definition at line 84 of file Analyzer_minbias.h.

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

TH1F* cms::Analyzer_minbias::hfSignalE

private

Definition at line 85 of file Analyzer_minbias.h.

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

TH2F* cms::Analyzer_minbias::hHBHEsize_vs_run

private

Definition at line 87 of file Analyzer_minbias.h.

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

TH2F* cms::Analyzer_minbias::hHFsize_vs_run

private

Definition at line 88 of file Analyzer_minbias.h.

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

TFile* cms::Analyzer_minbias::hOutputFile

private

Definition at line 76 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

int cms::Analyzer_minbias::ieta

private

Definition at line 91 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

int cms::Analyzer_minbias::iphi

private

Definition at line 91 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

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

private

Definition at line 99 of file Analyzer_minbias.h.

Referenced by beginJob().

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

private

Definition at line 99 of file Analyzer_minbias.h.

Referenced by beginJob().

float cms::Analyzer_minbias::mom0_Diff

private

Definition at line 95 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom0_MB

private

Definition at line 93 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom0_Noise

private

Definition at line 94 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom1_Diff

private

Definition at line 95 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom1_MB

private

Definition at line 93 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom1_Noise

private

Definition at line 94 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom2_Diff

private

Definition at line 95 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom2_MB

private

Definition at line 93 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom2_Noise

private

Definition at line 94 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom3_Diff

private

Definition at line 95 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom3_MB

private

Definition at line 93 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom3_Noise

private

Definition at line 94 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom4_Diff

private

Definition at line 95 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom4_MB

private

Definition at line 93 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

float cms::Analyzer_minbias::mom4_Noise

private

Definition at line 94 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

int cms::Analyzer_minbias::mydet

private

Definition at line 91 of file Analyzer_minbias.h.

Referenced by beginJob().

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

private

Definition at line 70 of file Analyzer_minbias.h.

Referenced by beginJob().

int cms::Analyzer_minbias::mysubd

private

Definition at line 91 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

TTree* cms::Analyzer_minbias::myTree

private

Definition at line 77 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

double cms::Analyzer_minbias::nevent

private

Definition at line 104 of file Analyzer_minbias.h.

Referenced by analyze(), and beginJob().

int cms::Analyzer_minbias::nevent_run

private

Definition at line 90 of file Analyzer_minbias.h.

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

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

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 100 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

float cms::Analyzer_minbias::occup

private

Definition at line 93 of file Analyzer_minbias.h.

Referenced by beginJob().

float cms::Analyzer_minbias::phi

private

Definition at line 92 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 128 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

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

private

Definition at line 129 of file Analyzer_minbias.h.

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

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

private

Definition at line 130 of file Analyzer_minbias.h.

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

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

private

Definition at line 125 of file Analyzer_minbias.h.

Referenced by beginJob(), and endJob().

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

private

Definition at line 126 of file Analyzer_minbias.h.

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

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

private

Definition at line 127 of file Analyzer_minbias.h.

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

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

private

Definition at line 110 of file Analyzer_minbias.h.

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

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

private

Definition at line 111 of file Analyzer_minbias.h.

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

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

private

Definition at line 112 of file Analyzer_minbias.h.

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

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

private

Definition at line 113 of file Analyzer_minbias.h.

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

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

private

Definition at line 105 of file Analyzer_minbias.h.

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

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

private

Definition at line 106 of file Analyzer_minbias.h.

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

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

private

Definition at line 107 of file Analyzer_minbias.h.

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

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

private

Definition at line 108 of file Analyzer_minbias.h.

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

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

private

Definition at line 120 of file Analyzer_minbias.h.

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

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

private

Definition at line 121 of file Analyzer_minbias.h.

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

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

private

Definition at line 122 of file Analyzer_minbias.h.

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

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

private

Definition at line 123 of file Analyzer_minbias.h.

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

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

private

Definition at line 115 of file Analyzer_minbias.h.

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

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

private

Definition at line 116 of file Analyzer_minbias.h.

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

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

private

Definition at line 117 of file Analyzer_minbias.h.

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

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

private

Definition at line 118 of file Analyzer_minbias.h.

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

bool cms::Analyzer_minbias::theRecalib

private

Definition at line 145 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 72 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 142 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 133 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 137 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 143 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 135 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 139 of file Analyzer_minbias.h.

Referenced by analyze(), and Analyzer_minbias().

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

private

Definition at line 134 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

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

private

Definition at line 138 of file Analyzer_minbias.h.

Referenced by Analyzer_minbias().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation