#include <EcalCosmicsHists.h>

Inheritance diagram for EcalCosmicsHists:

Public Member Functions
	EcalCosmicsHists (const edm::ParameterSet &)

	~EcalCosmicsHists () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

std::vector< bool >	determineTriggers (const edm::Event &, const edm::EventSetup &eventSetup)

void	endJob () override

void	initHists (int)

std::string	intToString (int num)

Private Attributes
TH1F *	allFedsE2Hist_

TH2F *	allFedsE2vsE1Hist_

TH1F *	allFedsenergyHighHist_

TH1F *	allFedsenergyHist_

TH1F *	allFedsenergyOnlyHighHist_

TH2F *	allFedsenergyvsE1Hist_

TH1F *	allFedsFreqTimeHist_

TH2F *	allFedsFreqTimeVsEtaHist_

TH2F *	allFedsFreqTimeVsEtaTTHist_

TH2F *	allFedsFreqTimeVsPhiHist_

TH2F *	allFedsFreqTimeVsPhiTTHist_

TH1F *	allFedsFrequencyHist_

TH1F *	allFedsHist_

TH1F *	allFedsiEtaProfileHist_

TH1F *	allFedsiPhiProfileHist_

TH1F *	allFedsNumXtalsInClusterHist_

TH3F *	allFedsOccupancyHighEnergyHist_

TH1F *	allFedsTimingEbmBottomHist_

TH1F *	allFedsTimingEbmHist_

TH1F *	allFedsTimingEbmTopHist_

TH1F *	allFedsTimingEbpBottomHist_

TH1F *	allFedsTimingEbpHist_

TH1F *	allFedsTimingEbpTopHist_

TH2F *	allFedsTimingEtaHist_

TH1F *	allFedsTimingHist_

TH1F *	allFedsTimingHistCSC_

TH1F *	allFedsTimingHistDT_

TH1F *	allFedsTimingHistECAL_

TH1F *	allFedsTimingHistEcalMuon_

TH1F *	allFedsTimingHistHCAL_

TH1F *	allFedsTimingHistRPC_

TH2F *	allFedsTimingLMHist_

TH2F *	allFedsTimingLMHistCSC_

TH2F *	allFedsTimingLMHistDT_

TH2F *	allFedsTimingLMHistECAL_

TH2F *	allFedsTimingLMHistHCAL_

TH2F *	allFedsTimingLMHistRPC_

TH2F *	allFedsTimingPhiEbmHist_

TH2F *	allFedsTimingPhiEbpHist_

TH3F *	allFedsTimingPhiEtaHist_

TH3F *	allFedsTimingPhiEtaHistCSC_

TH3F *	allFedsTimingPhiEtaHistDT_

TH3F *	allFedsTimingPhiEtaHistECAL_

TH3F *	allFedsTimingPhiEtaHistHCAL_

TH3F *	allFedsTimingPhiEtaHistRPC_

TH2F *	allFedsTimingPhiHist_

TH3F *	allFedsTimingTTHist_

TH3F *	allFedsTimingTTHistCSC_

TH3F *	allFedsTimingTTHistDT_

TH3F *	allFedsTimingTTHistECAL_

TH3F *	allFedsTimingTTHistHCAL_

TH3F *	allFedsTimingTTHistRPC_

TH2F *	allFedsTimingVsAmpHist_

TH2F *	allFedsTimingVsFreqHist_

TH2F *	allOccupancy_

TH2F *	allOccupancyCoarse_

TH2F *	allOccupancyCoarseCSC_

TH2F *	allOccupancyCoarseDT_

TH2F *	allOccupancyCoarseECAL_

TH2F *	allOccupancyCoarseExclusiveCSC_

TH2F *	allOccupancyCoarseExclusiveDT_

TH2F *	allOccupancyCoarseExclusiveECAL_

TH2F *	allOccupancyCoarseExclusiveHCAL_

TH2F *	allOccupancyCoarseExclusiveRPC_

TH2F *	allOccupancyCoarseHCAL_

TH2F *	allOccupancyCoarseRPC_

TH2F *	allOccupancyCSC_

TH2F *	allOccupancyDT_

TH2F *	allOccupancyECAL_

TH2F *	allOccupancyExclusiveCSC_

TH2F *	allOccupancyExclusiveDT_

TH2F *	allOccupancyExclusiveECAL_

TH2F *	allOccupancyExclusiveHCAL_

TH2F *	allOccupancyExclusiveRPC_

TH2F *	allOccupancyHCAL_

TH2F *	allOccupancyHighEnergy_

TH2F *	allOccupancyHighEnergyCoarse_

TH2F *	allOccupancyRPC_

TH2F *	allOccupancySingleXtal_

edm::InputTag	barrelClusterCollection_

int	cosmicCounter_

int	cosmicCounterEB_

int	cosmicCounterEEM_

int	cosmicCounterEEP_

int	cosmicCounterTopBottom_

TH1F *	dccBXErrorByFEDHist_

TH2F *	dccErrorVsBxHist_

TH2F *	dccEventVsBxHist_

TH1F *	dccOrbitErrorByFEDHist_

TH1F *	dccRuntypeErrorByFEDHist_

TH1F *	dccRuntypeHist_

TH2F *	deltaEtaDeltaPhiHist_

TH1F *	deltaEtaHist_

TH1F *	deltaPhiHist_

TH1F *	deltaRHist_

const EcalElectronicsMapping *	ecalElectronicsMap_

edm::InputTag	ecalRawDataColl_

edm::InputTag	ecalRecHitCollectionEB_

edm::InputTag	ecalRecHitCollectionEE_

TH2F *	EEM_AllOccupancy_

TH2F *	EEM_AllOccupancyCoarse_

TH1F *	EEM_energySingleXtalHist_

TH1F *	EEM_FedsE2Hist_

TH2F *	EEM_FedsE2vsE1Hist_

TH1F *	EEM_FedsenergyHighHist_

TH1F *	EEM_FedsenergyHist_

TH1F *	EEM_FedsenergyOnlyHighHist_

TH2F *	EEM_FedsenergyvsE1Hist_

TH1F *	EEM_FedsNumXtalsInClusterHist_

TH1F *	EEM_FedsSeedEnergyHist_

TH1F *	EEM_FedsTimingHist_

TH1F *	EEM_FedsTimingHistCSC_

TH1F *	EEM_FedsTimingHistDT_

TH1F *	EEM_FedsTimingHistECAL_

TH1F *	EEM_FedsTimingHistHCAL_

TH1F *	EEM_FedsTimingHistRPC_

TH3F *	EEM_FedsTimingTTHist_

TH3F *	EEM_FedsTimingTTHistCSC_

TH3F *	EEM_FedsTimingTTHistDT_

TH3F *	EEM_FedsTimingTTHistECAL_

TH3F *	EEM_FedsTimingTTHistHCAL_

TH3F *	EEM_FedsTimingTTHistRPC_

TH2F *	EEM_FedsTimingVsAmpHist_

TH1F *	EEM_numberofBCinSC_

TH1F *	EEM_numberofCosmicsHist_

TH1F *	EEM_NumXtalsInClusterHist_

TH2F *	EEM_numxtalsVsEnergy_

TH2F *	EEM_numxtalsVsHighEnergy_

TH2F *	EEM_OccupancyCoarseCSC_

TH2F *	EEM_OccupancyCoarseDT_

TH2F *	EEM_OccupancyCoarseECAL_

TH2F *	EEM_OccupancyCoarseExclusiveCSC_

TH2F *	EEM_OccupancyCoarseExclusiveDT_

TH2F *	EEM_OccupancyCoarseExclusiveECAL_

TH2F *	EEM_OccupancyCoarseExclusiveHCAL_

TH2F *	EEM_OccupancyCoarseExclusiveRPC_

TH2F *	EEM_OccupancyCoarseHCAL_

TH2F *	EEM_OccupancyCoarseRPC_

TH2F *	EEM_OccupancyCSC_

TH2F *	EEM_OccupancyDT_

TH2F *	EEM_OccupancyECAL_

TH2F *	EEM_OccupancyExclusiveCSC_

TH2F *	EEM_OccupancyExclusiveDT_

TH2F *	EEM_OccupancyExclusiveECAL_

TH2F *	EEM_OccupancyExclusiveHCAL_

TH2F *	EEM_OccupancyExclusiveRPC_

TH2F *	EEM_OccupancyHCAL_

TH2F *	EEM_OccupancyHighEnergy_

TH2F *	EEM_OccupancyHighEnergyCoarse_

TH2F *	EEM_OccupancyRPC_

TH2F *	EEM_OccupancySingleXtal_

TH1F *	EEM_triggerExclusiveHist_

TH1F *	EEM_triggerHist_

TH2F *	EEP_AllOccupancy_

TH2F *	EEP_AllOccupancyCoarse_

TH1F *	EEP_energySingleXtalHist_

TH1F *	EEP_FedsE2Hist_

TH2F *	EEP_FedsE2vsE1Hist_

TH1F *	EEP_FedsenergyHighHist_

TH1F *	EEP_FedsenergyHist_

TH1F *	EEP_FedsenergyOnlyHighHist_

TH2F *	EEP_FedsenergyvsE1Hist_

TH1F *	EEP_FedsNumXtalsInClusterHist_

TH1F *	EEP_FedsSeedEnergyHist_

TH1F *	EEP_FedsTimingHist_

TH1F *	EEP_FedsTimingHistCSC_

TH1F *	EEP_FedsTimingHistDT_

TH1F *	EEP_FedsTimingHistECAL_

TH1F *	EEP_FedsTimingHistHCAL_

TH1F *	EEP_FedsTimingHistRPC_

TH3F *	EEP_FedsTimingTTHist_

TH3F *	EEP_FedsTimingTTHistCSC_

TH3F *	EEP_FedsTimingTTHistDT_

TH3F *	EEP_FedsTimingTTHistECAL_

TH3F *	EEP_FedsTimingTTHistHCAL_

TH3F *	EEP_FedsTimingTTHistRPC_

TH2F *	EEP_FedsTimingVsAmpHist_

TH1F *	EEP_numberofBCinSC_

TH1F *	EEP_numberofCosmicsHist_

TH1F *	EEP_NumXtalsInClusterHist_

TH2F *	EEP_numxtalsVsEnergy_

TH2F *	EEP_numxtalsVsHighEnergy_

TH2F *	EEP_OccupancyCoarseCSC_

TH2F *	EEP_OccupancyCoarseDT_

TH2F *	EEP_OccupancyCoarseECAL_

TH2F *	EEP_OccupancyCoarseExclusiveCSC_

TH2F *	EEP_OccupancyCoarseExclusiveDT_

TH2F *	EEP_OccupancyCoarseExclusiveECAL_

TH2F *	EEP_OccupancyCoarseExclusiveHCAL_

TH2F *	EEP_OccupancyCoarseExclusiveRPC_

TH2F *	EEP_OccupancyCoarseHCAL_

TH2F *	EEP_OccupancyCoarseRPC_

TH2F *	EEP_OccupancyCSC_

TH2F *	EEP_OccupancyDT_

TH2F *	EEP_OccupancyECAL_

TH2F *	EEP_OccupancyExclusiveCSC_

TH2F *	EEP_OccupancyExclusiveDT_

TH2F *	EEP_OccupancyExclusiveECAL_

TH2F *	EEP_OccupancyExclusiveHCAL_

TH2F *	EEP_OccupancyExclusiveRPC_

TH2F *	EEP_OccupancyHCAL_

TH2F *	EEP_OccupancyHighEnergy_

TH2F *	EEP_OccupancyHighEnergyCoarse_

TH2F *	EEP_OccupancyRPC_

TH2F *	EEP_OccupancySingleXtal_

TH1F *	EEP_triggerExclusiveHist_

TH1F *	EEP_triggerHist_

edm::InputTag	endcapClusterCollection_

TH1F *	energySingleXtalHist_

EcalFedMap *	fedMap_

std::map< int, TH2F * >	FEDsAndDCCRuntypeVsBxHists_

std::map< int, TH1F * >	FEDsAndE2Hists_

std::map< int, TH2F * >	FEDsAndE2vsE1Hists_

std::map< int, TH1F * >	FEDsAndenergyHists_

std::map< int, TH2F * >	FEDsAndenergyvsE1Hists_

std::map< int, TH1F * >	FEDsAndFrequencyHists_

std::map< int, TH1F * >	FEDsAndHists_

std::map< int, TH1F * >	FEDsAndiEtaProfileHists_

std::map< int, TH1F * >	FEDsAndiPhiProfileHists_

std::map< int, TH1F * >	FEDsAndNumXtalsInClusterHists_

std::map< int, TH2F * >	FEDsAndOccupancyHists_

std::map< int, TH1F * >	FEDsAndTimingHists_

std::map< int, TH2F * >	FEDsAndTimingVsAmpHists_

std::map< int, TH2F * >	FEDsAndTimingVsEtaHists_

std::map< int, TH2F * >	FEDsAndTimingVsFreqHists_

std::map< int, TH2F * >	FEDsAndTimingVsModuleHists_

std::map< int, TH2F * >	FEDsAndTimingVsPhiHists_

TFile *	file

std::string	fileName_

TH1F *	hcalEnergy_HBHE_

TH1F *	hcalEnergy_HF_

TH1F *	hcalEnergy_HO_

TH2F *	hcalHEHBecalEB_

TH3F *	HighEnergy_0tracks_occu3D

TH3F *	HighEnergy_0tracks_occu3DXtal

TH3F *	HighEnergy_100GeV_occu3D

TH2F *	HighEnergy_100GeV_occuCoarse

TH3F *	HighEnergy_1tracks_occu3D

TH3F *	HighEnergy_1tracks_occu3DXtal

TH3F *	HighEnergy_2GeV_occu3D

TH2F *	HighEnergy_2GeV_occuCoarse

TH3F *	HighEnergy_2tracks_occu3D

TH3F *	HighEnergy_2tracks_occu3DXtal

TH1F *	HighEnergy_bestSeed

TH2F *	HighEnergy_bestSeedOccupancy

TH3F *	HighEnergy_energy3D

TH2F *	HighEnergy_energyNumXtal

TH1F *	HighEnergy_numClusHighEn

TH1F *	HighEnergy_numRecoTrackBarrel

TH1F *	HighEnergy_NumXtal

TH2F *	HighEnergy_NumXtalFedId

TH2F *	HighEnergy_NumXtaliphi

TH1F *	HighEnergy_ratioClusters

TH1F *	HighEnergy_TracksAngle

TH1F *	HighEnergy_TracksAngleTopBottom

double	histRangeMax_

double	histRangeMin_

std::vector< int >	l1Accepts_

edm::InputTag	l1GMTReadoutRecTag_

edm::InputTag	l1GTReadoutRecTag_

std::vector< std::string >	l1Names_

double	minHighEnergy_

double	minRecHitAmpEB_

double	minRecHitAmpEE_

double	minTimingAmpEB_

double	minTimingAmpEE_

double *	modEtaBins

int	naiveEvtNum_

TH1F *	numberofBCinSC_

TH2F *	numberofBCinSCphi_

TH1F *	numberofCosmicsHist_

TH1F *	numberofCosmicsHistEB_

TH1F *	numberofCosmicsTopBottomHist_

TH1F *	numberofCosmicsWTrackHist_

TH1F *	numberofCrossedEcalIdsHist_

TH1F *	numberofGoodEvtFreq_

int	numTimingBins_

TH1F *	NumXtalsInClusterHist_

TH2F *	numxtalsVsEnergy_

TH2F *	numxtalsVsHighEnergy_

TH1F *	ratioAssocClustersHist_

TH1F *	ratioAssocTracksHist_

bool	runInFileName_

int	runNum_

TH1F *	runNumberHist_

double	runTimeLength_

TH2F *	seedTrackEtaHist_

TH2F *	seedTrackPhiHist_

double	startTime_

TH2F *	trackAssoc_muonsEcal_

TrackDetectorAssociator	trackAssociator_

TrackAssociatorParameters	trackParameters_

TH1F *	triggerExclusiveHist_

TH1F *	triggerHist_

TH2F *	TrueBCOccupancy_

TH2F *	TrueBCOccupancyCoarse_

TH2F *	TrueOccupancy_

TH2F *	TrueOccupancyCoarse_

double *	ttEtaBins

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 &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

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

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

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 68 of file EcalCosmicsHists.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 56 of file EcalCosmicsHists.cc.

