#include <memory>
#include <string>
#include <vector>
#include <sstream>
#include <fstream>
#include <iostream>
#include <TH1F.h>
#include <TROOT.h>
#include <TFile.h>
#include <TSystem.h>
#include "DataFormats/FWLite/interface/Handle.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "DataFormats/PatCandidates/interface/CompositeCandidate.h"
#include "FWCore/FWLite/interface/FWLiteEnabler.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 20 of file JPsiFWLiteAnalyzer.cc.

References fwlite::Event::atEnd(), gather_cfg::cout, FWLiteEnabler::enable(), edmPickEvents::event, fwlite::Handle< T >::getByLabel(), mps_fire::i, iEvent, and fwlite::Event::toBegin().

 {
   // ----------------------------------------------------------------------
   // First Part: 
   //
   //  * enable FWLite 
   //  * book the histograms of interest 
   //  * open the input file
   // ----------------------------------------------------------------------
 
   // load framework libraries
   gSystem->Load( "libFWCoreFWLite" );
   FWLiteEnabler::enable();
     
   // open input file (can be located on castor)
   TFile* inFile = TFile::Open( "file:jpsi.root" );
 
   // ----------------------------------------------------------------------
   // 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
   // ----------------------------------------------------------------------
 
   // loop the events
   unsigned int iEvent=0;
   fwlite::Event event(inFile);
   for(event.toBegin(); !event.atEnd(); ++event, ++iEvent){
     // break loop after end of file is reached 
     // or after 1000 events have been processed
     if( iEvent==1000 ) break;
     
     // simple event counter
     if(iEvent>0 && iEvent%1==0){
       std::cout << "  processing event: " << iEvent << std::endl;
     }
 
     // fwlite::Handle to to jpsi collection
     fwlite::Handle<std::vector<pat::CompositeCandidate> > jpsis;
     jpsis.getByLabel(event, "patJPsiCandidates");
     
     // loop jpsi collection and fill histograms
     for(unsigned i=0; i<jpsis->size(); ++i){
       cout << "jpsi " << i << ", mass = " << jpsis->at(i).mass() << ", dR = " << jpsis->at(i).userFloat("dR") << endl;
     }
   }  
   // close input file
   inFile->Close();
   
   // that's it!
   return 0;
 }

Functions

Function Documentation