hcalCalib Class Reference

#include <hcalCalib.h>

Inheritance diagram for hcalCalib:

Public Member Functions
void	Begin (TTree *tree) override
void	GetCoefFromMtrxInvOfAve ()
Int_t	GetEntry (Long64_t entry, Int_t getall=0) override
TList *	GetOutputList () const override
	hcalCalib (TTree *=nullptr)
void	Init (TTree *tree) override
void	makeTextFile ()
Bool_t	Notify () override
Bool_t	Process (Long64_t entry) override
Bool_t	ReadPhiSymCor ()
void	SetApplyPhiSymCorFlag (Bool_t b)
void	SetCalibAbsIEtaMax (Int_t i)
void	SetCalibAbsIEtaMin (Int_t i)
void	SetCalibMethod (const TString &s)
void	SetCalibType (const TString &s)
void	SetCaloGeometry (const CaloGeometry *g, const HcalTopology *topo)
void	SetCombinePhiFlag (Bool_t b)
void	SetConeMaxDist (Float_t d)
void	SetHbClusterSize (Int_t i)
void	SetHeClusterSize (Int_t i)
void	SetHistoFileName (const TString &filename)
void	SetInputList (TList *input) override
void	SetMaxEOverP (Float_t e)
void	SetMaxEtThirdJet (Float_t et)
void	SetMaxProbeJetEmFrac (Float_t f)
void	SetMaxTagJetAbsEta (Float_t e)
void	SetMaxTagJetEmFrac (Float_t f)
void	SetMaxTargetE (Float_t e)
void	SetMaxTrkEmE (Float_t e)
void	SetMinCellE (Float_t e)
void	SetMinDPhiDiJets (Float_t dphi)
void	SetMinEOverP (Float_t e)
void	SetMinProbeJetAbsEta (Float_t e)
void	SetMinTagJetEt (Float_t e)
void	SetMinTargetE (Float_t e)
void	SetObject (TObject *obj) override
void	SetOption (const char *option) override
void	SetOutputCorCoefFileName (const TString &filename)
void	SetPhiSymCorFileName (const TString &filename)
void	SetSumDepthsFlag (Bool_t b)
void	SetSumSmallDepthsFlag (Bool_t b)
void	SetUseConeClustering (Bool_t b)
void	Terminate () override
Int_t	Version () const override
	~hcalCalib () override

Public Attributes
Bool_t	APPLY_PHI_SYM_COR_FLAG
TBranch *	b_cells
TBranch *	b_emEnergy
TBranch *	b_etVetoJet
TBranch *	b_eventNumber
TBranch *	b_iEtaHit
TBranch *	b_iPhiHit
TBranch *	b_probeJetEmFrac
TBranch *	b_probeJetP4
TBranch *	b_runNumber
TBranch *	b_tagJetEmFrac
TBranch *	b_tagJetP4
TBranch *	b_targetE
TBranch *	b_xTrkEcal
TBranch *	b_xTrkHcal
TBranch *	b_yTrkEcal
TBranch *	b_yTrkHcal
TBranch *	b_zTrkEcal
TBranch *	b_zTrkHcal
Int_t	CALIB_ABS_IETA_MAX
Int_t	CALIB_ABS_IETA_MIN
TString	CALIB_METHOD
TString	CALIB_TYPE
std::vector< std::vector< Float_t > >	cellEnergies
std::vector< std::vector< UInt_t > >	cellIds
TClonesArray *	cells
Bool_t	COMBINE_PHI
Float_t	emEnergy
Float_t	etVetoJet
UInt_t	eventNumber
	pointer to the analyzed TTree or TChain More...
TTree *	fChain
TH1F *	h1_allTrkP
TH1F *	h1_corResp
TH1F *	h1_corRespBarrel
TH1F *	h1_corRespEndcap
TH1F *	h1_corRespIEta [48]
TH1F *	h1_numEventsTwrIEta
TH1F *	h1_rawResp
TH1F *	h1_rawRespBarrel
TH1F *	h1_rawRespEndcap
TH1F *	h1_rawSumE
TH1F *	h1_selTrkP_iEta10
TH1F *	h1_trkP
TH2F *	h2_dHitRefBarrel
TH2F *	h2_dHitRefEndcap
Int_t	HB_CLUSTER_SIZE
Int_t	HE_CLUSTER_SIZE
TString	HISTO_FILENAME
TFile *	histoFile
std::map< Int_t, Float_t >	iEtaCoefMap
Int_t	iEtaHit
UInt_t	iPhiHit
Float_t	MAX_CONE_DIST
Float_t	MAX_EOVERP
Float_t	MAX_ET_THIRD_JET
Float_t	MAX_PROBEJET_EMFRAC
Float_t	MAX_TAGJET_ABSETA
Float_t	MAX_TAGJET_EMFRAC
Float_t	MAX_TARGET_E
Float_t	MAX_TRK_EME
Float_t	MIN_CELL_E
Float_t	MIN_DPHI_DIJETS
Float_t	MIN_EOVERP
Float_t	MIN_PROBEJET_ABSETA
Float_t	MIN_TAGJET_ET
Float_t	MIN_TARGET_E
UInt_t	nEvents
TString	OUTPUT_COR_COEF_FILENAME
TString	PHI_SYM_COR_FILENAME
std::map< UInt_t, Float_t >	phiSymCor
Float_t	probeJetEmFrac
TLorentzVector *	probeJetP4
std::vector< std::pair< Int_t, UInt_t > >	refIEtaIPhi
UInt_t	runNumber
std::map< UInt_t, Float_t >	solution
Bool_t	SUM_DEPTHS
Bool_t	SUM_SMALL_DEPTHS
Float_t	tagJetEmFrac
TLorentzVector *	tagJetP4
Float_t	targetE
std::vector< Float_t >	targetEnergies
const CaloGeometry *	theCaloGeometry
const HcalTopology *	topo_
Bool_t	USE_CONE_CLUSTERING
Float_t	xTrkEcal
Float_t	xTrkHcal
Float_t	yTrkEcal
Float_t	yTrkHcal
Float_t	zTrkEcal
Float_t	zTrkHcal

Detailed Description

Definition at line 46 of file hcalCalib.h.

Constructor & Destructor Documentation

hcalCalib()

hcalCalib::hcalCalib ( TTree *  = nullptr )

inline

Definition at line 121 of file hcalCalib.h.

{}

~hcalCalib()

hcalCalib::~hcalCalib ( )

inlineoverride

Definition at line 122 of file hcalCalib.h.

{}

Member Function Documentation

Begin()

void hcalCalib::Begin ( TTree *  tree )

override

Definition at line 28 of file hcalCalib.cc.

