#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

	~EcalSimple2007H4TBAnalyzer () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > 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 &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
std::string	digiCollection_

std::string	digiProducer_

std::string	eventHeaderCollection_

std::string	eventHeaderProducer_

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

std::string	hitCollection_

std::string	hitProducer_

std::string	hodoRecInfoCollection_

std::string	hodoRecInfoProducer_

std::string	rootfile_

std::string	tdcRecInfoCollection_

std::string	tdcRecInfoProducer_

const CaloGeometry *	theTBGeometry_

EEDetId	xtalInBeam_

EBDetId	xtalInBeamTmp

EEDetId	Xtals5x5 [25]

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	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

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

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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

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

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

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

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

Detailed Description

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

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

Definition at line 40 of file EcalSimple2007H4TBAnalyzer.h.

Constructor & Destructor Documentation

◆ EcalSimple2007H4TBAnalyzer()

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

explicit

Definition at line 49 of file EcalSimple2007H4TBAnalyzer.cc.

     : xtalInBeam_(0)
 //========================================================================
 {
   //now do what ever initialization is needed
   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");
  
   std::cout << "EcalSimple2007H4TBAnalyzer: fetching hitCollection: " << hitCollection_.c_str() << " produced by "
             << hitProducer_.c_str() << std::endl;
 }

References gather_cfg::cout, digiCollection_, digiProducer_, eventHeaderCollection_, eventHeaderProducer_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hitCollection_, hitProducer_, hodoRecInfoCollection_, hodoRecInfoProducer_, rootfile_, AlCaHLTBitMon_QueryRunRegistry::string, tdcRecInfoCollection_, and tdcRecInfoProducer_.

◆ ~EcalSimple2007H4TBAnalyzer()

EcalSimple2007H4TBAnalyzer::~EcalSimple2007H4TBAnalyzer ( )

override

Definition at line 71 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::EDAnalyzer.

Definition at line 221 of file EcalSimple2007H4TBAnalyzer.cc.

                                                                                             {
   //========================================================================
  
   using namespace edm;
   using namespace cms;
  
   Handle<EEDigiCollection> pdigis;
   const EEDigiCollection* digis = nullptr;
   //std::cout << "EcalSimple2007H4TBAnalyzer::analyze getting product with label: " << digiProducer_.c_str()<< " prodname: " << digiCollection_.c_str() << endl;
   iEvent.getByLabel(digiProducer_, digiCollection_, pdigis);
   if (pdigis.isValid()) {
     digis = pdigis.product();  // get a ptr to the product
     //iEvent.getByLabel( hitProducer_, phits);
   } 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;
   //std::cout << "EcalSimple2007H4TBAnalyzer::analyze getting product with label: " << digiProducer_.c_str()<< " prodname: " << digiCollection_.c_str() << endl;
   iEvent.getByLabel(hitProducer_, hitCollection_, phits);
   if (phits.isValid()) {
     hits = phits.product();  // get a ptr to the product
     //iEvent.getByLabel( hitProducer_, phits);
   } else {
     edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << hitCollection_;
   }
  
   Handle<EcalTBHodoscopeRecInfo> pHodo;
   const EcalTBHodoscopeRecInfo* recHodo = nullptr;
   //std::cout << "EcalSimple2007H4TBAnalyzer::analyze getting product with label: " << digiProducer_.c_str()<< " prodname: " << digiCollection_.c_str() << endl;
   iEvent.getByLabel(hodoRecInfoProducer_, hodoRecInfoCollection_, 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;
   //std::cout << "EcalSimple2007H4TBAnalyzer::analyze getting product with label: " << digiProducer_.c_str()<< " prodname: " << digiCollection_.c_str() << endl;
   iEvent.getByLabel(tdcRecInfoProducer_, tdcRecInfoCollection_, 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;
   //std::cout << "EcalSimple2007H4TBAnalyzer::analyze getting product with label: " << digiProducer_.c_str()<< " prodname: " << digiCollection_.c_str() << endl;
   iEvent.getByLabel(eventHeaderProducer_, 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) {
     //       std::cout << "1" << std::endl;
     return;
   }
  
   if (!recTDC) {
     //       std::cout << "2" << std::endl;
     return;
   }
  
   if (!recHodo) {
     //       std::cout << "3" << std::endl;
     return;
   }
  
   if (!evtHeader) {
     //       std::cout << "4" << std::endl;
     return;
   }
  
   if (hits->empty()) {
     //       std::cout << "5" << std::endl;
     return;
   }
  
   //Accessing various event information
   if (evtHeader->tableIsMoving())
     h_tableIsMoving->Fill(evtHeader->eventNumber());
  
   //    std::cout << "event " << evtHeader->eventNumber()
   //         << "\trun number " << evtHeader->runNumber()
   //         << "\tburst number " << evtHeader->burstNumber()
   //         << "\tbeginLV1A " << evtHeader->begBurstLV1A()
   //         << "\tendLV1A " << evtHeader->endBurstLV1A()
   //         << "\ttime " << evtHeader->date()
   //         << "\thas errors " << int(evtHeader->syncError())
   //         << std::endl;
  
   //    std::cout << "scaler";
   //    for (int iscaler=0;iscaler < 36;iscaler++)
   //      std::cout << "\t#" << iscaler << " " <<  evtHeader->scaler(iscaler);
   //    std::cout<<std::endl;
  
   //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);
     std::cout << "Xtal In Beam is " << xtalInBeam_.ic() << xtalInBeam_ << std::endl;
     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;
         std::cout << "** Xtal in the matrix **** row " << row << ", column " << column << ", xtal " << Xtals5x5[icry]
                   << " Position " << cell->getPosition(0.) << std::endl;
       }
     }
   } 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()) {
     std::cout << "Table is moving" << std::endl;
     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)) {
     std::cout << "No maxHit found" << std::endl;
     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;
       //  std::cout << analogSample << " ";
       if (eMax < analogSample) {
         eMax = analogSample;
       }
     }
     // std::cout << std::endl;
   }
  
   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);
   }
 }

