CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
ZMuPtScaleAnalyzer Class Reference
Inheritance diagram for ZMuPtScaleAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 ZMuPtScaleAnalyzer (const edm::ParameterSet &pset)
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Member Functions

virtual void analyze (const edm::Event &event, const edm::EventSetup &setup) override
 
virtual void endJob () override
 

Private Attributes

double accEtaMax_
 
double accEtaMin_
 
double accMassMax_
 
double accMassMaxDen_
 
double accMassMin_
 
double accMassMinDen_
 
double accPtMin_
 
double angMax_
 
edm::InputTag gen_
 
TH1F * h_etaZ_
 
TH1F * h_etaZMC_
 
TH1F * h_mZ_
 
TH1F * h_mZMC_
 
TH1F * h_nZ_
 
TH1F * h_phiZ_
 
TH1F * h_phiZMC_
 
TH1F * h_ptZ_
 
TH1F * h_ptZMC_
 
TH1F * h_rapidityZ_
 
TH1F * h_rapidityZMC_
 
TH1F * h_thetaZ_
 
TH1F * h_thetaZMC_
 
TH1F * hardeta
 
TH1F * hardpt
 
double massMax_
 
int muPdgStatus_
 
unsigned int nAcc_
 
unsigned int nAccPtScaleN_
 
unsigned int nAccPtScaleP_
 
unsigned int nAccPtScaleSmearedFlat_
 
unsigned int nAccPtScaleSmearedGaus_
 
unsigned int nbinsAng_
 
unsigned int nbinsMass_
 
unsigned int nbinsPt_
 
unsigned int nBothMuHasZHasGrandMa_
 
double ptMax_
 
double ptScale_
 
TH1F * softeta
 
TH1F * softpt
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- 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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 7 of file ZMuPtScaleAnalyzer.cc.

Constructor & Destructor Documentation

ZMuPtScaleAnalyzer::ZMuPtScaleAnalyzer ( const edm::ParameterSet &  pset)

Definition at line 45 of file ZMuPtScaleAnalyzer.cc.

References angMax_, gather_cfg::cout, h_etaZMC_, h_mZMC_, h_phiZMC_, h_ptZMC_, h_rapidityZMC_, h_thetaZMC_, hardeta, hardpt, TFileDirectory::make(), massMax_, TFileService::mkdir(), nAcc_, nAccPtScaleN_, nAccPtScaleP_, nAccPtScaleSmearedFlat_, nAccPtScaleSmearedGaus_, nbinsAng_, nbinsMass_, nbinsPt_, nBothMuHasZHasGrandMa_, ptMax_, softeta, and softpt.

45  :
46  gen_(pset.getParameter<InputTag>("genParticles")),
47  nbinsMass_(pset.getUntrackedParameter<unsigned int>("nbinsMass")),
48  nbinsPt_(pset.getUntrackedParameter<unsigned int>("nbinsPt")),
49  nbinsAng_(pset.getUntrackedParameter<unsigned int>("nbinsAng")),
50  massMax_(pset.getUntrackedParameter<double>("massMax")),
51  ptMax_(pset.getUntrackedParameter<double>("ptMax")),
52  angMax_(pset.getUntrackedParameter<double>("angMax")),
53  accPtMin_(pset.getUntrackedParameter<double>("accPtMin")),
54  accMassMin_(pset.getUntrackedParameter<double>("accMassMin")),
55  accMassMax_(pset.getUntrackedParameter<double>("accMassMax")),
56  accMassMinDen_(pset.getUntrackedParameter<double>("accMassMinDen")),
57  accMassMaxDen_(pset.getUntrackedParameter<double>("accMassMaxDen")),
58  accEtaMin_(pset.getUntrackedParameter<double>("accEtaMin")),
59  accEtaMax_(pset.getUntrackedParameter<double>("accEtaMax")),
60  ptScale_(pset.getUntrackedParameter<double>("ptScale")),
61  muPdgStatus_(pset.getUntrackedParameter<int>("muPdgStatus")){
62 
63  cout << ">>> Z Histogrammer constructor" << endl;
64  Service<TFileService> fs;
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 
73  h_phiZMC_ = ZMCHisto.make<TH1F>("ZMCPhi", "Z MC #phi", nbinsAng_, -angMax_, angMax_);
74  h_thetaZMC_ = ZMCHisto.make<TH1F>("ZMCTheta", "Z MC #theta", nbinsAng_, 0, angMax_);
75  h_etaZMC_ = ZMCHisto.make<TH1F>("ZMCEta", "Z MC #eta", nbinsAng_, -angMax_, angMax_);
76  h_rapidityZMC_ = ZMCHisto.make<TH1F>("ZMCRapidity", "Z MC y", nbinsAng_, -angMax_, angMax_);
77 
78  hardeta = ZMCHisto.make<TH1F>("hard muon eta", "hard muon #eta", nbinsAng_, -angMax_, angMax_);
79  softeta = ZMCHisto.make<TH1F>("soft muon eta", "soft muon #eta", nbinsAng_, -angMax_, angMax_);
80  nAcc_=0;
81  nAccPtScaleP_=0;
82  nAccPtScaleN_=0;
83  nAccPtScaleSmearedFlat_=0;
84  nAccPtScaleSmearedGaus_=0;
85  nBothMuHasZHasGrandMa_=0;}
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
ZMuPtScaleAnalyzer::nAccPtScaleP_
unsigned int nAccPtScaleP_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::nBothMuHasZHasGrandMa_
unsigned int nBothMuHasZHasGrandMa_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::nAccPtScaleN_
unsigned int nAccPtScaleN_
Definition: ZMuPtScaleAnalyzer.cc:20
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:59
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
ZMuPtScaleAnalyzer::nAccPtScaleSmearedFlat_
unsigned int nAccPtScaleSmearedFlat_
Definition: ZMuPtScaleAnalyzer.cc:20
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
ZMuPtScaleAnalyzer::nAccPtScaleSmearedGaus_
unsigned int nAccPtScaleSmearedGaus_
Definition: ZMuPtScaleAnalyzer.cc:20

