#include <TH1F.h>
#include <TROOT.h>
#include <TFile.h>
#include <TSystem.h>
#include "DataFormats/FWLite/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/FWLite/interface/FWLiteEnabler.h"
#include "DataFormats/FWLite/interface/InputSource.h"
#include "DataFormats/FWLite/interface/OutputFiles.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSetReader/interface/ParameterSetReader.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "PhysicsTools/FWLite/interface/TFileService.h"

Functions
int	main (int argc, char *argv[])

Function Documentation

int main	(	int	argc,
		char *	argv[]
	)

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Definition at line 19 of file FWLiteWithPythonConfig.cc.

References fwlite::Event::atEnd(), gather_cfg::cout, DeadROC_duringRun::dir, FWLiteEnabler::enable(), ev, beamvalidation::exit(), fwlite::OutputFiles::file(), edm::ParameterSet::getParameter(), TFileDirectory::make(), ResonanceBuilder::mass, TFileDirectory::mkdir(), dumpRecoGeometry_cfg::Muon, patZpeak::muons, LaserDQM_cfg::process, edm::readPSetsFrom(), TrackerOfflineValidation_Standalone_cff::TFileService, and fwlite::Event::toBegin().

                                  {
   // define what muon you are using; this is necessary as FWLite is not
   // capable of reading edm::Views
   using reco::Muon;
 
   // ----------------------------------------------------------------------
   // First Part:
   //
   //  * enable FWLite
   //  * book the histograms of interest
   //  * open the input file
   // ----------------------------------------------------------------------
 
   // load framework libraries
   gSystem->Load("libFWCoreFWLite");
   FWLiteEnabler::enable();
 
   // parse arguments
   if (argc < 2) {
     std::cout << "Usage : " << argv[0] << " [parameters.py]" << std::endl;
     return 0;
   }
 
   if (!edm::readPSetsFrom(argv[1])->existsAs<edm::ParameterSet>("process")) {
     std::cout << " ERROR: ParametersSet 'process' is missing in your configuration file" << std::endl;
     exit(0);
   }
   // get the python configuration
   const edm::ParameterSet& process = edm::readPSetsFrom(argv[1])->getParameter<edm::ParameterSet>("process");
   fwlite::InputSource inputHandler_(process);
   fwlite::OutputFiles outputHandler_(process);
 
   // now get each parameter
   const edm::ParameterSet& ana = process.getParameter<edm::ParameterSet>("muonAnalyzer");
   edm::InputTag muons_(ana.getParameter<edm::InputTag>("muons"));
 
   // book a set of histograms
   fwlite::TFileService fs = fwlite::TFileService(outputHandler_.file());
   TFileDirectory dir = fs.mkdir("analyzeBasicPat");
   TH1F* muonPt_ = dir.make<TH1F>("muonPt", "pt", 100, 0., 300.);
   TH1F* muonEta_ = dir.make<TH1F>("muonEta", "eta", 100, -3., 3.);
   TH1F* muonPhi_ = dir.make<TH1F>("muonPhi", "phi", 100, -5., 5.);
   TH1F* mumuMass_ = dir.make<TH1F>("mumuMass", "mass", 90, 30., 120.);
 
   // loop the events
   int ievt = 0;
   int maxEvents_(inputHandler_.maxEvents());
   for (unsigned int iFile = 0; iFile < inputHandler_.files().size(); ++iFile) {
     // open input file (can be located on castor)
     TFile* inFile = TFile::Open(inputHandler_.files()[iFile].c_str());
     if (inFile) {
       // ----------------------------------------------------------------------
       // Second Part:
       //
       //  * loop the events in the input file
       //  * receive the collections of interest via fwlite::Handle
       //  * fill the histograms
       //  * after the loop close the input file
       // ----------------------------------------------------------------------
       fwlite::Event ev(inFile);
       for (ev.toBegin(); !ev.atEnd(); ++ev, ++ievt) {
         edm::EventBase const& event = ev;
         // break loop if maximal number of events is reached
         if (maxEvents_ > 0 ? ievt + 1 > maxEvents_ : false)
           break;
         // simple event counter
         if (inputHandler_.reportAfter() != 0 ? (ievt > 0 && ievt % inputHandler_.reportAfter() == 0) : false)
           std::cout << "  processing event: " << ievt << std::endl;
 
         // Handle to the muon collection
         edm::Handle<std::vector<Muon> > muons;
         event.getByLabel(muons_, muons);
 
         // loop muon collection and fill histograms
         for (std::vector<Muon>::const_iterator mu1 = muons->begin(); mu1 != muons->end(); ++mu1) {
           muonPt_->Fill(mu1->pt());
           muonEta_->Fill(mu1->eta());
           muonPhi_->Fill(mu1->phi());
           if (mu1->pt() > 20 && fabs(mu1->eta()) < 2.1) {
             for (std::vector<Muon>::const_iterator mu2 = muons->begin(); mu2 != muons->end(); ++mu2) {
               if (mu2 > mu1) {                            // prevent double conting
                 if (mu1->charge() * mu2->charge() < 0) {  // check only muon pairs of unequal charge
                   if (mu2->pt() > 20 && fabs(mu2->eta()) < 2.1) {
                     mumuMass_->Fill((mu1->p4() + mu2->p4()).mass());
                   }
                 }
               }
             }
           }
         }
       }
       // close input file
       inFile->Close();
     }
     // break loop if maximal number of events is reached:
     // this has to be done twice to stop the file loop as well
     if (maxEvents_ > 0 ? ievt + 1 > maxEvents_ : false)
       break;
   }
   return 0;
 }

Functions

Function Documentation