References EcalMGPASample::adc(), CustomPhysics_cfi::amplitude, edm::DataFrameContainer::begin(), gather_cfg::cout, EcalTBEventHeader::crystalInBeam(), digiCollection_, digiProducer_, 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_, hitProducer_, hfClusterShapes_cfi::hits, hodoRecInfoCollection_, hodoRecInfoProducer_, mps_fire::i, EBDetId::ic(), EEDetId::ic(), iEvent, createfilelist::int, edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), MTVHistoProducerAlgoForTrackerBlock_cfi::maxHit, DetId::null(), EcalTBTDCRecInfo::offset(), phits, EcalTBHodoscopeRecInfo::posX(), EcalTBHodoscopeRecInfo::posY(), edm::Handle< T >::product(), EcalTBHodoscopeRecInfo::qualX(), EcalTBHodoscopeRecInfo::qualY(), EcalTBEventHeader::S6ADC(), simplePhotonAnalyzer_cfi::sample, EcalDataFrame::sample(), EcalDataFrame::size(), EcalTBHodoscopeRecInfo::slopeX(), EcalTBHodoscopeRecInfo::slopeY(), EBDetId::SMCRYSTALMODE, EcalTBEventHeader::tableIsMoving(), tdcRecInfoCollection_, tdcRecInfoProducer_, theTBGeometry_, x, xtalInBeam_, xtalInBeamTmp, Xtals5x5, y, and EEDetId::zside().

◆ beginRun()

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

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 99 of file EcalSimple2007H4TBAnalyzer.cc.

                                                                                     {
   //========================================================================
  
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
  
   theTBGeometry_ = pG.product();
   //  const std::vector<DetId>& validIds=theTBGeometry_->getValidDetIds(DetId::Ecal,EcalEndcap);
   //   std::cout << "Found " << validIds.size() << " channels in the geometry" << std::endl;
   //   for (unsigned int i=0;i<validIds.size();++i)
   //     std::cout << EEDetId(validIds[i]) << std::endl;
  
   // 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);
 }

References MillePedeFileConverter_cfg::e, edm::EventSetup::get(), get, 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, edm::ESHandle< T >::product(), and theTBGeometry_.

◆ endJob()

void EcalSimple2007H4TBAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 164 of file EcalSimple2007H4TBAnalyzer.cc.

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

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_.

Member Data Documentation

◆ digiCollection_

std::string EcalSimple2007H4TBAnalyzer::digiCollection_

private

Definition at line 51 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ digiProducer_

std::string EcalSimple2007H4TBAnalyzer::digiProducer_

private

Definition at line 52 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ eventHeaderCollection_

std::string EcalSimple2007H4TBAnalyzer::eventHeaderCollection_

private

Definition at line 59 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by EcalSimple2007H4TBAnalyzer().

◆ eventHeaderProducer_

std::string EcalSimple2007H4TBAnalyzer::eventHeaderProducer_

private

Definition at line 60 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ h_ampltdc

TH2F* EcalSimple2007H4TBAnalyzer::h_ampltdc

private

Definition at line 63 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_bprofx

TH1F* EcalSimple2007H4TBAnalyzer::h_bprofx

private

Definition at line 82 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_bprofy

TH1F* EcalSimple2007H4TBAnalyzer::h_bprofy

private

Definition at line 83 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e25

TH1F* EcalSimple2007H4TBAnalyzer::h_e1e25

private

Definition at line 78 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e25_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e25_mapx

private

Definition at line 97 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e25_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e25_mapy

private

Definition at line 98 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e9

TH1F* EcalSimple2007H4TBAnalyzer::h_e1e9

private

Definition at line 77 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e9_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e9_mapx

private

Definition at line 94 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1e9_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e9_mapy

private

Definition at line 95 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1x1

TH1F* EcalSimple2007H4TBAnalyzer::h_e1x1

private

Definition at line 69 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e1x1_center

TH1F* EcalSimple2007H4TBAnalyzer::h_e1x1_center

private

Definition at line 73 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e3x3

TH1F* EcalSimple2007H4TBAnalyzer::h_e3x3

private

Definition at line 70 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e3x3_center

TH1F* EcalSimple2007H4TBAnalyzer::h_e3x3_center

private

Definition at line 74 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e5x5

TH1F* EcalSimple2007H4TBAnalyzer::h_e5x5

private

Definition at line 71 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e5x5_center

TH1F* EcalSimple2007H4TBAnalyzer::h_e5x5_center

private

Definition at line 75 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e9e25

TH1F* EcalSimple2007H4TBAnalyzer::h_e9e25

private

Definition at line 79 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e9e25_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_e9e25_mapx

private

Definition at line 100 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_e9e25_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_e9e25_mapy

private

Definition at line 101 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_mapx

TH2F* EcalSimple2007H4TBAnalyzer::h_mapx[25]

private

Definition at line 91 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_mapy

TH2F* EcalSimple2007H4TBAnalyzer::h_mapy[25]

private

Definition at line 92 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_qualx

TH1F* EcalSimple2007H4TBAnalyzer::h_qualx

private

Definition at line 85 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_qualy

TH1F* EcalSimple2007H4TBAnalyzer::h_qualy

private

Definition at line 86 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_S6

TH1F* EcalSimple2007H4TBAnalyzer::h_S6

private

Definition at line 81 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_Shape_

TH2F* EcalSimple2007H4TBAnalyzer::h_Shape_

private

Definition at line 65 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_slopex

TH1F* EcalSimple2007H4TBAnalyzer::h_slopex

private

Definition at line 88 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_slopey

TH1F* EcalSimple2007H4TBAnalyzer::h_slopey

private

Definition at line 89 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ h_tableIsMoving

TH1F* EcalSimple2007H4TBAnalyzer::h_tableIsMoving

private

Definition at line 68 of file EcalSimple2007H4TBAnalyzer.h.

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

◆ hitCollection_

std::string EcalSimple2007H4TBAnalyzer::hitCollection_

private

Definition at line 53 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ hitProducer_

std::string EcalSimple2007H4TBAnalyzer::hitProducer_

private

Definition at line 54 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ hodoRecInfoCollection_

std::string EcalSimple2007H4TBAnalyzer::hodoRecInfoCollection_

private

Definition at line 55 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ hodoRecInfoProducer_

std::string EcalSimple2007H4TBAnalyzer::hodoRecInfoProducer_

private

Definition at line 56 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ rootfile_

std::string EcalSimple2007H4TBAnalyzer::rootfile_

private

Definition at line 50 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by EcalSimple2007H4TBAnalyzer(), and endJob().

◆ tdcRecInfoCollection_

std::string EcalSimple2007H4TBAnalyzer::tdcRecInfoCollection_

private

Definition at line 57 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ tdcRecInfoProducer_

std::string EcalSimple2007H4TBAnalyzer::tdcRecInfoProducer_

private

Definition at line 58 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

◆ theTBGeometry_

const CaloGeometry* EcalSimple2007H4TBAnalyzer::theTBGeometry_

private

Definition at line 107 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and beginRun().

◆ xtalInBeam_

EEDetId EcalSimple2007H4TBAnalyzer::xtalInBeam_

private

Definition at line 103 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ xtalInBeamTmp

EBDetId EcalSimple2007H4TBAnalyzer::xtalInBeamTmp

private

Definition at line 104 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

◆ Xtals5x5

EEDetId EcalSimple2007H4TBAnalyzer::Xtals5x5[25]

private

Definition at line 105 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()

Member Data Documentation

◆ digiCollection_

◆ digiProducer_

◆ eventHeaderCollection_

◆ eventHeaderProducer_

◆ 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_

◆ tdcRecInfoCollection_

◆ tdcRecInfoProducer_

◆ theTBGeometry_

◆ xtalInBeam_

◆ xtalInBeamTmp

◆ Xtals5x5