Member Function Documentation

void ZMuPtScaleAnalyzer::analyze ( const edm::Event &  event,
const edm::EventSetup &  setup 
)
overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 89 of file ZMuPtScaleAnalyzer.cc.

References funct::abs(), accEtaMax_, accEtaMin_, accMassMax_, accMassMaxDen_, accMassMin_, accMassMinDen_, accPtMin_, gather_cfg::cout, reco::LeafCandidate::eta(), f, alignCSCRings::ff, HcalObjRepresent::Fill(), relval_steps::gen(), gen_, h_etaZMC_, h_mZMC_, h_phiZMC_, h_ptZMC_, h_rapidityZMC_, h_thetaZMC_, hardeta, hardpt, i, reco::LeafCandidate::mass(), reco::CompositeRefCandidateT< D >::mother(), reco::Candidate::mother(), patZpeak::muons, muPdgStatus_, nAcc_, nAccPtScaleN_, nAccPtScaleP_, nAccPtScaleSmearedFlat_, nAccPtScaleSmearedGaus_, nBothMuHasZHasGrandMa_, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::CompositeRefCandidateT< D >::numberOfMothers(), reco::Candidate::numberOfMothers(), p4, reco::Candidate::pdgId(), reco::LeafCandidate::pdgId(), ptScale_, softeta, softpt, mathSSE::sqrt(), reco::LeafCandidate::status(), and create_public_pileup_plots::weights.

