#include <EcalSimple2007H4TBAnalyzer.h>

Inheritance diagram for EcalSimple2007H4TBAnalyzer:

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

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

	EcalSimple2007H4TBAnalyzer (const edm::ParameterSet &)

void	endJob () override

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

	~EcalSimple2007H4TBAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const std::string	digiCollection_

const std::string	digiProducer_

const edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_

const edm::EDGetTokenT< EEUncalibratedRecHitCollection >	eeUncalibratedRecHitToken_

const std::string	eventHeaderCollection_

const std::string	eventHeaderProducer_

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geometryToken_

TH2F *	h_ampltdc

TH1F *	h_bprofx

TH1F *	h_bprofy

TH1F *	h_e1e25

TH2F *	h_e1e25_mapx

TH2F *	h_e1e25_mapy

TH1F *	h_e1e9

TH2F *	h_e1e9_mapx

TH2F *	h_e1e9_mapy

TH1F *	h_e1x1

TH1F *	h_e1x1_center

TH1F *	h_e3x3

TH1F *	h_e3x3_center

TH1F *	h_e5x5

TH1F *	h_e5x5_center

TH1F *	h_e9e25

TH2F *	h_e9e25_mapx

TH2F *	h_e9e25_mapy

TH2F *	h_mapx [25]

TH2F *	h_mapy [25]

TH1F *	h_qualx

TH1F *	h_qualy

TH1F *	h_S6

TH2F *	h_Shape_

TH1F *	h_slopex

TH1F *	h_slopey

TH1F *	h_tableIsMoving

const std::string	hitCollection_

const std::string	hitProducer_

const std::string	hodoRecInfoCollection_

const std::string	hodoRecInfoProducer_

const std::string	rootfile_

const edm::EDGetTokenT< EcalTBEventHeader >	tbEventHeaderToken_

const edm::EDGetTokenT< EcalTBHodoscopeRecInfo >	tbHodoscopeRecInfoToken_

const edm::EDGetTokenT< EcalTBTDCRecInfo >	tbTDCRecInfoToken_

const std::string	tdcRecInfoCollection_

const std::string	tdcRecInfoProducer_

const CaloGeometry *	theTBGeometry_

EEDetId	xtalInBeam_

EBDetId	xtalInBeamTmp

EEDetId	Xtals5x5 [25]

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

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

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

Definition at line 43 of file EcalSimple2007H4TBAnalyzer.h.

Constructor & Destructor Documentation

◆ EcalSimple2007H4TBAnalyzer()

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

explicit

Definition at line 36 of file EcalSimple2007H4TBAnalyzer.cc.

References hitCollection_, and hitProducer_.

     : rootfile_(iConfig.getUntrackedParameter<std::string>("rootfile", "ecalSimpleTBanalysis.root")),
       digiCollection_(iConfig.getParameter<std::string>("digiCollection")),
       digiProducer_(iConfig.getParameter<std::string>("digiProducer")),
       hitCollection_(iConfig.getParameter<std::string>("hitCollection")),
       hitProducer_(iConfig.getParameter<std::string>("hitProducer")),
       hodoRecInfoCollection_(iConfig.getParameter<std::string>("hodoRecInfoCollection")),
       hodoRecInfoProducer_(iConfig.getParameter<std::string>("hodoRecInfoProducer")),
       tdcRecInfoCollection_(iConfig.getParameter<std::string>("tdcRecInfoCollection")),
       tdcRecInfoProducer_(iConfig.getParameter<std::string>("tdcRecInfoProducer")),
       eventHeaderCollection_(iConfig.getParameter<std::string>("eventHeaderCollection")),
       eventHeaderProducer_(iConfig.getParameter<std::string>("eventHeaderProducer")),
       eeDigiToken_(consumes<>(edm::InputTag(digiProducer_, digiCollection_))),
       eeUncalibratedRecHitToken_(consumes<>(edm::InputTag(hitProducer_, hitCollection_))),
       tbHodoscopeRecInfoToken_(consumes<>(edm::InputTag(hodoRecInfoProducer_, hodoRecInfoCollection_))),
       tbTDCRecInfoToken_(consumes<>(edm::InputTag(tdcRecInfoProducer_, tdcRecInfoCollection_))),
       tbEventHeaderToken_(consumes<>(edm::InputTag(eventHeaderProducer_))),
       geometryToken_(esConsumes<edm::Transition::BeginRun>()),
       xtalInBeam_(0)
 //========================================================================
 {
   //now do what ever initialization is needed
   edm::LogVerbatim("EcalSimple2007H4TBAnalyzer")
       << "EcalSimple2007H4TBAnalyzer: fetching hitCollection: " << hitCollection_.c_str() << " produced by "
       << hitProducer_.c_str();
 }