References APPLY_PHI_SYM_COR_FLAG, CALIB_TYPE, beamvalidation::exit(), h1_allTrkP, h1_corResp, h1_corRespBarrel, h1_corRespEndcap, h1_corRespIEta, h1_numEventsTwrIEta, h1_rawResp, h1_rawRespBarrel, h1_rawRespEndcap, h1_rawSumE, h1_selTrkP_iEta10, h1_trkP, h2_dHitRefBarrel, h2_dHitRefEndcap, HISTO_FILENAME, histoFile, HltBtagPostValidation_cff::histoName, mps_fire::i, l1tTowerCalibrationProducer_cfi::iEta, nEvents, fileinputsource_cfi::option, and ReadPhiSymCor().

                              {
  TString option = GetOption();

  nEvents = 0;

  if (APPLY_PHI_SYM_COR_FLAG && !ReadPhiSymCor()) {
    edm::LogError("HcalCalib") << "\nERROR: Failed to read the phi symmetry corrections.\n"
                               << "Check if the filename is correct. If the corrections are not needed, set the "
                                  "corresponding flag to \"false\"\n"
                               << "\nThe program will be terminated\n";

    exit(1);
  }

  //  cellEnergies.reserve(1000000);
  //  cellIds.reserve(1000000);
  //  targetEnergies.reserve(1000000);

  histoFile = new TFile(HISTO_FILENAME.Data(), "RECREATE");

  h1_trkP = new TH1F("h1_trkP", "Track momenta; p_{trk} (GeV); Number of tracks", 100, 0, 200);
  h1_allTrkP = new TH1F("h1_allTrkP", "Track momenta - all tracks; p_{trk} (GeV); Number of tracks", 100, 0, 200);

  h1_selTrkP_iEta10 = new TH1F(
      "h1_selTrkP_iEta10", "Track momenta - tracks with |iEta|<10; p_{trk} (GeV); Number of tracks", 100, 0, 200);

  if (CALIB_TYPE == "ISO_TRACK")
    h1_rawSumE = new TH1F("h1_rawSumE", "Cluster Energy; E_{cl} (GeV); Number of tracks", 100, 0, 200);
  else
    h1_rawSumE = new TH1F("h1_rawSumE", "Cluster Energy; E_{cl} (GeV); Number of tracks", 1000, 0, 2000);

  h1_rawResp = new TH1F("h1_rawResp", "Uncorrected response: |iEta|<24;  E_{had}/p; Number of tracks", 300, 0, 3);
  h1_corResp = new TH1F("h1_corResp", "Corrected response: |iEta|<24; E_{had}/p; Number of tracks", 300, 0, 3);

  h1_rawRespBarrel =
      new TH1F("h1_rawRespBarrel", "Uncorrected response: |iEta|<15;  E_{had}/p; Number of tracks", 300, 0, 3);
  h1_corRespBarrel =
      new TH1F("h1_corRespBarrel", "Corrected response: |iEta|<15; E_{had}/p; Number of tracks", 300, 0, 3);

  h1_rawRespEndcap =
      new TH1F("h1_rawRespEndcap", "Uncorrected response:  17<|iEta|<24;  E_{had}/p; Number of tracks", 300, 0, 3);
  h1_corRespEndcap =
      new TH1F("h1_corRespEndcap", "Corrected response: 17<|iEta|<24; E_{had}/p; Number of tracks", 300, 0, 3);

  h1_numEventsTwrIEta = new TH1F("h1_numEventsTwrIEta", "h1_numEventsTwrIEta", 80, -40, 40);

  h2_dHitRefBarrel = new TH2F("h2_dHitRefBarrel",
                              "{#Delta}i{#phi} vs {#Delta}i{#eta} of hit and most energetic "
                              "tower(|i{#eta}|<16);{#Delta}i{#eta}; {#Delta}i{#phi}",
                              10,
                              -5,
                              5,
                              10,
                              -5,
                              5);
  h2_dHitRefEndcap = new TH2F("h2_dHitRefEndcap",
                              "{#Delta}i{#phi} vs {#Delta}i{#eta} of hit and most energetic tower (16<|i{#eta}|<25) "
                              ";{#Delta}i{#eta}; {#Delta}i{#phi}",
                              10,
                              -5,
                              5,
                              10,
                              -5,
                              5);

  TString histoName = "isoTrack_";

  for (Int_t i = 0; i < 48; ++i) {
    Long_t iEta;
    if (i < 24)
      iEta = i - 24;
    else
      iEta = i - 23;
    TString hn = histoName + iEta;
    h1_corRespIEta[i] = new TH1F(hn, hn, 300, 0, 3.0);
  }

}  // end of Begin()

GetCoefFromMtrxInvOfAve()

void hcalCalib::GetCoefFromMtrxInvOfAve

(

)

Definition at line 471 of file hcalCalib.cc.

References A, funct::abs(), b, CALIB_ABS_IETA_MAX, CALIB_ABS_IETA_MIN, CALIB_METHOD, cellEnergies, cellIds, mps_fire::end, spr::find(), mps_fire::i, l1tTowerCalibrationProducer_cfi::iEta, HcalDetId::ieta(), iEtaCoefMap, dqmiolumiharvest::j, refIEtaIPhi, findQualityFiles::size, solution, targetEnergies, and x.

Referenced by Terminate().

                                        {
  // these maps are keyed by iEta
  std::map<Int_t, Float_t> aveTargetE;
  std::map<Int_t, Int_t> nEntries;  // count hits

  //  iEtaRef  iEtaCell, energy
  std::map<Int_t, std::map<Int_t, Float_t> > aveHitE;  // add energies in the loop, normalize after that

  for (unsigned int i = 0; i < cellEnergies.size(); ++i) {
    Int_t iEtaRef = refIEtaIPhi[i].first;
    aveTargetE[iEtaRef] += targetEnergies[i];
    nEntries[iEtaRef]++;

    // for hybrid method: matrix inv of averages preceeded by L3
    if (CALIB_METHOD == "L3_AND_MTRX_INV") {
      for (unsigned int j = 0; j < (cellEnergies[i]).size(); ++j) {
        aveHitE[iEtaRef][HcalDetId(cellIds[i][j]).ieta()] += (solution[cellIds[i][j]] * cellEnergies[i][j]);
      }
    } else if (CALIB_METHOD == "MATRIX_INV_OF_ETA_AVE") {
      for (unsigned int j = 0; j < (cellEnergies[i]).size(); ++j) {
        aveHitE[iEtaRef][HcalDetId(cellIds[i][j]).ieta()] += cellEnergies[i][j];
      }
    }

  }  // end of loop of entries

  // scale by number of entries to get the averages
  Float_t norm = 1.0;
  for (std::map<Int_t, Float_t>::iterator m_it = aveTargetE.begin(); m_it != aveTargetE.end(); ++m_it) {
    Int_t iEta = m_it->first;
    norm = (nEntries[iEta] > 0) ? 1.0 / (nEntries[iEta]) : 1.0;
    aveTargetE[iEta] *= norm;

    std::map<Int_t, Float_t>::iterator n_it = (aveHitE[iEta]).begin();

    for (; n_it != (aveHitE[iEta]).end(); ++n_it) {
      (n_it->second) *= norm;
    }

  }  // end of scaling by number of entries

  Int_t ONE_SIDE_IETA_RANGE = CALIB_ABS_IETA_MAX - CALIB_ABS_IETA_MIN + 1;

  // conversion from iEta to index for the linear system
  // contains elements only in the valid range for *matrix inversion*
  std::vector<Int_t> iEtaList;

  for (Int_t i = -CALIB_ABS_IETA_MAX; i <= CALIB_ABS_IETA_MAX; ++i) {
    if (abs(i) < CALIB_ABS_IETA_MIN)
      continue;
    iEtaList.push_back(i);
  }

  TMatrixD A(2 * ONE_SIDE_IETA_RANGE, 2 * ONE_SIDE_IETA_RANGE);
  TMatrixD b(2 * ONE_SIDE_IETA_RANGE, 1);
  TMatrixD x(2 * ONE_SIDE_IETA_RANGE, 1);

  for (Int_t i = 0; i < 2 * ONE_SIDE_IETA_RANGE; ++i) {
    for (Int_t j = 0; j < 2 * ONE_SIDE_IETA_RANGE; ++j) {
      A(i, j) = 0.0;
    }
  }

  for (UInt_t i = 0; i < iEtaList.size(); ++i) {
    b(i, 0) = aveTargetE[iEtaList[i]];

    std::map<Int_t, Float_t>::iterator n_it = aveHitE[iEtaList[i]].begin();
    for (; n_it != aveHitE[iEtaList[i]].end(); ++n_it) {
      if (fabs(n_it->first) > CALIB_ABS_IETA_MAX || fabs(n_it->first) < CALIB_ABS_IETA_MIN)
        continue;
      Int_t j = Int_t(find(iEtaList.begin(), iEtaList.end(), n_it->first) - iEtaList.begin());
      A(i, j) = n_it->second;
    }
  }

  TMatrixD coef = b;
  TDecompQRH qrh(A);
  Bool_t hasSolution = qrh.MultiSolve(coef);

  if (hasSolution && (CALIB_METHOD == "L3_AND_MTRX_INV" || CALIB_METHOD == "MATRIX_INV_OF_ETA_AVE")) {
    for (UInt_t i = 0; i < iEtaList.size(); ++i) {
      iEtaCoefMap[iEtaList[i]] = coef(i, 0);
    }
  }
}

GetEntry()

Int_t hcalCalib::GetEntry ( Long64_t  entry, Int_t  getall = 0  )

inlineoverride

Definition at line 129 of file hcalCalib.h.

References mps_splice::entry, and fChain.

Referenced by Process().

                                                            {
    return fChain ? fChain->GetTree()->GetEntry(entry, getall) : 0;
  }

GetOutputList()

TList* hcalCalib::GetOutputList ( ) const

inlineoverride

Definition at line 135 of file hcalCalib.h.

{ return fOutput; }

Init()

void hcalCalib::Init ( TTree *  tree )

inlineoverride

Definition at line 695 of file hcalCalib.cc.

