MaterialBudgetHcalHistos Class Reference

#include <MaterialBudgetHcalHistos.h>

Public Member Functions
void	fillBeginJob (const DDCompactView &)
void	fillEndTrack ()
void	fillPerStep (const G4Step *)
void	fillStartTrack (const G4Track *)
	MaterialBudgetHcalHistos (const edm::ParameterSet &p)
virtual	~MaterialBudgetHcalHistos ()

Private Member Functions
void	book ()
void	fillHisto (int ii)
void	fillLayer ()
std::vector< double >	getDDDArray (const std::string &str, const DDsvalues_type &sv)
std::vector< std::string >	getNames (DDFilteredView &fv)
void	hend ()
bool	isItEC (std::string)
bool	isItHF (const G4VTouchable *)
bool	isSensitive (std::string)

Private Attributes
int	binEta
int	binPhi
double	eta
double	etaHigh
double	etaLow
bool	fillHistos
std::vector< int >	hfLevels
std::vector< std::string >	hfNames
int	id
double	intLen
std::vector< double >	intLength
int	layer
std::vector< std::string >	matList
double	maxEta
TProfile *	me100 [maxSet]
TProfile2D *	me1000 [maxSet]
TProfile2D *	me1100 [maxSet]
TH2F *	me1200 [maxSet]
TH1F *	me1300 [maxSet2]
TH2F *	me1400 [maxSet2]
TProfile *	me1500 [maxSet2]
TProfile *	me200 [maxSet]
TProfile *	me300 [maxSet]
TH1F *	me400 [maxSet]
TProfile *	me500 [maxSet]
TProfile *	me600 [maxSet]
TProfile *	me700 [maxSet]
TH1F *	me800 [maxSet]
TProfile2D *	me900 [maxSet]
int	nlayHB
int	nlayHE
int	nlayHF
int	nlayHO
double	phi
bool	printSum
double	radLen
std::vector< double >	radLength
std::vector< std::string >	sensitiveEC
std::vector< std::string >	sensitives
double	stepLen
std::vector< double >	stepLength
int	steps

Static Private Attributes
static const int	maxSet = 25
static const int	maxSet2 = 9

Detailed Description

Definition at line 19 of file MaterialBudgetHcalHistos.h.

Constructor & Destructor Documentation

MaterialBudgetHcalHistos()

MaterialBudgetHcalHistos::MaterialBudgetHcalHistos

(

const edm::ParameterSet &

p

)

Definition at line 19 of file MaterialBudgetHcalHistos.cc.

                                                                           {
  binEta = p.getUntrackedParameter<int>("NBinEta", 260);
  binPhi = p.getUntrackedParameter<int>("NBinPhi", 180);
  maxEta = p.getUntrackedParameter<double>("MaxEta", 5.2);
  etaLow = p.getUntrackedParameter<double>("EtaLow", -5.2);
  etaHigh = p.getUntrackedParameter<double>("EtaHigh", 5.2);
  fillHistos = p.getUntrackedParameter<bool>("FillHisto", true);
  printSum = p.getUntrackedParameter<bool>("PrintSummary", false);
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: FillHisto : " << fillHistos << " PrintSummary "
                                     << printSum << " == Eta plot: NX " << binEta << " Range " << -maxEta << ":"
                                     << maxEta << " Phi plot: NX " << binPhi << " Range " << -pi << ":" << pi
                                     << " (Eta limit " << etaLow << ":" << etaHigh << ")";
  if (fillHistos)
    book();
}

References binEta, binPhi, book(), etaHigh, etaLow, fillHistos, maxEta, AlCaHLTBitMon_ParallelJobs::p, pi, and printSum.

~MaterialBudgetHcalHistos()

virtual MaterialBudgetHcalHistos::~MaterialBudgetHcalHistos ( )

inlinevirtual

Definition at line 22 of file MaterialBudgetHcalHistos.h.

{ hend(); }

References hend().

Member Function Documentation

book()

void MaterialBudgetHcalHistos::book ( )

private