◆ ~EcalSimple2007H4TBAnalyzer()

EcalSimple2007H4TBAnalyzer::~EcalSimple2007H4TBAnalyzer ( )

override

Definition at line 64 of file EcalSimple2007H4TBAnalyzer.cc.

 {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   // Amplitude vs TDC offset
   //   if (h_ampltdc)
   //   delete h_ampltdc;
 
   //   // Reconstructed energies
   //   delete h_e1x1;
   //   delete h_e3x3;
   //   delete h_e5x5;
 
   //   delete h_bprofx;
   //   delete h_bprofy;
 
   //   delete h_qualx;
   //   delete h_qualy;
 
   //   delete h_slopex;
   //   delete h_slopey;
 
   //   delete h_mapx;
   //   delete h_mapy;
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 211 of file EcalSimple2007H4TBAnalyzer.cc.

References EcalMGPASample::adc(), CustomPhysics_cfi::amplitude, edm::DataFrameContainer::begin(), EcalTBEventHeader::crystalInBeam(), digiCollection_, eeDigiToken_, eeUncalibratedRecHitToken_, cosmicPhotonAnalyzer_cfi::eMax, edm::DataFrameContainer::end(), eventHeaderProducer_, EcalTBEventHeader::eventNumber(), CaloGeometry::getGeometry(), h_ampltdc, h_bprofx, h_bprofy, h_e1e25, h_e1e25_mapx, h_e1e25_mapy, h_e1e9, h_e1e9_mapx, h_e1e9_mapy, h_e1x1, h_e1x1_center, h_e3x3, h_e3x3_center, h_e5x5, h_e5x5_center, h_e9e25, h_e9e25_mapx, h_e9e25_mapy, h_mapx, h_mapy, h_qualx, h_qualy, h_S6, h_Shape_, h_slopex, h_slopey, h_tableIsMoving, hitCollection_, hfClusterShapes_cfi::hits, hodoRecInfoCollection_, mps_fire::i, EBDetId::ic(), EEDetId::ic(), iEvent, createfilelist::int, edm::HandleBase::isValid(), EEDetId::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), EEDetId::iy(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), MTVHistoProducerAlgoForTrackerBlock_cfi::maxHit, DetId::null(), EcalTBTDCRecInfo::offset(), phits, EcalTBHodoscopeRecInfo::posX(), EcalTBHodoscopeRecInfo::posY(), edm::Handle< T >::product(), EcalTBHodoscopeRecInfo::qualX(), EcalTBHodoscopeRecInfo::qualY(), EcalTBEventHeader::S6ADC(), EcalDataFrame::sample(), ecalGpuTask_cfi::sample, EcalDataFrame::size(), EcalTBHodoscopeRecInfo::slopeX(), EcalTBHodoscopeRecInfo::slopeY(), EBDetId::SMCRYSTALMODE, EcalTBEventHeader::tableIsMoving(), tbEventHeaderToken_, tbHodoscopeRecInfoToken_, tbTDCRecInfoToken_, tdcRecInfoCollection_, theTBGeometry_, x, xtalInBeam_, xtalInBeamTmp, Xtals5x5, y, and EEDetId::zside().

                                                                                             {
   //========================================================================
 
   using namespace edm;
   using namespace cms;
 
   Handle<EEDigiCollection> pdigis;
   const EEDigiCollection* digis = nullptr;
   iEvent.getByToken(eeDigiToken_, pdigis);
   if (pdigis.isValid()) {
     digis = pdigis.product();  // get a ptr to the product
   } else {
     edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << digiCollection_;
   }
 
   // fetch the digis and compute signal amplitude
   Handle<EEUncalibratedRecHitCollection> phits;
   const EEUncalibratedRecHitCollection* hits = nullptr;
   iEvent.getByToken(eeUncalibratedRecHitToken_, phits);
   if (phits.isValid()) {
     hits = phits.product();  // get a ptr to the product
   } else {
     edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << hitCollection_;
   }
 
   Handle<EcalTBHodoscopeRecInfo> pHodo;
   const EcalTBHodoscopeRecInfo* recHodo = nullptr;
   iEvent.getByToken(tbHodoscopeRecInfoToken_, pHodo);
   if (pHodo.isValid()) {
     recHodo = pHodo.product();  // get a ptr to the product
   } else {
     edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << hodoRecInfoCollection_;
   }
 
   Handle<EcalTBTDCRecInfo> pTDC;
   const EcalTBTDCRecInfo* recTDC = nullptr;
   iEvent.getByToken(tbTDCRecInfoToken_, pTDC);
   if (pTDC.isValid()) {
     recTDC = pTDC.product();  // get a ptr to the product
   } else {
     edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << tdcRecInfoCollection_;
   }
 
   Handle<EcalTBEventHeader> pEventHeader;
   const EcalTBEventHeader* evtHeader = nullptr;
   iEvent.getByToken(tbEventHeaderToken_, pEventHeader);
   if (pEventHeader.isValid()) {
     evtHeader = pEventHeader.product();  // get a ptr to the product
   } else {
     edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << eventHeaderProducer_;
   }
 
   if (!hits) {
     return;
   }
 
   if (!recTDC) {
     return;
   }
 
   if (!recHodo) {
     return;
   }
 
   if (!evtHeader) {
     return;
   }
 
   if (hits->empty()) {
     return;
   }
 
   //Accessing various event information
   if (evtHeader->tableIsMoving())
     h_tableIsMoving->Fill(evtHeader->eventNumber());
 
   //S6 beam scintillator
   h_S6->Fill(evtHeader->S6ADC());
 
   if (xtalInBeamTmp.null()) {
     xtalInBeamTmp = EBDetId(1, evtHeader->crystalInBeam(), EBDetId::SMCRYSTALMODE);
     xtalInBeam_ = EEDetId(35 - ((xtalInBeamTmp.ic() - 1) % 20), int(int(xtalInBeamTmp.ic()) / int(20)) + 51, -1);
     edm::LogVerbatim("EcalSimple2007H4TBAnalyzer") << "Xtal In Beam is " << xtalInBeam_.ic() << xtalInBeam_;
     for (unsigned int icry = 0; icry < 25; icry++) {
       unsigned int row = icry / 5;
       unsigned int column = icry % 5;
       int ix = xtalInBeam_.ix() + row - 2;
       int iy = xtalInBeam_.iy() + column - 2;
       EEDetId tempId(ix, iy, xtalInBeam_.zside());
       //Selecting matrix of xtals used in 2007H4TB
       if (tempId.ix() < 16 || tempId.ix() > 35 || tempId.iy() < 51 || tempId.iy() > 75)
         Xtals5x5[icry] = EEDetId(0);
       else {
         Xtals5x5[icry] = tempId;
         auto cell = theTBGeometry_->getGeometry(Xtals5x5[icry]);
         if (!cell)
           continue;
         edm::LogVerbatim("EcalSimple2007H4TBAnalyzer")
             << "** Xtal in the matrix **** row " << row << ", column " << column << ", xtal " << Xtals5x5[icry]
             << " Position " << cell->getPosition(0.);
       }
     }
   } else if (xtalInBeamTmp !=
              EBDetId(1, evtHeader->crystalInBeam(), EBDetId::SMCRYSTALMODE))  //run analysis only on first xtal in beam
     return;
 
   //Avoid moving table events
   if (evtHeader->tableIsMoving()) {
     edm::LogVerbatim("EcalSimple2007H4TBAnalyzer") << "Table is moving";
     return;
   }
 
   // Searching for max amplitude xtal alternative to use xtalInBeam_
   EEDetId maxHitId(0);
   float maxHit = -999999.;
   for (EEUncalibratedRecHitCollection::const_iterator ithit = hits->begin(); ithit != hits->end(); ++ithit) {
     if (ithit->amplitude() >= maxHit) {
       maxHit = ithit->amplitude();
       maxHitId = ithit->id();
     }
   }
   if (maxHitId == EEDetId(0)) {
     edm::LogVerbatim("EcalSimple2007H4TBAnalyzer") << "No maxHit found";
     return;
   }
 
   //Filling the digis shape for the xtalInBeam
   double samples_save[10];
   for (int i = 0; i < 10; ++i)
     samples_save[i] = 0.0;
 
   double eMax = 0.;
   for (EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
     if (EEDetId((*digiItr).id()) != xtalInBeam_)
       continue;
 
     EEDataFrame myDigi = (*digiItr);
     for (int sample = 0; sample < myDigi.size(); ++sample) {
       double analogSample = myDigi.sample(sample).adc();
       samples_save[sample] = analogSample;
       if (eMax < analogSample) {
         eMax = analogSample;
       }
     }
   }
 
   for (int i = 0; i < 10; ++i)
     h_Shape_->Fill(double(i) + recTDC->offset(), samples_save[i]);
 
   // Taking amplitudes in 5x5
   double amplitude[25];
   double amplitude3x3 = 0;
   double amplitude5x5 = 0;
   for (unsigned int icry = 0; icry < 25; icry++) {
     if (!Xtals5x5[icry].null()) {
       amplitude[icry] = (hits->find(Xtals5x5[icry]))->amplitude();
       amplitude5x5 += amplitude[icry];
       // Is in 3x3?
       if (icry == 6 || icry == 7 || icry == 8 || icry == 11 || icry == 12 || icry == 13 || icry == 16 || icry == 17 ||
           icry == 18) {
         amplitude3x3 += amplitude[icry];
       }
     }
   }
 
   //Filling amplitudes
   h_e1x1->Fill(amplitude[12]);
   h_e3x3->Fill(amplitude3x3);
   h_e5x5->Fill(amplitude5x5);
 
   h_e1e9->Fill(amplitude[12] / amplitude3x3);
   h_e1e25->Fill(amplitude[12] / amplitude5x5);
   h_e9e25->Fill(amplitude3x3 / amplitude5x5);
 
   //Checking stability of amplitude vs TDC
   if (recTDC)
     h_ampltdc->Fill(recTDC->offset(), amplitude[12]);
 
   //Various amplitudes as a function of hodoscope coordinates
   if (recHodo) {
     float x = recHodo->posX();
     float y = recHodo->posY();
     float xslope = recHodo->slopeX();
     float yslope = recHodo->slopeY();
     float xqual = recHodo->qualX();
     float yqual = recHodo->qualY();
 
     //Filling beam profiles
     h_bprofx->Fill(x);
     h_bprofy->Fill(y);
     h_qualx->Fill(xqual);
     h_qualy->Fill(yqual);
     h_slopex->Fill(xslope);
     h_slopey->Fill(yslope);
 
     //Fill central events
 
     if (fabs(x + 2.5) < 2.5 && fabs(y + 0.5) < 2.5) {
       h_e1x1_center->Fill(amplitude[12]);
       h_e3x3_center->Fill(amplitude3x3);
       h_e5x5_center->Fill(amplitude5x5);
     }
 
     for (unsigned int icry = 0; icry < 25; icry++) {
       h_mapx[icry]->Fill(x, amplitude[icry]);
       h_mapy[icry]->Fill(y, amplitude[icry]);
     }
 
     h_e1e9_mapx->Fill(x, amplitude[12] / amplitude3x3);
     h_e1e9_mapy->Fill(y, amplitude[12] / amplitude3x3);
 
     h_e1e25_mapx->Fill(x, amplitude[12] / amplitude5x5);
     h_e1e25_mapy->Fill(y, amplitude[12] / amplitude5x5);
 
     h_e9e25_mapx->Fill(x, amplitude3x3 / amplitude5x5);
     h_e9e25_mapy->Fill(y, amplitude3x3 / amplitude5x5);
   }
 }

