Inheritance diagram for FactorizedJetCorrectorDemo:

Public Types
typedef reco::Particle::LorentzVector	LorentzVector

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

	~FactorizedJetCorrectorDemo ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

virtual void	endJob ()

Private Attributes
edm::Service< TFileService >	fs

bool	mDebug

double	mEtaMax

double	mEtaMin

TH2F *	mJECvsEta

TH2F *	mJECvsPt

std::string	mJetCorService

std::vector< std::string >	mLevels

int	mNGraphPoints

int	mNHistoPoints

std::string	mPayloadName

double	mPtMax

double	mPtMin

TRandom *	mRandom

TGraph *	mUncCorPt [100]

std::string	mUncertaintyFile

std::string	mUncertaintyTag

TGraph *	mUncEta [100]

bool	mUseCondDB

std::vector< double >	mVEta

TGraphErrors *	mVGraphCorPt [100]

TGraphErrors *	mVGraphEta [100]

TGraphErrors *	mVGraphPt [100]

std::vector< double >	mVPt

double	veta [100][1000]

double	vex_eta [100][1000]

double	vex_pt [100][1000]

double	vjec_eta [100][1000]

double	vjec_pt [100][1000]

double	vjecUnc_eta [100][1000]

double	vjecUnc_pt [100][1000]

double	vpt [100][1000]

double	vptcor [100][1000]

double	vUnc_eta [100][1000]

double	vUnc_pt [100][1000]

Additional Inherited Members
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::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 29 of file FactorizedJetCorrectorDemo.cc.

Member Typedef Documentation

typedef reco::Particle::LorentzVector FactorizedJetCorrectorDemo::LorentzVector

Definition at line 33 of file FactorizedJetCorrectorDemo.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 58 of file FactorizedJetCorrectorDemo.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mDebug, mEtaMax, mEtaMin, mLevels, mNGraphPoints, mNHistoPoints, mPayloadName, mPtMax, mPtMin, mUncertaintyFile, mUncertaintyTag, mVEta, and mVPt.

 {
   mLevels            = iConfig.getParameter<std::vector<std::string> > ("levels");
   mPayloadName       = iConfig.getParameter<std::string>          ("PayloadName");
   mUncertaintyTag    = iConfig.getParameter<std::string>          ("UncertaintyTag");
   mUncertaintyFile   = iConfig.getParameter<std::string>          ("UncertaintyFile");
   mNHistoPoints      = iConfig.getParameter<int>                  ("NHistoPoints");
   mNGraphPoints      = iConfig.getParameter<int>                  ("NGraphPoints");
   mEtaMin            = iConfig.getParameter<double>               ("EtaMin");
   mEtaMax            = iConfig.getParameter<double>               ("EtaMax");
   mPtMin             = iConfig.getParameter<double>               ("PtMin");
   mPtMax             = iConfig.getParameter<double>               ("PtMax");
   mVEta              = iConfig.getParameter<std::vector<double> > ("VEta");
   mVPt               = iConfig.getParameter<std::vector<double> > ("VPt");
   mDebug             = iConfig.getUntrackedParameter<bool>        ("Debug",false);
 }

FactorizedJetCorrectorDemo::~FactorizedJetCorrectorDemo ( )

Definition at line 75 of file FactorizedJetCorrectorDemo.cc.

 {
   
 }

Member Function Documentation

void FactorizedJetCorrectorDemo::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 80 of file FactorizedJetCorrectorDemo.cc.