References allFedsHist_, edm::EDConsumerBase::consumesCollector(), cosmicCounter_, cosmicCounterEB_, cosmicCounterEEM_, cosmicCounterEEP_, cosmicCounterTopBottom_, fedMap_, edm::ParameterSet::getParameter(), histRangeMax_, histRangeMin_, TrackAssociatorParameters::loadParameters(), naiveEvtNum_, trackAssociator_, trackParameters_, and TrackDetectorAssociator::useDefaultPropagator().

     : ecalRawDataColl_(iConfig.getParameter<edm::InputTag>("ecalRawDataColl")),
       ecalRecHitCollectionEB_(iConfig.getParameter<edm::InputTag>("ecalRecHitCollectionEB")),
       ecalRecHitCollectionEE_(iConfig.getParameter<edm::InputTag>("ecalRecHitCollectionEE")),
       barrelClusterCollection_(iConfig.getParameter<edm::InputTag>("barrelClusterCollection")),
       endcapClusterCollection_(iConfig.getParameter<edm::InputTag>("endcapClusterCollection")),
       l1GTReadoutRecTag_(iConfig.getUntrackedParameter<std::string>("L1GlobalReadoutRecord", "gtDigis")),
       l1GMTReadoutRecTag_(iConfig.getUntrackedParameter<std::string>("L1GlobalMuonReadoutRecord", "gtDigis")),
       runNum_(-1),
       histRangeMax_(iConfig.getUntrackedParameter<double>("histogramMaxRange", 1.8)),
       histRangeMin_(iConfig.getUntrackedParameter<double>("histogramMinRange", 0.0)),
       minTimingAmpEB_(iConfig.getUntrackedParameter<double>("MinTimingAmpEB", 0.100)),
       minTimingAmpEE_(iConfig.getUntrackedParameter<double>("MinTimingAmpEE", 0.100)),
       minRecHitAmpEB_(iConfig.getUntrackedParameter<double>("MinRecHitAmpEB", 0.027)),
       minRecHitAmpEE_(iConfig.getUntrackedParameter<double>("MinRecHitAmpEE", 0.18)),
       minHighEnergy_(iConfig.getUntrackedParameter<double>("MinHighEnergy", 2.0)),
       fileName_(iConfig.getUntrackedParameter<std::string>("fileName", std::string("ecalCosmicHists"))),
       runInFileName_(iConfig.getUntrackedParameter<bool>("runInFileName", true)),
       startTime_(iConfig.getUntrackedParameter<double>("TimeStampStart", 1215107133.)),
       runTimeLength_(iConfig.getUntrackedParameter<double>("TimeStampLength", 3.)),
       numTimingBins_(iConfig.getUntrackedParameter<int>("TimeStampBins", 1800)) {
   naiveEvtNum_ = 0;
   cosmicCounter_ = 0;
   cosmicCounterEB_ = 0;
   cosmicCounterEEM_ = 0;
   cosmicCounterEEP_ = 0;
   cosmicCounterTopBottom_ = 0;
 
   // TrackAssociator parameters
   edm::ParameterSet trkParameters = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
   edm::ConsumesCollector iC = consumesCollector();
   trackParameters_.loadParameters(trkParameters, iC);
   trackAssociator_.useDefaultPropagator();
 
   string title1 = "Seed Energy for All Feds; Seed Energy (GeV)";
   string name1 = "SeedEnergyAllFEDs";
   int numBins = 200;  //(int)round(histRangeMax_-histRangeMin_)+1;
   allFedsHist_ = new TH1F(name1.c_str(), title1.c_str(), numBins, histRangeMin_, histRangeMax_);
 
   fedMap_ = new EcalFedMap();
 }

EcalCosmicsHists::~EcalCosmicsHists ( )

override

Definition at line 98 of file EcalCosmicsHists.cc.

