Inheritance diagram for DTTPAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

	DTTPAnalyzer (const edm::ParameterSet &)

void	endJob () override

virtual	~DTTPAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::string	getHistoName (const DTLayerId &)

Private Attributes
edm::InputTag	digiLabel_

edm::ESHandle< DTGeometry >	dtGeom_

std::map< DTWireId, int >	nDigisPerWire_

TFile *	rootFile_

bool	subtractT0_

std::map< DTWireId, double >	sumW2PerWire_

std::map< DTWireId, double >	sumWPerWire_

DTTTrigBaseSync *	tTrigSync_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Author: A. Vilela Pereira

Definition at line 20 of file DTTPAnalyzer.cc.

Constructor & Destructor Documentation

DTTPAnalyzer::DTTPAnalyzer ( const edm::ParameterSet & pset )

Definition at line 68 of file DTTPAnalyzer.cc.

References reco::get(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), rootFile_, DTAnalyzerDetailed_cfi::rootFileName, AlCaHLTBitMon_QueryRunRegistry::string, subtractT0_, and tTrigSync_.

                                                      :
   subtractT0_(pset.getParameter<bool>("subtractT0")),
   digiLabel_(pset.getParameter<edm::InputTag>("digiLabel")),
   tTrigSync_(0) {
 
   std::string rootFileName = pset.getUntrackedParameter<std::string>("rootFileName");
   rootFile_ = new TFile(rootFileName.c_str(), "RECREATE");
   rootFile_->cd();
 
   if(subtractT0_) 
      tTrigSync_ = DTTTrigSyncFactory::get()->create(pset.getParameter<std::string>("tTrigMode"),
                                                     pset.getParameter<edm::ParameterSet>("tTrigModeConfig"));
 
 }

DTTPAnalyzer::~DTTPAnalyzer ( )

virtual

Definition at line 83 of file DTTPAnalyzer.cc.

References rootFile_.

                            {  
   rootFile_->Close();
 }

Member Function Documentation

void DTTPAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

override

Definition at line 99 of file DTTPAnalyzer.cc.

References digiLabel_, nDigisPerWire_, DTTTrigBaseSync::offset(), PFRecoTauDiscriminationByIsolation_cfi::offset, subtractT0_, sumW2PerWire_, sumWPerWire_, cscNeutronWriter_cfi::t0, and tTrigSync_.

                                                                              {
 
   // Get the digis from the event
   edm::Handle<DTDigiCollection> digis; 
   event.getByLabel(digiLabel_, digis);
 
   // Iterate through all digi collections ordered by LayerId   
   DTDigiCollection::DigiRangeIterator dtLayerIt;
   for (dtLayerIt = digis->begin();
        dtLayerIt != digis->end();
        ++dtLayerIt){
     // Get the iterators over the digis associated with this LayerId
     const DTDigiCollection::Range& digiRange = (*dtLayerIt).second;
   
     // Get the layerId
     const DTLayerId layerId = (*dtLayerIt).first; //FIXME: check to be in the right sector
 
     // Loop over all digis in the given layer
     for (DTDigiCollection::const_iterator digi = digiRange.first;
      digi != digiRange.second;
      ++digi) {
        const DTWireId wireId( layerId, (*digi).wire() );
 
        double t0 = (*digi).countsTDC();
 
        //FIXME: Reject digis not coming from TP
 
        if(subtractT0_) {
           const DTLayer* layer = 0; //fake
       const GlobalPoint glPt; //fake
       double offset = tTrigSync_->offset(layer, wireId, glPt);
           t0 -= offset;
        }
 
        if(nDigisPerWire_.find(wireId) == nDigisPerWire_.end()){
           nDigisPerWire_[wireId] = 0;
           sumWPerWire_[wireId] = 0.;
           sumW2PerWire_[wireId] = 0.;  
        }
 
        ++nDigisPerWire_[wireId]; 
        sumWPerWire_[wireId] += t0;
        sumW2PerWire_[wireId] += t0*t0;
     }
 
   }
 }

void DTTPAnalyzer::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

override

Definition at line 87 of file DTTPAnalyzer.cc.

References dtGeom_, edm::EventSetup::get(), DTTTrigBaseSync::setES(), subtractT0_, and tTrigSync_.

                                                                          {
   // Get the t0 map from the DB
   if(subtractT0_){ 
      /*ESHandle<DTT0> t0;
      setup.get<DTT0Rcd>().get(t0);
      tZeroMap_ = &*t0;*/
      tTrigSync_->setES(setup);
   }
   // Get the DT Geometry  
   setup.get<MuonGeometryRecord>().get(dtGeom_);
 }

void DTTPAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 147 of file DTTPAnalyzer.cc.

References DTTopology::channels(), dtGeom_, DTTopology::firstChannel(), getHistoName(), crabWrapper::key, DTGeometry::layer(), DTWireId::layerId(), RecoTauDiscriminantConfiguration::mean, nDigisPerWire_, RecoTauDiscriminantConfiguration::rms, rootFile_, DTLayer::specificTopology(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, sumW2PerWire_, sumWPerWire_, and DTWireId::wire().

Referenced by o2olib.O2ORunMgr::executeJob().

                           {
   rootFile_->cd();
   std::map<DTLayerId, TH1F*> meanHistoMap;
   std::map<DTLayerId, TH1F*> sigmaHistoMap; 
   for(std::map<DTWireId, int>::const_iterator wireIdIt = nDigisPerWire_.begin();
                                               wireIdIt != nDigisPerWire_.end(); ++wireIdIt){
      DTWireId wireId((*wireIdIt).first);
 
      int nDigis = nDigisPerWire_[wireId];
      double sumW = sumWPerWire_[wireId];
      double sumW2 = sumW2PerWire_[wireId]; 
 
      double mean = sumW/nDigis;
      double rms = sumW2/nDigis - mean*mean;
      rms = sqrt(rms);
 
      DTLayerId layerId = wireId.layerId();
      if(meanHistoMap.find(layerId) == meanHistoMap.end()) {
         std::string histoName = getHistoName(layerId);
         const int firstChannel = dtGeom_->layer(layerId)->specificTopology().firstChannel();
         const int nWires = dtGeom_->layer(layerId)->specificTopology().channels();
         TH1F* meanHistoTP = new TH1F((histoName + "_tpMean").c_str(),"mean from test pulses by channel", 
                                       nWires,firstChannel,(firstChannel + nWires));
         TH1F* sigmaHistoTP = new TH1F((histoName + "_tpSigma").c_str(),"sigma from test pulses by channel",
                                       nWires,firstChannel,(firstChannel + nWires));
         meanHistoMap[layerId] = meanHistoTP;
         sigmaHistoMap[layerId] = sigmaHistoTP;
      }
      // Fill the histograms
      int nBin = meanHistoMap[layerId]->GetXaxis()->FindFixBin(wireId.wire());
      meanHistoMap[layerId]->SetBinContent(nBin,mean);
      sigmaHistoMap[layerId]->SetBinContent(nBin,rms);
   }
 
   for(std::map<DTLayerId, TH1F*>::const_iterator key = meanHistoMap.begin();
                                                  key != meanHistoMap.end(); ++key){
      meanHistoMap[(*key).first]->Write();
      sigmaHistoMap[(*key).first]->Write(); 
   } 
 
 }

std::string DTTPAnalyzer::getHistoName ( const DTLayerId & lId )

private

Definition at line 189 of file DTTPAnalyzer.cc.

References DEFINE_FWK_MODULE, DTLayerId::layer(), DTChamberId::sector(), DTChamberId::station(), AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

Referenced by endJob().

                                                          {
   std::string histoName;
   std::stringstream theStream;
   theStream << "Ch_" << lId.wheel() << "_" << lId.station() << "_" << lId.sector()
         << "_SL" << lId.superlayer() << "_L" << lId.layer();
   theStream >> histoName;
   return histoName;
 }

Member Data Documentation

edm::InputTag DTTPAnalyzer::digiLabel_

private

Definition at line 34 of file DTTPAnalyzer.cc.

Referenced by analyze().

edm::ESHandle<DTGeometry> DTTPAnalyzer::dtGeom_

private

Definition at line 38 of file DTTPAnalyzer.cc.

Referenced by beginRun(), and endJob().

std::map<DTWireId, int> DTTPAnalyzer::nDigisPerWire_

private

Definition at line 42 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

TFile* DTTPAnalyzer::rootFile_

private

Definition at line 36 of file DTTPAnalyzer.cc.

Referenced by DTTPAnalyzer(), endJob(), and ~DTTPAnalyzer().

bool DTTPAnalyzer::subtractT0_

private

Definition at line 33 of file DTTPAnalyzer.cc.

Referenced by analyze(), beginRun(), and DTTPAnalyzer().

std::map<DTWireId, double> DTTPAnalyzer::sumW2PerWire_

private

Definition at line 44 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

std::map<DTWireId, double> DTTPAnalyzer::sumWPerWire_

private

Definition at line 43 of file DTTPAnalyzer.cc.

Referenced by analyze(), and endJob().

DTTTrigBaseSync* DTTPAnalyzer::tTrigSync_

private

Definition at line 39 of file DTTPAnalyzer.cc.

Referenced by analyze(), beginRun(), and DTTPAnalyzer().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation