CMS 3D CMS Logo

MaterialBudgetEcalHistos Class Reference

#include <Validation/Geometry/interface/MaterialBudgetEcalHistos.h>

Inheritance diagram for MaterialBudgetEcalHistos:

MaterialBudgetFormat

List of all members.

Public Member Functions

virtual void fillEndTrack ()
virtual void fillPerStep ()
virtual void fillStartTrack ()
 MaterialBudgetEcalHistos (MaterialBudgetData *data, TestHistoMgr *mgr, const std::string &fileName)
virtual ~MaterialBudgetEcalHistos ()

Private Member Functions

virtual void book ()
virtual void hend ()

Private Attributes

TestHistoMgr * hmgr
int MAXNUMBERSTEPS
double * theDmb
double * theMateId
double * theVoluId
double * theX
double * theY
double * theZ

Detailed Description

Definition at line 9 of file MaterialBudgetEcalHistos.h.

Constructor & Destructor Documentation

MaterialBudgetEcalHistos::MaterialBudgetEcalHistos ( MaterialBudgetData *  data,
TestHistoMgr *  mgr,
const std::string &  fileName 
)

Definition at line 7 of file MaterialBudgetEcalHistos.cc.

References book(), and MaterialBudgetFormat::theFileName.

00009                                                                               : MaterialBudgetFormat( data ), hmgr(mgr)
00010 {
00011   theFileName = fileName;
00012   book();
00013 
00014 }

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

Definition at line 16 of file MaterialBudgetEcalHistos.h.

References hend().

00016 { hend(); }

Member Function Documentation

void MaterialBudgetEcalHistos::book (  )  [private, virtual]

Definition at line 17 of file MaterialBudgetEcalHistos.cc.

References TestHistoMgr::addHisto1(), TestHistoMgr::addHisto2(), TestHistoMgr::addHistoProf1(), TestHistoMgr::addHistoProf2(), GenMuonPlsPt100GeV_cfg::cout, lat::endl(), and hmgr.

Referenced by MaterialBudgetEcalHistos().

00018 {
00019   std::cout << "=== booking user histos ===" << std::endl;
00020   
00021   // total X0
00022   hmgr->addHistoProf1( new TProfile("10", "MB prof Eta ", 250, -5., 5. ) );
00023   hmgr->addHisto1( new TH1F("11", "Eta " , 501, -5., 5. ) );
00024   hmgr->addHistoProf1( new TProfile("20", "MB prof Phi ", 180, -3.1416, 3.1416 ) );
00025   hmgr->addHisto1( new TH1F("21", "Phi " , 180, -3.1416, 3.1416 ) );
00026   hmgr->addHistoProf2( new TProfile2D("30", "MB prof Eta  Phi ", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00027   hmgr->addHisto2( new TH2F("31", "Eta vs Phi " , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00028   
00029   // rr
00030   
00031   // Support
00032   hmgr->addHistoProf1( new TProfile("110", "MB prof Eta [Support]", 250, -5.0, 5.0 ) );
00033   hmgr->addHisto1( new TH1F("111", "Eta [Support]" , 501, -5., 5. ) );
00034   hmgr->addHistoProf1( new TProfile("120", "MB prof Phi [Support]", 180, -3.1416, 3.1416 ) );
00035   hmgr->addHisto1( new TH1F("121", "Phi [Support]" , 180, -3.1416, 3.1416 ) );
00036   hmgr->addHistoProf2( new TProfile2D("130", "MB prof Eta  Phi [Support]", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00037   hmgr->addHisto2( new TH2F("131", "Eta vs Phi [Support]" , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00038   // Sensitive
00039   hmgr->addHistoProf1( new TProfile("210", "MB prof Eta [Sensitive]", 250, -5.0, 5.0 ) );
00040   hmgr->addHisto1( new TH1F("211", "Eta [Sensitive]" , 501, -5., 5. ) );
00041   hmgr->addHistoProf1( new TProfile("220", "MB prof Phi [Sensitive]", 180, -3.1416, 3.1416 ) );
00042   hmgr->addHisto1( new TH1F("221", "Phi [Sensitive]" , 180, -3.1416, 3.1416 ) );
00043   hmgr->addHistoProf2( new TProfile2D("230", "MB prof Eta  Phi [Sensitive]", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00044   hmgr->addHisto2( new TH2F("231", "Eta vs Phi [Sensitive]" , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00045   // Cables
00046   hmgr->addHistoProf1( new TProfile("310", "MB prof Eta [Cables]", 250, -5.0, 5.0 ) );
00047   hmgr->addHisto1( new TH1F("311", "Eta [Cables]" , 501, -5., 5. ) );
00048   hmgr->addHistoProf1( new TProfile("320", "MB prof Phi [Cables]", 180, -3.1416, 3.1416 ) );
00049   hmgr->addHisto1( new TH1F("321", "Phi [Cables]" , 180, -3.1416, 3.1416 ) );
00050   hmgr->addHistoProf2( new TProfile2D("330", "MB prof Eta  Phi [Cables]", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00051   hmgr->addHisto2( new TH2F("331", "Eta vs Phi [Cables]" , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00052   // Cooling
00053   hmgr->addHistoProf1( new TProfile("410", "MB prof Eta [Cooling]", 250, -5.0, 5.0 ) );
00054   hmgr->addHisto1( new TH1F("411", "Eta [Cooling]" , 501, -5., 5. ) );
00055   hmgr->addHistoProf1( new TProfile("420", "MB prof Phi [Cooling]", 180, -3.1416, 3.1416 ) );
00056   hmgr->addHisto1( new TH1F("421", "Phi [Cooling]" , 180, -3.1416, 3.1416 ) );
00057   hmgr->addHistoProf2( new TProfile2D("430", "MB prof Eta  Phi [Cooling]", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00058   hmgr->addHisto2( new TH2F("431", "Eta vs Phi [Cooling]" , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00059   // Electronics
00060   hmgr->addHistoProf1( new TProfile("510", "MB prof Eta [Electronics]", 250, -5.0, 5.0 ) );
00061   hmgr->addHisto1( new TH1F("511", "Eta [Electronics]" , 501, -5., 5. ) );
00062   hmgr->addHistoProf1( new TProfile("520", "MB prof Phi [Electronics]", 180, -3.1416, 3.1416 ) );
00063   hmgr->addHisto1( new TH1F("521", "Phi [Electronics]" , 180, -3.1416, 3.1416 ) );
00064   hmgr->addHistoProf2( new TProfile2D("530", "MB prof Eta  Phi [Electronics]", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00065   hmgr->addHisto2( new TH2F("531", "Eta vs Phi [Electronics]" , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00066   // Other
00067   hmgr->addHistoProf1( new TProfile("610", "MB prof Eta [Other]", 250, -5.0, 5.0 ) );
00068   hmgr->addHisto1( new TH1F("611", "Eta [Other]" , 501, -5., 5. ) );
00069   hmgr->addHistoProf1( new TProfile("620", "MB prof Phi [Other]", 180, -3.1416, 3.1416 ) );
00070   hmgr->addHisto1( new TH1F("621", "Phi [Other]" , 180, -3.1416, 3.1416 ) );
00071   hmgr->addHistoProf2( new TProfile2D("630", "MB prof Eta  Phi [Other]", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00072   hmgr->addHisto2( new TH2F("631", "Eta vs Phi [Other]" , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00073   // Air
00074   hmgr->addHistoProf1( new TProfile("710", "MB prof Eta [Air]", 250, -5.0, 5.0 ) );
00075   hmgr->addHisto1( new TH1F("711", "Eta [Air]" , 501, -5., 5. ) );
00076   hmgr->addHistoProf1( new TProfile("720", "MB prof Phi [Air]", 180, -3.1416, 3.1416 ) );
00077   hmgr->addHisto1( new TH1F("721", "Phi [Air]" , 180, -3.1416, 3.1416 ) );
00078   hmgr->addHistoProf2( new TProfile2D("730", "MB prof Eta  Phi [Air]", 250, -5., 5., 180, -3.1416, 3.1416 ) );
00079   hmgr->addHisto2( new TH2F("731", "Eta vs Phi [Air]" , 501, -5., 5., 180, -3.1416, 3.1416 ) );
00080   //
00081 
00082   // ECAL specific
00083 
00084   hmgr->addHistoProf1( new TProfile("1001", "MB prof Eta ECAL Barrel", 340, -1.5, 1.5 ) );
00085   hmgr->addHistoProf1( new TProfile("1002", "MB prof Phi ECAL Barrel", 180, -3.1416, 3.1416 ) );
00086   hmgr->addHistoProf1( new TProfile("1003", "MB prof Phi ECAL Barrel SM", 20, 0., 20. ) );
00087   hmgr->addHistoProf1( new TProfile("2003", "MB prof Phi ECAL Barrel SM", 10, 0., 20. ) );
00088   hmgr->addHistoProf1( new TProfile("1004", "MB prof Phi ECAL Barrel SM module 1", 20, 0., 20. ) );
00089   hmgr->addHistoProf1( new TProfile("1005", "MB prof Phi ECAL Barrel SM module 2", 20, 0., 20. ) );
00090   hmgr->addHistoProf1( new TProfile("1006", "MB prof Phi ECAL Barrel SM module 3", 20, 0., 20. ) );
00091   hmgr->addHistoProf1( new TProfile("1007", "MB prof Phi ECAL Barrel SM module 4", 20, 0., 20. ) );
00092 
00093   hmgr->addHistoProf1( new TProfile("1011", "MB prof Eta ECAL Preshower +", 100, 1.65, 2.6 ) );
00094   hmgr->addHistoProf1( new TProfile("1012", "MB prof Phi ECAL Preshower +", 180, -3.1416, 3.1416 ) );
00095   
00096   hmgr->addHistoProf1( new TProfile("1013", "MB prof Eta ECAL Preshower -", 100, -2.6, -1.65 ) );
00097   hmgr->addHistoProf1( new TProfile("1014", "MB prof Phi ECAL Preshower -", 180, -3.1416, 3.1416 ) );
00098   
00099   
00100   // rr
00101   
00102   std::cout << "=== booking user histos done ===" << std::endl;
00103 
00104 }

void MaterialBudgetEcalHistos::fillEndTrack (  )  [virtual]

Implements MaterialBudgetFormat.

Definition at line 119 of file MaterialBudgetEcalHistos.cc.

References MaterialBudgetData::getEta(), TestHistoMgr::getHisto1(), TestHistoMgr::getHisto2(), TestHistoMgr::getHistoProf1(), TestHistoMgr::getHistoProf2(), MaterialBudgetData::getPhi(), MaterialBudgetData::getTotalMB(), hmgr, phi, and MaterialBudgetFormat::theData.

00120 {
00121   // Total X0
00122   hmgr->getHisto1(11)->Fill(theData->getEta());
00123   hmgr->getHisto1(21)->Fill(theData->getPhi());
00124   hmgr->getHisto2(31)->Fill(theData->getEta(),theData->getPhi());
00125   
00126   hmgr->getHistoProf1(10)->Fill(theData->getEta(),theData->getTotalMB());
00127   hmgr->getHistoProf1(20)->Fill(theData->getPhi(),theData->getTotalMB());
00128   hmgr->getHistoProf2(30)->Fill(theData->getEta(),theData->getPhi(),theData->getTotalMB());
00129   
00130   
00131   // ECAL specific
00132 
00133   if (fabs(theData->getEta()) <= 1.479 ) {
00134     static const double twenty ( 20.*degree ) ;
00135     const double phi ( theData->getPhi()+M_PI ) ;
00136     const double phiModTwenty (( phi - floor(phi/twenty)*twenty )/degree) ;
00137     hmgr->getHistoProf1(1001)->Fill(theData->getEta(),theData->getTotalMB());
00138     hmgr->getHistoProf1(1002)->Fill(theData->getPhi(),theData->getTotalMB());
00139     hmgr->getHistoProf1(1003)->Fill(phiModTwenty,theData->getTotalMB());
00140     hmgr->getHistoProf1(2003)->Fill(phiModTwenty,theData->getTotalMB());
00141     if (fabs(theData->getEta()) >= 0. && fabs(theData->getEta()) < 0.435 ) {
00142       hmgr->getHistoProf1(1004)->Fill(phiModTwenty,theData->getTotalMB());
00143     }
00144     if (fabs(theData->getEta()) >= 0.435 && fabs(theData->getEta()) < 0.783 ) {
00145       hmgr->getHistoProf1(1005)->Fill(phiModTwenty,theData->getTotalMB());
00146     }
00147     if (fabs(theData->getEta()) > 0.783 && fabs(theData->getEta()) <= 1.131 ) {
00148       hmgr->getHistoProf1(1006)->Fill(phiModTwenty,theData->getTotalMB());
00149     }
00150     if (fabs(theData->getEta()) > 1.131 && fabs(theData->getEta()) <= 1.479 ) {
00151       hmgr->getHistoProf1(1007)->Fill(phiModTwenty,theData->getTotalMB());
00152     }
00153     
00154   }
00155 
00156   if (theData->getEta() >= 1.653 && theData->getEta() <= 2.6 ) {
00157     hmgr->getHistoProf1(1011)->Fill(theData->getEta(),theData->getTotalMB());
00158     hmgr->getHistoProf1(1012)->Fill(theData->getPhi(),theData->getTotalMB());
00159   }
00160 
00161   if (theData->getEta() >= -2.6 && theData->getEta() <= -1.653 ) {
00162     hmgr->getHistoProf1(1013)->Fill(theData->getEta(),theData->getTotalMB());
00163     hmgr->getHistoProf1(1014)->Fill(theData->getPhi(),theData->getTotalMB());
00164   }
00165 
00166 
00167   // rr
00168 }

void MaterialBudgetEcalHistos::fillPerStep (  )  [virtual]

Implements MaterialBudgetFormat.

Definition at line 113 of file MaterialBudgetEcalHistos.cc.

00114 {
00115 
00116 }

void MaterialBudgetEcalHistos::fillStartTrack (  )  [virtual]

Implements MaterialBudgetFormat.

Definition at line 107 of file MaterialBudgetEcalHistos.cc.

00108 {
00109 
00110 }

void MaterialBudgetEcalHistos::hend (  )  [private, virtual]

Definition at line 171 of file MaterialBudgetEcalHistos.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), hmgr, TestHistoMgr::save(), and MaterialBudgetFormat::theFileName.

Referenced by ~MaterialBudgetEcalHistos().

00172 {
00173   std::cout << "=== save user histos ===" << std::endl;
00174   hmgr->save( theFileName );
00175 
00176 }

Member Data Documentation

TestHistoMgr* MaterialBudgetEcalHistos::hmgr [private]

Definition at line 37 of file MaterialBudgetEcalHistos.h.

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

int MaterialBudgetEcalHistos::MAXNUMBERSTEPS [private]

Definition at line 29 of file MaterialBudgetEcalHistos.h.

double* MaterialBudgetEcalHistos::theDmb [private]

Definition at line 30 of file MaterialBudgetEcalHistos.h.

double* MaterialBudgetEcalHistos::theMateId [private]

Definition at line 35 of file MaterialBudgetEcalHistos.h.

double* MaterialBudgetEcalHistos::theVoluId [private]

Definition at line 34 of file MaterialBudgetEcalHistos.h.

double* MaterialBudgetEcalHistos::theX [private]

Definition at line 31 of file MaterialBudgetEcalHistos.h.

double* MaterialBudgetEcalHistos::theY [private]

Definition at line 32 of file MaterialBudgetEcalHistos.h.

double* MaterialBudgetEcalHistos::theZ [private]

Definition at line 33 of file MaterialBudgetEcalHistos.h.

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:28:03 2009 for CMSSW by  doxygen 1.5.4