References b_cells, b_emEnergy, b_etVetoJet, b_eventNumber, b_iEtaHit, b_iPhiHit, b_probeJetEmFrac, b_probeJetP4, b_runNumber, b_tagJetEmFrac, b_tagJetP4, b_targetE, b_xTrkEcal, b_xTrkHcal, b_yTrkEcal, b_yTrkHcal, b_zTrkEcal, b_zTrkHcal, cells, emEnergy, etVetoJet, eventNumber, fChain, iEtaHit, iPhiHit, probeJetEmFrac, probeJetP4, runNumber, tagJetEmFrac, tagJetP4, targetE, compare::tree, xTrkEcal, xTrkHcal, yTrkEcal, yTrkHcal, zTrkEcal, and zTrkHcal.

                                       {
  // The Init() function is called when the selector needs to initialize
  // a new tree or chain. Typically here the branch addresses and branch
  // pointers of the tree will be set.
  // It is normaly not necessary to make changes to the generated
  // code, but the routine can be extended by the user if needed.
  // Init() will be called many times when running on PROOF
  // (once per file to be processed).

  // Set object pointer
  cells = nullptr;
  tagJetP4 = nullptr;
  probeJetP4 = nullptr;

  // Set branch addresses and branch pointers
  if (!tree)
    return;
  fChain = tree;

  //      fChain->SetMakeClass(1);

  fChain->SetBranchAddress("eventNumber", &eventNumber, &b_eventNumber);
  fChain->SetBranchAddress("runNumber", &runNumber, &b_runNumber);
  fChain->SetBranchAddress("iEtaHit", &iEtaHit, &b_iEtaHit);
  fChain->SetBranchAddress("iPhiHit", &iPhiHit, &b_iPhiHit);
  fChain->SetBranchAddress("cells", &cells, &b_cells);
  fChain->SetBranchAddress("emEnergy", &emEnergy, &b_emEnergy);
  fChain->SetBranchAddress("targetE", &targetE, &b_targetE);
  fChain->SetBranchAddress("etVetoJet", &etVetoJet, &b_etVetoJet);

  fChain->SetBranchAddress("xTrkHcal", &xTrkHcal, &b_xTrkHcal);
  fChain->SetBranchAddress("yTrkHcal", &yTrkHcal, &b_yTrkHcal);
  fChain->SetBranchAddress("zTrkHcal", &zTrkHcal, &b_zTrkHcal);
  fChain->SetBranchAddress("xTrkEcal", &xTrkEcal, &b_xTrkEcal);
  fChain->SetBranchAddress("yTrkEcal", &yTrkEcal, &b_yTrkEcal);
  fChain->SetBranchAddress("zTrkEcal", &zTrkEcal, &b_zTrkEcal);

  fChain->SetBranchAddress("tagJetEmFrac", &tagJetEmFrac, &b_tagJetEmFrac);
  fChain->SetBranchAddress("probeJetEmFrac", &probeJetEmFrac, &b_probeJetEmFrac);

  fChain->SetBranchAddress("tagJetP4", &tagJetP4, &b_tagJetP4);
  fChain->SetBranchAddress("probeJetP4", &probeJetP4, &b_probeJetP4);
}

makeTextFile()

void hcalCalib::makeTextFile

(

)

Definition at line 620 of file hcalCalib.cc.

References funct::abs(), CALIB_ABS_IETA_MAX, CALIB_ABS_IETA_MIN, CALIB_METHOD, COMBINE_PHI, ztail::d, MillePedeFileConverter_cfg::e, PVValHelper::eta, HcalEndcap, HcalOuter, l1ctLayer2EG_cff::id, iEtaCoefMap, OUTPUT_COR_COEF_FILENAME, phi, nano_mu_digi_cff::rawId, heavyFlavorDQMFirstStep_cff::sdName, solution, SUM_DEPTHS, topo_, and HcalTopology::valid().

Referenced by Terminate().

                             {
  //{ HcalEmpty=0, HcalBarrel=1, HcalEndcap=2, HcalOuter=3, HcalForward=4, HcalTriggerTower=5, HcalOther=7 };

  TString sdName[8] = {"EMPTY", "HB", "HE", "HO", "HF", "TRITWR", "UNDEF", "OTHER"};

  FILE* constFile = fopen(OUTPUT_COR_COEF_FILENAME.Data(), "w+");

  // header of the constants file
  fprintf(constFile, "%1s%16s%16s%16s%16s%9s%11s\n", "#", "eta", "phi", "depth", "det", "value", "DetId");

  // Order loops to produce sequence of constants as for phi symmetry

  for (Int_t sd = 1; sd <= 4; sd++) {
    for (Int_t e = 1; e <= 50; e++) {
      Int_t eta;

      for (Int_t side = 0; side < 2; side++) {
        eta = (side == 0) ? -e : e;  //ta *= (-1);

        for (Int_t phi = 1; phi <= 72; phi++) {
          for (Int_t d = 1; d < 5; d++) {
            if (!topo_->valid(HcalDetId(HcalSubdetector(sd), eta, phi, d)))
              continue;
            HcalDetId id(HcalSubdetector(sd), eta, phi, d);
            Float_t corrFactor = 1.0;

            if (abs(eta) >= CALIB_ABS_IETA_MIN && abs(eta) <= CALIB_ABS_IETA_MAX && HcalSubdetector(sd) != HcalOuter) {
              //        if (abs(eta)>=CALIB_ABS_IETA_MIN && abs(eta)<=22 && HcalSubdetector(sd)!=HcalOuter) {

              // need some care when depths were summed for iEta=16 =>
              // the coeficients are saved in depth 1 of HB: affects
              Int_t subdetInd = sd;

              if (abs(eta) == 16 && HcalSubdetector(sd) == HcalEndcap && SUM_DEPTHS) {
                subdetInd = 1;
              }

              if (CALIB_METHOD == "L3" || CALIB_METHOD == "L3_AND_MTRX_INV") {
                if (SUM_DEPTHS && COMBINE_PHI)
                  corrFactor = solution[HcalDetId(HcalSubdetector(subdetInd), eta, 1, 1)];
                else if (SUM_DEPTHS)
                  corrFactor = solution[HcalDetId(HcalSubdetector(subdetInd), eta, phi, 1)];
                else if (COMBINE_PHI)
                  corrFactor = solution[HcalDetId(HcalSubdetector(sd), eta, 1, d)];
                else
                  corrFactor = solution[HcalDetId(HcalSubdetector(sd), eta, phi, d)];

                // Remark: a new case was added (sumSmallDepths) where the first two depths in towers 15,16
                // are summed and stored in  depth 1.
                // For now we create the correction coef for depth 2 (set equal to depth 1)
                // after the call to the L3 minimizer so that this case is also handled without modifying the
                // logic above. Probably it is better to move it here?

              }  // L3

              else if (CALIB_METHOD == "MATRIX_INV_OF_ETA_AVE") {
                corrFactor = iEtaCoefMap[eta];
              }

              else if (CALIB_METHOD == "ISO_TRK_PHI_SYM") {
                corrFactor = solution[HcalDetId(HcalSubdetector(sd), eta, phi, d)];
              }

            }  // if within the calibration range

            fprintf(constFile, "%17i%16i%16i%16s%9.5f%11X\n", eta, phi, d, sdName[sd].Data(), corrFactor, id.rawId());
          }
        }
      }
    }
  }

  return;
}

Notify()

Bool_t hcalCalib::Notify ( )

inlineoverride

Definition at line 739 of file hcalCalib.cc.

                                {
  // The Notify() function is called when a new file is opened. This
  // can be either for a new TTree in a TChain or when when a new TTree
  // is started when using PROOF. It is normaly not necessary to make changes
  // to the generated code, but the routine can be extended by the
  // user if needed. The return value is currently not used.

  return kTRUE;
}

Process()

Bool_t hcalCalib::Process ( Long64_t  entry )

override

Definition at line 111 of file hcalCalib.cc.