Definition at line 246 of file MaterialBudgetHcalHistos.cc.

                                    {
  // Book histograms
  edm::Service<TFileService> tfile;
 
  if (!tfile.isAvailable())
    throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                      << "please add it to config file";
 
  double maxPhi = pi;
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Booking user histos === with " << binEta
                                     << " bins in eta from " << -maxEta << " to " << maxEta << " and " << binPhi
                                     << " bins in phi from " << -maxPhi << " to " << maxPhi;
 
  std::string iter;
  // total X0
  for (int i = 0; i < maxSet; i++) {
    iter = std::to_string(i);
    me100[i] = tfile->make<TProfile>(
        std::to_string(i + 100).c_str(), ("MB(X0) prof Eta in region " + iter).c_str(), binEta, -maxEta, maxEta);
    me200[i] = tfile->make<TProfile>(
        std::to_string(i + 200).c_str(), ("MB(L0) prof Eta in region " + iter).c_str(), binEta, -maxEta, maxEta);
    me300[i] = tfile->make<TProfile>(
        std::to_string(i + 300).c_str(), ("MB(Step) prof Eta in region " + iter).c_str(), binEta, -maxEta, maxEta);
    me400[i] =
        tfile->make<TH1F>(std::to_string(i + 400).c_str(), ("Eta in region " + iter).c_str(), binEta, -maxEta, maxEta);
    me500[i] = tfile->make<TProfile>(
        std::to_string(i + 500).c_str(), ("MB(X0) prof Ph in region " + iter).c_str(), binPhi, -maxPhi, maxPhi);
    me600[i] = tfile->make<TProfile>(
        std::to_string(i + 600).c_str(), ("MB(L0) prof Ph in region " + iter).c_str(), binPhi, -maxPhi, maxPhi);
    me700[i] = tfile->make<TProfile>(
        std::to_string(i + 700).c_str(), ("MB(Step) prof Ph in region " + iter).c_str(), binPhi, -maxPhi, maxPhi);
    me800[i] =
        tfile->make<TH1F>(std::to_string(i + 800).c_str(), ("Phi in region " + iter).c_str(), binPhi, -maxPhi, maxPhi);
    me900[i] = tfile->make<TProfile2D>(std::to_string(i + 900).c_str(),
                                       ("MB(X0) prof Eta Phi in region " + iter).c_str(),
                                       binEta / 2,
                                       -maxEta,
                                       maxEta,
                                       binPhi / 2,
                                       -maxPhi,
                                       maxPhi);
    me1000[i] = tfile->make<TProfile2D>(std::to_string(i + 1000).c_str(),
                                        ("MB(L0) prof Eta Phi in region " + iter).c_str(),
                                        binEta / 2,
                                        -maxEta,
                                        maxEta,
                                        binPhi / 2,
                                        -maxPhi,
                                        maxPhi);
    me1100[i] = tfile->make<TProfile2D>(std::to_string(i + 1100).c_str(),
                                        ("MB(Step) prof Eta Phi in region " + iter).c_str(),
                                        binEta / 2,
                                        -maxEta,
                                        maxEta,
                                        binPhi / 2,
                                        -maxPhi,
                                        maxPhi);
    me1200[i] = tfile->make<TH2F>(std::to_string(i + 1200).c_str(),
                                  ("Eta vs Phi in region " + iter).c_str(),
                                  binEta / 2,
                                  -maxEta,
                                  maxEta,
                                  binPhi / 2,
                                  -maxPhi,
                                  maxPhi);
  }
  for (int i = 0; i < maxSet2; i++) {
    iter = std::to_string(i);
    me1300[i] = tfile->make<TH1F>(std::to_string(i + 1300).c_str(),
                                  ("Events with layers Hit (0 all, 1 HB, ..) for " + iter).c_str(),
                                  binEta,
                                  -maxEta,
                                  maxEta);
    me1400[i] = tfile->make<TH2F>(std::to_string(i + 1400).c_str(),
                                  ("Eta vs Phi for layers hit in " + iter).c_str(),
                                  binEta / 2,
                                  -maxEta,
                                  maxEta,
                                  binPhi / 2,
                                  -maxPhi,
                                  maxPhi);
    me1500[i] = tfile->make<TProfile>(std::to_string(i + 1500).c_str(),
                                      ("Number of layers crossed (0 all, 1 HB, ..) for " + iter).c_str(),
                                      binEta,
                                      -maxEta,
                                      maxEta);
  }
 
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Booking user histos done ===";
}