89  {
90  cout << ">>> Z HistogrammerZLONLOHistogrammer.cc analyze" << endl;
91 
92  Handle<GenParticleCollection> gen;
93  Handle<double> weights;
94 
95  event.getByLabel(gen_, gen);
96 
97 
98 
99 
100  // get weight and fill it to histogram
101 
102 
103 
104 
105  std::vector<GenParticle> muons;
106 
107  double mZGen = -100;
108 
109  for(unsigned int i = 0; i < gen->size(); ++i){
110  const GenParticle & muMC = (*gen)[i];
111  // filliScaledPng only muons coming form Z
112  if (abs(muMC.pdgId())==13 && muMC.status()==muPdgStatus_ && muMC.numberOfMothers()>0) {
113 
114  cout << "I'm getting a muon \n"
115  << "with " << "muMC.numberOfMothers() " << muMC.numberOfMothers() << "\n the first mother has pdgId " << muMC.mother()->pdgId()
116  << "with " << "muMC.mother()->numberOfMothers() " << muMC.mother()->numberOfMothers()<< "\n the first grandma has pdgId " << muMC.mother()->mother()->pdgId()<<endl;
117  cout << "with muMC.eta() " << muMC.eta()<<endl;
118  muons.push_back(muMC);
119  }
120  // introducing here the gen mass cut......................
121  /*
122  if (muPdgStatus_ ==1) {
123  mZGen = muMC.mother()->mother()->mass();
124  if (muMC.mother()->mother()->pdgId() ==23 && mZGen>accMassMinDen_ && mZGen<accMassMaxDen_ ) muons.push_back(muMC);}
125  }
126  if (muPdgStatus_ ==3) {
127  mZGen = muMC.mother()->mass();
128  if (muMC.mother()->pdgId() ==23 && mZGen>accMassMinDen_ && mZGen<accMassMaxDen_ ) muons.push_back(muMC);}
129  }
130 */
131 
132 
133  const GenParticle & zMC = (*gen)[i];
134  if (zMC.pdgId()==23 && zMC.status()==3 &&zMC.numberOfDaughters()>1 ) {
135  mZGen = zMC.mass();
136  cout << "I'm selecting a Z MC with mass " << mZGen << endl;
137  if(mZGen>accMassMinDen_ && mZGen<accMassMaxDen_ ) h_mZMC_->Fill(mZGen);
138 
139 
140 
141  }
142  }
143 
144 
145  cout << "finally I selected " << muons.size() << " muons" << endl;
146 
147 
148 
149 
150 
151 // if there are at least two muons,
152  // calculate invarant mass of first two and fill it into histogram
153 
154 
155 
156  double inv_mass = 0.0;
157  double Zpt_ = 0.0;
158  double Zeta_ = 0.0;
159  double Ztheta_ = 0.0;
160  double Zphi_ = 0.0;
161  double Zrapidity_ = 0.0;
162 
163  if(muons.size()>1) {
164 
165 
166  if (muons[0].mother()->mother()->pdgId()==23 && muons[1].mother()->mother()->pdgId()==23) nBothMuHasZHasGrandMa_ ++;
167  math::XYZTLorentzVector tot_momentum(muons[0].p4());
168  math::XYZTLorentzVector mom2(muons[1].p4());
169  tot_momentum += mom2;
170  inv_mass = sqrt(tot_momentum.mass2());
171  Zpt_=tot_momentum.pt();
172  Zeta_ = tot_momentum.eta();
173  Ztheta_ = tot_momentum.theta();
174  Zphi_ = tot_momentum.phi();
175  Zrapidity_ = tot_momentum.Rapidity();
176 
177 
178 
179 
180  double weight_sign = 1. ;
181 
182  //h_mZMC_->Fill(inv_mass);
183  h_ptZMC_->Fill(Zpt_,weight_sign);
184  h_etaZMC_->Fill(Zeta_,weight_sign);
185  h_thetaZMC_->Fill(Ztheta_,weight_sign);
186  h_phiZMC_->Fill(Zphi_,weight_sign);
187  h_rapidityZMC_-> Fill (Zrapidity_,weight_sign );
188 
189  double pt1 = muons[0].pt();
190  double pt2 = muons[1].pt();
191  double eta1 = muons[0].eta();
192  double eta2 = muons[1].eta();
193 
194 
195 
196 
197 
198  if(pt1>pt2) {
199  hardpt->Fill(pt1,weight_sign);
200  softpt->Fill(pt2,weight_sign);
201  hardeta->Fill(eta1,weight_sign);
202  softeta->Fill(eta2,weight_sign);
203  } else {
204  hardpt->Fill(pt2,weight_sign);
205  softpt->Fill(pt1,weight_sign);
206  hardeta->Fill(eta2,weight_sign);
207  softeta->Fill(eta1,weight_sign);
208  }
209 
210  //evaluating the geometric acceptance
211  if ( pt1 >= accPtMin_ && pt2 >= accPtMin_ && fabs(eta1)>= accEtaMin_ && fabs(eta2) >= accEtaMin_ && fabs(eta1)<= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass>= accMassMin_ && inv_mass<= accMassMax_) nAcc_++;
212 
213 
214  cout << "pt1" << pt1 << endl;
215 
216  // scaling the muon pt
217  double pt1ScaledP = pt1* ( 1. + ptScale_);
218  cout << "pt1 ScaledP of " << ( 1. + ptScale_) << endl;
219  cout << "pt1ScaledP" << pt1ScaledP << endl;
220 
221  double pt2ScaledP = pt2 * ( 1. + ptScale_);
222 
223 
224  //evaluating the geometric acceptance
225  if ( pt1ScaledP >= accPtMin_ && pt2ScaledP >= accPtMin_ && fabs(eta1)>= accEtaMin_ && fabs(eta2) >= accEtaMin_ && fabs(eta1)<= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass>= accMassMin_ && inv_mass<= accMassMax_)
226 nAccPtScaleP_++;
227 
228 
229 
230 
231  // scaling the muon pt
232  double pt1ScaledN = pt1* ( 1. - ptScale_);
233  double pt2ScaledN = pt2 * ( 1. - ptScale_);
234 
235 
236  //evaluating the geometric acceptance
237  if ( pt1ScaledN >= accPtMin_ && pt2ScaledN >= accPtMin_ && fabs(eta1)>= accEtaMin_ && fabs(eta2) >= accEtaMin_ && fabs(eta1)<= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass>= accMassMin_ && inv_mass<= accMassMax_)
238 nAccPtScaleN_++;
239 
240 
241  // scaling the muon pt
242  TRandom3 f;
243  f.SetSeed(123456789);
244  double pt1SmearedFlat = pt1* ( 1. + ptScale_ * f.Uniform() );
245  double pt2SmearedFlat = pt2 * ( 1. + ptScale_ * f.Uniform() ) ;
246 
247 
248  //evaluating the geometric acceptance
249  if ( pt1SmearedFlat >= accPtMin_ && pt2SmearedFlat >= accPtMin_ && fabs(eta1)>= accEtaMin_ && fabs(eta2) >= accEtaMin_ && fabs(eta1)<= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass>= accMassMin_ && inv_mass<= accMassMax_)
250  nAccPtScaleSmearedFlat_++;
251 
252 
253 // scaling the muon pt
254  TRandom3 ff;
255  ff.SetSeed(123456789);
256  double pt1SmearedGaus = pt1* ( 1. + ptScale_ * f.Gaus() );
257  double pt2SmearedGaus = pt2 * ( 1. + ptScale_ * f.Gaus() ) ;
258 
259 
260  //evaluating the geometric acceptance
261  if ( pt1SmearedGaus >= accPtMin_ && pt2SmearedGaus >= accPtMin_ && fabs(eta1)>= accEtaMin_ && fabs(eta2) >= accEtaMin_ && fabs(eta1)<= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass>= accMassMin_ && inv_mass<= accMassMax_)
262  nAccPtScaleSmearedGaus_++;
263 
264 
265 
266 
267  }}
i
int i
Definition: DBlmapReader.cc:9
reco::Candidate::mother
virtual const Candidate * mother(size_type i=0) const =0
return pointer to mother
reco::LeafCandidate::pdgId
virtual int pdgId() const GCC11_FINAL
PDG identifier.
Definition: LeafCandidate.h:211
ZMuPtScaleAnalyzer::nAccPtScaleP_
unsigned int nAccPtScaleP_
Definition: ZMuPtScaleAnalyzer.cc:20
reco::Candidate::numberOfMothers
virtual size_type numberOfMothers() const =0
number of mothers (zero or one in most of but not all the cases)
relval_steps.gen
def gen
Definition: relval_steps.py:246
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
ZMuPtScaleAnalyzer::nBothMuHasZHasGrandMa_
unsigned int nBothMuHasZHasGrandMa_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::nAccPtScaleN_
unsigned int nAccPtScaleN_
Definition: ZMuPtScaleAnalyzer.cc:20
reco::LeafCandidate::status
virtual int status() const GCC11_FINAL
status word
Definition: LeafCandidate.h:215
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
reco::CompositeRefCandidateT::numberOfMothers
virtual size_t numberOfMothers() const
number of mothers
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
p4
double p4[4]
Definition: TauolaWrapper.h:92
reco::CompositeRefCandidateT::numberOfDaughters
virtual size_t numberOfDaughters() const
number of daughters
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
f
double f[11][100]
Definition: MuScleFitUtils.cc:77
reco::LeafCandidate::eta
virtual float eta() const GCC11_FINAL
momentum pseudorapidity
Definition: LeafCandidate.h:159
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
reco::LeafCandidate::mass
virtual float mass() const GCC11_FINAL
mass
Definition: LeafCandidate.h:138
create_public_pileup_plots.weights
list weights
Definition: create_public_pileup_plots.py:205
ZMuPtScaleAnalyzer::nAccPtScaleSmearedFlat_
unsigned int nAccPtScaleSmearedFlat_
Definition: ZMuPtScaleAnalyzer.cc:20
patZpeak.muons
tuple muons
Definition: patZpeak.py:38
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
ZMuPtScaleAnalyzer::nAccPtScaleSmearedGaus_
unsigned int nAccPtScaleSmearedGaus_
Definition: ZMuPtScaleAnalyzer.cc:20
reco::CompositeRefCandidateT::mother
virtual const Candidate * mother(size_type=0) const
return mother at a given position, i = 0, ... numberOfMothers() - 1 (read only mode) ...
void ZMuPtScaleAnalyzer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 272 of file ZMuPtScaleAnalyzer.cc.

References gather_cfg::cout, h_mZMC_, nAcc_, nAccPtScaleN_, nAccPtScaleP_, nAccPtScaleSmearedFlat_, nAccPtScaleSmearedGaus_, nBothMuHasZHasGrandMa_, and mathSSE::sqrt().

272  {
273  cout << " number of events accepted :" << nAcc_ << endl;
274  cout << " number of total events :" << h_mZMC_->GetEntries() << endl;
275  cout << " number of cases in which BothMuHasZHasGrandMa :" << nBothMuHasZHasGrandMa_ << endl;
276  cout << " number of events pt scaled positively accepted :" << nAccPtScaleP_ << endl;
277 
278  cout << " number of events pt scaled negatively accepted :" << nAccPtScaleN_ << endl;
279 
280 cout << " number of events pt scaled smeared flattely accepted :" << nAccPtScaleSmearedFlat_ << endl;
281 
282 cout << " number of events pt scaled smeared gaussianely accepted :" << nAccPtScaleSmearedGaus_ << endl;
283 
284 
285  double eff = (double)nAcc_ / (double) h_mZMC_->GetEntries();
286  double err = sqrt( eff * (1. - eff) / (double) h_mZMC_->GetEntries() );
287  cout << " geometric acceptance: " << eff << "+/-" << err << endl;
288 
289  double effScaledP = (double)nAccPtScaleP_ / (double) h_mZMC_->GetEntries();
290  double errScaledP = sqrt( effScaledP * (1. - effScaledP) / (double) h_mZMC_->GetEntries() );
291  cout << " geometric acceptance when pt muon is positively scaled: " << effScaledP << "+/-" << errScaledP << endl;
292 
293  double effScaledN = (double)nAccPtScaleN_ / (double) h_mZMC_->GetEntries();
294  double errScaledN = sqrt( effScaledN * (1. - effScaledN) / (double) h_mZMC_->GetEntries() );
295  cout << " geometric acceptance when pt muon is negatively scaled: " << effScaledN << "+/-" << errScaledN << endl;
296 
297  double effSmearedFlat = (double) nAccPtScaleSmearedFlat_ / (double) h_mZMC_->GetEntries();
298  double errSmearedFlat = sqrt( effSmearedFlat * (1. - effSmearedFlat) / (double) h_mZMC_->GetEntries() );
299  cout << " geometric acceptance when pt muon is scaled with a flat smaering: " << effSmearedFlat << "+/-" << errSmearedFlat << endl;
300 
301  double effSmearedGaus = (double) nAccPtScaleSmearedGaus_ / (double) h_mZMC_->GetEntries();
302  double errSmearedGaus = sqrt( effSmearedGaus * (1. - effSmearedGaus) / (double) h_mZMC_->GetEntries() );
303  cout << " geometric acceptance when pt muon is scaled with a gaussian smearing: " << effSmearedGaus << "+/-" << errSmearedGaus << endl;
304 
305 
306 }
ZMuPtScaleAnalyzer::nAccPtScaleP_
unsigned int nAccPtScaleP_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::nBothMuHasZHasGrandMa_
unsigned int nBothMuHasZHasGrandMa_
Definition: ZMuPtScaleAnalyzer.cc:20
ZMuPtScaleAnalyzer::nAccPtScaleN_
unsigned int nAccPtScaleN_
Definition: ZMuPtScaleAnalyzer.cc:20
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
ZMuPtScaleAnalyzer::nAccPtScaleSmearedFlat_
unsigned int nAccPtScaleSmearedFlat_
Definition: ZMuPtScaleAnalyzer.cc:20
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
ZMuPtScaleAnalyzer::nAccPtScaleSmearedGaus_
unsigned int nAccPtScaleSmearedGaus_
Definition: ZMuPtScaleAnalyzer.cc:20

Member Data Documentation

double ZMuPtScaleAnalyzer::accEtaMax_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

double ZMuPtScaleAnalyzer::accEtaMin_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

double ZMuPtScaleAnalyzer::accMassMax_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

double ZMuPtScaleAnalyzer::accMassMaxDen_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

double ZMuPtScaleAnalyzer::accMassMin_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

double ZMuPtScaleAnalyzer::accMassMinDen_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

double ZMuPtScaleAnalyzer::accPtMin_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

double ZMuPtScaleAnalyzer::angMax_
private

Definition at line 15 of file ZMuPtScaleAnalyzer.cc.

Referenced by ZMuPtScaleAnalyzer().

edm::InputTag ZMuPtScaleAnalyzer::gen_
private

Definition at line 13 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

TH1F * ZMuPtScaleAnalyzer::h_etaZ_
private

Definition at line 17 of file ZMuPtScaleAnalyzer.cc.

TH1F * ZMuPtScaleAnalyzer::h_etaZMC_
private

Definition at line 18 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

TH1F * ZMuPtScaleAnalyzer::h_mZ_
private

Definition at line 17 of file ZMuPtScaleAnalyzer.cc.

TH1F* ZMuPtScaleAnalyzer::h_mZMC_
private

Definition at line 18 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), endJob(), and ZMuPtScaleAnalyzer().

TH1F* ZMuPtScaleAnalyzer::h_nZ_
private

Definition at line 17 of file ZMuPtScaleAnalyzer.cc.

TH1F * ZMuPtScaleAnalyzer::h_phiZ_
private

Definition at line 17 of file ZMuPtScaleAnalyzer.cc.

TH1F * ZMuPtScaleAnalyzer::h_phiZMC_
private

Definition at line 18 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

TH1F * ZMuPtScaleAnalyzer::h_ptZ_
private

Definition at line 17 of file ZMuPtScaleAnalyzer.cc.

TH1F * ZMuPtScaleAnalyzer::h_ptZMC_
private

Definition at line 18 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

TH1F * ZMuPtScaleAnalyzer::h_rapidityZ_
private

Definition at line 17 of file ZMuPtScaleAnalyzer.cc.

TH1F * ZMuPtScaleAnalyzer::h_rapidityZMC_
private

Definition at line 18 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

TH1F * ZMuPtScaleAnalyzer::h_thetaZ_
private

Definition at line 17 of file ZMuPtScaleAnalyzer.cc.

TH1F * ZMuPtScaleAnalyzer::h_thetaZMC_
private

Definition at line 18 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

TH1F * ZMuPtScaleAnalyzer::hardeta
private

Definition at line 19 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

TH1F* ZMuPtScaleAnalyzer::hardpt
private

Definition at line 19 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

double ZMuPtScaleAnalyzer::massMax_
private

Definition at line 15 of file ZMuPtScaleAnalyzer.cc.

Referenced by ZMuPtScaleAnalyzer().

int ZMuPtScaleAnalyzer::muPdgStatus_
private

Definition at line 21 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

unsigned int ZMuPtScaleAnalyzer::nAcc_
private

Definition at line 20 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), endJob(), and ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nAccPtScaleN_
private

Definition at line 20 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), endJob(), and ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nAccPtScaleP_
private

Definition at line 20 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), endJob(), and ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nAccPtScaleSmearedFlat_
private

Definition at line 20 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), endJob(), and ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nAccPtScaleSmearedGaus_
private

Definition at line 20 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), endJob(), and ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nbinsAng_
private

Definition at line 14 of file ZMuPtScaleAnalyzer.cc.

Referenced by ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nbinsMass_
private

Definition at line 14 of file ZMuPtScaleAnalyzer.cc.

Referenced by ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nbinsPt_
private

Definition at line 14 of file ZMuPtScaleAnalyzer.cc.

Referenced by ZMuPtScaleAnalyzer().

unsigned int ZMuPtScaleAnalyzer::nBothMuHasZHasGrandMa_
private

Definition at line 20 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), endJob(), and ZMuPtScaleAnalyzer().

double ZMuPtScaleAnalyzer::ptMax_
private

Definition at line 15 of file ZMuPtScaleAnalyzer.cc.

Referenced by ZMuPtScaleAnalyzer().

double ZMuPtScaleAnalyzer::ptScale_
private

Definition at line 16 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze().

TH1F * ZMuPtScaleAnalyzer::softeta
private

Definition at line 19 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

TH1F * ZMuPtScaleAnalyzer::softpt
private

Definition at line 19 of file ZMuPtScaleAnalyzer.cc.

Referenced by analyze(), and ZMuPtScaleAnalyzer().

Generated for CMSSW Reference Manual by doxygen 1.8.5