#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	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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::one::EDAnalyzer< edm::one::SharedResources, edm::one::WatchRuns, edm::one::WatchLuminosityBlocks >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	getKfactors ()

Private Attributes
const double	ap_

const double	b_

std::vector< DetId >	barrelCells

const std::string	barrelHits_

const edm::ESGetToken < EcalChannelStatus, EcalChannelStatusRcd >	channelStatusToken_

EcalGeomPhiSymHelper	e_

std::vector< TH1F * >	e_spectrum_b_histos

std::vector< TH1F * >	e_spectrum_e_histos

const edm::EDGetTokenT < EBRecHitCollection >	ebRecHitToken_

const std::string	ecalHitsProducer_

const double	eCut_barl_

const edm::EDGetTokenT < EERecHitCollection >	eeRecHitToken_

std::vector< DetId >	endcapCells

const 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]

const int	eventSet_

int	eventsinlb_

int	eventsinrun_

const edm::ESGetToken < CaloGeometry, CaloGeometryRecord >	geometryToken_

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]

const 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]

const bool	reiteration_

bool	spectra

const 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::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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

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

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

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

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

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

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

Definition at line 39 of file PhiSymmetryCalibration.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 42 of file PhiSymmetryCalibration.cc.

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.

     : ecalHitsProducer_(iConfig.getParameter<std::string>("ecalRecHitsProducer")),
       barrelHits_(iConfig.getParameter<std::string>("barrelHitCollection")),
       endcapHits_(iConfig.getParameter<std::string>("endcapHitCollection")),
       ebRecHitToken_(consumes<EBRecHitCollection>(edm::InputTag(ecalHitsProducer_, barrelHits_))),
       eeRecHitToken_(consumes<EERecHitCollection>(edm::InputTag(ecalHitsProducer_, endcapHits_))),
       channelStatusToken_(esConsumes()),
       geometryToken_(esConsumes()),
       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;
 
   // because ROOT draws something
   usesResource();
 }

PhiSymmetryCalibration::~PhiSymmetryCalibration ( )

override

Destructor.

Definition at line 77 of file PhiSymmetryCalibration.cc.

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

                                                 {
   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];
   }
 }

Member Function Documentation

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

overridevirtual

Called at each event.

Implements edm::one::EDAnalyzerBase.

Definition at line 228 of file PhiSymmetryCalibration.cc.

References funct::abs(), ap_, b_, EcalGeomPhiSymHelper::cellPos_, alignCSCRings::e, e_, e_spectrum_b_histos, e_spectrum_e_histos, ebRecHitToken_, DetId::Ecal, EcalBarrel, EcalEndcap, eCut_barl_, EcalGeomPhiSymHelper::endcapRing_, 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(), relativeConstraints::geometry, geometryToken_, edm::EventSetup::getData(), EcalGeomPhiSymHelper::goodCell_barl, EcalGeomPhiSymHelper::goodCell_endc, EBDetId::ieta(), EBDetId::iphi(), isfirstpass_, EEDetId::ix(), EEDetId::iy(), kEndcEtaRings, kNMiscalBinsEB, kNMiscalBinsEE, miscalEB_, miscalEE_, nevents_, nhits_barl_, nhits_endc_, EcalGeomPhiSymHelper::nRing_, oldCalibs_, reiteration_, relativeConstraints::ring, setUp(), jetcorrextractor::sign(), spectra, and EEDetId::zside().

                                                                                       {
   using namespace edm;
   using namespace std;
 
   if (isfirstpass_) {
     setUp(setup);
     isfirstpass_ = false;
   }
 
   Handle<EBRecHitCollection> barrelRecHitsHandle;
   Handle<EERecHitCollection> endcapRecHitsHandle;
 
   event.getByToken(ebRecHitToken_, barrelRecHitsHandle);
   if (!barrelRecHitsHandle.isValid()) {
     LogError("") << "[PhiSymmetryCalibration] Error! Can't get product!" << std::endl;
   }
 
   event.getByToken(ebRecHitToken_, endcapRecHitsHandle);
   if (!endcapRecHitsHandle.isValid()) {
     LogError("") << "[PhiSymmetryCalibration] Error! Can't get product!" << std::endl;
   }
 
   // get the ecal geometry
   const auto& geometry = setup.getData(geometryToken_);
   const auto* barrelGeometry = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
   const auto* endcapGeometry = geometry.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_++;
   }
 }

void PhiSymmetryCalibration::beginJob ( void )

overridevirtual

Called at beginning of job.

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 91 of file PhiSymmetryCalibration.cc.

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, kBarlRings, kBarlWedges, kEndcEtaRings, kEndcWedgesX, kEndcWedgesY, kMiscalRangeEB, kMiscalRangeEE, kNMiscalBinsEB, kNMiscalBinsEE, kSides, miscalEB_, miscalEE_, nhits_barl_, nhits_endc_, relativeConstraints::ring, jetcorrextractor::sign(), spectra, and submitPVValidationJobs::t.

                                       {
   // 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
 }

