CMS 3D CMS Logo

/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/DPGAnalysis/SiStripTools/src/DigiInvestigatorHistogramMaker.cc

Go to the documentation of this file.
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   //  char runname[100];
00115   //  sprintf(runname,"run_%d",nrun);
00116 
00117   edm::Service<TFileService> tfserv;
00118 
00119   //  currdir = &(*tfserv);
00120   //  _rhm.beginRun(nrun,*currdir);
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     //    char name[200];
00130     //    char title[500];
00131 
00132     //    TFileDirectory subd =_subdirs[i]->mkdir(runname);
00133 
00134     //    sprintf(name,"n%sdigivsorbrun",slab.c_str());
00135     //    sprintf(title,"%s %s mean multiplicity vs orbit",slab.c_str(),_hitname.c_str());
00136     //    _nmultvsorbrun[i] = subd.make<TProfile>(name,title,_norbbin,0.5,11223*_norbbin+0.5);
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