L1RCTTPGProvider Class Reference

#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 ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

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
typedef WorkerT< EDProducer >	WorkerType
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::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

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

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

Definition at line 37 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.

                         {
}

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, and useHcalCosmicTiming.

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

Referenced by produce().

int L1RCTTPGProvider::hbShift

private

Definition at line 54 of file L1RCTTPGProvider.h.

Referenced by produce().

edm::InputTag L1RCTTPGProvider::hcalTPG_

private

Definition at line 48 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::hfShift

private

Definition at line 53 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::postSamples

private

Definition at line 52 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

int L1RCTTPGProvider::preSamples

private

Definition at line 51 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

bool L1RCTTPGProvider::useEcalCosmicTiming

private

Definition at line 50 of file L1RCTTPGProvider.h.

Referenced by produce().

bool L1RCTTPGProvider::useHcalCosmicTiming

private

Definition at line 49 of file L1RCTTPGProvider.h.

Referenced by produce().