References binEta, binPhi, mps_fire::i, maxEta, HLT_2018_cff::maxPhi, maxSet, maxSet2, me100, me1000, me1100, me1200, me1300, me1400, me1500, me200, me300, me400, me500, me600, me700, me800, me900, pi, AlCaHLTBitMon_QueryRunRegistry::string, and compare::tfile.

Referenced by MaterialBudgetHcalHistos().

fillBeginJob()

void MaterialBudgetHcalHistos::fillBeginJob

(

const DDCompactView &

cpv

)

Definition at line 35 of file MaterialBudgetHcalHistos.cc.

                                                                    {
  if (fillHistos) {
    std::string attribute = "ReadOutName";
    std::string value = "HcalHits";
    DDSpecificsMatchesValueFilter filter1{DDValue(attribute, value, 0)};
    DDFilteredView fv1(cpv, filter1);
    sensitives = getNames(fv1);
    edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Names to be tested for " << attribute << " = "
                                       << value << " has " << sensitives.size() << " elements";
    for (unsigned int i = 0; i < sensitives.size(); i++)
      edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: sensitives[" << i << "] = " << sensitives[i];
 
    attribute = "Volume";
    value = "HF";
    DDSpecificsMatchesValueFilter filter2{DDValue(attribute, value, 0)};
    DDFilteredView fv2(cpv, filter2);
    hfNames = getNames(fv2);
    fv2.firstChild();
    DDsvalues_type sv(fv2.mergedSpecifics());
    std::vector<double> temp = getDDDArray("Levels", sv);
    edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Names to be tested for " << attribute << " = "
                                       << value << " has " << hfNames.size() << " elements";
    for (unsigned int i = 0; i < hfNames.size(); i++) {
      int level = static_cast<int>(temp[i]);
      hfLevels.push_back(level);
      edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos:  HF[" << i << "] = " << hfNames[i]
                                         << " at level " << hfLevels[i];
    }
 
    std::string ecalRO[2] = {"EcalHitsEB", "EcalHitsEE"};
    attribute = "ReadOutName";
    for (int k = 0; k < 2; k++) {
      value = ecalRO[k];
      DDSpecificsMatchesValueFilter filter3{DDValue(attribute, value, 0)};
      DDFilteredView fv3(cpv, filter3);
      std::vector<std::string> senstmp = getNames(fv3);
      edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Names to be tested for " << attribute << " = "
                                         << value << " has " << senstmp.size() << " elements";
      for (unsigned int i = 0; i < senstmp.size(); i++)
        sensitiveEC.push_back(senstmp[i]);
    }
    for (unsigned int i = 0; i < sensitiveEC.size(); i++)
      edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos:sensitiveEC[" << i << "] = " << sensitiveEC[i];
  }
}

References fillHistos, DDFilteredView::firstChild(), getDDDArray(), getNames(), hfLevels, hfNames, mps_fire::i, dqmdumpme::k, personalPlayback::level, DDFilteredView::mergedSpecifics(), sensitiveEC, sensitives, AlCaHLTBitMon_QueryRunRegistry::string, pfDeepBoostedJetPreprocessParams_cfi::sv, and groupFilesInBlocks::temp.

Referenced by MaterialBudgetHcal::update().

fillEndTrack()

void MaterialBudgetHcalHistos::fillEndTrack

(

)

Definition at line 231 of file MaterialBudgetHcalHistos.cc.

                                            {
  edm::LogVerbatim("MaterialBudget") << "Number of layers hit in HB:" << nlayHB << " HE:" << nlayHE << " HO:" << nlayHO
                                     << " HF:" << nlayHF;
  if (fillHistos) {
    fillHisto(maxSet - 1);
    fillLayer();
  }
  if (printSum) {
    for (unsigned int ii = 0; ii < matList.size(); ii++) {
      edm::LogVerbatim("MaterialBudget") << matList[ii] << "\t" << stepLength[ii] << "\t" << radLength[ii] << "\t"
                                         << intLength[ii];
    }
  }
}

