#include <CastorTTRecord.h>

Inheritance diagram for CastorTTRecord:

Public Member Functions
	CastorTTRecord (const edm::ParameterSet &ps)

void	getEnergy_fC (double energy[16][14], edm::Handle< CastorDigiCollection > &CastorDigiColl, edm::Event &e, const edm::EventSetup &c)

void	getTriggerDecisions (std::vector< bool > &decision, double energy[16][14]) const

void	getTriggerDecisionsPerOctant (std::vector< bool > tdps[16], double energy[16][14]) const

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

virtual	~CastorTTRecord ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

	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 Attributes
edm::EDGetTokenT < CastorDigiCollection >	CastorDigiColl_

unsigned int	CastorSignalTS_

double	reweighted_gain

std::vector< std::string >	TrigNames_

std::vector< double >	TrigThresholds_

std::vector< unsigned int >	ttpBits_

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

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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::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

Definition at line 11 of file CastorTTRecord.h.

Constructor & Destructor Documentation

CastorTTRecord::CastorTTRecord ( const edm::ParameterSet & ps )

explicit

Definition at line 14 of file CastorTTRecord.cc.

References CastorDigiColl_, CastorSignalTS_, edm::ParameterSet::getParameter(), reweighted_gain, TrigNames_, TrigThresholds_, and ttpBits_.

 {
     CastorDigiColl_ = consumes<CastorDigiCollection>(ps.getParameter<edm::InputTag>("CastorDigiCollection")) ;
     CastorSignalTS_ = ps.getParameter< unsigned int >("CastorSignalTS") ;
 
     ttpBits_        = ps.getParameter< std::vector<unsigned int> >("ttpBits");
     TrigNames_      = ps.getParameter< std::vector<std::string> >("TriggerBitNames");
     TrigThresholds_ = ps.getParameter< std::vector<double> >("TriggerThresholds");
 
     reweighted_gain = 1.0;
 
     produces<L1GtTechnicalTriggerRecord>();
 }

CastorTTRecord::~CastorTTRecord ( )

virtual

Definition at line 29 of file CastorTTRecord.cc.

29 {

30 }

Member Function Documentation

void CastorTTRecord::getEnergy_fC	(	double	energy[16][14],
		edm::Handle< CastorDigiCollection > &	CastorDigiColl,
		edm::Event &	e,
		const edm::EventSetup &	c
	)

Definition at line 69 of file CastorTTRecord.cc.

References CastorCoderDb::adc2fC(), CastorSignalTS_, CastorCalibrations::gain(), edm::EventSetup::get(), CastorDataFrame::id(), HcalCastorDetId::module(), CastorCalibrations::pedestal(), reweighted_gain, and HcalCastorDetId::sector().

Referenced by produce().

 {
     // std::cerr << "**** RUNNING THROUGH CastorTTRecord::getEnergy_fC" << std::endl;
 
     // Get Conditions
     edm::ESHandle<CastorDbService> conditions ;
     eventSetup.get<CastorDbRecord>().get(conditions) ;
     const CastorQIEShape* shape = conditions->getCastorShape () ; // this one is generic
 
     for(int isec=0; isec<16; isec++) for(int imod=0; imod<14; imod++) energy[isec][imod] = 0;
 
     // Loop over digis
     CastorDigiCollection::const_iterator idigi ;
     for (idigi=CastorDigiColl->begin(); idigi!=CastorDigiColl->end(); idigi++) {
         const CastorDataFrame & digi = (*idigi) ;
         HcalCastorDetId cell = digi.id() ;
 
         // Get Castor Coder
         const CastorQIECoder* channelCoder = conditions->getCastorCoder(cell);
         CastorCoderDb coder (*channelCoder, *shape);
 
         // Get Castor Calibration
         const CastorCalibrations& calibrations=conditions->getCastorCalibrations(cell);
 
         // convert adc to fC
         CaloSamples tool ;
         coder.adc2fC(digi,tool) ;
 
         // pedestal substraction
         int capid=digi[CastorSignalTS_].capid();
         double fC = tool[CastorSignalTS_] - calibrations.pedestal(capid);
 
         // to correct threshold levels in fC for different gains
         reweighted_gain = calibrations.gain(capid) / 0.015;
         
         energy[digi.id().sector()-1][digi.id().module()-1] = fC;
     }
 }

