|
|
[Rivet-svn] r3747 - trunk/src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Fri Jun 8 15:05:41 BST 2012
Author: hoeth
Date: Fri Jun 8 15:05:41 2012
New Revision: 3747
Log:
whitespace cleanup to follow Rivet conventions
Modified:
trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc
Modified: trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc
==============================================================================
--- trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc Fri Jun 8 15:04:50 2012 (r3746)
+++ trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc Fri Jun 8 15:05:41 2012 (r3747)
@@ -117,40 +117,40 @@
// candidate electrons
ParticleVector cand_e;
foreach ( const Particle & e,
- applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt() ) {
- double eta = e.momentum().eta();
- // remove electrons with pT<15 in old veto region
- if( fabs(eta)>1.37 && fabs(eta) < 1.52 && e.momentum().perp()< 15.*GeV)
- continue;
+ applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt() ) {
+ double eta = e.momentum().eta();
+ // remove electrons with pT<15 in old veto region
+ if( fabs(eta)>1.37 && fabs(eta) < 1.52 && e.momentum().perp()< 15.*GeV)
+ continue;
double pTinCone = -e.momentum().perp();
foreach ( const Particle & track, chg_tracks ) {
if ( deltaR(e.momentum(),track.momentum()) <= 0.2 )
pTinCone += track.momentum().pT();
}
- if (pTinCone/e.momentum().perp()<0.1) {
+ if (pTinCone/e.momentum().perp()<0.1) {
cand_e.push_back(e);
- }
+ }
}
// resolve jet/lepton ambiguity
Jets recon_jets;
foreach ( const Jet& jet, cand_jets ) {
- bool away_from_e = true;
- foreach ( const Particle & e, cand_e ) {
- if ( deltaR(e.momentum(),jet.momentum()) <= 0.2 ) {
- away_from_e = false;
- break;
- }
- }
- if ( away_from_e )
- recon_jets.push_back( jet );
+ bool away_from_e = true;
+ foreach ( const Particle & e, cand_e ) {
+ if ( deltaR(e.momentum(),jet.momentum()) <= 0.2 ) {
+ away_from_e = false;
+ break;
+ }
+ }
+ if ( away_from_e )
+ recon_jets.push_back( jet );
}
// only keep electrons more than R=0.4 from jets
ParticleVector recon_e;
for(unsigned int ie=0;ie<cand_e.size();++ie) {
- const Particle & e = cand_e[ie];
- // at least 0.4 from any jets
+ const Particle & e = cand_e[ie];
+ // at least 0.4 from any jets
bool away = true;
foreach ( const Jet& jet, recon_jets ) {
if ( deltaR(e.momentum(),jet.momentum()) < 0.4 ) {
@@ -158,24 +158,24 @@
break;
}
}
- // and 0.1 from any muons
- if ( away ) {
- foreach ( const Particle & mu, cand_mu ) {
- if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
- away = false;
- break;
- }
- }
- }
- // and 0.1 from electrons
- for(unsigned int ie2=0;ie2<cand_e.size();++ie2) {
- if(ie==ie2) continue;
- if ( deltaR(e.momentum(),cand_e[ie2].momentum()) < 0.1 ) {
- away = false;
- break;
- }
- }
- // if isolated keep it
+ // and 0.1 from any muons
+ if ( away ) {
+ foreach ( const Particle & mu, cand_mu ) {
+ if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
+ }
+ // and 0.1 from electrons
+ for(unsigned int ie2=0;ie2<cand_e.size();++ie2) {
+ if(ie==ie2) continue;
+ if ( deltaR(e.momentum(),cand_e[ie2].momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
+ // if isolated keep it
if ( away )
recon_e.push_back( e );
}
@@ -183,32 +183,32 @@
// only keep muons more than R=0.4 from jets
ParticleVector recon_mu;
for(unsigned int imu=0;imu<cand_mu.size();++imu) {
- const Particle & mu = cand_mu[imu];
+ const Particle & mu = cand_mu[imu];
bool away = true;
- // at least 0.4 from any jets
+ // at least 0.4 from any jets
foreach ( const Jet& jet, recon_jets ) {
if ( deltaR(mu.momentum(),jet.momentum()) < 0.4 ) {
away = false;
break;
}
}
- // and 0.1 from any electrona
- if ( away ) {
- foreach ( const Particle & e, cand_e ) {
- if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
- away = false;
- break;
- }
- }
- }
- // and 0.1 from muons
- for(unsigned int imu2=0;imu2<cand_mu.size();++imu2) {
- if(imu==imu2) continue;
- if ( deltaR(mu.momentum(),cand_mu[imu2].momentum()) < 0.1 ) {
- away = false;
- break;
- }
- }
+ // and 0.1 from any electrona
+ if ( away ) {
+ foreach ( const Particle & e, cand_e ) {
+ if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
+ }
+ // and 0.1 from muons
+ for(unsigned int imu2=0;imu2<cand_mu.size();++imu2) {
+ if(imu==imu2) continue;
+ if ( deltaR(mu.momentum(),cand_mu[imu2].momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
if ( away )
recon_mu.push_back( mu );
}
@@ -232,89 +232,89 @@
// ATLAS calo problem
if(rand()/static_cast<double>(RAND_MAX)<=0.42) {
- foreach ( const Particle & e, recon_e ) {
- double eta = e.momentum().eta();
- double phi = e.momentum().azimuthalAngle(MINUSPI_PLUSPI);
- if(eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
- vetoEvent;
- }
- foreach ( const Jet & jet, recon_jets ) {
- double eta = jet.momentum().rapidity();
- double phi = jet.momentum().azimuthalAngle(MINUSPI_PLUSPI);
- if(jet.momentum().perp()>40 && eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
- vetoEvent;
- }
+ foreach ( const Particle & e, recon_e ) {
+ double eta = e.momentum().eta();
+ double phi = e.momentum().azimuthalAngle(MINUSPI_PLUSPI);
+ if(eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
+ vetoEvent;
+ }
+ foreach ( const Jet & jet, recon_jets ) {
+ double eta = jet.momentum().rapidity();
+ double phi = jet.momentum().azimuthalAngle(MINUSPI_PLUSPI);
+ if(jet.momentum().perp()>40 && eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
+ vetoEvent;
+ }
}
// check at least one e/mu passing trigger
if( !( !recon_e .empty() && recon_e[0] .momentum().perp()>25.) &&
- !( !recon_mu.empty() && recon_mu[0].momentum().perp()>20.) ) {
- MSG_DEBUG("Hardest lepton fails trigger");
- vetoEvent;
+ !( !recon_mu.empty() && recon_mu[0].momentum().perp()>20.) ) {
+ MSG_DEBUG("Hardest lepton fails trigger");
+ vetoEvent;
}
// calculate meff
double meff = eTmiss;
- foreach ( const Particle & e , recon_e )
- meff += e.momentum().perp();
- foreach ( const Particle & mu, recon_mu )
- meff += mu.momentum().perp();
+ foreach ( const Particle & e , recon_e )
+ meff += e.momentum().perp();
+ foreach ( const Particle & mu, recon_mu )
+ meff += mu.momentum().perp();
foreach ( const Jet & jet, recon_jets ) {
- double pT = jet.momentum().perp();
- if(pT>40.) meff += pT;
+ double pT = jet.momentum().perp();
+ if(pT>40.) meff += pT;
}
double mSFOS=1e30, mdiff=1e30,mMin=1e30;
// mass of SFOS pairs closest to the Z mass
for(unsigned int ix=0;ix<recon_e.size();++ix) {
- for(unsigned int iy=ix+1;iy<recon_e.size();++iy) {
- if(recon_e[ix].pdgId()*recon_e[iy].pdgId()>0) continue;
- double mtest = (recon_e[ix].momentum()+recon_e[iy].momentum()).mass();
- if(fabs(mtest-90.)<mdiff) {
- mSFOS = mtest;
- mdiff = fabs(mtest-90.);
- }
- else if(mMin>mtest) {
- mMin = mtest;
- }
- }
+ for(unsigned int iy=ix+1;iy<recon_e.size();++iy) {
+ if(recon_e[ix].pdgId()*recon_e[iy].pdgId()>0) continue;
+ double mtest = (recon_e[ix].momentum()+recon_e[iy].momentum()).mass();
+ if(fabs(mtest-90.)<mdiff) {
+ mSFOS = mtest;
+ mdiff = fabs(mtest-90.);
+ }
+ else if(mMin>mtest) {
+ mMin = mtest;
+ }
+ }
}
for(unsigned int ix=0;ix<recon_mu.size();++ix) {
- for(unsigned int iy=ix+1;iy<recon_mu.size();++iy) {
- if(recon_mu[ix].pdgId()*recon_mu[iy].pdgId()>0) continue;
- double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass();
- if(fabs(mtest-91.118)<mdiff) {
- mSFOS = mtest;
- mdiff = fabs(mtest-91.118);
- }
- else if(mMin>mtest) {
- mMin = mtest;
- }
- }
+ for(unsigned int iy=ix+1;iy<recon_mu.size();++iy) {
+ if(recon_mu[ix].pdgId()*recon_mu[iy].pdgId()>0) continue;
+ double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass();
+ if(fabs(mtest-91.118)<mdiff) {
+ mSFOS = mtest;
+ mdiff = fabs(mtest-91.118);
+ }
+ else if(mMin>mtest) {
+ mMin = mtest;
+ }
+ }
}
// cut to reject low mass Drell-Yan
if(mMin<=20.) {
- vetoEvent;
+ vetoEvent;
}
// make the control plots
// lepton pT
unsigned int ie=0,imu=0;
for(unsigned int ix=0;ix<4;++ix) {
- double pTe = ie <recon_e .size() ?
- recon_e [ie ].momentum().perp() : -1*GeV;
- double pTmu = imu<recon_mu.size() ?
- recon_mu[imu].momentum().perp() : -1*GeV;
- if(pTe>pTmu) {
- _hist_leptonpT [ix]->fill(pTe ,weight);
- _hist_leptonpT_MC[ix]->fill(pTe ,weight);
- ++ie;
- }
- else {
- _hist_leptonpT [ix]->fill(pTmu,weight);
- _hist_leptonpT_MC[ix]->fill(pTmu,weight);
- ++imu;
- }
+ double pTe = ie <recon_e .size() ?
+ recon_e [ie ].momentum().perp() : -1*GeV;
+ double pTmu = imu<recon_mu.size() ?
+ recon_mu[imu].momentum().perp() : -1*GeV;
+ if(pTe>pTmu) {
+ _hist_leptonpT [ix]->fill(pTe ,weight);
+ _hist_leptonpT_MC[ix]->fill(pTe ,weight);
+ ++ie;
+ }
+ else {
+ _hist_leptonpT [ix]->fill(pTmu,weight);
+ _hist_leptonpT_MC[ix]->fill(pTmu,weight);
+ ++imu;
+ }
}
// njet
_hist_njet ->fill(recon_jets.size(),weight);
@@ -323,16 +323,16 @@
_hist_etmiss ->fill(eTmiss,weight);
_hist_etmiss_MC->fill(eTmiss,weight);
if(mSFOS<1e30) {
- _hist_mSFOS ->fill(mSFOS,weight);
- _hist_mSFOS_MC->fill(mSFOS,weight);
+ _hist_mSFOS ->fill(mSFOS,weight);
+ _hist_mSFOS_MC->fill(mSFOS,weight);
}
_hist_meff ->fill(meff,weight);
_hist_meff_MC->fill(meff,weight);
// finally the counts
if(eTmiss>50.) {
- _count_SR1->fill(0.5,weight);
- if(mdiff>10.) _count_SR2->fill(0.5,weight);
+ _count_SR1->fill(0.5,weight);
+ if(mdiff>10.) _count_SR2->fill(0.5,weight);
}
}
More information about the Rivet-svn
mailing list
|