CMS 3D CMS Logo

/data/refman/pasoursint/CMSSW_4_1_8_patch13/src/EventFilter/SiPixelRawToDigi/plugins/SiPixelRawToDigi.cc

Go to the documentation of this file.
00001 // Include parameter driven interface to SiPixelQuality for study purposes
00002 // exclude ROC(raw) based on bad ROC list in SiPixelQuality
00003 // enabled by: process.siPixelDigis.UseQualityInfo = True
00004 // 20-10-2010 Andrew York (Tennessee)
00005 
00006 #include "SiPixelRawToDigi.h"
00007 
00008 #include "DataFormats/Common/interface/Handle.h"
00009 #include "FWCore/Framework/interface/ESHandle.h"
00010 #include "FWCore/Framework/interface/ESTransientHandle.h"
00011 
00012 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00013 
00014 #include "DataFormats/SiPixelDigi/interface/PixelDigi.h"
00015 
00016 #include "DataFormats/SiPixelRawData/interface/SiPixelRawDataError.h"
00017 
00018 #include "DataFormats/Common/interface/DetSetVector.h"
00019 #include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
00020 #include "DataFormats/FEDRawData/interface/FEDRawData.h"
00021 
00022 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
00023 
00024 #include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingMap.h"
00025 #include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingTree.h"
00026 #include "EventFilter/SiPixelRawToDigi/interface/PixelDataFormatter.h"
00027 
00028 #include "CondFormats/SiPixelObjects/interface/SiPixelQuality.h"
00029 
00030 #include "EventFilter/SiPixelRawToDigi/interface/R2DTimerObserver.h"
00031 
00032 #include "TH1D.h"
00033 #include "TFile.h"
00034 
00035 using namespace std;
00036 
00037 // -----------------------------------------------------------------------------
00038 SiPixelRawToDigi::SiPixelRawToDigi( const edm::ParameterSet& conf ) 
00039   : config_(conf), 
00040     cabling_(0), 
00041     badPixelInfo_(0), 
00042     hCPU(0), hDigi(0), theTimer(0)
00043 {
00044 
00045   includeErrors = config_.getParameter<bool>("IncludeErrors");
00046   useQuality = config_.getParameter<bool>("UseQualityInfo");
00047   useCablingTree_ = config_.getUntrackedParameter<bool>("UseCablingTree",true);
00048 
00049   //start counters
00050   ndigis = 0;
00051   nwords = 0;
00052 
00053   // Products
00054   produces< edm::DetSetVector<PixelDigi> >();
00055   if(includeErrors) produces< edm::DetSetVector<SiPixelRawDataError> >();
00056 
00057   // Timing
00058   bool timing = config_.getUntrackedParameter<bool>("Timing",false);
00059   if (timing) {
00060     theTimer = new R2DTimerObserver("**** MY TIMING REPORT ***");
00061     hCPU = new TH1D ("hCPU","hCPU",100,0.,0.050);
00062     hDigi = new TH1D("hDigi","hDigi",50,0.,15000.);
00063   }
00064 }
00065 
00066 
00067 // -----------------------------------------------------------------------------
00068 SiPixelRawToDigi::~SiPixelRawToDigi() {
00069   edm::LogInfo("SiPixelRawToDigi")  << " HERE ** SiPixelRawToDigi destructor!";
00070 
00071   if(useCablingTree_) delete cabling_;
00072 
00073   if (theTimer) {
00074     TFile rootFile("analysis.root", "RECREATE", "my histograms");
00075     hCPU->Write();
00076     hDigi->Write();
00077     delete theTimer;
00078   }
00079 
00080 }
00081 
00082 
00083 // -----------------------------------------------------------------------------
00084 
00085 
00086 // -----------------------------------------------------------------------------
00087 void SiPixelRawToDigi::produce( edm::Event& ev,
00088                               const edm::EventSetup& es) 
00089 {
00090   const uint32_t dummydetid = 0xffffffff;
00091   debug = edm::MessageDrop::instance()->debugEnabled;
00092 
00093 // initialize cabling map or update if necessary
00094   if (recordWatcher.check( es )) {
00095     // cabling map, which maps online address (fed->link->ROC->local pixel) to offline (DetId->global pixel)
00096     edm::ESTransientHandle<SiPixelFedCablingMap> cablingMap;
00097     es.get<SiPixelFedCablingMapRcd>().get( cablingMap );
00098     fedList = cablingMap->fedIds();
00099     if (useCablingTree_ && cabling_) delete cabling_; 
00100     if (useCablingTree_) cabling_ = cablingMap->cablingTree();
00101     else cabling_ = cablingMap.product();
00102     LogDebug("map version:")<< cabling_->version();
00103   }
00104 // initialize quality record or update if necessary
00105   if (qualityWatcher.check( es )&&useQuality) {
00106     // quality info for dead pixel modules or ROCs
00107     edm::ESHandle<SiPixelQuality> qualityInfo;
00108     es.get<SiPixelQualityRcd>().get( qualityInfo );
00109     if (badPixelInfo_) delete badPixelInfo_;
00110     badPixelInfo_ = qualityInfo.product();
00111     if (!badPixelInfo_) {
00112       edm::LogError("**SiPixelRawToDigi**")<<" Configured to use SiPixelQuality, but SiPixelQuality not present"<<endl;
00113     }
00114   }
00115 
00116   edm::Handle<FEDRawDataCollection> buffers;
00117   label = config_.getParameter<edm::InputTag>("InputLabel");
00118   ev.getByLabel( label, buffers);
00119 
00120 // create product (digis & errors)
00121   std::auto_ptr< edm::DetSetVector<PixelDigi> > collection( new edm::DetSetVector<PixelDigi> );
00122   std::auto_ptr< edm::DetSetVector<SiPixelRawDataError> > errorcollection( new edm::DetSetVector<SiPixelRawDataError> );
00123 
00124   PixelDataFormatter formatter(cabling_);
00125   formatter.setErrorStatus(includeErrors);
00126   if (useQuality) formatter.setQualityStatus(useQuality, badPixelInfo_);
00127 
00128   if (theTimer) theTimer->start();
00129   bool errorsInEvent = false;
00130   PixelDataFormatter::DetErrors nodeterrors;
00131 
00132   typedef std::vector<unsigned int>::const_iterator IF;
00133   for (IF aFed = fedList.begin(); aFed != fedList.end(); ++aFed) {
00134     int fedId = *aFed;
00135     if(debug) LogDebug("SiPixelRawToDigi")<< " PRODUCE DIGI FOR FED: " <<  fedId << endl;
00136     PixelDataFormatter::Digis digis;
00137     PixelDataFormatter::Errors errors;
00138 
00139     //get event data for this fed
00140     const FEDRawData& fedRawData = buffers->FEDData( fedId );
00141 
00142     //convert data to digi and strip off errors
00143     formatter.interpretRawData( errorsInEvent, fedId, fedRawData, digis, errors);
00144 
00145     //pack digi into collection
00146     typedef PixelDataFormatter::Digis::iterator ID;
00147     for (ID it = digis.begin(); it != digis.end(); it++) {
00148       uint32_t detid = it->first;
00149       edm::DetSet<PixelDigi>& detSet = collection->find_or_insert(detid);
00150       detSet.data = it->second;
00151     }
00152 
00153     //pack errors into collection
00154     if(includeErrors) {
00155       typedef PixelDataFormatter::Errors::iterator IE;
00156       for (IE is = errors.begin(); is != errors.end(); is++) {
00157         uint32_t errordetid = is->first;
00158         if (errordetid==dummydetid) {           // errors given dummy detId must be sorted by Fed
00159           nodeterrors.insert( nodeterrors.end(), errors[errordetid].begin(), errors[errordetid].end() );
00160         } else {
00161           edm::DetSet<SiPixelRawDataError>& errorDetSet = errorcollection->find_or_insert(errordetid);
00162           errorDetSet.data = is->second;
00163         }
00164       }
00165     }
00166   }
00167 
00168   if(includeErrors) {
00169     edm::DetSet<SiPixelRawDataError>& errorDetSet = errorcollection->find_or_insert(dummydetid);
00170     errorDetSet.data = nodeterrors;
00171   }
00172   if (errorsInEvent) LogDebug("SiPixelRawToDigi") << "Error words were stored in this event";
00173 
00174   if (theTimer) {
00175     theTimer->stop();
00176     LogDebug("SiPixelRawToDigi") << "TIMING IS: (real)" << theTimer->lastMeasurement().real() ;
00177     ndigis += formatter.nDigis();
00178     nwords += formatter.nWords();
00179     LogDebug("SiPixelRawToDigi") << " (Words/Digis) this ev: "
00180          <<formatter.nWords()<<"/"<<formatter.nDigis() << "--- all :"<<nwords<<"/"<<ndigis;
00181     hCPU->Fill( theTimer->lastMeasurement().real() ); 
00182     hDigi->Fill(formatter.nDigis());
00183   }
00184 
00185   //send digis and errors back to framework 
00186   ev.put( collection );
00187   if(includeErrors) ev.put( errorcollection );
00188 }