ElectronCalibrationUniv Class Reference

#include <Calibration/EcalCalibAlgos/src/ElectronCalibrationUniv.cc>

Inheritance diagram for ElectronCalibrationUniv:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	beginRun (edm::Run const &, edm::EventSetup const &) override
	ElectronCalibrationUniv (const edm::ParameterSet &)
void	endJob () override
void	endRun (edm::Run const &, edm::EventSetup const &) override
	~ElectronCalibrationUniv () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const 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
DetId	findMaxHit (const std::vector< DetId > &v1, const EBRecHitCollection *EBhits, const EERecHitCollection *EEhits)
bool	TestEEvalidDetId (int crystal_ix, int crystal_iy, int iz)

Private Attributes
const std::string	calibAlgo_
int	calibClusterSize
TH1F *	calibinter
TH1F *	calibinterEndCapMinus
TH1F *	calibinterEndCapPlus
TH1F *	calibinterNoCuts
TH1F *	calibs
TH1F *	calibsEndCapMinus
TH1F *	calibsEndCapPlus
TH1F *	calibsNoCuts
const int	clusterSize_
const double	cut1_
const double	cut2_
const double	cut3_
const double	cutCalo1_
const double	cutCalo2_
const double	cutEPCalo1_
const double	cutEPCalo2_
const double	cutEPin1_
const double	cutEPin2_
const double	cutESeed_
TH1F *	e25
TH1F *	e25NoCuts
TH1F *	E25oP
TH1F *	E25oPNoCuts
TH2F *	E25oPvsEta
TH2F *	E25oPvsEtaEndCapMinus
TH2F *	E25oPvsEtaEndCapPlus
TH1F *	e25OverScE
TH1F *	e25OverScENoCuts
TH1F *	e9
TH1F *	e9NoCuts
TH1F *	e9Overe25
TH1F *	e9Overe25NoCuts
const edm::InputTag	ebRecHitLabel_
const edm::EDGetTokenT < EBRecHitCollection >	ebRecHitToken_
const edm::InputTag	eeRecHitLabel_
const edm::EDGetTokenT < EERecHitCollection >	eeRecHitToken_
const int	elecclass_
const edm::InputTag	electronLabel_
const double	elePt_
std::vector< float >	EnergyVector
std::vector< float >	EnergyVectorNoCuts
TH1F *	EoP
TH1F *	EoP_all
TH1F *	EoPNoCuts
TH1F *	Error1
TH1F *	Error1NoCuts
TH1F *	Error2
TH1F *	Error2NoCuts
TH1F *	Error3
TH1F *	Error3NoCuts
TH1F *	eSeedOverPout
TH1F *	eSeedOverPout2
TH1F *	eSeedOverPout2ESeed
TH1F *	eSeedOverPout2NoCuts
TH1F *	eSeedOverPoutNoCuts
int	etaMax
int	etaMin
std::vector< std::vector< float > >	EventMatrix
std::vector< std::vector< float > >	EventMatrixNoCuts
TH1F *	EventsAfterCuts
TFile *	f
TH2F *	GeneralMap
TH2F *	GeneralMapBeforePt
TH2F *	GeneralMapEndCapMinus
TH2F *	GeneralMapEndCapMinusBeforePt
TH2F *	GeneralMapEndCapPlus
TH2F *	GeneralMapEndCapPlusBeforePt
const edm::EDGetTokenT < reco::GsfElectronCollection >	gsfElectronToken_
TH1F *	hadOverEm
TH1F *	hadOverEmESeed
TH1F *	hadOverEmNoCuts
int	keventweight
const int	keventweight_
TH2F *	Map
TH2F *	Map3Dcalib
TH2F *	Map3DcalibEndCapMinus
TH2F *	Map3DcalibEndCapPlus
TH2F *	Map3DcalibNoCuts
TH2F *	MapCor1
TH2F *	MapCor10
TH2F *	MapCor10ESeed
TH2F *	MapCor10NoCuts
TH2F *	MapCor11
TH2F *	MapCor11ESeed
TH2F *	MapCor11NoCuts
TH2F *	MapCor1ESeed
TH2F *	MapCor1NoCuts
TH2F *	MapCor2
TH2F *	MapCor2ESeed
TH2F *	MapCor2NoCuts
TH2F *	MapCor3
TH2F *	MapCor3ESeed
TH2F *	MapCor3NoCuts
TH2F *	MapCor4
TH2F *	MapCor4ESeed
TH2F *	MapCor4NoCuts
TH2F *	MapCor5
TH2F *	MapCor5ESeed
TH2F *	MapCor5NoCuts
TH2F *	MapCor6
TH2F *	MapCor6ESeed
TH2F *	MapCor6NoCuts
TH2F *	MapCor7
TH2F *	MapCor7ESeed
TH2F *	MapCor7NoCuts
TH2F *	MapCor8
TH2F *	MapCor8ESeed
TH2F *	MapCor8NoCuts
TH2F *	MapCor9
TH2F *	MapCor9ESeed
TH2F *	MapCor9NoCuts
TH2F *	MapEndCapMinus
TH2F *	MapEndCapPlus
std::vector< int >	MaxCCeta
std::vector< int >	MaxCCetaNoCuts
std::vector< int >	MaxCCphi
std::vector< int >	MaxCCphiNoCuts
const int	maxeta_
const int	maxphi_
const int	mineta_
const int	minphi_
const std::string	miscalibfile_
const std::string	miscalibfileEndCap_
HouseholderDecomposition *	MyHH
MinL3Algorithm *	MyL3Algo1
const int	numevent_
int	phiMax
int	phiMin
TH1F *	PinMinPout
TH1F *	PinMinPoutNoCuts
TH1F *	PinOverPout
TH1F *	PinOverPoutNoCuts
int	read_events
const std::string	rootfile_
TH1F *	scE
TH1F *	scENoCuts
std::vector< float >	solution
std::vector< float >	solutionNoCuts
const CaloTopology *	theCaloTopology_
const edm::ESGetToken < CaloTopology, CaloTopologyRecord >	topologyToken_
TH1F *	trP
TH1F *	trPNoCuts
std::vector< std::vector< DetId > >	UnivEventIds
MinL3AlgoUniv< DetId > *	UnivL3
std::map< DetId, float >	Univsolution
std::vector< float >	WeightVector
std::vector< float >	WeightVectorNoCuts

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 ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: Perform single electron calibration (tested on TB data only).

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

Definition at line 58 of file ElectronCalibrationUniv.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 20 of file ElectronCalibrationUniv.cc.

    : ebRecHitLabel_(iConfig.getParameter<edm::InputTag>("ebRecHitsLabel")),
      eeRecHitLabel_(iConfig.getParameter<edm::InputTag>("eeRecHitsLabel")),
      electronLabel_(iConfig.getParameter<edm::InputTag>("electronLabel")),
      rootfile_(iConfig.getParameter<std::string>("rootfile")),
      calibAlgo_(iConfig.getParameter<std::string>("CALIBRATION_ALGO")),
      miscalibfile_(iConfig.getParameter<std::string>("miscalibfile")),
      miscalibfileEndCap_(iConfig.getParameter<std::string>("miscalibfileEndCap")),
      keventweight_(iConfig.getParameter<int>("keventweight")),
      elePt_(iConfig.getParameter<double>("ElePt")),
      maxeta_(iConfig.getParameter<int>("maxeta")),
      mineta_(iConfig.getParameter<int>("mineta")),
      maxphi_(iConfig.getParameter<int>("maxphi")),
      minphi_(iConfig.getParameter<int>("minphi")),
      cut1_(iConfig.getParameter<double>("cut1")),
      cut2_(iConfig.getParameter<double>("cut2")),
      cut3_(iConfig.getParameter<double>("cut3")),
      numevent_(iConfig.getParameter<int>("numevent")),
      cutEPCalo1_(iConfig.getParameter<double>("cutEPCaloMin")),
      cutEPCalo2_(iConfig.getParameter<double>("cutEPCaloMax")),
      cutEPin1_(iConfig.getParameter<double>("cutEPinMin")),
      cutEPin2_(iConfig.getParameter<double>("cutEPinMax")),
      cutCalo1_(iConfig.getParameter<double>("cutCaloMin")),
      cutCalo2_(iConfig.getParameter<double>("cutCaloMax")),
      cutESeed_(iConfig.getParameter<double>("cutESeed")),
      clusterSize_(iConfig.getParameter<int>("Clustersize")),
      elecclass_(iConfig.getParameter<int>("elecclass")),
      ebRecHitToken_(consumes<EBRecHitCollection>(ebRecHitLabel_)),
      eeRecHitToken_(consumes<EERecHitCollection>(eeRecHitLabel_)),
      gsfElectronToken_(consumes<reco::GsfElectronCollection>(electronLabel_)),
      topologyToken_(esConsumes<edm::Transition::BeginRun>()) {}

ElectronCalibrationUniv::~ElectronCalibrationUniv ( )

override

Definition at line 52 of file ElectronCalibrationUniv.cc.

{}

Member Function Documentation

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

overridevirtual

START HERE....

Implements edm::one::EDAnalyzerBase.

Definition at line 609 of file ElectronCalibrationUniv.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, PixelTestBeamValidation_cfi::ClusterSize, clusterSize_, gather_cfg::cout, cut1_, cut2_, cut3_, cutEPin1_, cutEPin2_, e25, e25NoCuts, E25oP, E25oPNoCuts, E25oPvsEta, E25oPvsEtaEndCapMinus, E25oPvsEtaEndCapPlus, e25OverScE, e25OverScENoCuts, e9, e9NoCuts, e9Overe25, e9Overe25NoCuts, ebRecHitToken_, DetId::Ecal, EcalBarrel, eeRecHitToken_, elecclass_, elePt_, relval_parameters_module::energy, EnergyVector, EoP, EoPNoCuts, Error1, Error1NoCuts, Error2, Error2NoCuts, Error3, Error3NoCuts, eSeedOverPout, eSeedOverPout2, eSeedOverPout2NoCuts, eSeedOverPoutNoCuts, EventMatrix, EventsAfterCuts, findMaxHit(), newFWLiteAna::found, GeneralMap, GeneralMapBeforePt, GeneralMapEndCapMinus, GeneralMapEndCapMinusBeforePt, GeneralMapEndCapPlus, GeneralMapEndCapPlusBeforePt, edm::Event::getByToken(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), gsfElectronToken_, hadOverEm, hadOverEmNoCuts, edm::isNotFinite(), Map, MapCor1, MapCor10, MapCor10NoCuts, MapCor11, MapCor11NoCuts, MapCor1NoCuts, MapCor2, MapCor2NoCuts, MapCor3, MapCor3NoCuts, MapCor4, MapCor4NoCuts, MapCor5, MapCor5NoCuts, MapCor6, MapCor6NoCuts, MapCor7, MapCor7NoCuts, MapCor8, MapCor8NoCuts, MapCor9, MapCor9NoCuts, MapEndCapMinus, MapEndCapPlus, PinMinPout, PinMinPoutNoCuts, PinOverPout, PinOverPoutNoCuts, funct::pow(), read_events, scE, scENoCuts, mathSSE::sqrt(), theCaloTopology_, trP, trPNoCuts, groupFilesInBlocks::tt, UnivEventIds, and trackerHitRTTI::vector.

                                                                                         {
  //=================================================================================
  using namespace edm;

  // Get EBRecHits
  edm::Handle<EBRecHitCollection> EBphits;
  iEvent.getByToken(ebRecHitToken_, EBphits);
  if (!EBphits.isValid()) {
    std::cerr << "Error! can't get the product EBRecHitCollection: " << std::endl;
  }
  const EBRecHitCollection* EBhits = EBphits.product();  // get a ptr to the product

  // Get EERecHits
  edm::Handle<EERecHitCollection> EEphits;

  iEvent.getByToken(eeRecHitToken_, EEphits);
  if (!EEphits.isValid()) {
    std::cerr << "Error! can't get the product EERecHitCollection: " << std::endl;
  }
  const EERecHitCollection* EEhits = EEphits.product();  // get a ptr to the product

  // Get pixelElectrons
  edm::Handle<reco::GsfElectronCollection> pElectrons;
  iEvent.getByToken(gsfElectronToken_, pElectrons);
  if (!pElectrons.isValid()) {
    std::cerr << "Error! can't get the product ElectronCollection: " << std::endl;
  }
  const reco::GsfElectronCollection* electronCollection = pElectrons.product();
  read_events++;
  if (read_events % 1000 == 0)
    std::cout << "read_events = " << read_events << std::endl;

  EventsAfterCuts->Fill(1);
  if (!EBhits || !EEhits)
    return;
  EventsAfterCuts->Fill(2);
  if (EBhits->empty() && EEhits->empty())
    return;
  EventsAfterCuts->Fill(3);
  if (!electronCollection)
    return;
  EventsAfterCuts->Fill(4);
  if (electronCollection->empty())
    return;

  //    ////////////////Need to recalibrate the events (copy code from EcalRecHitRecalib):


  reco::GsfElectronCollection::const_iterator eleIt = electronCollection->begin();

  reco::GsfElectron highPtElectron;

  float highestElePt = 0.;
  bool found = false;
  for (eleIt = electronCollection->begin(); eleIt != electronCollection->end(); eleIt++) {
    if (fabs(eleIt->eta()) > 2.4)
      continue;
    //     if(eleIt->eta()<0.0) continue;

    if (eleIt->pt() > highestElePt) {
      highestElePt = eleIt->pt();
      highPtElectron = *eleIt;
      found = true;
      //       std::cout<<" eleIt->pt( "<<eleIt->pt()<<" eleIt->eta() "<<eleIt->eta()<<std::endl;
    }
  }
  EventsAfterCuts->Fill(5);
  if (!found)
    return;

  const reco::SuperCluster& sc = *(highPtElectron.superCluster());
  //  if(fabs(sc.eta())>1.479){std::cout<<" SC not in Barrel "<<sc.eta()<<std::endl;;}
  //  const std::vector<DetId> & v1 = sc.getHitsByDetId();

  std::vector<DetId> v1;
  //Loop to fill the vector of DetIds
  for (std::vector<std::pair<DetId, float> >::const_iterator idsIt = sc.hitsAndFractions().begin();
       idsIt != sc.hitsAndFractions().end();
       ++idsIt) {
    v1.push_back(idsIt->first);
  }

  DetId maxHitId;

  maxHitId = findMaxHit(v1, (EBhits), (EEhits));
  //maxHitId = findMaxHit(v1,EBhits,EEhits);

  EventsAfterCuts->Fill(6);
  if (maxHitId.null()) {
    std::cout << " Null " << std::endl;
    return;
  }

  int maxCC_Eta = 0;
  int maxCC_Phi = 0;
  int Zside = 0;
  if (maxHitId.subdetId() != 1) {
    maxCC_Eta = ((EEDetId)maxHitId).ix();
    maxCC_Phi = ((EEDetId)maxHitId).iy();
    Zside = ((EEDetId)maxHitId).zside();
    //    std::cout<<" ++++++++ Zside "<<Zside<<std::endl;
  } else {
    maxCC_Eta = ((EBDetId)maxHitId).ieta();
    maxCC_Phi = ((EBDetId)maxHitId).iphi();
  }

  //   if(maxCC_Eta>maxeta_ ) ;
  //   if(maxCC_Eta<mineta_ )  ;

  // number of events per crystal is set
  //   eventcrystal[maxCC_Eta][maxCC_Phi]+=1;
  //   if(eventcrystal[maxCC_Eta][maxCC_Phi] > numevent_) ;

  // fill cluster energy
  std::vector<float> energy;
  float energy3x3 = 0.;
  float energy5x5 = 0.;
  //Should be moved to cfg file!
  int ClusterSize = clusterSize_;

  const CaloSubdetectorTopology* topology = theCaloTopology_->getSubdetectorTopology(DetId::Ecal, maxHitId.subdetId());
  std::vector<DetId> NxNaroundMax = topology->getWindow(maxHitId, ClusterSize, ClusterSize);
  //ToCompute 3x3
  std::vector<DetId> S9aroundMax = topology->getWindow(maxHitId, 3, 3);

  EventsAfterCuts->Fill(7);
  if ((int)NxNaroundMax.size() != ClusterSize * ClusterSize)
    return;
  EventsAfterCuts->Fill(8);
  if (S9aroundMax.size() != 9)
    return;

  //   std::cout<<" ******** New Event "<<std::endl;

  EventsAfterCuts->Fill(9);
  for (int icry = 0; icry < ClusterSize * ClusterSize; icry++) {
    if (NxNaroundMax[icry].subdetId() == EcalBarrel) {
      EBRecHitCollection::const_iterator itrechit;
      itrechit = EBhits->find(NxNaroundMax[icry]);
      if (itrechit == EBhits->end()) {
        //  std::cout << "EB DetId not in e25" << std::endl;
        energy.push_back(0.);
        energy5x5 += 0.;
        continue;
      }

      if (edm::isNotFinite(itrechit->energy())) {
        std::cout << " nan energy " << std::endl;
        return;
      }
      energy.push_back(itrechit->energy());
      energy5x5 += itrechit->energy();

      //Summing in 3x3 to cut later on:
      for (int tt = 0; tt < 9; tt++) {
        if (NxNaroundMax[icry] == S9aroundMax[tt])
          energy3x3 += itrechit->energy();
      }
    } else {
      EERecHitCollection::const_iterator itrechit;

      itrechit = EEhits->find(NxNaroundMax[icry]);

      if (itrechit == EEhits->end()) {
        //  std::cout << "EE DetId not in e25" << std::endl;
        //  std::cout<<" ******** putting 0 "<<std::endl;
        energy.push_back(0.);
        energy5x5 += 0.;
        continue;
      }

      if (edm::isNotFinite(itrechit->energy())) {
        std::cout << " nan energy " << std::endl;
        return;
      }
      energy.push_back(itrechit->energy());
      energy5x5 += itrechit->energy();

      //Summing in 3x3 to cut later on:
      for (int tt = 0; tt < 9; tt++) {
        if (NxNaroundMax[icry] == S9aroundMax[tt])
          energy3x3 += itrechit->energy();
      }
    }
  }
  //  if((read_events-50)%10000 ==0)cout << "++++++++++++ENERGY 5x5 " <<  energy5x5 << std::endl;
  EventsAfterCuts->Fill(10);
  //  std::cout<<" ******** NxNaroundMax.size() "<<NxNaroundMax.size()<<std::endl;
  //  std::cout<<" ******** energy.size() "<<energy.size()<<std::endl;
  if ((int)energy.size() != ClusterSize * ClusterSize)
    return;

  if (maxHitId.subdetId() == EcalBarrel) {
    GeneralMapBeforePt->Fill(maxCC_Eta, maxCC_Phi);
  } else {
    if (Zside < 0) {
      GeneralMapEndCapMinusBeforePt->Fill(maxCC_Eta, maxCC_Phi);
    } else {
      GeneralMapEndCapPlusBeforePt->Fill(maxCC_Eta, maxCC_Phi);
    }
  }

  EventsAfterCuts->Fill(11);
  if (highestElePt < elePt_)
    return;

  if (maxHitId.subdetId() == EcalBarrel) {
    GeneralMap->Fill(maxCC_Eta, maxCC_Phi);
  } else {
    if (Zside < 0) {
      GeneralMapEndCapMinus->Fill(maxCC_Eta, maxCC_Phi);
    } else {
      GeneralMapEndCapPlus->Fill(maxCC_Eta, maxCC_Phi);
    }
  }

  EventsAfterCuts->Fill(12);
  if (highPtElectron.classification() != elecclass_ && elecclass_ != -1)
    return;

  float Ptrack_in =
      sqrt(pow(highPtElectron.trackMomentumAtVtx().X(), 2) + pow(highPtElectron.trackMomentumAtVtx().Y(), 2) +
           pow(highPtElectron.trackMomentumAtVtx().Z(), 2));

  float UncorrectedPatCalo =
      sqrt(pow(highPtElectron.trackMomentumAtCalo().X(), 2) + pow(highPtElectron.trackMomentumAtCalo().Y(), 2) +
           pow(highPtElectron.trackMomentumAtCalo().Z(), 2));

  float Ptrack_out =
      sqrt(pow(highPtElectron.trackMomentumOut().X(), 2) + pow(highPtElectron.trackMomentumOut().Y(), 2) +
           pow(highPtElectron.trackMomentumOut().Z(), 2));

  e9NoCuts->Fill(energy3x3);
  e25NoCuts->Fill(energy5x5);
  e9Overe25NoCuts->Fill(energy3x3 / energy5x5);
  scENoCuts->Fill(sc.energy());

  trPNoCuts->Fill(UncorrectedPatCalo);

  EoPNoCuts->Fill(highPtElectron.eSuperClusterOverP());
  e25OverScENoCuts->Fill(energy5x5 / sc.energy());

  E25oPNoCuts->Fill(energy5x5 / UncorrectedPatCalo);

  PinOverPoutNoCuts->Fill(
      sqrt(pow(highPtElectron.trackMomentumAtVtx().X(), 2) + pow(highPtElectron.trackMomentumAtVtx().Y(), 2) +
           pow(highPtElectron.trackMomentumAtVtx().Z(), 2)) /
      sqrt(pow(highPtElectron.trackMomentumOut().X(), 2) + pow(highPtElectron.trackMomentumOut().Y(), 2) +
           pow(highPtElectron.trackMomentumOut().Z(), 2)));
  eSeedOverPoutNoCuts->Fill(highPtElectron.eSuperClusterOverP());

  MapCor1NoCuts->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / Ptrack_in);
  MapCor2NoCuts->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSuperClusterOverP());
  MapCor3NoCuts->Fill(energy5x5 / UncorrectedPatCalo, Ptrack_out / Ptrack_in);
  MapCor4NoCuts->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / highPtElectron.p());
  MapCor5NoCuts->Fill(energy5x5 / UncorrectedPatCalo, UncorrectedPatCalo / Ptrack_out);
  MapCor6NoCuts->Fill(Ptrack_out / Ptrack_in, energy5x5 / Ptrack_in);
  MapCor7NoCuts->Fill(Ptrack_out / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
  MapCor8NoCuts->Fill(energy5x5 / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
  MapCor9NoCuts->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSeedClusterOverPout());
  MapCor10NoCuts->Fill(highPtElectron.eSeedClusterOverPout(), Ptrack_out / Ptrack_in);
  MapCor11NoCuts->Fill(highPtElectron.eSeedClusterOverPout(), energy5x5 / Ptrack_in);

  PinMinPoutNoCuts->Fill((Ptrack_in - Ptrack_out) / Ptrack_in);

  Error1NoCuts->Fill(highPtElectron.trackMomentumError() / Ptrack_in);
  Error2NoCuts->Fill(highPtElectron.trackMomentumError() / Ptrack_out);
  Error3NoCuts->Fill(highPtElectron.trackMomentumError() / UncorrectedPatCalo);
  eSeedOverPout2NoCuts->Fill(highPtElectron.eSeedClusterOverPout());

  hadOverEmNoCuts->Fill(highPtElectron.hadronicOverEm());

  //Cuts!
  if ((energy3x3 / energy5x5) < cut1_)
    return;
  if ((Ptrack_out / Ptrack_in) < cut2_ || (Ptrack_out / Ptrack_in) > cut3_)
    return;
  if ((energy5x5 / Ptrack_in) < cutEPin1_ || (energy5x5 / Ptrack_in) > cutEPin2_)
    return;
  //    if(!highPtElectron.ecalDriven())return;
  //    if(!highPtElectron.passingCutBasedPreselection())return;

  // //  Apply Pietro cuts:
  //    EventsAfterCuts->Fill(13);
  //    //Module 1
  //    if(maxHitId.subdetId() == EcalBarrel){
  //      //Module 1
  //      if(maxCC_Eta <= 25){
  //        if(highPtElectron.eSuperClusterOverP()>1.05 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //        if(highPtElectron.eSeedClusterOverPout()>1.4 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //        if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.2)return ;
  //      }else{
  //        //Module 2
  //        if( maxCC_Eta > 25&& maxCC_Eta <= 45){
  //          if(highPtElectron.eSuperClusterOverP()>1.05 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //          if(highPtElectron.eSeedClusterOverPout()>1.25 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //          if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.2)return ;
  //        }else{
  //        //Module 3
  //          if( maxCC_Eta > 45&& maxCC_Eta <= 65){
  //        if(highPtElectron.eSuperClusterOverP()>1.05 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //        if(highPtElectron.eSeedClusterOverPout()>1.15 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //        if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.15)return ;
  //          }else{
  //          if( maxCC_Eta > 65&& maxCC_Eta <= 85){
  //        if(highPtElectron.eSuperClusterOverP()>1.05 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //        if(highPtElectron.eSeedClusterOverPout()>1.15 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //        if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.15)return ;
  //          }else{
  //        return;
  //          }
  //          }
  //        }
  //      }
  //    }else{
  //      //EndCapMinus Side:
  //      //EndCapPlus Side:
  //      int iR = sqrt((maxCC_Eta-50)*(maxCC_Eta-50) + (maxCC_Phi-50)*(maxCC_Phi-50));
  //      if( iR >= 22&& iR < 27){
  //        if(highPtElectron.eSuperClusterOverP()>1.05 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //        if(highPtElectron.eSeedClusterOverPout()>1.15 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //        if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.2)return ;
  //      }else{
  //        if( iR >= 27&& iR < 32){
  //          if(highPtElectron.eSuperClusterOverP()>1.1 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //          if(highPtElectron.eSeedClusterOverPout()>1.25 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //          if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.2)return ;
  //        }else{
  //          if( iR >= 32&& iR < 37){
  //        if(highPtElectron.eSuperClusterOverP()>1.05 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //        if(highPtElectron.eSeedClusterOverPout()>1.15 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //        if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.2)return ;
  //          }else{
  //        if( iR >= 37&& iR < 42){
  //          if(highPtElectron.eSuperClusterOverP()>1.1 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //          if(highPtElectron.eSeedClusterOverPout()>1.15 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //          if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.15)return ;
  //        }else{
  //          if( iR >= 42){
  //            if(highPtElectron.eSuperClusterOverP()>1.05 || highPtElectron.eSuperClusterOverP()<0.95)return ;
  //            if(highPtElectron.eSeedClusterOverPout()>1.15 || highPtElectron.eSeedClusterOverPout()<0.90)return ;
  //          if((Ptrack_in- Ptrack_out) / Ptrack_in <-0.05 || (Ptrack_in- Ptrack_out) / Ptrack_in >0.15)return ;
  //          }
  //        }
  //          }
  //        }
  //      }
  //    }

  if (maxHitId.subdetId() == EcalBarrel) {
    E25oPvsEta->Fill(maxCC_Eta, energy5x5 / UncorrectedPatCalo);
  } else {
    float Radius = sqrt((maxCC_Eta) * (maxCC_Eta) + (maxCC_Phi) * (maxCC_Phi));
    if (Zside < 0) {
      E25oPvsEtaEndCapMinus->Fill(Radius, energy5x5 / UncorrectedPatCalo);
    } else {
      E25oPvsEtaEndCapPlus->Fill(Radius, energy5x5 / UncorrectedPatCalo);
    }
  }
  e9->Fill(energy3x3);
  e25->Fill(energy5x5);
  e9Overe25->Fill(energy3x3 / energy5x5);
  scE->Fill(sc.energy());
  trP->Fill(UncorrectedPatCalo);

  EoP->Fill(highPtElectron.eSuperClusterOverP());
  e25OverScE->Fill(energy5x5 / sc.energy());

  E25oP->Fill(energy5x5 / UncorrectedPatCalo);

  if (maxHitId.subdetId() == EcalBarrel) {
    Map->Fill(maxCC_Eta, maxCC_Phi);
  } else {
    if (Zside < 0) {
      MapEndCapMinus->Fill(maxCC_Eta, maxCC_Phi);
    } else {
      MapEndCapPlus->Fill(maxCC_Eta, maxCC_Phi);
    }
  }

  PinOverPout->Fill(sqrt(pow(highPtElectron.trackMomentumAtVtx().X(), 2) +
                         pow(highPtElectron.trackMomentumAtVtx().Y(), 2) +
                         pow(highPtElectron.trackMomentumAtVtx().Z(), 2)) /
                    sqrt(pow(highPtElectron.trackMomentumOut().X(), 2) + pow(highPtElectron.trackMomentumOut().Y(), 2) +
                         pow(highPtElectron.trackMomentumOut().Z(), 2)));
  eSeedOverPout->Fill(highPtElectron.eSuperClusterOverP());

  MapCor1->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / Ptrack_in);
  MapCor2->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSuperClusterOverP());
  MapCor3->Fill(energy5x5 / UncorrectedPatCalo, Ptrack_out / Ptrack_in);
  MapCor4->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / highPtElectron.p());
  MapCor5->Fill(energy5x5 / UncorrectedPatCalo, UncorrectedPatCalo / Ptrack_out);
  MapCor6->Fill(Ptrack_out / Ptrack_in, energy5x5 / Ptrack_in);
  MapCor7->Fill(Ptrack_out / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
  MapCor8->Fill(energy5x5 / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
  MapCor9->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSeedClusterOverPout());
  MapCor10->Fill(highPtElectron.eSeedClusterOverPout(), Ptrack_out / Ptrack_in);
  MapCor11->Fill(highPtElectron.eSeedClusterOverPout(), energy5x5 / Ptrack_in);

  PinMinPout->Fill((Ptrack_in - Ptrack_out) / Ptrack_in);

  Error1->Fill(highPtElectron.trackMomentumError() / Ptrack_in);
  Error2->Fill(highPtElectron.trackMomentumError() / Ptrack_out);
  Error3->Fill(highPtElectron.trackMomentumError() / UncorrectedPatCalo);

  eSeedOverPout2->Fill(highPtElectron.eSeedClusterOverPout());
  hadOverEm->Fill(highPtElectron.hadronicOverEm());

  UnivEventIds.push_back(NxNaroundMax);
  EventMatrix.push_back(energy);
  EnergyVector.push_back(UncorrectedPatCalo);

  EventsAfterCuts->Fill(14);

  if (!highPtElectron.ecalDrivenSeed())
    EventsAfterCuts->Fill(15);

  return;
}

void ElectronCalibrationUniv::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 55 of file ElectronCalibrationUniv.cc.

References calibAlgo_, calibClusterSize, calibinter, calibinterEndCapMinus, calibinterEndCapPlus, calibinterNoCuts, calibs, calibsEndCapMinus, calibsEndCapPlus, calibsNoCuts, clusterSize_, gather_cfg::cout, e25, e25NoCuts, E25oP, E25oPNoCuts, E25oPvsEta, E25oPvsEtaEndCapMinus, E25oPvsEtaEndCapPlus, e25OverScE, e25OverScENoCuts, e9, e9NoCuts, e9Overe25, e9Overe25NoCuts, EoP, EoP_all, EoPNoCuts, Error1, Error1NoCuts, Error2, Error2NoCuts, Error3, Error3NoCuts, eSeedOverPout, eSeedOverPout2, eSeedOverPout2ESeed, eSeedOverPout2NoCuts, eSeedOverPoutNoCuts, etaMax, etaMin, EventsAfterCuts, f, GeneralMap, GeneralMapBeforePt, GeneralMapEndCapMinus, GeneralMapEndCapMinusBeforePt, GeneralMapEndCapPlus, GeneralMapEndCapPlusBeforePt, hadOverEm, hadOverEmESeed, hadOverEmNoCuts, keventweight_, Map, Map3Dcalib, Map3DcalibEndCapMinus, Map3DcalibEndCapPlus, Map3DcalibNoCuts, MapCor1, MapCor10, MapCor10ESeed, MapCor10NoCuts, MapCor11, MapCor11ESeed, MapCor11NoCuts, MapCor1ESeed, MapCor1NoCuts, MapCor2, MapCor2ESeed, MapCor2NoCuts, MapCor3, MapCor3ESeed, MapCor3NoCuts, MapCor4, MapCor4ESeed, MapCor4NoCuts, MapCor5, MapCor5ESeed, MapCor5NoCuts, MapCor6, MapCor6ESeed, MapCor6NoCuts, MapCor7, MapCor7ESeed, MapCor7NoCuts, MapCor8, MapCor8ESeed, MapCor8NoCuts, MapCor9, MapCor9ESeed, MapCor9NoCuts, MapEndCapMinus, MapEndCapPlus, maxeta_, maxphi_, mineta_, minphi_, MyHH, MyL3Algo1, phiMax, phiMin, PinMinPout, PinMinPoutNoCuts, PinOverPout, PinOverPoutNoCuts, read_events, rootfile_, scE, scENoCuts, trP, trPNoCuts, and UnivL3.

                                       {
  //========================================================================
  f = new TFile(rootfile_.c_str(), "RECREATE");
  f->cd();
  EventsAfterCuts = new TH1F("EventsAfterCuts", "Events After Cuts", 30, 0, 30);

  // Book histograms
  e9 = new TH1F("e9", "E9 energy", 300, 0., 150.);
  e25 = new TH1F("e25", "E25 energy", 300, 0., 150.);
  scE = new TH1F("scE", "SC energy", 300, 0., 150.);
  trP = new TH1F("trP", "Trk momentum", 300, 0., 150.);
  EoP = new TH1F("EoP", "EoP", 600, 0., 3.);
  EoP_all = new TH1F("EoP_all", "EoP_all", 600, 0., 3.);

  calibs = new TH1F("calib", "Calibration constants", 800, 0.5, 2.);
  calibsEndCapMinus = new TH1F("calibEndCapMinus", "Calibration constants EE-", 800, 0.5, 2.);
  calibsEndCapPlus = new TH1F("calibEndCapPlus", "Calibration constants EE+", 800, 0.5, 2.);

  e25OverScE = new TH1F("e25OverscE", "E25 / SC energy", 400, 0., 2.);
  E25oP = new TH1F("E25oP", "E25 / P", 750, 0., 1.5);

  Map = new TH2F("Map", "Nb Events in Crystal", 173, -86, 86, 362, 0, 361);
  e9Overe25 = new TH1F("e9Overe25", "E9 / E25", 400, 0., 2.);
  Map3Dcalib = new TH2F("3Dcalib", "3Dcalib", 173, -86, 86, 362, 0, 361);
  Map3DcalibEndCapMinus = new TH2F("3DcalibEndCapMinus", "3Dcalib EE-", 100, 0, 100, 100, 0, 100);
  Map3DcalibEndCapPlus = new TH2F("3DcalibEndCapPlus", "3Dcalib EE+", 100, 0, 100, 100, 0, 100);

  MapCor1 = new TH2F("MapCor1", "Correlation E25/Pcalo versus E25/Pin", 100, 0., 5., 100, 0., 5.);
  MapCor2 = new TH2F("MapCor2", "Correlation E25/Pcalo versus E/P", 100, 0., 5., 100, 0., 5.);
  MapCor3 = new TH2F("MapCor3", "Correlation E25/Pcalo versus Pout/Pin", 100, 0., 5., 100, 0., 5.);
  MapCor4 = new TH2F("MapCor4", "Correlation E25/Pcalo versus E25/highestP", 100, 0., 5., 100, 0., 5.);
  MapCor5 = new TH2F("MapCor5", "Correlation E25/Pcalo versus Pcalo/Pout", 100, 0., 5., 100, 0., 5.);
  MapCor6 = new TH2F("MapCor6", "Correlation Pout/Pin versus E25/Pin", 100, 0., 5., 100, 0., 5.);
  MapCor7 = new TH2F("MapCor7", "Correlation Pout/Pin versus Pcalo/Pout", 100, 0., 5., 100, 0., 5.);
  MapCor8 = new TH2F("MapCor8", "Correlation E25/Pin versus Pcalo/Pout", 100, 0., 5., 100, 0., 5.);
  MapCor9 = new TH2F("MapCor9", "Correlation  E25/Pcalo versus Eseed/Pout", 100, 0., 5., 100, 0., 5.);
  MapCor10 = new TH2F("MapCor10", "Correlation Eseed/Pout versus Pout/Pin", 100, 0., 5., 100, 0., 5.);
  MapCor11 = new TH2F("MapCor11", "Correlation Eseed/Pout versus E25/Pin", 100, 0., 5., 100, 0., 5.);
  //   MapCorCalib = new TH2F ("MapCorCalib", "Correlation Miscalibration versus Calibration constants", 500, 0.5,1.5, 500, 0.5, 1.5);

  E25oPvsEta = new TH2F("E25oPvsEta", "E/P vs Eta", 173, -86, 86, 600, 0.7, 1.3);
  E25oPvsEtaEndCapMinus = new TH2F("E25oPvsEtaEndCapMinus", "E/P vs R EE-", 100, 0, 100, 600, 0.7, 1.3);
  E25oPvsEtaEndCapPlus = new TH2F("E25oPvsEtaEndCapPlus", "E/P vs R EE+", 100, 0, 100, 600, 0.7, 1.3);

  PinMinPout = new TH1F("PinMinPout", "(Pin - Pout)/Pin", 600, -2.0, 2.0);

  calibinter = new TH1F("calibinter", "internal calibration constants", 800, 0.5, 2.);
  PinOverPout = new TH1F("PinOverPout", "pinOverpout", 600, 0., 3.);
  eSeedOverPout = new TH1F("eSeedOverPout", "eSeedOverpout ", 600, 0., 3.);
  //   MisCalibs = new TH1F("MisCalibs","Miscalibration constants",800,0.5,2.);
  //   RatioCalibs = new TH1F("RatioCalibs","Ratio in Calibration Constants", 800, 0.5, 2.0);
  //   DiffCalibs = new TH1F("DiffCalibs", "Difference in Calibration constants", 800, -1.0,1.0);
  calibinterEndCapMinus = new TH1F("calibinterEndCapMinus", "internal calibration constants", 800, 0.5, 2.);
  calibinterEndCapPlus = new TH1F("calibinterEndCapPlus", "internal calibration constants", 800, 0.5, 2.);
  //   MisCalibsEndCapMinus = new TH1F("MisCalibsEndCapMinus","Miscalibration constants",800,0.5,2.);
  //   MisCalibsEndCapPlus = new TH1F("MisCalibsEndCapPlus","Miscalibration constants",800,0.5,2.);
  //   RatioCalibsEndCapMinus = new TH1F("RatioCalibsEndCapMinus","Ratio in Calibration Constants", 800, 0.5, 2.0);
  //   RatioCalibsEndCapPlus = new TH1F("RatioCalibsEndCapPlus","Ratio in Calibration Constants", 800, 0.5, 2.0);
  //   DiffCalibsEndCapMinus = new TH1F("DiffCalibsEndCapMinus", "Difference in Calibration constants", 800, -1.0,1.0);
  //   DiffCalibsEndCapPlus = new TH1F("DiffCalibsEndCapPlus", "Difference in Calibration constants", 800, -1.0,1.0);
  Error1 = new TH1F("Error1", "DeltaP/Pin", 800, -1.0, 1.0);
  Error2 = new TH1F("Error2", "DeltaP/Pout", 800, -1.0, 1.0);
  Error3 = new TH1F("Error3", "DeltaP/Pcalo", 800, -1.0, 1.0);
  eSeedOverPout2 = new TH1F("eSeedOverPout2", "eSeedOverpout (No Supercluster)", 600, 0., 4.);
  hadOverEm = new TH1F("hadOverEm", "Had/EM distribution", 600, -2., 2.);

  // Book histograms
  Map3DcalibNoCuts = new TH2F("3DcalibNoCuts", "3Dcalib (Before Cuts)", 173, -86, 86, 362, 0, 361);
  e9NoCuts = new TH1F("e9NoCuts", "E9 energy (Before Cuts)", 300, 0., 150.);
  e25NoCuts = new TH1F("e25NoCuts", "E25 energy (Before Cuts)", 300, 0., 150.);
  scENoCuts = new TH1F("scENoCuts", "SC energy (Before Cuts)", 300, 0., 150.);
  trPNoCuts = new TH1F("trPNoCuts", "Trk momentum (Before Cuts)", 300, 0., 150.);
  EoPNoCuts = new TH1F("EoPNoCuts", "EoP (Before Cuts)", 600, 0., 3.);
  calibsNoCuts = new TH1F("calibNoCuts", "Calibration constants (Before Cuts)", 800, 0., 2.);
  e25OverScENoCuts = new TH1F("e25OverscENoCuts", "E25 / SC energy (Before Cuts)", 400, 0., 2.);
  E25oPNoCuts = new TH1F("E25oPNoCuts", "E25 / P (Before Cuts)", 750, 0., 1.5);
  MapEndCapMinus = new TH2F("MapEndCapMinus", "Nb Events in Crystal (EndCap)", 100, 0, 100, 100, 0, 100);
  MapEndCapPlus = new TH2F("MapEndCapPlus", "Nb Events in Crystal (EndCap)", 100, 0, 100, 100, 0, 100);
  e9Overe25NoCuts = new TH1F("e9Overe25NoCuts", "E9 / E25 (Before Cuts)", 400, 0., 2.);
  PinOverPoutNoCuts = new TH1F("PinOverPoutNoCuts", "pinOverpout (Before Cuts)", 600, 0., 3.);
  eSeedOverPoutNoCuts = new TH1F(" eSeedOverPoutNoCuts", "eSeedOverpout (Before Cuts) ", 600, 0., 4.);
  PinMinPoutNoCuts = new TH1F("PinMinPoutNoCuts", "(Pin - Pout)/Pin (Before Cuts)", 600, -2.0, 2.0);

  //   RatioCalibsNoCuts = new TH1F("RatioCalibsNoCuts","Ratio in Calibration Constants (Before Cuts)", 800, 0.5, 2.0);
  //   DiffCalibsNoCuts = new TH1F("DiffCalibsNoCuts", "Difference in Calibration constants (Before Cuts)", 800, -1.0,1.0);
  calibinterNoCuts = new TH1F("calibinterNoCuts", "internal calibration constants", 2000, 0.5, 2.);

  MapCor1NoCuts =
      new TH2F("MapCor1NoCuts", "Correlation E25/PatCalo versus E25/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor2NoCuts =
      new TH2F("MapCor2NoCuts", "Correlation E25/PatCalo versus E/P (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor3NoCuts =
      new TH2F("MapCor3NoCuts", "Correlation E25/PatCalo versus Pout/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor4NoCuts =
      new TH2F("MapCor4NoCuts", "Correlation E25/PatCalo versus E25/highestP (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor5NoCuts =
      new TH2F("MapCor5NoCuts", "Correlation E25/Pcalo versus Pcalo/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor6NoCuts =
      new TH2F("MapCor6NoCuts", "Correlation Pout/Pin versus E25/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor7NoCuts =
      new TH2F("MapCor7NoCuts", "Correlation Pout/Pin versus Pcalo/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor8NoCuts =
      new TH2F("MapCor8NoCuts", "Correlation E25/Pin versus Pcalo/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor9NoCuts =
      new TH2F("MapCor9NoCuts", "Correlation  E25/Pcalo versus Eseed/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor10NoCuts =
      new TH2F("MapCor10NoCuts", "Correlation Eseed/Pout versus Pout/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  MapCor11NoCuts =
      new TH2F("MapCor11NoCuts", "Correlation Eseed/Pout versus E25/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
  //   MapCorCalibEndCapMinus = new TH2F ("MapCorCalibEndCapMinus", "Correlation Miscalibration versus Calibration constants (EndCap)",  500, 0.5,1.5, 500, 0.5, 1.5);
  //   MapCorCalibEndCapPlus = new TH2F ("MapCorCalibEndCapPlus", "Correlation Miscalibration versus Calibration constants (EndCap)",  500, 0.5,1.5, 500, 0.5, 1.5);

  Error1NoCuts = new TH1F("Eror1NoCuts", "DeltaP/Pin (Before Cuts)", 800, -1.0, 1.0);
  Error2NoCuts = new TH1F("Error2NoCuts", "DeltaP/Pout (Before Cuts)", 800, -1.0, 1.0);
  Error3NoCuts = new TH1F("Error3NoCuts", "DeltaP/Pcalo (Before Cuts)", 800, -1.0, 1.0);
  eSeedOverPout2NoCuts = new TH1F("eSeedOverPout2NoCuts", "eSeedOverpout (No Supercluster, Before Cuts)", 600, 0., 4.);
  hadOverEmNoCuts = new TH1F("hadOverEmNoCuts", "Had/EM distribution (Before Cuts)", 600, -2., 2.);

  //Book histograms after ESeed cut
  MapCor1ESeed =
      new TH2F("MapCor1ESeed", "Correlation E25/Pcalo versus E25/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor2ESeed =
      new TH2F("MapCor2ESeed", "Correlation E25/Pcalo versus E/P (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor3ESeed = new TH2F(
      "MapCor3ESeed", "Correlation E25/Pcalo versus Pout/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor4ESeed = new TH2F(
      "MapCor4ESeed", "Correlation E25/Pcalo versus E25/highestP (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor5ESeed = new TH2F(
      "MapCor5ESeed", "Correlation E25/Pcalo versus Pcalo/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor6ESeed =
      new TH2F("MapCor6ESeed", "Correlation Pout/Pin versus E25/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor7ESeed = new TH2F(
      "MapCor7ESeed", "Correlation Pout/Pin versus Pcalo/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor8ESeed = new TH2F(
      "MapCor8ESeed", "Correlation E25/Pin versus Pcalo/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor9ESeed = new TH2F(
      "MapCor9ESeed", "Correlation  E25/Pcalo versus Eseed/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor10ESeed = new TH2F(
      "MapCor10ESeed", "Correlation Eseed/Pout versus Pout/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
  MapCor11ESeed = new TH2F(
      "MapCor11ESeed", "Correlation Eseed/Pout versus E25/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);

  eSeedOverPout2ESeed =
      new TH1F("eSeedOverPout2ESeed", "eSeedOverpout (No Supercluster, after Eseed/Pout cut)", 600, 0., 4.);

  hadOverEmESeed = new TH1F("hadOverEmESeed", "Had/EM distribution (after Eseed/Pout cut)", 600, -2., 2.);

  //Book histograms without any cut
  GeneralMap = new TH2F("GeneralMap", "Map without any cuts", 173, -86, 86, 362, 0, 361);
  GeneralMapEndCapMinus = new TH2F("GeneralMapEndCapMinus", "Map without any cuts", 100, 0, 100, 100, 0, 100);
  GeneralMapEndCapPlus = new TH2F("GeneralMapEndCapPlus", "Map without any cuts", 100, 0, 100, 100, 0, 100);
  GeneralMapBeforePt = new TH2F("GeneralMapBeforePt", "Map without any cuts", 173, -86, 86, 362, 0, 361);
  GeneralMapEndCapMinusBeforePt =
      new TH2F("GeneralMapEndCapMinusBeforePt", "Map without any cuts", 100, 0, 100, 100, 0, 100);
  GeneralMapEndCapPlusBeforePt =
      new TH2F("GeneralMapEndCapPlusBeforePt", "Map without any cuts", 100, 0, 100, 100, 0, 100);

  calibClusterSize = clusterSize_;
  etaMin = int(mineta_);
  etaMax = int(maxeta_);
  phiMin = int(minphi_);
  phiMax = int(maxphi_);
  if (calibAlgo_ == "L3") {
    MyL3Algo1 = new MinL3Algorithm(keventweight_, calibClusterSize, etaMin, etaMax, phiMin, phiMax);
  } else {
    if (calibAlgo_ == "L3Univ") {
      UnivL3 = new MinL3AlgoUniv<DetId>(keventweight_);
    } else {
      if (calibAlgo_ == "HH" || calibAlgo_ == "HHReg") {
        MyHH = new HouseholderDecomposition(calibClusterSize, etaMin, etaMax, phiMin, phiMax);
      } else {
        std::cout << " Name of Algorithm is not recognize " << calibAlgo_
                  << " Should be either L3, HH or HHReg. Abort! " << std::endl;
      }
    }
  }
  read_events = 0;
}

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

override

Definition at line 235 of file ElectronCalibrationUniv.cc.

References edm::EventSetup::getData(), theCaloTopology_, and topologyToken_.

                                                                                 {
  //========================================================================

  //To Deal with Geometry:
  theCaloTopology_ = &iSetup.getData(topologyToken_);
}

void ElectronCalibrationUniv::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 248 of file ElectronCalibrationUniv.cc.

References calibAlgo_, calibinter, calibinterEndCapMinus, calibinterEndCapPlus, calibinterNoCuts, calibs, calibsEndCapMinus, calibsEndCapPlus, calibsNoCuts, gather_cfg::cout, e25, e25NoCuts, E25oP, E25oPNoCuts, E25oPvsEta, E25oPvsEtaEndCapMinus, E25oPvsEtaEndCapPlus, e25OverScE, e25OverScENoCuts, e9, e9NoCuts, e9Overe25, e9Overe25NoCuts, dataset::end, EnergyVector, EoP, EoP_all, EoPNoCuts, Error1, Error1NoCuts, Error2, Error2NoCuts, Error3, Error3NoCuts, eSeedOverPout, eSeedOverPout2, eSeedOverPout2ESeed, eSeedOverPout2NoCuts, eSeedOverPoutNoCuts, EventMatrix, EventsAfterCuts, f, GeneralMap, GeneralMapBeforePt, GeneralMapEndCapMinus, GeneralMapEndCapMinusBeforePt, GeneralMapEndCapPlus, GeneralMapEndCapPlusBeforePt, hadOverEm, hadOverEmESeed, hadOverEmNoCuts, EBDetId::ieta(), EBDetId::iphi(), MinL3Algorithm::iterate(), MinL3AlgoUniv< IDdet >::iterate(), HouseholderDecomposition::iterate(), EEDetId::ix(), EEDetId::iy(), Map, Map3Dcalib, Map3DcalibEndCapMinus, Map3DcalibEndCapPlus, Map3DcalibNoCuts, MapCor1, MapCor10, MapCor10ESeed, MapCor10NoCuts, MapCor11, MapCor11ESeed, MapCor11NoCuts, MapCor1ESeed, MapCor1NoCuts, MapCor2, MapCor2ESeed, MapCor2NoCuts, MapCor3, MapCor3ESeed, MapCor3NoCuts, MapCor4, MapCor4ESeed, MapCor4NoCuts, MapCor5, MapCor5ESeed, MapCor5NoCuts, MapCor6, MapCor6ESeed, MapCor6NoCuts, MapCor7, MapCor7ESeed, MapCor7NoCuts, MapCor8, MapCor8ESeed, MapCor8NoCuts, MapCor9, MapCor9ESeed, MapCor9NoCuts, MapEndCapMinus, MapEndCapPlus, MaxCCeta, MaxCCphi, maxeta_, maxphi_, mineta_, minphi_, MyHH, MyL3Algo1, PinMinPout, PinMinPoutNoCuts, PinOverPout, PinOverPoutNoCuts, read_events, HouseholderDecomposition::runRegional(), scE, scENoCuts, solution, command_line::start, DetId::subdetId(), trP, trPNoCuts, UnivEventIds, UnivL3, Univsolution, calibXMLwriter::writeLine(), and EEDetId::zside().

Referenced by o2olib.O2ORunMgr::executeJob().

                                     {
  //========================================================================

  f->cd();
  time_t start, end;
  time_t cpu_time_used;
  start = time(nullptr);

  //In order to do only one loop to use properly looper properties, ask only for 1 iterations!
  int nIterations = 10;
  if (calibAlgo_ == "L3") {
    solution = MyL3Algo1->iterate(EventMatrix, MaxCCeta, MaxCCphi, EnergyVector, nIterations);
  } else {
    if (calibAlgo_ == "L3Univ") {
      //Univsolution= UnivL3->getSolution();
      //     std::cout<<" Should derive solution "<<EnergyVector.size()<<std::endl;
      Univsolution = UnivL3->iterate(EventMatrix, UnivEventIds, EnergyVector, nIterations);
      //std::cout<<" solution size "<<Univsolution.size()<<std::endl;
    } else {
      if (calibAlgo_ == "HH") {
        solution = MyHH->iterate(EventMatrix, MaxCCeta, MaxCCphi, EnergyVector, 1, false);
      } else {
        if (calibAlgo_ == "HHReg") {
          solution = MyHH->runRegional(EventMatrix, MaxCCeta, MaxCCphi, EnergyVector, 2);
        } else {
          std::cout << " Calibration not run due to problem in Algo Choice..." << std::endl;
          return;
        }
      }
    }
  }
  end = time(nullptr);
  cpu_time_used = end - start;
  //     std::cout<<"222 solution size "<<Univsolution.size()<<std::endl;

  calibXMLwriter write_calibrations;

  //   FILE* MisCalib;
  //   //char* calibfile="miscalibfile";
  //   MisCalib = fopen(miscalibfile_.c_str(),"r");

  //   int fileStatus=0;
  //   int eta=-1;
  //   int phi=-1;
  //   float coeff=-1;

  std::map<EBDetId, float> OldCoeff;

  //  while(fileStatus != EOF) {
  //    fileStatus = fscanf(MisCalib,"%d %d %f\n",  &eta,&phi,&coeff);
  //    if(eta!=-1&&phi!=-1&& coeff!=-1){
  //      //     std::cout<<" We have read correctly the coefficient " << coeff << " corresponding to eta "<<eta<<" and  phi "<<phi<<std::endl;
  //      OldCoeff.insert(std::make_pair(EBDetId(eta,phi,EBDetId::ETAPHIMODE),coeff ));
  //    }
  //  }

  //  fclose(MisCalib);
  //   FILE* MisCalibEndCap;
  //   //char* calibfile="miscalibfile";
  //   MisCalibEndCap = fopen(miscalibfileEndCap_.c_str(),"r");

  //   int fileStatus2=0;
  //   int X=-1;
  //   int Y=-1;
  //   float coeff2=-1;
  std::map<EEDetId, float> OldCoeffEndCap;

  //  while(fileStatus2 != EOF) {
  //    fileStatus2 = fscanf(MisCalibEndCap,"%d %d %f\n",  &X,&Y,&coeff2);
  //    if(X!=-1&&Y!=-1&& coeff2!=-1){
  //      //     std::cout<<" We have read correctly the coefficient " << coeff << " corresponding to eta "<<eta<<" and  phi "<<phi<<std::endl;
  //      if(TestEEvalidDetId(X,Y,1)){
  //        OldCoeffEndCap.insert(std::make_pair(EEDetId(X,Y,1,EEDetId::XYMODE),coeff2 ));
  //      }
  //    }
  //  }

  // fclose(MisCalibEndCap);
  std::map<DetId, float>::const_iterator itmap;
  for (itmap = Univsolution.begin(); itmap != Univsolution.end(); itmap++) {
    const DetId Id(itmap->first);
    if (Id.subdetId() == 1) {
      const EBDetId IChannelDetId(itmap->first);
      if (IChannelDetId.ieta() < mineta_) {
        continue;
      }
      if (IChannelDetId.ieta() > maxeta_) {
        continue;
      }
      if (IChannelDetId.iphi() < minphi_) {
        continue;
      }
      if (IChannelDetId.iphi() > maxphi_) {
        continue;
      }
      //      float Compare=1;
      //      std::map<EBDetId,float>::iterator iter = OldCoeff.find(itmap->first);
      //      if( iter != OldCoeff.end() )Compare = iter->second;
      Map3Dcalib->Fill(IChannelDetId.ieta(), IChannelDetId.iphi(), itmap->second);
      calibs->Fill(itmap->second);
      //DiffCalibs->Fill(newCalibs[icry]-miscalib[IChannelDetId.ieta()-1][IChannelDetId.iphi()-21]);
      //RatioCalibs->Fill(newCalibs[icry]/miscalib[IChannelDetId.ieta()-1][IChannelDetId.iphi()-21]);
      if (IChannelDetId.ieta() < mineta_ + 2) {
        continue;
      }
      if (IChannelDetId.ieta() > maxeta_ - 2) {
        continue;
      }
      if (IChannelDetId.iphi() < minphi_ + 2) {
        continue;
      }
      if (IChannelDetId.iphi() > maxphi_ - 2) {
        continue;
      }
      write_calibrations.writeLine(IChannelDetId, itmap->second);
      calibinter->Fill(itmap->second);
      //       MapCorCalib->Fill(itmap->second,Compare);
      //       DiffCalibs->Fill(itmap->second-Compare);
      //       RatioCalibs->Fill(itmap->second*Compare);
    } else {
      const EEDetId IChannelDetId(itmap->first);
      //       if (IChannelDetId.ix()<0 ){continue;}
      //       if (IChannelDetId.ix()>100 ){continue;}
      //       if (IChannelDetId.iy()<0 ){continue;}
      //       if (IChannelDetId.iy()>100 ){continue;}
      //     std::map<EEDetId,float>::iterator iter = OldCoeffEndCap.find(itmap->first);
      //      float Compare=1;
      //      if( iter != OldCoeffEndCap.end() )Compare = iter->second;
      if (IChannelDetId.zside() < 0) {
        Map3DcalibEndCapMinus->Fill(IChannelDetId.ix(), IChannelDetId.iy(), itmap->second);
        calibsEndCapMinus->Fill(itmap->second);
        calibinterEndCapMinus->Fill(itmap->second);
        //  DiffCalibsEndCapMinus->Fill(itmap->second-Compare);
        //  RatioCalibsEndCapMinus->Fill(itmap->second*Compare);
        //  MapCorCalibEndCapMinus->Fill(itmap->second,Compare);
      } else {
        Map3DcalibEndCapPlus->Fill(IChannelDetId.ix(), IChannelDetId.iy(), itmap->second);
        calibsEndCapPlus->Fill(itmap->second);
        calibinterEndCapPlus->Fill(itmap->second);
        //  DiffCalibsEndCapPlus->Fill(itmap->second-Compare);
        //  RatioCalibsEndCapPlus->Fill(itmap->second*Compare);
        //  MapCorCalibEndCapPlus->Fill(itmap->second,Compare);
      }
      write_calibrations.writeLine(IChannelDetId, itmap->second);
    }
  }
  EventsAfterCuts->Write();

  // Book histograms
  e25->Write();
  e9->Write();
  scE->Write();
  trP->Write();
  EoP->Write();
  EoP_all->Write();
  calibs->Write();
  calibsEndCapMinus->Write();
  calibsEndCapPlus->Write();
  e9Overe25->Write();
  e25OverScE->Write();
  Map->Write();
  E25oP->Write();

  PinOverPout->Write();
  eSeedOverPout->Write();
  //   MisCalibs->Write();
  //   RatioCalibs->Write();
  //   DiffCalibs->Write();
  //   RatioCalibsNoCuts->Write();
  //   DiffCalibsNoCuts->Write();
  //   MisCalibsEndCapMinus->Write();
  //   MisCalibsEndCapPlus->Write();
  //   RatioCalibsEndCapMinus->Write();
  //   RatioCalibsEndCapPlus->Write();
  //   DiffCalibsEndCapMinus->Write();
  //   DiffCalibsEndCapPlus->Write();

  e25NoCuts->Write();
  e9NoCuts->Write();
  scENoCuts->Write();
  trPNoCuts->Write();
  EoPNoCuts->Write();
  calibsNoCuts->Write();
  e9Overe25NoCuts->Write();
  e25OverScENoCuts->Write();
  MapEndCapMinus->Write();
  MapEndCapPlus->Write();
  E25oPNoCuts->Write();
  Map3Dcalib->Write();
  Map3DcalibEndCapMinus->Write();
  Map3DcalibEndCapPlus->Write();
  Map3DcalibNoCuts->Write();
  calibinter->Write();
  calibinterEndCapMinus->Write();
  calibinterEndCapPlus->Write();
  calibinterNoCuts->Write();
  PinOverPoutNoCuts->Write();
  eSeedOverPoutNoCuts->Write();

  GeneralMap->Write();
  GeneralMapEndCapMinus->Write();
  GeneralMapEndCapPlus->Write();
  GeneralMapBeforePt->Write();
  GeneralMapEndCapMinusBeforePt->Write();
  GeneralMapEndCapPlusBeforePt->Write();

  MapCor1->Write();
  MapCor2->Write();
  MapCor3->Write();
  MapCor4->Write();
  MapCor5->Write();
  MapCor6->Write();
  MapCor7->Write();
  MapCor8->Write();
  MapCor9->Write();
  MapCor10->Write();
  MapCor11->Write();
  //  MapCorCalib->Write();

  MapCor1NoCuts->Write();
  MapCor2NoCuts->Write();
  MapCor3NoCuts->Write();
  MapCor4NoCuts->Write();
  MapCor5NoCuts->Write();
  MapCor6NoCuts->Write();
  MapCor7NoCuts->Write();
  MapCor8NoCuts->Write();
  MapCor9NoCuts->Write();
  MapCor10NoCuts->Write();
  MapCor11NoCuts->Write();
  //   MapCorCalibEndCapMinus->Write();
  //   MapCorCalibEndCapPlus->Write();

  MapCor1ESeed->Write();
  MapCor2ESeed->Write();
  MapCor3ESeed->Write();
  MapCor4ESeed->Write();
  MapCor5ESeed->Write();
  MapCor6ESeed->Write();
  MapCor7ESeed->Write();
  MapCor8ESeed->Write();
  MapCor9ESeed->Write();
  MapCor10ESeed->Write();
  MapCor11ESeed->Write();

  E25oPvsEta->Write();
  E25oPvsEtaEndCapMinus->Write();
  E25oPvsEtaEndCapPlus->Write();

  PinMinPout->Write();
  PinMinPoutNoCuts->Write();

  Error1->Write();
  Error2->Write();
  Error3->Write();
  Error1NoCuts->Write();
  Error2NoCuts->Write();
  Error3NoCuts->Write();

  eSeedOverPout2->Write();
  eSeedOverPout2NoCuts->Write();
  eSeedOverPout2ESeed->Write();

  hadOverEm->Write();
  hadOverEmNoCuts->Write();
  hadOverEmESeed->Write();

  f->Write();

  f->Close();
  //   if(MyL3Algo1)delete MyL3Algo1;
  //   if(UnivL3)delete UnivL3;
  //   if(MyHH)delete MyHH;
  //  delete f;

  std::cout << " " << std::endl;
  std::cout << "************* STATISTICS **************" << std::endl;
  std::cout << " Events Studied " << read_events << std::endl;
  std::cout << "Timing info:" << std::endl;
  std::cout << "CPU time usage  -- calibrating: " << cpu_time_used << " sec." << std::endl;

}

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

override

Definition at line 243 of file ElectronCalibrationUniv.cc.

{}

DetId ElectronCalibrationUniv::findMaxHit ( const std::vector< DetId > &  v1, const EBRecHitCollection *  EBhits, const EERecHitCollection *  EEhits  )

private

Definition at line 533 of file ElectronCalibrationUniv.cc.

References gather_cfg::cout, edm::SortedCollection< T, SORT >::end(), and edm::SortedCollection< T, SORT >::find().

Referenced by analyze().

                                                                            {
  //=================================================================================

  double currEnergy = 0.;
  DetId maxHit;

  for (std::vector<DetId>::const_iterator idsIt = v1.begin(); idsIt != v1.end(); ++idsIt) {
    if (idsIt->subdetId() == 1) {
      EBRecHitCollection::const_iterator itrechit;
      itrechit = EBhits->find(*idsIt);
      if (itrechit == EBhits->end()) {
        std::cout << "ElectronCalibration::findMaxHit: rechit not found! " << (EBDetId)(*idsIt) << std::endl;
        continue;
      }
      if (itrechit->energy() > currEnergy) {
        currEnergy = itrechit->energy();
        maxHit = *idsIt;
      }
    } else {
      EERecHitCollection::const_iterator itrechit;
      itrechit = EEhits->find(*idsIt);
      if (itrechit == EEhits->end()) {
        std::cout << "ElectronCalibration::findMaxHit: rechit not found! idsIt = " << (EEDetId)(*idsIt) << std::endl;
        continue;
      }

      if (itrechit->energy() > currEnergy) {
        currEnergy = itrechit->energy();
        maxHit = *idsIt;
      }
    }
  }

  return maxHit;
}

bool ElectronCalibrationUniv::TestEEvalidDetId ( int  crystal_ix, int  crystal_iy, int  iz  )

private

Definition at line 571 of file ElectronCalibrationUniv.cc.

References funct::abs().

                                                                                     {
  bool valid = false;
  if (crystal_ix < 1 || crystal_ix > 100 || crystal_iy < 1 || crystal_iy > 100 || abs(iz) != 1) {
    return valid;
  }
  if ((crystal_ix >= 1 && crystal_ix <= 3 && (crystal_iy <= 40 || crystal_iy > 60)) ||
      (crystal_ix >= 4 && crystal_ix <= 5 && (crystal_iy <= 35 || crystal_iy > 65)) ||
      (crystal_ix >= 6 && crystal_ix <= 8 && (crystal_iy <= 25 || crystal_iy > 75)) ||
      (crystal_ix >= 9 && crystal_ix <= 13 && (crystal_iy <= 20 || crystal_iy > 80)) ||
      (crystal_ix >= 14 && crystal_ix <= 15 && (crystal_iy <= 15 || crystal_iy > 85)) ||
      (crystal_ix >= 16 && crystal_ix <= 20 && (crystal_iy <= 13 || crystal_iy > 87)) ||
      (crystal_ix >= 21 && crystal_ix <= 25 && (crystal_iy <= 8 || crystal_iy > 92)) ||
      (crystal_ix >= 26 && crystal_ix <= 35 && (crystal_iy <= 5 || crystal_iy > 95)) ||
      (crystal_ix >= 36 && crystal_ix <= 39 && (crystal_iy <= 3 || crystal_iy > 97)) ||
      (crystal_ix >= 98 && crystal_ix <= 100 && (crystal_iy <= 40 || crystal_iy > 60)) ||
      (crystal_ix >= 96 && crystal_ix <= 97 && (crystal_iy <= 35 || crystal_iy > 65)) ||
      (crystal_ix >= 93 && crystal_ix <= 95 && (crystal_iy <= 25 || crystal_iy > 75)) ||
      (crystal_ix >= 88 && crystal_ix <= 92 && (crystal_iy <= 20 || crystal_iy > 80)) ||
      (crystal_ix >= 86 && crystal_ix <= 87 && (crystal_iy <= 15 || crystal_iy > 85)) ||
      (crystal_ix >= 81 && crystal_ix <= 85 && (crystal_iy <= 13 || crystal_iy > 87)) ||
      (crystal_ix >= 76 && crystal_ix <= 80 && (crystal_iy <= 8 || crystal_iy > 92)) ||
      (crystal_ix >= 66 && crystal_ix <= 75 && (crystal_iy <= 5 || crystal_iy > 95)) ||
      (crystal_ix >= 62 && crystal_ix <= 65 && (crystal_iy <= 3 || crystal_iy > 97)) ||
      ((crystal_ix == 40 || crystal_ix == 61) &&
       ((crystal_iy >= 46 && crystal_iy <= 55) || crystal_iy <= 3 || crystal_iy > 97)) ||
      ((crystal_ix == 41 || crystal_ix == 60) && crystal_iy >= 44 && crystal_iy <= 57) ||
      ((crystal_ix == 42 || crystal_ix == 59) && crystal_iy >= 43 && crystal_iy <= 58) ||
      ((crystal_ix == 43 || crystal_ix == 58) && crystal_iy >= 42 && crystal_iy <= 59) ||
      ((crystal_ix == 44 || crystal_ix == 45 || crystal_ix == 57 || crystal_ix == 56) && crystal_iy >= 41 &&
       crystal_iy <= 60) ||
      (crystal_ix >= 46 && crystal_ix <= 55 && crystal_iy >= 40 && crystal_iy <= 61)) {
    return valid;
  }
  valid = true;
  return valid;
}

Member Data Documentation

const std::string ElectronCalibrationUniv::calibAlgo_

private

Definition at line 79 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

int ElectronCalibrationUniv::calibClusterSize

private

Definition at line 104 of file ElectronCalibrationUniv.h.

Referenced by beginJob().

TH1F* ElectronCalibrationUniv::calibinter

private

Definition at line 181 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::calibinterEndCapMinus

private

Definition at line 182 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::calibinterEndCapPlus

private

Definition at line 183 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::calibinterNoCuts

private

Definition at line 184 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::calibs

private

Definition at line 144 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::calibsEndCapMinus

private

Definition at line 145 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::calibsEndCapPlus

private

Definition at line 146 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::calibsNoCuts

private

Definition at line 171 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

const int ElectronCalibrationUniv::clusterSize_

private

Definition at line 99 of file ElectronCalibrationUniv.h.

Referenced by analyze(), and beginJob().

const double ElectronCalibrationUniv::cut1_

private

Definition at line 88 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const double ElectronCalibrationUniv::cut2_

private

Definition at line 89 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const double ElectronCalibrationUniv::cut3_

private

Definition at line 90 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const double ElectronCalibrationUniv::cutCalo1_

private

Definition at line 96 of file ElectronCalibrationUniv.h.

const double ElectronCalibrationUniv::cutCalo2_

private

Definition at line 97 of file ElectronCalibrationUniv.h.

const double ElectronCalibrationUniv::cutEPCalo1_

private

Definition at line 92 of file ElectronCalibrationUniv.h.

const double ElectronCalibrationUniv::cutEPCalo2_

private

Definition at line 93 of file ElectronCalibrationUniv.h.

const double ElectronCalibrationUniv::cutEPin1_

private

Definition at line 94 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const double ElectronCalibrationUniv::cutEPin2_

private

Definition at line 95 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const double ElectronCalibrationUniv::cutESeed_

private

Definition at line 98 of file ElectronCalibrationUniv.h.

TH1F* ElectronCalibrationUniv::e25

private

Definition at line 138 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::e25NoCuts

private

Definition at line 166 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::E25oP

private

Definition at line 150 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::E25oPNoCuts

private

Definition at line 176 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::E25oPvsEta

private

Definition at line 235 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::E25oPvsEtaEndCapMinus

private

Definition at line 236 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::E25oPvsEtaEndCapPlus

private

Definition at line 237 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::e25OverScE

private

Definition at line 148 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::e25OverScENoCuts

private

Definition at line 173 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::e9

private

Definition at line 139 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::e9NoCuts

private

Definition at line 167 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::e9Overe25

private

Definition at line 147 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::e9Overe25NoCuts

private

Definition at line 172 of file ElectronCalibrationUniv.h.

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

const edm::InputTag ElectronCalibrationUniv::ebRecHitLabel_

private

Definition at line 75 of file ElectronCalibrationUniv.h.

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

private

Definition at line 119 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const edm::InputTag ElectronCalibrationUniv::eeRecHitLabel_

private

Definition at line 76 of file ElectronCalibrationUniv.h.

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

private

Definition at line 120 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const int ElectronCalibrationUniv::elecclass_

private

Definition at line 100 of file ElectronCalibrationUniv.h.

Referenced by analyze().

const edm::InputTag ElectronCalibrationUniv::electronLabel_

private

Definition at line 77 of file ElectronCalibrationUniv.h.

const double ElectronCalibrationUniv::elePt_

private

Definition at line 83 of file ElectronCalibrationUniv.h.

Referenced by analyze().

std::vector<float> ElectronCalibrationUniv::EnergyVector

private

Definition at line 106 of file ElectronCalibrationUniv.h.

Referenced by analyze(), and endJob().

std::vector<float> ElectronCalibrationUniv::EnergyVectorNoCuts

private

Definition at line 111 of file ElectronCalibrationUniv.h.

TH1F* ElectronCalibrationUniv::EoP

private

Definition at line 142 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::EoP_all

private

Definition at line 143 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::EoPNoCuts

private

Definition at line 170 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::Error1

private

Definition at line 242 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::Error1NoCuts

private

Definition at line 245 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::Error2

private

Definition at line 243 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::Error2NoCuts

private

Definition at line 246 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::Error3

private

Definition at line 244 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::Error3NoCuts

private

Definition at line 247 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::eSeedOverPout

private

Definition at line 153 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::eSeedOverPout2

private

Definition at line 249 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::eSeedOverPout2ESeed

private

Definition at line 251 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::eSeedOverPout2NoCuts

private

Definition at line 250 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::eSeedOverPoutNoCuts

private

Definition at line 187 of file ElectronCalibrationUniv.h.

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

int ElectronCalibrationUniv::etaMax

private

Definition at line 105 of file ElectronCalibrationUniv.h.

Referenced by beginJob().

int ElectronCalibrationUniv::etaMin

private

Definition at line 105 of file ElectronCalibrationUniv.h.

Referenced by beginJob().

std::vector<std::vector<float> > ElectronCalibrationUniv::EventMatrix

private

Definition at line 108 of file ElectronCalibrationUniv.h.

Referenced by analyze(), and endJob().

std::vector<std::vector<float> > ElectronCalibrationUniv::EventMatrixNoCuts

private

Definition at line 113 of file ElectronCalibrationUniv.h.

TH1F* ElectronCalibrationUniv::EventsAfterCuts

private

Definition at line 136 of file ElectronCalibrationUniv.h.

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

TFile* ElectronCalibrationUniv::f

private

Definition at line 134 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), endJob(), and ztail.Decoder::initial_synchronize().

TH2F* ElectronCalibrationUniv::GeneralMap

private

Definition at line 189 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::GeneralMapBeforePt

private

Definition at line 192 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::GeneralMapEndCapMinus

private

Definition at line 190 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::GeneralMapEndCapMinusBeforePt

private

Definition at line 193 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::GeneralMapEndCapPlus

private

Definition at line 191 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::GeneralMapEndCapPlusBeforePt

private

Definition at line 194 of file ElectronCalibrationUniv.h.

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

const edm::EDGetTokenT<reco::GsfElectronCollection> ElectronCalibrationUniv::gsfElectronToken_

private

Definition at line 121 of file ElectronCalibrationUniv.h.

Referenced by analyze().

TH1F* ElectronCalibrationUniv::hadOverEm

private

Definition at line 253 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::hadOverEmESeed

private

Definition at line 255 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH1F* ElectronCalibrationUniv::hadOverEmNoCuts

private

Definition at line 254 of file ElectronCalibrationUniv.h.

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

int ElectronCalibrationUniv::keventweight

private

Definition at line 104 of file ElectronCalibrationUniv.h.

const int ElectronCalibrationUniv::keventweight_

private

Definition at line 82 of file ElectronCalibrationUniv.h.

Referenced by beginJob().

TH2F* ElectronCalibrationUniv::Map

private

Definition at line 149 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::Map3Dcalib

private

Definition at line 177 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::Map3DcalibEndCapMinus

private

Definition at line 178 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::Map3DcalibEndCapPlus

private

Definition at line 179 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::Map3DcalibNoCuts

private

Definition at line 180 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor1

private

Definition at line 196 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor10

private

Definition at line 205 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor10ESeed

private

Definition at line 232 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor10NoCuts

private

Definition at line 218 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor11

private

Definition at line 206 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor11ESeed

private

Definition at line 233 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor11NoCuts

private

Definition at line 219 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor1ESeed

private

Definition at line 223 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor1NoCuts

private

Definition at line 209 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor2

private

Definition at line 197 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor2ESeed

private

Definition at line 224 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor2NoCuts

private

Definition at line 210 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor3

private

Definition at line 198 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor3ESeed

private

Definition at line 225 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor3NoCuts

private

Definition at line 211 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor4

private

Definition at line 199 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor4ESeed

private

Definition at line 226 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor4NoCuts

private

Definition at line 212 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor5

private

Definition at line 200 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor5ESeed

private

Definition at line 227 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor5NoCuts

private

Definition at line 213 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor6

private

Definition at line 201 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor6ESeed

private

Definition at line 228 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor6NoCuts

private

Definition at line 214 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor7

private

Definition at line 202 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor7ESeed

private

Definition at line 229 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor7NoCuts

private

Definition at line 215 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor8

private

Definition at line 203 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor8ESeed

private

Definition at line 230 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor8NoCuts

private

Definition at line 216 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor9

private

Definition at line 204 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapCor9ESeed

private

Definition at line 231 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

TH2F* ElectronCalibrationUniv::MapCor9NoCuts

private

Definition at line 217 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapEndCapMinus

private

Definition at line 174 of file ElectronCalibrationUniv.h.

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

TH2F* ElectronCalibrationUniv::MapEndCapPlus

private

Definition at line 175 of file ElectronCalibrationUniv.h.

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

std::vector<int> ElectronCalibrationUniv::MaxCCeta

private

Definition at line 109 of file ElectronCalibrationUniv.h.

Referenced by endJob().

std::vector<int> ElectronCalibrationUniv::MaxCCetaNoCuts

private

Definition at line 114 of file ElectronCalibrationUniv.h.

std::vector<int> ElectronCalibrationUniv::MaxCCphi

private

Definition at line 110 of file ElectronCalibrationUniv.h.

Referenced by endJob().

std::vector<int> ElectronCalibrationUniv::MaxCCphiNoCuts

private

Definition at line 115 of file ElectronCalibrationUniv.h.

const int ElectronCalibrationUniv::maxeta_

private

Definition at line 84 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

const int ElectronCalibrationUniv::maxphi_

private

Definition at line 86 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

const int ElectronCalibrationUniv::mineta_

private

Definition at line 85 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

const int ElectronCalibrationUniv::minphi_

private

Definition at line 87 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

const std::string ElectronCalibrationUniv::miscalibfile_

private

Definition at line 80 of file ElectronCalibrationUniv.h.

const std::string ElectronCalibrationUniv::miscalibfileEndCap_

private

Definition at line 81 of file ElectronCalibrationUniv.h.

HouseholderDecomposition* ElectronCalibrationUniv::MyHH

private

Definition at line 185 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

MinL3Algorithm* ElectronCalibrationUniv::MyL3Algo1

private

Definition at line 116 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

const int ElectronCalibrationUniv::numevent_

private

Definition at line 91 of file ElectronCalibrationUniv.h.

int ElectronCalibrationUniv::phiMax

private

Definition at line 105 of file ElectronCalibrationUniv.h.

Referenced by beginJob().

int ElectronCalibrationUniv::phiMin

private

Definition at line 105 of file ElectronCalibrationUniv.h.

Referenced by beginJob().

TH1F* ElectronCalibrationUniv::PinMinPout

private

Definition at line 239 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::PinMinPoutNoCuts

private

Definition at line 240 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::PinOverPout

private

Definition at line 152 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::PinOverPoutNoCuts

private

Definition at line 186 of file ElectronCalibrationUniv.h.

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

int ElectronCalibrationUniv::read_events

private

Definition at line 102 of file ElectronCalibrationUniv.h.

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

const std::string ElectronCalibrationUniv::rootfile_

private

Definition at line 78 of file ElectronCalibrationUniv.h.

Referenced by beginJob().

TH1F* ElectronCalibrationUniv::scE

private

Definition at line 140 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::scENoCuts

private

Definition at line 168 of file ElectronCalibrationUniv.h.

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

std::vector<float> ElectronCalibrationUniv::solution

private

Definition at line 125 of file ElectronCalibrationUniv.h.

Referenced by endJob().

std::vector<float> ElectronCalibrationUniv::solutionNoCuts

private

Definition at line 126 of file ElectronCalibrationUniv.h.

const CaloTopology* ElectronCalibrationUniv::theCaloTopology_

private

Definition at line 123 of file ElectronCalibrationUniv.h.

Referenced by analyze(), and beginRun().

const edm::ESGetToken<CaloTopology, CaloTopologyRecord> ElectronCalibrationUniv::topologyToken_

private

Definition at line 122 of file ElectronCalibrationUniv.h.

Referenced by beginRun().

TH1F* ElectronCalibrationUniv::trP

private

Definition at line 141 of file ElectronCalibrationUniv.h.

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

TH1F* ElectronCalibrationUniv::trPNoCuts

private

Definition at line 169 of file ElectronCalibrationUniv.h.

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

std::vector<std::vector<DetId> > ElectronCalibrationUniv::UnivEventIds

private

Definition at line 129 of file ElectronCalibrationUniv.h.

Referenced by analyze(), and endJob().

MinL3AlgoUniv<DetId>* ElectronCalibrationUniv::UnivL3

private

Definition at line 117 of file ElectronCalibrationUniv.h.

Referenced by beginJob(), and endJob().

std::map<DetId, float> ElectronCalibrationUniv::Univsolution

private

Definition at line 130 of file ElectronCalibrationUniv.h.

Referenced by endJob().

std::vector<float> ElectronCalibrationUniv::WeightVector

private

Definition at line 107 of file ElectronCalibrationUniv.h.

std::vector<float> ElectronCalibrationUniv::WeightVectorNoCuts

private

Definition at line 112 of file ElectronCalibrationUniv.h.