◆ beginRun()

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

override

Definition at line 92 of file EcalSimple2007H4TBAnalyzer.cc.

References MillePedeFileConverter_cfg::e, geometryToken_, edm::EventSetup::getData(), h_ampltdc, h_bprofx, h_bprofy, h_e1e25, h_e1e25_mapx, h_e1e25_mapy, h_e1e9, h_e1e9_mapx, h_e1e9_mapy, h_e1x1, h_e1x1_center, h_e3x3, h_e3x3_center, h_e5x5, h_e5x5_center, h_e9e25, h_e9e25_mapx, h_e9e25_mapy, h_mapx, h_mapy, h_qualx, h_qualy, h_S6, h_Shape_, h_slopex, h_slopey, h_tableIsMoving, and theTBGeometry_.

                                                                                     {
   //========================================================================
 
   theTBGeometry_ = &iSetup.getData(geometryToken_);
 
   // Amplitude vs TDC offset
   h_ampltdc = new TH2F("h_ampltdc", "Max Amplitude vs TDC offset", 100, 0., 1., 1000, 0., 4000.);
 
   // Reconstructed energies
   h_tableIsMoving = new TH1F("h_tableIsMoving", "TableIsMoving", 100000, 0., 100000.);
 
   h_e1x1 = new TH1F("h_e1x1", "E1x1 energy", 1000, 0., 4000.);
   h_e3x3 = new TH1F("h_e3x3", "E3x3 energy", 1000, 0., 4000.);
   h_e5x5 = new TH1F("h_e5x5", "E5x5 energy", 1000, 0., 4000.);
 
   h_e1x1_center = new TH1F("h_e1x1_center", "E1x1 energy", 1000, 0., 4000.);
   h_e3x3_center = new TH1F("h_e3x3_center", "E3x3 energy", 1000, 0., 4000.);
   h_e5x5_center = new TH1F("h_e5x5_center", "E5x5 energy", 1000, 0., 4000.);
 
   h_e1e9 = new TH1F("h_e1e9", "E1/E9 ratio", 600, 0., 1.2);
   h_e1e25 = new TH1F("h_e1e25", "E1/E25 ratio", 600, 0., 1.2);
   h_e9e25 = new TH1F("h_e9e25", "E9/E25 ratio", 600, 0., 1.2);
 
   h_S6 = new TH1F("h_S6", "Amplitude S6", 1000, 0., 4000.);
 
   h_bprofx = new TH1F("h_bprofx", "Beam Profile X", 100, -20., 20.);
   h_bprofy = new TH1F("h_bprofy", "Beam Profile Y", 100, -20., 20.);
 
   h_qualx = new TH1F("h_qualx", "Beam Quality X", 5000, 0., 5.);
   h_qualy = new TH1F("h_qualy", "Beam Quality X", 5000, 0., 5.);
 
   h_slopex = new TH1F("h_slopex", "Beam Slope X", 500, -5e-4, 5e-4);
   h_slopey = new TH1F("h_slopey", "Beam Slope Y", 500, -5e-4, 5e-4);
 
   char hname[50];
   char htitle[50];
   for (unsigned int icry = 0; icry < 25; icry++) {
     sprintf(hname, "h_mapx_%d", icry);
     sprintf(htitle, "Max Amplitude vs X %d", icry);
     h_mapx[icry] = new TH2F(hname, htitle, 80, -20, 20, 1000, 0., 4000.);
     sprintf(hname, "h_mapy_%d", icry);
     sprintf(htitle, "Max Amplitude vs Y %d", icry);
     h_mapy[icry] = new TH2F(hname, htitle, 80, -20, 20, 1000, 0., 4000.);
   }
 
   h_e1e9_mapx = new TH2F("h_e1e9_mapx", "E1/E9 vs X", 80, -20, 20, 600, 0., 1.2);
   h_e1e9_mapy = new TH2F("h_e1e9_mapy", "E1/E9 vs Y", 80, -20, 20, 600, 0., 1.2);
 
   h_e1e25_mapx = new TH2F("h_e1e25_mapx", "E1/E25 vs X", 80, -20, 20, 600, 0., 1.2);
   h_e1e25_mapy = new TH2F("h_e1e25_mapy", "E1/E25 vs Y", 80, -20, 20, 600, 0., 1.2);
 
   h_e9e25_mapx = new TH2F("h_e9e25_mapx", "E9/E25 vs X", 80, -20, 20, 600, 0., 1.2);
   h_e9e25_mapy = new TH2F("h_e9e25_mapy", "E9/E25 vs Y", 80, -20, 20, 600, 0., 1.2);
 
   h_Shape_ = new TH2F("h_Shape_", "Xtal in Beam Shape", 250, 0, 10, 350, 0, 3500);
 }