References fillHisto(), fillHistos, fillLayer(), cuy::ii, intLength, matList, maxSet, nlayHB, nlayHE, nlayHF, nlayHO, printSum, radLength, and stepLength.

Referenced by MaterialBudgetHcal::update().

fillHisto()

void MaterialBudgetHcalHistos::fillHisto ( int  ii )

private

Definition at line 337 of file MaterialBudgetHcalHistos.cc.

                                               {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos:FillHisto called with index " << ii
                                     << " integrated  step " << stepLen << " X0 " << radLen << " Lamda " << intLen;
#endif
  if (ii >= 0 && ii < maxSet) {
    me100[ii]->Fill(eta, radLen);
    me200[ii]->Fill(eta, intLen);
    me300[ii]->Fill(eta, stepLen);
    me400[ii]->Fill(eta);
 
    if (eta >= etaLow && eta <= etaHigh) {
      me500[ii]->Fill(phi, radLen);
      me600[ii]->Fill(phi, intLen);
      me700[ii]->Fill(phi, stepLen);
      me800[ii]->Fill(phi);
    }
 
    me900[ii]->Fill(eta, phi, radLen);
    me1000[ii]->Fill(eta, phi, intLen);
    me1100[ii]->Fill(eta, phi, stepLen);
    me1200[ii]->Fill(eta, phi);
  }
}

References eta, etaHigh, etaLow, cuy::ii, intLen, maxSet, me100, me1000, me1100, me1200, me200, me300, me400, me500, me600, me700, me800, me900, phi, radLen, and stepLen.

Referenced by fillEndTrack(), and fillPerStep().

fillLayer()

void MaterialBudgetHcalHistos::fillLayer ( )

private

Definition at line 362 of file MaterialBudgetHcalHistos.cc.

                                         {
  me1300[0]->Fill(eta);
  me1400[0]->Fill(eta, phi);
  if (nlayHB > 0) {
    me1300[1]->Fill(eta);
    me1400[1]->Fill(eta, phi);
  }
  if (nlayHB >= 16) {
    me1300[2]->Fill(eta);
    me1400[2]->Fill(eta, phi);
  }
  if (nlayHE > 0) {
    me1300[3]->Fill(eta);
    me1400[3]->Fill(eta, phi);
  }
  if (nlayHE >= 16) {
    me1300[4]->Fill(eta);
    me1400[4]->Fill(eta, phi);
  }
  if (nlayHO > 0) {
    me1300[5]->Fill(eta);
    me1400[5]->Fill(eta, phi);
  }
  if (nlayHO >= 2) {
    me1300[6]->Fill(eta);
    me1400[6]->Fill(eta, phi);
  }
  if (nlayHF > 0) {
    me1300[7]->Fill(eta);
    me1400[7]->Fill(eta, phi);
  }
  if (nlayHB > 0 || nlayHE > 0 || (nlayHF > 0 && std::abs(eta) > 3.0)) {
    me1300[8]->Fill(eta);
    me1400[8]->Fill(eta, phi);
  }
  me1500[0]->Fill(eta, (double)(nlayHB + nlayHO + nlayHE + nlayHF));
  me1500[1]->Fill(eta, (double)(nlayHB));
  me1500[2]->Fill(eta, (double)(nlayHE));
  me1500[4]->Fill(eta, (double)(nlayHF));
}

References funct::abs(), eta, me1300, me1400, me1500, nlayHB, nlayHE, nlayHF, nlayHO, and phi.

Referenced by fillEndTrack().

fillPerStep()

void MaterialBudgetHcalHistos::fillPerStep

(

const G4Step *

aStep

)