References gather_cfg::cout, alignCSCRings::e, eta(), edm::EventSetup::get(), i, patTestJEC_cfi::jec, testEve_cfg::level, mDebug, mEtaMax, mEtaMin, mJECvsEta, mJECvsPt, mLevels, mNGraphPoints, mNHistoPoints, mPayloadName, mPtMax, mPtMin, mRandom, mVEta, mVPt, Parameters::parameters, funct::pow(), tmp, veta, vex_eta, vjec_eta, vjec_pt, vpt, and vptcor.

 {
 
   if ( mDebug ) 
     std::cout << "Hello from FactorizedJetCorrectorDemo" << std::endl;
   // retreive parameters from the DB this still need a proper configurable 
   // payloadName like: JetCorrectorParametersCollection_Spring10_AK5Calo.
   edm::ESHandle<JetCorrectorParametersCollection> parameters;
   iSetup.get<JetCorrectionsRecord>().get(mPayloadName, parameters); 
 
   std::vector<JetCorrectorParameters> params;
   for(std::vector<std::string>::const_iterator level=mLevels.begin(); level!=mLevels.end(); ++level){ 
     const JetCorrectorParameters& ip = (*parameters)[*level]; //ip.printScreen();
     if ( mDebug ) 
       std::cout << "Adding level " << *level << std::endl;
     params.push_back(ip); 
   } 
 
   boost::shared_ptr<FactorizedJetCorrector> corrector ( new FactorizedJetCorrector(params));
 
 
   double jec,rawPt,corPt,eta;
   TLorentzVector P4;
   double dEta = (mEtaMax-mEtaMin)/mNGraphPoints;
   if ( mDebug )
     std::cout << "Making JEC vs Eta and pT" << std::endl;
   for(int i=0;i<mNHistoPoints;i++)
     {
       rawPt  = mRandom->Uniform(mPtMin,mPtMax);
       eta = mRandom->Uniform(mEtaMin,mEtaMax);
       P4.SetPtEtaPhiE(rawPt,eta,0,0); 
       corrector->setJetEta( eta );
       corrector->setJetPt( rawPt );
       corrector->setJetE( P4.E() );
       jec = corrector->getCorrection();
       mJECvsEta->Fill(eta,jec);
       mJECvsPt->Fill(rawPt,jec);
     }
   if ( mDebug )
     std::cout << "Making JEC vs pT for different etas" << std::endl;
   //--------- Pt Graphs ------------------
   for(unsigned ieta=0;ieta<mVEta.size();ieta++)
     {
       double rPt  = pow((3500./TMath::CosH(mVEta[ieta]))/mPtMin,1./mNGraphPoints);
       for(int i=0;i<mNGraphPoints;i++) 
         {
           rawPt  = mPtMin*pow(rPt,i);
           eta = mVEta[ieta];
           vpt[ieta][i] = rawPt;
           P4.SetPtEtaPhiE(rawPt,eta,0,0); 
       corrector->setJetEta( eta );
       corrector->setJetPt( rawPt );
       corrector->setJetE( P4.E() );
       jec = corrector->getCorrection();
           vjec_eta[ieta][i] = jec;
           vptcor[ieta][i] = rawPt*jec;
           vex_eta[ieta][i] = 0.0;
           if (mDebug)
             std::cout<<rawPt<<" "<<eta<<" "<<jec<<" "<<rawPt*jec<<std::endl;
         }
     }
   if ( mDebug )
     std::cout << "Making JEC vs eta for different pTs" << std::endl;
   //--------- Eta Graphs -------------------
   for(unsigned ipt=0;ipt<mVPt.size();ipt++)
     {
       for(int i=0;i<mNGraphPoints;i++)
         {
           eta = mEtaMin + i*dEta;
           corPt  = mVPt[ipt];
           veta[ipt][i] = eta;
           //---------- find the raw pt -----------
           double e = 1.0;
           int nLoop(0);
           rawPt = corPt; 
           while(e > 0.0001 && nLoop < 10) 
              {
                P4.SetPtEtaPhiE(rawPt,eta,0,0); 
                LorentzVector rawP4(P4.Px(),P4.Py(),P4.Pz(),P4.E()); 
            corrector->setJetEta( eta );
            corrector->setJetPt( rawPt );
            corrector->setJetE( P4.E() );
            jec = corrector->getCorrection();
                double tmp = rawPt * jec;
                e = fabs(tmp-corPt)/corPt;
                if (jec > 0)
                  rawPt = corPt/jec;
                nLoop++;
              } 
           //--------- calculate the jec for the rawPt --------
           P4.SetPtEtaPhiE(rawPt,eta,0,0); 
           LorentzVector rawP4(P4.Px(),P4.Py(),P4.Pz(),P4.E());
       corrector->setJetEta( eta );
       corrector->setJetPt( rawPt );
       corrector->setJetE( P4.E() );
       jec = corrector->getCorrection();
           vjec_pt[ipt][i] = jec;
           if (mDebug)
             std::cout<<rawPt<<" "<<eta<<" "<<jec<<" "<<rawPt*jec<<std::endl;
         }
     }
   if ( mDebug )
     std::cout << "See ya!" << std::endl;
 
 }