◆ endJob()

void EcalSimple2007H4TBAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 154 of file EcalSimple2007H4TBAnalyzer.cc.

References f, h_ampltdc, h_bprofx, h_bprofy, h_e1e25, h_e1e25_mapx, h_e1e25_mapy, h_e1e9, h_e1e9_mapx, h_e1e9_mapy, h_e1x1, h_e1x1_center, h_e3x3, h_e3x3_center, h_e5x5, h_e5x5_center, h_e9e25, h_e9e25_mapx, h_e9e25_mapy, h_mapx, h_mapy, h_qualx, h_qualy, h_S6, h_Shape_, h_slopex, h_slopey, h_tableIsMoving, and rootfile_.

                                         {
   //========================================================================
 
   TFile f(rootfile_.c_str(), "RECREATE");
 
   // Amplitude vs TDC offset
   h_ampltdc->Write();
 
   // Reconstructed energies
   h_e1x1->Write();
   h_e3x3->Write();
   h_e5x5->Write();
 
   h_e1x1_center->Write();
   h_e3x3_center->Write();
   h_e5x5_center->Write();
 
   h_e1e9->Write();
   h_e1e25->Write();
   h_e9e25->Write();
 
   h_S6->Write();
   h_bprofx->Write();
   h_bprofy->Write();
 
   h_qualx->Write();
   h_qualy->Write();
 
   h_slopex->Write();
   h_slopey->Write();
 
   h_Shape_->Write();
 
   for (unsigned int icry = 0; icry < 25; icry++) {
     h_mapx[icry]->Write();
     h_mapy[icry]->Write();
   }
 
   h_e1e9_mapx->Write();
   h_e1e9_mapy->Write();
 
   h_e1e25_mapx->Write();
   h_e1e25_mapy->Write();
 
   h_e9e25_mapx->Write();
   h_e9e25_mapy->Write();
 
   h_tableIsMoving->Write();
 
   f.Close();
 }

