Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-03-17 11:27:27

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