void PhiSymmetryCalibration::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

override

Definition at line 501 of file PhiSymmetryCalibration.cc.

501 {}

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

override

Definition at line 377 of file PhiSymmetryCalibration.cc.

377 {}

void PhiSymmetryCalibration::endJob ( void )

overridevirtual

Called at end of job.

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 151 of file PhiSymmetryCalibration.cc.

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_, validate-o2o-wbm::f, getKfactors(), mps_fire::i, k_barl_, k_endc_, kBarlRings, kBarlWedges, kEndcEtaRings, kEndcWedgesX, kEndcWedgesY, kSides, nevents_, nhits_barl_, nhits_endc_, submitPVResolutionJobs::out, relativeConstraints::ring, jetcorrextractor::sign(), and spectra.

Referenced by o2olib.O2ORunMgr::executeJob().

                                     {
   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;
 }

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

override

Definition at line 503 of file PhiSymmetryCalibration.cc.

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

                                                                                                   {
   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;
 }

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

override

Definition at line 379 of file PhiSymmetryCalibration.cc.

References edm::RunBase::beginTime(), gather_cfg::cout, edm::RunBase::endTime(), eventsinrun_, edm::RunBase::run(), and edm::Timestamp::value().

                                                                            {
   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;
 }

void PhiSymmetryCalibration::getKfactors ( )

private

Definition at line 392 of file PhiSymmetryCalibration.cc.

References gather_cfg::cout, etsum_barl_miscal_, etsum_endc_miscal_, validate-o2o-wbm::f, k_barl_, k_endc_, kBarlRings, kEndcEtaRings, kMiscalRangeEB, kMiscalRangeEE, kNMiscalBinsEB, kNMiscalBinsEE, miscalEB_, miscalEE_, relativeConstraints::ring, and submitPVValidationJobs::t.

Referenced by endJob().

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

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

Definition at line 479 of file PhiSymmetryCalibration.cc.

References channelStatusToken_, e_, edm::FileInPath::fullPath(), relativeConstraints::geometry, geometryToken_, edm::EventSetup::getData(), h, oldcalibfile_, oldCalibs_, EcalFloatCondObjectContainerXMLTranslator::readXML(), reiteration_, runTheMatrix::ret, EcalGeomPhiSymHelper::setup(), and statusThreshold_.

Referenced by analyze().

                                                              {
   const auto& chStatus = setup.getData(channelStatusToken_);
 
   const auto& geometry = setup.getData(geometryToken_);
 
   e_.setup(&geometry, &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;
   }
 }

Member Data Documentation

const double PhiSymmetryCalibration::ap_

private

Definition at line 135 of file PhiSymmetryCalibration.h.

Referenced by analyze().

const double PhiSymmetryCalibration::b_

private

Definition at line 136 of file PhiSymmetryCalibration.h.

Referenced by analyze().

std::vector<DetId> PhiSymmetryCalibration::barrelCells

private

Definition at line 101 of file PhiSymmetryCalibration.h.

const std::string PhiSymmetryCalibration::barrelHits_

private

Definition at line 123 of file PhiSymmetryCalibration.h.

const edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> PhiSymmetryCalibration::channelStatusToken_

private

Definition at line 128 of file PhiSymmetryCalibration.h.

Referenced by setUp().

EcalGeomPhiSymHelper PhiSymmetryCalibration::e_

private

Definition at line 77 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 156 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 158 of file PhiSymmetryCalibration.h.

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

const edm::EDGetTokenT<EBRecHitCollection> PhiSymmetryCalibration::ebRecHitToken_

private

Definition at line 126 of file PhiSymmetryCalibration.h.

Referenced by analyze().

const std::string PhiSymmetryCalibration::ecalHitsProducer_

private

Definition at line 122 of file PhiSymmetryCalibration.h.

const double PhiSymmetryCalibration::eCut_barl_

private

Definition at line 132 of file PhiSymmetryCalibration.h.

Referenced by analyze().

const edm::EDGetTokenT<EERecHitCollection> PhiSymmetryCalibration::eeRecHitToken_

private

Definition at line 127 of file PhiSymmetryCalibration.h.

std::vector<DetId> PhiSymmetryCalibration::endcapCells

private

Definition at line 102 of file PhiSymmetryCalibration.h.

const std::string PhiSymmetryCalibration::endcapHits_

private

Definition at line 124 of file PhiSymmetryCalibration.h.

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

private

Definition at line 113 of file PhiSymmetryCalibration.h.

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

private

Definition at line 114 of file PhiSymmetryCalibration.h.

double PhiSymmetryCalibration::esumMean_barl_[kBarlRings]

private

Definition at line 92 of file PhiSymmetryCalibration.h.

double PhiSymmetryCalibration::esumMean_endc_[kEndcEtaRings]

private

Definition at line 93 of file PhiSymmetryCalibration.h.

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

private

Definition at line 155 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 157 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 80 of file PhiSymmetryCalibration.h.

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

double PhiSymmetryCalibration::etsum_barl_miscal_[kNMiscalBinsEB][kBarlRings]

private

Definition at line 89 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 81 of file PhiSymmetryCalibration.h.

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

double PhiSymmetryCalibration::etsum_endc_miscal_[kNMiscalBinsEE][kEndcEtaRings]

private

Definition at line 90 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 82 of file PhiSymmetryCalibration.h.

double PhiSymmetryCalibration::etsumMean_barl_[kBarlRings]

private

Definition at line 83 of file PhiSymmetryCalibration.h.

double PhiSymmetryCalibration::etsumMean_endc_[kEndcEtaRings]

private

Definition at line 84 of file PhiSymmetryCalibration.h.

const int PhiSymmetryCalibration::eventSet_

private

Definition at line 138 of file PhiSymmetryCalibration.h.

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

int PhiSymmetryCalibration::eventsinlb_

private

Definition at line 163 of file PhiSymmetryCalibration.h.

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

int PhiSymmetryCalibration::eventsinrun_

private

Definition at line 162 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 129 of file PhiSymmetryCalibration.h.

Referenced by analyze(), and setUp().

bool PhiSymmetryCalibration::isfirstpass_

private

Definition at line 152 of file PhiSymmetryCalibration.h.

Referenced by analyze(), and PhiSymmetryCalibration().

double PhiSymmetryCalibration::k_barl_[kBarlRings]

private

Definition at line 96 of file PhiSymmetryCalibration.h.

Referenced by endJob(), and getKfactors().

double PhiSymmetryCalibration::k_endc_[kEndcEtaRings]

private

Definition at line 97 of file PhiSymmetryCalibration.h.

Referenced by endJob(), and getKfactors().

const int PhiSymmetryCalibration::kMaxEndciPhi = 360

staticprivate

Definition at line 142 of file PhiSymmetryCalibration.h.

const float PhiSymmetryCalibration::kMiscalRangeEB = .05

staticprivate

Definition at line 72 of file PhiSymmetryCalibration.h.

Referenced by beginJob(), and getKfactors().

const float PhiSymmetryCalibration::kMiscalRangeEE = .10

staticprivate

Definition at line 75 of file PhiSymmetryCalibration.h.

Referenced by beginJob(), and getKfactors().

const int PhiSymmetryCalibration::kNMiscalBinsEB = 21

staticprivate

Definition at line 71 of file PhiSymmetryCalibration.h.

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

const int PhiSymmetryCalibration::kNMiscalBinsEE = 41

staticprivate

Definition at line 74 of file PhiSymmetryCalibration.h.

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

double PhiSymmetryCalibration::miscalEB_[kNMiscalBinsEB]

private

Definition at line 98 of file PhiSymmetryCalibration.h.

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

double PhiSymmetryCalibration::miscalEE_[kNMiscalBinsEE]

private

Definition at line 99 of file PhiSymmetryCalibration.h.

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

int PhiSymmetryCalibration::nevents_

private

Definition at line 161 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 109 of file PhiSymmetryCalibration.h.

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

private

Definition at line 110 of file PhiSymmetryCalibration.h.

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

private

Definition at line 86 of file PhiSymmetryCalibration.h.

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

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

private

Definition at line 87 of file PhiSymmetryCalibration.h.

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

const std::string PhiSymmetryCalibration::oldcalibfile_

private

Definition at line 147 of file PhiSymmetryCalibration.h.

Referenced by setUp().

EcalIntercalibConstants PhiSymmetryCalibration::oldCalibs_

private

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

Definition at line 150 of file PhiSymmetryCalibration.h.

Referenced by analyze(), and setUp().

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

private

Definition at line 105 of file PhiSymmetryCalibration.h.

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

private

Definition at line 106 of file PhiSymmetryCalibration.h.

float PhiSymmetryCalibration::phi_endc[kMaxEndciPhi][kEndcEtaRings]

private

Definition at line 144 of file PhiSymmetryCalibration.h.

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

private

Definition at line 117 of file PhiSymmetryCalibration.h.

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

private

Definition at line 118 of file PhiSymmetryCalibration.h.

const bool PhiSymmetryCalibration::reiteration_

private

Definition at line 146 of file PhiSymmetryCalibration.h.

Referenced by analyze(), and setUp().

bool PhiSymmetryCalibration::spectra

private

Definition at line 160 of file PhiSymmetryCalibration.h.

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

const int PhiSymmetryCalibration::statusThreshold_

private

threshold in channel status beyond which channel is marked bad

Definition at line 140 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

Member Function Documentation

Member Data Documentation