#include <EcalPulseShapeGrapher.h>

Inheritance diagram for EcalPulseShapeGrapher:

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

	~EcalPulseShapeGrapher () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	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 &	itemsToGetFrom (BranchType iType) 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

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

void	endJob () override

std::string	intToString (int)

Private Attributes
int	abscissa [10]

int	ampCut_

std::map< int, TH1F * >	ampHistMap_

std::map< int, TH1F * >	cutAmpHistMap_

edm::InputTag	EBDigis_

edm::InputTag	EBUncalibratedRecHitCollection_

edm::InputTag	EEDigis_

edm::InputTag	EEUncalibratedRecHitCollection_

EcalFedMap *	fedMap_

TFile *	file_

std::map< int, TH1F * >	firstSampleHistMap_

std::vector< int >	listChannels_

int	ordinate [10]

std::map< int, TH2F * >	pulseShapeHistMap_

std::map< int, TH2F * >	rawPulseShapeHistMap_

std::string	rootFilename_

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 &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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

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

Definition at line 46 of file EcalPulseShapeGrapher.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 32 of file EcalPulseShapeGrapher.cc.

References abscissa, ampHistMap_, cutAmpHistMap_, fedMap_, firstSampleHistMap_, edm::ParameterSet::getUntrackedParameter(), mps_fire::i, intToString(), listChannels_, dataset::name, pulseShapeHistMap_, rawPulseShapeHistMap_, AlCaHLTBitMon_QueryRunRegistry::string, and fftjetcommon_cfi::title.

                                                                            :
 EBUncalibratedRecHitCollection_ (iConfig.getParameter<edm::InputTag>("EBUncalibratedRecHitCollection")),
 EBDigis_ (iConfig.getParameter<edm::InputTag>("EBDigiCollection")),
 EEUncalibratedRecHitCollection_ (iConfig.getParameter<edm::InputTag>("EEUncalibratedRecHitCollection")),
 EEDigis_ (iConfig.getParameter<edm::InputTag>("EEDigiCollection")),
 ampCut_ (iConfig.getUntrackedParameter<int>("AmplitudeCutADC", 13)),
 rootFilename_ (iConfig.getUntrackedParameter<std::string>("rootFilename","pulseShapeGrapher"))
 {
    //now do what ever initialization is needed
 
    std::vector<int> listDefaults;
    listDefaults.push_back(-1);
    listChannels_ = iConfig.getUntrackedParameter<std::vector<int> >("listChannels", listDefaults);
 
    for(std::vector<int>::const_iterator itr = listChannels_.begin(); itr != listChannels_.end(); ++itr)
    {
      std::string title = "Amplitude of cry "+intToString(*itr);
      std::string name = "ampOfCry"+intToString(*itr);
      ampHistMap_[*itr] = new TH1F(name.c_str(),title.c_str(),100,0,100);
      ampHistMap_[*itr]->GetXaxis()->SetTitle("ADC");
 
      title = "Amplitude (over 13 ADC) of cry "+intToString(*itr);
      name = "cutAmpOfCry"+intToString(*itr);
      cutAmpHistMap_[*itr] = new TH1F(name.c_str(),title.c_str(),100,0,100);
      cutAmpHistMap_[*itr]->GetXaxis()->SetTitle("ADC");
 
      title = "Pulse shape of cry "+intToString(*itr);
      name = "PulseShapeCry"+intToString(*itr);
      pulseShapeHistMap_[*itr] = new TH2F(name.c_str(),title.c_str(),10,0,10,220,-20,2);
      pulseShapeHistMap_[*itr]->GetXaxis()->SetTitle("sample");
 
      title = "Raw Pulse shape of cry "+intToString(*itr);
      name = "RawPulseShapeCry"+intToString(*itr);
      rawPulseShapeHistMap_[*itr] = new TH2F(name.c_str(),title.c_str(),10,0,10,500,0,500);
      rawPulseShapeHistMap_[*itr]->GetXaxis()->SetTitle("sample");
      
      title = "Amplitude of first sample, cry "+intToString(*itr);
      name = "AmpOfFirstSampleCry"+intToString(*itr);
      firstSampleHistMap_[*itr] = new TH1F(name.c_str(),title.c_str(),300,100,400);
      firstSampleHistMap_[*itr]->GetXaxis()->SetTitle("ADC");
    }
    
    fedMap_ = new EcalFedMap();
    
    for (int i=0; i<10; i++)
      abscissa[i] = i;
 }