References funct::abs(), APPLY_PHI_SYM_COR_FLAG, CALIB_TYPE, cellEnergies, cellIds, cells, COMBINE_PHI, combinePhi(), TCell::e(), emEnergy, mps_splice::entry, etVetoJet, filterCells3x3(), filterCells5x5(), filterCellsInCone(), GetEntry(), getIEtaIPhiForHighestE(), h1_allTrkP, h1_rawSumE, h1_selTrkP_iEta10, h2_dHitRefBarrel, h2_dHitRefEndcap, HB_CLUSTER_SIZE, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HE_CLUSTER_SIZE, mps_fire::i, TCell::id(), iEtaHit, iPhiHit, M_PI, MAX_CONE_DIST, MAX_EOVERP, MAX_ET_THIRD_JET, MAX_PROBEJET_EMFRAC, MAX_TAGJET_ABSETA, MAX_TARGET_E, MAX_TRK_EME, MIN_CELL_E, MIN_DPHI_DIJETS, MIN_EOVERP, MIN_PROBEJET_ABSETA, MIN_TAGJET_ET, nEvents, phiSymCor, probeJetEmFrac, probeJetP4, refIEtaIPhi, edm::second(), TCell::SetE(), SUM_DEPTHS, SUM_SMALL_DEPTHS, sumDepths(), sumSmallDepths(), tagJetP4, targetE, targetEnergies, theCaloGeometry, USE_CONE_CLUSTERING, xTrkHcal, yTrkHcal, and zTrkHcal.

                                        {
  //  fChain->GetTree()->GetEntry(entry);
  GetEntry(entry);

  std::set<UInt_t> uniqueIds;  // for testing: check if there are duplicate cells   (AA)

  Bool_t acceptEvent = kTRUE;

  ++nEvents;

  if (!(nEvents % 100000))
    edm::LogVerbatim("HcalCalib") << "event: " << nEvents;

  h1_allTrkP->Fill(targetE);

  if (targetE < MIN_TARGET_E || targetE > MAX_TARGET_E)
    return kFALSE;
  ;

  // make local copy as the cells may be modified  due to phi/depth sum, phi corrections etc
  std::vector<TCell> selectCells;

  if (cells->GetSize() == 0)
    return kFALSE;

  if (CALIB_TYPE == "DI_JET" && probeJetEmFrac > 0.999)
    return kTRUE;

  for (Int_t i = 0; i < cells->GetSize(); ++i) {
    TCell* thisCell = (TCell*)cells->At(i);

    if (HcalDetId(thisCell->id()).subdet() == HcalOuter)
      continue;  // reject HO, make a switch!

    if (HcalDetId(thisCell->id()).subdet() != HcalBarrel && HcalDetId(thisCell->id()).subdet() != HcalEndcap &&
        HcalDetId(thisCell->id()).subdet() != HcalForward) {
      edm::LogWarning("HcalCalib") << "Unknown or wrong hcal subdetector: " << HcalDetId(thisCell->id()).subdet();
    }

    // Apply phi symmetry corrections if the flag is set
    if (APPLY_PHI_SYM_COR_FLAG)
      thisCell->SetE(phiSymCor[thisCell->id()] * thisCell->e());

    if (thisCell->e() > MIN_CELL_E)
      selectCells.push_back(*thisCell);
  }

  if (selectCells.empty()) {
    edm::LogWarning("HcalCalib") << "NO CELLS ABOVE THRESHOLD FOUND FOR TARGET!!!";
  }

  if (SUM_DEPTHS)
    sumDepths(selectCells);
  else if (SUM_SMALL_DEPTHS)
    sumSmallDepths(selectCells);  // depth 1,2 in twrs 15,16

  // most energetic tower (IsoTracks) or centroid of probe jet (DiJets)
  std::pair<Int_t, UInt_t> refPos;

  Int_t dEtaHitRef = 999;
  Int_t dPhiHitRef = 999;

  if (CALIB_TYPE == "ISO_TRACK") {
    Int_t iEtaMaxE;   // filled by reference in getIEtaIPhiForHighestE
    UInt_t iPhiMaxE;  //

    getIEtaIPhiForHighestE(selectCells, iEtaMaxE, iPhiMaxE);

    dEtaHitRef = iEtaMaxE - iEtaHit;
    dPhiHitRef = iPhiMaxE - iPhiHit;

    if (dPhiHitRef < -36)
      dPhiHitRef += 72;
    if (dPhiHitRef > 36)
      dPhiHitRef -= 72;

    if (iEtaHit * iEtaMaxE < 0) {
      if (dEtaHitRef < 0)
        dEtaHitRef += 1;
      if (dEtaHitRef > 0)
        dEtaHitRef -= 1;
    }

    if (abs(iEtaHit) < 16)
      h2_dHitRefBarrel->Fill(dEtaHitRef, dPhiHitRef);
    if (abs(iEtaHit) > 16 && abs(iEtaHit) < 25)
      h2_dHitRefEndcap->Fill(dEtaHitRef, dPhiHitRef);

    // --------------------------------------------------
    // Choice of cluster definition
    //
    // fixed size NxN clusters as specified in to config file
    if (!USE_CONE_CLUSTERING) {
      if (abs(iEtaMaxE) < 16 && HB_CLUSTER_SIZE == 3)
        filterCells3x3(selectCells, iEtaMaxE, iPhiMaxE);
      if (abs(iEtaMaxE) > 15 && HE_CLUSTER_SIZE == 3)
        filterCells3x3(selectCells, iEtaMaxE, iPhiMaxE);

      if (abs(iEtaMaxE) < 16 && HB_CLUSTER_SIZE == 5)
        filterCells5x5(selectCells, iEtaMaxE, iPhiMaxE);
      if (abs(iEtaMaxE) > 15 && HE_CLUSTER_SIZE == 5)
        filterCells5x5(selectCells, iEtaMaxE, iPhiMaxE);
    } else {
      //  calculate distance at hcal surface
      const GlobalPoint hitPositionHcal(xTrkHcal, yTrkHcal, zTrkHcal);
      filterCellsInCone(selectCells, hitPositionHcal, MAX_CONE_DIST, theCaloGeometry);
    }

    refPos.first = iEtaMaxE;
    refPos.second = iPhiMaxE;

  } else if (CALIB_TYPE == "DI_JET") {  // Apply selection cuts on DiJet events here

    if (etVetoJet > MAX_ET_THIRD_JET)
      acceptEvent = kFALSE;

    Float_t jetsDPhi = probeJetP4->DeltaPhi(*tagJetP4);
    if (fabs(jetsDPhi * 180.0 / M_PI) < MIN_DPHI_DIJETS)
      acceptEvent = kFALSE;

    if (probeJetEmFrac > MAX_PROBEJET_EMFRAC)
      acceptEvent = kFALSE;
    if (fabs(probeJetP4->Eta()) < MIN_PROBEJET_ABSETA)
      acceptEvent = kFALSE;
    if (fabs(tagJetP4->Eta()) > MAX_TAGJET_ABSETA)
      acceptEvent = kFALSE;
    if (fabs(tagJetP4->Et()) < MIN_TAGJET_ET)
      acceptEvent = kFALSE;

    if (acceptEvent) {
      Int_t iEtaMaxE;   // filled by reference in getIEtaIPhiForHighestE
      UInt_t iPhiMaxE;  //

      getIEtaIPhiForHighestE(selectCells, iEtaMaxE, iPhiMaxE);

      // The ref position for the jet is not used in the minimization at this time.
      // It will be needed if we attempt to do matrix inversion: then the question is
      // which value is better suited: the centroid of the jet or the hottest tower...

      //    refPos.first  = iEtaHit;
      //    refPos.second = iPhiHit;

      refPos.first = iEtaMaxE;
      refPos.second = iPhiMaxE;

      if (abs(iEtaMaxE) > 40)
        acceptEvent = kFALSE;  // for testing :  set as parameter (AA)
    }
  }

  if (COMBINE_PHI)
    combinePhi(selectCells);

  // fill the containers for the minimization prcedures
  std::vector<Float_t> energies;
  std::vector<UInt_t> ids;

  for (std::vector<TCell>::iterator i_it = selectCells.begin(); i_it != selectCells.end(); ++i_it) {
    // for testing : fill only unique id's

    if (uniqueIds.insert(i_it->id()).second) {
      energies.push_back(i_it->e());
      ids.push_back(i_it->id());
    }
  }

  if (CALIB_TYPE == "ISO_TRACK") {
    if (accumulate(energies.begin(), energies.end(), 0.0) / targetE < MIN_EOVERP)
      acceptEvent = kFALSE;
    if (accumulate(energies.begin(), energies.end(), 0.0) / targetE > MAX_EOVERP)
      acceptEvent = kFALSE;

    if (emEnergy > MAX_TRK_EME)
      acceptEvent = kFALSE;

    if (abs(dEtaHitRef) > 1 || abs(dPhiHitRef) > 1)
      acceptEvent = kFALSE;

    // Have to check if for |iEta|>20 (and neighboring region) the dPhiHitRef
    // should be relaxed to 2. The neighboring towers have dPhi=2...
  }

  h1_rawSumE->Fill(accumulate(energies.begin(), energies.end(), 0.0));

  // here we fill the information for the minimization procedure
  if (acceptEvent) {
    cellEnergies.push_back(energies);
    cellIds.push_back(ids);
    targetEnergies.push_back(targetE);
    refIEtaIPhi.push_back(refPos);

    if (abs(refPos.first) <= 10)
      h1_selTrkP_iEta10->Fill(targetE);
  }

  // Clean up
  energies.clear();
  ids.clear();
  selectCells.clear();

  return kTRUE;
}

ReadPhiSymCor()

Bool_t hcalCalib::ReadPhiSymCor

(

)

Definition at line 557 of file hcalCalib.cc.

References hcalRecHitTable_cff::depth, hcalRecHitTable_cff::detId, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, l1tTowerCalibrationProducer_cfi::iEta, mps_splice::line, PHI_SYM_COR_FILENAME, phiSymCor, heavyFlavorDQMFirstStep_cff::sdName, AlCaHLTBitMon_QueryRunRegistry::string, topo_, HcalTopology::valid(), and relativeConstraints::value.

Referenced by Begin().

                                {
  std::ifstream phiSymFile(PHI_SYM_COR_FILENAME.Data());

  if (!phiSymFile) {
    edm::LogWarning("HcalCalib") << "\nERROR: Can not find file with phi symmetry constants \""
                                 << PHI_SYM_COR_FILENAME.Data() << "\"";
    return kFALSE;
  }

  // assumes the format used in CSA08, first line is a comment

  Int_t iEta;
  UInt_t iPhi;
  UInt_t depth;
  TString sdName;
  UInt_t detId;

  Float_t value;
  HcalSubdetector sd;

  std::string line;

  while (getline(phiSymFile, line)) {
    if (line.empty() || line[0] == '#')
      continue;

    std::istringstream linestream(line);
    linestream >> iEta >> iPhi >> depth >> sdName >> value >> std::hex >> detId;

    if (sdName == "HB")
      sd = HcalBarrel;
    else if (sdName == "HE")
      sd = HcalEndcap;
    else if (sdName == "HO")
      sd = HcalOuter;
    else if (sdName == "HF")
      sd = HcalForward;
    else {
      edm::LogWarning("HcalCalib") << "\nInvalid detector name in phi symmetry constants file: " << sdName.Data()
                                   << "\nCheck file and rerun!\n";
      return kFALSE;
    }

    // check if the data is consistent

    if (HcalDetId(sd, iEta, iPhi, depth) != HcalDetId(detId)) {
      edm::LogWarning("HcalCalib")
          << "\nInconsistent info in phi symmetry file: subdet, iEta, iPhi, depth do not match rawId!\n";
      return kFALSE;
    }
    HcalDetId hId(detId);
    if (!topo_->valid(hId)) {
      edm::LogWarning("HcalCalib") << "\nInvalid DetId from: iEta=" << iEta << " iPhi=" << iPhi << " depth=" << depth
                                   << " subdet=" << sdName.Data() << " detId=" << detId << "\n";
      return kFALSE;
    }

    phiSymCor[HcalDetId(sd, iEta, iPhi, depth)] = value;
  }

  return kTRUE;
}

SetApplyPhiSymCorFlag()

void hcalCalib::SetApplyPhiSymCorFlag ( Bool_t  b )

inline

Definition at line 204 of file hcalCalib.h.

References APPLY_PHI_SYM_COR_FLAG, and b.

Referenced by HcalCalibrator::endJob().

{ APPLY_PHI_SYM_COR_FLAG = b; }

SetCalibAbsIEtaMax()

void hcalCalib::SetCalibAbsIEtaMax ( Int_t  i )

inline

Definition at line 200 of file hcalCalib.h.

References CALIB_ABS_IETA_MAX, and mps_fire::i.

Referenced by HcalCalibrator::endJob().

{ CALIB_ABS_IETA_MAX = i; }

SetCalibAbsIEtaMin()

void hcalCalib::SetCalibAbsIEtaMin ( Int_t  i )

inline

Definition at line 201 of file hcalCalib.h.

References CALIB_ABS_IETA_MIN, and mps_fire::i.

Referenced by HcalCalibrator::endJob().

{ CALIB_ABS_IETA_MIN = i; }

SetCalibMethod()

void hcalCalib::SetCalibMethod ( const TString &  s )

inline

Definition at line 193 of file hcalCalib.h.

References CALIB_METHOD, and alignCSCRings::s.

Referenced by HcalCalibrator::endJob().

{ CALIB_METHOD = s; }

SetCalibType()

void hcalCalib::SetCalibType ( const TString &  s )

inline

Definition at line 192 of file hcalCalib.h.

References CALIB_TYPE, and alignCSCRings::s.

Referenced by HcalCalibrator::endJob().

{ CALIB_TYPE = s; }

SetCaloGeometry()

void hcalCalib::SetCaloGeometry ( const CaloGeometry *  g, const HcalTopology *  topo  )

inline

Definition at line 214 of file hcalCalib.h.

References g, theCaloGeometry, and topo_.

Referenced by HcalCalibrator::endJob().

                                                                        {
    theCaloGeometry = g;
    topo_ = topo;
  }

SetCombinePhiFlag()

void hcalCalib::SetCombinePhiFlag ( Bool_t  b )

inline

Definition at line 187 of file hcalCalib.h.

References b, and COMBINE_PHI.

Referenced by HcalCalibrator::endJob().

{ COMBINE_PHI = b; }

SetConeMaxDist()

void hcalCalib::SetConeMaxDist ( Float_t  d )

inline

Definition at line 198 of file hcalCalib.h.

References ztail::d, and MAX_CONE_DIST.

Referenced by HcalCalibrator::endJob().

{ MAX_CONE_DIST = d; }

SetHbClusterSize()

void hcalCalib::SetHbClusterSize ( Int_t  i )

inline

Definition at line 194 of file hcalCalib.h.

References HB_CLUSTER_SIZE, and mps_fire::i.

Referenced by HcalCalibrator::endJob().

{ HB_CLUSTER_SIZE = i; }

SetHeClusterSize()

void hcalCalib::SetHeClusterSize ( Int_t  i )

inline

Definition at line 195 of file hcalCalib.h.

References HE_CLUSTER_SIZE, and mps_fire::i.

Referenced by HcalCalibrator::endJob().

{ HE_CLUSTER_SIZE = i; }

SetHistoFileName()

void hcalCalib::SetHistoFileName ( const TString &  filename )

inline

Definition at line 212 of file hcalCalib.h.

References corrVsCorr::filename, and HISTO_FILENAME.

Referenced by HcalCalibrator::endJob().

{ HISTO_FILENAME = filename; }

SetInputList()

void hcalCalib::SetInputList ( TList *  input )

inlineoverride

Definition at line 134 of file hcalCalib.h.

References input.

{ fInput = input; }

SetMaxEOverP()

void hcalCalib::SetMaxEOverP ( Float_t  e )

inline

Definition at line 190 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MAX_EOVERP.

Referenced by HcalCalibrator::endJob().

{ MAX_EOVERP = e; }

SetMaxEtThirdJet()

void hcalCalib::SetMaxEtThirdJet ( Float_t  et )

inline

Definition at line 202 of file hcalCalib.h.

References EgHLTOffHistBins_cfi::et, and MAX_ET_THIRD_JET.

Referenced by HcalCalibrator::endJob().

{ MAX_ET_THIRD_JET = et; }

SetMaxProbeJetEmFrac()

void hcalCalib::SetMaxProbeJetEmFrac ( Float_t  f )

inline

Definition at line 206 of file hcalCalib.h.

References f, and MAX_PROBEJET_EMFRAC.

Referenced by HcalCalibrator::endJob().

{ MAX_PROBEJET_EMFRAC = f; }

SetMaxTagJetAbsEta()

void hcalCalib::SetMaxTagJetAbsEta ( Float_t  e )

inline

Definition at line 208 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MAX_TAGJET_ABSETA.

Referenced by HcalCalibrator::endJob().

{ MAX_TAGJET_ABSETA = e; }

SetMaxTagJetEmFrac()

void hcalCalib::SetMaxTagJetEmFrac ( Float_t  f )

inline

Definition at line 207 of file hcalCalib.h.

References f, and MAX_TAGJET_EMFRAC.

Referenced by HcalCalibrator::endJob().

{ MAX_TAGJET_EMFRAC = f; }

SetMaxTargetE()

void hcalCalib::SetMaxTargetE ( Float_t  e )

inline

Definition at line 184 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MAX_TARGET_E.

Referenced by HcalCalibrator::endJob().

{ MAX_TARGET_E = e; }

SetMaxTrkEmE()

void hcalCalib::SetMaxTrkEmE ( Float_t  e )

inline

Definition at line 191 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MAX_TRK_EME.

Referenced by HcalCalibrator::endJob().

{ MAX_TRK_EME = e; }

SetMinCellE()

void hcalCalib::SetMinCellE ( Float_t  e )

inline

Definition at line 188 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MIN_CELL_E.

Referenced by HcalCalibrator::endJob().

{ MIN_CELL_E = e; }

SetMinDPhiDiJets()

void hcalCalib::SetMinDPhiDiJets ( Float_t  dphi )

inline

Definition at line 203 of file hcalCalib.h.

References MIN_DPHI_DIJETS.

Referenced by HcalCalibrator::endJob().

{ MIN_DPHI_DIJETS = dphi; }

SetMinEOverP()

void hcalCalib::SetMinEOverP ( Float_t  e )

inline

Definition at line 189 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MIN_EOVERP.

Referenced by HcalCalibrator::endJob().

{ MIN_EOVERP = e; }

SetMinProbeJetAbsEta()

void hcalCalib::SetMinProbeJetAbsEta ( Float_t  e )

inline

Definition at line 210 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MIN_PROBEJET_ABSETA.

Referenced by HcalCalibrator::endJob().

{ MIN_PROBEJET_ABSETA = e; }

SetMinTagJetEt()

void hcalCalib::SetMinTagJetEt ( Float_t  e )

inline

Definition at line 209 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MIN_TAGJET_ET.

Referenced by HcalCalibrator::endJob().

{ MIN_TAGJET_ET = e; }

SetMinTargetE()

void hcalCalib::SetMinTargetE ( Float_t  e )

inline

Definition at line 183 of file hcalCalib.h.

References MillePedeFileConverter_cfg::e, and MIN_TARGET_E.

Referenced by HcalCalibrator::endJob().

{ MIN_TARGET_E = e; }

SetObject()

void hcalCalib::SetObject ( TObject *  obj )

inlineoverride

Definition at line 133 of file hcalCalib.h.

References getGTfromDQMFile::obj.

{ fObject = obj; }

SetOption()

void hcalCalib::SetOption ( const char *  option )

inlineoverride

Definition at line 132 of file hcalCalib.h.

References fileinputsource_cfi::option.

{ fOption = option; }

SetOutputCorCoefFileName()

void hcalCalib::SetOutputCorCoefFileName ( const TString &  filename )

inline

Definition at line 211 of file hcalCalib.h.

References corrVsCorr::filename, and OUTPUT_COR_COEF_FILENAME.

Referenced by HcalCalibrator::endJob().

{ OUTPUT_COR_COEF_FILENAME = filename; }

SetPhiSymCorFileName()

void hcalCalib::SetPhiSymCorFileName ( const TString &  filename )

inline

Definition at line 205 of file hcalCalib.h.

References corrVsCorr::filename, and PHI_SYM_COR_FILENAME.

Referenced by HcalCalibrator::endJob().

{ PHI_SYM_COR_FILENAME = filename; }

SetSumDepthsFlag()

void hcalCalib::SetSumDepthsFlag ( Bool_t  b )

inline

Definition at line 185 of file hcalCalib.h.

References b, and SUM_DEPTHS.

Referenced by HcalCalibrator::endJob().

{ SUM_DEPTHS = b; }

SetSumSmallDepthsFlag()

void hcalCalib::SetSumSmallDepthsFlag ( Bool_t  b )

inline

Definition at line 186 of file hcalCalib.h.

References b, and SUM_SMALL_DEPTHS.

Referenced by HcalCalibrator::endJob().

{ SUM_SMALL_DEPTHS = b; }

SetUseConeClustering()

void hcalCalib::SetUseConeClustering ( Bool_t  b )

inline

Definition at line 197 of file hcalCalib.h.

References b, and USE_CONE_CLUSTERING.

Referenced by HcalCalibrator::endJob().

{ USE_CONE_CLUSTERING = b; }

Terminate()

void hcalCalib::Terminate ( )

override

Definition at line 316 of file hcalCalib.cc.

References funct::abs(), CALIB_ABS_IETA_MAX, CALIB_ABS_IETA_MIN, CALIB_METHOD, CALIB_TYPE, cellEnergies, cellIds, HcalDetId::depth(), GetCoefFromMtrxInvOfAve(), h1_allTrkP, h1_corResp, h1_corRespBarrel, h1_corRespEndcap, h1_corRespIEta, h1_numEventsTwrIEta, h1_rawResp, h1_rawRespBarrel, h1_rawRespEndcap, h1_rawSumE, h1_selTrkP_iEta10, h1_trkP, h2_dHitRefBarrel, h2_dHitRefEndcap, HcalBarrel, histoFile, mps_fire::i, l1tTowerCalibrationProducer_cfi::iEta, HcalDetId::ieta(), HcalDetId::ietaAbs(), iEtaCoefMap, HcalDetId::iphi(), MinL3AlgoUniv< IDdet >::iterate(), dqmiolumiharvest::j, makeTextFile(), refIEtaIPhi, findQualityFiles::size, solution, HcalDetId::subdet(), SUM_DEPTHS, SUM_SMALL_DEPTHS, targetEnergies, and trackerHitRTTI::vector.

                          {
  edm::LogVerbatim("HcalCalib") << "\n\nFinished reading the events.\n"
                                << "Number of input objects: " << cellIds.size()
                                << "\nPerforming minimization: depending on selected method can take some time...\n\n";

  for (std::vector<std::pair<Int_t, UInt_t> >::iterator it_rp = refIEtaIPhi.begin(); it_rp != refIEtaIPhi.end();
       ++it_rp) {
    Float_t weight = (abs(it_rp->first) < 21) ? 1.0 / 72.0 : 1.0 / 36.0;
    h1_numEventsTwrIEta->Fill(it_rp->first, weight);
  }

  if (CALIB_METHOD == "MATRIX_INV_OF_ETA_AVE") {
    GetCoefFromMtrxInvOfAve();
  } else if (CALIB_METHOD == "L3" || CALIB_METHOD == "L3_AND_MTRX_INV") {
    int eventWeight = 2;  // 2 is the default (try at some point 0,1,2,3)
    MinL3AlgoUniv<UInt_t>* thisL3Algo = new MinL3AlgoUniv<UInt_t>(eventWeight);
    int numIterations = 10;  // default 10

    solution = thisL3Algo->iterate(cellEnergies, cellIds, targetEnergies, numIterations);

    // in order to handle the case where sumDepths="false", but the flag to sum depths 1,2 in HB towers 15, 16
    // is set (sumSmallDepths) we create entries in "solution" to create equal correction here
    // for each encountered coef in depth one.

    if (!SUM_DEPTHS && SUM_SMALL_DEPTHS) {
      std::vector<UInt_t> idForSummedCells;

      for (std::map<UInt_t, Float_t>::iterator m_it = solution.begin(); m_it != solution.end(); ++m_it) {
        if (HcalDetId(m_it->first).ietaAbs() != 15 && HcalDetId(m_it->first).ietaAbs() != 16)
          continue;
        if (HcalDetId(m_it->first).subdet() != HcalBarrel)
          continue;
        if (HcalDetId(m_it->first).depth() == 1)
          idForSummedCells.push_back(HcalDetId(m_it->first));
      }

      for (std::vector<UInt_t>::iterator v_it = idForSummedCells.begin(); v_it != idForSummedCells.end(); ++v_it) {
        UInt_t addCoefId = HcalDetId(HcalBarrel, HcalDetId(*v_it).ieta(), HcalDetId(*v_it).iphi(), 2);
        solution[addCoefId] = solution[*v_it];
      }

    }  // end of special treatment for "sumSmallDepths" mode

    if (CALIB_METHOD == "L3_AND_MTRX_INV") {
      GetCoefFromMtrxInvOfAve();

      // loop over the solution from L3 and multiply by the additional factor from
      // the matrix inversion. Set coef outside of the valid calibration region =1.
      for (std::map<UInt_t, Float_t>::iterator it_s = solution.begin(); it_s != solution.end(); ++it_s) {
        Int_t iEtaSol = HcalDetId(it_s->first).ieta();
        if (abs(iEtaSol) < CALIB_ABS_IETA_MIN || abs(iEtaSol) > CALIB_ABS_IETA_MAX)
          it_s->second = 1.0;
        else
          it_s->second *= iEtaCoefMap[iEtaSol];
      }

    }  // if CALIB_METHOD=="L3_AND_MTRX_INV"

  }  // end of L3 or L3 + mtrx inversion minimization

  // done getting the constants -> write the formatted file
  makeTextFile();

  // fill some histograms
  Float_t rawResp = 0;
  Float_t corResp = 0;
  Int_t maxIEta = 0;
  Int_t minIEta = 999;

  for (unsigned int i = 0; i < cellEnergies.size(); ++i) {
    Int_t iEta;
    for (unsigned int j = 0; j < (cellEnergies[i]).size(); ++j) {
      iEta = HcalDetId(cellIds[i][j]).ieta();
      rawResp += (cellEnergies[i])[j];

      if (CALIB_METHOD == "L3_AND_MTRX_INV") {
        corResp += (solution[cellIds[i][j]] * cellEnergies[i][j]);
      } else if (CALIB_METHOD == "L3") {
        corResp += (solution[cellIds[i][j]] * (cellEnergies[i][j]));
      } else if (CALIB_METHOD == "MATRIX_INV_OF_ETA_AVE") {
        corResp += (iEtaCoefMap[iEta] * cellEnergies[i][j]);
      }

      if (maxIEta < abs(iEta))
        maxIEta = abs(iEta);
      if (minIEta > abs(iEta))
        minIEta = abs(iEta);
    }

    rawResp /= targetEnergies[i];
    corResp /= targetEnergies[i];

    // fill histograms based on iEta on ref point of the cluster (for now: hottest tower)
    // expect range |iEta|<=24 (to do: add flexibility for arbitrary range)

    if (CALIB_TYPE == "ISO_TRACK") {
      Int_t ind = refIEtaIPhi[i].first;
      (ind < 0) ? (ind += 24) : (ind += 23);
      if (ind >= 0 && ind < 48) {
        h1_corRespIEta[ind]->Fill(corResp);
      }

      // fill histograms for cases where all towers are in the barrel or endcap
      if (maxIEta < 25) {
        h1_rawResp->Fill(rawResp);
        h1_corResp->Fill(corResp);
      }
      if (maxIEta < 15) {
        h1_rawRespBarrel->Fill(rawResp);
        h1_corRespBarrel->Fill(corResp);
      } else if (maxIEta < 25 && minIEta > 16) {
        h1_rawRespEndcap->Fill(rawResp);
        h1_corRespEndcap->Fill(corResp);
      }
    }  // histograms for isotrack calibration

    else {
      // put jet plots here
    }

    rawResp = 0;
    corResp = 0;
    maxIEta = 0;
    minIEta = 999;
  }

  // save the histograms
  h1_trkP->Write();
  h1_allTrkP->Write();

  h1_selTrkP_iEta10->Write();

  h1_rawSumE->Write();
  h1_rawResp->Write();
  h1_corResp->Write();
  h1_rawRespBarrel->Write();
  h1_corRespBarrel->Write();
  h1_rawRespEndcap->Write();
  h1_corRespEndcap->Write();
  h1_numEventsTwrIEta->Write();
  h2_dHitRefBarrel->Write();
  h2_dHitRefEndcap->Write();
  for (Int_t i = 0; i < 48; ++i) {
    h1_corRespIEta[i]->Write();
  }

  histoFile->Write();
  histoFile->Close();

  edm::LogVerbatim("HcalCalib") << "\n Finished calibration.\n ";

}  // end of Terminate()

Version()

Int_t hcalCalib::Version ( ) const

inlineoverride

Definition at line 123 of file hcalCalib.h.

{ return 2; }

Member Data Documentation

APPLY_PHI_SYM_COR_FLAG

Bool_t hcalCalib::APPLY_PHI_SYM_COR_FLAG

Definition at line 175 of file hcalCalib.h.

Referenced by Begin(), Process(), and SetApplyPhiSymCorFlag().

b_cells

TBranch* hcalCalib::b_cells

Definition at line 77 of file hcalCalib.h.

Referenced by Init().

b_emEnergy

TBranch* hcalCalib::b_emEnergy

Definition at line 78 of file hcalCalib.h.

Referenced by Init().

b_etVetoJet

TBranch* hcalCalib::b_etVetoJet

Definition at line 80 of file hcalCalib.h.

Referenced by Init().

b_eventNumber

TBranch* hcalCalib::b_eventNumber

Definition at line 73 of file hcalCalib.h.

Referenced by Init().

b_iEtaHit

TBranch* hcalCalib::b_iEtaHit

Definition at line 75 of file hcalCalib.h.

Referenced by Init().

b_iPhiHit

TBranch* hcalCalib::b_iPhiHit

Definition at line 76 of file hcalCalib.h.

Referenced by Init().

b_probeJetEmFrac

TBranch* hcalCalib::b_probeJetEmFrac

Definition at line 90 of file hcalCalib.h.

Referenced by Init().

b_probeJetP4

TBranch* hcalCalib::b_probeJetP4

Definition at line 93 of file hcalCalib.h.

Referenced by Init().

b_runNumber

TBranch* hcalCalib::b_runNumber

Definition at line 74 of file hcalCalib.h.

Referenced by Init().

b_tagJetEmFrac

TBranch* hcalCalib::b_tagJetEmFrac

Definition at line 89 of file hcalCalib.h.

Referenced by Init().

b_tagJetP4

TBranch* hcalCalib::b_tagJetP4

Definition at line 92 of file hcalCalib.h.

Referenced by Init().

b_targetE

TBranch* hcalCalib::b_targetE

Definition at line 79 of file hcalCalib.h.

Referenced by Init().

b_xTrkEcal

TBranch* hcalCalib::b_xTrkEcal

Definition at line 85 of file hcalCalib.h.

Referenced by Init().

b_xTrkHcal

TBranch* hcalCalib::b_xTrkHcal

Definition at line 82 of file hcalCalib.h.

Referenced by Init().

b_yTrkEcal

TBranch* hcalCalib::b_yTrkEcal

Definition at line 86 of file hcalCalib.h.

Referenced by Init().

b_yTrkHcal

TBranch* hcalCalib::b_yTrkHcal

Definition at line 83 of file hcalCalib.h.

Referenced by Init().

b_zTrkEcal

TBranch* hcalCalib::b_zTrkEcal

Definition at line 87 of file hcalCalib.h.

Referenced by Init().

b_zTrkHcal

TBranch* hcalCalib::b_zTrkHcal

Definition at line 84 of file hcalCalib.h.

Referenced by Init().

CALIB_ABS_IETA_MAX

Int_t hcalCalib::CALIB_ABS_IETA_MAX

Definition at line 161 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), makeTextFile(), SetCalibAbsIEtaMax(), and Terminate().

CALIB_ABS_IETA_MIN

Int_t hcalCalib::CALIB_ABS_IETA_MIN

Definition at line 162 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), makeTextFile(), SetCalibAbsIEtaMin(), and Terminate().

CALIB_METHOD

TString hcalCalib::CALIB_METHOD

Definition at line 172 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), makeTextFile(), SetCalibMethod(), and Terminate().

CALIB_TYPE

TString hcalCalib::CALIB_TYPE

Definition at line 171 of file hcalCalib.h.

Referenced by Begin(), Process(), SetCalibType(), and Terminate().

cellEnergies

std::vector<std::vector<Float_t> > hcalCalib::cellEnergies

Definition at line 226 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), Process(), and Terminate().

cellIds

std::vector<std::vector<UInt_t> > hcalCalib::cellIds

Definition at line 227 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), Process(), and Terminate().

cells

TClonesArray* hcalCalib::cells

Definition at line 54 of file hcalCalib.h.

Referenced by Init(), and Process().

COMBINE_PHI

Bool_t hcalCalib::COMBINE_PHI

Definition at line 153 of file hcalCalib.h.

Referenced by makeTextFile(), Process(), and SetCombinePhiFlag().

emEnergy

Float_t hcalCalib::emEnergy

Definition at line 55 of file hcalCalib.h.

Referenced by Init(), and Process().

etVetoJet

Float_t hcalCalib::etVetoJet

Definition at line 57 of file hcalCalib.h.

Referenced by Init(), and Process().

eventNumber

UInt_t hcalCalib::eventNumber

pointer to the analyzed TTree or TChain

Definition at line 50 of file hcalCalib.h.

Referenced by Init().

fChain

TTree* hcalCalib::fChain

Definition at line 48 of file hcalCalib.h.

Referenced by GetEntry(), and Init().

h1_allTrkP

TH1F* hcalCalib::h1_allTrkP

Definition at line 101 of file hcalCalib.h.

Referenced by Begin(), Process(), and Terminate().

h1_corResp

TH1F* hcalCalib::h1_corResp

Definition at line 107 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_corRespBarrel

TH1F* hcalCalib::h1_corRespBarrel

Definition at line 109 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_corRespEndcap

TH1F* hcalCalib::h1_corRespEndcap

Definition at line 111 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_corRespIEta

TH1F* hcalCalib::h1_corRespIEta[48]

Definition at line 119 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_numEventsTwrIEta

TH1F* hcalCalib::h1_numEventsTwrIEta

Definition at line 112 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_rawResp

TH1F* hcalCalib::h1_rawResp

Definition at line 106 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_rawRespBarrel

TH1F* hcalCalib::h1_rawRespBarrel

Definition at line 108 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_rawRespEndcap

TH1F* hcalCalib::h1_rawRespEndcap

Definition at line 110 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h1_rawSumE

TH1F* hcalCalib::h1_rawSumE

Definition at line 105 of file hcalCalib.h.

Referenced by Begin(), Process(), and Terminate().

h1_selTrkP_iEta10

TH1F* hcalCalib::h1_selTrkP_iEta10

Definition at line 103 of file hcalCalib.h.

Referenced by Begin(), Process(), and Terminate().

h1_trkP

TH1F* hcalCalib::h1_trkP

Definition at line 100 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

h2_dHitRefBarrel

TH2F* hcalCalib::h2_dHitRefBarrel

Definition at line 114 of file hcalCalib.h.

Referenced by Begin(), Process(), and Terminate().

h2_dHitRefEndcap

TH2F* hcalCalib::h2_dHitRefEndcap

Definition at line 115 of file hcalCalib.h.

Referenced by Begin(), Process(), and Terminate().

HB_CLUSTER_SIZE

Int_t hcalCalib::HB_CLUSTER_SIZE

Definition at line 155 of file hcalCalib.h.

Referenced by Process(), and SetHbClusterSize().

HE_CLUSTER_SIZE

Int_t hcalCalib::HE_CLUSTER_SIZE

Definition at line 156 of file hcalCalib.h.

Referenced by Process(), and SetHeClusterSize().

HISTO_FILENAME

TString hcalCalib::HISTO_FILENAME

Definition at line 178 of file hcalCalib.h.

Referenced by Begin(), and SetHistoFileName().

histoFile

TFile* hcalCalib::histoFile

Definition at line 97 of file hcalCalib.h.

Referenced by Begin(), and Terminate().

iEtaCoefMap

std::map<Int_t, Float_t> hcalCalib::iEtaCoefMap

Definition at line 235 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), makeTextFile(), and Terminate().

iEtaHit

Int_t hcalCalib::iEtaHit

Definition at line 52 of file hcalCalib.h.

Referenced by Init(), and Process().

iPhiHit

UInt_t hcalCalib::iPhiHit

Definition at line 53 of file hcalCalib.h.

Referenced by Init(), and Process().

MAX_CONE_DIST

Float_t hcalCalib::MAX_CONE_DIST

Definition at line 159 of file hcalCalib.h.

Referenced by Process(), and SetConeMaxDist().

MAX_EOVERP

Float_t hcalCalib::MAX_EOVERP

Definition at line 145 of file hcalCalib.h.

Referenced by Process(), and SetMaxEOverP().

MAX_ET_THIRD_JET

Float_t hcalCalib::MAX_ET_THIRD_JET

Definition at line 148 of file hcalCalib.h.

Referenced by Process(), and SetMaxEtThirdJet().

MAX_PROBEJET_EMFRAC

Float_t hcalCalib::MAX_PROBEJET_EMFRAC

Definition at line 164 of file hcalCalib.h.

Referenced by Process(), and SetMaxProbeJetEmFrac().

MAX_TAGJET_ABSETA

Float_t hcalCalib::MAX_TAGJET_ABSETA

Definition at line 166 of file hcalCalib.h.

Referenced by Process(), and SetMaxTagJetAbsEta().

MAX_TAGJET_EMFRAC

Float_t hcalCalib::MAX_TAGJET_EMFRAC

Definition at line 165 of file hcalCalib.h.

Referenced by SetMaxTagJetEmFrac().

MAX_TARGET_E

Float_t hcalCalib::MAX_TARGET_E

Definition at line 141 of file hcalCalib.h.

Referenced by Process(), and SetMaxTargetE().

MAX_TRK_EME

Float_t hcalCalib::MAX_TRK_EME

Definition at line 146 of file hcalCalib.h.

Referenced by Process(), and SetMaxTrkEmE().

MIN_CELL_E

Float_t hcalCalib::MIN_CELL_E

Definition at line 143 of file hcalCalib.h.

Referenced by Process(), and SetMinCellE().

MIN_DPHI_DIJETS

Float_t hcalCalib::MIN_DPHI_DIJETS

Definition at line 149 of file hcalCalib.h.

Referenced by Process(), and SetMinDPhiDiJets().

MIN_EOVERP

Float_t hcalCalib::MIN_EOVERP

Definition at line 144 of file hcalCalib.h.

Referenced by Process(), and SetMinEOverP().

MIN_PROBEJET_ABSETA

Float_t hcalCalib::MIN_PROBEJET_ABSETA

Definition at line 169 of file hcalCalib.h.

Referenced by Process(), and SetMinProbeJetAbsEta().

MIN_TAGJET_ET

Float_t hcalCalib::MIN_TAGJET_ET

Definition at line 167 of file hcalCalib.h.

Referenced by Process(), and SetMinTagJetEt().

MIN_TARGET_E

Float_t hcalCalib::MIN_TARGET_E

Definition at line 140 of file hcalCalib.h.

Referenced by SetMinTargetE().

nEvents

UInt_t hcalCalib::nEvents

Definition at line 95 of file hcalCalib.h.

Referenced by Begin(), looper.Looper::loop(), and Process().

OUTPUT_COR_COEF_FILENAME

TString hcalCalib::OUTPUT_COR_COEF_FILENAME

Definition at line 177 of file hcalCalib.h.

Referenced by makeTextFile(), and SetOutputCorCoefFileName().

PHI_SYM_COR_FILENAME

TString hcalCalib::PHI_SYM_COR_FILENAME

Definition at line 174 of file hcalCalib.h.

Referenced by ReadPhiSymCor(), and SetPhiSymCorFileName().

phiSymCor

std::map<UInt_t, Float_t> hcalCalib::phiSymCor

Definition at line 231 of file hcalCalib.h.

Referenced by Process(), and ReadPhiSymCor().

probeJetEmFrac

Float_t hcalCalib::probeJetEmFrac

Definition at line 70 of file hcalCalib.h.

Referenced by Init(), and Process().

probeJetP4

TLorentzVector* hcalCalib::probeJetP4

Definition at line 67 of file hcalCalib.h.

Referenced by Init(), and Process().

refIEtaIPhi

std::vector<std::pair<Int_t, UInt_t> > hcalCalib::refIEtaIPhi

Definition at line 228 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), Process(), and Terminate().

runNumber

UInt_t hcalCalib::runNumber

Definition at line 51 of file hcalCalib.h.

Referenced by Init(), and TH2PolyOfflineMaps.TH2PolyOfflineMaps::PrintTrackerMaps().

solution

std::map<UInt_t, Float_t> hcalCalib::solution

Definition at line 233 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), makeTextFile(), and Terminate().

SUM_DEPTHS

Bool_t hcalCalib::SUM_DEPTHS

Definition at line 151 of file hcalCalib.h.

Referenced by makeTextFile(), Process(), SetSumDepthsFlag(), and Terminate().

SUM_SMALL_DEPTHS

Bool_t hcalCalib::SUM_SMALL_DEPTHS

Definition at line 152 of file hcalCalib.h.

Referenced by Process(), SetSumSmallDepthsFlag(), and Terminate().

tagJetEmFrac

Float_t hcalCalib::tagJetEmFrac

Definition at line 69 of file hcalCalib.h.

Referenced by Init().

tagJetP4

TLorentzVector* hcalCalib::tagJetP4

Definition at line 66 of file hcalCalib.h.

Referenced by Init(), and Process().

targetE

Float_t hcalCalib::targetE

Definition at line 56 of file hcalCalib.h.

Referenced by Init(), and Process().

targetEnergies

std::vector<Float_t> hcalCalib::targetEnergies

Definition at line 229 of file hcalCalib.h.

Referenced by GetCoefFromMtrxInvOfAve(), Process(), and Terminate().

theCaloGeometry

const CaloGeometry* hcalCalib::theCaloGeometry

Definition at line 180 of file hcalCalib.h.

Referenced by Process(), and SetCaloGeometry().

topo_

const HcalTopology* hcalCalib::topo_

Definition at line 181 of file hcalCalib.h.

Referenced by makeTextFile(), ReadPhiSymCor(), and SetCaloGeometry().

USE_CONE_CLUSTERING

Bool_t hcalCalib::USE_CONE_CLUSTERING

Definition at line 158 of file hcalCalib.h.

Referenced by Process(), and SetUseConeClustering().

xTrkEcal

Float_t hcalCalib::xTrkEcal

Definition at line 62 of file hcalCalib.h.

Referenced by Init().

xTrkHcal

Float_t hcalCalib::xTrkHcal

Definition at line 59 of file hcalCalib.h.

Referenced by Init(), and Process().

yTrkEcal

Float_t hcalCalib::yTrkEcal

Definition at line 63 of file hcalCalib.h.

Referenced by Init().

yTrkHcal

Float_t hcalCalib::yTrkHcal

Definition at line 60 of file hcalCalib.h.

Referenced by Init(), and Process().

zTrkEcal

Float_t hcalCalib::zTrkEcal

Definition at line 64 of file hcalCalib.h.

Referenced by Init().

zTrkHcal

Float_t hcalCalib::zTrkHcal

Definition at line 61 of file hcalCalib.h.

Referenced by Init(), and Process().