void CastorTTRecord::getTriggerDecisions	(	std::vector< bool > &	decision,
		double	energy[16][14]
	)		const

Definition at line 109 of file CastorTTRecord.cc.

References getTriggerDecisionsPerOctant().

Referenced by produce().

 {
     // std::cerr << "**** RUNNING THROUGH CastorTTRecord::getTriggerDecisions" << std::endl;
 
     // check if number of bits is at least four
     if( decision.size() < 4 ) return;
 
     std::vector<bool> tdpo[8]; // TriggerDecisionsPerOctant
     getTriggerDecisionsPerOctant(tdpo,energy);
 
 
     // preset trigger decisions
     decision.at(0) = true;
     decision.at(1) = false;
     decision.at(2) = false;
     decision.at(3) = false;
 
     // canceld for low pt jet
     // bool EM_decision = false;
     // bool HAD_decision = false;
     // loop over castor octants
     for(int ioct=0; ioct<8; ioct++) {
         int next_oct = (ioct+1)%8;
         int prev_oct = (ioct+8-1)%8;
 
         // gap Trigger
         if( !tdpo[ioct].at(0) ) decision.at(0) = false;
         if( !tdpo[ioct].at(1) ) decision.at(0) = false;
 
         // jet Trigger
         if( tdpo[ioct].at(2) ) decision.at(1) = true;
 
         // electron
         // canceld for low pt jet
         // if( tdpo[ioct].at(3) ) EM_decision = true;
         // if( tdpo[ioct].at(4) ) HAD_decision = true;
 
         // iso muon
         if( tdpo[ioct].at(5) ) {
             // was one of the other sectors 
             // in the octant empty ?
             if( tdpo[ioct].at(0) ) {
                 if( tdpo[prev_oct].at(1) && 
                     tdpo[next_oct].at(0) && 
                     tdpo[next_oct].at(1) )
                     decision.at(3) = true;
             }
             else if( tdpo[ioct].at(1) ) {
                 if( tdpo[prev_oct].at(0) && 
                     tdpo[prev_oct].at(1) && 
                     tdpo[next_oct].at(0) )
                     decision.at(3) = true;
             }
             // when not no iso muon
         }
 
         // low pt jet Trigger
         if( tdpo[ioct].at(6) ) decision.at(2) = true;
     }
 
     // for EM Trigger whole castor not hadronic and somewhere EM
     // canceld for low pt jet
     // decision.at(2) = EM_decision && !HAD_decision;
 }

void CastorTTRecord::getTriggerDecisionsPerOctant	(	std::vector< bool >	tdps[16],
		double	energy[16][14]
	)		const

Definition at line 174 of file CastorTTRecord.cc.

References reweighted_gain, and TrigThresholds_.

Referenced by getTriggerDecisions().

 {
     // std::cerr << "**** RUNNING THROUGH CastorTTRecord::getTriggerDecisionsPerOctant" << std::endl;
 
     // loop over octatants
     for(int ioct=0; ioct<8; ioct++)
     {
         // six bits from HTR card
         // 0. first sector empty
         // 1. second sector empty
         // 2. jet any sector
         // 3. EM any sector
         // 4. HAD any sector
         // 5. muon any sector
         // add instead of EM Trigger (not bit 6 in real)
         // 6. low pt jet any sector
         tdpo[ioct].resize(7);
 
         for(int ibit=0; ibit<7; ibit++)
             tdpo[ioct].at(ibit) = false;
 
         // loop over castor sectors in octant
         for(int ioctsec=0; ioctsec<2; ioctsec++)
         {
             // absolute sector number
             int isec = 2*ioct+ioctsec;
 
             // init module sums for every sector
             double fCsum_mod = 0;
             double fCsum_em = 0, fCsum_ha = 0;
             double fCsum_jet_had = 0;
             double fCsum_col[3] = { 0, 0, 0 };
             
             // loop over modules
             for(int imod=0; imod<14; imod++) { 
                 // total sum
                 fCsum_mod += energy[isec][imod];
                 
                 // EM & HAD sum
                 if( imod < 2 ) fCsum_em += energy[isec][imod];
                 if( imod > 2 && imod < 12 ) fCsum_ha += energy[isec][imod];
                 
                 // sum over three sector parts
                 if( imod < 4 )       fCsum_col[0] += energy[isec][imod];
                 else if( imod < 8 )  fCsum_col[1] += energy[isec][imod];
                 else if( imod < 12 ) fCsum_col[2] += energy[isec][imod];
 
                 // HAD sum for jet trigger v2
                 if( imod > 1 && imod < 5 ) fCsum_jet_had += energy[isec][imod];
             }
 
             // gap Trigger
             if( fCsum_mod < TrigThresholds_.at(0) ) {
                 if( ioctsec == 0 ) tdpo[ioct].at(0) = true;
                 else if( ioctsec == 1 ) tdpo[ioct].at(1) = true;
             }
 
             // jet Trigger 
             // with gain correction
             /* old version of jet trigger ( deprecated because of saturation )
             if( fCsum_mod > TrigThresholds_.at(1) / reweighted_gain ) 
                 tdpo[ioct].at(2) = true; 
             */
             if( fCsum_jet_had > TrigThresholds_.at(1) / reweighted_gain )
             // additional high threshold near saturation for EM part
                 if( energy[isec][0] > 26000 / reweighted_gain && energy[isec][1] > 26000 / reweighted_gain )
                     tdpo[ioct].at(2) = true;
 
             // low pt jet Trigger
             if( fCsum_mod > TrigThresholds_.at(5) / reweighted_gain ) 
                 tdpo[ioct].at(6) = true; 
 
             // egamma Trigger
             // with gain correction only in the EM threshold
             if( fCsum_em > TrigThresholds_.at(2) / reweighted_gain )
                 tdpo[ioct].at(3) = true;
             if( fCsum_ha > TrigThresholds_.at(3) ) 
                 tdpo[ioct].at(4) = true; 
 
             // muon Trigger
             int countColumns = 0;
             for( int icol=0; icol<3; icol++ ) 
                 if( fCsum_col[icol] > TrigThresholds_.at(4) ) 
                     countColumns++;
             if( countColumns >= 2 )
                 tdpo[ioct].at(5) = true;
         }
     }
 }

void CastorTTRecord::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

Implements edm::EDProducer.

Definition at line 32 of file CastorTTRecord.cc.

References CastorDigiColl_, edm::Event::getByToken(), getEnergy_fC(), getTriggerDecisions(), i, convertSQLitetoXML_cfg::output, edm::Event::put(), TrigNames_, and ttpBits_.

                                                                          {
 
     // std::cerr << "**** RUNNING THROUGH CastorTTRecord::produce" << std::endl;
 
     std::vector<L1GtTechnicalTrigger> vecTT(ttpBits_.size()) ;
 
     // Get Inputs
     edm::Handle<CastorDigiCollection> CastorDigiColl ; 
     e.getByToken(CastorDigiColl_,CastorDigiColl) ;
 
     if ( !CastorDigiColl.failedToGet() ) { 
 
         double cas_efC[16][14];
         getEnergy_fC(cas_efC,CastorDigiColl,e,eventSetup);
 
 
         std::vector<bool> decision(ttpBits_.size());
 
         getTriggerDecisions(decision,cas_efC);
 
         for(unsigned int i=0; i<ttpBits_.size(); i++) {
             // if( decision.at(i) ) std::cerr << "**** Something Triggered" << std::endl;
             // std::cout << "Run CastorTTRecord::produce. TriggerBit = " << ttpBits_.at(i) << "; TriggerName = " << TrigNames_.at(i) << "; Decision = " << decision[i] << std::endl;
             vecTT.at(i) = L1GtTechnicalTrigger(TrigNames_.at(i), ttpBits_.at(i), 0, decision.at(i)) ;
         }
 
     } else {
         vecTT.clear() ;
     }
 
     // Put output into event
     std::auto_ptr<L1GtTechnicalTriggerRecord> output(new L1GtTechnicalTriggerRecord()) ;
     output->setGtTechnicalTrigger(vecTT) ;    
     e.put(output) ;
 }