Directory Validation/Commissioning_dimuon_2010

Author: Lukas9318

README.md

@@ -13,22 +13,24 @@ This step is only needed the first time.
 cmsrel CMSSW_4_2_8
 ```
 
-## Cloning and copying the 2010-commissioning-dimuon repository from Github
+## Cloning the 2010-commissioning-dimuon repository from Github
 ```
-git clone https://github.com/cms-opendata-validation/2010-commissioning-dimuon
-cp -R 2010-commissioning-dimuon/Demo CMSSW_4_2_8/src
+git clone https://github.com/cms-opendata-validation/2010-commissioning-dimuon  Validation/Commissioning_dimuon_2010
+
 ```
 
 ## Setting up the CMS environment
 ```
-cd CMSSW_4_2_8/src/Demo/DemoAnalyzer
+cd CMSSW_4_2_8/src/Validation/Commissioning_dimuon_2010
 cmsenv
+
 ```
 
 ## Compiling and Running
 ```
 scram b
 cmsRun demoanalyzer_cfg.py
+
 ```
 
 After analysis, Commissioning00val.root file should be created. 

README.md~

@@ -0,0 +1,59 @@
+# 2010-commissioning-dimuon
+                        
+Validation code for 2010 Commissioning dataset, based on dimuon mass spectrum.
+You need to work in a Virtual Machine properly contextualized for CMS.
+
+In order to run the demoanalyzer_cfg.py to create the Commissioning ROOT files, 
+you need to create a working area and set up a proper CMS environment.
+
+## Creating the Working Area
+
+This step is only needed the first time.
+```
+cmsrel CMSSW_4_2_8
+```
+
+## Cloning and copying the 2010-commissioning-dimuon repository from Github
+```
+git clone https://github.com/cms-opendata-validation/2010-commissioning-dimuon
+cp -R 2010-commissioning-dimuon/Demo CMSSW_4_2_8/src
+```
+
+## Setting up the CMS environment
+```
+cd CMSSW_4_2_8/src/Demo/DemoAnalyzer
+cmsenv
+```
+
+## Compiling and Running
+```
+scram b
+cmsRun demoanalyzer_cfg.py
+```
+
+After analysis, Commissioning00val.root file should be created. 
+Then, is necessary to change input and also output files in demoanalyzer_cfg.py, to save these changes and to rerun program. 
+In the  repository from github are index files with names:
+ CMS_Run2010B_Commissioning_AOD_Apr21ReReco-v1_0000_file_index.txt
+ CMS_Run2010B_Commissioning_AOD_Apr21ReReco-v1_0002_file_index.txt
+ CMS_Run2010B_Commissioning_AOD_Apr21ReReco-v1_0003_file_index.txt
+ CMS_Run2010B_Commissioning_AOD_Apr21ReReco-v1_0004_file_index.txt   
+
+ 
+When you rerun for all four index files you should have four root files with names:
+ Commissioning00val.root
+ Commissioning02val.root
+ Commissioning03val.root
+ Commissioning04val.root
+
+The last thing what you should do is to merge these four root files into one root file.
+You do this as follow. In the downloaded repository is also file with name mergeCommissioning.C .
+Press command for opening ROOT program: root
+In this programe just type command: .x mergeCommissioning.C .
+This merging create a root file called CommissioningAllval.root.
+
+To look at this output, write down command: 
+```
+new TBrowser
+```
+and navigate to relevant file.

demoanalyzer_cfg.py

@@ -34,14 +34,14 @@
 # ****************************************************************
 
 # *** Commissioning data set ***
-files2010data = FileUtils.loadListFromFile ('/home/cms-opendata/CMSSW_4_2_8/src/Validation/Commissioning/datasets/CMS_Run2010B_Commissioning_AOD_Apr21ReReco-v1_0000_file_index.txt')
+files2010data = FileUtils.loadListFromFile ('/home/cms-opendata/CMSSW_4_2_8/src/Validation/Commissioning-dimuon-2010/datasets/CMS_Run2010B_Commissioning_AOD_Apr21ReReco-v1_0000_file_index.txt')
 process.source = cms.Source("PoolSource",fileNames = cms.untracked.vstring(*files2010data))
 
 # *************************************************
 # number of events to be skipped (0 by default)   *
 # *************************************************
 process.source.skipEvents = cms.untracked.uint32(0)
-process.demo = cms.EDAnalyzer('Commissioning')
+process.demo = cms.EDAnalyzer('Commissioning-dimuon-2010')
 # ***********************************************************
 # output file name                                          *
 # change this according to input file, or according to wish *

demoanalyzer_cfg.py~

@@ -0,0 +1,50 @@
+import FWCore.ParameterSet.Config as cms
+from RecoMuon.TrackingTools.MuonServiceProxy_cff import *
+import PhysicsTools.PythonAnalysis.LumiList as LumiList
+import FWCore.ParameterSet.Types as CfgTypes
+#import FWCore.PythonUtilities.LumiList as LumiList
+process = cms.Process("Validation")
+
+# intialize MessageLogger and output report
+process.load("FWCore.MessageLogger.MessageLogger_cfi")
+process.MessageLogger.cerr.threshold = 'INFO'
+process.MessageLogger.categories.append('Validation')
+process.MessageLogger.cerr.INFO = cms.untracked.PSet(limit = cms.untracked.int32(-1))
+process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
+# **********************************************************************
+# set the maximum number of events to be processed                     *
+#    this number (argument of int32) is to be modified by the user     *
+#    according to need and wish                                        *
+#    default is preset to 10000 events                                 *
+# **********************************************************************
+process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(20000000) )
+
+# ****************************************************************************
+# define the input data set here by inserting the appropriate .txt file list *
+# ****************************************************************************
+import FWCore.Utilities.FileUtils as FileUtils
+#
+# ****************************************************************
+# exactly one of the following 'files2010data =' statements      *
+# should be uncommented.                                         *
+# For MUO-10-004, the default is the Mu data set.                *
+# To run over all sets, replace '0000' by '0001' etc.            *
+# consecutively (make sure you save the output before rerunning) *
+# and add up the histograms using root tools.                    *
+# ****************************************************************
+
+# *** Commissioning data set ***
+files2010data = FileUtils.loadListFromFile ('/home/cms-opendata/CMSSW_4_2_8/src/Validation/Commissioning/datasets/CMS_Run2010B_Commissioning_AOD_Apr21ReReco-v1_0000_file_index.txt')
+process.source = cms.Source("PoolSource",fileNames = cms.untracked.vstring(*files2010data))
+
+# *************************************************
+# number of events to be skipped (0 by default)   *
+# *************************************************
+process.source.skipEvents = cms.untracked.uint32(0)
+process.demo = cms.EDAnalyzer('Commissioning')
+# ***********************************************************
+# output file name                                          *
+# change this according to input file, or according to wish *
+# ***********************************************************
+process.TFileService = cms.Service("TFileService",fileName = cms.string('Commissioning00val.root'))
+process.p = cms.Path(process.demo)

src/DemoAnalyzer.cc

@@ -1,9 +1,9 @@
 // -*- C++ -*-
 //
-// Package:    Commissioning
-// Class:      Commissioning
+// Package:    Commissioning_dimuon_2010
+// Class:      Commissioning_dimuon_2010
 // 
-/**\class Commissioning DemoAnalyzer.cc Validation/Commissioning/src/DemoAnalyzer.cc
+/**\class Commissioning DemoAnalyzer.cc Validation/Commissioning_dimuon_2010/src/DemoAnalyzer.cc
 
  Description: [one line class summary]
 
@@ -64,10 +64,10 @@
 // class declaration
 //
 
-class Commissioning: public edm::EDAnalyzer {
+class Commissioning_dimuon_2010: public edm::EDAnalyzer {
 public:
-	explicit Commissioning(const edm::ParameterSet&);
-	~Commissioning();
+	explicit Commissioning_dimuon_2010(const edm::ParameterSet&);
+	~Commissioning_dimuon_2010();
 
 private:
 	virtual void beginJob();
@@ -117,7 +117,7 @@ TH1D *h203;
 // constructors and destructor
 //
 
-Commissioning::Commissioning(const edm::ParameterSet& iConfig) {
+Commissioning::Commissioning_dimuon_2010(const edm::ParameterSet& iConfig) {
 
 // *****************************************************************
 // This is the main analysis routine
@@ -225,7 +225,7 @@ h203 = fs->make<TH1D>("Lumi section", "Lumi section", 300, 0, 3000);
 }
 
 
-Commissioning::~Commissioning() {
+Commissioning_dimuon_2010::~Commissioning_dimuon_2010() {
 	// do anything here that needs to be done at destruction time
 	// (e.g. close files, deallocate resources etc.)
 }
@@ -237,7 +237,7 @@ Commissioning::~Commissioning() {
 //
 
 // ------------ method called for each event  ------------
-void Commissioning::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
+void Commissioning_dimuon_2010::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
 
 // **********************************************
 // here each relevant event will get analyzed 
@@ -459,7 +459,7 @@ if (!providesGoodLumisection(iEvent)) {
 
 
 // ---- method called for each event to check good quality lumi section ---- //
-bool Commissioning::providesGoodLumisection(const edm::Event& iEvent) {
+bool Commissioning_dimuon_2010::providesGoodLumisection(const edm::Event& iEvent) {
 
 // check JSON "by hand"
 // This is a 'primitive' check which has the advantage that it also works 
@@ -775,15 +775,15 @@ bool Commissioning::providesGoodLumisection(const edm::Event& iEvent) {
 
 
 // ------------ method called once each job just before starting event loop  ------------
-void Commissioning::beginJob() {
+void Commissioning_dimuon_2010::beginJob() {
 
 }
 
 // ------------ method called once each job just after ending the event loop  ------------
-void Commissioning::endJob() {
+void Commissioning_dimuon_2010::endJob() {
 }
 
 //define this as a plug-in
-DEFINE_FWK_MODULE(Commissioning);
+DEFINE_FWK_MODULE(Commissioning_dimuon_2010);
 
 

src/DemoAnalyzer.cc~

@@ -1,9 +1,9 @@
 // -*- C++ -*-
 //
-// Package:    DemoAnalyzer
-// Class:      DemoAnalyzer
+// Package:    Commissioning
+// Class:      Commissioning
 // 
-/**\class DemoAnalyzer DemoAnalyzer.cc Demo/DemoAnalyzer/src/DemoAnalyzer.cc
+/**\class Commissioning DemoAnalyzer.cc Validation/Commissioning/src/DemoAnalyzer.cc
 
  Description: [one line class summary]
 
@@ -64,10 +64,10 @@
 // class declaration
 //
 
-class DemoAnalyzer: public edm::EDAnalyzer {
+class Commissioning: public edm::EDAnalyzer {
 public:
-	explicit DemoAnalyzer(const edm::ParameterSet&);
-	~DemoAnalyzer();
+	explicit Commissioning(const edm::ParameterSet&);
+	~Commissioning();
 
 private:
 	virtual void beginJob();
@@ -117,7 +117,7 @@ TH1D *h203;
 // constructors and destructor
 //
 
-DemoAnalyzer::DemoAnalyzer(const edm::ParameterSet& iConfig) {
+Commissioning::Commissioning(const edm::ParameterSet& iConfig) {
 
 // *****************************************************************
 // This is the main analysis routine
@@ -225,7 +225,7 @@ h203 = fs->make<TH1D>("Lumi section", "Lumi section", 300, 0, 3000);
 }
 
 
-DemoAnalyzer::~DemoAnalyzer() {
+Commissioning::~Commissioning() {
 	// do anything here that needs to be done at destruction time
 	// (e.g. close files, deallocate resources etc.)
 }
@@ -237,7 +237,7 @@ DemoAnalyzer::~DemoAnalyzer() {
 //
 
 // ------------ method called for each event  ------------
-void DemoAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
+void Commissioning::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
 
 // **********************************************
 // here each relevant event will get analyzed 
@@ -261,7 +261,7 @@ using namespace std;
 
 //--------------Fill basic Event-information histograms---------------------//
 // debug print commented out
-// LogInfo("Demo")<<"Event id: "<<iEvent.id()<<" Run number: "<<iEvent.run();
+// LogInfo("Validation")<<"Event id: "<<iEvent.id()<<" Run number: "<<iEvent.run();
 
 h200->Fill(iEvent.run());   // Run number which the analyzed Event belongs to
 h202->Fill((iEvent.id()).event());	// Event's EventNumber
@@ -282,12 +282,12 @@ if (!providesGoodLumisection(iEvent)) {
 
 // event is not in JSON (see above)
   h201->Fill(iEvent.run());
-  LogInfo("Demo") << "not in JSON, Event id (" << iEvent.id() << ") --> analysing anyway";
+  LogInfo("Validation") << "not in JSON, Event id (" << iEvent.id() << ") --> analysing anyway";
 }
 // analyze in any case! (use this only for Commissioning!)
 // else { // Event is to be analyzed
 
-  LogInfo("Demo")
+  LogInfo("Validation")
   << "Starting to analyze \n"
   << "Event number: " << (iEvent.id()).event()
   << ", Run number: " << iEvent.run()
@@ -344,8 +344,8 @@ if (!providesGoodLumisection(iEvent)) {
     h55->Fill(it->normalizedChi2());
 
 // the following can be uncommented if more log information is wished    
-//   LogInfo("Demo")  <<"global  muon track pointer "<<it;
-//   LogInfo("Demo")<<"global muon track p"<<it->p()<<"  global muon track pos"<<it->referencePoint()<<" global muon track vertex"<<it->vertex();
+//   LogInfo("Validation")  <<"global  muon track pointer "<<it;
+//   LogInfo("Validation")<<"global muon track p"<<it->p()<<"  global muon track pos"<<it->referencePoint()<<" global muon track vertex"<<it->vertex();
 
 //-----------------prepare variables to determine quality cuts---------------//
 // WHAT: 1) Find out the number of Hits in the current globalMuon-Track
@@ -459,7 +459,7 @@ if (!providesGoodLumisection(iEvent)) {
 
 
 // ---- method called for each event to check good quality lumi section ---- //
-bool DemoAnalyzer::providesGoodLumisection(const edm::Event& iEvent) {
+bool Commissioning::providesGoodLumisection(const edm::Event& iEvent) {
 
 // check JSON "by hand"
 // This is a 'primitive' check which has the advantage that it also works 
@@ -775,15 +775,15 @@ bool DemoAnalyzer::providesGoodLumisection(const edm::Event& iEvent) {
 
 
 // ------------ method called once each job just before starting event loop  ------------
-void DemoAnalyzer::beginJob() {
+void Commissioning::beginJob() {
 
 }
 
 // ------------ method called once each job just after ending the event loop  ------------
-void DemoAnalyzer::endJob() {
+void Commissioning::endJob() {
 }
 
 //define this as a plug-in
-DEFINE_FWK_MODULE(DemoAnalyzer);
+DEFINE_FWK_MODULE(Commissioning);