Definition at line 108 of file MaterialBudgetHcalHistos.cc.

                                                              {
  G4Material* material = aStep->GetPreStepPoint()->GetMaterial();
  double step = aStep->GetStepLength();
  double radl = material->GetRadlen();
  double intl = material->GetNuclearInterLength();
  double density = material->GetDensity() / (g / cm3);
 
  int idOld = id;
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  std::string name = touch->GetVolume(0)->GetName();
  const std::string& matName = material->GetName();
  if (printSum) {
    bool found = false;
    for (unsigned int ii = 0; ii < matList.size(); ii++) {
      if (matList[ii] == matName) {
        stepLength[ii] += step;
        radLength[ii] += (step / radl);
        intLength[ii] += (step / intl);
        found = true;
        break;
      }
    }
    if (!found) {
      matList.push_back(matName);
      stepLength.push_back(step);
      radLength.push_back(step / radl);
      intLength.push_back(step / intl);
    }
    edm::LogVerbatim("MaterialBudget") << "Volume " << name << " id " << id << ":" << idOld << " Step " << step
                                       << " Material " << matName << " Old Length " << stepLen << " X0 " << step / radl
                                       << ":" << radLen << " Lambda " << step / intl << ":" << intLen;
  } else {
    edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Step at " << name << " id " << id << ":" << idOld
                                       << " Length " << step << " in " << matName << " of density " << density
                                       << " g/cc; Radiation Length " << radl << " mm; Interaction Length " << intl
                                       << " mm\n"
                                       << "                          Position "
                                       << aStep->GetPreStepPoint()->GetPosition() << " Cylindrical R "
                                       << aStep->GetPreStepPoint()->GetPosition().perp() << " Length (so far) "
                                       << stepLen << " L/X0 " << step / radl << "/" << radLen << " L/Lambda "
                                       << step / intl << "/" << intLen;
  }
 
  int det = 0, lay = 0;
  if (fillHistos) {
    if (isItEC(name)) {
      det = 1;
      lay = 1;
    } else {
      if (isSensitive(name)) {
        if (isItHF(touch)) {
          det = 5;
          lay = 21;
          if (lay != layer)
            nlayHF++;
        } else {
          det = (touch->GetReplicaNumber(1)) / 1000;
          lay = (touch->GetReplicaNumber(0) / 10) % 100 + 3;
          if (det == 4) {
            double abeta = std::abs(eta);
            if (abeta < 1.479)
              lay = layer + 1;
            else
              lay--;
            if (lay < 3)
              lay = 3;
            if (lay == layer)
              lay++;
            if (lay > 20)
              lay = 20;
            if (lay != layer)
              nlayHE++;
          } else if (lay != layer) {
            if (lay >= 20)
              nlayHO++;
            else
              nlayHB++;
          }
        }
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Det " << det << " Layer " << lay << " Eta "
                                           << eta << " Phi " << phi / CLHEP::deg;
#endif
      } else if (layer == 1) {
        det = -1;
        lay = 2;
      }
    }
    if (det != 0) {
      if (lay != layer) {
        id = lay;
        layer = lay;
      }
    }
 
    if (id > idOld) {
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Step at " << name << " calls filHisto with "
                                         << (id - 1);
#endif
      fillHisto(id - 1);
    }
  }
 
  stepLen += step;
  radLen += (step / radl);
  intLen += (step / intl);
  if (fillHistos) {
    if (id == 21) {
      if (!isItHF(aStep->GetPostStepPoint()->GetTouchable())) {
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("MaterialBudget")
            << "MaterialBudgetHcalHistos: After HF in " << name << ":"
            << aStep->GetPostStepPoint()->GetTouchable()->GetVolume(0)->GetName() << " calls fillHisto with " << id;
#endif
        fillHisto(idOld);
        id++;
        layer = 0;
      }
    }
  }
}

References funct::abs(), fastSimProducer_cff::density, eta, fillHisto(), fillHistos, newFWLiteAna::found, g, id, cuy::ii, intLen, intLength, isItEC(), isItHF(), isSensitive(), layer, matList, Skims_PA_cff::name, nlayHB, nlayHE, nlayHF, nlayHO, phi, printSum, radLen, radLength, stepLen, stepLength, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MaterialBudgetHcal::update().

fillStartTrack()

void MaterialBudgetHcalHistos::fillStartTrack

(

const G4Track *

aTrack

)

Definition at line 81 of file MaterialBudgetHcalHistos.cc.

                                                                   {
  id = layer = steps = 0;
  radLen = intLen = stepLen = 0;
  nlayHB = nlayHE = nlayHF = nlayHO = 0;
 
  const G4ThreeVector& dir = aTrack->GetMomentum();
  if (dir.theta() != 0) {
    eta = dir.eta();
  } else {
    eta = -99;
  }
  phi = dir.phi();
  double theEnergy = aTrack->GetTotalEnergy();
  int theID = (int)(aTrack->GetDefinition()->GetPDGEncoding());
 
  if (printSum) {
    matList.clear();
    stepLength.clear();
    radLength.clear();
    intLength.clear();
  }
 
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Track " << aTrack->GetTrackID() << " Code " << theID
                                     << " Energy " << theEnergy / CLHEP::GeV << " GeV; Eta " << eta << " Phi "
                                     << phi / CLHEP::deg << " PT " << dir.perp() / CLHEP::GeV << " GeV *****";
}

References DeadROC_duringRun::dir, eta, GeV, createfilelist::int, intLen, intLength, layer, matList, nlayHB, nlayHE, nlayHF, nlayHO, phi, printSum, radLen, radLength, stepLen, stepLength, and steps.

Referenced by MaterialBudgetHcal::update().

getDDDArray()

std::vector< double > MaterialBudgetHcalHistos::getDDDArray ( const std::string &  str, const DDsvalues_type &  sv  )

private

Definition at line 424 of file MaterialBudgetHcalHistos.cc.

                                                                                                      {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos:getDDDArray called for " << str;
#endif
  DDValue value(str);
  if (DDfetch(&sv, value)) {
    LogDebug("MaterialBudget") << value;
    const std::vector<double>& fvec = value.doubles();
    int nval = fvec.size();
    if (nval < 1) {
      throw cms::Exception("MaterialBudgetHcalHistos") << "nval = " << nval << " < 1 for array " << str << "\n";
    }
 
    return fvec;
  } else {
    throw cms::Exception("MaterialBudgetHcalHistos") << "cannot get array " << str << "\n";
  }
}

References DDfetch(), Exception, LogDebug, str, pfDeepBoostedJetPreprocessParams_cfi::sv, and relativeConstraints::value.

Referenced by fillBeginJob().

getNames()

std::vector< std::string > MaterialBudgetHcalHistos::getNames ( DDFilteredView &  fv )

private

Definition at line 407 of file MaterialBudgetHcalHistos.cc.

                                                                          {
  std::vector<std::string> tmp;
  bool dodet = fv.firstChild();
  while (dodet) {
    const DDLogicalPart& log = fv.logicalPart();
    std::string namx = log.name().name();
    bool ok = true;
    for (unsigned int i = 0; i < tmp.size(); i++)
      if (namx == tmp[i])
        ok = false;
    if (ok)
      tmp.push_back(namx);
    dodet = fv.next();
  }
  return tmp;
}

References DDFilteredView::firstChild(), mps_fire::i, dqm-mbProfile::log, DDFilteredView::logicalPart(), DDFilteredView::next(), convertSQLiteXML::ok, AlCaHLTBitMon_QueryRunRegistry::string, and createJobs::tmp.

Referenced by fillBeginJob().

hend()

void MaterialBudgetHcalHistos::hend ( )

private

Definition at line 403 of file MaterialBudgetHcalHistos.cc.

                                    {
  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalHistos: Save user histos ===";
}

Referenced by ~MaterialBudgetHcalHistos().

isItEC()

bool MaterialBudgetHcalHistos::isItEC ( std::string  name )

private

Definition at line 465 of file MaterialBudgetHcalHistos.cc.

                                                    {
  std::vector<std::string>::const_iterator it = sensitiveEC.begin();
  std::vector<std::string>::const_iterator itEnd = sensitiveEC.end();
  for (; it != itEnd; ++it)
    if (name == *it)
      return true;
  return false;
}

References Skims_PA_cff::name, and sensitiveEC.

Referenced by fillPerStep().

isItHF()

bool MaterialBudgetHcalHistos::isItHF ( const G4VTouchable *  touch )

private

Definition at line 452 of file MaterialBudgetHcalHistos.cc.

                                                               {
  int levels = ((touch->GetHistoryDepth()) + 1);
  for (unsigned int it = 0; it < hfNames.size(); it++) {
    if (levels >= hfLevels[it]) {
      std::string name = touch->GetVolume(levels - hfLevels[it])->GetName();
      if (name == hfNames[it]) {
        return true;
      }
    }
  }
  return false;
}

References hfLevels, hfNames, jets_cff::levels, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by fillPerStep().

isSensitive()

bool MaterialBudgetHcalHistos::isSensitive ( std::string  name )

private

Definition at line 443 of file MaterialBudgetHcalHistos.cc.

                                                         {
  std::vector<std::string>::const_iterator it = sensitives.begin();
  std::vector<std::string>::const_iterator itEnd = sensitives.end();
  for (; it != itEnd; ++it)
    if (name == *it)
      return true;
  return false;
}

References Skims_PA_cff::name, and sensitives.

Referenced by fillPerStep().

Member Data Documentation

binEta

int MaterialBudgetHcalHistos::binEta

private

Definition at line 45 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and MaterialBudgetHcalHistos().

binPhi

int MaterialBudgetHcalHistos::binPhi

private

Definition at line 45 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and MaterialBudgetHcalHistos().

eta

double MaterialBudgetHcalHistos::eta

private

Definition at line 57 of file MaterialBudgetHcalHistos.h.