◆ endRun()

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

override

Definition at line 150 of file EcalSimple2007H4TBAnalyzer.cc.

150 {}

Member Data Documentation

◆ digiCollection_

const std::string EcalSimple2007H4TBAnalyzer::digiCollection_

private

Definition at line 55 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ digiProducer_

const std::string EcalSimple2007H4TBAnalyzer::digiProducer_

private

Definition at line 56 of file EcalSimple2007H4TBAnalyzer.h.

◆ eeDigiToken_

const edm::EDGetTokenT<EEDigiCollection> EcalSimple2007H4TBAnalyzer::eeDigiToken_

private

Definition at line 66 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ eeUncalibratedRecHitToken_

const edm::EDGetTokenT<EEUncalibratedRecHitCollection> EcalSimple2007H4TBAnalyzer::eeUncalibratedRecHitToken_

private

Definition at line 67 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ eventHeaderCollection_

const std::string EcalSimple2007H4TBAnalyzer::eventHeaderCollection_

private

Definition at line 63 of file EcalSimple2007H4TBAnalyzer.h.

◆ eventHeaderProducer_

const std::string EcalSimple2007H4TBAnalyzer::eventHeaderProducer_

private

Definition at line 64 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ geometryToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> EcalSimple2007H4TBAnalyzer::geometryToken_

private

Definition at line 71 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by beginRun().

◆ h_ampltdc

TH2F* EcalSimple2007H4TBAnalyzer::h_ampltdc

private

Definition at line 74 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_bprofx

TH1F* EcalSimple2007H4TBAnalyzer::h_bprofx

private

Definition at line 93 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_bprofy

TH1F* EcalSimple2007H4TBAnalyzer::h_bprofy

private

Definition at line 94 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e25

TH1F* EcalSimple2007H4TBAnalyzer::h_e1e25

private

Definition at line 89 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e25_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e25_mapx

private

Definition at line 108 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e25_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e25_mapy

private

Definition at line 109 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e9

TH1F* EcalSimple2007H4TBAnalyzer::h_e1e9

private

Definition at line 88 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e9_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e9_mapx

private

Definition at line 105 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e9_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e9_mapy

private

Definition at line 106 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1x1

TH1F* EcalSimple2007H4TBAnalyzer::h_e1x1

private

Definition at line 80 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1x1_center

TH1F* EcalSimple2007H4TBAnalyzer::h_e1x1_center

private

Definition at line 84 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e3x3

TH1F* EcalSimple2007H4TBAnalyzer::h_e3x3

private

Definition at line 81 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e3x3_center

TH1F* EcalSimple2007H4TBAnalyzer::h_e3x3_center

private

Definition at line 85 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e5x5

TH1F* EcalSimple2007H4TBAnalyzer::h_e5x5

private

Definition at line 82 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e5x5_center

TH1F* EcalSimple2007H4TBAnalyzer::h_e5x5_center

private

Definition at line 86 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e9e25

TH1F* EcalSimple2007H4TBAnalyzer::h_e9e25

private

Definition at line 90 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e9e25_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_e9e25_mapx

private

Definition at line 111 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e9e25_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_e9e25_mapy

private

Definition at line 112 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_mapx[25]

private

Definition at line 102 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_mapy[25]

private

Definition at line 103 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_qualx

TH1F* EcalSimple2007H4TBAnalyzer::h_qualx

private

Definition at line 96 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_qualy

TH1F* EcalSimple2007H4TBAnalyzer::h_qualy

private

Definition at line 97 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_S6

TH1F* EcalSimple2007H4TBAnalyzer::h_S6

private

Definition at line 92 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_Shape_

TH2F* EcalSimple2007H4TBAnalyzer::h_Shape_

private

Definition at line 76 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_slopex

TH1F* EcalSimple2007H4TBAnalyzer::h_slopex

private

Definition at line 99 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_slopey

TH1F* EcalSimple2007H4TBAnalyzer::h_slopey

private

Definition at line 100 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_tableIsMoving

TH1F* EcalSimple2007H4TBAnalyzer::h_tableIsMoving

private

Definition at line 79 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ hitCollection_

const std::string EcalSimple2007H4TBAnalyzer::hitCollection_

private

Definition at line 57 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ hitProducer_

const std::string EcalSimple2007H4TBAnalyzer::hitProducer_

private

Definition at line 58 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by EcalSimple2007H4TBAnalyzer().

◆ hodoRecInfoCollection_

const std::string EcalSimple2007H4TBAnalyzer::hodoRecInfoCollection_

private

Definition at line 59 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ hodoRecInfoProducer_

const std::string EcalSimple2007H4TBAnalyzer::hodoRecInfoProducer_

private

Definition at line 60 of file EcalSimple2007H4TBAnalyzer.h.

◆ rootfile_

const std::string EcalSimple2007H4TBAnalyzer::rootfile_

private

Definition at line 54 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by endJob().

◆ tbEventHeaderToken_

const edm::EDGetTokenT<EcalTBEventHeader> EcalSimple2007H4TBAnalyzer::tbEventHeaderToken_

private

Definition at line 70 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ tbHodoscopeRecInfoToken_

const edm::EDGetTokenT<EcalTBHodoscopeRecInfo> EcalSimple2007H4TBAnalyzer::tbHodoscopeRecInfoToken_

private

Definition at line 68 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ tbTDCRecInfoToken_

const edm::EDGetTokenT<EcalTBTDCRecInfo> EcalSimple2007H4TBAnalyzer::tbTDCRecInfoToken_

private

Definition at line 69 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ tdcRecInfoCollection_

const std::string EcalSimple2007H4TBAnalyzer::tdcRecInfoCollection_

private

Definition at line 61 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ tdcRecInfoProducer_

const std::string EcalSimple2007H4TBAnalyzer::tdcRecInfoProducer_

private

Definition at line 62 of file EcalSimple2007H4TBAnalyzer.h.

◆ theTBGeometry_

const CaloGeometry* EcalSimple2007H4TBAnalyzer::theTBGeometry_

private

Definition at line 118 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and beginRun().

◆ xtalInBeam_

EEDetId EcalSimple2007H4TBAnalyzer::xtalInBeam_

private

Definition at line 114 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ xtalInBeamTmp

EBDetId EcalSimple2007H4TBAnalyzer::xtalInBeamTmp

private

Definition at line 115 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ Xtals5x5

EEDetId EcalSimple2007H4TBAnalyzer::Xtals5x5[25]

private

Definition at line 116 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalSimple2007H4TBAnalyzer()

◆ ~EcalSimple2007H4TBAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ endRun()

Member Data Documentation

◆ digiCollection_

◆ digiProducer_

◆ eeDigiToken_

◆ eeUncalibratedRecHitToken_

◆ eventHeaderCollection_

◆ eventHeaderProducer_

◆ geometryToken_

◆ h_ampltdc

◆ h_bprofx

◆ h_bprofy

◆ h_e1e25

◆ h_e1e25_mapx

◆ h_e1e25_mapy

◆ h_e1e9

◆ h_e1e9_mapx

◆ h_e1e9_mapy

◆ h_e1x1

◆ h_e1x1_center

◆ h_e3x3

◆ h_e3x3_center

◆ h_e5x5

◆ h_e5x5_center

◆ h_e9e25

◆ h_e9e25_mapx

◆ h_e9e25_mapy

◆ h_mapx

◆ h_mapy

◆ h_qualx

◆ h_qualy

◆ h_S6

◆ h_Shape_

◆ h_slopex

◆ h_slopey

◆ h_tableIsMoving

◆ hitCollection_

◆ hitProducer_

◆ hodoRecInfoCollection_

◆ hodoRecInfoProducer_

◆ rootfile_

◆ tbEventHeaderToken_

◆ tbHodoscopeRecInfoToken_

◆ tbTDCRecInfoToken_

◆ tdcRecInfoCollection_

◆ tdcRecInfoProducer_

◆ theTBGeometry_

◆ xtalInBeam_

◆ xtalInBeamTmp

◆ Xtals5x5