#include <L1Trigger/L1RCTTPGProvider/src/L1RCTTPGProvider.cc>

Inheritance diagram for L1RCTTPGProvider:

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

	~L1RCTTPGProvider ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	endJob ()

virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
edm::InputTag	ecalTPG_

int	hbShift

edm::InputTag	hcalTPG_

int	hfShift

int	postSamples

int	preSamples

bool	useEcalCosmicTiming

bool	useHcalCosmicTiming

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

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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 36 of file L1RCTTPGProvider.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 5 of file L1RCTTPGProvider.cc.

References i, postSamples, and preSamples.

                                                                 :
   ecalTPG_(iConfig.getParameter<edm::InputTag>("ecalTPGs")),
   hcalTPG_(iConfig.getParameter<edm::InputTag>("hcalTPGs")),
   useHcalCosmicTiming(iConfig.getParameter<bool>("useHCALCosmicTiming")),
   useEcalCosmicTiming(iConfig.getParameter<bool>("useECALCosmicTiming")),
   preSamples(iConfig.getParameter<int>("preSamples")),
   postSamples(iConfig.getParameter<int>("postSamples")),
   hfShift(iConfig.getParameter<int>("HFShift")),
   hbShift(iConfig.getParameter<int>("HBShift"))
 {
 
   //Output :The new manipulated TPGs
   //make it smart - to name the collections
   //correctly
     char ecal_label[200];
     char hcal_label[200];
 
 
 
     for(int i=preSamples;i>0;--i)
       {
         sprintf(ecal_label,"ECALBxminus%d",i);
         sprintf(hcal_label,"HCALBxminus%d",i);
         produces<EcalTrigPrimDigiCollection>(ecal_label);
         produces<HcalTrigPrimDigiCollection>(hcal_label);
       }
     
     produces<EcalTrigPrimDigiCollection>("ECALBx0");
     produces<HcalTrigPrimDigiCollection>("HCALBx0");
 
     for(int i=0;i<postSamples;++i)
       {
        sprintf(ecal_label,"ECALBxplus%d",i+1);
        sprintf(hcal_label,"HCALBxplus%d",i+1);
        produces<EcalTrigPrimDigiCollection>(ecal_label);
        produces<HcalTrigPrimDigiCollection>(hcal_label);
       }
 }

L1RCTTPGProvider::~L1RCTTPGProvider ( )

Definition at line 44 of file L1RCTTPGProvider.cc.

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

Member Function Documentation

void L1RCTTPGProvider::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 404 of file L1RCTTPGProvider.cc.

404 {

405 }

void L1RCTTPGProvider::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 57 of file L1RCTTPGProvider.cc.