98 {}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 105 of file EcalCosmicsHists.cc.

                                                                                   {
   bool hasEndcapClusters = true;
   int ievt = iEvent.id().event();
 
   edm::Handle<reco::SuperClusterCollection> bscHandle;
   edm::Handle<reco::SuperClusterCollection> escHandle;
 
   naiveEvtNum_++;
 
   //LogDebug("EcalCosmicsHists")<< "  My Event: " << naiveEvtNum_ << " " << iEvent.id().run() << " " << iEvent.id().event() << " " << iEvent.time().value();
   //LogDebug("EcalCosmicsHists")<< "Timestamp: " << iEvent.id().run() << " " << iEvent.id().event() << " " << iEvent.time().value();
 
   // check DB payload
   edm::ESHandle<EcalADCToGeVConstant> pAgc;
   iSetup.get<EcalADCToGeVConstantRcd>().get(pAgc);
   const EcalADCToGeVConstant* agc = pAgc.product();
   if (naiveEvtNum_ <= 1) {
     LogWarning("EcalCosmicsHists") << "Global EB ADC->GeV scale: " << agc->getEBValue() << " GeV/ADC count";
     LogWarning("EcalCosmicsHists") << "Global EE ADC->GeV scale: " << agc->getEEValue() << " GeV/ADC count";
   }
   //float adcEBconst = agc->getEBValue();
   //float adcEEconst = agc->getEEValue();
 
   //
 
   //===================TIMESTAMP INFORMTION=================================
   // Code added to get the time in seconds
   unsigned int timeStampLow = (0xFFFFFFFF & iEvent.time().value());
   unsigned int timeStampHigh = (iEvent.time().value() >> 32);
   double rawEventTime = (double)(timeStampHigh) + ((double)(timeStampLow) / 1000000.);
   double eventTime = rawEventTime - startTime_;  //Notice the subtraction of the "starttime"
   //std::cout << "Event Time " << eventTime << " High " <<timeStampHigh<< " low"<<timeStampLow <<" value " <<iEvent.time().value() << std::endl;
   //========================================================================
 
   iEvent.getByLabel(barrelClusterCollection_, bscHandle);
   if (!(bscHandle.isValid())) {
     LogWarning("EcalCosmicsHists") << barrelClusterCollection_ << " not available";
     return;
   }
   LogDebug("EcalCosmicsHists") << "event " << ievt;
 
   iEvent.getByLabel(endcapClusterCollection_, escHandle);
   if (!(escHandle.isValid())) {
     LogWarning("EcalCosmicsHists") << endcapClusterCollection_ << " not available";
     hasEndcapClusters = false;
     //return;
   }
 
   Handle<EcalRecHitCollection> hits;
   iEvent.getByLabel(ecalRecHitCollectionEB_, hits);
   if (!(hits.isValid())) {
     LogWarning("EcalCosmicsHists") << ecalRecHitCollectionEB_ << " not available";
     return;
   }
   Handle<EcalRecHitCollection> hitsEE;
   iEvent.getByLabel(ecalRecHitCollectionEE_, hitsEE);
   if (!(hitsEE.isValid())) {
     LogWarning("EcalCosmicsHists") << ecalRecHitCollectionEE_ << " not available";
     return;
   }
 
   Handle<EcalRawDataCollection> DCCHeaders;
   iEvent.getByLabel(ecalRawDataColl_, DCCHeaders);
   if (!DCCHeaders.isValid())
     LogWarning("EcalCosmicsHists") << "DCC headers not available";
 
   //make the bx histos right here
   //TODO: Right now we are filling histos for errors...
   int orbit = -100;
   int bx = -100;
   int runType = -100;
 
   for (EcalRawDataCollection::const_iterator headerItr = DCCHeaders->begin(); headerItr != DCCHeaders->end();
        ++headerItr) {
     headerItr->getEventSettings();
     int myorbit = headerItr->getOrbit();
     int mybx = headerItr->getBX();
     int myRunType = headerItr->getRunType();
     int FEDid = headerItr->fedId();
     TH2F* dccRuntypeHist = FEDsAndDCCRuntypeVsBxHists_[FEDid];
     if (dccRuntypeHist == nullptr) {
       initHists(FEDid);
       dccRuntypeHist = FEDsAndDCCRuntypeVsBxHists_[FEDid];
     }
     dccRuntypeHist->Fill(mybx, myRunType);
 
     if (bx == -100) {
       bx = mybx;
     } else if (bx != mybx) {
       LogWarning("EcalCosmicsHists") << "This header has a conflicting bx OTHERS were " << bx << " here " << mybx;
       dccBXErrorByFEDHist_->Fill(headerItr->fedId());
       if (bx != -100) {
         dccErrorVsBxHist_->Fill(bx, 0);
       }
     }
 
     if (runType == -100) {
       runType = myRunType;
     } else if (runType != myRunType) {
       LogWarning("EcalCosmicsHists") << "This header has a conflicting runType OTHERS were " << bx << " here " << mybx;
       dccRuntypeErrorByFEDHist_->Fill(headerItr->fedId());
       if (bx != -100)
         dccErrorVsBxHist_->Fill(bx, 2);
     }
 
     if (orbit == -100) {
       orbit = myorbit;
     } else if (orbit != myorbit) {
       LogWarning("EcalCosmicsHists") << "This header has a conflicting orbit; OTHERS were " << orbit << " here "
                                      << myorbit;
       dccOrbitErrorByFEDHist_->Fill(headerItr->fedId());
       if (bx != -100)
         dccErrorVsBxHist_->Fill(bx, 1);
     }
   }
   dccEventVsBxHist_->Fill(bx, runType);
   dccRuntypeHist_->Fill(runType);
 
   std::vector<bool> l1Triggers = determineTriggers(iEvent, iSetup);
   bool isEcalL1 = l1Triggers[4];
   bool isHCALL1 = l1Triggers[3];
   bool isRPCL1 = l1Triggers[2];
   bool isCSCL1 = l1Triggers[1];
   bool isDTL1 = l1Triggers[0];
 
   if (runNum_ == -1) {
     runNum_ = iEvent.id().run();
   }
 
   int numberOfCosmics = 0;
   int numberOfCosmicsEB = 0;
   int numberOfCosmicsEEP = 0;
   int numberOfCosmicsEEM = 0;
   int numberOfCosmicsTop = 0;
   int numberOfCosmicsBottom = 0;
   int numberOfHighEClusters = 0;
   //int eventnum = iEvent.id().event();
   std::vector<EBDetId> seeds;
 
   //++++++++++++++++++BEGIN LOOP OVER EB SUPERCLUSTERS+++++++++++++++++++++++++//
 
   const reco::SuperClusterCollection* clusterCollection_p = bscHandle.product();
   for (reco::SuperClusterCollection::const_iterator clus = clusterCollection_p->begin();
        clus != clusterCollection_p->end();
        ++clus) {
     double energy = clus->energy();
     double phi = clus->phi();
     double eta = clus->eta();
     double time = -1000.0;
     double ampli = 0.;
     double secondMin = 0.;
     double secondTime = -1000.;
     int numXtalsinCluster = 0;
 
     EBDetId maxDet;
     EBDetId secDet;
 
     numberofBCinSC_->Fill(clus->clustersSize());
     numberofBCinSCphi_->Fill(phi, clus->clustersSize());
 
     for (reco::CaloCluster_iterator bclus = (clus->clustersBegin()); bclus != (clus->clustersEnd()); ++bclus) {
       double cphi = (*bclus)->phi();
       double ceta = (*bclus)->eta();
       TrueBCOccupancy_->Fill(cphi, ceta);
       TrueBCOccupancyCoarse_->Fill(cphi, ceta);
     }
 
     std::vector<std::pair<DetId, float> > clusterDetIds = clus->hitsAndFractions();  //get these from the cluster
     for (std::vector<std::pair<DetId, float> >::const_iterator detitr = clusterDetIds.begin();
          detitr != clusterDetIds.end();
          ++detitr) {
       //Here I use the "find" on a digi collection... I have been warned...
       if ((*detitr).first.det() != DetId::Ecal) {
         std::cout << " det is " << (*detitr).first.det() << std::endl;
         continue;
       }
       if ((*detitr).first.subdetId() != EcalBarrel) {
         std::cout << " subdet is " << (*detitr).first.subdetId() << std::endl;
         continue;
       }
       EcalRecHitCollection::const_iterator thishit = hits->find((*detitr).first);
       if (thishit == hits->end())
         continue;
       //The checking above should no longer be needed...... as only those in the cluster would already have rechits..
 
       EcalRecHit myhit = (*thishit);
 
       double thisamp = myhit.energy();
       if (thisamp > minRecHitAmpEB_) {
         numXtalsinCluster++;
       }
       if (thisamp > secondMin) {
         secondMin = thisamp;
         secondTime = myhit.time();
         secDet = (EBDetId)(*detitr).first;
       }
       if (secondMin > ampli) {
         std::swap(ampli, secondMin);
         std::swap(time, secondTime);
         std::swap(maxDet, secDet);
       }
     }
 
     double fullnumXtalsinCluster = clusterDetIds.size();
 
     float E2 = ampli + secondMin;
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId((EBDetId)maxDet);
     int FEDid = 600 + elecId.dccId();
 
     numberOfCosmics++;
     numberOfCosmicsEB++;
 
     //Set some more values
     seeds.push_back(maxDet);
     int ieta = maxDet.ieta();
     int iphi = maxDet.iphi();
     int ietaSM = maxDet.ietaSM();
     int iphiSM = maxDet.iphiSM();
 
     int LM = ecalElectronicsMap_->getLMNumber(maxDet);  //FIX ME
 
     // top and bottom clusters
     if (iphi > 0 && iphi < 180) {
       numberOfCosmicsTop++;
     } else {
       numberOfCosmicsBottom++;
     }
 
     // fill the proper hist
     TH1F* uRecHist = FEDsAndHists_[FEDid];
     TH1F* E2uRecHist = FEDsAndE2Hists_[FEDid];
     TH1F* energyuRecHist = FEDsAndenergyHists_[FEDid];
     TH1F* timingHist = FEDsAndTimingHists_[FEDid];
     TH1F* freqHist = FEDsAndFrequencyHists_[FEDid];
     TH1F* iphiProfileHist = FEDsAndiPhiProfileHists_[FEDid];
     TH1F* ietaProfileHist = FEDsAndiEtaProfileHists_[FEDid];
     TH2F* timingHistVsFreq = FEDsAndTimingVsFreqHists_[FEDid];
     TH2F* timingHistVsAmp = FEDsAndTimingVsAmpHists_[FEDid];
     TH2F* E2vsE1uRecHist = FEDsAndE2vsE1Hists_[FEDid];
     TH2F* energyvsE1uRecHist = FEDsAndenergyvsE1Hists_[FEDid];
     TH2F* occupHist = FEDsAndOccupancyHists_[FEDid];
     TH2F* timingHistVsPhi = FEDsAndTimingVsPhiHists_[FEDid];
     TH2F* timingHistVsModule = FEDsAndTimingVsModuleHists_[FEDid];
 
     if (uRecHist == nullptr) {
       initHists(FEDid);
       uRecHist = FEDsAndHists_[FEDid];
       E2uRecHist = FEDsAndE2Hists_[FEDid];
       energyuRecHist = FEDsAndenergyHists_[FEDid];
       timingHist = FEDsAndTimingHists_[FEDid];
       freqHist = FEDsAndFrequencyHists_[FEDid];
       timingHistVsFreq = FEDsAndTimingVsFreqHists_[FEDid];
       timingHistVsAmp = FEDsAndTimingVsAmpHists_[FEDid];
       iphiProfileHist = FEDsAndiPhiProfileHists_[FEDid];
       ietaProfileHist = FEDsAndiEtaProfileHists_[FEDid];
       E2vsE1uRecHist = FEDsAndE2vsE1Hists_[FEDid];
       energyvsE1uRecHist = FEDsAndenergyvsE1Hists_[FEDid];
       occupHist = FEDsAndOccupancyHists_[FEDid];
       timingHistVsPhi = FEDsAndTimingVsPhiHists_[FEDid];
       timingHistVsModule = FEDsAndTimingVsModuleHists_[FEDid];
     }
 
     uRecHist->Fill(ampli);
     E2uRecHist->Fill(E2);
     E2vsE1uRecHist->Fill(ampli, E2);
     energyuRecHist->Fill(energy);
     energyvsE1uRecHist->Fill(ampli, energy);
     allFedsHist_->Fill(ampli);
     allFedsE2Hist_->Fill(E2);
     allFedsenergyHist_->Fill(energy);
     allFedsenergyHighHist_->Fill(energy);
     allFedsE2vsE1Hist_->Fill(ampli, E2);
     allFedsenergyvsE1Hist_->Fill(ampli, energy);
     freqHist->Fill(naiveEvtNum_);
     iphiProfileHist->Fill(iphi);
     ietaProfileHist->Fill(ieta);
 
     allFedsFrequencyHist_->Fill(naiveEvtNum_);
     allFedsiPhiProfileHist_->Fill(iphi);
     allFedsiEtaProfileHist_->Fill(ieta);
     allOccupancy_->Fill(iphi, ieta);
     TrueOccupancy_->Fill(phi, eta);
     allOccupancyCoarse_->Fill(iphi, ieta);
     TrueOccupancyCoarse_->Fill(phi, eta);
     allFedsNumXtalsInClusterHist_->Fill(numXtalsinCluster);
     NumXtalsInClusterHist_->Fill(fullnumXtalsinCluster);
     numxtalsVsEnergy_->Fill(energy, numXtalsinCluster);
     numxtalsVsHighEnergy_->Fill(energy, numXtalsinCluster);
 
     //Fill the hists for the time stamp information
     allFedsFreqTimeVsPhiHist_->Fill(iphi, eventTime);
     allFedsFreqTimeVsPhiTTHist_->Fill(iphi, eventTime);
     allFedsFreqTimeVsEtaHist_->Fill(eventTime, ieta);
     allFedsFreqTimeVsEtaTTHist_->Fill(eventTime, ieta);
     //end time stamp hists
 
     occupHist->Fill(ietaSM, iphiSM);
     if (fullnumXtalsinCluster == 1) {
       allOccupancySingleXtal_->Fill(iphi, ieta);
       energySingleXtalHist_->Fill(energy);
     }
 
     // Exclusive trigger plots
 
     if (isEcalL1 && !isDTL1 && !isRPCL1 && !isCSCL1 && !isHCALL1) {
       allOccupancyExclusiveECAL_->Fill(iphi, ieta);
       allOccupancyCoarseExclusiveECAL_->Fill(iphi, ieta);
       if (ampli > minTimingAmpEB_) {
         allFedsTimingHistECAL_->Fill(time);
         allFedsTimingPhiEtaHistECAL_->Fill(iphi, ieta, time);
         allFedsTimingTTHistECAL_->Fill(iphi, ieta, time);
         allFedsTimingLMHistECAL_->Fill(LM, time);
       }
       triggerExclusiveHist_->Fill(0);
     }
 
     if (!isEcalL1 && !isDTL1 && !isRPCL1 && !isCSCL1 && isHCALL1) {
       allOccupancyExclusiveHCAL_->Fill(iphi, ieta);
       allOccupancyCoarseExclusiveHCAL_->Fill(iphi, ieta);
       if (ampli > minTimingAmpEB_) {
         allFedsTimingHistHCAL_->Fill(time);
         allFedsTimingPhiEtaHistHCAL_->Fill(iphi, ieta, time);
         allFedsTimingTTHistHCAL_->Fill(iphi, ieta, time);
         allFedsTimingLMHistHCAL_->Fill(LM, time);
       }
       triggerExclusiveHist_->Fill(1);
     }
 
     if (!isEcalL1 && isDTL1 && !isRPCL1 && !isCSCL1 && !isHCALL1) {
       allOccupancyExclusiveDT_->Fill(iphi, ieta);
       allOccupancyCoarseExclusiveDT_->Fill(iphi, ieta);
       if (ampli > minTimingAmpEB_) {
         allFedsTimingHistDT_->Fill(time);
         allFedsTimingPhiEtaHistDT_->Fill(iphi, ieta, time);
         allFedsTimingTTHistDT_->Fill(iphi, ieta, time);
         allFedsTimingLMHistDT_->Fill(LM, time);
       }
       triggerExclusiveHist_->Fill(2);
     }
 
     if (!isEcalL1 && !isDTL1 && isRPCL1 && !isCSCL1 && !isHCALL1) {
       allOccupancyExclusiveRPC_->Fill(iphi, ieta);
       allOccupancyCoarseExclusiveRPC_->Fill(iphi, ieta);
       if (ampli > minTimingAmpEB_) {
         allFedsTimingHistRPC_->Fill(time);
         allFedsTimingPhiEtaHistRPC_->Fill(iphi, ieta, time);
         allFedsTimingTTHistRPC_->Fill(iphi, ieta, time);
         allFedsTimingLMHistRPC_->Fill(LM, time);
       }
       triggerExclusiveHist_->Fill(3);
     }
 
     if (!isEcalL1 && !isDTL1 && !isRPCL1 && isCSCL1 && !isHCALL1) {
       allOccupancyExclusiveCSC_->Fill(iphi, ieta);
       allOccupancyCoarseExclusiveCSC_->Fill(iphi, ieta);
       if (ampli > minTimingAmpEB_) {
         allFedsTimingHistCSC_->Fill(time);
         allFedsTimingPhiEtaHistCSC_->Fill(iphi, ieta, time);
         allFedsTimingTTHistCSC_->Fill(iphi, ieta, time);
         allFedsTimingLMHistCSC_->Fill(LM, time);
       }
       triggerExclusiveHist_->Fill(4);
     }
 
     // Inclusive trigger plots
 
     if (isEcalL1) {
       triggerHist_->Fill(0);
       allOccupancyECAL_->Fill(iphi, ieta);
       allOccupancyCoarseECAL_->Fill(iphi, ieta);
     }
     if (isHCALL1) {
       triggerHist_->Fill(1);
       allOccupancyHCAL_->Fill(iphi, ieta);
       allOccupancyCoarseHCAL_->Fill(iphi, ieta);
     }
     if (isDTL1) {
       triggerHist_->Fill(2);
       allOccupancyDT_->Fill(iphi, ieta);
       allOccupancyCoarseDT_->Fill(iphi, ieta);
     }
     if (isRPCL1) {
       triggerHist_->Fill(3);
       allOccupancyRPC_->Fill(iphi, ieta);
       allOccupancyCoarseRPC_->Fill(iphi, ieta);
     }
     if (isCSCL1) {
       triggerHist_->Fill(4);
       allOccupancyCSC_->Fill(iphi, ieta);
       allOccupancyCoarseCSC_->Fill(iphi, ieta);
     }
 
     // Fill histo for Ecal+muon coincidence
     if (isEcalL1 && (isCSCL1 || isRPCL1 || isDTL1) && !isHCALL1)
       allFedsTimingHistEcalMuon_->Fill(time);
 
     if (ampli > minTimingAmpEB_) {
       timingHist->Fill(time);
       timingHistVsFreq->Fill(time, naiveEvtNum_);
       timingHistVsAmp->Fill(time, ampli);
       allFedsTimingHist_->Fill(time);
       allFedsTimingVsAmpHist_->Fill(time, ampli);
       allFedsTimingVsFreqHist_->Fill(time, naiveEvtNum_);
       timingHistVsPhi->Fill(time, iphiSM);
       timingHistVsModule->Fill(time, ietaSM);
       allFedsTimingPhiHist_->Fill(iphi, time);
       allFedsTimingPhiEtaHist_->Fill(iphi, ieta, time);
       allFedsTimingTTHist_->Fill(iphi, ieta, time);
       allFedsTimingLMHist_->Fill(LM, time);
       if (FEDid >= 610 && FEDid <= 627)
         allFedsTimingPhiEbmHist_->Fill(iphi, time);
       if (FEDid >= 628 && FEDid <= 645)
         allFedsTimingPhiEbpHist_->Fill(iphi, time);
 
       if (FEDid >= 610 && FEDid <= 627)
         allFedsTimingEbmHist_->Fill(time);
       if (FEDid >= 628 && FEDid <= 645)
         allFedsTimingEbpHist_->Fill(time);
       if (FEDid >= 613 && FEDid <= 616)
         allFedsTimingEbmTopHist_->Fill(time);
       if (FEDid >= 631 && FEDid <= 634)
         allFedsTimingEbpTopHist_->Fill(time);
       if (FEDid >= 622 && FEDid <= 625)
         allFedsTimingEbmBottomHist_->Fill(time);
       if (FEDid >= 640 && FEDid <= 643)
         allFedsTimingEbpBottomHist_->Fill(time);
     }
 
     // *** High Energy Clusters Analysis ** //
 
     if (energy > minHighEnergy_) {
       LogInfo("EcalCosmicsHists") << "High energy event " << iEvent.id().run() << " : " << iEvent.id().event() << " "
                                   << naiveEvtNum_ << " : " << energy << " " << numXtalsinCluster << " : " << iphi << " "
                                   << ieta << " : " << isEcalL1 << isHCALL1 << isDTL1 << isRPCL1 << isCSCL1;
 
       numberOfHighEClusters++;
       allOccupancyHighEnergy_->Fill(iphi, ieta);
       allOccupancyHighEnergyCoarse_->Fill(iphi, ieta);
       allFedsOccupancyHighEnergyHist_->Fill(iphi, ieta, energy);
       allFedsenergyOnlyHighHist_->Fill(energy);
 
       HighEnergy_2GeV_occuCoarse->Fill(iphi, ieta);
       HighEnergy_2GeV_occu3D->Fill(iphi, ieta, energy);
 
       HighEnergy_NumXtal->Fill(fullnumXtalsinCluster);
       HighEnergy_NumXtalFedId->Fill(FEDid, fullnumXtalsinCluster);
       HighEnergy_NumXtaliphi->Fill(iphi, fullnumXtalsinCluster);
       HighEnergy_energy3D->Fill(iphi, ieta, energy);
       HighEnergy_energyNumXtal->Fill(fullnumXtalsinCluster, energy);
 
       if (energy > 100.0) {
         LogInfo("EcalCosmicsHists") << "Very high energy event " << iEvent.id().run() << " : " << iEvent.id().event()
                                     << " " << naiveEvtNum_ << " : " << energy << " " << numXtalsinCluster << " : "
                                     << iphi << " " << ieta << " : " << isEcalL1 << isHCALL1 << isDTL1 << isRPCL1
                                     << isCSCL1;
 
         HighEnergy_100GeV_occuCoarse->Fill(iphi, ieta);
         HighEnergy_100GeV_occu3D->Fill(iphi, ieta, energy);
       }
     }
 
     // *** end of High Energy Clusters analysis *** //
 
   }  //++++++++++++++++++END LOOP OVER EB SUPERCLUSTERS+++++++++++++++++++++++//
 
   //+++++++++++++++++++LOOP OVER ENDCAP EE CLUSTERS++++++++++++++++++++//
 
   if (hasEndcapClusters) {
     clusterCollection_p = escHandle.product();
     for (reco::SuperClusterCollection::const_iterator clus = clusterCollection_p->begin();
          clus != clusterCollection_p->end();
          ++clus) {
       double energy = clus->energy();
       //double phi    = clus->phi();
       //double eta    = clus->eta();
       double time = -1000.0;
       double ampli = 0.;
       double secondMin = 0.;
       double secondTime = -1000.;
       int clusSize = clus->clustersSize();
       int numXtalsinCluster = 0;
 
       EEDetId maxDet;
       EEDetId secDet;
       //LogInfo("EcalCosmicsHists") << "Here is what we initialized the maxDet to: " << maxDet;
 
       //      for (reco::basicCluster_iterator bclus = (clus->clustersBegin()); bclus != (clus->clustersEnd()); ++bclus) {
       //    //double cphi = (*bclus)->phi();
       //    //double ceta = (*bclus)->eta();
       //         //TODO: extend histos to EE
       //         //TrueBCOccupancy_->Fill(cphi,ceta);
       //         //TrueBCOccupancyCoarse_->Fill(cphi,ceta);
       //       }
 
       std::vector<std::pair<DetId, float> > clusterDetIds = clus->hitsAndFractions();  //get these from the cluster
       for (std::vector<std::pair<DetId, float> >::const_iterator detitr = clusterDetIds.begin();
            detitr != clusterDetIds.end();
            ++detitr) {
         //LogInfo("EcalCosmicsHists") << " Here is the DetId inside the cluster: " << (EEDetId)(*detitr);
         //Here I use the "find" on a digi collection... I have been warned...
 
         if ((*detitr).first.det() != DetId::Ecal) {
           LogError("EcalCosmicsHists") << " det is " << (*detitr).first.det();
           continue;
         }
         if ((*detitr).first.subdetId() != EcalEndcap) {
           LogError("EcalCosmicsHists") << " subdet is " << (*detitr).first.subdetId();
           continue;
         }
 
         EcalRecHitCollection::const_iterator thishit = hitsEE->find((*detitr).first);
 
         if (thishit == hitsEE->end()) {
           LogInfo("EcalCosmicsHists") << " WARNING: EEDetId not found in the RecHit collection!";
           continue;
         }
         // The checking above should no longer be needed......
         // as only those in the cluster would already have rechits..
 
         EcalRecHit myhit = (*thishit);
 
         //LogInfo("EcalCosmicsHists") << " Found hit for DetId: " << (EEDetId)(*detitr);
         double thisamp = myhit.energy();
         if (thisamp > minRecHitAmpEE_) {
           numXtalsinCluster++;
         }
         if (thisamp > secondMin) {
           secondMin = thisamp;
           secondTime = myhit.time();
           secDet = (EEDetId)(*detitr).first;
         }
         if (secondMin > ampli) {
           std::swap(ampli, secondMin);
           std::swap(time, secondTime);
           std::swap(maxDet, secDet);
         }
 
         //LogInfo("EcalCosmicsHists") << "maxDetId is now: " << (EEDetId)(maxDet);
       }
 
       double fullnumXtalsinCluster = clusterDetIds.size();
 
       float E2 = ampli + secondMin;
 
       ecalElectronicsMap_->getElectronicsId((EEDetId)maxDet);
       //int FEDid = 600+elecId.dccId();
 
       //Set some more values
       //TODO: need to fix the seeds vector to be DetId or have another one for EE
       //seeds.push_back(maxDet);
 
       int ix = maxDet.ix();
       int iy = maxDet.iy();
       int iz = maxDet.zside();
 
       //      LogWarning("EcalCosmicsHists") << "EE cluster (x,y,z) : ( "
       //                     << ix << " , " << iy << " , " << iz
       //                     << " ) " << std::endl;
 
       if (!EEDetId::validDetId(ix, iy, iz)) {
         LogWarning("EcalCosmicsHists") << "INVALID EE DetId !!!" << endl;
         return;
       }
 
       numberOfCosmics++;
       if (iz < 0) {
         numberOfCosmicsEEM++;
       } else {
         numberOfCosmicsEEP++;
       }
 
       //int LM = ecalElectronicsMap_->getLMNumber(maxDet) ;//FIX ME
 
       //TODO: extend histos to EE
       //TH1F* uRecHist = FEDsAndHists_[FEDid];
       //TH1F* E2uRecHist = FEDsAndE2Hists_[FEDid];
       //TH1F* energyuRecHist = FEDsAndenergyHists_[FEDid];
       //TH1F* timingHist = FEDsAndTimingHists_[FEDid];
       //TH1F* freqHist = FEDsAndFrequencyHists_[FEDid];
       //TH1F* iphiProfileHist = FEDsAndiPhiProfileHists_[FEDid];
       //TH1F* ietaProfileHist = FEDsAndiEtaProfileHists_[FEDid];
       //TH2F* timingHistVsFreq = FEDsAndTimingVsFreqHists_[FEDid];
       //TH2F* timingHistVsAmp = FEDsAndTimingVsAmpHists_[FEDid];
       //TH2F* E2vsE1uRecHist = FEDsAndE2vsE1Hists_[FEDid];
       //TH2F* energyvsE1uRecHist = FEDsAndenergyvsE1Hists_[FEDid];
       //TH1F* numXtalInClusterHist = FEDsAndNumXtalsInClusterHists_[FEDid];
       //TH2F* occupHist = FEDsAndOccupancyHists_[FEDid];
       //TH2F* timingHistVsPhi = FEDsAndTimingVsPhiHists_[FEDid];
       //TH2F* timingHistVsModule = FEDsAndTimingVsModuleHists_[FEDid];
       //if(uRecHist==0)
       //{
       //  initHists(FEDid);
       //  uRecHist = FEDsAndHists_[FEDid];
       //  E2uRecHist = FEDsAndE2Hists_[FEDid];
       //  energyuRecHist = FEDsAndenergyHists_[FEDid];
       //  timingHist = FEDsAndTimingHists_[FEDid];
       //  freqHist = FEDsAndFrequencyHists_[FEDid];
       //  timingHistVsFreq = FEDsAndTimingVsFreqHists_[FEDid];
       //  timingHistVsAmp = FEDsAndTimingVsAmpHists_[FEDid];
       //  iphiProfileHist = FEDsAndiPhiProfileHists_[FEDid];
       //  ietaProfileHist = FEDsAndiEtaProfileHists_[FEDid];
       //  E2vsE1uRecHist = FEDsAndE2vsE1Hists_[FEDid];
       //  energyvsE1uRecHist = FEDsAndenergyvsE1Hists_[FEDid];
       //  numXtalInClusterHist = FEDsAndNumXtalsInClusterHists_[FEDid];
       //  occupHist = FEDsAndOccupancyHists_[FEDid];
       //  timingHistVsPhi = FEDsAndTimingVsPhiHists_[FEDid];
       //  timingHistVsModule = FEDsAndTimingVsModuleHists_[FEDid];
       //}
       //uRecHist->Fill(ampli);
       //E2uRecHist->Fill(E2);
       //E2vsE1uRecHist->Fill(ampli,E2);
       //energyuRecHist->Fill(energy);
       //energyvsE1uRecHist->Fill(ampli,energy);
       //allFedsHist_->Fill(ampli);
 
       if (iz < 0) {
         EEM_FedsSeedEnergyHist_->Fill(ampli);
         EEM_FedsenergyHist_->Fill(energy);
         EEM_FedsenergyHighHist_->Fill(energy);
         EEM_FedsE2Hist_->Fill(E2);
         EEM_FedsE2vsE1Hist_->Fill(ampli, E2);
         EEM_FedsenergyvsE1Hist_->Fill(ampli, energy);
         EEM_AllOccupancyCoarse_->Fill(ix - 0.5, iy - 0.5);
         EEM_AllOccupancy_->Fill(ix - 0.5, iy - 0.5);
 
         EEM_FedsNumXtalsInClusterHist_->Fill(numXtalsinCluster);
         EEM_NumXtalsInClusterHist_->Fill(fullnumXtalsinCluster);
         EEM_numxtalsVsEnergy_->Fill(energy, numXtalsinCluster);
         EEM_numxtalsVsHighEnergy_->Fill(energy, numXtalsinCluster);
         EEM_numberofBCinSC_->Fill(clusSize);
 
         if (fullnumXtalsinCluster == 1) {
           EEM_OccupancySingleXtal_->Fill(ix - 0.5, iy - 0.5);
           EEM_energySingleXtalHist_->Fill(energy);
         }
 
         if (ampli > minTimingAmpEE_) {
           EEM_FedsTimingHist_->Fill(time);
           EEM_FedsTimingVsAmpHist_->Fill(time, ampli);
           EEM_FedsTimingTTHist_->Fill(ix - 0.5, iy - 0.5, time);
         }
 
         // Exclusive trigger plots
 
         if (isEcalL1 && !isDTL1 && !isRPCL1 && !isCSCL1 && !isHCALL1) {
           EEM_OccupancyExclusiveECAL_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseExclusiveECAL_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEM_FedsTimingHistECAL_->Fill(time);
             EEM_FedsTimingTTHistECAL_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEM_triggerExclusiveHist_->Fill(0);
         }
 
         if (!isEcalL1 && !isDTL1 && !isRPCL1 && !isCSCL1 && isHCALL1) {
           EEM_OccupancyExclusiveHCAL_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseExclusiveHCAL_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEM_FedsTimingHistHCAL_->Fill(time);
             EEM_FedsTimingTTHistHCAL_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEM_triggerExclusiveHist_->Fill(1);
         }
 
         if (!isEcalL1 && isDTL1 && !isRPCL1 && !isCSCL1 && !isHCALL1) {
           EEM_OccupancyExclusiveDT_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseExclusiveDT_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEM_FedsTimingHistDT_->Fill(time);
             EEM_FedsTimingTTHistDT_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEM_triggerExclusiveHist_->Fill(2);
         }
 
         if (!isEcalL1 && !isDTL1 && isRPCL1 && !isCSCL1 && !isHCALL1) {
           EEM_OccupancyExclusiveRPC_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseExclusiveRPC_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEM_FedsTimingHistRPC_->Fill(time);
             EEM_FedsTimingTTHistRPC_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEM_triggerExclusiveHist_->Fill(3);
         }
 
         if (!isEcalL1 && !isDTL1 && !isRPCL1 && isCSCL1 && !isHCALL1) {
           EEM_OccupancyExclusiveCSC_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseExclusiveCSC_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEM_FedsTimingHistCSC_->Fill(time);
             EEM_FedsTimingTTHistCSC_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEM_triggerExclusiveHist_->Fill(4);
         }
 
         // Inclusive trigger plots
 
         if (isEcalL1) {
           EEM_triggerHist_->Fill(0);
           EEM_OccupancyECAL_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseECAL_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isHCALL1) {
           EEM_triggerHist_->Fill(1);
           EEM_OccupancyHCAL_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseHCAL_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isDTL1) {
           EEM_triggerHist_->Fill(2);
           EEM_OccupancyDT_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseDT_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isRPCL1) {
           EEM_triggerHist_->Fill(3);
           EEM_OccupancyRPC_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseRPC_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isCSCL1) {
           EEM_triggerHist_->Fill(4);
           EEM_OccupancyCSC_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyCoarseCSC_->Fill(ix - 0.5, iy - 0.5);
         }
 
       } else {
         EEP_FedsSeedEnergyHist_->Fill(ampli);
         EEP_FedsenergyHist_->Fill(energy);
         EEP_FedsenergyHighHist_->Fill(energy);
         EEP_FedsE2Hist_->Fill(E2);
         EEP_FedsE2vsE1Hist_->Fill(ampli, E2);
         EEP_FedsenergyvsE1Hist_->Fill(ampli, energy);
         EEP_AllOccupancyCoarse_->Fill(ix - 0.5, iy - 0.5);
         EEP_AllOccupancy_->Fill(ix - 0.5, iy - 0.5);
 
         EEP_FedsNumXtalsInClusterHist_->Fill(numXtalsinCluster);
         EEP_NumXtalsInClusterHist_->Fill(fullnumXtalsinCluster);
         EEP_numxtalsVsEnergy_->Fill(energy, numXtalsinCluster);
         EEP_numxtalsVsHighEnergy_->Fill(energy, numXtalsinCluster);
         EEP_numberofBCinSC_->Fill(clusSize);
 
         if (fullnumXtalsinCluster == 1) {
           EEP_OccupancySingleXtal_->Fill(ix - 0.5, iy - 0.5);
           EEP_energySingleXtalHist_->Fill(energy);
         }
 
         if (ampli > minTimingAmpEE_) {
           EEP_FedsTimingHist_->Fill(time);
           EEP_FedsTimingVsAmpHist_->Fill(time, ampli);
           EEP_FedsTimingTTHist_->Fill(ix - 0.5, iy - 0.5, time);
         }
 
         // Exclusive trigger plots
 
         if (isEcalL1 && !isDTL1 && !isRPCL1 && !isCSCL1 && !isHCALL1) {
           EEP_OccupancyExclusiveECAL_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseExclusiveECAL_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEP_FedsTimingHistECAL_->Fill(time);
             EEP_FedsTimingTTHistECAL_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEP_triggerExclusiveHist_->Fill(0);
         }
 
         if (!isEcalL1 && !isDTL1 && !isRPCL1 && !isCSCL1 && isHCALL1) {
           EEP_OccupancyExclusiveHCAL_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseExclusiveHCAL_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEP_FedsTimingHistHCAL_->Fill(time);
             EEP_FedsTimingTTHistHCAL_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEP_triggerExclusiveHist_->Fill(1);
         }
 
         if (!isEcalL1 && isDTL1 && !isRPCL1 && !isCSCL1 && !isHCALL1) {
           EEP_OccupancyExclusiveDT_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseExclusiveDT_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEP_FedsTimingHistDT_->Fill(time);
             EEP_FedsTimingTTHistDT_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEP_triggerExclusiveHist_->Fill(2);
         }
 
         if (!isEcalL1 && !isDTL1 && isRPCL1 && !isCSCL1 && !isHCALL1) {
           EEP_OccupancyExclusiveRPC_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseExclusiveRPC_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEP_FedsTimingHistRPC_->Fill(time);
             EEP_FedsTimingTTHistRPC_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEP_triggerExclusiveHist_->Fill(3);
         }
 
         if (!isEcalL1 && !isDTL1 && !isRPCL1 && isCSCL1 && !isHCALL1) {
           EEP_OccupancyExclusiveCSC_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseExclusiveCSC_->Fill(ix - 0.5, iy - 0.5);
           if (ampli > minTimingAmpEE_) {
             EEP_FedsTimingHistCSC_->Fill(time);
             EEP_FedsTimingTTHistCSC_->Fill(ix - 0.5, iy - 0.5, time);
           }
           EEP_triggerExclusiveHist_->Fill(4);
         }
 
         // Inclusive trigger plots
 
         if (isEcalL1) {
           EEP_triggerHist_->Fill(0);
           EEP_OccupancyECAL_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseECAL_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isHCALL1) {
           EEP_triggerHist_->Fill(1);
           EEP_OccupancyHCAL_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseHCAL_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isDTL1) {
           EEP_triggerHist_->Fill(2);
           EEP_OccupancyDT_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseDT_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isRPCL1) {
           EEP_triggerHist_->Fill(3);
           EEP_OccupancyRPC_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseRPC_->Fill(ix - 0.5, iy - 0.5);
         }
         if (isCSCL1) {
           EEP_triggerHist_->Fill(4);
           EEP_OccupancyCSC_->Fill(ix - 0.5, iy - 0.5);
           EEP_OccupancyCoarseCSC_->Fill(ix - 0.5, iy - 0.5);
         }
       }
 
       // *** High Energy Clusters Analysis ** //
 
       if (energy > minHighEnergy_) {
         LogInfo("EcalCosmicsHists") << "High energy event in EE " << iEvent.id().run() << " : " << iEvent.id().event()
                                     << " " << naiveEvtNum_ << " : " << energy << " " << numXtalsinCluster << " : " << ix
                                     << " " << iy << " : " << isEcalL1 << isHCALL1 << isDTL1 << isRPCL1 << isCSCL1;
 
         //  numberOfHighEClusters++;
         //  allOccupancyHighEnergy_->Fill(iphi, ieta);
         //  allOccupancyHighEnergyCoarse_->Fill(iphi, ieta);
         //  allFedsOccupancyHighEnergyHist_->Fill(iphi,ieta,energy);
         if (iz < 0) {
           EEM_FedsenergyOnlyHighHist_->Fill(energy);
           EEM_OccupancyHighEnergy_->Fill(ix - 0.5, iy - 0.5);
           EEM_OccupancyHighEnergyCoarse_->Fill(ix - 0.5, iy - 0.5);
         } else {
           EEP_FedsenergyOnlyHighHist_->Fill(energy);
         }
         //  HighEnergy_2GeV_occuCoarse->Fill(iphi,ieta);
         //  HighEnergy_2GeV_occu3D->Fill(iphi,ieta,energy);
 
         //  HighEnergy_NumXtal->Fill(fullnumXtalsinCluster);
         //  HighEnergy_NumXtalFedId->Fill(FEDid,fullnumXtalsinCluster);
         //  HighEnergy_NumXtaliphi->Fill(iphi,fullnumXtalsinCluster);
         //  HighEnergy_energy3D->Fill(iphi,ieta,energy);
         //  HighEnergy_energyNumXtal->Fill(fullnumXtalsinCluster,energy);
 
         if (energy > 100.0) {
           LogInfo("EcalCosmicsHists") << "Very high energy event in EE " << iEvent.id().run() << " : "
                                       << iEvent.id().event() << " " << naiveEvtNum_ << " : " << energy << " "
                                       << numXtalsinCluster << " : " << ix << " " << iy << " : " << isEcalL1 << isHCALL1
                                       << isDTL1 << isRPCL1 << isCSCL1;
           //    HighEnergy_100GeV_occuCoarse->Fill(iphi,ieta);
           //    HighEnergy_100GeV_occu3D->Fill(iphi,ieta,energy);
         }
       }
 
       // *** end of High Energy Clusters analysis *** //
 
     }  //++++++++++++++++++END LOOP OVER EE SUPERCLUSTERS+++++++++++++++++++++++++++++++++++
   }
 
   HighEnergy_numClusHighEn->Fill(numberOfHighEClusters);
   HighEnergy_ratioClusters->Fill((double(numberOfHighEClusters)) / (double(numberOfCosmics)));
 
   numberofCosmicsHist_->Fill(numberOfCosmics);
   EEP_numberofCosmicsHist_->Fill(numberOfCosmicsEEP);
   EEM_numberofCosmicsHist_->Fill(numberOfCosmicsEEM);
   numberofCosmicsHistEB_->Fill(numberOfCosmicsEB);
 
   if (numberOfCosmics > 0) {
     cosmicCounter_++;
     numberofGoodEvtFreq_->Fill(naiveEvtNum_);
     allFedsFreqTimeHist_->Fill(eventTime);
     //This line will work in 21X!!
     //std::cout << " Orbit " <<   iEvent.orbitNumber() << " BX " << iEvent.bunchCrossing()<< std::endl;
     //std::cout << " BX " << iEvent.experimentType() << std::endl;
   }
 
   if (numberOfCosmicsEB > 0)
     cosmicCounterEB_++;
   if (numberOfCosmicsEEP > 0)
     cosmicCounterEEP_++;
   if (numberOfCosmicsEEM > 0)
     cosmicCounterEEM_++;
 
   if (numberOfCosmicsTop && numberOfCosmicsBottom) {
     cosmicCounterTopBottom_++;
     numberofCosmicsTopBottomHist_->Fill(numberOfCosmicsTop + numberOfCosmicsBottom);
   }
 
   // *** TrackAssociator *** //
 
   // get reco tracks
   edm::Handle<reco::TrackCollection> recoTracks;
   iEvent.getByLabel("cosmicMuons", recoTracks);
 
   if (recoTracks.isValid()) {
     //    LogWarning("EcalCosmicsHists") << "... Valid TrackAssociator recoTracks !!! " << recoTracks.product()->size();
     std::map<int, std::vector<DetId> > trackDetIdMap;
     int tracks = 0;
     for (reco::TrackCollection::const_iterator recoTrack = recoTracks->begin(); recoTrack != recoTracks->end();
          ++recoTrack) {
       if (fabs(recoTrack->d0()) > 70 || fabs(recoTrack->dz()) > 70)
         continue;
       if (recoTrack->numberOfValidHits() < 20)
         continue;
 
       //if (recoTrack->pt() < 2) continue; // skip low Pt tracks
 
       TrackDetMatchInfo info = trackAssociator_.associate(iEvent, iSetup, *recoTrack, trackParameters_);
 
       //       edm::LogVerbatim("TrackAssociator") << "\n-------------------------------------------------------\n Track (pt,eta,phi): " <<
       //    recoTrack->pt() << " , " << recoTrack->eta() << " , " << recoTrack->phi() ;
       //       edm::LogVerbatim("TrackAssociator") << "Ecal energy in crossed crystals based on RecHits: " <<
       //    info.crossedEnergy(TrackDetMatchInfo::EcalRecHits);
       //       edm::LogVerbatim("TrackAssociator") << "Ecal energy in 3x3 crystals based on RecHits: " <<
       //    info.nXnEnergy(TrackDetMatchInfo::EcalRecHits, 1);
       //       edm::LogVerbatim("TrackAssociator") << "Hcal energy in crossed towers based on RecHits: " <<
       //    info.crossedEnergy(TrackDetMatchInfo::HcalRecHits);
       //       edm::LogVerbatim("TrackAssociator") << "Hcal energy in 3x3 towers based on RecHits: " <<
       //    info.nXnEnergy(TrackDetMatchInfo::HcalRecHits, 1);
       //       edm::LogVerbatim("TrackAssociator") << "Number of muon segment matches: " << info.numberOfSegments();
 
       //       std::cout << "\n-------------------------------------------------------\n Track (pt,eta,phi): " <<
       //    recoTrack->pt() << " , " << recoTrack->eta() << " , " << recoTrack->phi() << std::endl;
       //       std::cout << "Ecal energy in crossed crystals based on RecHits: " <<
       //    info.crossedEnergy(TrackDetMatchInfo::EcalRecHits) << std::endl;
       //       std::cout << "Ecal energy in 3x3 crystals based on RecHits: " <<
       //    info.nXnEnergy(TrackDetMatchInfo::EcalRecHits, 1) << std::endl;
       //       std::cout << "Hcal energy in crossed towers based on RecHits: " <<
       //    info.crossedEnergy(TrackDetMatchInfo::HcalRecHits) << std::endl;
       //       std::cout << "Hcal energy in 3x3 towers based on RecHits: " <<
       //    info.nXnEnergy(TrackDetMatchInfo::HcalRecHits, 1) << std::endl;
 
       for (unsigned int i = 0; i < info.crossedEcalIds.size(); i++) {
         // only checks for barrel
         if (info.crossedEcalIds[i].det() == DetId::Ecal && info.crossedEcalIds[i].subdetId() == 1) {
           EBDetId ebDetId(info.crossedEcalIds[i]);
           trackAssoc_muonsEcal_->Fill(ebDetId.iphi(), ebDetId.ieta());
           //std::cout << "Crossed iphi: " << ebDetId.iphi()
           //    << " ieta: " << ebDetId.ieta() << " : nCross " << info.crossedEcalIds.size() << std::endl;
 
           EcalRecHitCollection::const_iterator thishit = hits->find(ebDetId);
           if (thishit == hits->end())
             continue;
 
           //EcalRecHit myhit = (*thishit);
           //double thisamp = myhit.energy();
           //std::cout << " Crossed energy: " << thisamp << " : nCross " << info.crossedEcalIds.size() << std::endl;
         }
       }
 
       //edm::LogVerbatim("TrackAssociator") << " crossedEcalIds size: " << info.crossedEcalIds.size()
       //                  << " crossedEcalRecHits size: " << info.crossedEcalRecHits.size();
       numberofCrossedEcalIdsHist_->Fill(info.crossedEcalIds.size());
       tracks++;
       if (!info.crossedEcalIds.empty())
         trackDetIdMap.insert(std::pair<int, std::vector<DetId> >(tracks, info.crossedEcalIds));
     }
 
     // Now to match recoTracks with cosmic clusters
 
     int numAssocTracks = 0;
     int numAssocClusters = 0;
     edm::LogVerbatim("TrackAssociator") << "Matching cosmic clusters to tracks...";
     int numSeeds = seeds.size();
     int numTracks = trackDetIdMap.size();
     while (!seeds.empty() && !trackDetIdMap.empty()) {
       double bestDr = 1000;
       double bestDPhi = 1000;
       double bestDEta = 1000;
       double bestEtaTrack = 1000;
       double bestEtaSeed = 1000;
       double bestPhiTrack = 1000;
       double bestPhiSeed = 1000;
       EBDetId bestTrackDet;
       EBDetId bestSeed;
       int bestTrack = -1;
       std::map<EBDetId, EBDetId> trackDetIdToSeedMap;
 
       //edm::LogVerbatim("TrackAssociator") << "NumTracks:" << trackDetIdMap.size() << " numClusters:" << seeds.size();
 
       for (std::vector<EBDetId>::const_iterator seedItr = seeds.begin(); seedItr != seeds.end(); ++seedItr) {
         for (std::map<int, std::vector<DetId> >::const_iterator mapItr = trackDetIdMap.begin();
              mapItr != trackDetIdMap.end();
              ++mapItr) {
           for (unsigned int i = 0; i < mapItr->second.size(); i++) {
             // only checks for barrel
             if (mapItr->second[i].det() == DetId::Ecal && mapItr->second[i].subdetId() == 1) {
               EBDetId ebDet = (mapItr->second[i]);
               double seedEta = seedItr->ieta();
               double deta = ebDet.ieta() - seedEta;
               if (seedEta * ebDet.ieta() < 0)
                 deta > 0 ? (deta = deta - 1.) : (deta = deta + 1.);
               double dR;
               double dphi = ebDet.iphi() - seedItr->iphi();
               if (abs(dphi) > 180)
                 dphi > 0 ? (dphi = 360 - dphi) : (dphi = -360 - dphi);
               dR = sqrt(deta * deta + dphi * dphi);
               if (dR < bestDr) {
                 bestDr = dR;
                 bestDPhi = dphi;
                 bestDEta = deta;
                 bestTrackDet = mapItr->second[i];
                 bestTrack = mapItr->first;
                 bestSeed = (*seedItr);
                 bestEtaTrack = ebDet.ieta();
                 bestEtaSeed = seedEta;
                 bestPhiTrack = ebDet.iphi();
                 bestPhiSeed = seedItr->iphi();
               }
             }
           }
         }
       }
       if (bestDr < 1000) {
         //edm::LogVerbatim("TrackAssociator") << "Best deltaR from matched DetId's to cluster:" << bestDr;
         deltaRHist_->Fill(bestDr);
         deltaPhiHist_->Fill(bestDPhi);
         deltaEtaHist_->Fill(bestDEta);
         deltaEtaDeltaPhiHist_->Fill(bestDEta, bestDPhi);
         seedTrackEtaHist_->Fill(bestEtaSeed, bestEtaTrack);
         seedTrackPhiHist_->Fill(bestPhiSeed, bestPhiTrack);
         seeds.erase(find(seeds.begin(), seeds.end(), bestSeed));
         trackDetIdMap.erase(trackDetIdMap.find(bestTrack));
         trackDetIdToSeedMap[bestTrackDet] = bestSeed;
         numAssocTracks++;
         numAssocClusters++;
 
         // for high energy analysis
         if (bestDPhi < 1.5 && bestDEta < 1.8) {
           //edm::LogVerbatim("TrackAssociator") << "Best seed ieta and iphi: ("
           //                    << bestSeed.ieta() << ", " << bestSeed.iphi() << ") ";
           //check if the bestSeed is a High Energy one
           EcalRecHitCollection::const_iterator Ecalhit = hits->find(bestSeed);
           if (Ecalhit == hits->end()) {
             continue;
           }
           double EcalhitEnergy = Ecalhit->energy();
           //edm::LogVerbatim("TrackAssociator") << "Best seed energy: " << EcalhitEnergy ;
           HighEnergy_bestSeed->Fill(EcalhitEnergy);
           HighEnergy_bestSeedOccupancy->Fill(bestSeed.iphi(), bestSeed.ieta());
         }
       } else {
         edm::LogVerbatim("TrackAssociator") << "could not match cluster seed to track." << bestDr;
         break;  // no match found
       }
     }
     if (numSeeds > 0 && numTracks > 0) {
       ratioAssocClustersHist_->AddBinContent(1, numAssocClusters);
       ratioAssocClustersHist_->AddBinContent(2, numSeeds);
     }
     if (numTracks > 0) {
       ratioAssocTracksHist_->AddBinContent(1, numAssocTracks);
       ratioAssocTracksHist_->AddBinContent(2, numTracks);
       numberofCosmicsWTrackHist_->Fill(numSeeds);
     }
   } else {
     //    LogWarning("EcalCosmicsHists") << "!!! No TrackAssociator recoTracks !!!";
   }
 
   // Study on Tracks for High Energy events
   edm::Handle<reco::TrackCollection> recoTracksBarrel;
   iEvent.getByLabel("cosmicMuonsBarrelOnly", recoTracksBarrel);
   if (!recoTracksBarrel.isValid()) {
     //edm::LogWarning("EcalCosmicsHists") << "RecoTracksBarrel not valid!! " ;
   } else {
     //edm::LogWarning("EcalCosmicsHists") << "Number of barrel reco tracks found in the event: " << recoTracksBarrel->size() ;
     HighEnergy_numRecoTrackBarrel->Fill(recoTracksBarrel->size());
     for (reco::TrackCollection::const_iterator recoTrack = recoTracksBarrel->begin();
          recoTrack != recoTracksBarrel->end();
          ++recoTrack) {
       //edm::LogWarning("EcalCosmicsHists") << "Track (pt,eta,phi): " << recoTrack->pt() << " , " << recoTrack->eta() << " , " << recoTrack->phi() ;
       //edm::LogWarning("EcalCosmicsHists") << "Track innermost hit: " << recoTrack->innerPosition().phi() ;
     }
 
     for (reco::SuperClusterCollection::const_iterator clus = clusterCollection_p->begin();
          clus != clusterCollection_p->end();
          ++clus) {
       double energy = clus->energy();
       double phi = clus->phi();
       double eta = clus->eta();
 
       if (recoTracksBarrel->empty())
         HighEnergy_0tracks_occu3D->Fill(phi, eta, energy);
       if (recoTracksBarrel->size() == 1)
         HighEnergy_1tracks_occu3D->Fill(phi, eta, energy);
       if (recoTracksBarrel->size() == 2)
         HighEnergy_2tracks_occu3D->Fill(phi, eta, energy);
 
       std::vector<std::pair<DetId, float> > clusterDetIds = clus->hitsAndFractions();  //get these from the cluster
       for (std::vector<std::pair<DetId, float> >::const_iterator detitr = clusterDetIds.begin();
            detitr != clusterDetIds.end();
            ++detitr) {
         if ((*detitr).first.det() != DetId::Ecal) {
           continue;
         }
         if ((*detitr).first.subdetId() != EcalBarrel) {
           continue;
         }
         EcalRecHitCollection::const_iterator thishit = hits->find((*detitr).first);
         if (thishit == hits->end()) {
           continue;
         }
 
         double rechitenergy = thishit->energy();
         int ieta = ((EBDetId)(*detitr).first).ieta();
         int iphi = ((EBDetId)(*detitr).first).iphi();
         if (rechitenergy > minRecHitAmpEB_) {
           if (recoTracksBarrel->empty())
             HighEnergy_0tracks_occu3DXtal->Fill(iphi, ieta, rechitenergy);
           if (recoTracksBarrel->size() == 1)
             HighEnergy_1tracks_occu3DXtal->Fill(iphi, ieta, rechitenergy);
           if (recoTracksBarrel->size() == 2)
             HighEnergy_2tracks_occu3DXtal->Fill(iphi, ieta, rechitenergy);
 
           if (rechitenergy > 10)
             edm::LogWarning("EcalCosmicsHists") << "!!!!! Crystal with energy " << rechitenergy << " at (ieta,iphi) ("
                                                 << ieta << " ," << iphi << "); Id: " << (*detitr).first.det();
         }
       }
     }
 
     // look at angle between 2 recoTracks
     if (recoTracksBarrel->size() == 2) {
       reco::TrackCollection::const_iterator Track1 = recoTracksBarrel->begin();
       reco::TrackCollection::const_iterator Track2 = (recoTracksBarrel->begin()) + 1;
       float angle = (acos(sin(Track1->theta()) * cos(Track1->phi()) * sin(Track2->theta()) * cos(Track2->phi()) +
                           sin(Track1->theta()) * sin(Track1->phi()) * sin(Track2->theta()) * sin(Track2->phi()) +
                           cos(Track1->theta()) * cos(Track2->theta()))) /
                     0.017453292519943;
       //edm::LogWarning("EcalCosmicsHists") << "Tracks angle: " << angle;
       HighEnergy_TracksAngle->Fill(angle);
       if ((Track1->innerPosition().phi()) > 0 && (Track2->innerPosition().phi()) < 0) {
         //edm::LogWarning("EcalCosmicsHists") << "Top-bottom tracks";
         HighEnergy_TracksAngleTopBottom->Fill(angle);
       }
     }
   }
 
   // *** end of TrackAssociator code *** //
 
   // *** HCAL RecHits code *** //
 
   //hcal rechits
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByLabel("hbhereco", hbhe);
 
   edm::Handle<HFRecHitCollection> hfrh;
   iEvent.getByLabel("hfreco", hfrh);
 
   edm::Handle<HORecHitCollection> horh;
   iEvent.getByLabel("horeco", horh);
 
   if (hbhe.isValid()) {
     //    LogInfo("EcalCosmicHists") << "event " << ievt << " HBHE RecHits collection size " << hbhe->size();
     const HBHERecHitCollection hbheHit = *(hbhe.product());
     for (HBHERecHitCollection::const_iterator hhit = hbheHit.begin(); hhit != hbheHit.end(); hhit++) {
       //      if (hhit->energy() > 0.6){
       hcalEnergy_HBHE_->Fill(hhit->energy());
       //      }
     }
   } else {
     //    LogWarning("EcalCosmicHists") << " HBHE RecHits **NOT** VALID!! " << endl;
   }
 
   if (hfrh.isValid()) {
     //    LogInfo("EcalCosmicHists") << "event " << ievt << " HF RecHits collection size " << hfrh->size();
     const HFRecHitCollection hfHit = *(hfrh.product());
     for (HFRecHitCollection::const_iterator hhit = hfHit.begin(); hhit != hfHit.end(); hhit++) {
       hcalEnergy_HF_->Fill(hhit->energy());
     }
   } else {
     //    LogWarning("EcalCosmicHists") << " HF RecHits **NOT** VALID!! " << endl;
   }
 
   if (horh.isValid()) {
     //    LogInfo("EcalCosmicHists") << "event " << ievt << " HO RecHits collection size " << horh->size();
     const HORecHitCollection hoHit = *(horh.product());
     for (HORecHitCollection::const_iterator hhit = hoHit.begin(); hhit != hoHit.end(); hhit++) {
       //     if (hhit->energy() > 0.6){
       hcalEnergy_HO_->Fill(hhit->energy());
       //     }
     }
   } else {
     //    LogWarning("EcalCosmicHists") << " HO RecHits **NOT** VALID!! " << endl;
   }
 
   // *** end of HCAL code *** //
 }

void EcalCosmicsHists::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	eventSetup
	)

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1569 of file EcalCosmicsHists.cc.

                                                                               {
   edm::ESHandle<EcalElectronicsMapping> handle;
   eventSetup.get<EcalMappingRcd>().get(handle);
   ecalElectronicsMap_ = handle.product();
 
   //Here I will init some of the specific histograms
   int numBins = 200;  //(int)round(histRangeMax_-histRangeMin_)+1;
 
   //=============Special Bins for TT and Modules borders=============================
   double ttEtaBins[36] = {-85, -80, -75, -70, -65, -60, -55, -50, -45, -40, -35, -30, -25, -20, -15, -10, -5, 0,
                           1,   6,   11,  16,  21,  26,  31,  36,  41,  46,  51,  56,  61,  66,  71,  76,  81, 86};
   double modEtaBins[10] = {-85, -65, -45, -25, 0, 1, 26, 46, 66, 86};
   double ttPhiBins[73];
   double modPhiBins[19];
   double timingBins[79];
   double highEBins[11];
   for (int i = 0; i < 79; ++i) {
     timingBins[i] = -7. + double(i) * 14. / 78.;
     if (i < 73) {
       ttPhiBins[i] = 1 + 5 * i;
       if (i < 19) {
         modPhiBins[i] = 1 + 20 * i;
         if (i < 11) {
           highEBins[i] = 10. + double(i) * 20.;
         }
       }
     }
   }
   //=============END Special Bins for TT and Modules borders===========================
 
   //====================Frequency for Event timing information=========================
   double timingEndInSeconds = runTimeLength_ * 3600.;
   double timingBinWidth = timingEndInSeconds / double(numTimingBins_);
   //====================END Frequency for Event timing information=====================
 
   allFedsenergyHist_ =
       new TH1F("energy_AllClusters", "energy_AllClusters;Cluster Energy (GeV)", numBins, histRangeMin_, histRangeMax_);
   allFedsenergyHighHist_ =
       new TH1F("energyHigh_AllClusters", "energyHigh_AllClusters;Cluster Energy (GeV)", numBins, histRangeMin_, 200.0);
   allFedsenergyOnlyHighHist_ = new TH1F("energyHigh_HighEnergyClusters",
                                         "energy of High Energy Clusters;Cluster Energy (GeV)",
                                         numBins,
                                         histRangeMin_,
                                         200.0);
   allFedsE2Hist_ = new TH1F(
       "E2_AllClusters", "E2_AllClusters;Seed+highest neighbor energy (GeV)", numBins, histRangeMin_, histRangeMax_);
   allFedsE2vsE1Hist_ = new TH2F("E2vsE1_AllClusters",
                                 "E2vsE1_AllClusters;Seed Energy (GeV);Seed+highest neighbor energy (GeV)",
                                 numBins,
                                 histRangeMin_,
                                 histRangeMax_,
                                 numBins,
                                 histRangeMin_,
                                 histRangeMax_);
   allFedsenergyvsE1Hist_ = new TH2F("energyvsE1_AllClusters",
                                     "energyvsE1_AllClusters;Seed Energy (GeV);Energy(GeV)",
                                     numBins,
                                     histRangeMin_,
                                     histRangeMax_,
                                     numBins,
                                     histRangeMin_,
                                     histRangeMax_);
   allFedsTimingHist_ = new TH1F("timeForAllFeds", "timeForAllFeds;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingVsFreqHist_ = new TH2F("timeVsFreqAllEvent",
                                       "time Vs Freq All events;Relative Time (1 clock = 25ns);Event Number",
                                       78,
                                       -7,
                                       7,
                                       2000,
                                       0.,
                                       200000);
   allFedsTimingVsAmpHist_ = new TH2F("timeVsAmpAllEvents",
                                      "time Vs Amp All Events;Relative Time (1 clock = 25ns);Amplitude (GeV)",
                                      78,
                                      -7,
                                      7,
                                      numBins,
                                      histRangeMin_,
                                      histRangeMax_);
   allFedsFrequencyHist_ = new TH1F("FrequencyAllEvent", "Frequency for All events;Event Number", 2000, 0., 200000);
 
   //--------Special Hists for times stamp info----------------------------------------------
   allFedsFreqTimeHist_ = new TH1F("FrequencyAllEventsInTime",
                                   Form("Time of Cosmic Events; Time (s);Passing Event rate/%5g s", timingBinWidth),
                                   numTimingBins_,
                                   0.,
                                   timingEndInSeconds);
   allFedsFreqTimeVsPhiHist_ =
       new TH2F("FrequencyAllEventsInTimeVsPhi",
                Form("Time of Cosmic Events vs iPhi; iPhi;Time (s)/%5g s", timingBinWidth * 360.),
                360,
                1.,
                361.,
                numTimingBins_ / 360,
                0.,
                timingEndInSeconds);
   allFedsFreqTimeVsPhiTTHist_ =
       new TH2F("FrequencyAllEventsInTimeVsTTPhi",
                Form("Time of Cosmic Events vs iPhi (TT bins); iPhi;Time (s)/%5g s", timingBinWidth * 72.),
                72,
                1.,
                361.,
                numTimingBins_ / 72,
                0.,
                timingEndInSeconds);
   allFedsFreqTimeVsEtaHist_ =
       new TH2F("FrequencyAllEventsInTimeVsEta",
                Form("Time of Cosmic Events vs iEta; Time (s)/%5g s; iEta", timingBinWidth * 172.),
                numTimingBins_ / 172,
                0.,
                timingEndInSeconds,
                172,
                -86.,
                86.);
   allFedsFreqTimeVsEtaTTHist_ =
       new TH2F("FrequencyAllEventsInTimeVsTTEta",
                Form("Time of Cosmic Events vs Eta (TT bins);Time (s)/%5g s; iEta", timingBinWidth * 35.),
                numTimingBins_ / 35,
                0.,
                timingEndInSeconds,
                35,
                ttEtaBins);
   //--------END Special Hists for times stamp info------------------------------------------
 
   allFedsiPhiProfileHist_ = new TH1F("iPhiProfileAllEvents", "iPhi Profile all events;i#phi", 360, 1., 361.);
   allFedsiEtaProfileHist_ = new TH1F("iEtaProfileAllEvents", "iEta Profile all events;i#eta", 172, -86, 86);
 
   allOccupancy_ = new TH2F("OccupancyAllEvents", "Occupancy all events;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   TrueOccupancy_ =
       new TH2F("TrueOccupancyAllEvents", "True Occupancy all events;#phi;#eta", 360, -3.14159, 3.14159, 172, -1.5, 1.5);
   allOccupancyCoarse_ =
       new TH2F("OccupancyAllEventsCoarse", "Occupancy all events Coarse;i#phi;i#eta", 360 / 5, 1, 361., 35, ttEtaBins);
   TrueOccupancyCoarse_ = new TH2F("TrueOccupancyAllEventsCoarse",
                                   "True Occupancy all events Coarse;#phi;#eta",
                                   360 / 5,
                                   -3.14159,
                                   3.14159,
                                   34,
                                   -1.5,
                                   1.5);
 
   // single xtal cluster occupancy
   allOccupancySingleXtal_ =
       new TH2F("OccupancySingleXtal", "Occupancy single xtal clusters;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   energySingleXtalHist_ = new TH1F(
       "energy_SingleXtalClusters", "Energy single xtal clusters;Cluster Energy (GeV)", numBins, histRangeMin_, 200.0);
 
   allFedsTimingPhiHist_ = new TH2F("timePhiAllFEDs",
                                    "time vs Phi for all FEDs (TT binning);i#phi;Relative Time (1 clock = 25ns)",
                                    72,
                                    1,
                                    361,
                                    78,
                                    -7,
                                    7);
   allFedsTimingPhiEbpHist_ = new TH2F("timePhiEBP",
                                       "time vs Phi for FEDs in EB+ (TT binning) ;i#phi;Relative Time (1 clock = 25ns)",
                                       72,
                                       1,
                                       361,
                                       78,
                                       -7,
                                       7);
   allFedsTimingPhiEbmHist_ = new TH2F("timePhiEBM",
                                       "time vs Phi for FEDs in EB- (TT binning);i#phi;Relative Time (1 clock = 25ns)",
                                       72,
                                       1,
                                       361,
                                       78,
                                       -7,
                                       7);
   allFedsTimingPhiEtaHist_ =
       new TH3F("timePhiEtaAllFEDs",
                "(Phi,Eta,time) for all FEDs (SM,M binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                18,
                modPhiBins,
                9,
                modEtaBins,
                78,
                timingBins);
   allFedsTimingTTHist_ =
       new TH3F("timeTTAllFEDs",
                "(Phi,Eta,time) for all FEDs (SM,TT binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                72,
                ttPhiBins,
                35,
                ttEtaBins,
                78,
                timingBins);
   allFedsTimingLMHist_ = new TH2F(
       "timeLMAllFEDs", "(LM,time) for all FEDs (SM,LM binning);LM;Relative Time (1 clock = 25ns)", 92, 1, 92, 78, -7, 7);
 
   allFedsTimingEbpHist_ = new TH1F("timeEBP", "time for FEDs in EB+;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingEbmHist_ = new TH1F("timeEBM", "time for FEDs in EB-;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingEbpTopHist_ =
       new TH1F("timeEBPTop", "time for FEDs in EB+ Top;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingEbmTopHist_ =
       new TH1F("timeEBMTop", "time for FEDs in EB- Top;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingEbpBottomHist_ =
       new TH1F("timeEBPBottom", "time for FEDs in EB+ Bottom;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingEbmBottomHist_ =
       new TH1F("timeEBMBottom", "time for FEDs in EB- Bottom;Relative Time (1 clock = 25ns)", 78, -7, 7);
 
   numberofCosmicsHist_ =
       new TH1F("numberofCosmicsPerEvent", "Number of cosmics per event;Number of Cosmics", 30, 0, 30);
   numberofCosmicsHistEB_ =
       new TH1F("numberofCosmicsPerEvent_EB", "Number of cosmics per event EB;Number of Cosmics", 30, 0, 30);
 
   numberofCosmicsWTrackHist_ =
       new TH1F("numberofCosmicsWTrackPerEvent", "Number of cosmics with track per event", 30, 0, 30);
   numberofCosmicsTopBottomHist_ = new TH1F(
       "numberofCosmicsTopBottomPerEvent", "Number of top bottom cosmics per event;Number of Cosmics", 30, 0, 30);
   numberofGoodEvtFreq_ = new TH1F("frequencyOfGoodEvents",
                                   "Number of events with cosmic vs Event;Event Number;Number of Good Events/100 Events",
                                   2000,
                                   0,
                                   200000);
 
   numberofCrossedEcalIdsHist_ = new TH1F("numberofCrossedEcalCosmicsPerEvent",
                                          "Number of crossed ECAL cosmics per event;Number of Crossed Cosmics",
                                          10,
                                          0,
                                          10);
 
   allOccupancyExclusiveECAL_ = new TH2F("OccupancyAllEvents_ExclusiveECAL",
                                         "Occupancy all events Exclusive ECAL ;i#phi;i#eta",
                                         360,
                                         1.,
                                         361.,
                                         172,
                                         -86,
                                         86);
   allOccupancyCoarseExclusiveECAL_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveECAL",
                                               "Occupancy all events Coarse Exclusive ECAL;i#phi;i#eta",
                                               360 / 5,
                                               ttPhiBins,
                                               35,
                                               ttEtaBins);
   allOccupancyECAL_ =
       new TH2F("OccupancyAllEvents_ECAL", "Occupancy all events ECAL;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseECAL_ = new TH2F("OccupancyAllEventsCoarse_ECAL",
                                      "Occupancy all events Coarse ECAL;i#phi;i#eta",
                                      360 / 5,
                                      ttPhiBins,
                                      35,
                                      ttEtaBins);
   allFedsTimingHistECAL_ =
       new TH1F("timeForAllFeds_ECAL", "timeForAllFeds ECAL;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingPhiEtaHistECAL_ =
       new TH3F("timePhiEtaAllFEDs_ECAL",
                "ECAL (Phi,Eta,time) for all FEDs (SM,M binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                18,
                modPhiBins,
                9,
                modEtaBins,
                78,
                timingBins);
   allFedsTimingTTHistECAL_ =
       new TH3F("timeTTAllFEDs_ECAL",
                "ECAL (Phi,Eta,time) for all FEDs (SM,TT binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                72,
                ttPhiBins,
                35,
                ttEtaBins,
                78,
                timingBins);
   allFedsTimingLMHistECAL_ = new TH2F("timeLMAllFEDs_ECAL",
                                       "ECAL (LM,time) for all FEDs (SM,LM binning);LM;Relative Time (1 clock = 25ns)",
                                       92,
                                       1,
                                       92,
                                       78,
                                       -7,
                                       7);
 
   allOccupancyExclusiveDT_ = new TH2F(
       "OccupancyAllEvents_ExclusiveDT", "Occupancy all events Exclusive DT;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseExclusiveDT_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveDT",
                                             "Occupancy all events Coarse Exclusive DT;i#phi;i#eta",
                                             360 / 5,
                                             1,
                                             361.,
                                             35,
                                             ttEtaBins);
   allOccupancyDT_ =
       new TH2F("OccupancyAllEvents_DT", "Occupancy all events DT;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseDT_ = new TH2F(
       "OccupancyAllEventsCoarse_DT", "Occupancy all events Coarse DT;i#phi;i#eta", 360 / 5, 1, 361., 35, ttEtaBins);
   allFedsTimingHistDT_ = new TH1F("timeForAllFeds_DT", "timeForAllFeds DT;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingPhiEtaHistDT_ =
       new TH3F("timePhiEtaAllFEDs_DT",
                "DT (Phi,Eta,time) for all FEDs (SM,M binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                18,
                modPhiBins,
                9,
                modEtaBins,
                78,
                timingBins);
   allFedsTimingTTHistDT_ =
       new TH3F("timeTTAllFEDs_DT",
                "DT (Phi,Eta,time) for all FEDs (SM,TT binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                72,
                ttPhiBins,
                35,
                ttEtaBins,
                78,
                timingBins);
   allFedsTimingLMHistDT_ = new TH2F("timeLMAllFEDs_DT",
                                     "DT (LM,time) for all FEDs (SM,LM binning);LM;Relative Time (1 clock = 25ns)",
                                     92,
                                     1,
                                     92,
                                     78,
                                     -7,
                                     7);
 
   allOccupancyExclusiveRPC_ = new TH2F(
       "OccupancyAllEvents_ExclusiveRPC", "Occupancy all events Exclusive RPC;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseExclusiveRPC_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveRPC",
                                              "Occupancy all events Coarse Exclusive RPC;i#phi;i#eta",
                                              360 / 5,
                                              1,
                                              361.,
                                              35,
                                              ttEtaBins);
   allOccupancyRPC_ =
       new TH2F("OccupancyAllEvents_RPC", "Occupancy all events RPC;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseRPC_ = new TH2F(
       "OccupancyAllEventsCoarse_RPC", "Occupancy all events Coarse RPC;i#phi;i#eta", 360 / 5, 1, 361., 35, ttEtaBins);
   allFedsTimingHistRPC_ =
       new TH1F("timeForAllFeds_RPC", "timeForAllFeds RPC;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingPhiEtaHistRPC_ =
       new TH3F("timePhiEtaAllFEDs_RPC",
                "RPC (Phi,Eta,time) for all FEDs (SM,M binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                18,
                modPhiBins,
                9,
                modEtaBins,
                78,
                timingBins);
   allFedsTimingTTHistRPC_ =
       new TH3F("timeTTAllFEDs_RPC",
                "RPC (Phi,Eta,time) for all FEDs (SM,TT binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                72,
                ttPhiBins,
                35,
                ttEtaBins,
                78,
                timingBins);
   allFedsTimingLMHistRPC_ = new TH2F("timeLMAllFEDs_RPC",
                                      "RPC (LM,time) for all FEDs (SM,LM binning);LM;Relative Time (1 clock = 25ns)",
                                      92,
                                      1,
                                      92,
                                      78,
                                      -7,
                                      7);
 
   allOccupancyExclusiveCSC_ = new TH2F(
       "OccupancyAllEvents_ExclusiveCSC", "Occupancy all events Exclusive CSC;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseExclusiveCSC_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveCSC",
                                              "Occupancy all events Coarse Exclusive CSC;i#phi;i#eta",
                                              360 / 5,
                                              1,
                                              361.,
                                              35,
                                              ttEtaBins);
   allOccupancyCSC_ =
       new TH2F("OccupancyAllEvents_CSC", "Occupancy all events CSC;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseCSC_ = new TH2F(
       "OccupancyAllEventsCoarse_CSC", "Occupancy all events Coarse CSC;i#phi;i#eta", 360 / 5, 1, 361., 35, ttEtaBins);
   allFedsTimingHistCSC_ =
       new TH1F("timeForAllFeds_CSC", "timeForAllFeds CSC;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingPhiEtaHistCSC_ =
       new TH3F("timePhiEtaAllFEDs_CSC",
                "CSC (Phi,Eta,time) for all FEDs (SM,M binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                18,
                modPhiBins,
                9,
                modEtaBins,
                78,
                timingBins);
   allFedsTimingTTHistCSC_ =
       new TH3F("timeTTAllFEDs_CSC",
                "CSC (Phi,Eta,time) for all FEDs (SM,TT binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                72,
                ttPhiBins,
                35,
                ttEtaBins,
                78,
                timingBins);
   allFedsTimingLMHistCSC_ = new TH2F("timeLMAllFEDs_CSC",
                                      "CSC (LM,time) for all FEDs (SM,LM binning);LM;Relative Time (1 clock = 25ns)",
                                      92,
                                      1,
                                      62,
                                      78,
                                      -7,
                                      7);
 
   allOccupancyExclusiveHCAL_ = new TH2F("OccupancyAllEvents_ExclusiveHCAL",
                                         "Occupancy all events Exclusive HCAL;i#phi;i#eta",
                                         360,
                                         1.,
                                         361.,
                                         172,
                                         -86,
                                         86);
   allOccupancyCoarseExclusiveHCAL_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveHCAL",
                                               "Occupancy all events Coarse Exclusive HCAL;i#phi;i#eta",
                                               360 / 5,
                                               1,
                                               361.,
                                               35,
                                               ttEtaBins);
   allOccupancyHCAL_ =
       new TH2F("OccupancyAllEvents_HCAL", "Occupancy all events HCAL;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyCoarseHCAL_ = new TH2F(
       "OccupancyAllEventsCoarse_HCAL", "Occupancy all events Coarse HCAL;i#phi;i#eta", 360 / 5, 1, 361., 35, ttEtaBins);
   allFedsTimingHistHCAL_ =
       new TH1F("timeForAllFeds_HCAL", "timeForAllFeds HCAL;Relative Time (1 clock = 25ns)", 78, -7, 7);
   allFedsTimingPhiEtaHistHCAL_ =
       new TH3F("timePhiEtaAllFEDs_HCAL",
                "HCAL (Phi,Eta,time) for all FEDs (SM,M binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                18,
                modPhiBins,
                9,
                modEtaBins,
                78,
                timingBins);
   allFedsTimingTTHistHCAL_ =
       new TH3F("timeTTAllFEDs_HCAL",
                "HCAL (Phi,Eta,time) for all FEDs (SM,TT binning);i#phi;i#eta;Relative Time (1 clock = 25ns)",
                72,
                ttPhiBins,
                35,
                ttEtaBins,
                78,
                timingBins);
   allFedsTimingLMHistHCAL_ = new TH2F("timeLMAllFEDs_HCAL",
                                       "HCAL (LM,time) for all FEDs (SM,LM binning);LM;Relative Time (1 clock = 25ns)",
                                       92,
                                       1,
                                       92,
                                       78,
                                       -7,
                                       7);
 
   TrueBCOccupancy_ = new TH2F(
       "BCTrueOccupancyAllEvents", "True SB Occupancy all events;#phi;#eta", 360, -3.14159, 3.14159, 172, -1.5, 1.5);
   TrueBCOccupancyCoarse_ = new TH2F("BCTrueOccupancyAllEventsCoarse",
                                     "True BC Occupancy all events Coarse;#phi;#eta",
                                     360 / 5,
                                     -3.14159,
                                     3.14159,
                                     34,
                                     -1.5,
                                     1.5);
 
   numberofBCinSC_ =
       new TH1F("numberofBCinSC", "Number of Basic Clusters in Super Cluster;Num Basic Clusters", 20, 0, 20);  //SC
   numberofBCinSCphi_ = new TH2F("numberofBCinSCphi",
                                 "Number of Basic Clusters in Super Cluster;phi;Num Basic Clusters",
                                 360 / 5,
                                 -3.14159,
                                 3.14159,
                                 20,
                                 0,
                                 20);  //SC
 
   allFedsTimingHistEcalMuon_ =
       new TH1F("timeForAllFeds_EcalMuon", "timeForAllFeds Ecal+Muon;Relative Time (1 clock = 25ns)", 78, -7, 7);
 
   triggerHist_ = new TH1F("triggerHist", "Trigger Number", 5, 0, 5);
   triggerHist_->GetXaxis()->SetBinLabel(1, "ECAL");
   triggerHist_->GetXaxis()->SetBinLabel(2, "HCAL");
   triggerHist_->GetXaxis()->SetBinLabel(3, "DT");
   triggerHist_->GetXaxis()->SetBinLabel(4, "RPC");
   triggerHist_->GetXaxis()->SetBinLabel(5, "CSC");
 
   triggerExclusiveHist_ = new TH1F("triggerExclusiveHist", "Trigger Number (Mutually Exclusive)", 5, 0, 5);
   triggerExclusiveHist_->GetXaxis()->SetBinLabel(1, "ECAL");
   triggerExclusiveHist_->GetXaxis()->SetBinLabel(2, "HCAL");
   triggerExclusiveHist_->GetXaxis()->SetBinLabel(3, "DT");
   triggerExclusiveHist_->GetXaxis()->SetBinLabel(4, "RPC");
   triggerExclusiveHist_->GetXaxis()->SetBinLabel(5, "CSC");
 
   runNumberHist_ = new TH1F("runNumberHist", "Run Number", 1, 0, 1);
 
   deltaRHist_ = new TH1F("deltaRHist", "deltaR", 500, -0.5, 499.5);
   deltaEtaHist_ = new TH1F("deltaIEtaHist", "deltaIEta", 170, -85.5, 84.5);
   deltaPhiHist_ = new TH1F("deltaIPhiHist", "deltaIPhi", 720, -360.5, 359.5);
   ratioAssocTracksHist_ = new TH1F("ratioAssocTracks", "num assoc. tracks/tracks through Ecal", 11, 0, 1.1);
   ratioAssocClustersHist_ = new TH1F("ratioAssocClusters", "num assoc. clusters/total clusters", 11, 0, 1.1);
   trackAssoc_muonsEcal_ = new TH2F(
       "trackAssoc_muonsEcal", "Map of muon hits in Ecal", 360, 1., 361., 172, -86, 86);  //360, 0 , 360, 170,-85 ,85);
   deltaEtaDeltaPhiHist_ =
       new TH2F("deltaEtaDeltaPhi", "Delta ieta vs. delta iphi", 170, -85.5, 84.5, 720, -360.5, 359.5);
   seedTrackEtaHist_ = new TH2F("seedTrackEta", "track ieta vs. seed ieta", 170, -85.5, 84.5, 170, -85.5, 84.5);
   seedTrackPhiHist_ = new TH2F("seedTrackPhi", "track iphi vs. seed iphi", 720, -360.5, 359.5, 720, -360.5, 359.5);
 
   dccEventVsBxHist_ = new TH2F("dccEventVsBx", "DCC Runtype vs. bunch crossing", 3600, 0, 3600, 24, 0, 24);
   dccBXErrorByFEDHist_ = new TH1F("dccBXErrorByFED", "Incorrect BX number by FED", 54, 601, 655);
   dccOrbitErrorByFEDHist_ = new TH1F("dccOrbitErrorByFED", "Incorrect orbit number by FED", 54, 601, 655);
   dccRuntypeErrorByFEDHist_ = new TH1F("dccRuntypeErrorByFED", "Incorrect DCC Runtype by FED", 54, 601, 655);
   dccRuntypeHist_ = new TH1F("dccRuntype", "DCC Runtype frequency", 24, 0, 24);
   dccErrorVsBxHist_ = new TH2F("dccErrorVsBX", "DCC Errors vs. BX", 3600, 0, 3600, 3, 0, 3);
 
   hcalEnergy_HBHE_ = new TH1F("hcalEnergy_HBHE", "RecHit Energy HBHE", 440, -10, 100);
   hcalEnergy_HF_ = new TH1F("hcalEnergy_HF", "RecHit Energy HF", 440, -10, 100);
   hcalEnergy_HO_ = new TH1F("hcalEnergy_HO", "RecHit Energy HO", 440, -10, 100);
   hcalHEHBecalEB_ =
       new TH2F("hcalHEHBecalEB", "HCAL HBHE RecHit energy vs ECAL EB energy", numBins, histRangeMin_, 300.0, 40, -5, 5);
 
   NumXtalsInClusterHist_ = new TH1F("NumXtalsInClusterAllHist", "Number of Xtals in Cluster;NumXtals", 150, 0, 150);
   numxtalsVsEnergy_ = new TH2F("NumXtalsVsEnergy",
                                "Number of Xtals in Cluster vs Energy;Energy (GeV);Number of Xtals in Cluster",
                                numBins,
                                histRangeMin_,
                                histRangeMax_,
                                150,
                                0,
                                150);
   numxtalsVsHighEnergy_ = new TH2F("NumXtalsVsHighEnergy",
                                    "Number of Xtals in Cluster vs Energy;Energy (GeV);Number of Xtals in Cluster",
                                    numBins,
                                    histRangeMin_,
                                    200.,
                                    150,
                                    0,
                                    150);
 
   // high energy analysis
   allOccupancyHighEnergy_ =
       new TH2F("OccupancyHighEnergyEvents", "Occupancy high energy events;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
   allOccupancyHighEnergyCoarse_ = new TH2F("OccupancyHighEnergyEventsCoarse",
                                            "Occupancy high energy events Coarse;i#phi;i#eta",
                                            72,
                                            ttPhiBins,
                                            35,
                                            ttEtaBins);
   allFedsOccupancyHighEnergyHist_ = new TH3F("OccupancyHighEnergyEvents3D",
                                              "(Phi,Eta,energy) for all high energy events;i#phi;i#eta;energy (GeV)",
                                              18,
                                              modPhiBins,
                                              9,
                                              modEtaBins,
                                              10,
                                              highEBins);
   allFedsNumXtalsInClusterHist_ =
       new TH1F("NumActiveXtalsInClusterAllHist", "Number of active Xtals in Cluster;NumXtals", 100, 0, 100);
 
   HighEnergy_NumXtal = new TH1F("HighEnergy_NumXtal", "Num crystals in high E clusters;num crystals", 150, 0, 150);
   HighEnergy_NumXtalFedId = new TH2F(
       "HighEnergy_NumXtalFedId", "Num crystals in cluster vs FedId;FedId;num crystals", 36, 610., 645., 150, 0, 150);
   HighEnergy_NumXtaliphi = new TH2F(
       "HighEnergy_NumXtaliphi", "Num crystals in cluster vs iphi;i#phi;num crystals", 360, 1., 361., 150, 0, 150);
   HighEnergy_energy3D = new TH3F("HighEnergy_energy3D",
                                  "(Phi,Eta,energy) for all high energy events;i#phi;i#eta;energy (GeV)",
                                  72,
                                  ttPhiBins,
                                  35,
                                  ttEtaBins,
                                  10,
                                  highEBins);
 
   HighEnergy_energyNumXtal = new TH2F("HighEnergy_energyNumXtal",
                                       "Energy in cluster vs Num crystals in cluster;num crystals;energy",
                                       150,
                                       0,
                                       150,
                                       200,
                                       0.,
                                       200.);
 
   HighEnergy_bestSeed = new TH1F("HighEnergy_bestSeedEnergy", "BestSeed Energy from TrackAss", 200, 0., 200.);
   HighEnergy_bestSeedOccupancy = new TH2F(
       "HighEnergy_bestSeedOccupancy", "Occupancy HighEn events from TrackAss;i#phi;i#eta", 360, 1., 361., 172, -86, 86);
 
   HighEnergy_numClusHighEn = new TH1F("HighEnergy_numClusHighEn", "Num High Energy Clusters", 7, 0, 7);
   HighEnergy_ratioClusters =
       new TH1F("HighEnergy_ratioClusters", "Num High Energy Clusters/Num tot Clusters", 100, 0., 1.1);
 
   HighEnergy_numRecoTrackBarrel = new TH1F("HighEnergy_numRecoTracksBarrel", "Num BarrelRecoTracks", 10, 0, 10);
   HighEnergy_TracksAngle = new TH1F("HighEnergy_TracksAngle", "Angle between tracks", 720, 0., 180.);
   HighEnergy_TracksAngleTopBottom =
       new TH1F("HighEnergy_TopBottomTracksAngle", "Angle between top-bottom tracks", 720, 0., 180.);
 
   HighEnergy_2GeV_occuCoarse = new TH2F("HighEnergy_occu2GeV_Coarse",
                                         "Occupancy high energy events with more than 2 GeV;i#phi;i#eta",
                                         72,
                                         ttPhiBins,
                                         35,
                                         ttEtaBins);
   HighEnergy_2GeV_occu3D = new TH3F("HighEnergy_2GeV_energy3D",
                                     "(iphi,ieta,energy) for all high energy events w > 10 GeV;i#phi;i#eta;energy (GeV)",
                                     72,
                                     ttPhiBins,
                                     35,
                                     ttEtaBins,
                                     10,
                                     highEBins);
   HighEnergy_100GeV_occuCoarse = new TH2F("HighEnergy_occu100GeV_Coarse",
                                           "Occupancy high energy events with more than 100 GeV;i#phi;i#eta",
                                           72,
                                           ttPhiBins,
                                           35,
                                           ttEtaBins);
   HighEnergy_100GeV_occu3D =
       new TH3F("HighEnergy_100GeV_energy3D",
                "(iphi,ieta,energy) for all high energy events more than 100 GeV;i#phi;i#eta;energy (GeV)",
                72,
                ttPhiBins,
                35,
                ttEtaBins,
                10,
                highEBins);
   HighEnergy_0tracks_occu3D = new TH3F("HighEnergy_0Tracks_energy3D",
                                        "(iphi,ieta,energy) for all events with 0 tracks;i#phi;i#eta;energy (GeV)",
                                        72,
                                        ttPhiBins,
                                        35,
                                        ttEtaBins,
                                        10,
                                        highEBins);
   HighEnergy_1tracks_occu3D = new TH3F("HighEnergy_1Tracks_energy3D",
                                        "(iphi,ieta,energy) for all events with 1 tracks;i#phi;i#eta;energy (GeV)",
                                        72,
                                        ttPhiBins,
                                        35,
                                        ttEtaBins,
                                        10,
                                        highEBins);
   HighEnergy_2tracks_occu3D = new TH3F("HighEnergy_2Tracks_energy3D",
                                        "(iphi,ieta,energy) for all events with 2 tracks;i#phi;i#eta;energy (GeV)",
                                        72,
                                        ttPhiBins,
                                        35,
                                        ttEtaBins,
                                        10,
                                        highEBins);
   HighEnergy_0tracks_occu3DXtal = new TH3F("HighEnergy_0Tracks_energy3DXtal",
                                            "(iphi,ieta,energy) for all events with 0 tracks;i#phi;i#eta;energy (GeV)",
                                            360,
                                            1.,
                                            361.,
                                            172,
                                            -86,
                                            86,
                                            200,
                                            0.,
                                            200.);
   HighEnergy_1tracks_occu3DXtal = new TH3F("HighEnergy_1Tracks_energy3DXtal",
                                            "(iphi,ieta,energy) for all events with 1 tracks;i#phi;i#eta;energy (GeV)",
                                            360,
                                            1.,
                                            361.,
                                            172,
                                            -86,
                                            86,
                                            200,
                                            0.,
                                            200.);
   HighEnergy_2tracks_occu3DXtal = new TH3F("HighEnergy_2Tracks_energy3DXtal",
                                            "(iphi,ieta,energy) for all events with 2 tracks;i#phi;i#eta;energy (GeV)",
                                            360,
                                            1.,
                                            361.,
                                            172,
                                            -86,
                                            86,
                                            200,
                                            0.,
                                            200.);
 
   //EE histograms
 
   // EE-
   EEM_FedsSeedEnergyHist_ =
       new TH1F("SeedEnergyAllFEDs", "Seed Energy for EEM Feds; Seed Energy (GeV)", 200, histRangeMin_, 10.0);
 
   EEM_AllOccupancyCoarse_ =
       new TH2F("OccupancyAllEventsCoarse", "Occupancy all events Coarse EEM;ix;iy", 20, 0, 100, 20, 0, 100);
   EEM_AllOccupancy_ = new TH2F("OccupancyAllEvents", "Occupancy all events EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_FedsenergyHist_ =
       new TH1F("energy_AllClusters", "energy_AllClusters_EEM;Cluster Energy (GeV)", numBins, histRangeMin_, 10.0);
   EEM_FedsenergyHighHist_ = new TH1F(
       "energyHigh_AllClusters", "energyHigh_AllClusters in EEM;Cluster Energy (GeV)", numBins, histRangeMin_, 200.0);
   EEM_FedsenergyOnlyHighHist_ = new TH1F("energyHigh_HighEnergyClusters",
                                          "energy of High Energy Clusters in EEM;Cluster Energy (GeV)",
                                          numBins,
                                          histRangeMin_,
                                          200.0);
   EEM_FedsE2Hist_ =
       new TH1F("E2_AllClusters", "E2_AllClusters_EEM;Seed+highest neighbor energy (GeV)", numBins, histRangeMin_, 10.0);
   EEM_FedsE2vsE1Hist_ = new TH2F("E2vsE1_AllClusters",
                                  "E2vsE1_AllClusters_EEM;Seed Energy (GeV);Seed+highest neighbor energy (GeV)",
                                  numBins,
                                  histRangeMin_,
                                  10.0,
                                  numBins,
                                  histRangeMin_,
                                  10.0);
   EEM_FedsenergyvsE1Hist_ = new TH2F("energyvsE1_AllClusters",
                                      "energyvsE1_AllClusters_EEM;Seed Energy (GeV);Energy(GeV)",
                                      numBins,
                                      histRangeMin_,
                                      10.0,
                                      numBins,
                                      histRangeMin_,
                                      10.0);
   EEM_FedsTimingHist_ = new TH1F("timeForAllFeds", "timeForAllFeds_EEM;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEM_numberofCosmicsHist_ =
       new TH1F("numberofCosmicsPerEvent", "Number of cosmics per event EEM;Number of Cosmics", 30, 0, 30);
 
   EEM_FedsTimingVsAmpHist_ = new TH2F("timeVsAmpAllEvents",
                                       "time Vs Amp All Events EEM;Relative Time (1 clock = 25ns);Amplitude (GeV)",
                                       78,
                                       -7,
                                       7,
                                       numBins,
                                       histRangeMin_,
                                       10.0);
   EEM_FedsTimingTTHist_ = new TH3F("timeTTAllFEDs",
                                    "(ix,iy,time) for all FEDs (SM,TT binning) EEM;ix;iy;Relative Time (1 clock = 25ns)",
                                    20,
                                    0,
                                    100,
                                    20,
                                    0,
                                    100,
                                    78,
                                    -7,
                                    7);
 
   EEM_OccupancySingleXtal_ =
       new TH2F("OccupancySingleXtal", "Occupancy single xtal clusters EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_energySingleXtalHist_ = new TH1F("energy_SingleXtalClusters",
                                        "Energy single xtal clusters EEM;Cluster Energy (GeV)",
                                        numBins,
                                        histRangeMin_,
                                        200.0);
 
   EEM_OccupancyExclusiveECAL_ = new TH2F(
       "OccupancyAllEvents_ExclusiveECAL", "Occupancy all events Exclusive ECAL  EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseExclusiveECAL_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveECAL",
                                                "Occupancy all events Coarse Exclusive ECAL EEM;ix;iy",
                                                20,
                                                0,
                                                100,
                                                20,
                                                0,
                                                100);
   EEM_OccupancyECAL_ =
       new TH2F("OccupancyAllEvents_ECAL", "Occupancy all events ECAL EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseECAL_ =
       new TH2F("OccupancyAllEventsCoarse_ECAL", "Occupancy all events Coarse ECAL EEM;ix;iy", 20, 0, 100, 20, 0, 100);
   EEM_FedsTimingHistECAL_ =
       new TH1F("timeForAllFeds_ECAL", "timeForAllFeds ECAL EEM;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEM_FedsTimingTTHistECAL_ =
       new TH3F("timeTTAllFEDs_ECAL",
                "(ix,iy,time) for all FEDs (SM,TT binning) ECAL EEM;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEM_OccupancyExclusiveDT_ = new TH2F(
       "OccupancyAllEvents_ExclusiveDT", "Occupancy all events Exclusive DT EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseExclusiveDT_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveDT",
                                              "Occupancy all events Coarse Exclusive DT EEM;ix;iy",
                                              20,
                                              0,
                                              100,
                                              20,
                                              0,
                                              100);
   EEM_OccupancyDT_ = new TH2F("OccupancyAllEvents_DT", "Occupancy all events DT EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseDT_ =
       new TH2F("OccupancyAllEventsCoarse_DT", "Occupancy all events Coarse DT EEM;ix;iy", 20, 0, 100, 20, 0, 100);
   EEM_FedsTimingHistDT_ =
       new TH1F("timeForAllFeds_DT", "timeForAllFeds DT EEM;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEM_FedsTimingTTHistDT_ =
       new TH3F("timeTTAllFEDs_DT",
                "(ix,iy,time) for all FEDs (SM,TT binning) DT EEM;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEM_OccupancyExclusiveRPC_ = new TH2F(
       "OccupancyAllEvents_ExclusiveRPC", "Occupancy all events Exclusive RPC EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseExclusiveRPC_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveRPC",
                                               "Occupancy all events Coarse Exclusive RPC EEM;ix;iy",
                                               20,
                                               0,
                                               100,
                                               20,
                                               0,
                                               100);
   EEM_OccupancyRPC_ =
       new TH2F("OccupancyAllEvents_RPC", "Occupancy all events RPC EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseRPC_ =
       new TH2F("OccupancyAllEventsCoarse_RPC", "Occupancy all events Coarse RPC EEM;ix;iy", 20, 0, 100, 20, 0, 100);
   EEM_FedsTimingHistRPC_ =
       new TH1F("timeForAllFeds_RPC", "timeForAllFeds RPC EEM;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEM_FedsTimingTTHistRPC_ =
       new TH3F("timeTTAllFEDs_RPC",
                "(ix,iy,time) for all FEDs (SM,TT binning) RPC EEM;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEM_OccupancyExclusiveCSC_ = new TH2F(
       "OccupancyAllEvents_ExclusiveCSC", "Occupancy all events Exclusive CSC EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseExclusiveCSC_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveCSC",
                                               "Occupancy all events Coarse Exclusive CSC EEM;ix;iy",
                                               20,
                                               0,
                                               100,
                                               20,
                                               0,
                                               100);
   EEM_OccupancyCSC_ =
       new TH2F("OccupancyAllEvents_CSC", "Occupancy all events CSC EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseCSC_ =
       new TH2F("OccupancyAllEventsCoarse_CSC", "Occupancy all events Coarse CSC EEM;ix;iy", 20, 0, 100, 20, 0, 100);
   EEM_FedsTimingHistCSC_ =
       new TH1F("timeForAllFeds_CSC", "timeForAllFeds CSC EEM;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEM_FedsTimingTTHistCSC_ =
       new TH3F("timeTTAllFEDs_CSC",
                "(ix,iy,time) for all FEDs (SM,TT binning) CSC EEM;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEM_OccupancyExclusiveHCAL_ = new TH2F(
       "OccupancyAllEvents_ExclusiveHCAL", "Occupancy all events Exclusive HCAL EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseExclusiveHCAL_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveHCAL",
                                                "Occupancy all events Coarse Exclusive HCAL EEM;ix;iy",
                                                20,
                                                0,
                                                100,
                                                20,
                                                0,
                                                100);
   EEM_OccupancyHCAL_ =
       new TH2F("OccupancyAllEvents_HCAL", "Occupancy all events HCAL EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyCoarseHCAL_ =
       new TH2F("OccupancyAllEventsCoarse_HCAL", "Occupancy all events Coarse HCAL EEM;ix;iy", 20, 0, 100, 20, 0, 100);
   EEM_FedsTimingHistHCAL_ =
       new TH1F("timeForAllFeds_HCAL", "timeForAllFeds HCAL EEM;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEM_FedsTimingTTHistHCAL_ =
       new TH3F("timeTTAllFEDs_HCAL",
                "(ix,iy,time) for all FEDs (SM,TT binning) HCAL EEM;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEM_numberofBCinSC_ =
       new TH1F("numberofBCinSC", "Number of Basic Clusters in Super Cluster EEM;Num Basic Clusters", 20, 0, 20);
 
   EEM_triggerHist_ = new TH1F("triggerHist", "Trigger Number EEM", 5, 0, 5);
   EEM_triggerHist_->GetXaxis()->SetBinLabel(1, "ECAL");
   EEM_triggerHist_->GetXaxis()->SetBinLabel(2, "HCAL");
   EEM_triggerHist_->GetXaxis()->SetBinLabel(3, "DT");
   EEM_triggerHist_->GetXaxis()->SetBinLabel(4, "RPC");
   EEM_triggerHist_->GetXaxis()->SetBinLabel(5, "CSC");
 
   EEM_triggerExclusiveHist_ = new TH1F("triggerExclusiveHist", "Trigger Number (Mutually Exclusive) EEM", 5, 0, 5);
   triggerExclusiveHist_->GetXaxis()->SetBinLabel(1, "ECAL");
   EEM_triggerExclusiveHist_->GetXaxis()->SetBinLabel(2, "HCAL");
   EEM_triggerExclusiveHist_->GetXaxis()->SetBinLabel(3, "DT");
   EEM_triggerExclusiveHist_->GetXaxis()->SetBinLabel(4, "RPC");
   EEM_triggerExclusiveHist_->GetXaxis()->SetBinLabel(5, "CSC");
 
   EEM_NumXtalsInClusterHist_ =
       new TH1F("NumXtalsInClusterAllHist", "Number of Xtals in Cluster EEM;NumXtals", 150, 0, 150);
   EEM_numxtalsVsEnergy_ = new TH2F("NumXtalsVsEnergy",
                                    "Number of Xtals in Cluster vs Energy EEM;Energy (GeV);Number of Xtals in Cluster",
                                    numBins,
                                    histRangeMin_,
                                    10.0,
                                    150,
                                    0,
                                    150);
   EEM_numxtalsVsHighEnergy_ =
       new TH2F("NumXtalsVsHighEnergy",
                "Number of Xtals in Cluster vs Energy EEM;Energy (GeV);Number of Xtals in Cluster",
                numBins,
                histRangeMin_,
                200.,
                150,
                0,
                150);
 
   EEM_OccupancyHighEnergy_ =
       new TH2F("OccupancyHighEnergyEvents", "Occupancy high energy events EEM;ix;iy", 100, 0, 100, 100, 0, 100);
   EEM_OccupancyHighEnergyCoarse_ = new TH2F(
       "OccupancyHighEnergyEventsCoarse", "Occupancy high energy events Coarse EEM;ix;iy", 20, 0, 100, 20, 0, 100);
 
   EEM_FedsNumXtalsInClusterHist_ =
       new TH1F("NumActiveXtalsInClusterAllHist", "Number of active Xtals in Cluster EEM;NumXtals", 100, 0, 100);
 
   // EE+
   EEP_FedsSeedEnergyHist_ =
       new TH1F("SeedEnergyAllFEDs", "Seed Energy for EEP Feds; Seed Energy (GeV)", 200, histRangeMin_, 10.0);
 
   EEP_AllOccupancyCoarse_ =
       new TH2F("OccupancyAllEventsCoarse", "Occupancy all events Coarse EEP;ix;iy", 20, 0, 100, 20, 0, 100);
   EEP_AllOccupancy_ = new TH2F("OccupancyAllEvents", "Occupancy all events EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_FedsenergyHist_ =
       new TH1F("energy_AllClusters", "energy_AllClusters_EEP;Cluster Energy (GeV)", numBins, histRangeMin_, 10.0);
   EEP_FedsenergyHighHist_ = new TH1F(
       "energyHigh_AllClusters", "energyHigh_AllClusters in EEP;Cluster Energy (GeV)", numBins, histRangeMin_, 200.0);
   EEP_FedsenergyOnlyHighHist_ = new TH1F("energyHigh_HighEnergyClusters",
                                          "energy of High Energy Clusters in EEP;Cluster Energy (GeV)",
                                          numBins,
                                          histRangeMin_,
                                          200.0);
   EEP_FedsE2Hist_ =
       new TH1F("E2_AllClusters", "E2_AllClusters_EEP;Seed+highest neighbor energy (GeV)", numBins, histRangeMin_, 10.0);
   EEP_FedsE2vsE1Hist_ = new TH2F("E2vsE1_AllClusters",
                                  "E2vsE1_AllClusters_EEP;Seed Energy (GeV);Seed+highest neighbor energy (GeV)",
                                  numBins,
                                  histRangeMin_,
                                  10.0,
                                  numBins,
                                  histRangeMin_,
                                  10.0);
   EEP_FedsenergyvsE1Hist_ = new TH2F("energyvsE1_AllClusters",
                                      "energyvsE1_AllClusters_EEP;Seed Energy (GeV);Energy(GeV)",
                                      numBins,
                                      histRangeMin_,
                                      10.0,
                                      numBins,
                                      histRangeMin_,
                                      10.0);
   EEP_FedsTimingHist_ = new TH1F("timeForAllFeds", "timeForAllFeds_EEP;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEP_numberofCosmicsHist_ =
       new TH1F("numberofCosmicsPerEvent", "Number of cosmics per event EEP;Number of Cosmics", 30, 0, 30);
 
   EEP_FedsTimingVsAmpHist_ = new TH2F("timeVsAmpAllEvents",
                                       "time Vs Amp All Events EEP;Relative Time (1 clock = 25ns);Amplitude (GeV)",
                                       78,
                                       -7,
                                       7,
                                       numBins,
                                       histRangeMin_,
                                       10.0);
   EEP_FedsTimingTTHist_ = new TH3F("timeTTAllFEDs",
                                    "(ix,iy,time) for all FEDs (SM,TT binning) EEP;ix;iy;Relative Time (1 clock = 25ns)",
                                    20,
                                    0,
                                    100,
                                    20,
                                    0,
                                    100,
                                    78,
                                    -7,
                                    7);
 
   EEP_OccupancySingleXtal_ =
       new TH2F("OccupancySingleXtal", "Occupancy single xtal clusters EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_energySingleXtalHist_ = new TH1F("energy_SingleXtalClusters",
                                        "Energy single xtal clusters EEP;Cluster Energy (GeV)",
                                        numBins,
                                        histRangeMin_,
                                        200.0);
 
   EEP_OccupancyExclusiveECAL_ = new TH2F(
       "OccupancyAllEvents_ExclusiveECAL", "Occupancy all events Exclusive ECAL  EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseExclusiveECAL_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveECAL",
                                                "Occupancy all events Coarse Exclusive ECAL EEP;ix;iy",
                                                20,
                                                0,
                                                100,
                                                20,
                                                0,
                                                100);
   EEP_OccupancyECAL_ =
       new TH2F("OccupancyAllEvents_ECAL", "Occupancy all events ECAL EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseECAL_ =
       new TH2F("OccupancyAllEventsCoarse_ECAL", "Occupancy all events Coarse ECAL EEP;ix;iy", 20, 0, 100, 20, 0, 100);
   EEP_FedsTimingHistECAL_ =
       new TH1F("timeForAllFeds_ECAL", "timeForAllFeds ECAL EEP;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEP_FedsTimingTTHistECAL_ =
       new TH3F("timeTTAllFEDs_ECAL",
                "(ix,iy,time) for all FEDs (SM,TT binning) ECAL EEP;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEP_OccupancyExclusiveDT_ = new TH2F(
       "OccupancyAllEvents_ExclusiveDT", "Occupancy all events Exclusive DT EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseExclusiveDT_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveDT",
                                              "Occupancy all events Coarse Exclusive DT EEP;ix;iy",
                                              20,
                                              0,
                                              100,
                                              20,
                                              0,
                                              100);
   EEP_OccupancyDT_ = new TH2F("OccupancyAllEvents_DT", "Occupancy all events DT EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseDT_ =
       new TH2F("OccupancyAllEventsCoarse_DT", "Occupancy all events Coarse DT EEP;ix;iy", 20, 0, 100, 20, 0, 100);
   EEP_FedsTimingHistDT_ =
       new TH1F("timeForAllFeds_DT", "timeForAllFeds DT EEP;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEP_FedsTimingTTHistDT_ =
       new TH3F("timeTTAllFEDs_DT",
                "(ix,iy,time) for all FEDs (SM,TT binning) DT EEP;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEP_OccupancyExclusiveRPC_ = new TH2F(
       "OccupancyAllEvents_ExclusiveRPC", "Occupancy all events Exclusive RPC EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseExclusiveRPC_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveRPC",
                                               "Occupancy all events Coarse Exclusive RPC EEP;ix;iy",
                                               20,
                                               0,
                                               100,
                                               20,
                                               0,
                                               100);
   EEP_OccupancyRPC_ =
       new TH2F("OccupancyAllEvents_RPC", "Occupancy all events RPC EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseRPC_ =
       new TH2F("OccupancyAllEventsCoarse_RPC", "Occupancy all events Coarse RPC EEP;ix;iy", 20, 0, 100, 20, 0, 100);
   EEP_FedsTimingHistRPC_ =
       new TH1F("timeForAllFeds_RPC", "timeForAllFeds RPC EEP;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEP_FedsTimingTTHistRPC_ =
       new TH3F("timeTTAllFEDs_RPC",
                "(ix,iy,time) for all FEDs (SM,TT binning) RPC EEP;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEP_OccupancyExclusiveCSC_ = new TH2F(
       "OccupancyAllEvents_ExclusiveCSC", "Occupancy all events Exclusive CSC EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseExclusiveCSC_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveCSC",
                                               "Occupancy all events Coarse Exclusive CSC EEP;ix;iy",
                                               20,
                                               0,
                                               100,
                                               20,
                                               0,
                                               100);
   EEP_OccupancyCSC_ =
       new TH2F("OccupancyAllEvents_CSC", "Occupancy all events CSC EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseCSC_ =
       new TH2F("OccupancyAllEventsCoarse_CSC", "Occupancy all events Coarse CSC EEP;ix;iy", 20, 0, 100, 20, 0, 100);
   EEP_FedsTimingHistCSC_ =
       new TH1F("timeForAllFeds_CSC", "timeForAllFeds CSC EEP;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEP_FedsTimingTTHistCSC_ =
       new TH3F("timeTTAllFEDs_CSC",
                "(ix,iy,time) for all FEDs (SM,TT binning) CSC EEP;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEP_OccupancyExclusiveHCAL_ = new TH2F(
       "OccupancyAllEvents_ExclusiveHCAL", "Occupancy all events Exclusive HCAL EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseExclusiveHCAL_ = new TH2F("OccupancyAllEventsCoarse_ExclusiveHCAL",
                                                "Occupancy all events Coarse Exclusive HCAL EEP;ix;iy",
                                                20,
                                                0,
                                                100,
                                                20,
                                                0,
                                                100);
   EEP_OccupancyHCAL_ =
       new TH2F("OccupancyAllEvents_HCAL", "Occupancy all events HCAL EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyCoarseHCAL_ =
       new TH2F("OccupancyAllEventsCoarse_HCAL", "Occupancy all events Coarse HCAL EEP;ix;iy", 20, 0, 100, 20, 0, 100);
   EEP_FedsTimingHistHCAL_ =
       new TH1F("timeForAllFeds_HCAL", "timeForAllFeds HCAL EEP;Relative Time (1 clock = 25ns)", 78, -7, 7);
   EEP_FedsTimingTTHistHCAL_ =
       new TH3F("timeTTAllFEDs_HCAL",
                "(ix,iy,time) for all FEDs (SM,TT binning) HCAL EEP;ix;iy;Relative Time (1 clock = 25ns)",
                20,
                0,
                100,
                20,
                0,
                100,
                78,
                -7,
                7);
 
   EEP_numberofBCinSC_ =
       new TH1F("numberofBCinSC", "Number of Basic Clusters in Super Cluster EEP;Num Basic Clusters", 20, 0, 20);
 
   EEP_triggerHist_ = new TH1F("triggerHist", "Trigger Number EEP", 5, 0, 5);
   EEP_triggerHist_->GetXaxis()->SetBinLabel(1, "ECAL");
   EEP_triggerHist_->GetXaxis()->SetBinLabel(2, "HCAL");
   EEP_triggerHist_->GetXaxis()->SetBinLabel(3, "DT");
   EEP_triggerHist_->GetXaxis()->SetBinLabel(4, "RPC");
   EEP_triggerHist_->GetXaxis()->SetBinLabel(5, "CSC");
 
   EEP_triggerExclusiveHist_ = new TH1F("triggerExclusiveHist", "Trigger Number (Mutually Exclusive) EEP", 5, 0, 5);
   EEP_triggerExclusiveHist_->GetXaxis()->SetBinLabel(1, "ECAL");
   EEP_triggerExclusiveHist_->GetXaxis()->SetBinLabel(2, "HCAL");
   EEP_triggerExclusiveHist_->GetXaxis()->SetBinLabel(3, "DT");
   EEP_triggerExclusiveHist_->GetXaxis()->SetBinLabel(4, "RPC");
   EEP_triggerExclusiveHist_->GetXaxis()->SetBinLabel(5, "CSC");
 
   EEP_NumXtalsInClusterHist_ =
       new TH1F("NumXtalsInClusterAllHist", "Number of Xtals in Cluster EEP;NumXtals", 150, 0, 150);
   EEP_numxtalsVsEnergy_ = new TH2F("NumXtalsVsEnergy",
                                    "Number of Xtals in Cluster vs Energy EEP;Energy (GeV);Number of Xtals in Cluster",
                                    numBins,
                                    histRangeMin_,
                                    10.0,
                                    150,
                                    0,
                                    150);
   EEP_numxtalsVsHighEnergy_ =
       new TH2F("NumXtalsVsHighEnergy",
                "Number of Xtals in Cluster vs Energy EEP;Energy (GeV);Number of Xtals in Cluster",
                numBins,
                histRangeMin_,
                200.,
                150,
                0,
                150);
 
   EEP_OccupancyHighEnergy_ =
       new TH2F("OccupancyHighEnergyEvents", "Occupancy high energy events EEP;ix;iy", 100, 0, 100, 100, 0, 100);
   EEP_OccupancyHighEnergyCoarse_ = new TH2F(
       "OccupancyHighEnergyEventsCoarse", "Occupancy high energy events Coarse EEP;ix;iy", 20, 0, 100, 20, 0, 100);
 
   EEP_FedsNumXtalsInClusterHist_ =
       new TH1F("NumActiveXtalsInClusterAllHist", "Number of active Xtals in Cluster EEP;NumXtals", 100, 0, 100);
 }

std::vector< bool > EcalCosmicsHists::determineTriggers	(	const edm::Event &	iEvent,
		const edm::EventSetup &	eventSetup
	)

private

Definition at line 1307 of file EcalCosmicsHists.cc.

References L1GtPsbWord::aData(), L1GtPsbWord::bData(), L1GlobalTriggerReadoutRecord::decisionWord(), edm::EventSetup::get(), edm::Event::getByLabel(), L1MuGMTReadoutCollection::getRecords(), L1GtTriggerMenu::gtAlgorithmResult(), L1GlobalTriggerReadoutRecord::gtPsbWord(), mps_fire::i, edm::HandleBase::isValid(), l1GMTReadoutRecTag_, l1GTReadoutRecTag_, edm::Handle< T >::product(), and edm::ESHandle< T >::product().

Referenced by analyze().

                                                                                                            {
   std::vector<bool> l1Triggers;  //DT,CSC,RPC,HCAL,ECAL
                                  //0 , 1 , 2 , 3  , 4
   for (int i = 0; i < 5; i++)
     l1Triggers.push_back(false);
 
   // get the GMTReadoutCollection
   Handle<L1MuGMTReadoutCollection> gmtrc_handle;
   iEvent.getByLabel(l1GMTReadoutRecTag_, gmtrc_handle);
   L1MuGMTReadoutCollection const* gmtrc = gmtrc_handle.product();
   if (!(gmtrc_handle.isValid())) {
     LogWarning("EcalCosmicsHists") << "l1MuGMTReadoutCollection"
                                    << " not available";
     return l1Triggers;
   }
   // get hold of L1GlobalReadoutRecord
   Handle<L1GlobalTriggerReadoutRecord> L1GTRR;
   iEvent.getByLabel(l1GTReadoutRecTag_, L1GTRR);
 
   //Ecal
   edm::ESHandle<L1GtTriggerMenu> menuRcd;
   eventSetup.get<L1GtTriggerMenuRcd>().get(menuRcd);
   const L1GtTriggerMenu* menu = menuRcd.product();
   edm::Handle<L1GlobalTriggerReadoutRecord> gtRecord;
   iEvent.getByLabel(edm::InputTag("gtDigis"), gtRecord);
   // Get dWord after masking disabled bits
   const DecisionWord dWord = gtRecord->decisionWord();
 
   bool l1SingleEG1 = menu->gtAlgorithmResult("L1_SingleEG1", dWord);
   bool l1SingleEG5 = menu->gtAlgorithmResult("L1_SingleEG5", dWord);
   bool l1SingleEG8 = menu->gtAlgorithmResult("L1_SingleEG8", dWord);
   bool l1SingleEG10 = menu->gtAlgorithmResult("L1_SingleEG10", dWord);
   bool l1SingleEG12 = menu->gtAlgorithmResult("L1_SingleEG12", dWord);
   bool l1SingleEG15 = menu->gtAlgorithmResult("L1_SingleEG15", dWord);
   bool l1SingleEG20 = menu->gtAlgorithmResult("L1_SingleEG20", dWord);
   bool l1SingleEG25 = menu->gtAlgorithmResult("L1_SingleEG25", dWord);
   bool l1DoubleNoIsoEGBTBtight = menu->gtAlgorithmResult("L1_DoubleNoIsoEG_BTB_tight", dWord);
   bool l1DoubleNoIsoEGBTBloose = menu->gtAlgorithmResult("L1_DoubleNoIsoEG_BTB_loose ", dWord);
   bool l1DoubleNoIsoEGTopBottom = menu->gtAlgorithmResult("L1_DoubleNoIsoEGTopBottom", dWord);
   bool l1DoubleNoIsoEGTopBottomCen = menu->gtAlgorithmResult("L1_DoubleNoIsoEGTopBottomCen", dWord);
   bool l1DoubleNoIsoEGTopBottomCen2 = menu->gtAlgorithmResult("L1_DoubleNoIsoEGTopBottomCen2", dWord);
   bool l1DoubleNoIsoEGTopBottomCenVert = menu->gtAlgorithmResult("L1_DoubleNoIsoEGTopBottomCenVert", dWord);
 
   l1Triggers[4] = l1SingleEG1 || l1SingleEG5 || l1SingleEG8 || l1SingleEG10 || l1SingleEG12 || l1SingleEG15 ||
                   l1SingleEG20 || l1SingleEG25 || l1DoubleNoIsoEGBTBtight || l1DoubleNoIsoEGBTBloose ||
                   l1DoubleNoIsoEGTopBottom || l1DoubleNoIsoEGTopBottomCen || l1DoubleNoIsoEGTopBottomCen2 ||
                   l1DoubleNoIsoEGTopBottomCenVert;
 
   std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
   std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr;
   for (igmtrr = gmt_records.begin(); igmtrr != gmt_records.end(); igmtrr++) {
     std::vector<L1MuRegionalCand>::const_iterator iter1;
     std::vector<L1MuRegionalCand> rmc;
 
     //DT triggers
     int idt = 0;
     rmc = igmtrr->getDTBXCands();
     for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
       if (!(*iter1).empty()) {
         idt++;
       }
     }
     //if(idt>0) std::cout << "Found " << idt << " valid DT candidates in bx wrt. L1A = "
     //  << igmtrr->getBxInEvent() << std::endl;
     if (igmtrr->getBxInEvent() == 0 && idt > 0)
       l1Triggers[0] = true;
 
     //RPC triggers
     int irpcb = 0;
     rmc = igmtrr->getBrlRPCCands();
     for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
       if (!(*iter1).empty()) {
         irpcb++;
       }
     }
     //if(irpcb>0) std::cout << "Found " << irpcb << " valid RPC candidates in bx wrt. L1A = "
     //  << igmtrr->getBxInEvent() << std::endl;
     if (igmtrr->getBxInEvent() == 0 && irpcb > 0)
       l1Triggers[2] = true;
 
     //CSC Triggers
     int icsc = 0;
     rmc = igmtrr->getCSCCands();
     for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
       if (!(*iter1).empty()) {
         icsc++;
       }
     }
     //if(icsc>0) std::cout << "Found " << icsc << " valid CSC candidates in bx wrt. L1A = "
     //  << igmtrr->getBxInEvent() << std::endl;
     if (igmtrr->getBxInEvent() == 0 && icsc > 0)
       l1Triggers[1] = true;
   }
 
   L1GlobalTriggerReadoutRecord const* gtrr = L1GTRR.product();
 
   for (int ibx = -1; ibx <= 1; ibx++) {
     bool hcal_top = false;
     bool hcal_bot = false;
     const L1GtPsbWord psb = gtrr->gtPsbWord(0xbb0d, ibx);
     std::vector<int> valid_phi;
     if ((psb.aData(4) & 0x3f) >= 1) {
       valid_phi.push_back((psb.aData(4) >> 10) & 0x1f);
     }
     if ((psb.bData(4) & 0x3f) >= 1) {
       valid_phi.push_back((psb.bData(4) >> 10) & 0x1f);
     }
     if ((psb.aData(5) & 0x3f) >= 1) {
       valid_phi.push_back((psb.aData(5) >> 10) & 0x1f);
     }
     if ((psb.bData(5) & 0x3f) >= 1) {
       valid_phi.push_back((psb.bData(5) >> 10) & 0x1f);
     }
     std::vector<int>::const_iterator iphi;
     for (iphi = valid_phi.begin(); iphi != valid_phi.end(); iphi++) {
       //std::cout << "Found HCAL mip with phi=" << *iphi << " in bx wrt. L1A = " << ibx << std::endl;
       if (*iphi < 9)
         hcal_top = true;
       if (*iphi > 8)
         hcal_bot = true;
     }
     if (ibx == 0 && hcal_top && hcal_bot)
       l1Triggers[3] = true;
   }
 
   edm::LogInfo("EcalCosmicsHists") << "**** Trigger SourceSource ****";
   if (l1Triggers[0])
     edm::LogInfo("EcalCosmicsHists") << "DT";
   if (l1Triggers[2])
     edm::LogInfo("EcalCosmicsHists") << "RPC";
   if (l1Triggers[1])
     edm::LogInfo("EcalCosmicsHists") << "CSC";
   if (l1Triggers[3])
     edm::LogInfo("EcalCosmicsHists") << "HCAL";
   if (l1Triggers[4])
     edm::LogInfo("EcalCosmicsHists") << "ECAL";
   edm::LogInfo("EcalCosmicsHists") << "************************";
 
   return l1Triggers;
 }

void EcalCosmicsHists::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 2758 of file EcalCosmicsHists.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

                               {
   using namespace std;
   if (runInFileName_) {
     fileName_ += "-" + intToString(runNum_) + ".graph.root";
   } else {
     fileName_ += ".root";
   }
 
   TFile root_file_(fileName_.c_str(), "RECREATE");
 
   for (map<int, TH1F*>::const_iterator itr = FEDsAndHists_.begin(); itr != FEDsAndHists_.end(); ++itr) {
     string dir = fedMap_->getSliceFromFed(itr->first);
     TDirectory* FEDdir = gDirectory->mkdir(dir.c_str());
     FEDdir->cd();
 
     TH1F* hist = itr->second;
     if (hist != nullptr)
       hist->Write();
     else {
       cerr << "EcalCosmicsHists: Error: This shouldn't happen!" << endl;
     }
     // Write out timing hist
     hist = FEDsAndTimingHists_[itr->first];
     if<

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation