#include <PhiSymmetryCalibration.h>

Inheritance diagram for PhiSymmetryCalibration:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Called at each event. More...

void	beginJob () override
	Called at beginning of job. More...

void	endJob () override
	Called at end of job. More...

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

	PhiSymmetryCalibration (const edm::ParameterSet &iConfig)
	Constructor. More...

void	setUp (const edm::EventSetup &setup)

	~PhiSymmetryCalibration () override
	Destructor. More...

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 Member Functions
void	getKfactors ()

Private Attributes
double	ap_

double	b_

std::vector< DetId >	barrelCells

std::string	barrelHits_

EcalGeomPhiSymHelper	e_

std::vector< TH1F * >	e_spectrum_b_histos

std::vector< TH1F * >	e_spectrum_e_histos

std::string	ecalHitsProducer_

double	eCut_barl_

std::vector< DetId >	endcapCells

std::string	endcapHits_

float	epsilon_M_barl [kBarlRings][kBarlWedges][kSides]

float	epsilon_M_endc [kEndcWedgesX][kEndcWedgesY][kSides]

double	esumMean_barl_ [kBarlRings]

double	esumMean_endc_ [kEndcEtaRings]

std::vector< TH1F * >	et_spectrum_b_histos

std::vector< TH1F * >	et_spectrum_e_histos

double	etsum_barl_ [kBarlRings][kBarlWedges][kSides]

double	etsum_barl_miscal_ [kNMiscalBinsEB][kBarlRings]

double	etsum_endc_ [kEndcWedgesX][kEndcWedgesX][kSides]

double	etsum_endc_miscal_ [kNMiscalBinsEE][kEndcEtaRings]

double	etsum_endc_uncorr [kEndcWedgesX][kEndcWedgesX][kSides]

double	etsumMean_barl_ [kBarlRings]

double	etsumMean_endc_ [kEndcEtaRings]

int	eventSet_

int	eventsinlb_

int	eventsinrun_

bool	isfirstpass_

double	k_barl_ [kBarlRings]

double	k_endc_ [kEndcEtaRings]

double	miscalEB_ [kNMiscalBinsEB]

double	miscalEE_ [kNMiscalBinsEE]

int	nevents_

double	newCalibs_barl [kBarlRings][kBarlWedges][kSides]

double	newCalibs_endc [kEndcWedgesX][kEndcWedgesX][kSides]

unsigned int	nhits_barl_ [kBarlRings][kBarlWedges][kSides]

unsigned int	nhits_endc_ [kEndcWedgesX][kEndcWedgesX][kSides]

std::string	oldcalibfile_

EcalIntercalibConstants	oldCalibs_
	the old calibration constants (when reiterating, the last ones derived) More...

double	oldCalibs_barl [kBarlRings][kBarlWedges][kSides]

double	oldCalibs_endc [kEndcWedgesX][kEndcWedgesY][kSides]

float	phi_endc [kMaxEndciPhi][kEndcEtaRings]

float	rawconst_barl [kBarlRings][kBarlWedges][kSides]

float	rawconst_endc [kEndcWedgesX][kEndcWedgesX][kSides]

bool	reiteration_

bool	spectra

int	statusThreshold_
	threshold in channel status beyond which channel is marked bad More...

Static Private Attributes
static const int	kMaxEndciPhi = 360

static const float	kMiscalRangeEB = .05

static const float	kMiscalRangeEE = .10

static const int	kNMiscalBinsEB = 21

static const int	kNMiscalBinsEE = 41

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

Definition at line 34 of file PhiSymmetryCalibration.h.

Constructor & Destructor Documentation

◆ PhiSymmetryCalibration()

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

Constructor.

Definition at line 49 of file PhiSymmetryCalibration.cc.

     :
  
       ecalHitsProducer_(iConfig.getParameter<std::string>("ecalRecHitsProducer")),
       barrelHits_(iConfig.getParameter<std::string>("barrelHitCollection")),
       endcapHits_(iConfig.getParameter<std::string>("endcapHitCollection")),
       eCut_barl_(iConfig.getParameter<double>("eCut_barrel")),
       ap_(iConfig.getParameter<double>("ap")),
       b_(iConfig.getParameter<double>("b")),
       eventSet_(iConfig.getParameter<int>("eventSet")),
       statusThreshold_(iConfig.getUntrackedParameter<int>("statusThreshold", 3)),
       reiteration_(iConfig.getUntrackedParameter<bool>("reiteration", false)),
       oldcalibfile_(iConfig.getUntrackedParameter<std::string>("oldcalibfile", "EcalintercalibConstants.xml")) {
   isfirstpass_ = true;
  
   et_spectrum_b_histos.resize(kBarlRings);
   e_spectrum_b_histos.resize(kBarlRings);
   et_spectrum_e_histos.resize(kEndcEtaRings);
   e_spectrum_e_histos.resize(kEndcEtaRings);
  
   spectra = true;
  
   nevents_ = 0;
   eventsinrun_ = 0;
   eventsinlb_ = 0;
 }

References e_spectrum_b_histos, e_spectrum_e_histos, et_spectrum_b_histos, et_spectrum_e_histos, eventsinlb_, eventsinrun_, isfirstpass_, kBarlRings, kEndcEtaRings, nevents_, and spectra.

◆ ~PhiSymmetryCalibration()

PhiSymmetryCalibration::~PhiSymmetryCalibration ( )

override

Destructor.

Definition at line 79 of file PhiSymmetryCalibration.cc.

                                                 {
   for (Int_t i = 0; i < kBarlRings; i++) {
     delete et_spectrum_b_histos[i];
     delete e_spectrum_b_histos[i];
   }
   for (Int_t i = 0; i < kEndcEtaRings; i++) {
     delete et_spectrum_e_histos[i];
     delete e_spectrum_e_histos[i];
   }
 }

References e_spectrum_b_histos, e_spectrum_e_histos, et_spectrum_b_histos, et_spectrum_e_histos, mps_fire::i, kBarlRings, and kEndcEtaRings.

Member Function Documentation

◆ analyze()

void PhiSymmetryCalibration::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Called at each event.

Implements edm::EDAnalyzer.

Definition at line 230 of file PhiSymmetryCalibration.cc.

                                                                                       {
   using namespace edm;
   using namespace std;
  
   if (isfirstpass_) {
     setUp(setup);
     isfirstpass_ = false;
   }
  
   Handle<EBRecHitCollection> barrelRecHitsHandle;
   Handle<EERecHitCollection> endcapRecHitsHandle;
  
   event.getByLabel(ecalHitsProducer_, barrelHits_, barrelRecHitsHandle);
   if (!barrelRecHitsHandle.isValid()) {
     LogError("") << "[PhiSymmetryCalibration] Error! Can't get product!" << std::endl;
   }
  
   event.getByLabel(ecalHitsProducer_, endcapHits_, endcapRecHitsHandle);
   if (!endcapRecHitsHandle.isValid()) {
     LogError("") << "[PhiSymmetryCalibration] Error! Can't get product!" << std::endl;
   }
  
   // get the ecal geometry
   edm::ESHandle<CaloGeometry> geoHandle;
   setup.get<CaloGeometryRecord>().get(geoHandle);
   const CaloSubdetectorGeometry* barrelGeometry = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
   const CaloSubdetectorGeometry* endcapGeometry = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
  
   bool pass = false;
   // select interesting EcalRecHits (barrel)
   EBRecHitCollection::const_iterator itb;
   for (itb = barrelRecHitsHandle->begin(); itb != barrelRecHitsHandle->end(); itb++) {
     EBDetId hit = EBDetId(itb->id());
     float eta = barrelGeometry->getGeometry(hit)->getPosition().eta();
     float et = itb->energy() / cosh(eta);
     float e = itb->energy();
  
     // if iterating, correct by the previous calib constants found,
     // which are supplied in the form of correction
     if (reiteration_) {
       et = et * oldCalibs_[hit];
       e = e * oldCalibs_[hit];
     }
  
     float et_thr = eCut_barl_ / cosh(eta) + 1.;
  
     int sign = hit.ieta() > 0 ? 1 : 0;
  
     if (e > eCut_barl_ && et < et_thr && e_.goodCell_barl[abs(hit.ieta()) - 1][hit.iphi() - 1][sign]) {
       etsum_barl_[abs(hit.ieta()) - 1][hit.iphi() - 1][sign] += et;
       nhits_barl_[abs(hit.ieta()) - 1][hit.iphi() - 1][sign]++;
       pass = true;
     }  //if energy
  
     if (eventSet_ == 1) {
       // apply a miscalibration to all crystals and increment the
       // ET sum, combined for all crystals
       for (int imiscal = 0; imiscal < kNMiscalBinsEB; imiscal++) {
         if (miscalEB_[imiscal] * e > eCut_barl_ && miscalEB_[imiscal] * et < et_thr &&
             e_.goodCell_barl[abs(hit.ieta()) - 1][hit.iphi() - 1][sign]) {
           etsum_barl_miscal_[imiscal][abs(hit.ieta()) - 1] += miscalEB_[imiscal] * et;
         }
       }
  
       // spectra stuff
       if (spectra && hit.ieta() > 0)  //POSITIVE!!!
       //      if(spectra && hit.ieta()<0) //NEGATIVE!!!
       {
         et_spectrum_b_histos[abs(hit.ieta()) - 1]->Fill(et * 1000.);
         e_spectrum_b_histos[abs(hit.ieta()) - 1]->Fill(e * 1000.);
       }  //if spectra
  
     }  //if eventSet_==1
   }    //for barl
  
   // select interesting EcalRecHits (endcaps)
   EERecHitCollection::const_iterator ite;
   for (ite = endcapRecHitsHandle->begin(); ite != endcapRecHitsHandle->end(); ite++) {
     EEDetId hit = EEDetId(ite->id());
     float eta = abs(endcapGeometry->getGeometry(hit)->getPosition().eta());
     //float phi = endcapGeometry->getGeometry(hit)->getPosition().phi();
  
     float et = ite->energy() / cosh(eta);
     float e = ite->energy();
  
     // if iterating, multiply by the previous correction factor
     if (reiteration_) {
       et = et * oldCalibs_[hit];
       e = e * oldCalibs_[hit];
     }
  
     int sign = hit.zside() > 0 ? 1 : 0;
  
     // changes of eCut_endc_ -> variable linearthr
     // e_cut = ap + eta_ring*b
  
     double eCut_endc = 0;
     for (int ring = 0; ring < kEndcEtaRings; ring++) {
       if (eta > e_.etaBoundary_[ring] && eta < e_.etaBoundary_[ring + 1]) {
         float eta_ring = abs(e_.cellPos_[ring][50].eta());
         eCut_endc = ap_ + eta_ring * b_;
       }
     }
  
     float et_thr = eCut_endc / cosh(eta) + 1.;
  
     if (e > eCut_endc && et < et_thr && e_.goodCell_endc[hit.ix() - 1][hit.iy() - 1][sign]) {
       etsum_endc_[hit.ix() - 1][hit.iy() - 1][sign] += et;
       nhits_endc_[hit.ix() - 1][hit.iy() - 1][sign]++;
       pass = true;
     }
  
     if (eventSet_ == 1) {
       // apply a miscalibration to all crystals and increment the
       // ET sum, combined for all crystals
       for (int imiscal = 0; imiscal < kNMiscalBinsEE; imiscal++) {
         if (miscalEE_[imiscal] * e > eCut_endc && et * miscalEE_[imiscal] < et_thr &&
             e_.goodCell_endc[hit.ix() - 1][hit.iy() - 1][sign]) {
           int ring = e_.endcapRing_[hit.ix() - 1][hit.iy() - 1];
           etsum_endc_miscal_[imiscal][ring] += miscalEE_[imiscal] * et;
         }
       }
  
       // spectra stuff
       if (spectra && hit.zside() > 0)  //POSITIVE!!!
  
       {
         int ring = e_.endcapRing_[hit.ix() - 1][hit.iy() - 1];
  
         et_spectrum_e_histos[ring]->Fill(et * 1000.);
         e_spectrum_e_histos[ring]->Fill(e * 1000.);
  
         if (ring == 16) {
           //int iphi_endc = 0;
           for (int ip = 0; ip < e_.nRing_[ring]; ip++) {
             //if (phi==e_.phi_endc_[ip][ring]) iphi_endc=ip;
           }
         }
       }  //if spectra
  
     }  //if eventSet_==1
   }    //for endc
  
   if (pass) {
     nevents_++;
     eventsinrun_++;
     eventsinlb_++;
   }
 }

