#include <EcalTrivialObjectAnalyzer.h>

Inheritance diagram for EcalTrivialObjectAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &e, const edm::EventSetup &c)

	EcalTrivialObjectAnalyzer (edm::ParameterSet const &p)

	EcalTrivialObjectAnalyzer (int i)

virtual	~EcalTrivialObjectAnalyzer ()

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

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

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

virtual	~EDConsumerBase () noexcept(false)

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

Definition at line 10 of file EcalTrivialObjectAnalyzer.h.

Constructor & Destructor Documentation

EcalTrivialObjectAnalyzer::EcalTrivialObjectAnalyzer ( edm::ParameterSet const & p )

inlineexplicit

Definition at line 13 of file EcalTrivialObjectAnalyzer.h.

14 { }

EcalTrivialObjectAnalyzer::EcalTrivialObjectAnalyzer ( int i )

inlineexplicit

Definition at line 15 of file EcalTrivialObjectAnalyzer.h.

16 { }

virtual EcalTrivialObjectAnalyzer::~EcalTrivialObjectAnalyzer ( )

inlinevirtual

Definition at line 17 of file EcalTrivialObjectAnalyzer.h.

References analyze(), EnergyCorrector::c, and MillePedeFileConverter_cfg::e.

17 { }

Member Function Documentation

void EcalTrivialObjectAnalyzer::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

Definition at line 80 of file EcalTrivialObjectAnalyzer.cc.

References gather_cfg::cout, EcalCondObjectContainer< T >::end(), edm::EventID::event(), EcalCondObjectContainer< T >::find(), EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), edm::EventSetup::get(), EcalTimeOffsetConstant::getEBValue(), EcalADCToGeVConstant::getEBValue(), EcalSampleMask::getEcalSampleMaskRecordEB(), EcalSampleMask::getEcalSampleMaskRecordEE(), EcalADCToGeVConstant::getEEValue(), EcalTimeOffsetConstant::getEEValue(), EcalLaserAPDPNRatios::getLaserMap(), EcalTBWeights::getMap(), EcalCondObjectContainer< T >::getMap(), EcalChannelStatusCode::getStatusCode(), EcalLaserAPDPNRatios::getTimeMap(), mps_fire::i, EcalXtalGroupId::id(), edm::EventBase::id(), createfilelist::int, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p1, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p2, EcalFunParams::params(), edm::ESHandle< T >::product(), diffTwoXMLs::r1, rand(), DetId::rawId(), edm::EventID::run(), EcalLaserAPDPNRatios::EcalLaserTimeStamp::t1, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t2, and edm::Timestamp::value().

Referenced by ~EcalTrivialObjectAnalyzer().

 {
     using namespace edm::eventsetup;
     // Context is not used.
     std::cout <<">>> EcalTrivialObjectAnalyzer: processing run "<<e.id().run() << " event: " << e.id().event() << std::endl;
 
     edm::ESHandle<EcalPedestals> pPeds;
     context.get<EcalPedestalsRcd>().get(pPeds);
 
     // ADC -> GeV Scale
     edm::ESHandle<EcalADCToGeVConstant> pAgc;
     context.get<EcalADCToGeVConstantRcd>().get(pAgc);
     const EcalADCToGeVConstant* agc = pAgc.product();
     std::cout << "Global ADC->GeV scale: EB " << std::setprecision(6) << agc->getEBValue() << " GeV/ADC count" 
           " EE " << std::setprecision(6) << agc->getEEValue() << " GeV/ADC count" << std::endl;
 
     // use a channel to fetch values from DB
     double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
     int ieta =  int( 1 + r1*85 );
     r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
     int iphi =  int( 1 + r1*20 );
 
     EBDetId ebid(ieta,iphi); //eta,phi
     std::cout << "EcalTrivialObjectAnalyzer: using EBDetId: " << ebid << std::endl;
 
     const EcalPedestals* myped=pPeds.product();
     EcalPedestals::const_iterator it = myped->find(ebid.rawId());
     if( it!=myped->end() ){
       std::cout << "EcalPedestal: "
                 << "  mean_x1:  " << std::setprecision(8) << (*it).mean_x1 << " rms_x1: " << (*it).rms_x1
                 << "  mean_x6:  " <<(*it).mean_x6 << " rms_x6: " << (*it).rms_x6
                 << "  mean_x12: " <<(*it).mean_x12 << " rms_x12: " << (*it).rms_x12
                 << std::endl;
     } else {
      std::cout << "No pedestal found for this xtal! something wrong with EcalPedestals in your DB? "
                << std::endl;
     }
 
     // fetch map of groups of xtals
     edm::ESHandle<EcalWeightXtalGroups> pGrp;
     context.get<EcalWeightXtalGroupsRcd>().get(pGrp);
     const EcalWeightXtalGroups* grp = pGrp.product();
 
     EcalXtalGroupsMap::const_iterator git = grp->getMap().find( ebid.rawId() );
     EcalXtalGroupId gid;
     if( git != grp->getMap().end() ) {
       std::cout << "XtalGroupId.id() = " << std::setprecision(3) << (*git).id() << std:: endl;
       gid = (*git);
     } else {
       std::cout << "No group id found for this crystal. something wrong with EcalWeightXtalGroups in your DB?"
                 << std::endl;
    }
 
     // Gain Ratios
     edm::ESHandle<EcalGainRatios> pRatio;
     context.get<EcalGainRatiosRcd>().get(pRatio);
     const EcalGainRatios* gr = pRatio.product();
 
     EcalGainRatioMap::const_iterator grit=gr->getMap().find(ebid.rawId());
     EcalMGPAGainRatio mgpa;
     if( grit!=gr->getMap().end() ){
       mgpa = (*grit);
 
       std::cout << "EcalMGPAGainRatio: "
                 << "gain 12/6 :  " << std::setprecision(4) << mgpa.gain12Over6() << " gain 6/1: " << mgpa.gain6Over1()
                 << std::endl;
     } else {
      std::cout << "No MGPA Gain Ratio found for this xtal! something wrong with EcalGainRatios in your DB? "
                << std::endl;
     }
 
     // Intercalib constants
     edm::ESHandle<EcalIntercalibConstants> pIcal;
     context.get<EcalIntercalibConstantsRcd>().get(pIcal);
     const EcalIntercalibConstants* ical = pIcal.product();
 
     EcalIntercalibConstantMap::const_iterator icalit=ical->getMap().find(ebid.rawId());
     EcalIntercalibConstant icalconst;
     if( icalit!=ical->getMap().end() ){
       icalconst = (*icalit);
 
       std::cout << "EcalIntercalibConstant: "
                 <<std::setprecision(6)
                 << icalconst
                 << std::endl;
     } else {
      std::cout << "No intercalib const found for this xtal! something wrong with EcalIntercalibConstants in your DB? "
                << std::endl;
     }
 
     // Intercalib errors
     edm::ESHandle<EcalIntercalibErrors> pIcalErr;
     context.get<EcalIntercalibErrorsRcd>().get(pIcalErr);
     const EcalIntercalibErrors* icalErr = pIcalErr.product();
 
     EcalIntercalibErrorMap::const_iterator icalitErr=icalErr->getMap().find(ebid.rawId());
     EcalIntercalibError icalconstErr;
     if( icalitErr!=icalErr->getMap().end() ){
       icalconstErr = (*icalitErr);
 
       std::cout << "EcalIntercalibError: "
                 <<std::setprecision(6)
                 << icalconstErr
                 << std::endl;
     } else {
      std::cout << "No intercalib const found for this xtal! something wrong with EcalIntercalibErrors in your DB? "
                << std::endl;
     }
 
 
 
     { // quick and dirty for cut and paste ;) it is a test program...
             // TimeCalib constants
             edm::ESHandle<EcalTimeCalibConstants> pIcal;
             context.get<EcalTimeCalibConstantsRcd>().get(pIcal);
             const EcalTimeCalibConstants* ical = pIcal.product();
 
             EcalTimeCalibConstantMap::const_iterator icalit=ical->getMap().find(ebid.rawId());
             EcalTimeCalibConstant icalconst;
             if( icalit!=ical->getMap().end() ){
                     icalconst = (*icalit);
 
                     std::cout << "EcalTimeCalibConstant: "
                             <<std::setprecision(6)
                             << icalconst
                             << std::endl;
             } else {
                     std::cout << "No intercalib const found for this xtal! something wrong with EcalTimeCalibConstants in your DB? "
                             << std::endl;
             }
 
             // TimeCalib errors
             edm::ESHandle<EcalTimeCalibErrors> pIcalErr;
             context.get<EcalTimeCalibErrorsRcd>().get(pIcalErr);
             const EcalTimeCalibErrors* icalErr = pIcalErr.product();
 
             EcalTimeCalibErrorMap::const_iterator icalitErr=icalErr->getMap().find(ebid.rawId());
             EcalTimeCalibError icalconstErr;
             if( icalitErr!=icalErr->getMap().end() ){
                     icalconstErr = (*icalitErr);
 
                     std::cout << "EcalTimeCalibError: "
                             <<std::setprecision(6)
                             << icalconstErr
                             << std::endl;
             } else {
                     std::cout << "No intercalib const found for this xtal! something wrong with EcalTimeCalibErrors in your DB? "
                             << std::endl;
             }
     }
 
    // fetch Time Offset
    //std::cout <<"Fetching TimeOffsetConstant from DB " << std::endl;
 
    // Time Offset constants
    edm::ESHandle<EcalTimeOffsetConstant> pTOff;
    context.get<EcalTimeOffsetConstantRcd>().get(pTOff);
    const EcalTimeOffsetConstant* TOff = pTOff.product();
 
    std::cout << "EcalTimeOffsetConstant: "
          << " EB " << TOff->getEBValue()
          << " EE " << TOff->getEEValue()
          << std::endl;
 
    // fetch TB weights
    std::cout <<"Fetching EcalTBWeights from DB " << std::endl;
    edm::ESHandle<EcalTBWeights> pWgts;
    context.get<EcalTBWeightsRcd>().get(pWgts);
    const EcalTBWeights* wgts = pWgts.product();
    std::cout << "EcalTBWeightMap.size(): " << std::setprecision(3) << wgts->getMap().size() << std::endl;
 
 
    // look up the correct weights for this  xtal
    //EcalXtalGroupId gid( (*git) );
    EcalTBWeights::EcalTDCId tdcid(1);
 
    std::cout << "Lookup EcalWeightSet for groupid: " << std::setprecision(3) 
              << gid.id() << " and TDC id " << tdcid << std::endl;
    EcalTBWeights::EcalTBWeightMap::const_iterator wit = wgts->getMap().find( std::make_pair(gid,tdcid) );
    EcalWeightSet  wset;
    if( wit != wgts->getMap().end() ) {
       wset = wit->second;
       std::cout << "check size of data members in EcalWeightSet" << std::endl;
       //wit->second.print(std::cout);
 
 
       //typedef std::vector< std::vector<EcalWeight> > EcalWeightMatrix;
       const EcalWeightSet::EcalWeightMatrix& mat1 = wit->second.getWeightsBeforeGainSwitch();
       const EcalWeightSet::EcalWeightMatrix& mat2 = wit->second.getWeightsAfterGainSwitch();
 
 //       std::cout << "WeightsBeforeGainSwitch.size: " << mat1.size()
 //                 << ", WeightsAfterGainSwitch.size: " << mat2.size() << std::endl;
 
 
       CLHEP::HepMatrix clmat1(3,10,0);
       CLHEP::HepMatrix clmat2(3,10,0);
       for(int irow=0; irow<3; irow++) {
        for(int icol=0; icol<10; icol++) {
          clmat1[irow][icol] = mat1(irow,icol);
          clmat2[irow][icol] = mat2(irow,icol);
        }
      }
      std::cout << "weight matrix before gain switch:" << std::endl;
      std::cout << clmat1 << std::endl;
      std::cout << "weight matrix after gain switch:" << std::endl;
      std::cout << clmat2 << std::endl;
 
    } else {
      std::cout << "No weights found for EcalGroupId: " << gid.id() << " and  EcalTDCId: " << tdcid << std::endl;
   }
 
    // cluster functions/corrections
    edm::ESHandle<EcalClusterLocalContCorrParameters> pLocalCont;
    context.get<EcalClusterLocalContCorrParametersRcd>().get(pLocalCont);
    const EcalClusterLocalContCorrParameters* paramLocalCont = pLocalCont.product();
    std::cout << "LocalContCorrParameters:";
    for ( EcalFunctionParameters::const_iterator it = paramLocalCont->params().begin(); it != paramLocalCont->params().end(); ++it ) {
            std::cout << " " << *it;
    }
    std::cout << "\n";
    edm::ESHandle<EcalClusterCrackCorrParameters> pCrack;
    context.get<EcalClusterCrackCorrParametersRcd>().get(pCrack);
    const EcalClusterCrackCorrParameters* paramCrack = pCrack.product();
    std::cout << "CrackCorrParameters:";
    for ( EcalFunctionParameters::const_iterator it = paramCrack->params().begin(); it != paramCrack->params().end(); ++it ) {
            std::cout << " " << *it;
    }
    std::cout << "\n";
    edm::ESHandle<EcalClusterEnergyCorrectionParameters> pEnergyCorrection;
    context.get<EcalClusterEnergyCorrectionParametersRcd>().get(pEnergyCorrection);
    const EcalClusterEnergyCorrectionParameters* paramEnergyCorrection = pEnergyCorrection.product();
    std::cout << "EnergyCorrectionParameters:";
    for ( EcalFunctionParameters::const_iterator it = paramEnergyCorrection->params().begin(); it != paramEnergyCorrection->params().end(); ++it ) {
            std::cout << " " << *it;
    }
    std::cout << "\n";
    edm::ESHandle<EcalClusterEnergyUncertaintyParameters> pEnergyUncertainty;
    context.get<EcalClusterEnergyUncertaintyParametersRcd>().get(pEnergyUncertainty);
    const EcalClusterEnergyUncertaintyParameters* paramEnergyUncertainty = pEnergyUncertainty.product();
    std::cout << "EnergyCorrectionParameters:";
    for ( EcalFunctionParameters::const_iterator it = paramEnergyUncertainty->params().begin(); it != paramEnergyUncertainty->params().end(); ++it ) {
            std::cout << " " << *it;
    }
    std::cout << "\n";
    edm::ESHandle<EcalClusterEnergyCorrectionObjectSpecificParameters> pEnergyCorrectionObjectSpecific;
    context.get<EcalClusterEnergyCorrectionObjectSpecificParametersRcd>().get(pEnergyCorrectionObjectSpecific);
    const EcalClusterEnergyCorrectionObjectSpecificParameters* paramEnergyCorrectionObjectSpecific = pEnergyCorrectionObjectSpecific.product();
    std::cout << "EnergyCorrectionObjectSpecificParameters:";
    for ( EcalFunctionParameters::const_iterator it = paramEnergyCorrectionObjectSpecific->params().begin(); it != paramEnergyCorrectionObjectSpecific->params().end(); ++it ) {
      std::cout << " " << *it;
    }
    std::cout << "\n";
 
 
    // laser correction
    
    // laser alphas
    edm::ESHandle<EcalLaserAlphas> pAlpha; 
    context.get<EcalLaserAlphasRcd>().get(pAlpha);
    const EcalLaserAlphas* lalpha = pAlpha.product();
    
    EcalLaserAlphaMap::const_iterator lalphait;
    lalphait = lalpha->getMap().find(ebid.rawId());
    if( lalphait!=lalpha->getMap().end() ){
      std::cout << "EcalLaserAlpha: "
            <<std::setprecision(6)
            << (*lalphait)
            << std::endl;
    } else {
      std::cout << "No laser alpha found for this xtal! something wrong with EcalLaserAlphas in your DB? "
                << std::endl;
    }
 
    // laser apdpnref
    edm::ESHandle<EcalLaserAPDPNRatiosRef> pAPDPNRatiosRef;          
    context.get<EcalLaserAPDPNRatiosRefRcd>().get(pAPDPNRatiosRef);  
    const EcalLaserAPDPNRatiosRef* lref = pAPDPNRatiosRef.product(); 
    
    EcalLaserAPDPNRatiosRef::const_iterator lrefit;
    lrefit = lref->getMap().find(ebid.rawId());
    if( lrefit!=lref->getMap().end() ){
      std::cout << "EcalLaserAPDPNRatiosRef: "
            <<std::setprecision(6)
            << (*lrefit)
            << std::endl;
    } else {
      std::cout << "No laser apd/pn ref found for this xtal! something wrong with EcalLaserAPDPNRatiosRef in your DB? "
            << std::endl;
    }
 
    // laser apdpn ratios 
    edm::ESHandle<EcalLaserAPDPNRatios> pAPDPNRatios;          
    context.get<EcalLaserAPDPNRatiosRcd>().get(pAPDPNRatios);  
    const EcalLaserAPDPNRatios* lratio = pAPDPNRatios.product();
 
    EcalLaserAPDPNRatios::EcalLaserAPDPNRatiosMap::const_iterator lratioit; 
    lratioit = lratio->getLaserMap().find(ebid.rawId());
    EcalLaserAPDPNRatios::EcalLaserAPDPNpair lratioconst;
    
    if( lratioit!=lratio->getLaserMap().end() ){
      lratioconst = (*lratioit);
      
      std::cout << "EcalLaserAPDPNRatios: "
        //          <<e.id().run() << " " << e.id().event() << " " 
                << std::setprecision(6)
                << lratioconst.p1 << " " << lratioconst.p2
                << std::endl;
    } else {
      std::cout << "No laser apd/pn ratio found for this xtal! something wrong with EcalLaserAPDPNRatios in your DB? "
            << std::endl;
    }
 
    // laser timestamps
    EcalLaserAPDPNRatios::EcalLaserTimeStamp ltimestamp;
    EcalLaserAPDPNRatios::EcalLaserTimeStampMap::const_iterator ltimeit;
    for (int i=1; i<=92; i++) {     
      ltimestamp = lratio->getTimeMap()[i];
      std::cout << "i = " << std::setprecision(6) << i
              << ltimestamp.t1.value() << " " << ltimestamp.t2.value() << " : " ;
    }
    std::cout << "Tests finished." << std::endl;
 
    // channel status
    edm::ESHandle<EcalChannelStatus> pChannelStatus;
    context.get<EcalChannelStatusRcd>().get(pChannelStatus);
    const EcalChannelStatus *ch_status = pChannelStatus.product();
 
    EcalChannelStatusMap::const_iterator chit;
    chit = ch_status->getMap().find(ebid.rawId());
    if( chit != ch_status->getMap().end() ){
            EcalChannelStatusCode ch_code = (*chit);
            std::cout << "EcalChannelStatus: "
                    <<std::setprecision(6)
                    << ch_code.getStatusCode()
                    << std::endl;
    } else {
            std::cout << "No channel status found for this xtal! something wrong with EcalChannelStatus in your DB? "
                    << std::endl;
    }
 
 
    // laser transparency correction
 
    // Mask to ignore sample
    edm::ESHandle<EcalSampleMask> pSMask;
    context.get<EcalSampleMaskRcd>().get(pSMask);
    const EcalSampleMask* smask = pSMask.product();
    std::cout << "Sample Mask EB " << std::hex << smask->getEcalSampleMaskRecordEB() 
          << " EE " << std::hex << smask->getEcalSampleMaskRecordEE() << std::endl;
    
 
 /*
     std::cout << "make CLHEP matrices from vector<vector<Ecalweight>>" << std::endl;
     CLHEP::HepMatrix clmat1(3,8,0);
     CLHEP::HepMatrix clmat2(3,8,0);
     for(int irow=0; irow<3; irow++) {
      for(int icol=0; icol<8; icol++) {
        clmat1[irow][icol] = (mat1[irow])[icol]();
        clmat2[irow][icol] = (mat2[irow])[icol]();
      }
    }
    std::cout << clmat1 << std::endl;
    std::cout << clmat2 << std::endl;
 */
 
 } //end of ::Analyze()

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation