Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 14:32:39

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