References patCandidatesForDimuonsSequences_cff::ecal, ecalTPG_, EcalTriggerTower, edm::Event::getByLabel(), hbShift, patCandidatesForDimuonsSequences_cff::hcal, hcalTPG_, hfShift, i, j, postSamples, preSamples, edm::Event::put(), compare_using_db::sample, HcalTriggerPrimitiveDigi::setPresamples(), EcalTriggerPrimitiveDigi::setSample(), HcalTriggerPrimitiveDigi::setSample(), EcalTriggerPrimitiveDigi::setSize(), HcalTriggerPrimitiveDigi::setSize(), useEcalCosmicTiming, useHcalCosmicTiming, and ecaldqm::zside().

 {
    using namespace edm;
    //Declare handles
    Handle<EcalTrigPrimDigiCollection> ecal;
    Handle<HcalTrigPrimDigiCollection> hcal;
 
    //Declare vector of collection to send for output !
 
    std::vector<EcalTrigPrimDigiCollection> ecalColl(preSamples+1+postSamples);
    std::vector<HcalTrigPrimDigiCollection> hcalColl(preSamples+1+postSamples);
 
   unsigned nSamples = preSamples + postSamples + 1;
 
 
   if(iEvent.getByLabel(ecalTPG_,ecal))
     if (ecal.isValid()) 
       { 
       // loop through all ecal digis
       for (EcalTrigPrimDigiCollection::const_iterator ecal_it = ecal->begin(); ecal_it != ecal->end(); ecal_it++)
     {
       short zside             = ecal_it->id().zside();
       unsigned short ietaAbs  = ecal_it->id().ietaAbs();
       short iphi              = ecal_it->id().iphi();
       unsigned short digiSize = ecal_it->size();
       unsigned short nSOI     = (unsigned short) ( ecal_it->sampleOfInterest() );
       if (digiSize < nSamples || nSOI < preSamples || ((int)(digiSize - nSOI) < (int)(nSamples - preSamples)))
         {
           unsigned short preLoopsZero = (unsigned short) (preSamples)- nSOI;
 
           // fill extra bx's at beginning with zeros
           for (int sample = 0; sample < preLoopsZero; sample++)
         {
           // fill first few with zeros
           EcalTriggerPrimitiveDigi
             ecalDigi(EcalTrigTowerDetId((int) zside, EcalTriggerTower,
                         (int) ietaAbs, (int) iphi));
           ecalDigi.setSize(1);
           ecalDigi.setSample(0, EcalTriggerPrimitiveSample(0,false,0));
           ecalColl[sample].push_back(ecalDigi);
         }
 
           // loop through existing data
           for (int sample = preLoopsZero;sample < (preLoopsZero + digiSize); sample++)
         {
           // go through data
           EcalTriggerPrimitiveDigi
             ecalDigi(EcalTrigTowerDetId((int) zside, EcalTriggerTower,
                         (int) ietaAbs, (int) iphi));
           ecalDigi.setSize(1);
           if (useEcalCosmicTiming && iphi >= 1 && iphi <= 36)
             {
               if (nSOI < (preSamples + 1))
             {
               edm::LogWarning ("TooLittleData")
                 << "ECAL data needs at least one presample "
                 << "more than the number requested "
                 << "to use ecal cosmic timing mod!  "
                 << "reverting to useEcalCosmicTiming = false "
                 << "for rest of job.";
               useEcalCosmicTiming = false;
             }
               else
             {
               ecalDigi.setSample(0, EcalTriggerPrimitiveSample
                          (ecal_it->sample(nSOI + sample - 
                                   preSamples - 
                                   1).raw()));
             }
             }
           //else
           if ((!useEcalCosmicTiming) || (iphi >=37 && iphi <= 72))
             {
               ecalDigi.setSample(0, EcalTriggerPrimitiveSample
                      (ecal_it->sample(nSOI + sample - 
                               preSamples).raw()));
             }
           ecalColl[sample].push_back(ecalDigi);
         }
           
           // fill extra bx's at end with zeros
           for (unsigned int sample = (preLoopsZero + digiSize); 
            sample < nSamples; sample++)
         {
           // fill zeros!
           EcalTriggerPrimitiveDigi
             ecalDigi(EcalTrigTowerDetId((int) zside, EcalTriggerTower,
                         (int) ietaAbs, (int) iphi));
           ecalDigi.setSize(1);
           ecalDigi.setSample(0, EcalTriggerPrimitiveSample(0,false,0));
           ecalColl[sample].push_back(ecalDigi);
         }
         }
       else
         {
           for (unsigned short sample = 0; sample < nSamples; sample++)
         {
           // put each time sample into its own digi
           short zside = ecal_it->id().zside();
           unsigned short ietaAbs = ecal_it->id().ietaAbs();
           short iphi = ecal_it->id().iphi();
           EcalTriggerPrimitiveDigi
             ecalDigi(EcalTrigTowerDetId((int) zside, EcalTriggerTower,
                         (int) ietaAbs, (int) iphi));
           ecalDigi.setSize(1);
 
           if (useEcalCosmicTiming && iphi >= 1 && iphi <=36)
             {
               if (nSOI < (preSamples + 1))
             {
               edm::LogWarning ("TooLittleData")
                 << "ECAL data needs at least one presample "
                 << "more than the number requested "
                 << "to use ecal cosmic timing mod!  "
                 << "reverting to useEcalCosmicTiming = false "
                 << "for rest of job.";
               useEcalCosmicTiming = false;
             }
               else
             {
               ecalDigi.setSample(0, EcalTriggerPrimitiveSample
                          (ecal_it->sample
                           (ecal_it->sampleOfInterest() + 
                            sample - preSamples - 
                            1).raw()));
             }
             }
           //else
           if ((!useEcalCosmicTiming) || (iphi >=37 && iphi <= 72))
             {
               ecalDigi.setSample(0, EcalTriggerPrimitiveSample
                      (ecal_it->sample
                       (ecal_it->sampleOfInterest() + 
                        sample - preSamples).raw()));
             }
           // push back each digi into correct "time sample" of coll
           ecalColl[sample].push_back(ecalDigi);
         }
         }
     }
     }
 
   if(iEvent.getByLabel(hcalTPG_,hcal))
     if (hcal.isValid()) 
     { 
       // loop through all hcal digis
       for (HcalTrigPrimDigiCollection::const_iterator hcal_it = hcal->begin(); hcal_it != hcal->end(); hcal_it++)
     {
       short ieta = hcal_it->id().ieta();
       short iphi = hcal_it->id().iphi();
       // loop through time samples for each digi
       unsigned short digiSize = hcal_it->size();
       // (size of each digi must be no less than nSamples)
       unsigned short nSOI = (unsigned short) (hcal_it->presamples());
       if (digiSize < nSamples || nSOI < preSamples
           || ((int)(digiSize - nSOI) < (int)(nSamples - preSamples)))
         {
           unsigned short preLoopsZero = (unsigned short) (preSamples) 
         - nSOI;
           // fill extra bx's at beginning with zeros
           for (int sample = 0; sample < preLoopsZero; sample++)
         {
           // fill first few with zeros
           HcalTriggerPrimitiveDigi hcalDigi(HcalTrigTowerDetId((int) ieta, (int) iphi));
           hcalDigi.setSize(1);
           hcalDigi.setPresamples(0);
           hcalDigi.setSample(0, HcalTriggerPrimitiveSample(0,false,0,0));
           hcalColl[sample].push_back(hcalDigi);
         }
 
           // loop through existing data
           for (int sample = preLoopsZero; 
            sample < (preLoopsZero + digiSize); sample++)
         {
           // go through data
           HcalTriggerPrimitiveDigi hcalDigi(HcalTrigTowerDetId((int) ieta, (int) iphi));
           hcalDigi.setSize(1);
           hcalDigi.setPresamples(0);
 
           if (useHcalCosmicTiming && iphi >= 1 && iphi <= 36)
             {
               if (nSOI < (preSamples + 1))
             {
               edm::LogWarning ("TooLittleData")
                 << "HCAL data needs at least one presample "
                 << "more than the number requested "
                 << "to use hcal cosmic timing mod!  "
                 << "reverting to useHcalCosmicTiming = false "
                 << "for rest of job.";
               useHcalCosmicTiming = false;
             }
               else
             {
               hcalDigi.setSample(0, HcalTriggerPrimitiveSample
                          (hcal_it->sample(hcal_it->
                                   presamples() + 
                                   sample - 
                                   preSamples - 
                                   1).raw()));
             }
             }
           //else
           if ((!useHcalCosmicTiming) || (iphi >= 37 && iphi <= 72))
             {
 
                   hcalDigi.setSample(0, HcalTriggerPrimitiveSample
                              (hcal_it->sample(hcal_it->
                              presamples() + 
                               sample - 
                               preSamples).raw()));
 
             }
           hcalColl[sample].push_back(hcalDigi);
         }
           
           // fill extra bx's at end with zeros
           for (unsigned int sample = (preLoopsZero + digiSize); 
            sample < nSamples; sample++)
         {
           // fill zeros!
           HcalTriggerPrimitiveDigi
             hcalDigi(HcalTrigTowerDetId((int) ieta, (int) iphi));
           hcalDigi.setSize(1);
           hcalDigi.setPresamples(0);
           hcalDigi.setSample(0, HcalTriggerPrimitiveSample(0,false,0,0));
           hcalColl[sample].push_back(hcalDigi);
         }
         }
       else
         {
           for (unsigned short sample = 0; sample < nSamples; sample++)
         {
           // put each (relevant) time sample into its own digi
           HcalTriggerPrimitiveDigi hcalDigi(HcalTrigTowerDetId((int) ieta, (int) iphi));
           hcalDigi.setSize(1);
           hcalDigi.setPresamples(0);
 
           if (useHcalCosmicTiming && iphi >= 1 && iphi <= 36)
             {
               if (nSOI < (preSamples + 1))
             {
               edm::LogWarning ("TooLittleData")
                 << "HCAL data needs at least one presample "
                 << "more than the number requested "
                 << "to use hcal cosmic timing mod!  "
                 << "reverting to useHcalCosmicTiming = false "
                 << "for rest of job.";
               useHcalCosmicTiming = false;
             }
               else
             {
               hcalDigi.setSample(0, HcalTriggerPrimitiveSample
                          (hcal_it->sample(hcal_it->
                                   presamples() + 
                                   sample - 
                                   preSamples - 
                                   1).raw()));
             }
             }
           //else
           if ((!useHcalCosmicTiming) || (iphi >= 37 && iphi <= 72))
             {
               if(ieta>-29 && ieta<29) 
                 hcalDigi.setSample(0, HcalTriggerPrimitiveSample
                            (hcal_it->sample(hcal_it->
                                   presamples() + 
                                 sample - 
                                 preSamples+hbShift).raw()));
               if(ieta<=-29 || ieta>=29)
                 hcalDigi.setSample(0, HcalTriggerPrimitiveSample
                            (hcal_it->sample(hcal_it->
                                   presamples() + 
                                 sample - 
                                 preSamples+hfShift).raw()));
             }
           hcalColl[sample].push_back(hcalDigi);  
         }
         }
     }
     }
 
 
    //Now put the events back to file
 
   for(int i=0;i<preSamples;++i)
     {
       char ecal_label[200];
       char hcal_label[200];
 
 
       sprintf(ecal_label,"ECALBxminus%d",preSamples-i);
       sprintf(hcal_label,"HCALBxminus%d",preSamples-i);
 
       std::auto_ptr<EcalTrigPrimDigiCollection> ecalIn(new EcalTrigPrimDigiCollection); 
       std::auto_ptr<HcalTrigPrimDigiCollection> hcalIn(new HcalTrigPrimDigiCollection);
       for(unsigned int j=0;j<ecalColl[i].size();++j)
         {
           ecalIn->push_back((ecalColl[i])[j]);
         }
       for(unsigned int j=0;j<hcalColl[i].size();++j)
         hcalIn->push_back((hcalColl[i])[j]);
 
          iEvent.put(ecalIn,ecal_label);
          iEvent.put(hcalIn,hcal_label);
 
     }
 
 
    std::auto_ptr<EcalTrigPrimDigiCollection> ecal0(new EcalTrigPrimDigiCollection); 
    std::auto_ptr<HcalTrigPrimDigiCollection> hcal0(new HcalTrigPrimDigiCollection);
    for(unsigned int j=0;j<ecalColl[preSamples].size();++j)
      ecal0->push_back((ecalColl[preSamples])[j]);
    for(unsigned int j=0;j<hcalColl[preSamples].size();++j)
      hcal0->push_back((hcalColl[preSamples])[j]);
 
    iEvent.put(ecal0,"ECALBx0");
    iEvent.put(hcal0,"HCALBx0");
 
 
   for(int i=preSamples+1;i<preSamples+postSamples+1;++i)
     {
       char ecal_label[200];
       char hcal_label[200];
 
       sprintf(ecal_label,"ECALBxplus%d",i-preSamples);
       sprintf(hcal_label,"HCALBxplus%d",i-preSamples);
 
       std::auto_ptr<EcalTrigPrimDigiCollection> ecalIn2(new EcalTrigPrimDigiCollection); 
       std::auto_ptr<HcalTrigPrimDigiCollection> hcalIn2(new HcalTrigPrimDigiCollection);
 
       for(unsigned int j=0;j<ecalColl[i].size();++j)
     ecalIn2->push_back((ecalColl[i])[j]);
 
       for(unsigned int j=0;j<hcalColl[i].size();++j)
     hcalIn2->push_back((hcalColl[i])[j]);
 
 
          iEvent.put(ecalIn2,ecal_label);
          iEvent.put(hcalIn2,hcal_label);
     }
 }

Member Data Documentation

edm::InputTag L1RCTTPGProvider::ecalTPG_

private

Definition at line 46 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::hbShift

private

Definition at line 53 of file L1RCTTPGProvider.h.

Referenced by produce().

edm::InputTag L1RCTTPGProvider::hcalTPG_

private

Definition at line 47 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::hfShift

private

Definition at line 52 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::postSamples

private

Definition at line 51 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

int L1RCTTPGProvider::preSamples

private

Definition at line 50 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

bool L1RCTTPGProvider::useEcalCosmicTiming

private

Definition at line 49 of file L1RCTTPGProvider.h.

Referenced by produce().

bool L1RCTTPGProvider::useHcalCosmicTiming

private

Definition at line 48 of file L1RCTTPGProvider.h.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation