00001 #include "DPGAnalysis/SiStripTools/interface/DigiInvestigatorHistogramMaker.h"
00002 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00003 #include "FWCore/ParameterSet/interface/ParameterSet.h"
00004 #include "FWCore/ServiceRegistry/interface/Service.h"
00005 #include "CommonTools/UtilAlgos/interface/TFileService.h"
00006 #include "TProfile.h"
00007 #include "TH1F.h"
00008
00009 #include "DPGAnalysis/SiStripTools/interface/SiStripTKNumbers.h"
00010
00011
00012 DigiInvestigatorHistogramMaker::DigiInvestigatorHistogramMaker():
00013 _hitname(), _nbins(500), m_maxLS(100), m_LSfrac(4), _scalefact(), _runHisto(true), _binmax(), _labels(), _nmultvsorbrun(), _nmult() { }
00014
00015 DigiInvestigatorHistogramMaker::DigiInvestigatorHistogramMaker(const edm::ParameterSet& iConfig):
00016 _hitname(iConfig.getUntrackedParameter<std::string>("hitName","digi")),
00017 _nbins(iConfig.getUntrackedParameter<int>("numberOfBins",500)),
00018 m_maxLS(iConfig.getUntrackedParameter<unsigned int>("maxLSBeforeRebin",100)),
00019 m_LSfrac(iConfig.getUntrackedParameter<unsigned int>("startingLSFraction",4)),
00020 _scalefact(iConfig.getUntrackedParameter<int>("scaleFactor",5)),
00021 _runHisto(iConfig.getUntrackedParameter<bool>("runHisto",true)),
00022 _labels(), _rhm(), _nmultvsorbrun(), _nmult(), _subdirs()
00023 {
00024
00025 std::vector<edm::ParameterSet>
00026 wantedsubds(iConfig.getUntrackedParameter<std::vector<edm::ParameterSet> >("wantedSubDets",std::vector<edm::ParameterSet>()));
00027
00028 for(std::vector<edm::ParameterSet>::iterator ps=wantedsubds.begin();ps!=wantedsubds.end();++ps) {
00029 _labels[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<std::string>("detLabel");
00030 _binmax[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<int>("binMax");
00031 }
00032
00033
00034 }
00035
00036
00037 DigiInvestigatorHistogramMaker::~DigiInvestigatorHistogramMaker() {
00038
00039 for(std::map<unsigned int,std::string>::const_iterator lab=_labels.begin();lab!=_labels.end();lab++) {
00040
00041 const unsigned int i = lab->first; const std::string slab = lab->second;
00042
00043 delete _subdirs[i];
00044 }
00045
00046 }
00047
00048
00049
00050 void DigiInvestigatorHistogramMaker::book(const std::string dirname, const std::map<unsigned int, std::string>& labels) {
00051
00052 _labels = labels;
00053 book(dirname);
00054
00055 }
00056
00057 void DigiInvestigatorHistogramMaker::book(const std::string dirname) {
00058
00059 edm::Service<TFileService> tfserv;
00060 TFileDirectory subev = tfserv->mkdir(dirname);
00061
00062 SiStripTKNumbers trnumb;
00063
00064 edm::LogInfo("NumberOfBins") << "Number of Bins: " << _nbins;
00065 edm::LogInfo("NumberOfMaxLS") << "Max number of LS before rebinning: " << m_maxLS;
00066 edm::LogInfo("StartingLSFrac") << "Fraction of LS in one bin before rebinning: " << m_LSfrac;
00067 edm::LogInfo("ScaleFactors") << "x-axis range scale factor: " << _scalefact;
00068 edm::LogInfo("BinMaxValue") << "Setting bin max values";
00069
00070 for(std::map<unsigned int,std::string>::const_iterator lab=_labels.begin();lab!=_labels.end();lab++) {
00071
00072 const unsigned int i = lab->first; const std::string slab = lab->second;
00073
00074 if(_binmax.find(i)==_binmax.end()) {
00075 edm::LogVerbatim("NotConfiguredBinMax") << "Bin max for " << lab->second
00076 << " not configured: " << trnumb.nstrips(i) << " used";
00077 _binmax[i] = trnumb.nstrips(i);
00078 }
00079
00080 edm::LogVerbatim("BinMaxValue") << "Bin max for " << lab->second << " is " << _binmax[i];
00081
00082 }
00083
00084 for(std::map<unsigned int,std::string>::const_iterator lab=_labels.begin();lab!=_labels.end();++lab) {
00085
00086 const int i = lab->first; const std::string slab = lab->second;
00087
00088 char name[200];
00089 char title[500];
00090
00091 _subdirs[i] = new TFileDirectory(subev.mkdir(slab.c_str()));
00092
00093 if(_subdirs[i]) {
00094 sprintf(name,"n%sdigi",slab.c_str());
00095 sprintf(title,"%s %s multiplicity",slab.c_str(),_hitname.c_str());
00096 _nmult[i] = _subdirs[i]->make<TH1F>(name,title,_nbins,0.,(1+_binmax[i]/(_scalefact*_nbins))*_nbins);
00097 _nmult[i]->GetXaxis()->SetTitle("Number of Hits"); _nmult[i]->GetYaxis()->SetTitle("Events");
00098
00099 if(_runHisto) {
00100 sprintf(name,"n%sdigivsorbrun",slab.c_str());
00101 sprintf(title,"%s %s mean multiplicity vs orbit",slab.c_str(),_hitname.c_str());
00102 _nmultvsorbrun[i] = _rhm.makeTProfile(name,title,m_LSfrac*m_maxLS,0,m_maxLS*262144);
00103 }
00104
00105 }
00106
00107 }
00108
00109
00110 }
00111
00112 void DigiInvestigatorHistogramMaker::beginRun(const unsigned int nrun) {
00113
00114
00115
00116
00117 edm::Service<TFileService> tfserv;
00118
00119
00120
00121
00122 _rhm.beginRun(nrun,*tfserv);
00123
00124
00125 for(std::map<unsigned int,std::string>::const_iterator lab=_labels.begin();lab!=_labels.end();++lab) {
00126
00127 const int i = lab->first; const std::string slab = lab->second;
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137 if(_runHisto) {
00138 (*_nmultvsorbrun[i])->GetXaxis()->SetTitle("time [orbit#]"); (*_nmultvsorbrun[i])->GetYaxis()->SetTitle("Hits");
00139 (*_nmultvsorbrun[i])->SetBit(TH1::kCanRebin);
00140 }
00141 }
00142
00143
00144 }
00145
00146 void DigiInvestigatorHistogramMaker::fill(const unsigned int orbit, const std::map<unsigned int,int>& ndigi) {
00147
00148 for(std::map<unsigned int,int>::const_iterator digi=ndigi.begin();digi!=ndigi.end();digi++) {
00149
00150 if(_labels.find(digi->first) != _labels.end()) {
00151
00152 const unsigned int i=digi->first;
00153
00154 _nmult[i]->Fill(digi->second);
00155 if(_runHisto) {
00156 if(_nmultvsorbrun[i] && *_nmultvsorbrun[i]) (*_nmultvsorbrun[i])->Fill(orbit,digi->second);
00157 }
00158 }
00159
00160 }
00161 }
00162