#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
static const int	maxSet2 = 9
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

Detailed Description

Definition at line 19 of file MaterialBudgetHcalHistos.h.

Constructor & Destructor Documentation

MaterialBudgetHcalHistos::MaterialBudgetHcalHistos ( const edm::ParameterSet & p )

Definition at line 15 of file MaterialBudgetHcalHistos.cc.

References binEta, binPhi, book(), etaHigh, etaLow, fillHistos, edm::ParameterSet::getUntrackedParameter(), maxEta, pi, and printSum.

                                                                          {

  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::LogInfo("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();

}

virtual MaterialBudgetHcalHistos::~MaterialBudgetHcalHistos ( ) [inline, virtual]

Definition at line 24 of file MaterialBudgetHcalHistos.h.

References hend().

{ hend(); }

Member Function Documentation

void MaterialBudgetHcalHistos::book ( ) [private]

Definition at line 260 of file MaterialBudgetHcalHistos.cc.

References binEta, binPhi, i, edm::Service< T >::isAvailable(), maxEta, maxSet, maxSet2, me100, me1000, me1100, me1200, me1300, me1400, me1500, me200, me300, me400, me500, me600, me700, me800, me900, mergeVDriftHistosByStation::name, pi, and indexGen::title.

Referenced by MaterialBudgetHcalHistos().

                                    {

  // 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::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: Booking user "
                                 << "histos === with " << binEta << " bins "
                                 << "in eta from " << -maxEta << " to "
                                 << maxEta << " and " << binPhi << " bins "
                                 << "in phi from " << -maxPhi << " to " 
                                 << maxPhi;
  
  char  name[10], title[40];
  // total X0
  for (int i=0; i<maxSet; i++) {
    sprintf(name, "%d", i+100);
    sprintf(title, "MB(X0) prof Eta in region %d", i);
    me100[i] =  tfile->make<TProfile>(name, title, binEta, -maxEta, maxEta);
    sprintf(name, "%d", i+200);
    sprintf(title, "MB(L0) prof Eta in region %d", i);
    me200[i] = tfile->make<TProfile>(name, title, binEta, -maxEta, maxEta);
    sprintf(name, "%d", i+300);
    sprintf(title, "MB(Step) prof Eta in region %d", i);
    me300[i] = tfile->make<TProfile>(name, title, binEta, -maxEta, maxEta);
    sprintf(name, "%d", i+400);
    sprintf(title, "Eta in region %d", i);
    me400[i] = tfile->make<TH1F>(name, title, binEta, -maxEta, maxEta);
    sprintf(name, "%d", i+500);
    sprintf(title, "MB(X0) prof Ph in region %d", i);
    me500[i] = tfile->make<TProfile>(name, title, binPhi, -maxPhi, maxPhi);
    sprintf(name, "%d", i+600);
    sprintf(title, "MB(L0) prof Ph in region %d", i);
    me600[i] = tfile->make<TProfile>(name, title, binPhi, -maxPhi, maxPhi);
    sprintf(name, "%d", i+700);
    sprintf(title, "MB(Step) prof Ph in region %d", i);
    me700[i] = tfile->make<TProfile>(name, title, binPhi, -maxPhi, maxPhi);
    sprintf(name, "%d", i+800);
    sprintf(title, "Phi in region %d", i);
    me800[i] = tfile->make<TH1F>(name, title, binPhi, -maxPhi, maxPhi);
    sprintf(name, "%d", i+900);
    sprintf(title, "MB(X0) prof Eta Phi in region %d", i);
    me900[i] = tfile->make<TProfile2D>(name, title, binEta/2, -maxEta, maxEta,
                                       binPhi/2, -maxPhi, maxPhi);
    sprintf(name, "%d", i+1000);
    sprintf(title, "MB(L0) prof Eta Phi in region %d", i);
    me1000[i]= tfile->make<TProfile2D>(name, title, binEta/2, -maxEta, maxEta,
                                       binPhi/2, -maxPhi, maxPhi);
    sprintf(name, "%d", i+1100);
    sprintf(title, "MB(Step) prof Eta Phi in region %d", i);
    me1100[i]= tfile->make<TProfile2D>(name, title, binEta/2, -maxEta, maxEta,
                                       binPhi/2, -maxPhi, maxPhi);
    sprintf(name, "%d", i+1200);
    sprintf(title, "Eta vs Phi in region %d", i);
    me1200[i]= tfile->make<TH2F>(name, title, binEta/2, -maxEta, maxEta, 
                                 binPhi/2, -maxPhi, maxPhi);
  }
  for (int i=0; i<maxSet2; i++) {
    sprintf(name, "%d", i+1300);
    sprintf(title, "Events with layers Hit (0 all, 1 HB, ..) for %d", i);
    me1300[i]= tfile->make<TH1F>(name, title, binEta, -maxEta, maxEta);
    sprintf(name, "%d", i+1400);
    sprintf(title, "Eta vs Phi for layers hit in %d", i);
    me1400[i]= tfile->make<TH2F>(name, title, binEta/2, -maxEta, maxEta, 
                                 binPhi/2, -maxPhi, maxPhi);
    sprintf(name, "%d", i+1500);
    sprintf(title, "Number of layers crossed (0 all, 1 HB, ..) for %d", i);
    me1500[i]= tfile->make<TProfile>(name, title, binEta, -maxEta, maxEta);
  }

  edm::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: Booking user "
                                 << "histos done ===";

}

void MaterialBudgetHcalHistos::fillBeginJob ( const DDCompactView & cpv )

Definition at line 35 of file MaterialBudgetHcalHistos.cc.

References DDFilteredView::addFilter(), DDSpecificsFilter::equals, fillHistos, DDFilteredView::firstChild(), getDDDArray(), getNames(), hfLevels, hfNames, i, gen::k, testEve_cfg::level, DDFilteredView::mergedSpecifics(), sensitiveEC, sensitives, DDSpecificsFilter::setCriteria(), AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, and relativeConstraints::value.

Referenced by MaterialBudgetHcal::update().

                                                                     {

  if (fillHistos) {
    std::string attribute = "ReadOutName";
    std::string value     = "HcalHits";
    DDSpecificsFilter filter1;
    DDValue           ddv1(attribute,value,0);
    filter1.setCriteria(ddv1,DDSpecificsFilter::equals);
    DDFilteredView fv1(cpv);
    fv1.addFilter(filter1);
    sensitives = getNames(fv1);
    edm::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: Names to be "
                                   << "tested for " << attribute << " = " 
                                   << value << " has " << sensitives.size()
                                   << " elements";
    for (unsigned int i=0; i<sensitives.size(); i++) 
      edm::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: sensitives["
                                     << i << "] = " << sensitives[i];

    attribute = "Volume";
    value     = "HF";
    DDSpecificsFilter filter2;
    DDValue           ddv2(attribute,value,0);
    filter2.setCriteria(ddv2,DDSpecificsFilter::equals);
    DDFilteredView fv2(cpv);
    fv2.addFilter(filter2);
    hfNames = getNames(fv2);
    fv2.firstChild();
    DDsvalues_type sv(fv2.mergedSpecifics());
    std::vector<double> temp =  getDDDArray("Levels",sv);
    edm::LogInfo("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::LogInfo("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];
      DDSpecificsFilter filter3;
      DDValue           ddv3(attribute,value,0);
      filter3.setCriteria(ddv3,DDSpecificsFilter::equals);
      DDFilteredView fv3(cpv);
      fv3.addFilter(filter3);
      std::vector<std::string> senstmp = getNames(fv3);
      edm::LogInfo("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::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos:sensitiveEC["
                                     << i << "] = " << sensitiveEC[i];
  }
}

void MaterialBudgetHcalHistos::fillEndTrack ( )

Definition at line 243 of file MaterialBudgetHcalHistos.cc.

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

Referenced by MaterialBudgetHcal::update().

                                            {
  edm::LogInfo("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::LogInfo("MaterialBudget") << matList[ii] << "\t" << stepLength[ii]
                                     << "\t" << radLength[ii] << "\t"
                                     << intLength[ii];
    }
  }
}

void MaterialBudgetHcalHistos::fillHisto ( int ii ) [private]

Definition at line 339 of file MaterialBudgetHcalHistos.cc.

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

Referenced by fillEndTrack(), and fillPerStep().

                                               {

  LogDebug("MaterialBudget") << "MaterialBudgetHcalHistos:FillHisto called "
                             << "with index " << ii << " integrated  step "
                             << stepLen << " X0 " << radLen << " Lamda " 
                             << intLen;
  
  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);
    
  }
}

void MaterialBudgetHcalHistos::fillLayer ( ) [private]

Definition at line 367 of file MaterialBudgetHcalHistos.cc.

References abs, eta, me1300, me1400, me1500, nlayHB, nlayHE, nlayHF, nlayHO, and phi.

Referenced by fillEndTrack().

                                         {

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

void MaterialBudgetHcalHistos::fillPerStep ( const G4Step * aStep )

Definition at line 131 of file MaterialBudgetHcalHistos.cc.

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

Referenced by MaterialBudgetHcal::update().

                                                              {

  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();
  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::LogInfo("MaterialBudget") << name << " " << step << " " << matName 
                                   << " " << stepLen << " " << step/radl << " " 
                                   << radLen << " " <<step/intl << " " <<intLen;
  } else {
    edm::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: Step at " 
                                   << name << " 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++;
          }
        }
        LogDebug("MaterialBudget") << "MaterialBudgetHcalHistos: Det " << det
                                   << " Layer " << lay << " Eta " << eta 
                                   << " Phi " << phi/deg;
      } else if (layer == 1) {
        det = -1;
        lay = 2;
      }
    }
    if (det != 0) {
      if (lay != layer) {
        id    = lay;
        layer = lay;
      }
    }

    if (id > idOld) {
      //    edm::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: Step at " << name;
      fillHisto(id-1);
    }
  }

  stepLen += step;
  radLen  += step/radl;
  intLen  += step/intl;
  if (fillHistos) {
    if (layer == 21 && det == 5) {
      if (!isItHF(aStep->GetPostStepPoint()->GetTouchable())) {
        LogDebug("MaterialBudget") << "MaterialBudgetHcalHistos: After HF in " 
                                   << aStep->GetPostStepPoint()->GetTouchable()->GetVolume(0)->GetName();
        fillHisto(id);
        id++;
        layer = 0;
      }
    }
  }
}

void MaterialBudgetHcalHistos::fillStartTrack ( const G4Track * aTrack )

Definition at line 100 of file MaterialBudgetHcalHistos.cc.

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

Referenced by MaterialBudgetHcal::update().

                                                                   {

  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::LogInfo("MaterialBudget") << "MaterialBudgetHcalHistos: Track " 
                                 << aTrack->GetTrackID() << " Code " << theID
                                 << " Energy " << theEnergy/GeV << " GeV; Eta "
                                 << eta << " Phi " << phi/deg << " PT "
                                 << dir.perp()/GeV << " GeV *****";
}

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

Definition at line 430 of file MaterialBudgetHcalHistos.cc.

References DDfetch(), DDValue::doubles(), Exception, LogDebug, and relativeConstraints::value.

Referenced by fillBeginJob().

                                                                                     {

  LogDebug("MaterialBudget") << "MaterialBudgetHcalHistos:getDDDArray called "
                             << "for " << str;
  DDValue value(str);
  if (DDfetch(&sv,value)) {
    LogDebug("MaterialBudget") << value;
    const std::vector<double> & fvec = value.doubles();
    int nval = fvec.size();
    if (nval < 1) {
      edm::LogError("MaterialBudget") << "MaterialBudgetHcalHistos : # of " 
                                      << str << " bins " << nval
                                      << " < 1 ==> illegal";
      throw cms::Exception("Unknown", "MaterialBudgetHcalHistos")
        << "nval < 1 for array " << str << "\n";
    }

    return fvec;
  } else {
    edm::LogError("MaterialBudget") << "MaterialBudgetHcalHistos : cannot get "
                                    << "array " << str;
    throw cms::Exception("Unknown", "MaterialBudgetHcalHistos")
      << "cannot get array " << str << "\n";
  }
}

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

Definition at line 414 of file MaterialBudgetHcalHistos.cc.

References DDFilteredView::firstChild(), i, create_public_lumi_plots::log, DDFilteredView::logicalPart(), DDName::name(), DDBase< N, C >::name(), DDFilteredView::next(), convertSQLiteXML::ok, AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

Referenced by fillBeginJob().

                                                                          {

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

void MaterialBudgetHcalHistos::hend ( ) [private]

Definition at line 409 of file MaterialBudgetHcalHistos.cc.

Referenced by ~MaterialBudgetHcalHistos().

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

bool MaterialBudgetHcalHistos::isItEC ( std::string name ) [private]

Definition at line 479 of file MaterialBudgetHcalHistos.cc.

References sensitiveEC.

Referenced by fillPerStep().

                                                     {

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

bool MaterialBudgetHcalHistos::isItHF ( const G4VTouchable * touch ) [private]

Definition at line 466 of file MaterialBudgetHcalHistos.cc.

References hfLevels, hfNames, mergeVDriftHistosByStation::name, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by fillPerStep().

                                                                {

  // std::vector<std::string>::const_iterator it = hfNames.begin();
  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;
}

bool MaterialBudgetHcalHistos::isSensitive ( std::string name ) [private]

Definition at line 457 of file MaterialBudgetHcalHistos.cc.

References sensitives.

Referenced by fillPerStep().

                                                          {

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