References funct::abs(), ap_, b_, barrelHits_, edm::SortedCollection< T, SORT >::begin(), EcalGeomPhiSymHelper::cellPos_, MillePedeFileConverter_cfg::e, e_, e_spectrum_b_histos, e_spectrum_e_histos, DetId::Ecal, EcalBarrel, EcalEndcap, ecalHitsProducer_, eCut_barl_, edm::SortedCollection< T, SORT >::end(), endcapHits_, EcalGeomPhiSymHelper::endcapRing_, EgHLTOffHistBins_cfi::et, et_spectrum_b_histos, et_spectrum_e_histos, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), EcalGeomPhiSymHelper::etaBoundary_, etsum_barl_, etsum_barl_miscal_, etsum_endc_, etsum_endc_miscal_, eventSet_, eventsinlb_, eventsinrun_, HcalObjRepresent::Fill(), get, CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), EcalGeomPhiSymHelper::goodCell_barl, EcalGeomPhiSymHelper::goodCell_endc, isfirstpass_, edm::HandleBase::isValid(), kEndcEtaRings, kNMiscalBinsEB, kNMiscalBinsEE, miscalEB_, miscalEE_, nevents_, nhits_barl_, nhits_endc_, EcalGeomPhiSymHelper::nRing_, oldCalibs_, reiteration_, relativeConstraints::ring, singleTopDQM_cfi::setup, setUp(), Validation_hcalonly_cfi::sign, and spectra.

◆ beginJob()

void PhiSymmetryCalibration::beginJob ( void )

overridevirtual

Called at beginning of job.

Reimplemented from edm::EDAnalyzer.

Definition at line 93 of file PhiSymmetryCalibration.cc.

                                       {
   // initialize arrays
   for (int sign = 0; sign < kSides; sign++) {
     for (int ieta = 0; ieta < kBarlRings; ieta++) {
       for (int iphi = 0; iphi < kBarlWedges; iphi++) {
         etsum_barl_[ieta][iphi][sign] = 0.;
         nhits_barl_[ieta][iphi][sign] = 0;
       }
     }
     for (int ix = 0; ix < kEndcWedgesX; ix++) {
       for (int iy = 0; iy < kEndcWedgesY; iy++) {
         etsum_endc_[ix][iy][sign] = 0.;
         nhits_endc_[ix][iy][sign] = 0;
       }
     }
   }
  
   for (int imiscal = 0; imiscal < kNMiscalBinsEB; imiscal++) {
     miscalEB_[imiscal] = (1 - kMiscalRangeEB) + float(imiscal) * (2 * kMiscalRangeEB / (kNMiscalBinsEB - 1));
     for (int ieta = 0; ieta < kBarlRings; ieta++)
       etsum_barl_miscal_[imiscal][ieta] = 0.;
   }
  
   for (int imiscal = 0; imiscal < kNMiscalBinsEE; imiscal++) {
     miscalEE_[imiscal] = (1 - kMiscalRangeEE) + float(imiscal) * (2 * kMiscalRangeEE / (kNMiscalBinsEE - 1));
     for (int ring = 0; ring < kEndcEtaRings; ring++)
       etsum_endc_miscal_[imiscal][ring] = 0.;
   }
  
   // start spectra stuff
   if (eventSet_ != 1)
     spectra = false;
  
   if (spectra) {
     ostringstream t;
     for (Int_t i = 0; i < kBarlRings; i++) {
       t << "et_spectrum_b_" << i + 1;
       et_spectrum_b_histos[i] = new TH1F(t.str().c_str(), ";E_{T} [MeV]", 50, 0., 500.);
       t.str("");
  
       t << "e_spectrum_b_" << i + 1;
       e_spectrum_b_histos[i] = new TH1F(t.str().c_str(), ";E [MeV]", 50, 0., 500.);
       t.str("");
     }
     for (Int_t i = 0; i < kEndcEtaRings; i++) {
       t << "et_spectrum_e_" << i + 1;
       et_spectrum_e_histos[i] = new TH1F(t.str().c_str(), ";E_{T} [MeV]", 75, 0., 1500.);
       t.str("");
  
       t << "e_spectrum_e_" << i + 1;
       e_spectrum_e_histos[i] = new TH1F(t.str().c_str(), ";E [MeV]", 75, 0., 1500.);
       t.str("");
     }
   }
   // end spectra stuff
 }

References e_spectrum_b_histos, e_spectrum_e_histos, et_spectrum_b_histos, et_spectrum_e_histos, etsum_barl_, etsum_barl_miscal_, etsum_endc_, etsum_endc_miscal_, eventSet_, mps_fire::i, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, kBarlRings, kBarlWedges, kEndcEtaRings, kEndcWedgesX, kEndcWedgesY, kMiscalRangeEB, kMiscalRangeEE, kNMiscalBinsEB, kNMiscalBinsEE, kSides, miscalEB_, miscalEE_, nhits_barl_, nhits_endc_, relativeConstraints::ring, Validation_hcalonly_cfi::sign, spectra, and submitPVValidationJobs::t.

◆ endJob()

void PhiSymmetryCalibration::endJob ( void )

overridevirtual

Called at end of job.

Reimplemented from edm::EDAnalyzer.

Definition at line 153 of file PhiSymmetryCalibration.cc.

                                     {
   edm::LogInfo("Calibration") << "[PhiSymmetryCalibration] At end of job";
  
   // start spectra stuff
   if (spectra) {
     TFile f("Espectra_plus.root", "recreate");
  
     for (int i = 0; i < kBarlRings; i++) {
       et_spectrum_b_histos[i]->Write();
       e_spectrum_b_histos[i]->Write();
     }
  
     for (int i = 0; i < kEndcEtaRings; i++) {
       et_spectrum_e_histos[i]->Write();
       e_spectrum_e_histos[i]->Write();
     }
  
     f.Close();
   }
  
   if (eventSet_ == 1) {
     // calculate factors to convert from fractional deviation of ET sum from
     // the mean to the estimate of the miscalibration factor
     getKfactors();
  
     std::ofstream k_barl_out("k_barl.dat", ios::out);
     for (int ieta = 0; ieta < kBarlRings; ieta++)
       k_barl_out << ieta << " " << k_barl_[ieta] << endl;
     k_barl_out.close();
  
     std::ofstream k_endc_out("k_endc.dat", ios::out);
     for (int ring = 0; ring < kEndcEtaRings; ring++)
       k_endc_out << ring << " " << k_endc_[ring] << endl;
     k_endc_out.close();
   }
  
   if (eventSet_ != 0) {
     // output ET sums
  
     stringstream etsum_file_barl;
     etsum_file_barl << "etsum_barl_" << eventSet_ << ".dat";
  
     std::ofstream etsum_barl_out(etsum_file_barl.str().c_str(), ios::out);
  
     for (int ieta = 0; ieta < kBarlRings; ieta++) {
       for (int iphi = 0; iphi < kBarlWedges; iphi++) {
         for (int sign = 0; sign < kSides; sign++) {
           etsum_barl_out << eventSet_ << " " << ieta << " " << iphi << " " << sign << " "
                          << etsum_barl_[ieta][iphi][sign] << " " << nhits_barl_[ieta][iphi][sign] << endl;
         }
       }
     }
     etsum_barl_out.close();
  
     stringstream etsum_file_endc;
     etsum_file_endc << "etsum_endc_" << eventSet_ << ".dat";
  
     std::ofstream etsum_endc_out(etsum_file_endc.str().c_str(), ios::out);
     for (int ix = 0; ix < kEndcWedgesX; ix++) {
       for (int iy = 0; iy < kEndcWedgesY; iy++) {
         int ring = e_.endcapRing_[ix][iy];
         if (ring != -1) {
           for (int sign = 0; sign < kSides; sign++) {
             etsum_endc_out << eventSet_ << " " << ix << " " << iy << " " << sign << " " << etsum_endc_[ix][iy][sign]
                            << " " << nhits_endc_[ix][iy][sign] << " " << e_.endcapRing_[ix][iy] << endl;
           }
         }
       }
     }
     etsum_endc_out.close();
   }
   cout << "Events processed " << nevents_ << endl;
 }

References gather_cfg::cout, e_, e_spectrum_b_histos, e_spectrum_e_histos, EcalGeomPhiSymHelper::endcapRing_, et_spectrum_b_histos, et_spectrum_e_histos, etsum_barl_, etsum_endc_, eventSet_, f, getKfactors(), mps_fire::i, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, k_barl_, k_endc_, kBarlRings, kBarlWedges, kEndcEtaRings, kEndcWedgesX, kEndcWedgesY, kSides, nevents_, nhits_barl_, nhits_endc_, MillePedeFileConverter_cfg::out, relativeConstraints::ring, Validation_hcalonly_cfi::sign, and spectra.

Referenced by o2olib.O2ORunMgr::executeJob().

◆ endLuminosityBlock()

void PhiSymmetryCalibration::endLuminosityBlock	(	edm::LuminosityBlock const &	lb,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 511 of file PhiSymmetryCalibration.cc.

                                                                                                   {
   if ((lb.endTime().value() >> 32) - (lb.beginTime().value() >> 32) < 60)
     return;
  
   std::cout << "PHILB : run " << lb.run() << " id " << lb.id() << " start " << (lb.beginTime().value() >> 32) << " end "
             << (lb.endTime().value() >> 32) << " dur " << (lb.endTime().value() >> 32) - (lb.beginTime().value() >> 32)
  
             << " npass " << eventsinlb_ << std::endl;
  
   eventsinlb_ = 0;
 }

References edm::LuminosityBlockBase::beginTime(), gather_cfg::cout, edm::LuminosityBlockBase::endTime(), eventsinlb_, edm::LuminosityBlockBase::id(), edm::LuminosityBlockBase::run(), and edm::Timestamp::value().

◆ endRun()

void PhiSymmetryCalibration::endRun	(	edm::Run const &	run,
		const edm::EventSetup &
	)

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 380 of file PhiSymmetryCalibration.cc.

                                                                            {
   std::cout << "PHIREPRT : run " << run.run() << " start " << (run.beginTime().value() >> 32) << " end "
             << (run.endTime().value() >> 32) << " dur "
             << (run.endTime().value() >> 32) - (run.beginTime().value() >> 32)
  
             << " npass " << eventsinrun_ << std::endl;
   eventsinrun_ = 0;
  
   return;
 }

References gather_cfg::cout, eventsinrun_, and writedatasetfile::run.

◆ getKfactors()

void PhiSymmetryCalibration::getKfactors ( )

private

Definition at line 393 of file PhiSymmetryCalibration.cc.

                                          {
   float epsilon_T_eb[kNMiscalBinsEB];
   float epsilon_M_eb[kNMiscalBinsEB];
  
   float epsilon_T_ee[kNMiscalBinsEE];
   float epsilon_M_ee[kNMiscalBinsEE];
  
   std::vector<TGraph*> k_barl_graph(kBarlRings);
   std::vector<TCanvas*> k_barl_plot(kBarlRings);
  
   //Create our own TF1 to avoid threading problems
   TF1 mypol1("mypol1", "pol1");
   for (int ieta = 0; ieta < kBarlRings; ieta++) {
     for (int imiscal = 0; imiscal < kNMiscalBinsEB; imiscal++) {
       int middlebin = int(kNMiscalBinsEB / 2);
       epsilon_T_eb[imiscal] = etsum_barl_miscal_[imiscal][ieta] / etsum_barl_miscal_[middlebin][ieta] - 1.;
       epsilon_M_eb[imiscal] = miscalEB_[imiscal] - 1.;
     }
     k_barl_graph[ieta] = new TGraph(kNMiscalBinsEB, epsilon_M_eb, epsilon_T_eb);
     k_barl_graph[ieta]->Fit(&mypol1);
  
     ostringstream t;
     t << "k_barl_" << ieta + 1;
     k_barl_plot[ieta] = new TCanvas(t.str().c_str(), "");
     k_barl_plot[ieta]->SetFillColor(10);
     k_barl_plot[ieta]->SetGrid();
     k_barl_graph[ieta]->SetMarkerSize(1.);
     k_barl_graph[ieta]->SetMarkerColor(4);
     k_barl_graph[ieta]->SetMarkerStyle(20);
     k_barl_graph[ieta]->GetXaxis()->SetLimits(-1. * kMiscalRangeEB, kMiscalRangeEB);
     k_barl_graph[ieta]->GetXaxis()->SetTitleSize(.05);
     k_barl_graph[ieta]->GetYaxis()->SetTitleSize(.05);
     k_barl_graph[ieta]->GetXaxis()->SetTitle("#epsilon_{M}");
     k_barl_graph[ieta]->GetYaxis()->SetTitle("#epsilon_{T}");
     k_barl_graph[ieta]->Draw("AP");
  
     k_barl_[ieta] = k_barl_graph[ieta]->GetFunction("pol1")->GetParameter(1);
     std::cout << "k_barl_[" << ieta << "]=" << k_barl_[ieta] << std::endl;
   }
  
   std::vector<TGraph*> k_endc_graph(kEndcEtaRings);
   std::vector<TCanvas*> k_endc_plot(kEndcEtaRings);
  
   for (int ring = 0; ring < kEndcEtaRings; ring++) {
     for (int imiscal = 0; imiscal < kNMiscalBinsEE; imiscal++) {
       int middlebin = int(kNMiscalBinsEE / 2);
       epsilon_T_ee[imiscal] = etsum_endc_miscal_[imiscal][ring] / etsum_endc_miscal_[middlebin][ring] - 1.;
       epsilon_M_ee[imiscal] = miscalEE_[imiscal] - 1.;
     }
     k_endc_graph[ring] = new TGraph(kNMiscalBinsEE, epsilon_M_ee, epsilon_T_ee);
     k_endc_graph[ring]->Fit(&mypol1);
  
     ostringstream t;
     t << "k_endc_" << ring + 1;
     k_endc_plot[ring] = new TCanvas(t.str().c_str(), "");
     k_endc_plot[ring]->SetFillColor(10);
     k_endc_plot[ring]->SetGrid();
     k_endc_graph[ring]->SetMarkerSize(1.);
     k_endc_graph[ring]->SetMarkerColor(4);
     k_endc_graph[ring]->SetMarkerStyle(20);
     k_endc_graph[ring]->GetXaxis()->SetLimits(-1 * kMiscalRangeEE, kMiscalRangeEE);
     k_endc_graph[ring]->GetXaxis()->SetTitleSize(.05);
     k_endc_graph[ring]->GetYaxis()->SetTitleSize(.05);
     k_endc_graph[ring]->GetXaxis()->SetTitle("#epsilon_{M}");
     k_endc_graph[ring]->GetYaxis()->SetTitle("#epsilon_{T}");
     k_endc_graph[ring]->Draw("AP");
  
     k_endc_[ring] = k_endc_graph[ring]->GetFunction("pol1")->GetParameter(1);
     std::cout << "k_endc_[" << ring << "]=" << k_endc_[ring] << std::endl;
   }
  
   TFile f("PhiSymmetryCalibration_kFactors.root", "recreate");
   for (int ieta = 0; ieta < kBarlRings; ieta++) {
     k_barl_plot[ieta]->Write();
     delete k_barl_plot[ieta];
     delete k_barl_graph[ieta];
   }
   for (int ring = 0; ring < kEndcEtaRings; ring++) {
     k_endc_plot[ring]->Write();
     delete k_endc_plot[ring];
     delete k_endc_graph[ring];
   }
   f.Close();
 }

References gather_cfg::cout, etsum_barl_miscal_, etsum_endc_miscal_, f, LEDCalibrationChannels::ieta, createfilelist::int, k_barl_, k_endc_, kBarlRings, kEndcEtaRings, kMiscalRangeEB, kMiscalRangeEE, kNMiscalBinsEB, kNMiscalBinsEE, miscalEB_, miscalEE_, relativeConstraints::ring, and submitPVValidationJobs::t.

Referenced by endJob().

◆ setUp()

void PhiSymmetryCalibration::setUp ( const edm::EventSetup & setup )

Definition at line 480 of file PhiSymmetryCalibration.cc.

                                                              {
   edm::ESHandle<EcalChannelStatus> chStatus;
   setup.get<EcalChannelStatusRcd>().get(chStatus);
  
   edm::ESHandle<CaloGeometry> geoHandle;
   setup.get<CaloGeometryRecord>().get(geoHandle);
  
   e_.setup(&(*geoHandle), &(*chStatus), statusThreshold_);
  
   if (reiteration_) {
     EcalCondHeader h;
     // namespace fs = boost::filesystem;
     //     fs::path p(oldcalibfile_.c_str(),fs::native);
     //     if (!fs::exists(p)) edm::LogError("PhiSym") << "File not found: "
     //                      << oldcalibfile_ <<endl;
  
     edm::FileInPath fip("Calibration/EcalCalibAlgos/data/" + oldcalibfile_);
  
     int ret = EcalIntercalibConstantsXMLTranslator::readXML(fip.fullPath(), h, oldCalibs_);
     if (ret)
       edm::LogError("PhiSym") << "Error reading XML files" << endl;
     ;
  
   } else {
     // in fact if not reiterating, oldCalibs_ will never be used
     edm::ESHandle<EcalIntercalibConstants> pIcal;
     setup.get<EcalIntercalibConstantsRcd>().get(pIcal);
     oldCalibs_ = *pIcal;
   }
 }

References e_, edm::FileInPath::fullPath(), get, h, oldcalibfile_, oldCalibs_, EcalFloatCondObjectContainerXMLTranslator::readXML(), reiteration_, runTheMatrix::ret, EcalGeomPhiSymHelper::setup(), singleTopDQM_cfi::setup, and statusThreshold_.

Referenced by analyze().

Member Data Documentation

◆ ap_

double PhiSymmetryCalibration::ap_

private

Definition at line 122 of file PhiSymmetryCalibration.h.

Referenced by analyze().

◆ b_

double PhiSymmetryCalibration::b_

private

Definition at line 123 of file PhiSymmetryCalibration.h.

Referenced by analyze().

◆ barrelCells

std::vector<DetId> PhiSymmetryCalibration::barrelCells

private

Definition at line 93 of file PhiSymmetryCalibration.h.

◆ barrelHits_

std::string PhiSymmetryCalibration::barrelHits_

private

Definition at line 115 of file PhiSymmetryCalibration.h.

Referenced by analyze().

◆ e_

EcalGeomPhiSymHelper PhiSymmetryCalibration::e_

private

Definition at line 69 of file PhiSymmetryCalibration.h.

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

◆ e_spectrum_b_histos

std::vector<TH1F*> PhiSymmetryCalibration::e_spectrum_b_histos

private

Definition at line 143 of file PhiSymmetryCalibration.h.

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

◆ e_spectrum_e_histos

std::vector<TH1F*> PhiSymmetryCalibration::e_spectrum_e_histos

private

Definition at line 145 of file PhiSymmetryCalibration.h.

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

◆ ecalHitsProducer_

std::string PhiSymmetryCalibration::ecalHitsProducer_

private

Definition at line 114 of file PhiSymmetryCalibration.h.

Referenced by analyze().

◆ eCut_barl_

double PhiSymmetryCalibration::eCut_barl_

private

Definition at line 119 of file PhiSymmetryCalibration.h.

Referenced by analyze().

◆ endcapCells

std::vector<DetId> PhiSymmetryCalibration::endcapCells

private

Definition at line 94 of file PhiSymmetryCalibration.h.

◆ endcapHits_

std::string PhiSymmetryCalibration::endcapHits_

private

Definition at line 116 of file PhiSymmetryCalibration.h.

Referenced by analyze().

◆ epsilon_M_barl

float PhiSymmetryCalibration::epsilon_M_barl[kBarlRings][kBarlWedges][kSides]

private

Definition at line 105 of file PhiSymmetryCalibration.h.

◆ epsilon_M_endc

float PhiSymmetryCalibration::epsilon_M_endc[kEndcWedgesX][kEndcWedgesY][kSides]

private

Definition at line 106 of file PhiSymmetryCalibration.h.

◆ esumMean_barl_

double PhiSymmetryCalibration::esumMean_barl_[kBarlRings]

private

Definition at line 84 of file PhiSymmetryCalibration.h.

◆ esumMean_endc_

double PhiSymmetryCalibration::esumMean_endc_[kEndcEtaRings]

private

Definition at line 85 of file PhiSymmetryCalibration.h.

◆ et_spectrum_b_histos

std::vector<TH1F*> PhiSymmetryCalibration::et_spectrum_b_histos

private

Definition at line 142 of file PhiSymmetryCalibration.h.

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

◆ et_spectrum_e_histos

std::vector<TH1F*> PhiSymmetryCalibration::et_spectrum_e_histos

private

Definition at line 144 of file PhiSymmetryCalibration.h.

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

◆ etsum_barl_

double PhiSymmetryCalibration::etsum_barl_[kBarlRings][kBarlWedges][kSides]

private

Definition at line 72 of file PhiSymmetryCalibration.h.

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

◆ etsum_barl_miscal_

double PhiSymmetryCalibration::etsum_barl_miscal_[kNMiscalBinsEB][kBarlRings]

private

Definition at line 81 of file PhiSymmetryCalibration.h.

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

◆ etsum_endc_

double PhiSymmetryCalibration::etsum_endc_[kEndcWedgesX][kEndcWedgesX][kSides]

private

Definition at line 73 of file PhiSymmetryCalibration.h.

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

◆ etsum_endc_miscal_

double PhiSymmetryCalibration::etsum_endc_miscal_[kNMiscalBinsEE][kEndcEtaRings]

private

Definition at line 82 of file PhiSymmetryCalibration.h.

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

◆ etsum_endc_uncorr

double PhiSymmetryCalibration::etsum_endc_uncorr[kEndcWedgesX][kEndcWedgesX][kSides]

private

Definition at line 74 of file PhiSymmetryCalibration.h.

◆ etsumMean_barl_

double PhiSymmetryCalibration::etsumMean_barl_[kBarlRings]

private

Definition at line 75 of file PhiSymmetryCalibration.h.

◆ etsumMean_endc_

double PhiSymmetryCalibration::etsumMean_endc_[kEndcEtaRings]

private

Definition at line 76 of file PhiSymmetryCalibration.h.

◆ eventSet_

int PhiSymmetryCalibration::eventSet_

private

Definition at line 125 of file PhiSymmetryCalibration.h.

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

◆ eventsinlb_

int PhiSymmetryCalibration::eventsinlb_

private

Definition at line 150 of file PhiSymmetryCalibration.h.

Referenced by analyze(), endLuminosityBlock(), and PhiSymmetryCalibration().

◆ eventsinrun_

int PhiSymmetryCalibration::eventsinrun_

private

Definition at line 149 of file PhiSymmetryCalibration.h.

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

◆ isfirstpass_

bool PhiSymmetryCalibration::isfirstpass_

private

Definition at line 139 of file PhiSymmetryCalibration.h.

Referenced by analyze(), and PhiSymmetryCalibration().

◆ k_barl_

double PhiSymmetryCalibration::k_barl_[kBarlRings]

private

Definition at line 88 of file PhiSymmetryCalibration.h.

Referenced by endJob(), and getKfactors().

◆ k_endc_

double PhiSymmetryCalibration::k_endc_[kEndcEtaRings]

private

Definition at line 89 of file PhiSymmetryCalibration.h.

Referenced by endJob(), and getKfactors().

◆ kMaxEndciPhi

const int PhiSymmetryCalibration::kMaxEndciPhi = 360

staticprivate

Definition at line 129 of file PhiSymmetryCalibration.h.

◆ kMiscalRangeEB

const float PhiSymmetryCalibration::kMiscalRangeEB = .05

staticprivate

Definition at line 64 of file PhiSymmetryCalibration.h.

Referenced by beginJob(), and getKfactors().

◆ kMiscalRangeEE

const float PhiSymmetryCalibration::kMiscalRangeEE = .10

staticprivate

Definition at line 67 of file PhiSymmetryCalibration.h.

Referenced by beginJob(), and getKfactors().

◆ kNMiscalBinsEB

const int PhiSymmetryCalibration::kNMiscalBinsEB = 21

staticprivate

Definition at line 63 of file PhiSymmetryCalibration.h.

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

◆ kNMiscalBinsEE

const int PhiSymmetryCalibration::kNMiscalBinsEE = 41

staticprivate

Definition at line 66 of file PhiSymmetryCalibration.h.

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

◆ miscalEB_

double PhiSymmetryCalibration::miscalEB_[kNMiscalBinsEB]

private

Definition at line 90 of file PhiSymmetryCalibration.h.

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

◆ miscalEE_

double PhiSymmetryCalibration::miscalEE_[kNMiscalBinsEE]

private

Definition at line 91 of file PhiSymmetryCalibration.h.

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

◆ nevents_

int PhiSymmetryCalibration::nevents_

private

Definition at line 148 of file PhiSymmetryCalibration.h.

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

◆ newCalibs_barl

double PhiSymmetryCalibration::newCalibs_barl[kBarlRings][kBarlWedges][kSides]

private

Definition at line 101 of file PhiSymmetryCalibration.h.

◆ newCalibs_endc

double PhiSymmetryCalibration::newCalibs_endc[kEndcWedgesX][kEndcWedgesX][kSides]

private

Definition at line 102 of file PhiSymmetryCalibration.h.

◆ nhits_barl_

unsigned int PhiSymmetryCalibration::nhits_barl_[kBarlRings][kBarlWedges][kSides]

private

Definition at line 78 of file PhiSymmetryCalibration.h.

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

◆ nhits_endc_

unsigned int PhiSymmetryCalibration::nhits_endc_[kEndcWedgesX][kEndcWedgesX][kSides]

private

Definition at line 79 of file PhiSymmetryCalibration.h.

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

◆ oldcalibfile_

std::string PhiSymmetryCalibration::oldcalibfile_

private

Definition at line 134 of file PhiSymmetryCalibration.h.

Referenced by setUp().

◆ oldCalibs_

EcalIntercalibConstants PhiSymmetryCalibration::oldCalibs_

private

the old calibration constants (when reiterating, the last ones derived)

Definition at line 137 of file PhiSymmetryCalibration.h.

Referenced by analyze(), and setUp().

◆ oldCalibs_barl

double PhiSymmetryCalibration::oldCalibs_barl[kBarlRings][kBarlWedges][kSides]

private

Definition at line 97 of file PhiSymmetryCalibration.h.

◆ oldCalibs_endc

double PhiSymmetryCalibration::oldCalibs_endc[kEndcWedgesX][kEndcWedgesY][kSides]

private

Definition at line 98 of file PhiSymmetryCalibration.h.

◆ phi_endc

float PhiSymmetryCalibration::phi_endc[kMaxEndciPhi][kEndcEtaRings]

private

Definition at line 131 of file PhiSymmetryCalibration.h.

◆ rawconst_barl

float PhiSymmetryCalibration::rawconst_barl[kBarlRings][kBarlWedges][kSides]

private

Definition at line 109 of file PhiSymmetryCalibration.h.

◆ rawconst_endc

float PhiSymmetryCalibration::rawconst_endc[kEndcWedgesX][kEndcWedgesX][kSides]

private

Definition at line 110 of file PhiSymmetryCalibration.h.

◆ reiteration_

bool PhiSymmetryCalibration::reiteration_

private

Definition at line 133 of file PhiSymmetryCalibration.h.

Referenced by analyze(), and setUp().

◆ spectra

bool PhiSymmetryCalibration::spectra

private

Definition at line 147 of file PhiSymmetryCalibration.h.

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

◆ statusThreshold_

int PhiSymmetryCalibration::statusThreshold_

private

threshold in channel status beyond which channel is marked bad

Definition at line 127 of file PhiSymmetryCalibration.h.

Referenced by setUp().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PhiSymmetryCalibration()

◆ ~PhiSymmetryCalibration()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ endLuminosityBlock()

◆ endRun()

◆ getKfactors()

◆ setUp()

Member Data Documentation

◆ ap_

◆ b_

◆ barrelCells

◆ barrelHits_

◆ e_

◆ e_spectrum_b_histos

◆ e_spectrum_e_histos

◆ ecalHitsProducer_

◆ eCut_barl_

◆ endcapCells

◆ endcapHits_

◆ epsilon_M_barl

◆ epsilon_M_endc

◆ esumMean_barl_

◆ esumMean_endc_

◆ et_spectrum_b_histos

◆ et_spectrum_e_histos

◆ etsum_barl_

◆ etsum_barl_miscal_

◆ etsum_endc_

◆ etsum_endc_miscal_

◆ etsum_endc_uncorr

◆ etsumMean_barl_

◆ etsumMean_endc_

◆ eventSet_

◆ eventsinlb_

◆ eventsinrun_

◆ isfirstpass_

◆ k_barl_

◆ k_endc_

◆ kMaxEndciPhi

◆ kMiscalRangeEB

◆ kMiscalRangeEE

◆ kNMiscalBinsEB

◆ kNMiscalBinsEE

◆ miscalEB_

◆ miscalEE_

◆ nevents_

◆ newCalibs_barl

◆ newCalibs_endc

◆ nhits_barl_

◆ nhits_endc_

◆ oldcalibfile_

◆ oldCalibs_

◆ oldCalibs_barl

◆ oldCalibs_endc

◆ phi_endc

◆ rawconst_barl

◆ rawconst_endc

◆ reiteration_

◆ spectra

◆ statusThreshold_