|
|
|
#include <EcalTrivialObjectAnalyzer.h>
Public Member Functions | |
| virtual void | analyze (const edm::Event &e, const edm::EventSetup &c) |
| EcalTrivialObjectAnalyzer (int i) | |
| EcalTrivialObjectAnalyzer (edm::ParameterSet const &p) | |
| virtual | ~EcalTrivialObjectAnalyzer () |
Definition at line 11 of file EcalTrivialObjectAnalyzer.h.
| EcalTrivialObjectAnalyzer::EcalTrivialObjectAnalyzer | ( | edm::ParameterSet const & | p | ) | [inline, explicit] |
Definition at line 14 of file EcalTrivialObjectAnalyzer.h.
{ }
| EcalTrivialObjectAnalyzer::EcalTrivialObjectAnalyzer | ( | int | i | ) | [inline, explicit] |
Definition at line 16 of file EcalTrivialObjectAnalyzer.h.
{ }
| virtual EcalTrivialObjectAnalyzer::~EcalTrivialObjectAnalyzer | ( | ) | [inline, virtual] |
Definition at line 18 of file EcalTrivialObjectAnalyzer.h.
{ }
| void EcalTrivialObjectAnalyzer::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 73 of file EcalTrivialObjectAnalyzer.cc.
References gather_cfg::cout, edm::EventID::event(), EcalCondObjectContainer< T >::find(), edm::EventSetup::get(), EcalLaserAPDPNRatios::getLaserMap(), EcalCondObjectContainer< T >::getMap(), EcalChannelStatusCode::getStatusCode(), i, edm::EventBase::id(), edm::ESHandle< T >::product(), DetId::rawId(), and edm::EventID::run().
{
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 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";
// 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 finihed." << 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
/*
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()
1.7.3