Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-10-25 10:06:22

0001 #include "Validation/Geometry/interface/MaterialBudgetMtdHistos.h"
0002 #include "Validation/Geometry/interface/MaterialBudgetData.h"
0003 
0004 #include "CLHEP/Units/GlobalSystemOfUnits.h"
0005 
0006 MaterialBudgetMtdHistos::MaterialBudgetMtdHistos(std::shared_ptr<MaterialBudgetData> data,
0007                                                  std::shared_ptr<TestHistoMgr> mgr,
0008                                                  const std::string& fileName)
0009     : MaterialBudgetFormat(data), hmgr(mgr) {
0010   theFileName = fileName;
0011   book();
0012 }
0013 
0014 void MaterialBudgetMtdHistos::book() {
0015   edm::LogInfo("MaterialBudget") << "MaterialBudgetMtdHistos: Booking user histos";
0016 
0017   static constexpr double minEta = -5.;
0018   static constexpr double maxEta = 5.;
0019   static constexpr double minPhi = -3.1416;
0020   static constexpr double maxPhi = 3.1416;
0021   static constexpr int nbinEta = 250;
0022   static constexpr int nbinPhi = 180;
0023 
0024   static constexpr double minEtaBTLZoom = 0.;
0025   static constexpr double maxEtaBTLZoom = 0.087;
0026   static constexpr double minPhiBTLZoom = 0.;
0027   static constexpr double maxPhiBTLZoom = 0.35;
0028   static constexpr int nbinEtaBTLZoom = 64;
0029   static constexpr int nbinPhiBTLZoom = 20;
0030 
0031   static constexpr double minMB = 0.;
0032   static constexpr double maxMBetl = 0.025;
0033   static constexpr double maxMBbtl = 0.5;
0034   static constexpr int nbinMB = 25;
0035 
0036   // Material budget: radiation length
0037   // total X0
0038   hmgr->addHistoProf1(new TProfile("10", "MB prof Eta [Total];#eta;x/X_{0} ", nbinEta, minEta, maxEta));
0039   hmgr->addHisto1(new TH1F("11", "Eta ", nbinEta, minEta, maxEta));
0040   hmgr->addHistoProf1(new TProfile("20", "MB prof Phi [Total];#varphi [rad];x/X_{0} ", nbinPhi, minPhi, maxPhi));
0041   hmgr->addHisto1(new TH1F("21", "Phi ", nbinPhi, minPhi, maxPhi));
0042   hmgr->addHistoProf2(new TProfile2D(
0043       "30", "MB prof Eta  Phi [Total];#eta;#varphi;x/X_{0} ", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0044   hmgr->addHisto2(new TH2F("31", "Eta vs Phi ", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0045 
0046   // Support
0047   hmgr->addHistoProf1(new TProfile("110", "MB prof Eta [Support];#eta;x/X_{0}", nbinEta, minEta, maxEta));
0048   hmgr->addHisto1(new TH1F("111", "Eta [Support]", nbinEta, minEta, maxEta));
0049   hmgr->addHistoProf1(new TProfile("120", "MB prof Phi [Support];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
0050   hmgr->addHisto1(new TH1F("121", "Phi [Support]", nbinPhi, minPhi, maxPhi));
0051   hmgr->addHistoProf2(new TProfile2D(
0052       "130", "MB prof Eta  Phi [Support];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0053   hmgr->addHisto2(new TH2F("131", "Eta vs Phi [Support]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0054 
0055   // Sensitive
0056   hmgr->addHistoProf1(new TProfile("210", "MB prof Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta));
0057   hmgr->addHisto1(new TH1F("211", "Eta [Sensitive]", nbinEta, minEta, maxEta));
0058   hmgr->addHistoProf1(new TProfile("220", "MB prof Phi [Sensitive];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
0059   hmgr->addHisto1(new TH1F("221", "Phi [Sensitive]", nbinPhi, minPhi, maxPhi));
0060   hmgr->addHistoProf2(new TProfile2D(
0061       "230", "MB prof Eta  Phi [Sensitive];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0062 
0063   hmgr->addHistoProf2(new TProfile2D("10230",
0064                                      "MB prof Eta  Phi [Sensitive];#eta;#varphi;x/X_{0}",
0065                                      nbinEtaBTLZoom,
0066                                      minEtaBTLZoom,
0067                                      maxEtaBTLZoom,
0068                                      nbinPhiBTLZoom,
0069                                      minPhiBTLZoom,
0070                                      maxPhiBTLZoom));
0071   hmgr->addHisto2(
0072       new TH2F("10234", "MB vs Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta, nbinMB, minMB, maxMBetl));
0073   hmgr->addHisto2(new TH2F(
0074       "20234", "MB along z vs Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta, nbinMB, minMB, maxMBetl));
0075   hmgr->addHisto2(
0076       new TH2F("10235", "MB vs Eta [Sensitive];#eta;x/X_{0}", nbinEta, minEta, maxEta, nbinMB, minMB, maxMBbtl));
0077 
0078   hmgr->addHisto2(new TH2F("231", "Eta vs Phi [Sensitive]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0079 
0080   // Cables
0081   hmgr->addHistoProf1(new TProfile("310", "MB prof Eta [Cables];#eta;x/X_{0}", nbinEta, minEta, maxEta));
0082   hmgr->addHisto1(new TH1F("311", "Eta [Cables]", nbinEta, minEta, maxEta));
0083   hmgr->addHistoProf1(new TProfile("320", "MB prof Phi [Cables];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
0084   hmgr->addHisto1(new TH1F("321", "Phi [Cables]", nbinPhi, minPhi, maxPhi));
0085   hmgr->addHistoProf2(new TProfile2D(
0086       "330", "MB prof Eta  Phi [Cables];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0087   hmgr->addHisto2(new TH2F("331", "Eta vs Phi [Cables]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0088 
0089   // Cooling
0090   hmgr->addHistoProf1(new TProfile("410", "MB prof Eta [Cooling];#eta;x/X_{0}", nbinEta, minEta, maxEta));
0091   hmgr->addHisto1(new TH1F("411", "Eta [Cooling]", nbinEta, minEta, maxEta));
0092   hmgr->addHistoProf1(new TProfile("420", "MB prof Phi [Cooling];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
0093   hmgr->addHisto1(new TH1F("421", "Phi [Cooling]", nbinPhi, minPhi, maxPhi));
0094   hmgr->addHistoProf2(new TProfile2D(
0095       "430", "MB prof Eta  Phi [Cooling];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0096   hmgr->addHisto2(new TH2F("431", "Eta vs Phi [Cooling]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0097 
0098   // Electronics
0099   hmgr->addHistoProf1(new TProfile("510", "MB prof Eta [Electronics];#eta;x/X_{0}", nbinEta, minEta, maxEta));
0100   hmgr->addHisto1(new TH1F("511", "Eta [Electronics]", nbinEta, minEta, maxEta));
0101   hmgr->addHistoProf1(new TProfile("520", "MB prof Phi [Electronics];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
0102   hmgr->addHisto1(new TH1F("521", "Phi [Electronics]", nbinPhi, minPhi, maxPhi));
0103   hmgr->addHistoProf2(new TProfile2D(
0104       "530", "MB prof Eta  Phi [Electronics];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0105   hmgr->addHisto2(new TH2F("531", "Eta vs Phi [Electronics]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0106 
0107   // Other
0108   hmgr->addHistoProf1(new TProfile("610", "MB prof Eta [Other];#eta;x/X_{0}", nbinEta, minEta, maxEta));
0109   hmgr->addHisto1(new TH1F("611", "Eta [Other]", nbinEta, minEta, maxEta));
0110   hmgr->addHistoProf1(new TProfile("620", "MB prof Phi [Other];#varphi [rad];x/X_{0}", nbinPhi, minPhi, maxPhi));
0111   hmgr->addHisto1(new TH1F("621", "Phi [Other]", nbinPhi, minPhi, maxPhi));
0112   hmgr->addHistoProf2(new TProfile2D(
0113       "630", "MB prof Eta  Phi [Other];#eta;#varphi;x/X_{0}", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0114   hmgr->addHisto2(new TH2F("631", "Eta vs Phi [Other]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0115 
0116   // Material budget: interaction length
0117   // total Lambda0
0118   hmgr->addHistoProf1(new TProfile("1010", "MB prof Eta [Total];#eta;#lambda/#lambda_{0} ", nbinEta, minEta, maxEta));
0119   hmgr->addHisto1(new TH1F("1011", "Eta ", nbinEta, minEta, maxEta));
0120   hmgr->addHistoProf1(
0121       new TProfile("1020", "MB prof Phi [Total];#varphi [rad];#lambda/#lambda_{0} ", nbinPhi, minPhi, maxPhi));
0122   hmgr->addHisto1(new TH1F("1021", "Phi ", nbinPhi, minPhi, maxPhi));
0123   hmgr->addHistoProf2(new TProfile2D("1030",
0124                                      "MB prof Eta  Phi [Total];#eta;#varphi;#lambda/#lambda_{0} ",
0125                                      nbinEta,
0126                                      minEta,
0127                                      maxEta,
0128                                      nbinPhi,
0129                                      minPhi,
0130                                      maxPhi));
0131   hmgr->addHisto2(new TH2F("1031", "Eta vs Phi ", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0132 
0133   // Support
0134   hmgr->addHistoProf1(new TProfile("1110", "MB prof Eta [Support];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
0135   hmgr->addHisto1(new TH1F("1111", "Eta [Support]", nbinEta, minEta, maxEta));
0136   hmgr->addHistoProf1(
0137       new TProfile("1120", "MB prof Phi [Support];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
0138   hmgr->addHisto1(new TH1F("1121", "Phi [Support]", nbinPhi, minPhi, maxPhi));
0139   hmgr->addHistoProf2(new TProfile2D("1130",
0140                                      "MB prof Eta  Phi [Support];#eta;#varphi;#lambda/#lambda_{0}",
0141                                      nbinEta,
0142                                      minEta,
0143                                      maxEta,
0144                                      nbinPhi,
0145                                      minPhi,
0146                                      maxPhi));
0147   hmgr->addHisto2(new TH2F("1131", "Eta vs Phi [Support]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0148 
0149   // Sensitive
0150   hmgr->addHistoProf1(
0151       new TProfile("1210", "MB prof Eta [Sensitive];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
0152   hmgr->addHisto1(new TH1F("1211", "Eta [Sensitive]", nbinEta, minEta, maxEta));
0153   hmgr->addHistoProf1(
0154       new TProfile("1220", "MB prof Phi [Sensitive];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
0155   hmgr->addHisto1(new TH1F("1221", "Phi [Sensitive]", nbinPhi, minPhi, maxPhi));
0156   hmgr->addHistoProf2(new TProfile2D("1230",
0157                                      "MB prof Eta  Phi [Sensitive];#eta;#varphi;#lambda/#lambda_{0}",
0158                                      nbinEta,
0159                                      minEta,
0160                                      maxEta,
0161                                      nbinPhi,
0162                                      minPhi,
0163                                      maxPhi));
0164   hmgr->addHisto2(new TH2F("1231", "Eta vs Phi [Sensitive]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0165 
0166   // Cables
0167   hmgr->addHistoProf1(new TProfile("1310", "MB prof Eta [Cables];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
0168   hmgr->addHisto1(new TH1F("1311", "Eta [Cables]", nbinEta, minEta, maxEta));
0169   hmgr->addHistoProf1(
0170       new TProfile("1320", "MB prof Phi [Cables];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
0171   hmgr->addHisto1(new TH1F("1321", "Phi [Cables]", nbinPhi, minPhi, maxPhi));
0172   hmgr->addHistoProf2(new TProfile2D("1330",
0173                                      "MB prof Eta  Phi [Cables];#eta;#varphi;#lambda/#lambda_{0}",
0174                                      nbinEta,
0175                                      minEta,
0176                                      maxEta,
0177                                      nbinPhi,
0178                                      minPhi,
0179                                      maxPhi));
0180   hmgr->addHisto2(new TH2F("1331", "Eta vs Phi [Cables]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0181 
0182   // Cooling
0183   hmgr->addHistoProf1(new TProfile("1410", "MB prof Eta [Cooling];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
0184   hmgr->addHisto1(new TH1F("1411", "Eta [Cooling]", nbinEta, minEta, maxEta));
0185   hmgr->addHistoProf1(
0186       new TProfile("1420", "MB prof Phi [Cooling];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
0187   hmgr->addHisto1(new TH1F("1421", "Phi [Cooling]", nbinPhi, minPhi, maxPhi));
0188   hmgr->addHistoProf2(new TProfile2D("1430",
0189                                      "MB prof Eta  Phi [Cooling];#eta;#varphi;#lambda/#lambda_{0}",
0190                                      nbinEta,
0191                                      minEta,
0192                                      maxEta,
0193                                      nbinPhi,
0194                                      minPhi,
0195                                      maxPhi));
0196   hmgr->addHisto2(new TH2F("1431", "Eta vs Phi [Cooling]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0197 
0198   // Electronics
0199   hmgr->addHistoProf1(
0200       new TProfile("1510", "MB prof Eta [Electronics];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
0201   hmgr->addHisto1(new TH1F("1511", "Eta [Electronics]", nbinEta, minEta, maxEta));
0202   hmgr->addHistoProf1(
0203       new TProfile("1520", "MB prof Phi [Electronics];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
0204   hmgr->addHisto1(new TH1F("1521", "Phi [Electronics]", nbinPhi, minPhi, maxPhi));
0205   hmgr->addHistoProf2(new TProfile2D("1530",
0206                                      "MB prof Eta  Phi [Electronics];#eta;#varphi;#lambda/#lambda_{0}",
0207                                      nbinEta,
0208                                      minEta,
0209                                      maxEta,
0210                                      nbinPhi,
0211                                      minPhi,
0212                                      maxPhi));
0213   hmgr->addHisto2(new TH2F("1531", "Eta vs Phi [Electronics]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0214 
0215   // Other
0216   hmgr->addHistoProf1(new TProfile("1610", "MB prof Eta [Other];#eta;#lambda/#lambda_{0}", nbinEta, minEta, maxEta));
0217   hmgr->addHisto1(new TH1F("1611", "Eta [Other]", nbinEta, minEta, maxEta));
0218   hmgr->addHistoProf1(
0219       new TProfile("1620", "MB prof Phi [Other];#varphi [rad];#lambda/#lambda_{0}", nbinPhi, minPhi, maxPhi));
0220   hmgr->addHisto1(new TH1F("1621", "Phi [Other]", nbinPhi, minPhi, maxPhi));
0221   hmgr->addHistoProf2(new TProfile2D("1630",
0222                                      "MB prof Eta  Phi [Other];#eta;#varphi;#lambda/#lambda_{0}",
0223                                      nbinEta,
0224                                      minEta,
0225                                      maxEta,
0226                                      nbinPhi,
0227                                      minPhi,
0228                                      maxPhi));
0229   hmgr->addHisto2(new TH2F("1631", "Eta vs Phi [Other]", nbinEta, minEta, maxEta, nbinPhi, minPhi, maxPhi));
0230 
0231   edm::LogInfo("MaterialBudget") << "MaterialBudgetMtdHistos: booking user histos done";
0232 }
0233 
0234 void MaterialBudgetMtdHistos::fillStartTrack() {}
0235 
0236 void MaterialBudgetMtdHistos::fillPerStep() {}
0237 
0238 void MaterialBudgetMtdHistos::fillEndTrack() {
0239   // Total X0
0240   hmgr->getHisto1(11)->Fill(theData->getEta());
0241   hmgr->getHisto1(21)->Fill(theData->getPhi());
0242   hmgr->getHisto2(31)->Fill(theData->getEta(), theData->getPhi());
0243 
0244   hmgr->getHistoProf1(10)->Fill(theData->getEta(), theData->getTotalMB());
0245   hmgr->getHistoProf1(20)->Fill(theData->getPhi(), theData->getTotalMB());
0246   hmgr->getHistoProf2(30)->Fill(theData->getEta(), theData->getPhi(), theData->getTotalMB());
0247 
0248   // Support
0249   hmgr->getHisto1(111)->Fill(theData->getEta());
0250   hmgr->getHisto1(121)->Fill(theData->getPhi());
0251   hmgr->getHisto2(131)->Fill(theData->getEta(), theData->getPhi());
0252 
0253   hmgr->getHistoProf1(110)->Fill(theData->getEta(), theData->getSupportMB());
0254   hmgr->getHistoProf1(120)->Fill(theData->getPhi(), theData->getSupportMB());
0255   hmgr->getHistoProf2(130)->Fill(theData->getEta(), theData->getPhi(), theData->getSupportMB());
0256 
0257   // Sensitive
0258   hmgr->getHisto1(211)->Fill(theData->getEta());
0259   hmgr->getHisto1(221)->Fill(theData->getPhi());
0260   hmgr->getHisto2(231)->Fill(theData->getEta(), theData->getPhi());
0261 
0262   hmgr->getHistoProf1(210)->Fill(theData->getEta(), theData->getSensitiveMB());
0263   hmgr->getHistoProf1(220)->Fill(theData->getPhi(), theData->getSensitiveMB());
0264   hmgr->getHistoProf2(230)->Fill(theData->getEta(), theData->getPhi(), theData->getSensitiveMB());
0265 
0266   static constexpr double bfTransitionEta = 1.55;
0267 
0268   if (std::abs(theData->getEta()) > bfTransitionEta) {
0269     hmgr->getHisto2(10234)->Fill(theData->getEta(), theData->getSensitiveMB());
0270     double norma = std::tanh(std::abs(theData->getEta()));
0271     hmgr->getHisto2(20234)->Fill(theData->getEta(), theData->getSensitiveMB() * norma);
0272   } else {
0273     hmgr->getHistoProf2(10230)->Fill(theData->getEta(), theData->getPhi(), theData->getSensitiveMB());
0274     hmgr->getHisto2(10235)->Fill(theData->getEta(), theData->getSensitiveMB());
0275   }
0276 
0277   // Cables
0278   hmgr->getHisto1(311)->Fill(theData->getEta());
0279   hmgr->getHisto1(321)->Fill(theData->getPhi());
0280   hmgr->getHisto2(331)->Fill(theData->getEta(), theData->getPhi());
0281 
0282   hmgr->getHistoProf1(310)->Fill(theData->getEta(), theData->getCablesMB());
0283   hmgr->getHistoProf1(320)->Fill(theData->getPhi(), theData->getCablesMB());
0284   hmgr->getHistoProf2(330)->Fill(theData->getEta(), theData->getPhi(), theData->getCablesMB());
0285 
0286   // Cooling
0287   hmgr->getHisto1(411)->Fill(theData->getEta());
0288   hmgr->getHisto1(421)->Fill(theData->getPhi());
0289   hmgr->getHisto2(431)->Fill(theData->getEta(), theData->getPhi());
0290 
0291   hmgr->getHistoProf1(410)->Fill(theData->getEta(), theData->getCoolingMB());
0292   hmgr->getHistoProf1(420)->Fill(theData->getPhi(), theData->getCoolingMB());
0293   hmgr->getHistoProf2(430)->Fill(theData->getEta(), theData->getPhi(), theData->getCoolingMB());
0294 
0295   // Electronics
0296   hmgr->getHisto1(511)->Fill(theData->getEta());
0297   hmgr->getHisto1(521)->Fill(theData->getPhi());
0298   hmgr->getHisto2(531)->Fill(theData->getEta(), theData->getPhi());
0299 
0300   hmgr->getHistoProf1(510)->Fill(theData->getEta(), theData->getElectronicsMB());
0301   hmgr->getHistoProf1(520)->Fill(theData->getPhi(), theData->getElectronicsMB());
0302   hmgr->getHistoProf2(530)->Fill(theData->getEta(), theData->getPhi(), theData->getElectronicsMB());
0303 
0304   // Other
0305   hmgr->getHisto1(611)->Fill(theData->getEta());
0306   hmgr->getHisto1(621)->Fill(theData->getPhi());
0307   hmgr->getHisto2(631)->Fill(theData->getEta(), theData->getPhi());
0308 
0309   hmgr->getHistoProf1(610)->Fill(theData->getEta(), theData->getOtherMB());
0310   hmgr->getHistoProf1(620)->Fill(theData->getPhi(), theData->getOtherMB());
0311   hmgr->getHistoProf2(630)->Fill(theData->getEta(), theData->getPhi(), theData->getOtherMB());
0312 
0313   // Total Lambda0
0314   hmgr->getHisto1(1011)->Fill(theData->getEta());
0315   hmgr->getHisto1(1021)->Fill(theData->getPhi());
0316   hmgr->getHisto2(1031)->Fill(theData->getEta(), theData->getPhi());
0317 
0318   hmgr->getHistoProf1(1010)->Fill(theData->getEta(), theData->getTotalIL());
0319   hmgr->getHistoProf1(1020)->Fill(theData->getPhi(), theData->getTotalIL());
0320   hmgr->getHistoProf2(1030)->Fill(theData->getEta(), theData->getPhi(), theData->getTotalIL());
0321 
0322   // Support
0323   hmgr->getHisto1(1111)->Fill(theData->getEta());
0324   hmgr->getHisto1(1121)->Fill(theData->getPhi());
0325   hmgr->getHisto2(1131)->Fill(theData->getEta(), theData->getPhi());
0326 
0327   hmgr->getHistoProf1(1110)->Fill(theData->getEta(), theData->getSupportIL());
0328   hmgr->getHistoProf1(1120)->Fill(theData->getPhi(), theData->getSupportIL());
0329   hmgr->getHistoProf2(1130)->Fill(theData->getEta(), theData->getPhi(), theData->getSupportIL());
0330 
0331   // Sensitive
0332   hmgr->getHisto1(1211)->Fill(theData->getEta());
0333   hmgr->getHisto1(1221)->Fill(theData->getPhi());
0334   hmgr->getHisto2(1231)->Fill(theData->getEta(), theData->getPhi());
0335 
0336   hmgr->getHistoProf1(1210)->Fill(theData->getEta(), theData->getSensitiveIL());
0337   hmgr->getHistoProf1(1220)->Fill(theData->getPhi(), theData->getSensitiveIL());
0338   hmgr->getHistoProf2(1230)->Fill(theData->getEta(), theData->getPhi(), theData->getSensitiveIL());
0339 
0340   // Cables
0341   hmgr->getHisto1(1311)->Fill(theData->getEta());
0342   hmgr->getHisto1(1321)->Fill(theData->getPhi());
0343   hmgr->getHisto2(1331)->Fill(theData->getEta(), theData->getPhi());
0344 
0345   hmgr->getHistoProf1(1310)->Fill(theData->getEta(), theData->getCablesIL());
0346   hmgr->getHistoProf1(1320)->Fill(theData->getPhi(), theData->getCablesIL());
0347   hmgr->getHistoProf2(1330)->Fill(theData->getEta(), theData->getPhi(), theData->getCablesIL());
0348 
0349   // Cooling
0350   hmgr->getHisto1(1411)->Fill(theData->getEta());
0351   hmgr->getHisto1(1421)->Fill(theData->getPhi());
0352   hmgr->getHisto2(1431)->Fill(theData->getEta(), theData->getPhi());
0353 
0354   hmgr->getHistoProf1(1410)->Fill(theData->getEta(), theData->getCoolingIL());
0355   hmgr->getHistoProf1(1420)->Fill(theData->getPhi(), theData->getCoolingIL());
0356   hmgr->getHistoProf2(1430)->Fill(theData->getEta(), theData->getPhi(), theData->getCoolingIL());
0357 
0358   // Electronics
0359   hmgr->getHisto1(1511)->Fill(theData->getEta());
0360   hmgr->getHisto1(1521)->Fill(theData->getPhi());
0361   hmgr->getHisto2(1531)->Fill(theData->getEta(), theData->getPhi());
0362 
0363   hmgr->getHistoProf1(1510)->Fill(theData->getEta(), theData->getElectronicsIL());
0364   hmgr->getHistoProf1(1520)->Fill(theData->getPhi(), theData->getElectronicsIL());
0365   hmgr->getHistoProf2(1530)->Fill(theData->getEta(), theData->getPhi(), theData->getElectronicsIL());
0366 
0367   // Other
0368   hmgr->getHisto1(1611)->Fill(theData->getEta());
0369   hmgr->getHisto1(1621)->Fill(theData->getPhi());
0370   hmgr->getHisto2(1631)->Fill(theData->getEta(), theData->getPhi());
0371 
0372   hmgr->getHistoProf1(1610)->Fill(theData->getEta(), theData->getOtherIL());
0373   hmgr->getHistoProf1(1620)->Fill(theData->getPhi(), theData->getOtherIL());
0374   hmgr->getHistoProf2(1630)->Fill(theData->getEta(), theData->getPhi(), theData->getOtherIL());
0375 }
0376 
0377 void MaterialBudgetMtdHistos::endOfRun() {
0378   edm::LogInfo("MaterialBudget") << "MaterialBudgetMtdHistos: Writing histos ROOT file to:" << theFileName;
0379   hmgr->save(theFileName);
0380 }