Referenced by Particle.Particle::__str__(), fillHisto(), fillLayer(), fillPerStep(), fillStartTrack(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

etaHigh

double MaterialBudgetHcalHistos::etaHigh

private

Definition at line 46 of file MaterialBudgetHcalHistos.h.

Referenced by fillHisto(), and MaterialBudgetHcalHistos().

etaLow

double MaterialBudgetHcalHistos::etaLow

private

Definition at line 46 of file MaterialBudgetHcalHistos.h.

Referenced by fillHisto(), and MaterialBudgetHcalHistos().

fillHistos

bool MaterialBudgetHcalHistos::fillHistos

private

Definition at line 44 of file MaterialBudgetHcalHistos.h.

Referenced by fillBeginJob(), fillEndTrack(), fillPerStep(), and MaterialBudgetHcalHistos().

hfLevels

std::vector<int> MaterialBudgetHcalHistos::hfLevels

private

Definition at line 43 of file MaterialBudgetHcalHistos.h.

Referenced by fillBeginJob(), and isItHF().

hfNames

std::vector<std::string> MaterialBudgetHcalHistos::hfNames

private

Definition at line 42 of file MaterialBudgetHcalHistos.h.

Referenced by fillBeginJob(), and isItHF().

id

int MaterialBudgetHcalHistos::id

private

Definition at line 55 of file MaterialBudgetHcalHistos.h.

Referenced by fillPerStep().

intLen

double MaterialBudgetHcalHistos::intLen

private

Definition at line 56 of file MaterialBudgetHcalHistos.h.

Referenced by fillHisto(), fillPerStep(), and fillStartTrack().

intLength

std::vector<double> MaterialBudgetHcalHistos::intLength

private

Definition at line 48 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillPerStep(), and fillStartTrack().

layer

int MaterialBudgetHcalHistos::layer

private

Definition at line 55 of file MaterialBudgetHcalHistos.h.

Referenced by fillPerStep(), fillStartTrack(), geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

matList

std::vector<std::string> MaterialBudgetHcalHistos::matList

private

Definition at line 47 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillPerStep(), and fillStartTrack().

maxEta

double MaterialBudgetHcalHistos::maxEta

private

Definition at line 46 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and MaterialBudgetHcalHistos().

maxSet

const int MaterialBudgetHcalHistos::maxSet = 25

staticprivate

Definition at line 41 of file MaterialBudgetHcalHistos.h.

Referenced by book(), fillEndTrack(), and fillHisto().

maxSet2

const int MaterialBudgetHcalHistos::maxSet2 = 9

staticprivate

Definition at line 41 of file MaterialBudgetHcalHistos.h.

Referenced by book().

me100

TProfile* MaterialBudgetHcalHistos::me100[maxSet]

private

Definition at line 51 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me1000

TProfile2D * MaterialBudgetHcalHistos::me1000[maxSet]

private

Definition at line 54 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me1100

TProfile2D * MaterialBudgetHcalHistos::me1100[maxSet]

private

Definition at line 54 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me1200

TH2F* MaterialBudgetHcalHistos::me1200[maxSet]

private

Definition at line 50 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me1300

TH1F * MaterialBudgetHcalHistos::me1300[maxSet2]

private

Definition at line 49 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillLayer().

me1400

TH2F * MaterialBudgetHcalHistos::me1400[maxSet2]

private

Definition at line 50 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillLayer().

me1500

TProfile* MaterialBudgetHcalHistos::me1500[maxSet2]

private

Definition at line 53 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillLayer().

me200

TProfile * MaterialBudgetHcalHistos::me200[maxSet]

private

Definition at line 51 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me300

TProfile * MaterialBudgetHcalHistos::me300[maxSet]

private

Definition at line 51 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me400

TH1F* MaterialBudgetHcalHistos::me400[maxSet]

private

Definition at line 49 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me500

TProfile* MaterialBudgetHcalHistos::me500[maxSet]

private

Definition at line 52 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me600

TProfile * MaterialBudgetHcalHistos::me600[maxSet]

private

Definition at line 52 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me700

TProfile * MaterialBudgetHcalHistos::me700[maxSet]

private

Definition at line 52 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me800

TH1F * MaterialBudgetHcalHistos::me800[maxSet]

private

Definition at line 49 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

me900

TProfile2D* MaterialBudgetHcalHistos::me900[maxSet]

private

Definition at line 54 of file MaterialBudgetHcalHistos.h.

Referenced by book(), and fillHisto().

nlayHB

int MaterialBudgetHcalHistos::nlayHB

private

Definition at line 58 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillLayer(), fillPerStep(), and fillStartTrack().

nlayHE

int MaterialBudgetHcalHistos::nlayHE

private

Definition at line 58 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillLayer(), fillPerStep(), and fillStartTrack().

nlayHF

int MaterialBudgetHcalHistos::nlayHF

private

Definition at line 58 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillLayer(), fillPerStep(), and fillStartTrack().

nlayHO

int MaterialBudgetHcalHistos::nlayHO

private

Definition at line 58 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillLayer(), fillPerStep(), and fillStartTrack().

phi

double MaterialBudgetHcalHistos::phi

private

Definition at line 57 of file MaterialBudgetHcalHistos.h.

Referenced by Particle.Particle::__str__(), fillHisto(), fillLayer(), fillPerStep(), fillStartTrack(), and ntupleDataFormat.Track::phiPull().

printSum

bool MaterialBudgetHcalHistos::printSum

private

Definition at line 44 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillPerStep(), fillStartTrack(), and MaterialBudgetHcalHistos().

radLen

double MaterialBudgetHcalHistos::radLen

private

Definition at line 56 of file MaterialBudgetHcalHistos.h.

Referenced by fillHisto(), fillPerStep(), and fillStartTrack().

radLength

std::vector<double> MaterialBudgetHcalHistos::radLength

private

Definition at line 48 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillPerStep(), and fillStartTrack().

sensitiveEC

std::vector<std::string> MaterialBudgetHcalHistos::sensitiveEC

private

Definition at line 42 of file MaterialBudgetHcalHistos.h.

Referenced by fillBeginJob(), and isItEC().

sensitives

std::vector<std::string> MaterialBudgetHcalHistos::sensitives

private

Definition at line 42 of file MaterialBudgetHcalHistos.h.

Referenced by fillBeginJob(), and isSensitive().

stepLen

double MaterialBudgetHcalHistos::stepLen

private

Definition at line 56 of file MaterialBudgetHcalHistos.h.

Referenced by fillHisto(), fillPerStep(), and fillStartTrack().

stepLength

std::vector<double> MaterialBudgetHcalHistos::stepLength

private

Definition at line 48 of file MaterialBudgetHcalHistos.h.

Referenced by fillEndTrack(), fillPerStep(), and fillStartTrack().

steps

int MaterialBudgetHcalHistos::steps

private

Definition at line 55 of file MaterialBudgetHcalHistos.h.

Referenced by fillStartTrack().