|
|
|
Public Types | |
| typedef math::XYZVector | Vector |
Public Member Functions | |
| ZToLLEdmNtupleDumper (const edm::ParameterSet &) | |
Private Member Functions | |
| void | produce (edm::Event &, const edm::EventSetup &) |
Private Attributes | |
| std::vector< double > | alpha_ |
| edm::InputTag | beamSpot_ |
| std::vector< double > | beta_ |
| int | counter |
| std::vector< double > | dREcal_ |
| std::vector< double > | dRHcal_ |
| std::vector< double > | dRTrk_ |
| std::vector< double > | dRVetoTrk_ |
| std::vector< double > | etEcalThreshold_ |
| std::vector< double > | etHcalThreshold_ |
| std::vector< string > | hltPath_ |
| edm::InputTag | primaryVertices_ |
| std::vector< double > | ptThreshold_ |
| std::vector< double > | relativeIsolation_ |
| std::vector< edm::InputTag > | z_ |
| std::vector< edm::InputTag > | zGenParticlesMatch_ |
| std::vector< std::string > | zName_ |
Definition at line 37 of file ZToLLEdmNtupleDumper.cc.
Definition at line 39 of file ZToLLEdmNtupleDumper.cc.
| ZToLLEdmNtupleDumper::ZToLLEdmNtupleDumper | ( | const edm::ParameterSet & | cfg | ) |
Definition at line 105 of file ZToLLEdmNtupleDumper.cc.
References ExpressReco_HICollisions_FallBack::alias, alpha, beta, ExpressReco_HICollisions_FallBack::dREcal, ExpressReco_HICollisions_FallBack::dRHcal, edm::ParameterSet::getParameter(), i, and z.
{
string alias;
vector<ParameterSet> psets = cfg.getParameter<vector<ParameterSet> > ( "zBlocks" );
for( std::vector<edm::ParameterSet>::const_iterator i = psets.begin(); i != psets.end(); ++ i ) {
string zName = i->getParameter<string>( "zName" );
InputTag z = i->getParameter<InputTag>( "z" );
InputTag zGenParticlesMatch = i->getParameter<InputTag>( "zGenParticlesMatch" );
beamSpot_ = i->getParameter<InputTag>( "beamSpot" );
primaryVertices_= i->getParameter<InputTag>( "primaryVertices" ) ;
double ptThreshold = i->getParameter<double>("ptThreshold");
double etEcalThreshold = i->getParameter<double>("etEcalThreshold");
double etHcalThreshold= i->getParameter<double>("etHcalThreshold");
double dRVetoTrk=i->getParameter<double>("deltaRVetoTrk");
double dRTrk=i->getParameter<double>("deltaRTrk");
double dREcal=i->getParameter<double>("deltaREcal");
double dRHcal=i->getParameter<double>("deltaRHcal");
double alpha=i->getParameter<double>("alpha");
double beta=i->getParameter<double>("beta");
bool relativeIsolation = i->getParameter<bool>("relativeIsolation");
string hltPath = i ->getParameter<std::string >("hltPath");
zName_.push_back( zName );
z_.push_back( z );
zGenParticlesMatch_.push_back( zGenParticlesMatch );
ptThreshold_.push_back( ptThreshold );
etEcalThreshold_.push_back(etEcalThreshold);
etHcalThreshold_.push_back(etHcalThreshold);
dRVetoTrk_.push_back(dRVetoTrk);
dRTrk_.push_back(dRTrk);
dREcal_.push_back(dREcal);
dRHcal_.push_back(dRHcal);
alpha_.push_back(alpha);
beta_.push_back(beta);
relativeIsolation_.push_back(relativeIsolation);
hltPath_.push_back(hltPath);
produces<vector<unsigned int> >( alias = zName + "EventNumber" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "RunNumber" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "LumiBlock" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Mass" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "MassSa" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Pt" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Eta" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Phi" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Y" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1Pt" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2Pt" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1SaPt" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2SaPt" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau1HLTBit" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau2HLTBit" ).setBranchAlias( alias );
produces<vector<int> >( alias = zName + "Dau1Q" ).setBranchAlias( alias );
produces<vector<int> >( alias = zName + "Dau2Q" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1Eta" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2Eta" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1SaEta" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2SaEta" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1Phi" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2Phi" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1SaPhi" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2SaPhi" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1Iso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2Iso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1TrkIso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2TrkIso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1EcalIso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2EcalIso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1HcalIso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2HcalIso" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1MuEnergyEm" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1MuEnergyHad" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2MuEnergyEm" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2MuEnergyHad" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "VtxNormChi2" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau1NofHit" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau2NofHit" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau1NofHitTk" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau1NofHitSta" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau2NofHitTk" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau2NofHitSta" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau1NofMuChambers" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau2NofMuChambers" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau1NofMuMatches" ).setBranchAlias( alias );
produces<vector<unsigned int> >( alias = zName + "Dau2NofMuMatches" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1Chi2" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2Chi2" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1TrkChi2" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2TrkChi2" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1dxyFromBS" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2dxyFromBS" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1dzFromBS" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2dzFromBS" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1dxyFromPV" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2dxyFromPV" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau1dzFromPV" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "Dau2dzFromPV" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "TrueMass" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "TruePt" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "TrueEta" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "TruePhi" ).setBranchAlias( alias );
produces<vector<float> >( alias = zName + "TrueY" ).setBranchAlias( alias );
}
}
| void ZToLLEdmNtupleDumper::produce | ( | edm::Event & | evt, |
| const edm::EventSetup & | |||
| ) | [private, virtual] |
only for ZGolden evaluated zMassSa for the mu+sta pdf, see zmumuSaMassHistogram.cc
I fill the dau1 with positive and dau2 with negatove values for the pt, in order to flag the muons used for building zMassSa
I fill the dau1 with negatove and dau2 with positive values for the pt
Implements edm::EDProducer.
Definition at line 210 of file ZToLLEdmNtupleDumper.cc.
References trackerHits::c, reco::Muon::calEnergy(), candIsolation(), reco::Candidate::charge(), cmsDriverOptions::counter, gather_cfg::cout, reco::Candidate::daughter(), reco::MuonEnergy::em, reco::Candidate::eta(), edm::EventID::event(), event(), Exception, reco::Candidate::get(), edm::Event::getByLabel(), reco::MuonEnergy::had, reco::Candidate::hasMasterClone(), i, edm::EventBase::id(), pat::Muon::innerTrack(), edm::errors::InvalidReference, reco::Muon::isGlobalMuon(), edm::Ref< C, T, F >::isNonnull(), reco::Muon::isStandAloneMuon(), reco::Muon::isTrackerMuon(), fjr2json::lumi, edm::EventBase::luminosityBlock(), reco::Candidate::mass(), reco::Candidate::masterClone(), pat::Muon::normChi2(), reco::Muon::numberOfChambers(), reco::Muon::numberOfMatches(), pat::Muon::numberOfValidHits(), pat::Muon::outerTrack(), reco::Candidate::phi(), reco::Candidate::polarP4(), ExpressReco_HICollisions_FallBack::pt, reco::Candidate::pt(), edm::Event::put(), reco::Candidate::rapidity(), edm::EventID::run(), CrabTask::run, findQualityFiles::size, pat::GenericParticle::track(), pat::PATObject< ObjectType >::triggerObjectMatchesByPath(), reco::Candidate::vertexNormalizedChi2(), and z.
{
Handle<reco::BeamSpot> beamSpotHandle;
if (!evt.getByLabel(beamSpot_, beamSpotHandle)) {
std::cout << ">>> No beam spot found !!!"<<std::endl;
}
Handle<reco::VertexCollection> primaryVertices; // Collection of primary Vertices
if (!evt.getByLabel(primaryVertices_, primaryVertices)){
std::cout << ">>> No primary verteces found !!!"<<std::endl;
}
unsigned int size = z_.size();
for( unsigned int c = 0; c < size; ++ c ) {
Handle<CandidateView> zColl;
evt.getByLabel( z_[c], zColl );
bool isMCMatchTrue = false;
//if (zGenParticlesMatch_[c] != "") isMCMatchTrue = true;
Handle<GenParticleMatch> zGenParticlesMatch;
if (evt.getByLabel( zGenParticlesMatch_[c], zGenParticlesMatch )) {
isMCMatchTrue=true;
}
unsigned int zSize = zColl->size();
auto_ptr<vector<unsigned int> > event( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > run( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > lumi( new vector<unsigned int > );
auto_ptr<vector<float> > zMass( new vector<float> );
auto_ptr<vector<float> > zMassSa( new vector<float> );
auto_ptr<vector<float> > zPt( new vector<float> );
auto_ptr<vector<float> > zEta( new vector<float> );
auto_ptr<vector<float> > zPhi( new vector<float> );
auto_ptr<vector<float> > zY( new vector<float> );
auto_ptr<vector<float> > zDau1Pt( new vector<float> );
auto_ptr<vector<float> > zDau2Pt( new vector<float> );
auto_ptr<vector<float> > zDau1SaPt( new vector<float> );
auto_ptr<vector<float> > zDau2SaPt( new vector<float> );
auto_ptr<vector<unsigned int> > zDau1HLTBit( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau2HLTBit( new vector<unsigned int> );
auto_ptr<vector<int> > zDau1Q( new vector<int> );
auto_ptr<vector<int> > zDau2Q( new vector<int> );
auto_ptr<vector<float> > zDau1Eta( new vector<float> );
auto_ptr<vector<float> > zDau2Eta( new vector<float> );
auto_ptr<vector<float> > zDau1SaEta( new vector<float> );
auto_ptr<vector<float> > zDau2SaEta( new vector<float> );
auto_ptr<vector<float> > zDau1Phi( new vector<float> );
auto_ptr<vector<float> > zDau2Phi( new vector<float> );
auto_ptr<vector<float> > zDau1SaPhi( new vector<float> );
auto_ptr<vector<float> > zDau2SaPhi( new vector<float> );
auto_ptr<vector<float> > zDau1Iso( new vector<float> );
auto_ptr<vector<float> > zDau2Iso( new vector<float> );
auto_ptr<vector<float> > zDau1TrkIso( new vector<float> );
auto_ptr<vector<float> > zDau2TrkIso( new vector<float> );
auto_ptr<vector<float> > zDau1EcalIso( new vector<float> );
auto_ptr<vector<float> > zDau2EcalIso( new vector<float> );
auto_ptr<vector<float> > zDau1HcalIso( new vector<float> );
auto_ptr<vector<float> > zDau2HcalIso( new vector<float> );
auto_ptr<vector<float> > zDau1MuEnergyEm( new vector<float> );
auto_ptr<vector<float> > zDau2MuEnergyEm( new vector<float> );
auto_ptr<vector<float> > zDau1MuEnergyHad( new vector<float> );
auto_ptr<vector<float> > zDau2MuEnergyHad( new vector<float> );
auto_ptr<vector<float> > vtxNormChi2( new vector<float> );
auto_ptr<vector<unsigned int> > zDau1NofHit( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau2NofHit( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau1NofHitTk( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau2NofHitTk( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau1NofHitSta( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau2NofHitSta( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau1NofMuChambers( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau2NofMuChambers( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau1NofMuMatches( new vector<unsigned int> );
auto_ptr<vector<unsigned int> > zDau2NofMuMatches( new vector<unsigned int> );
auto_ptr<vector<float> > zDau1Chi2( new vector<float> );
auto_ptr<vector<float> > zDau2Chi2( new vector<float> );
auto_ptr<vector<float> > zDau1TrkChi2( new vector<float> );
auto_ptr<vector<float> > zDau2TrkChi2( new vector<float> );
auto_ptr<vector<float> > zDau1dxyFromBS( new vector<float> );
auto_ptr<vector<float> > zDau2dxyFromBS( new vector<float> );
auto_ptr<vector<float> > zDau1dzFromBS( new vector<float> );
auto_ptr<vector<float> > zDau2dzFromBS( new vector<float> );
auto_ptr<vector<float> > zDau1dxyFromPV( new vector<float> );
auto_ptr<vector<float> > zDau2dxyFromPV( new vector<float> );
auto_ptr<vector<float> > zDau1dzFromPV( new vector<float> );
auto_ptr<vector<float> > zDau2dzFromPV( new vector<float> );
auto_ptr<vector<float> > trueZMass( new vector<float> );
auto_ptr<vector<float> > trueZPt( new vector<float> );
auto_ptr<vector<float> > trueZEta( new vector<float> );
auto_ptr<vector<float> > trueZPhi( new vector<float> );
auto_ptr<vector<float> > trueZY( new vector<float> );
event -> push_back(evt.id().event());
run -> push_back(evt.id().run());
lumi -> push_back(evt.luminosityBlock());
for( unsigned int i = 0; i < zSize; ++ i ) {
const Candidate & z = (*zColl)[ i ];
CandidateBaseRef zRef = zColl->refAt(i);
zMass->push_back( z.mass() );
zPt->push_back( z.pt() );
zEta->push_back( z.eta() );
zPhi->push_back( z.phi() );
zY->push_back( z.rapidity() );
vtxNormChi2->push_back(z.vertexNormalizedChi2() );
const Candidate * dau1 = z.daughter(0);
const Candidate * dau2 = z.daughter(1);
zDau1Pt->push_back( dau1->pt() );
zDau2Pt->push_back( dau2->pt() );
zDau1Q->push_back( dau1->charge() );
zDau2Q->push_back( dau2->charge() );
zDau1Eta->push_back( dau1->eta() );
zDau2Eta->push_back( dau2->eta() );
zDau1Phi->push_back( dau1->phi() );
zDau2Phi->push_back( dau2->phi() );
if(!(dau1->hasMasterClone()&&dau2->hasMasterClone()))
throw edm::Exception(edm::errors::InvalidReference)
<< "Candidate daughters have no master clone\n";
const CandidateBaseRef & mr1 = dau1->masterClone(), & mr2 = dau2->masterClone();
const Candidate * m1 = &*mr1, * m2 = &*mr2;
// isolation as defined by us into the analyzer
double iso1 = candIsolation(m1,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], alpha_[c], beta_[c], relativeIsolation_[c]);
double iso2 = candIsolation(m2,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], alpha_[c], beta_[c], relativeIsolation_[c] );
// tracker isolation : alpha =0
double trkIso1 = candIsolation(m1,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], 0.0, beta_[c], relativeIsolation_[c]);
double trkIso2 = candIsolation(m2,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], 0.0, beta_[c], relativeIsolation_[c] );
// ecal isolation : alpha =1, beta =1
double ecalIso1 = candIsolation(m1,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], 1.0, 1.0, relativeIsolation_[c]);
double ecalIso2 = candIsolation(m2,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], 1.0, 1.0, relativeIsolation_[c] );
// hcal isolation : alpha =1, beta =-1
double hcalIso1 = candIsolation(m1,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], 1.0, -1.0, relativeIsolation_[c]);
double hcalIso2 = candIsolation(m2,ptThreshold_[c], etEcalThreshold_[c], etHcalThreshold_[c] ,dRVetoTrk_[c], dRTrk_[c], dREcal_[c] , dRHcal_[c], 1.0, -1.0, relativeIsolation_[c] );
zDau1Iso->push_back( iso1 );
zDau2Iso->push_back( iso2 );
zDau1TrkIso->push_back( trkIso1 );
zDau2TrkIso->push_back( trkIso2 );
zDau1EcalIso->push_back( ecalIso1 );
zDau2EcalIso->push_back( ecalIso2 );
zDau1HcalIso->push_back( hcalIso1 );
zDau2HcalIso->push_back( hcalIso2 );
if (isMCMatchTrue){
GenParticleRef trueZRef = (*zGenParticlesMatch)[zRef];
//CandidateRef trueZRef = trueZIter->val;
if( trueZRef.isNonnull() ) {
const Candidate & z = * trueZRef;
trueZMass->push_back( z.mass() );
trueZPt->push_back( z.pt() );
trueZEta->push_back( z.eta() );
trueZPhi->push_back( z.phi() );
trueZY->push_back( z.rapidity() );
} else {
trueZMass->push_back( -100 );
trueZPt->push_back( -100 );
trueZEta->push_back( -100 );
trueZPhi->push_back( -100 );
trueZY->push_back( -100 );
}
}
// quality variables
const pat::Muon * mu1 = dynamic_cast<const pat::Muon*>(m1);
// protection for standalone and trackerMuon
if (mu1->isGlobalMuon() == true){
zDau1NofHit->push_back(mu1->numberOfValidHits());
zDau1NofHitTk->push_back(mu1->innerTrack()->numberOfValidHits());
zDau1NofHitSta->push_back(mu1->outerTrack()->numberOfValidHits());
zDau1Chi2->push_back(mu1->normChi2());
TrackRef mu1TrkRef = mu1->innerTrack();
zDau1TrkChi2->push_back( mu1TrkRef->normalizedChi2());
zDau1dxyFromBS->push_back(mu1TrkRef->dxy(beamSpotHandle->position()));
zDau1dzFromBS->push_back(mu1TrkRef->dz(beamSpotHandle->position()));
zDau1dxyFromPV->push_back(mu1TrkRef->dxy(primaryVertices->begin()->position() ));
zDau1dzFromPV->push_back(mu1TrkRef->dz(primaryVertices->begin()->position() ));
zDau1MuEnergyEm->push_back( mu1->calEnergy().em);
zDau1MuEnergyHad->push_back( mu1->calEnergy().had);
} else if (mu1->isStandAloneMuon() == true) {
// the muon is a standalone
TrackRef mu1StaRef = mu1->outerTrack();
zDau1NofHit->push_back(mu1StaRef->numberOfValidHits());
zDau1NofHitTk->push_back(0);
zDau1NofHitSta->push_back(mu1StaRef->numberOfValidHits());
zDau1Chi2->push_back(mu1StaRef->normalizedChi2());
zDau1TrkChi2->push_back(0);
zDau1dxyFromBS->push_back(mu1StaRef->dxy(beamSpotHandle->position()));
zDau1dzFromBS->push_back(mu1StaRef->dz(beamSpotHandle->position()));
zDau1dxyFromPV->push_back(mu1StaRef->dxy(primaryVertices->begin()->position() ));
zDau1dzFromPV->push_back(mu1StaRef->dz(primaryVertices->begin()->position() ));
zDau1MuEnergyEm->push_back( -1);
zDau1MuEnergyHad->push_back( -1);
} else if (mu1->isTrackerMuon() == true) {
// the muon is a trackerMuon
TrackRef mu1TrkRef = mu1->innerTrack();
zDau1NofHit->push_back(mu1TrkRef->numberOfValidHits());
zDau1NofHitTk->push_back(mu1TrkRef->numberOfValidHits());
zDau1NofHitSta->push_back(0);
zDau1Chi2->push_back(mu1TrkRef->normalizedChi2());
zDau1TrkChi2->push_back(mu1TrkRef->normalizedChi2());
zDau1dxyFromBS->push_back(mu1TrkRef->dxy(beamSpotHandle->position()));
zDau1dzFromBS->push_back(mu1TrkRef->dz(beamSpotHandle->position()));
zDau1dxyFromPV->push_back(mu1TrkRef->dxy(primaryVertices->begin()->position() ));
zDau1dzFromPV->push_back(mu1TrkRef->dz(primaryVertices->begin()->position() ));
zDau1MuEnergyEm->push_back( mu1->calEnergy().em);
zDau1MuEnergyHad->push_back( mu1->calEnergy().had);
}
zDau1NofMuChambers->push_back(mu1->numberOfChambers());
zDau1NofMuMatches->push_back(mu1->numberOfMatches());
// would we like to add another variables???
// HLT trigger bit
const pat::TriggerObjectStandAloneCollection mu1HLTMatches = mu1->triggerObjectMatchesByPath( hltPath_[c] );
int dimTrig1 = mu1HLTMatches.size();
if(dimTrig1 !=0 ){
zDau1HLTBit->push_back(1);
} else {
zDau1HLTBit->push_back(0);
}
const pat::Muon * mu2 = dynamic_cast<const pat::Muon*>(m2);
if (mu2!=0 ) {
if (mu2->isGlobalMuon() == true) {
zDau2NofHit->push_back(mu2->numberOfValidHits());
zDau2NofHitTk->push_back(mu2->innerTrack()->numberOfValidHits());
zDau2NofHitSta->push_back(mu2->outerTrack()->numberOfValidHits());
zDau2Chi2->push_back(mu2->normChi2());
TrackRef mu2TrkRef = mu2->innerTrack();
zDau1TrkChi2->push_back( mu2TrkRef->normalizedChi2());
zDau2dxyFromBS->push_back(mu2TrkRef->dxy(beamSpotHandle->position()));
zDau2dzFromBS->push_back(mu2TrkRef->dz(beamSpotHandle->position()));
zDau2dxyFromPV->push_back(mu2TrkRef->dxy(primaryVertices->begin()->position() ));
zDau2dzFromPV->push_back(mu2TrkRef->dz(primaryVertices->begin()->position() ));
zDau2MuEnergyEm->push_back( mu2->calEnergy().em);
zDau2MuEnergyHad->push_back( mu2->calEnergy().had);
} else if (mu2->isStandAloneMuon() == true){
// its' a standalone
zDau2HLTBit->push_back(0);
TrackRef mu2StaRef = mu2->outerTrack();
zDau2NofHit->push_back(mu2StaRef->numberOfValidHits());
zDau2NofHitTk->push_back(0);
zDau2NofHitSta->push_back(mu2StaRef->numberOfValidHits());
zDau2Chi2->push_back(mu2StaRef->normalizedChi2());
zDau2TrkChi2->push_back(0);
zDau2dxyFromBS->push_back(mu2StaRef->dxy(beamSpotHandle->position()));
zDau2dzFromBS->push_back(mu2StaRef->dz(beamSpotHandle->position()));
zDau2dxyFromPV->push_back(mu2StaRef->dxy(primaryVertices->begin()->position() ));
zDau2dzFromPV->push_back(mu2StaRef->dz(primaryVertices->begin()->position() ));
zDau1MuEnergyEm->push_back( -1);
zDau1MuEnergyHad->push_back( -1);
} else if (mu2->isTrackerMuon() == true) {
// the muon is a trackerMuon
TrackRef mu2TrkRef = mu2->innerTrack();
zDau2NofHit->push_back(mu2TrkRef->numberOfValidHits());
zDau2NofHitSta->push_back(0);
zDau2NofHitTk->push_back(mu2TrkRef->numberOfValidHits());
zDau2Chi2->push_back(mu2TrkRef->normalizedChi2());
zDau2TrkChi2->push_back(mu2TrkRef->normalizedChi2());
zDau2dxyFromBS->push_back(mu2TrkRef->dxy(beamSpotHandle->position()));
zDau2dzFromBS->push_back(mu2TrkRef->dz(beamSpotHandle->position()));
zDau2dxyFromPV->push_back(mu2TrkRef->dxy(primaryVertices->begin()->position() ));
zDau2dzFromPV->push_back(mu2TrkRef->dz(primaryVertices->begin()->position() ));
zDau2MuEnergyEm->push_back( mu2->calEnergy().em);
zDau2MuEnergyHad->push_back( mu2->calEnergy().had);
}
// HLT trigger bit
const pat::TriggerObjectStandAloneCollection mu2HLTMatches = mu2->triggerObjectMatchesByPath( hltPath_[c] );
int dimTrig2 = mu2HLTMatches.size();
if(dimTrig2 !=0 ){
zDau2HLTBit->push_back(1);
}
else {
zDau2HLTBit->push_back(0);
}
if ( mu1->isGlobalMuon() && mu2->isGlobalMuon() ) {
TrackRef stAloneTrack1;
TrackRef stAloneTrack2;
Vector momentum;
Candidate::PolarLorentzVector p4_1;
double mu_mass;
stAloneTrack1 = dau1->get<TrackRef,reco::StandAloneMuonTag>();
stAloneTrack2 = dau2->get<TrackRef,reco::StandAloneMuonTag>();
zDau1SaEta->push_back(stAloneTrack1->eta());
zDau2SaEta->push_back(stAloneTrack2->eta());
zDau1SaPhi->push_back(stAloneTrack1->phi());
zDau2SaPhi->push_back(stAloneTrack2->phi());
if(counter % 2 == 0) {
momentum = stAloneTrack1->momentum();
p4_1 = dau2->polarP4();
mu_mass = dau1->mass();
zDau1SaPt->push_back(stAloneTrack1 ->pt());
zDau2SaPt->push_back(- stAloneTrack2->pt());
}else{
momentum = stAloneTrack2->momentum();
p4_1= dau1->polarP4();
mu_mass = dau2->mass();
zDau1SaPt->push_back( - stAloneTrack1->pt());
zDau2SaPt->push_back( stAloneTrack2->pt());
}
Candidate::PolarLorentzVector p4_2(momentum.rho(), momentum.eta(),momentum.phi(), mu_mass);
double mass = (p4_1+p4_2).mass();
zMassSa->push_back(mass);
++counter;
}
zDau2NofMuChambers->push_back(mu2->numberOfChambers());
zDau2NofMuMatches->push_back(mu2->numberOfMatches());
} else{
// for ZMuTk case...
// it's a track......
const pat::GenericParticle * trk2 = dynamic_cast<const pat::GenericParticle*>(m2);
TrackRef mu2TrkRef = trk2->track();
zDau2NofHit->push_back(mu2TrkRef->numberOfValidHits());
zDau2NofHitTk->push_back( mu2TrkRef->numberOfValidHits());
zDau2NofHitSta->push_back( 0);
zDau2NofMuChambers->push_back(0);
zDau2NofMuMatches->push_back(0);
zDau2Chi2->push_back(mu2TrkRef->normalizedChi2());
zDau2dxyFromBS->push_back(mu2TrkRef->dxy(beamSpotHandle->position()));
zDau2dzFromBS->push_back(mu2TrkRef->dz(beamSpotHandle->position()));
zDau2dxyFromPV->push_back(mu2TrkRef->dxy(primaryVertices->begin()->position() ));
zDau2dzFromPV->push_back(mu2TrkRef->dz(primaryVertices->begin()->position() ));
zDau1MuEnergyEm->push_back( -1);
zDau1MuEnergyHad->push_back( -1);
}
}
const string & zName = zName_[c];
evt.put( event,zName + "EventNumber" );
evt.put( run, zName + "RunNumber" );
evt.put( lumi,zName + "LumiBlock" );
evt.put( zMass, zName + "Mass" );
evt.put( zMassSa, zName + "MassSa" );
evt.put( zPt, zName + "Pt" );
evt.put( zEta, zName + "Eta" );
evt.put( zPhi, zName + "Phi" );
evt.put( zY, zName + "Y" );
evt.put( zDau1Pt, zName + "Dau1Pt" );
evt.put( zDau2Pt, zName + "Dau2Pt" );
evt.put( zDau1SaPt, zName + "Dau1SaPt" );
evt.put( zDau2SaPt, zName + "Dau2SaPt" );
evt.put( zDau1HLTBit, zName + "Dau1HLTBit" );
evt.put( zDau2HLTBit, zName + "Dau2HLTBit" );
evt.put( zDau1Q, zName + "Dau1Q" );
evt.put( zDau2Q, zName + "Dau2Q" );
evt.put( zDau1Eta, zName + "Dau1Eta" );
evt.put( zDau2Eta, zName + "Dau2Eta" );
evt.put( zDau1SaEta, zName + "Dau1SaEta" );
evt.put( zDau2SaEta, zName + "Dau2SaEta" );
evt.put( zDau1Phi, zName + "Dau1Phi" );
evt.put( zDau2Phi, zName + "Dau2Phi" );
evt.put( zDau1SaPhi, zName + "Dau1SaPhi" );
evt.put( zDau2SaPhi, zName + "Dau2SaPhi" );
evt.put( zDau1Iso, zName + "Dau1Iso" );
evt.put( zDau2Iso, zName + "Dau2Iso" );
evt.put( zDau1TrkIso, zName + "Dau1TrkIso" );
evt.put( zDau2TrkIso, zName + "Dau2TrkIso" );
evt.put( zDau1EcalIso, zName + "Dau1EcalIso" );
evt.put( zDau2EcalIso, zName + "Dau2EcalIso" );
evt.put( zDau1HcalIso, zName + "Dau1HcalIso" );
evt.put( zDau2HcalIso, zName + "Dau2HcalIso" );
evt.put( zDau1MuEnergyEm, zName + "Dau1MuEnergyEm" );
evt.put( zDau2MuEnergyEm, zName + "Dau2MuEnergyEm" );
evt.put( zDau1MuEnergyHad, zName + "Dau1MuEnergyHad" );
evt.put( zDau2MuEnergyHad, zName + "Dau2MuEnergyHad" );
evt.put( vtxNormChi2, zName + "VtxNormChi2" );
evt.put( zDau1NofHit, zName + "Dau1NofHit" );
evt.put( zDau2NofHit, zName + "Dau2NofHit" );
evt.put( zDau1NofHitTk, zName + "Dau1NofHitTk" );
evt.put( zDau2NofHitTk, zName + "Dau2NofHitTk" );
evt.put( zDau1NofHitSta, zName + "Dau1NofHitSta" );
evt.put( zDau2NofHitSta, zName + "Dau2NofHitSta" );
evt.put( zDau1NofMuChambers, zName + "Dau1NofMuChambers" );
evt.put( zDau1NofMuMatches, zName + "Dau1NofMuMatches" );
evt.put( zDau2NofMuChambers, zName + "Dau2NofMuChambers" );
evt.put( zDau2NofMuMatches, zName + "Dau2NofMuMatches" );
evt.put( zDau1Chi2, zName + "Dau1Chi2" );
evt.put( zDau2Chi2, zName + "Dau2Chi2" );
evt.put( zDau1TrkChi2, zName + "Dau1TrkChi2" );
evt.put( zDau2TrkChi2, zName + "Dau2TrkChi2" );
evt.put( zDau1dxyFromBS, zName + "Dau1dxyFromBS" );
evt.put( zDau2dxyFromBS, zName + "Dau2dxyFromBS" );
evt.put( zDau1dxyFromPV, zName + "Dau1dxyFromPV" );
evt.put( zDau2dxyFromPV, zName + "Dau2dxyFromPV" );
evt.put( zDau1dzFromBS, zName + "Dau1dzFromBS" );
evt.put( zDau2dzFromBS, zName + "Dau2dzFromBS" );
evt.put( zDau1dzFromPV, zName + "Dau1dzFromPV" );
evt.put( zDau2dzFromPV, zName + "Dau2dzFromPV" );
evt.put( trueZMass, zName + "TrueMass" );
evt.put( trueZPt, zName + "TruePt" );
evt.put( trueZEta, zName + "TrueEta" );
evt.put( trueZPhi, zName + "TruePhi" );
evt.put( trueZY, zName + "TrueY" );
}
}
std::vector<double> ZToLLEdmNtupleDumper::alpha_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
edm::InputTag ZToLLEdmNtupleDumper::beamSpot_ [private] |
Definition at line 46 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::beta_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
int ZToLLEdmNtupleDumper::counter [private] |
Definition at line 51 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::dREcal_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::dRHcal_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::dRTrk_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::dRVetoTrk_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::etEcalThreshold_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::etHcalThreshold_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
std::vector<string> ZToLLEdmNtupleDumper::hltPath_ [private] |
Definition at line 50 of file ZToLLEdmNtupleDumper.cc.
Definition at line 46 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::ptThreshold_ [private] |
Definition at line 48 of file ZToLLEdmNtupleDumper.cc.
std::vector<double> ZToLLEdmNtupleDumper::relativeIsolation_ [private] |
Definition at line 49 of file ZToLLEdmNtupleDumper.cc.
std::vector<edm::InputTag> ZToLLEdmNtupleDumper::z_ [private] |
Definition at line 45 of file ZToLLEdmNtupleDumper.cc.
std::vector<edm::InputTag> ZToLLEdmNtupleDumper::zGenParticlesMatch_ [private] |
Definition at line 45 of file ZToLLEdmNtupleDumper.cc.
std::vector<std::string> ZToLLEdmNtupleDumper::zName_ [private] |
Definition at line 44 of file ZToLLEdmNtupleDumper.cc.
1.7.3