void FactorizedJetCorrectorDemo::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 186 of file FactorizedJetCorrectorDemo.cc.

References edm::hlt::Exception, fs, TFileDirectory::make(), mEtaMax, mEtaMin, mJECvsEta, mJECvsPt, mNGraphPoints, mPtMax, mPtMin, mRandom, mVEta, and mVPt.

 {
   if (mNGraphPoints > 1000)
     throw  cms::Exception("FactorizedJetCorrectorDemo","Too many graph points !!! Maximum is 1000 !!!");
   if (mVEta.size() > 100)
     throw  cms::Exception("FactorizedJetCorrectorDemo","Too many eta values !!! Maximum is 100 !!!");
   if (mVPt.size() > 100)
     throw  cms::Exception("FactorizedJetCorrectorDemo","Too many pt values !!! Maximum is 100 !!!");
   mJECvsEta = fs->make<TH2F>("JECvsEta","JECvsEta",200,mEtaMin,mEtaMax,100,0,5);
   mJECvsPt  = fs->make<TH2F>("JECvsPt","JECvsPt",200,mPtMin,mPtMax,100,0,5);
   mRandom   = new TRandom();
   mRandom->SetSeed(0);
 }

void FactorizedJetCorrectorDemo::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 200 of file FactorizedJetCorrectorDemo.cc.

References fs, TFileDirectory::make(), mNGraphPoints, mUncCorPt, mUncEta, mVEta, mVGraphCorPt, mVGraphEta, mVGraphPt, mVPt, mergeVDriftHistosByStation::name, veta, vex_eta, vex_pt, vjec_eta, vjec_pt, vjecUnc_eta, vjecUnc_pt, vpt, vptcor, vUnc_eta, and vUnc_pt.

 {
   char name[1000];
   for(unsigned ipt=0;ipt<mVPt.size();ipt++)
     {
       mVGraphEta[ipt] = fs->make<TGraphErrors>(mNGraphPoints,veta[ipt],vjec_pt[ipt],vex_pt[ipt],vjecUnc_pt[ipt]);
       sprintf(name,"JEC_vs_Eta_CorPt%1.1f",mVPt[ipt]);
       mVGraphEta[ipt]->SetName(name);
       mUncEta[ipt] = fs->make<TGraph>(mNGraphPoints,veta[ipt],vUnc_pt[ipt]);
       sprintf(name,"UNC_vs_Eta_CorPt%1.1f",mVPt[ipt]);
       mUncEta[ipt]->SetName(name);
     }
   for(unsigned ieta=0;ieta<mVEta.size();ieta++)
     {
       mVGraphPt[ieta] = fs->make<TGraphErrors>(mNGraphPoints,vpt[ieta],vjec_eta[ieta],vex_eta[ieta],vjecUnc_eta[ieta]);
       sprintf(name,"JEC_vs_RawPt_eta%1.1f",mVEta[ieta]);
       mVGraphPt[ieta]->SetName(name);
       mVGraphCorPt[ieta] = fs->make<TGraphErrors>(mNGraphPoints,vptcor[ieta],vjec_eta[ieta],vex_eta[ieta],vjecUnc_eta[ieta]);
       sprintf(name,"JEC_vs_CorPt_eta%1.1f",mVEta[ieta]);
       mVGraphCorPt[ieta]->SetName(name);
       mUncCorPt[ieta] = fs->make<TGraph>(mNGraphPoints,vptcor[ieta],vUnc_eta[ieta]);
       sprintf(name,"UNC_vs_CorPt_eta%1.1f",mVEta[ieta]);
       mUncCorPt[ieta]->SetName(name);
     }
 }