EcalPulseShapeGrapher::~EcalPulseShapeGrapher ( )

override

Definition at line 81 of file EcalPulseShapeGrapher.cc.

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

Member Function Documentation

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

overrideprivate

Definition at line 96 of file EcalPulseShapeGrapher.cc.

References EcalMGPASample::adc(), ampCut_, ampHistMap_, CustomPhysics_cfi::amplitude, EcalUncalibratedRecHit::amplitude(), edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), cutAmpHistMap_, EcalElectronicsId::dccId(), EBDigis_, EBUncalibratedRecHitCollection_, EEDigis_, EEUncalibratedRecHitCollection_, edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), fedMap_, firstSampleHistMap_, muonCSCDigis_cfi::gain, EcalMGPASample::gainId(), edm::Event::getByLabel(), EcalFedMap::getSliceFromFed(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), mps_fire::i, EcalUncalibratedRecHit::id(), createfilelist::int, listChannels_, muonCSCDigis_cfi::pedestal, pulseShapeHistMap_, rawPulseShapeHistMap_, EcalDataFrame::sample(), and findQualityFiles::size.

 {
    using namespace edm;
    using namespace std;
 
    int numHitsWithActivity = 0;
    
    //int eventNum = iEvent.id().event();
    //vector<EcalUncalibratedRecHit> sampleHitsPastCut;
    
    Handle<EcalUncalibratedRecHitCollection> EBHits;
    iEvent.getByLabel(EBUncalibratedRecHitCollection_, EBHits);
    Handle<EcalUncalibratedRecHitCollection> EEHits;
    iEvent.getByLabel(EEUncalibratedRecHitCollection_, EEHits);
    //cout << "event: " << eventNum << " sample hits collection size: " << sampleHits->size() << endl;
    //Handle<EcalUncalibratedRecHitCollection> fittedHits;
    //iEvent.getByLabel(FittedUncalibratedRecHitCollection_, fittedHits);
    Handle<EBDigiCollection>  EBdigis;
    iEvent.getByLabel(EBDigis_, EBdigis);
    Handle<EEDigiCollection>  EEdigis;
    iEvent.getByLabel(EEDigis_, EEdigis);
    //cout << "event: " << eventNum << " digi collection size: " << digis->size() << endl;
 
    auto_ptr<EcalElectronicsMapping> ecalElectronicsMap(new EcalElectronicsMapping);
 
    // Loop over the hits
    for(EcalUncalibratedRecHitCollection::const_iterator hitItr = EBHits->begin(); hitItr != EBHits->end(); ++hitItr)
    {
      EcalUncalibratedRecHit hit = (*hitItr);
      float amplitude = hit.amplitude();
      EBDetId hitDetId = hit.id();
 
      // Get the Fedid
      EcalElectronicsId elecId = ecalElectronicsMap->getElectronicsId(hitDetId);
      int FEDid = 600+elecId.dccId();
      string SMname = fedMap_->getSliceFromFed(FEDid);
      
      vector<int>::const_iterator itr = listChannels_.begin();
      while(itr != listChannels_.end() && (*itr) != hitDetId.hashedIndex())
      {
        itr++;
      }
      if(itr==listChannels_.end())
        continue;
 
      ampHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
      //cout << "Cry hash:" << hitDetId.hashedIndex() << " amplitude: " << amplitude << endl;
      if(amplitude < ampCut_)
        continue;
 
      cutAmpHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
      numHitsWithActivity++;
      EBDigiCollection::const_iterator digiItr= EBdigis->begin();
      while(digiItr != EBdigis->end() && digiItr->id() != hitItr->id())
      {
        digiItr++;
      }
      if(digiItr == EBdigis->end())
        continue;
 
      double sampleADC[10];
      EBDataFrame df(*digiItr);
      double pedestal = 200;
 
      if(df.sample(0).gainId()!=1 || df.sample(1).gainId()!=1) continue; //goes to the next digi
      else {
        sampleADC[0] = df.sample(0).adc();
        sampleADC[1] = df.sample(1).adc();
        pedestal = (double)(sampleADC[0]+sampleADC[1])/(double)2;
      } 
 
      for (int i=0; (unsigned int)i< digiItr->size(); ++i ) {
        EBDataFrame df(*digiItr);
        double gain = 12.;
        if(df.sample(i).gainId()==1)
          gain = 1.;
        else if(df.sample(i).gainId()==2)
          gain = 2.;
        sampleADC[i] = pedestal+(df.sample(i).adc()-pedestal)*gain;
      }
 
      //cout << "1) maxsample amp:" << maxSampleAmp << " maxSampleIndex:" << maxSampleIndex << endl;
      for(int i=0; i<10;++i)
      {
        //cout << "Filling hist for:" << hitDetId.hashedIndex() << " with sample:" << i 
          //<< " amp:" <<(float)(df.sample(i).adc()-baseline)/(maxSampleAmp-baseline) << endl;
        //cout << "ADC of sample:" << df.sample(i).adc() << " baseline:" << baseline << " maxSampleAmp:" 
          //<< maxSampleAmp << endl << endl;
        pulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i]-pedestal)/amplitude);
        rawPulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i]));
      }
      firstSampleHistMap_[hitDetId.hashedIndex()]->Fill(sampleADC[0]);
    }
 
    // Now do the same for the EE hits
    for(EcalUncalibratedRecHitCollection::const_iterator hitItr = EEHits->begin(); hitItr != EEHits->end(); ++hitItr)
    {
      EcalUncalibratedRecHit hit = (*hitItr);
      float amplitude = hit.amplitude();
      EEDetId hitDetId = hit.id();
 
      // Get the Fedid
      EcalElectronicsId elecId = ecalElectronicsMap->getElectronicsId(hitDetId);
      int FEDid = 600+elecId.dccId();
      string SMname = fedMap_->getSliceFromFed(FEDid);
      
      vector<int>::const_iterator itr = listChannels_.begin();
      while(itr != listChannels_.end() && (*itr) != hitDetId.hashedIndex())
      {
        itr++;
      }
      if(itr==listChannels_.end())
        continue;
 
      ampHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
      //cout << "Cry hash:" << hitDetId.hashedIndex() << " amplitude: " << amplitude << endl;
      if(amplitude < ampCut_)
        continue;
 
      cutAmpHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
      numHitsWithActivity++;
      EEDigiCollection::const_iterator digiItr= EEdigis->begin();
      while(digiItr != EEdigis->end() && digiItr->id() != hitItr->id())
      {
        digiItr++;
      }
      if(digiItr == EEdigis->end())
        continue;
 
      double sampleADC[10];
      EEDataFrame df(*digiItr);
      double pedestal = 200;
 
      if(df.sample(0).gainId()!=1 || df.sample(1).gainId()!=1) continue; //goes to the next digi
      else {
        sampleADC[0] = df.sample(0).adc();
        sampleADC[1] = df.sample(1).adc();
        pedestal = (double)(sampleADC[0]+sampleADC[1])/(double)2;
      } 
 
      for (int i=0; (unsigned int)i< digiItr->size(); ++i ) {
        EEDataFrame df(*digiItr);
        double gain = 12.;
        if(df.sample(i).gainId()==1)
          gain = 1.;
        else if(df.sample(i).gainId()==2)
          gain = 2.;
        sampleADC[i] = pedestal+(df.sample(i).adc()-pedestal)*gain;
      }
 
      //cout << "1) maxsample amp:" << maxSampleAmp << " maxSampleIndex:" << maxSampleIndex << endl;
      for(int i=0; i<10;++i)
      {
        //cout << "Filling hist for:" << hitDetId.hashedIndex() << " with sample:" << i 
          //<< " amp:" <<(float)(df.sample(i).adc()-baseline)/(maxSampleAmp-baseline) << endl;
        //cout << "ADC of sample:" << df.sample(i).adc() << " baseline:" << baseline << " maxSampleAmp:" 
          //<< maxSampleAmp << endl << endl;
        pulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i]-pedestal)/amplitude);
        rawPulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i]));
      }
      firstSampleHistMap_[hitDetId.hashedIndex()]->Fill(sampleADC[0]);
    }
 }

void EcalPulseShapeGrapher::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 263 of file EcalPulseShapeGrapher.cc.

References ampHistMap_, cutAmpHistMap_, file_, firstSampleHistMap_, listChannels_, rawPulseShapeHistMap_, and rootFilename_.

Referenced by o2olib.O2ORunMgr::executeJob().

                               {
   
    rootFilename_+=".root";
    file_ = new TFile(rootFilename_.c_str(),"RECREATE");
    TH1::AddDirectory(false);
   
    for(std::vector<int>::const_iterator itr = listChannels_.begin(); itr != listChannels_.end(); ++itr)
    {
      ampHistMap_[*itr]->Write();
      cutAmpHistMap_[*itr]->Write();
      firstSampleHistMap_[*itr]->Write();
 
      rawPulseShapeHistMap_[*itr]->Write();
      TProfile* t2 = (TProfile*) (rawPulseShapeHistMap_[*itr]->ProfileX());
      t2->Write();
      //TODO: fix the normalization so these are correct
      //pulseShapeHistMap_[*itr]->Write();
      //TProfile* t1 = (TProfile*) (pulseShapeHistMap_[*itr]->ProfileX());
      //t1->Write();
    }
 
    
   file_->Write();
   file_->Close();
 }

std::string EcalPulseShapeGrapher::intToString ( int num )

private

Definition at line 289 of file EcalPulseShapeGrapher.cc.

Referenced by EcalPulseShapeGrapher().

 {
   using namespace std;
   ostringstream myStream;
   myStream << num << flush;
   return(myStream.str()); //returns the string form of the stringstream object
 }

Member Data Documentation

int EcalPulseShapeGrapher::abscissa[10]

private

Definition at line 63 of file EcalPulseShapeGrapher.h.

Referenced by EcalPulseShapeGrapher().

int EcalPulseShapeGrapher::ampCut_

private

Definition at line 72 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

std::map<int,TH1F*> EcalPulseShapeGrapher::ampHistMap_

private

Definition at line 66 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

std::map<int,TH1F*> EcalPulseShapeGrapher::cutAmpHistMap_

private

Definition at line 70 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

edm::InputTag EcalPulseShapeGrapher::EBDigis_

private

Definition at line 59 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

edm::InputTag EcalPulseShapeGrapher::EBUncalibratedRecHitCollection_

private

Definition at line 58 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

edm::InputTag EcalPulseShapeGrapher::EEDigis_

private

Definition at line 61 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

edm::InputTag EcalPulseShapeGrapher::EEUncalibratedRecHitCollection_

private

Definition at line 60 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

EcalFedMap* EcalPulseShapeGrapher::fedMap_

private

Definition at line 77 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), and EcalPulseShapeGrapher().

TFile* EcalPulseShapeGrapher::file_

private

Definition at line 75 of file EcalPulseShapeGrapher.h.

Referenced by endJob().

std::map<int,TH1F*> EcalPulseShapeGrapher::firstSampleHistMap_

private

Definition at line 68 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

std::vector<int> EcalPulseShapeGrapher::listChannels_

private

Definition at line 65 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

int EcalPulseShapeGrapher::ordinate[10]

private

Definition at line 64 of file EcalPulseShapeGrapher.h.

std::map<int,TH2F*> EcalPulseShapeGrapher::pulseShapeHistMap_

private

Definition at line 67 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), and EcalPulseShapeGrapher().

std::map<int,TH2F*> EcalPulseShapeGrapher::rawPulseShapeHistMap_

private

Definition at line 69 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

std::string EcalPulseShapeGrapher::rootFilename_

private

Definition at line 73 of file EcalPulseShapeGrapher.h.

Referenced by endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation