CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
ZMuPtScaleAnalyzer.cc
Go to the documentation of this file.
5 #include "TH1.h"
6 #include "TRandom3.h"
7 
9 public:
11 
12 private:
13  void analyze(const edm::Event& event, const edm::EventSetup& setup) override;
14  void endJob() override;
16  unsigned int nbinsMass_, nbinsPt_, nbinsAng_;
25 };
26 
31 
32 #include "HepMC/WeightContainer.h"
33 #include "HepMC/GenEvent.h"
39 #include <cmath>
40 #include <iostream>
41 
42 using namespace std;
43 using namespace reco;
44 using namespace edm;
45 
47  : genToken_(consumes<GenParticleCollection>(pset.getParameter<InputTag>("genParticles"))),
48  nbinsMass_(pset.getUntrackedParameter<unsigned int>("nbinsMass")),
49  nbinsPt_(pset.getUntrackedParameter<unsigned int>("nbinsPt")),
50  nbinsAng_(pset.getUntrackedParameter<unsigned int>("nbinsAng")),
51  massMax_(pset.getUntrackedParameter<double>("massMax")),
52  ptMax_(pset.getUntrackedParameter<double>("ptMax")),
53  angMax_(pset.getUntrackedParameter<double>("angMax")),
54  accPtMin_(pset.getUntrackedParameter<double>("accPtMin")),
55  accMassMin_(pset.getUntrackedParameter<double>("accMassMin")),
56  accMassMax_(pset.getUntrackedParameter<double>("accMassMax")),
57  accMassMinDen_(pset.getUntrackedParameter<double>("accMassMinDen")),
58  accMassMaxDen_(pset.getUntrackedParameter<double>("accMassMaxDen")),
59  accEtaMin_(pset.getUntrackedParameter<double>("accEtaMin")),
60  accEtaMax_(pset.getUntrackedParameter<double>("accEtaMax")),
61  ptScale_(pset.getUntrackedParameter<double>("ptScale")),
62  muPdgStatus_(pset.getUntrackedParameter<int>("muPdgStatus")) {
63  cout << ">>> Z Histogrammer constructor" << endl;
65  TFileDirectory ZMCHisto = fs->mkdir("ZMCHisto");
66 
67  h_mZMC_ = ZMCHisto.make<TH1F>("ZMCMass", "Z MC mass (GeV/c^{2})", nbinsMass_, 0, massMax_);
68  h_ptZMC_ = ZMCHisto.make<TH1F>("ZMCPt", "Z MC p_{t} (GeV/c)", nbinsPt_, 0, ptMax_);
69  hardpt = ZMCHisto.make<TH1F>("hardpt", "hard muon p_{t} (GeV/c)", nbinsPt_, 0, ptMax_);
70  softpt = ZMCHisto.make<TH1F>("softpt", "soft muon p_{t} (GeV/c)", nbinsPt_, 0, ptMax_);
71 
72  h_phiZMC_ = ZMCHisto.make<TH1F>("ZMCPhi", "Z MC #phi", nbinsAng_, -angMax_, angMax_);
73  h_thetaZMC_ = ZMCHisto.make<TH1F>("ZMCTheta", "Z MC #theta", nbinsAng_, 0, angMax_);
74  h_etaZMC_ = ZMCHisto.make<TH1F>("ZMCEta", "Z MC #eta", nbinsAng_, -angMax_, angMax_);
75  h_rapidityZMC_ = ZMCHisto.make<TH1F>("ZMCRapidity", "Z MC y", nbinsAng_, -angMax_, angMax_);
76 
77  hardeta = ZMCHisto.make<TH1F>("hard muon eta", "hard muon #eta", nbinsAng_, -angMax_, angMax_);
78  softeta = ZMCHisto.make<TH1F>("soft muon eta", "soft muon #eta", nbinsAng_, -angMax_, angMax_);
79  nAcc_ = 0;
80  nAccPtScaleP_ = 0;
81  nAccPtScaleN_ = 0;
85 }
86 
88  cout << ">>> Z HistogrammerZLONLOHistogrammer.cc analyze" << endl;
89 
92 
93  event.getByToken(genToken_, gen);
94 
95  // get weight and fill it to histogram
96 
97  std::vector<GenParticle> muons;
98 
99  double mZGen = -100;
100 
101  for (unsigned int i = 0; i < gen->size(); ++i) {
102  const GenParticle& muMC = (*gen)[i];
103  // filliScaledPng only muons coming form Z
104  if (abs(muMC.pdgId()) == 13 && muMC.status() == muPdgStatus_ && muMC.numberOfMothers() > 0) {
105  cout << "I'm getting a muon \n"
106  << "with "
107  << "muMC.numberOfMothers() " << muMC.numberOfMothers() << "\n the first mother has pdgId "
108  << muMC.mother()->pdgId() << "with "
109  << "muMC.mother()->numberOfMothers() " << muMC.mother()->numberOfMothers()
110  << "\n the first grandma has pdgId " << muMC.mother()->mother()->pdgId() << endl;
111  cout << "with muMC.eta() " << muMC.eta() << endl;
112  muons.push_back(muMC);
113  }
114  // introducing here the gen mass cut......................
115  /*
116  if (muPdgStatus_ ==1) {
117  mZGen = muMC.mother()->mother()->mass();
118  if (muMC.mother()->mother()->pdgId() ==23 && mZGen>accMassMinDen_ && mZGen<accMassMaxDen_ ) muons.push_back(muMC);}
119  }
120  if (muPdgStatus_ ==3) {
121  mZGen = muMC.mother()->mass();
122  if (muMC.mother()->pdgId() ==23 && mZGen>accMassMinDen_ && mZGen<accMassMaxDen_ ) muons.push_back(muMC);}
123  }
124 */
125 
126  const GenParticle& zMC = (*gen)[i];
127  if (zMC.pdgId() == 23 && zMC.status() == 3 && zMC.numberOfDaughters() > 1) {
128  mZGen = zMC.mass();
129  cout << "I'm selecting a Z MC with mass " << mZGen << endl;
130  if (mZGen > accMassMinDen_ && mZGen < accMassMaxDen_)
131  h_mZMC_->Fill(mZGen);
132  }
133  }
134 
135  cout << "finally I selected " << muons.size() << " muons" << endl;
136 
137  // if there are at least two muons,
138  // calculate invarant mass of first two and fill it into histogram
139 
140  double inv_mass = 0.0;
141  double Zpt_ = 0.0;
142  double Zeta_ = 0.0;
143  double Ztheta_ = 0.0;
144  double Zphi_ = 0.0;
145  double Zrapidity_ = 0.0;
146 
147  if (muons.size() > 1) {
148  if (muons[0].mother()->mother()->pdgId() == 23 && muons[1].mother()->mother()->pdgId() == 23)
150  math::XYZTLorentzVector tot_momentum(muons[0].p4());
151  math::XYZTLorentzVector mom2(muons[1].p4());
152  tot_momentum += mom2;
153  inv_mass = sqrt(tot_momentum.mass2());
154  Zpt_ = tot_momentum.pt();
155  Zeta_ = tot_momentum.eta();
156  Ztheta_ = tot_momentum.theta();
157  Zphi_ = tot_momentum.phi();
158  Zrapidity_ = tot_momentum.Rapidity();
159 
160  double weight_sign = 1.;
161 
162  //h_mZMC_->Fill(inv_mass);
163  h_ptZMC_->Fill(Zpt_, weight_sign);
164  h_etaZMC_->Fill(Zeta_, weight_sign);
165  h_thetaZMC_->Fill(Ztheta_, weight_sign);
166  h_phiZMC_->Fill(Zphi_, weight_sign);
167  h_rapidityZMC_->Fill(Zrapidity_, weight_sign);
168 
169  double pt1 = muons[0].pt();
170  double pt2 = muons[1].pt();
171  double eta1 = muons[0].eta();
172  double eta2 = muons[1].eta();
173 
174  if (pt1 > pt2) {
175  hardpt->Fill(pt1, weight_sign);
176  softpt->Fill(pt2, weight_sign);
177  hardeta->Fill(eta1, weight_sign);
178  softeta->Fill(eta2, weight_sign);
179  } else {
180  hardpt->Fill(pt2, weight_sign);
181  softpt->Fill(pt1, weight_sign);
182  hardeta->Fill(eta2, weight_sign);
183  softeta->Fill(eta1, weight_sign);
184  }
185 
186  //evaluating the geometric acceptance
187  if (pt1 >= accPtMin_ && pt2 >= accPtMin_ && fabs(eta1) >= accEtaMin_ && fabs(eta2) >= accEtaMin_ &&
188  fabs(eta1) <= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass >= accMassMin_ && inv_mass <= accMassMax_)
189  nAcc_++;
190 
191  cout << "pt1" << pt1 << endl;
192 
193  // scaling the muon pt
194  double pt1ScaledP = pt1 * (1. + ptScale_);
195  cout << "pt1 ScaledP of " << (1. + ptScale_) << endl;
196  cout << "pt1ScaledP" << pt1ScaledP << endl;
197 
198  double pt2ScaledP = pt2 * (1. + ptScale_);
199 
200  //evaluating the geometric acceptance
201  if (pt1ScaledP >= accPtMin_ && pt2ScaledP >= accPtMin_ && fabs(eta1) >= accEtaMin_ && fabs(eta2) >= accEtaMin_ &&
202  fabs(eta1) <= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass >= accMassMin_ && inv_mass <= accMassMax_)
203  nAccPtScaleP_++;
204 
205  // scaling the muon pt
206  double pt1ScaledN = pt1 * (1. - ptScale_);
207  double pt2ScaledN = pt2 * (1. - ptScale_);
208 
209  //evaluating the geometric acceptance
210  if (pt1ScaledN >= accPtMin_ && pt2ScaledN >= accPtMin_ && fabs(eta1) >= accEtaMin_ && fabs(eta2) >= accEtaMin_ &&
211  fabs(eta1) <= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass >= accMassMin_ && inv_mass <= accMassMax_)
212  nAccPtScaleN_++;
213 
214  // scaling the muon pt
215  TRandom3 f;
216  f.SetSeed(123456789);
217  double pt1SmearedFlat = pt1 * (1. + ptScale_ * f.Uniform());
218  double pt2SmearedFlat = pt2 * (1. + ptScale_ * f.Uniform());
219 
220  //evaluating the geometric acceptance
221  if (pt1SmearedFlat >= accPtMin_ && pt2SmearedFlat >= accPtMin_ && fabs(eta1) >= accEtaMin_ &&
222  fabs(eta2) >= accEtaMin_ && fabs(eta1) <= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass >= accMassMin_ &&
223  inv_mass <= accMassMax_)
225 
226  // scaling the muon pt
227  TRandom3 ff;
228  ff.SetSeed(123456789);
229  double pt1SmearedGaus = pt1 * (1. + ptScale_ * f.Gaus());
230  double pt2SmearedGaus = pt2 * (1. + ptScale_ * f.Gaus());
231 
232  //evaluating the geometric acceptance
233  if (pt1SmearedGaus >= accPtMin_ && pt2SmearedGaus >= accPtMin_ && fabs(eta1) >= accEtaMin_ &&
234  fabs(eta2) >= accEtaMin_ && fabs(eta1) <= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass >= accMassMin_ &&
235  inv_mass <= accMassMax_)
237  }
238 }
239 
241  cout << " number of events accepted :" << nAcc_ << endl;
242  cout << " number of total events :" << h_mZMC_->GetEntries() << endl;
243  cout << " number of cases in which BothMuHasZHasGrandMa :" << nBothMuHasZHasGrandMa_ << endl;
244  cout << " number of events pt scaled positively accepted :" << nAccPtScaleP_ << endl;
245 
246  cout << " number of events pt scaled negatively accepted :" << nAccPtScaleN_ << endl;
247 
248  cout << " number of events pt scaled smeared flattely accepted :" << nAccPtScaleSmearedFlat_ << endl;
249 
250  cout << " number of events pt scaled smeared gaussianely accepted :" << nAccPtScaleSmearedGaus_ << endl;
251 
252  double eff = (double)nAcc_ / (double)h_mZMC_->GetEntries();
253  double err = sqrt(eff * (1. - eff) / (double)h_mZMC_->GetEntries());
254  cout << " geometric acceptance: " << eff << "+/-" << err << endl;
255 
256  double effScaledP = (double)nAccPtScaleP_ / (double)h_mZMC_->GetEntries();
257  double errScaledP = sqrt(effScaledP * (1. - effScaledP) / (double)h_mZMC_->GetEntries());
258  cout << " geometric acceptance when pt muon is positively scaled: " << effScaledP << "+/-" << errScaledP << endl;
259 
260  double effScaledN = (double)nAccPtScaleN_ / (double)h_mZMC_->GetEntries();
261  double errScaledN = sqrt(effScaledN * (1. - effScaledN) / (double)h_mZMC_->GetEntries());
262  cout << " geometric acceptance when pt muon is negatively scaled: " << effScaledN << "+/-" << errScaledN << endl;
263 
264  double effSmearedFlat = (double)nAccPtScaleSmearedFlat_ / (double)h_mZMC_->GetEntries();
265  double errSmearedFlat = sqrt(effSmearedFlat * (1. - effSmearedFlat) / (double)h_mZMC_->GetEntries());
266  cout << " geometric acceptance when pt muon is scaled with a flat smaering: " << effSmearedFlat << "+/-"
267  << errSmearedFlat << endl;
268 
269  double effSmearedGaus = (double)nAccPtScaleSmearedGaus_ / (double)h_mZMC_->GetEntries();
270  double errSmearedGaus = sqrt(effSmearedGaus * (1. - effSmearedGaus) / (double)h_mZMC_->GetEntries());
271  cout << " geometric acceptance when pt muon is scaled with a gaussian smearing: " << effSmearedGaus << "+/-"
272  << errSmearedGaus << endl;
273 }
275 
PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
ZMuPtScaleAnalyzer::accMassMaxDen_
double accMassMaxDen_
Definition: ZMuPtScaleAnalyzer.cc:18
ZMuPtScaleAnalyzer::endJob
void endJob() override
Definition: ZMuPtScaleAnalyzer.cc:240
Handle.h
ZMuPtScaleAnalyzer::h_ptZMC_
TH1F * h_ptZMC_
Definition: ZMuPtScaleAnalyzer.cc:20
mps_fire.i
i
Definition: mps_fire.py:428
reco::CompositeRefCandidateT::numberOfDaughters
size_t numberOfDaughters() const override
number of daughters
ZMuPtScaleAnalyzer::h_phiZ_
TH1F * h_phiZ_
Definition: ZMuPtScaleAnalyzer.cc:19
ZMuPtScaleAnalyzer::accMassMin_
double accMassMin_
Definition: ZMuPtScaleAnalyzer.cc:18
ZMuPtScaleAnalyzer::softpt
TH1F * softpt
Definition: ZMuPtScaleAnalyzer.cc:21
ZMuPtScaleAnalyzer::ptScale_
double ptScale_
Definition: ZMuPtScaleAnalyzer.cc:18
ZMuPtScaleAnalyzer::hardeta
TH1F * hardeta
Definition: ZMuPtScaleAnalyzer.cc:21
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
ZMuPtScaleAnalyzer::muPdgStatus_
int muPdgStatus_
Definition: ZMuPtScaleAnalyzer.cc:24
reco::LeafCandidate::status
int status() const final
status word
Definition: LeafCandidate.h:180
ZMuPtScaleAnalyzer::h_ptZ_
TH1F * h_ptZ_
Definition: ZMuPtScaleAnalyzer.cc:19
edm::EDGetTokenT< reco::GenParticleCollection >
reco::CompositeRefCandidateT::mother
const Candidate * mother(size_type=0) const override
return mother at a given position, i = 0, ... numberOfMothers() - 1 (read only mode)
TFileDirectory::make
T * make(const Args &... args) const
make new ROOT object
Definition: TFileDirectory.h:53
edm
HLT enums.
Definition: AlignableModifier.h:19
ZMuPtScaleAnalyzer::nbinsPt_
unsigned int nbinsPt_
Definition: ZMuPtScaleAnalyzer.cc:16
gather_cfg.cout
cout
Definition: gather_cfg.py:144
reco::GenParticleCollection
std::vector< GenParticle > GenParticleCollection
collection of GenParticles
Definition: GenParticleFwd.h:13
ZMuPtScaleAnalyzer::accEtaMin_
double accEtaMin_
Definition: ZMuPtScaleAnalyzer.cc:18
reco::CompositeRefCandidateT::numberOfMothers
size_t numberOfMothers() const override
number of mothers
EDAnalyzer.h
TFileDirectory
Definition: TFileDirectory.h:24
ZMuPtScaleAnalyzer::hardpt
TH1F * hardpt
Definition: ZMuPtScaleAnalyzer.cc:21
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
ZMuPtScaleAnalyzer::nAccPtScaleSmearedFlat_
unsigned int nAccPtScaleSmearedFlat_
Definition: ZMuPtScaleAnalyzer.cc:22
reco::Candidate::mother
virtual const Candidate * mother(size_type i=0) const =0
return pointer to mother
edm::Handle
Definition: AssociativeIterator.h:50
GenParticle
Definition: GenParticle.py:1
singleTopDQM_cfi.setup
setup
Definition: singleTopDQM_cfi.py:37
ZMuPtScaleAnalyzer::nAccPtScaleN_
unsigned int nAccPtScaleN_
Definition: ZMuPtScaleAnalyzer.cc:22
CandMatchMap.h
edm::EDAnalyzer
Definition: EDAnalyzer.h:28
GenParticle.h
CandidateFwd.h
ZMuPtScaleAnalyzer::nBothMuHasZHasGrandMa_
unsigned int nBothMuHasZHasGrandMa_
Definition: ZMuPtScaleAnalyzer.cc:22
MakerMacros.h
reco::Candidate::numberOfMothers
virtual size_type numberOfMothers() const =0
number of mothers (zero or one in most of but not all the cases)
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
ZMuPtScaleAnalyzer::h_etaZMC_
TH1F * h_etaZMC_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::accPtMin_
double accPtMin_
Definition: ZMuPtScaleAnalyzer.cc:18
Service.h
alignCSCRings.ff
ff
Definition: alignCSCRings.py:148
HLT_FULL_cff.weights
weights
Definition: HLT_FULL_cff.py:99178
ZMuPtScaleAnalyzer::genToken_
edm::EDGetTokenT< reco::GenParticleCollection > genToken_
Definition: ZMuPtScaleAnalyzer.cc:15
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
HLT_FULL_cff.eta2
eta2
Definition: HLT_FULL_cff.py:9542
GenParticleFwd.h
ZMuPtScaleAnalyzer::softeta
TH1F * softeta
Definition: ZMuPtScaleAnalyzer.cc:21
HLT_FULL_cff.pt1
pt1
Definition: HLT_FULL_cff.py:9870
gen
Definition: PythiaDecays.h:13
TFileService.h
HLT_FULL_cff.eta1
eta1
Definition: HLT_FULL_cff.py:9541
ZMuPtScaleAnalyzer::h_thetaZ_
TH1F * h_thetaZ_
Definition: ZMuPtScaleAnalyzer.cc:19
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
edm::ParameterSet
Definition: ParameterSet.h:47
ZMuPtScaleAnalyzer::h_mZ_
TH1F * h_mZ_
Definition: ZMuPtScaleAnalyzer.cc:19
ZMuPtScaleAnalyzer::h_thetaZMC_
TH1F * h_thetaZMC_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::accMassMax_
double accMassMax_
Definition: ZMuPtScaleAnalyzer.cc:18
reco::LeafCandidate::mass
double mass() const final
mass
Definition: LeafCandidate.h:131
Event.h
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:152
ZMuPtScaleAnalyzer::h_nZ_
TH1F * h_nZ_
Definition: ZMuPtScaleAnalyzer.cc:19
reco::LeafCandidate::pdgId
int pdgId() const final
PDG identifier.
Definition: LeafCandidate.h:176
ZMuPtScaleAnalyzer::analyze
void analyze(const edm::Event &event, const edm::EventSetup &setup) override
Definition: ZMuPtScaleAnalyzer.cc:87
ZMuPtScaleAnalyzer::ZMuPtScaleAnalyzer
ZMuPtScaleAnalyzer(const edm::ParameterSet &pset)
Definition: ZMuPtScaleAnalyzer.cc:46
ZMuPtScaleAnalyzer::ptMax_
double ptMax_
Definition: ZMuPtScaleAnalyzer.cc:17
edm::Service< TFileService >
createfilelist.int
int
Definition: createfilelist.py:10
ZMuPtScaleAnalyzer::accMassMinDen_
double accMassMinDen_
Definition: ZMuPtScaleAnalyzer.cc:18
p4
double p4[4]
Definition: TauolaWrapper.h:92
ZMuPtScaleAnalyzer::h_etaZ_
TH1F * h_etaZ_
Definition: ZMuPtScaleAnalyzer.cc:19
edm::EventSetup
Definition: EventSetup.h:57
HLT_FULL_cff.pt2
pt2
Definition: HLT_FULL_cff.py:9872
submitPVResolutionJobs.err
err
Definition: submitPVResolutionJobs.py:85
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
InputTag.h
ZMuPtScaleAnalyzer::h_mZMC_
TH1F * h_mZMC_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::h_rapidityZMC_
TH1F * h_rapidityZMC_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::nbinsMass_
unsigned int nbinsMass_
Definition: ZMuPtScaleAnalyzer.cc:16
std
Definition: JetResolutionObject.h:76
ZMuPtScaleAnalyzer::nAccPtScaleSmearedGaus_
unsigned int nAccPtScaleSmearedGaus_
Definition: ZMuPtScaleAnalyzer.cc:22
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
ZMuPtScaleAnalyzer::accEtaMax_
double accEtaMax_
Definition: ZMuPtScaleAnalyzer.cc:18
ZMuPtScaleAnalyzer::angMax_
double angMax_
Definition: ZMuPtScaleAnalyzer.cc:17
relval_steps.gen
def gen(fragment, howMuch)
Production test section ####.
Definition: relval_steps.py:509
Candidate.h
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ParameterSet.h
ZMuPtScaleAnalyzer::nAcc_
unsigned int nAcc_
Definition: ZMuPtScaleAnalyzer.cc:22
HepMCProduct.h
ZMuPtScaleAnalyzer::nAccPtScaleP_
unsigned int nAccPtScaleP_
Definition: ZMuPtScaleAnalyzer.cc:22
ZMuPtScaleAnalyzer::massMax_
double massMax_
Definition: ZMuPtScaleAnalyzer.cc:17
event
Definition: event.py:1
edm::Event
Definition: Event.h:73
ZMuPtScaleAnalyzer::h_rapidityZ_
TH1F * h_rapidityZ_
Definition: ZMuPtScaleAnalyzer.cc:19
edm::InputTag
Definition: InputTag.h:15
ZMuPtScaleAnalyzer::nbinsAng_
unsigned int nbinsAng_
Definition: ZMuPtScaleAnalyzer.cc:16
ZMuPtScaleAnalyzer
Definition: ZMuPtScaleAnalyzer.cc:8
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
ZMuPtScaleAnalyzer::h_phiZMC_
TH1F * h_phiZMC_
Definition: ZMuPtScaleAnalyzer.cc:20