#include <DPGAnalysis/SiStripTools/plugins/APVShotsAnalyzer.cc>

Inheritance diagram for APVShotsAnalyzer:

Public Member Functions
	APVShotsAnalyzer (const edm::ParameterSet &)

	~APVShotsAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

virtual void	beginRun (const edm::Run &, const edm::EventSetup &)

virtual void	endJob ()

virtual void	endRun (const edm::Run &, const edm::EventSetup &)

void	updateDetCabling (const edm::EventSetup &setup)

Private Attributes
uint32_t	_cacheIdDet
	DB cache ID used to establish if the cabling has changed during the run. More...

const SiStripDetCabling *	_detCabling
	The cabling object. More...

edm::InputTag	_digicollection

TH1F *	_fed

TH1F *	_median

TH2F *	_medianVsFED

int	_nevents

TH1F *	_nShots

TH2F *	_nShotsVsFED

TProfile *	_nShotsVsTime

TH1F *	_stripMult

TH1F *	_subDetector

std::string	_suffix

bool	_useCabling

TH1F *	_whichAPV

bool	_zs

TkHistoMap *	tkhisto

TkHistoMap *	tkhisto2

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 66 of file APVShotsAnalyzer.cc.

Constructor & Destructor Documentation

APVShotsAnalyzer::APVShotsAnalyzer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 120 of file APVShotsAnalyzer.cc.

References _fed, _median, _medianVsFED, _nShots, _nShotsVsFED, _nShotsVsTime, _stripMult, _subDetector, _suffix, _useCabling, _whichAPV, _zs, TFileDirectory::make(), tkhisto, and tkhisto2.

                                                                 :
   _digicollection(iConfig.getParameter<edm::InputTag>("digiCollection")),
   _zs(iConfig.getUntrackedParameter<bool>("zeroSuppressed",true)),
   _suffix(iConfig.getParameter<std::string>("mapSuffix")),
   _nevents(0),
   _useCabling(iConfig.getUntrackedParameter<bool>("useCabling",true)),
   _cacheIdDet(0),
   _detCabling(0)
 {
    //now do what ever initialization is needed
 
   if(!_zs) _suffix += "_notZS";
 
  edm::Service<TFileService> tfserv;
 
  _nShots = tfserv->make<TH1F>("nShots","Number of Shots per event",200,-0.5,199.5);
  _nShots->GetXaxis()->SetTitle("Shots");  _nShots->GetYaxis()->SetTitle("Events"); 
  _nShots->StatOverflows(kTRUE);
 
  _nShotsVsTime = tfserv->make<TProfile>("nShotsVsTime","Mean number of shots vs orbit number",3600,0.5,3600*11223+0.5);
  _nShotsVsTime->GetXaxis()->SetTitle("Orbit");  _nShotsVsTime->GetYaxis()->SetTitle("Number of Shots");
  _nShotsVsTime->SetBit(TH1::kCanRebin);
 
  _whichAPV = tfserv->make<TH1F>("whichAPV","APV with shots",6,-0.5,5.5);
  _whichAPV->GetXaxis()->SetTitle("APV");  _whichAPV->GetYaxis()->SetTitle("Shots"); 
 
  _stripMult = tfserv->make<TH1F>("stripMultiplicity","Shot Strip Multiplicity",129,-0.5,128.5);
  _stripMult->GetXaxis()->SetTitle("Number of Strips");  _stripMult->GetYaxis()->SetTitle("Shots");
 
  _median = tfserv->make<TH1F>("median","APV Shot charge median",256,-0.5,255.5);
  _median->GetXaxis()->SetTitle("Charge [ADC]");  _median->GetYaxis()->SetTitle("Shots");
 
  _subDetector = tfserv->make<TH1F>("subDets","SubDetector Shot distribution",10,-0.5,9.5);
  _subDetector->GetYaxis()->SetTitle("Shots");
 
  if (_useCabling) {
    _fed = tfserv->make<TH1F>("fed","FED Shot distribution",440,50,490);
    _fed->GetYaxis()->SetTitle("Shots");
 
    _nShotsVsFED = tfserv->make<TH2F>("nShotsVsFED","Number of Shots per event vs fedid",440,50,490,200,-0.5,199.5);
    _nShotsVsFED->GetXaxis()->SetTitle("fedId");  _nShots->GetYaxis()->SetTitle("Shots");  _nShots->GetZaxis()->SetTitle("Events");
    _nShotsVsFED->StatOverflows(kTRUE);
 
    _medianVsFED = tfserv->make<TH2F>("medianVsFED","APV Shot charge median vs fedid",440,50,490,256,-0.5,255.5);
    _medianVsFED->GetXaxis()->SetTitle("fedId");_medianVsFED->GetYaxis()->SetTitle("Charge [ADC]");  _median->GetZaxis()->SetTitle("Shots");
  }
 
  tkhisto      =new TkHistoMap("ShotMultiplicity","ShotMultiplicity",-1); 
  tkhisto2      =new TkHistoMap("StripMultiplicity","StripMultiplicity",-1); 
 }

APVShotsAnalyzer::~APVShotsAnalyzer ( )

Definition at line 172 of file APVShotsAnalyzer.cc.

References _detCabling.

 {
  
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
   if ( _detCabling ) _detCabling = 0;
 
 }

Member Function Documentation

void APVShotsAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 188 of file APVShotsAnalyzer.cc.

References _detCabling, _digicollection, _fed, _median, _medianVsFED, _nevents, _nShots, _nShotsVsFED, _nShotsVsTime, _stripMult, _subDetector, _useCabling, _whichAPV, _zs, TkHistoMap::add(), gather_cfg::cout, sistrip::FED_ID_MAX, sistrip::FED_ID_MIN, TkHistoMap::fill(), edm::Event::getByLabel(), SiStripDetCabling::getConnections(), APVShotFinder::getShots(), edm::EventBase::orbitNumber(), tkhisto, tkhisto2, and updateDetCabling().

 {
    using namespace edm;
 
    if (_useCabling){
      //retrieve cabling
      updateDetCabling( iSetup );
    }
 
    _nevents++;
 
    Handle<edm::DetSetVector<SiStripDigi> > digis;
    iEvent.getByLabel(_digicollection,digis);
 
    // loop on detector with digis
 
    int nshots=0;
    std::vector<int> nshotsperFed;
 
    const uint16_t lNumFeds = sistrip::FED_ID_MAX-sistrip::FED_ID_MIN+1;
    if (_useCabling){
      nshotsperFed.resize(lNumFeds,0);
    }
 
    APVShotFinder apvsf(*digis,_zs);
    const std::vector<APVShot>& shots = apvsf.getShots();
 
    for(std::vector<APVShot>::const_iterator shot=shots.begin();shot!=shots.end();++shot) {
      if(shot->isGenuine()) {
 
        //get the fedid from the detid
 
        uint32_t det=shot->detId();
        if (_useCabling){
      const std::vector<FedChannelConnection> & conns = _detCabling->getConnections( det );
 
      if (!(conns.size())) continue;
      uint16_t lFedId = 0;
      for (uint32_t ch = 0; ch<conns.size(); ch++) {
        lFedId = conns[ch].fedId();
        //uint16_t lFedCh = conns[ch].fedCh();
    
        if (lFedId < sistrip::FED_ID_MIN || lFedId > sistrip::FED_ID_MAX){
          std::cout << " -- Invalid fedid " << lFedId << " for detid " << det << " connection " << ch << std::endl;
          continue;
        }
        else break;
      }
 
      if (lFedId < sistrip::FED_ID_MIN || lFedId > sistrip::FED_ID_MAX){
        std::cout << " -- No valid fedid (=" << lFedId << ") found for detid " << det << std::endl;
        continue;
      }
      ++nshotsperFed[lFedId-sistrip::FED_ID_MIN];
      _fed->Fill(lFedId);
      _medianVsFED->Fill(lFedId,shot->median());
 
 
        }
 
        ++nshots;
 
 
        _whichAPV->Fill(shot->apvNumber());
        _median->Fill(shot->median());
        _stripMult->Fill(shot->nStrips());
        _subDetector->Fill(shot->subDet());
        tkhisto2->fill(det,shot->nStrips());;
        tkhisto->add(det,1);
 
 
 
      }
    }
 
    _nShots->Fill(nshots);
    if (_useCabling){
      for (uint16_t lFed(0); lFed<lNumFeds; lFed++){
        _nShotsVsFED->Fill(lFed+sistrip::FED_ID_MIN,nshotsperFed[lFed]);
      }
    }
 
    _nShotsVsTime->Fill(iEvent.orbitNumber(),nshots);
    
 
 }

void APVShotsAnalyzer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 287 of file APVShotsAnalyzer.cc.

 {
 
 }

void APVShotsAnalyzer::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 276 of file APVShotsAnalyzer.cc.

277 {}

void APVShotsAnalyzer::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 294 of file APVShotsAnalyzer.cc.

References _nevents, _suffix, TkHistoMap::dumpInTkMap(), TrackerMap::save(), TkHistoMap::save(), TrackerMap::setPalette(), tkhisto, and tkhisto2.

                          {
 
   edm::LogInfo("EndOfJob") << _nevents << " analyzed events";
 
 #include "CommonTools/TrackerMap/interface/TrackerMap.h"
   TrackerMap tkmap,tkmap2;
 
   tkmap.setPalette(1);
   tkmap2.setPalette(1);
   tkhisto->dumpInTkMap(&tkmap);
   tkhisto2->dumpInTkMap(&tkmap2);
   std::string tkshotmultmapname = "ShotMultiplicity_" + _suffix + ".png";
   tkmap.save(true,0,0,tkshotmultmapname);
   std::string tkstripmultmapname = "StripMultiplicity_" + _suffix + ".png";
   tkmap2.save(true,0,0,tkstripmultmapname);
 
   std::string rootmapname = "TKMap_"+_suffix+".root";
   tkhisto->save(rootmapname);
   tkhisto2->save(rootmapname);
 }

void APVShotsAnalyzer::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 280 of file APVShotsAnalyzer.cc.

281 {

282 }

void APVShotsAnalyzer::updateDetCabling ( const edm::EventSetup & setup )

private

Definition at line 316 of file APVShotsAnalyzer.cc.

References _cacheIdDet, _detCabling, _useCabling, trackerHits::c, edm::EventSetup::get(), and edm::ESHandle< class >::product().

Referenced by analyze().

 {
   if (_useCabling){
     uint32_t cache_id = setup.get<SiStripDetCablingRcd>().cacheIdentifier();//.get( cabling_ );
    
     if ( _cacheIdDet != cache_id ) { // If the cache ID has changed since the last update...
       // Update the cabling object
       edm::ESHandle<SiStripDetCabling> c;
       setup.get<SiStripDetCablingRcd>().get( c );
       _detCabling = c.product();
       _cacheIdDet = cache_id;
     } // end of new cache ID check
   }
 }