|
|
|
#include <HLTOniaSource.h>
Definition at line 28 of file HLTOniaSource.h.
| HLTOniaSource::HLTOniaSource | ( | const edm::ParameterSet & | pset | ) | [explicit] |
Definition at line 25 of file HLTOniaSource.cc.
References beamSpotTag_, edm::ParameterSet::getUntrackedParameter(), hltProcessName_, oniaMuonTag_, pixelTag_, pixelTagsAfterFilter_, subsystemFolder_, trackTag_, trackTagsAfterFilter_, triggerPath_, and triggerSummaryRAWTag_.
{
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Constructor";
//HLTrigger Path Names
std::vector<std::string> myTriggerPaths;
myTriggerPaths.push_back("HLT_Mu0_Track0_Jpsi");
myTriggerPaths.push_back("HLT_Mu3_Track0_Jpsi");
myTriggerPaths.push_back("HLT_Mu5_Track0_Jpsi");
triggerPath_ = pset.getUntrackedParameter<std::vector<std::string> >("TriggerPathNames",myTriggerPaths);
//Tag for Onia Muons
std::vector<edm::InputTag> myOniaMuonTags;
myOniaMuonTags.push_back(edm::InputTag("hltMu0TrackJpsiL3Filtered0", "", "HLT"));
myOniaMuonTags.push_back(edm::InputTag("hltMu3TrackJpsiL3Filtered3", "", "HLT"));
myOniaMuonTags.push_back(edm::InputTag("hltMu5TrackJpsiL3Filtered5", "", "HLT"));
oniaMuonTag_ = pset.getUntrackedParameter<std::vector<edm::InputTag> >("OniaMuonTag",myOniaMuonTags);
//Tag for Pixel tracks before Onia filter
pixelTag_ = pset.getUntrackedParameter<edm::InputTag>("PixelTag",edm::InputTag("hltPixelTracks", "", "HLT"));
//Tag for Tracker tracks before Onia filter
trackTag_ = pset.getUntrackedParameter<edm::InputTag>("TrackTag",edm::InputTag("hltMuTrackJpsiCtfTrackCands","", "HLT"));
beamSpotTag_ = pset.getUntrackedParameter<edm::InputTag>("BeamSpotTag",edm::InputTag("hltOfflineBeamSpot", "", "HLT"));
//Tag Trigger Summary
triggerSummaryRAWTag_ = pset.getUntrackedParameter<edm::InputTag>("TriggerSummaryTag",edm::InputTag("hltTriggerSummaryRAW", "", "HLT"));
hltProcessName_ = pset.getUntrackedParameter<std::string>("TriggerProcessName","HLT");
//Tag for Pixel tracks after Onia filter
std::vector<edm::InputTag> pxlTagsAfterFilter;
pxlTagsAfterFilter.push_back(edm::InputTag("hltMu0TrackJpsiPixelMassFiltered", "", "HLT"));
pxlTagsAfterFilter.push_back(edm::InputTag("hltMu3TrackJpsiPixelMassFiltered", "", "HLT"));
pxlTagsAfterFilter.push_back(edm::InputTag("hltMu5TrackJpsiPixelMassFiltered", "", "HLT"));
pixelTagsAfterFilter_= pset.getUntrackedParameter< std::vector<edm::InputTag> >("PixelTagAfterFilter",pxlTagsAfterFilter);
//Tag for Tracker tracks after Onia filter
std::vector<edm::InputTag> trxTagsAfterFilter;
trxTagsAfterFilter.push_back(edm::InputTag("hltMu0TrackJpsiTrackMassFiltered", "", "HLT"));
trxTagsAfterFilter.push_back(edm::InputTag("hltMu3TrackJpsiTrackMassFiltered", "", "HLT"));
trxTagsAfterFilter.push_back(edm::InputTag("hltMu5TrackJpsiTrackMassFiltered", "", "HLT"));
trackTagsAfterFilter_ = pset.getUntrackedParameter< std::vector<edm::InputTag> >("TrackTagAfterFilter",trxTagsAfterFilter);
//Foldering output
subsystemFolder_ = pset.getUntrackedParameter<std::string>("SubSystemFolder","HLT/HLTMonMuon/Onia");
}
| HLTOniaSource::~HLTOniaSource | ( | ) |
| void HLTOniaSource::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 141 of file HLTOniaSource.cc.
References beamSpotTag_, BSPosition_, fillInvariantMass(), fillOniaTriggerMEs(), edm::Event::getByLabel(), hltConfigInit_, i, edm::HandleBase::isValid(), label, edm::InputTag::label(), muonME_, oniaMuonTag_, pixelAfterFilterME_, pixelME_, pixelTag_, pixelTagsAfterFilter_, findQualityFiles::size, python::multivaluedict::sort(), trackAfterFilterME_, trackME_, trackTag_, trackTagsAfterFilter_, trigger::TriggerMuon, triggerSummaryRAWTag_, and trigger::TriggerTrack.
{
if(!hltConfigInit_) return;
//Get Pixel Tracks
edm::Handle<reco::TrackCollection> pixelCands;
iEvent.getByLabel(pixelTag_, pixelCands);
reco::TrackCollection mypixelCands; //This is needed for the sort!!!!
if (pixelCands.isValid()) {
mypixelCands = * pixelCands;
sort(mypixelCands.begin(), mypixelCands.end(),PtGreater());
this->fillOniaTriggerMEs(pixelCands , pixelTag_.label(), pixelME_ );
}else {
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Could not access pixel collection with tag "<<pixelTag_;
}
//Get Tracker Tracks
edm::Handle<reco::RecoChargedCandidateCollection> trackCands;
iEvent.getByLabel(trackTag_, trackCands);
reco::RecoChargedCandidateCollection mytrackCands; //This is needed for the sort!!!!
if(trackCands.isValid()) {
mytrackCands = * trackCands;
sort(mytrackCands.begin(),mytrackCands.end(),PtGreater());
this->fillOniaTriggerMEs(trackCands , trackTag_.label(), trackME_ );
}else {
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Could not access track collection with tag "<<trackTag_;
}
//Get Beamspot
edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
iEvent.getByLabel(beamSpotTag_, recoBeamSpotHandle);
if (recoBeamSpotHandle.isValid()) {
BSPosition_ = recoBeamSpotHandle->position();
}else {
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Could not access beam spot info with tag "<<beamSpotTag_;
}
//Get Trigger Summary RA
edm::Handle<trigger::TriggerEventWithRefs> rawTriggerEvent;
iEvent.getByLabel(triggerSummaryRAWTag_, rawTriggerEvent );
if( rawTriggerEvent.isValid() ){
for(size_t i=0; i<oniaMuonTag_.size(); i++){
std::vector<reco::RecoChargedCandidateRef> myMuonFilterCands;
std::vector<reco::RecoChargedCandidateRef> myPixelFilterCands;
std::vector<reco::RecoChargedCandidateRef> myTrackFilterCands;
//Get Onia Muons
size_t indexM = rawTriggerEvent->filterIndex(oniaMuonTag_[i]);
if ( indexM < rawTriggerEvent->size() ){
rawTriggerEvent->getObjects( indexM, trigger::TriggerMuon, myMuonFilterCands );
this->fillOniaTriggerMEs( myMuonFilterCands, oniaMuonTag_[i].label(), muonME_ );
}else{
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Did not find muons with tag "<<oniaMuonTag_[i];
}
//Get Onia Pixel
size_t indexP = rawTriggerEvent->filterIndex(pixelTagsAfterFilter_[i]);
if ( indexP < rawTriggerEvent->size() ){
rawTriggerEvent->getObjects( indexP, trigger::TriggerTrack , myPixelFilterCands );
this->fillOniaTriggerMEs( myPixelFilterCands,pixelTagsAfterFilter_[i].label(), pixelAfterFilterME_);
sort(myPixelFilterCands.begin(), myPixelFilterCands.end(),PtGreaterRef());
}else{
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Did not find pixel with tag "<<pixelTagsAfterFilter_[i];
}
//Get Onia Tracker Tracks
size_t indexT = rawTriggerEvent->filterIndex(trackTagsAfterFilter_[i]);
if ( indexT < rawTriggerEvent->size() ){
rawTriggerEvent->getObjects( indexT, trigger::TriggerTrack , myTrackFilterCands );
this->fillOniaTriggerMEs( myTrackFilterCands,trackTagsAfterFilter_[i].label(), trackAfterFilterME_ );
sort(myTrackFilterCands.begin(), myTrackFilterCands.end(),PtGreaterRef());
}else{
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Did not find tracks with tag "<<trackTagsAfterFilter_[i];
}
if( myMuonFilterCands.size() > 0){
if ( myPixelFilterCands.size() > 0 )this->fillInvariantMass( myMuonFilterCands, myPixelFilterCands ,oniaMuonTag_[i].label(),pixelTagsAfterFilter_[i].label());
if ( myTrackFilterCands.size() > 0 )this->fillInvariantMass( myMuonFilterCands, myTrackFilterCands ,oniaMuonTag_[i].label(),trackTagsAfterFilter_[i].label());
if (pixelCands.isValid()) this->fillInvariantMass( myMuonFilterCands, mypixelCands , oniaMuonTag_[i].label(), pixelTag_.label());
if (trackCands.isValid()) this->fillInvariantMass( myMuonFilterCands, mytrackCands , oniaMuonTag_[i].label(), trackTag_.label());
}
}//ADD INVARIANT MASSES
}else{
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Could not access trigger collection with tag "<<triggerSummaryRAWTag_;
}
}
| void HLTOniaSource::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 76 of file HLTOniaSource.cc.
References dbe_, and cmsCodeRules::cppFunctionSkipper::operator.
{
dbe_ = edm::Service<DQMStore>().operator->();
if( !dbe_ ) {
edm::LogVerbatim ("oniatriggermonitor") << "[HLTOniaSource]: Could not access DQM Store.";
return;
}
}
| void HLTOniaSource::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | setup | ||
| ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 85 of file HLTOniaSource.cc.
References bookOniaTriggerInvariantMassMEs(), bookOniaTriggerMEs(), checkHLTConfiguration(), dbe_, hltConfigInit_, hltProcessName_, i, label, edm::InputTag::label(), massME_, muonME_, oniaMuonTag_, pixelAfterFilterME_, pixelME_, pixelTag_, pixelTagsAfterFilter_, DQMStore::setCurrentFolder(), subsystemFolder_, trackAfterFilterME_, trackME_, trackTag_, and trackTagsAfterFilter_.
{
if (!dbe_) return;
// Book Pixel Histos
if (pixelTag_.label()!= ""){
dbe_->setCurrentFolder(subsystemFolder_+"/Pixel");
this->bookOniaTriggerMEs(pixelME_, pixelTag_.label());
}
// Book Track Histos
if (trackTag_.label()!= ""){
dbe_->setCurrentFolder(subsystemFolder_+"/Track");
this->bookOniaTriggerMEs(trackME_, trackTag_.label());
}
//Book Onia Histos
for (size_t i = 0 ; i<oniaMuonTag_.size() && i<pixelTagsAfterFilter_.size() && i<trackTagsAfterFilter_.size(); i++){
if (oniaMuonTag_[i].label()!= "") {
dbe_->setCurrentFolder(subsystemFolder_+"/MuonFilters");
this->bookOniaTriggerMEs(muonME_, oniaMuonTag_[i].label());
}
if (pixelTagsAfterFilter_[i].label() != ""){
dbe_->setCurrentFolder(subsystemFolder_+"/PixelFilters");
this->bookOniaTriggerMEs(pixelAfterFilterME_, pixelTagsAfterFilter_[i].label() );
}
if (trackTagsAfterFilter_[i].label() != ""){
dbe_->setCurrentFolder(subsystemFolder_+"/TrackFilters");
this->bookOniaTriggerMEs(trackAfterFilterME_, trackTagsAfterFilter_[i].label() );
}
if (oniaMuonTag_[i].label()!= "" && pixelTag_.label() != ""){
dbe_->setCurrentFolder(subsystemFolder_+"/Pixel");
this->bookOniaTriggerInvariantMassMEs( massME_, oniaMuonTag_[i].label(),pixelTag_.label() );
}
if (oniaMuonTag_[i].label() != "" && trackTag_.label()!= "" ){
dbe_->setCurrentFolder(subsystemFolder_+"/Track");
this->bookOniaTriggerInvariantMassMEs( massME_, oniaMuonTag_[i].label(), trackTag_.label() );
}
if (oniaMuonTag_[i].label()!= "" && pixelTagsAfterFilter_[i].label()!= ""){
dbe_->setCurrentFolder(subsystemFolder_+"/PixelFilters");
this->bookOniaTriggerInvariantMassMEs( massME_,oniaMuonTag_[i].label(),pixelTagsAfterFilter_[i].label() );
}
if (oniaMuonTag_[i].label()!= "" && trackTagsAfterFilter_[i].label()!= ""){
dbe_->setCurrentFolder(subsystemFolder_+"/TrackFilters");
this->bookOniaTriggerInvariantMassMEs( massME_, oniaMuonTag_[i].label(), trackTagsAfterFilter_[i].label() );
}
}
hltConfigInit_ = this->checkHLTConfiguration(run , setup , hltProcessName_);
}
| void HLTOniaSource::bookOniaTriggerInvariantMassMEs | ( | std::map< std::string, MonitorElement * > & | myMap, |
| std::string | label1, | ||
| std::string | label2 | ||
| ) | [private] |
Definition at line 326 of file HLTOniaSource.cc.
References DQMStore::book1D(), dbe_, and massME_.
Referenced by beginRun().
{
std::stringstream meName;
//Same charge
meName.str("");
meName<<label1<<"_"<<label2<<"_SameCharge_InvariantMass";
massME_[label1+label2+"same"]=dbe_->book1D(meName.str(), meName.str(),120, 0.0, 6.0);
massME_[label1+label2+"same"]->setAxisTitle("mass (GeV)", 1);
//Opposite charge
meName.str("");
meName<<label1<<"_"<<label2<<"_OppositeCharge_InvariantMass";
massME_[label1+label2+"opposite"]=dbe_->book1D(meName.str(), meName.str(),120, 0.0, 6.0);
massME_[label1+label2+"opposite"]->setAxisTitle("mass (GeV)", 1);
//Same charge Highest PT
meName.str("");
meName<<label1<<"_"<<label2<<"_SameCharge_HighestPT_InvariantMass";
massME_[label1+label2+"same"+"highestpt"]=dbe_->book1D(meName.str(), meName.str(),120, 0.0, 6.0);
massME_[label1+label2+"same"+"highestpt"]->setAxisTitle("mass (GeV)", 1);
//Opposite charge Highest PT
meName.str("");
meName<<label1<<"_"<<label2<<"_OppositeCharge_HighestPT_InvariantMass";
massME_[label1+label2+"opposite"+"highestpt"]=dbe_->book1D(meName.str(), meName.str(),120, 0.0, 6.0);
massME_[label1+label2+"opposite"+"highestpt"]->setAxisTitle("mass (GeV)", 1);
// Same Charge Deltaz Muon - Track Highest PT
meName.str("");
meName<<label1<<"_"<<label2<<"_SameCharge_HighestPT_MaxDzMuonTrack";
massME_[label1+label2+"same"+"highestpt"+"maxdzmuontrack"]=dbe_->book1D(meName.str(), meName.str(),50, 0.0, 5.0);
massME_[label1+label2+"same"+"highestpt"+"maxdzmuontrack"]->setAxisTitle("MaxDz Muon-Track", 1);
// Same Charge Deltaz Muon - Track
meName.str("");
meName<<label1<<"_"<<label2<<"_SameCharge_MaxDzMuonTrack";
massME_[label1+label2+"same"+"maxdzmuontrack"]=dbe_->book1D(meName.str(), meName.str(),50, 0.0, 5.0);
massME_[label1+label2+"same"+"maxdzmuontrack"]->setAxisTitle("MaxDz Muon-Track", 1);
// Opposite charge Deltaz Muon - Track Highest PT
meName.str("");
meName<<label1<<"_"<<label2<<"_OppositeCharge_HighestPT_MaxDzMuonTrack";
massME_[label1+label2+"opposite"+"highestpt"+"maxdzmuontrack"]=dbe_->book1D(meName.str(), meName.str(),50, 0.0, 5.0);
massME_[label1+label2+"opposite"+"highestpt"+"maxdzmuontrack"]->setAxisTitle("MaxDz Muon-Track", 1);
// Opposite charge Deltaz Muon - Track
meName.str("");
meName<<label1<<"_"<<label2<<"_OppositeCharge_MaxDzMuonTrack";
massME_[label1+label2+"opposite"+"maxdzmuontrack"]=dbe_->book1D(meName.str(), meName.str(),50, 0.0, 5.0);
massME_[label1+label2+"opposite"+"maxdzmuontrack"]->setAxisTitle("MaxDz Muon-Track", 1);
}
| void HLTOniaSource::bookOniaTriggerMEs | ( | std::map< std::string, MonitorElement * > & | myMap, |
| std::string | meName | ||
| ) | [private] |
Definition at line 239 of file HLTOniaSource.cc.
References DQMStore::book1D(), DQMStore::book2D(), dbe_, and MonitorElement::setAxisTitle().
Referenced by beginRun().
{
std::stringstream myMeName;
//PT
myMeName.str("");
myMeName<<meName<<"_Pt";
myMap[meName+"pt"]=dbe_->book1D(myMeName.str(), myMeName.str(),200, 0.0, 20.0);
myMap[meName+"pt"]->setAxisTitle("Pt (GeV)", 1);
//P
myMeName.str("");
myMeName<<meName<<"_P";
myMap[meName+"p"]=dbe_->book1D(myMeName.str(), myMeName.str(),250, 0.0, 50.0);
myMap[meName+"p"]->setAxisTitle("P (GeV)", 1);
//Eta
myMeName.str("");
myMeName<<meName<<"_Eta";
myMap[meName+"eta"]=dbe_->book1D(myMeName.str(), myMeName.str(), 50, -2.5, 2.5 );
myMap[meName+"eta"]->setAxisTitle("Eta", 1);
//Phi
myMeName.str("");
myMeName<<meName<<"_Phi";
myMap[meName+"phi"]=dbe_->book1D(myMeName.str(), myMeName.str(), 144, -3.1416, 3.1416 );
myMap[meName+"phi"]->setAxisTitle("Phi", 1);
//Phi
myMeName.str("");
myMeName<<meName<<"_Phi";
myMap[meName+"phi"]=dbe_->book1D(myMeName.str(), myMeName.str(), 144, -3.1416, 3.1416 );
myMap[meName+"phi"]->setAxisTitle("Phi", 1);
//Charge
myMeName.str("");
myMeName<<meName<<"_Charge";
myMap[meName+"charge"]=dbe_->book1D(myMeName.str(), myMeName.str(), 3, -1.5, 1.5 );
myMap[meName+"charge"]->setAxisTitle("charge", 1);
//Dz
myMeName.str("");
myMeName<<meName<<"_Dz";
myMap[meName+"dz"]=dbe_->book1D(myMeName.str(), myMeName.str(), 400, -20.0, 20.0 );
myMap[meName+"dz"]->setAxisTitle("dz", 1);
//Dxy
myMeName.str("");
myMeName<<meName<<"_Dxy";
myMap[meName+"dxy"]=dbe_->book1D(myMeName.str(), myMeName.str(), 100, -0.5, 0.5 );
myMap[meName+"dxy"]->setAxisTitle("dxy", 1);
//EtaVsPhi
myMeName.str("");
myMeName<<meName<<"_EtaPhi";
myMap[meName+"etaphi"]=dbe_->book2D(myMeName.str(), myMeName.str(), 144, -3.1416, 3.1416 ,50, -2.5, 2.5 );
myMap[meName+"etaphi"]->setAxisTitle("Phi", 1);
myMap[meName+"etaphi"]->setAxisTitle("Eta", 2);
//EtaVsPt
myMeName.str("");
myMeName<<meName<<"_EtaPt";
myMap[meName+"etapt"]=dbe_->book2D(myMeName.str(), myMeName.str(), 100, 0.0, 100.0, 50, -2.5, 2.5 );
myMap[meName+"etapt"]->setAxisTitle("Pt (GeV)", 1);
myMap[meName+"etapt"]->setAxisTitle("Eta", 2);
//ValidHits
myMeName.str("");
myMeName<<meName<<"_ValidHits";
myMap[meName+"validhits"]=dbe_->book1D(myMeName.str(), myMeName.str(), 50, 0.0, 50.0 );
myMap[meName+"validhits"]->setAxisTitle("ValidHits", 1);
//Norm Chi2
myMeName.str("");
myMeName<<meName<<"_NormChi2";
myMap[meName+"normchi"]=dbe_->book1D(myMeName.str(), myMeName.str(), 300, 0.0, 30.0 );
myMap[meName+"normchi"]->setAxisTitle("Normalized Chi2", 1);
//Number Of Candi
myMeName.str("");
myMeName<<meName<<"_NrCandidates";
myMap[meName+"nrcand"]=dbe_->book1D(myMeName.str(), myMeName.str(), 50, 0.0, 50.0 );
myMap[meName+"nrcand"]->setAxisTitle("Nr per Event", 1);
//return true;
}
| bool HLTOniaSource::checkHLTConfiguration | ( | const edm::Run & | run, |
| const edm::EventSetup & | setup, | ||
| std::string | triggerProcessName | ||
| ) | [private] |
Definition at line 597 of file HLTOniaSource.cc.
References spr::find(), HltComparatorCreateWorkflow::hltConfig, i, HLTConfigProvider::init(), j, HLTConfigProvider::moduleLabels(), AlCaRecoCosmics_cfg::name, HLTConfigProvider::triggerNames(), and triggerPath_.
Referenced by beginRun().
{
HLTConfigProvider hltConfig;
bool changed(false);
if(hltConfig.init(run , setup, triggerProcessName, changed)){
edm::LogVerbatim("hltoniasource") << "Successfully initialized HLTConfigProvider with process name: "<<triggerProcessName;
std::stringstream os;
std::vector<std::string> triggerNames = hltConfig.triggerNames();
for( size_t i = 0; i < triggerNames.size(); i++) {
if (find(triggerPath_.begin(), triggerPath_.end(), triggerNames[i]) == triggerPath_.end()) continue;
edm::LogVerbatim("hltoniasource") << "[HLTOniaSource]: Trigger Path: "<<triggerNames[i];
std::vector<std::string> moduleNames = hltConfig.moduleLabels( triggerNames[i] );
for( size_t j = 0; j < moduleNames.size(); j++) {
TString name = moduleNames[j];
edm::LogVerbatim("hltoniasource") << "\t Fliter Name: "<<moduleNames[j];
}
}
return true;
}else{
edm::LogVerbatim("hltoniasource") << "Could not initialize HLTConfigProvider with process name: "<<triggerProcessName;
return false;
}
return true;
}
| void HLTOniaSource::endJob | ( | void | ) | [private, virtual] |
| void HLTOniaSource::fillInvariantMass | ( | std::vector< reco::RecoChargedCandidateRef > & | cand1, |
| reco::RecoChargedCandidateCollection & | cand2, | ||
| std::string | cand1Label, | ||
| std::string | cand2Label | ||
| ) | [private] |
Definition at line 511 of file HLTOniaSource.cc.
References DeDxDiscriminatorTools::charge(), i, massME_, and p4.
Referenced by analyze().
{
typedef reco::RecoChargedCandidateCollection::const_iterator cand;
//Loop on collection to calculate invariate mass
for(size_t i = 0 ; i< cand1.size(); i++) {
//Highest PT
if(cand2.begin() != cand2.end() && cand2.begin()->p()>3) {
cand candItr = cand2.begin();
std::string chargeLabel = "same";
//Check relative charge
if(cand1[i]->charge() * candItr->charge() < 0) chargeLabel = "opposite";
if(massME_[cand1Label+cand2Label+chargeLabel+"highestpt"]){
massME_[cand1Label+cand2Label+chargeLabel+"highestpt"]->Fill((cand1[i]->p4()+candItr->p4()).mass());
}
if(massME_[cand1Label+cand2Label+chargeLabel+"highestpt"+"maxdzmuontrack"]){
reco::RecoChargedCandidate tk1 = (*cand1[i]);
reco::RecoChargedCandidate tk2 = (*candItr);
// massME_[cand1Label+cand2Label+chargeLabel+"highestpt"+"maxdzmuontrack"]->Fill(fabs(tk1->dz(BSPosition_) - tk2->dz(BSPosition_)));
}
}
for (cand candItr2= cand2.begin(); candItr2!=cand2.end(); candItr2++) {
if(candItr2->p() < 3) continue; //Check if momentum is greater than 3GeV.
std::string chargeLabel = "same";
//Check relative charge
if(cand1[i]->charge() * candItr2->charge() < 0) chargeLabel = "opposite";
if(massME_[cand1Label+cand2Label+chargeLabel]){
massME_[cand1Label+cand2Label+chargeLabel]->Fill((cand1[i]->p4()+candItr2->p4()).mass());
}
if(massME_[cand1Label+cand2Label+chargeLabel+"maxdzmuontrack"]){
reco::RecoChargedCandidate tk1 = (*cand1[i]);
reco::RecoChargedCandidate tk2 = (*candItr2);
// massME_[cand1Label+cand2Label+chargeLabel+"maxdzmuontrack"]->Fill(fabs(tk1->dz(BSPosition_) - tk2->dz(BSPosition_)));
}
}
}
}
| void HLTOniaSource::fillInvariantMass | ( | std::vector< reco::RecoChargedCandidateRef > & | cand1, |
| std::vector< reco::RecoChargedCandidateRef > & | cand2, | ||
| std::string | cand1Label, | ||
| std::string | cand2Label | ||
| ) | [private] |
Definition at line 470 of file HLTOniaSource.cc.
References DeDxDiscriminatorTools::charge(), i, j, massME_, L1TEmulatorMonitor_cff::p, and p4.
{
//Loop on collection to calculate invariate mass
for(size_t i = 0 ; i< cand1.size(); i++) {
//Highest PT
std::string chargeLabel = "same";
//Check relative charge
if(cand1[i]->charge() * cand2[0]->charge() < 0) chargeLabel = "opposite";
if(massME_[cand1Label+cand2Label+chargeLabel+"highestpt"]){
massME_[cand1Label+cand2Label+chargeLabel+"highestpt"]->Fill((cand1[i]->p4()+cand2[0]->p4()).mass());
}
if(massME_[cand1Label+cand2Label+chargeLabel+"highestpt"+"maxdzmuontrack"]){
reco::RecoChargedCandidate tk1 = (*cand1[i]);
reco::RecoChargedCandidate tk2 = (*cand2[0]);
// TrackRef tk1 = cand1[i]->get<TrackRef>();
// massME_[cand1Label+cand2Label+chargeLabel+"highestpt"+"maxdzmuontrack"]->Fill(fabs(tk1->dz(BSPosition_) - tk2->dz(BSPosition_)));
}
for (size_t j= 0; j< cand2.size(); j++) {
if(cand2[j]->p() < 3) continue; //Check if momentum is greater than 3GeV.
std::string chargeLabel = "same";
//Check relative charge
if(cand1[i]->charge() * cand2[j]->charge() < 0) chargeLabel = "opposite";
if(massME_[cand1Label+cand2Label+chargeLabel]){
massME_[cand1Label+cand2Label+chargeLabel]->Fill((cand1[i]->p4()+cand2[j]->p4()).mass());
}
if(massME_[cand1Label+cand2Label+chargeLabel+"maxdzmuontrack"]){
reco::RecoChargedCandidate tk1 = (*cand1[i]);
reco::RecoChargedCandidate tk2 = (*cand2[j]);
// massME_[cand1Label+cand2Label+chargeLabel+"maxdzmuontrack"]->Fill(fabs(tk1->dz(BSPosition_) - tk2->dz(BSPosition_)));
}
}
}
}
| void HLTOniaSource::fillInvariantMass | ( | std::vector< reco::RecoChargedCandidateRef > & | cand1, |
| reco::TrackCollection & | cand2, | ||
| std::string | cand1Label, | ||
| std::string | cand2Label | ||
| ) | [private] |
Definition at line 552 of file HLTOniaSource.cc.
References DeDxDiscriminatorTools::charge(), i, massME_, and p4.
{
typedef reco::TrackCollection::const_iterator cand;
//Loop on collection to calculate invariate mass
for(size_t i = 0 ; i< cand1.size(); i++) {
//Highest PT
if(cand2.begin() != cand2.end() && cand2.begin()->p()>3) {
cand candItr = cand2.begin();
math::PtEtaPhiMLorentzVector bestPtCandLVector(candItr->pt(), candItr->eta(), candItr->phi(), 1.056);
std::string chargeLabel = "same";
//Check relative charge
if(cand1[i]->charge() * candItr->charge() < 0) chargeLabel = "opposite";
if(massME_[cand1Label+cand2Label+chargeLabel+"highestpt"]){
massME_[cand1Label+cand2Label+chargeLabel+"highestpt"]->Fill((cand1[i]->p4()+bestPtCandLVector).mass());
}
if(massME_[cand1Label+cand2Label+chargeLabel+"highestpt"+"maxdzmuontrack"]){
reco::RecoChargedCandidate tk1 = (*cand1[i]);
// massME_[cand1Label+cand2Label+chargeLabel+"highestpt"+"maxdzmuontrack"]->Fill(fabs(tk->dz(BSPosition_) - candItr->dz(BSPosition_)));
}
}
for (cand candIter= cand2.begin(); candIter!=cand2.end(); candIter++) {
if(candIter->p() < 3) continue; //Check if momentum is greater than 3GeV.
math::PtEtaPhiMLorentzVector candLVector(candIter->pt(), candIter->eta(), candIter->phi(), 1.056);
std::string chargeLabel = "same";
//Check relative charge
if(cand1[i]->charge() * candIter->charge() < 0) chargeLabel = "opposite";
if(massME_[cand1Label+cand2Label+chargeLabel]){
massME_[cand1Label+cand2Label+chargeLabel]->Fill((cand1[i]->p4()+candLVector).mass());
}
if(massME_[cand1Label+cand2Label+chargeLabel+"maxdzmuontrack"]){
reco::RecoChargedCandidate tk1 = (*cand1[i]);
// massME_[cand1Label+cand2Label+chargeLabel+"maxdzmuontrack"]->Fill(fabs(tk->dz(BSPosition_) - candIter->dz(BSPosition_)));
}
}
}
}
| void HLTOniaSource::fillOniaTriggerMEs | ( | edm::Handle< reco::RecoChargedCandidateCollection > & | collection, |
| std::string | collectionLabel, | ||
| std::map< std::string, MonitorElement * > & | mapME | ||
| ) | [private] |
Definition at line 437 of file HLTOniaSource.cc.
References reco::LeafCandidate::charge(), runEdmFileComparison::collection, reco::LeafCandidate::eta(), i, edm::HandleBase::isValid(), reco::LeafCandidate::p(), reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().
Referenced by analyze().
{
reco::RecoChargedCandidateCollection myCollection;
if (collection.isValid()) {
myCollection = * collection;
// int nCollection= myCollection.size();
int num = 0;
typedef reco::RecoChargedCandidateCollection::const_iterator cand;
for (cand i=myCollection.begin(); i!=myCollection.end(); i++) {
reco::RecoChargedCandidate tk = (*i);
num++;
//Fill MEs
if(mapME[collectionLabel+"pt"]){ mapME[collectionLabel+"pt"]->Fill(tk.pt()); }
if(mapME[collectionLabel+"p"]) { mapME[collectionLabel+"p"]->Fill(tk.p()); }
if(mapME[collectionLabel+"eta"]) { mapME[collectionLabel+"eta"]->Fill(tk.eta()); }
if(mapME[collectionLabel+"phi"]) { mapME[collectionLabel+"phi"]->Fill(tk.phi()); }
if(mapME[collectionLabel+"etaphi"]){ mapME[collectionLabel+"etaphi"]->Fill(tk.phi(),tk.eta()); }
if(mapME[collectionLabel+"etapt"]){ mapME[collectionLabel+"etapt"]->Fill(tk.pt(),tk.eta()); }
if(mapME[collectionLabel+"charge"]){ mapME[collectionLabel+"charge"]->Fill(tk.charge()); }
// if(mapME[collectionLabel+"dxy"]){ mapME[collectionLabel+"dxy"]->Fill(tk.dxy(BSPosition_)); }
// if(mapME[collectionLabel+"dz"]){ mapME[collectionLabel+"dz"]->Fill(tk.dz(BSPosition_)); }
// if(mapME[collectionLabel+"validhits"]) { mapME[collectionLabel+"validhits"]->Fill(tk.numberOfValidHits()); }
// if(mapME[collectionLabel+"normchi"]){ mapME[collectionLabel+"normchi"]->Fill(tk.normalizedChi2()); }
}
if(mapME[collectionLabel+"nrcand"]){ mapME[collectionLabel+"nrcand"]->Fill(num);}
}
}
| void HLTOniaSource::fillOniaTriggerMEs | ( | std::vector< reco::RecoChargedCandidateRef > & | candidateVector, |
| std::string | collectionLabel, | ||
| std::map< std::string, MonitorElement * > & | mapME | ||
| ) | [private] |
Definition at line 412 of file HLTOniaSource.cc.
References reco::LeafCandidate::charge(), reco::LeafCandidate::eta(), i, reco::LeafCandidate::p(), reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().
{
for (unsigned int i=0; i!=candidateVector.size(); i++) {
reco::RecoChargedCandidate tk = (*candidateVector[i]);
//Fill MEs
if(mapME[collectionLabel+"pt"]){ mapME[collectionLabel+"pt"]->Fill(tk.pt()); }
if(mapME[collectionLabel+"p"]) { mapME[collectionLabel+"p"]->Fill(tk.p()); }
if(mapME[collectionLabel+"eta"]) { mapME[collectionLabel+"eta"]->Fill(tk.eta()); }
if(mapME[collectionLabel+"phi"]) { mapME[collectionLabel+"phi"]->Fill(tk.phi()); }
if(mapME[collectionLabel+"etaphi"]){ mapME[collectionLabel+"etaphi"]->Fill(tk.phi(),tk.eta()); }
if(mapME[collectionLabel+"etapt"]){ mapME[collectionLabel+"etapt"]->Fill(tk.pt(),tk.eta()); }
if(mapME[collectionLabel+"charge"]){ mapME[collectionLabel+"charge"]->Fill(tk.charge()); }
// if(mapME[collectionLabel+"dxy"]){ mapME[collectionLabel+"dxy"]->Fill(tk.dxy(BSPosition_)); }
// if(mapME[collectionLabel+"dz"]){ mapME[collectionLabel+"dz"]->Fill(tk.dz(BSPosition_)); }
// if(mapME[collectionLabel+"validhits"]) { mapME[collectionLabel+"validhits"]->Fill(tk.numberOfValidHits()); }
// if(mapME[collectionLabel+"normchi"]){ mapME[collectionLabel+"normchi"]->Fill(tk.normalizedChi2()); }
}
if(mapME[collectionLabel+"nrcand"]){ mapME[collectionLabel+"nrcand"]->Fill( candidateVector.size());}
}
| void HLTOniaSource::fillOniaTriggerMEs | ( | edm::Handle< reco::TrackCollection > & | collection, |
| std::string | collectionLabel, | ||
| std::map< std::string, MonitorElement * > & | mapME | ||
| ) | [private] |
Definition at line 378 of file HLTOniaSource.cc.
References runEdmFileComparison::collection, and edm::HandleBase::isValid().
{
// cout << "[HLTOniaSource]: fillOniaTriggerMEs " << collectionLabel << endl;
reco::TrackCollection myCollection;
if (collection.isValid()) {
myCollection = * collection;
// int nCollection= myCollection.size();
typedef reco::TrackCollection::const_iterator cand;
int num = 0;
for (cand tk=myCollection.begin(); tk!=myCollection.end(); tk++) {
num++;
//Fill MEs
if(mapME[collectionLabel+"pt"]){ mapME[collectionLabel+"pt"]->Fill(tk->pt()); }
if(mapME[collectionLabel+"p"]) { mapME[collectionLabel+"p"]->Fill(tk->p()); }
if(mapME[collectionLabel+"eta"]) { mapME[collectionLabel+"eta"]->Fill(tk->eta()); }
if(mapME[collectionLabel+"phi"]) { mapME[collectionLabel+"phi"]->Fill(tk->phi()); }
if(mapME[collectionLabel+"etaphi"]){ mapME[collectionLabel+"etaphi"]->Fill(tk->phi(),tk->eta()); }
if(mapME[collectionLabel+"etapt"]){ mapME[collectionLabel+"etapt"]->Fill(tk->pt(),tk->eta()); }
if(mapME[collectionLabel+"charge"]){ mapME[collectionLabel+"charge"]->Fill(tk->charge()); }
// if(mapME[collectionLabel+"dxy"]){ mapME[collectionLabel+"dxy"]->Fill(tk->dxy(BSPosition_)); }
// if(mapME[collectionLabel+"dz"]){ mapME[collectionLabel+"dz"]->Fill(tk->dz(BSPosition_)); }
// if(mapME[collectionLabel+"validhits"]) { mapME[collectionLabel+"validhits"]->Fill(tk->numberOfValidHits()); }
// if(mapME[collectionLabel+"normchi"]){ mapME[collectionLabel+"normchi"]->Fill(tk->normalizedChi2()); }
}
if(mapME[collectionLabel+"nrcand"]){ mapME[collectionLabel+"nrcand"]->Fill(num);}
}
}
edm::InputTag HLTOniaSource::beamSpotTag_ [private] |
Definition at line 64 of file HLTOniaSource.h.
Referenced by analyze(), and HLTOniaSource().
Definition at line 41 of file HLTOniaSource.h.
Referenced by analyze().
DQMStore* HLTOniaSource::dbe_ [private] |
Definition at line 66 of file HLTOniaSource.h.
Referenced by beginJob(), beginRun(), bookOniaTriggerInvariantMassMEs(), bookOniaTriggerMEs(), and ~HLTOniaSource().
bool HLTOniaSource::hltConfigInit_ [private] |
Definition at line 73 of file HLTOniaSource.h.
Referenced by analyze(), and beginRun().
std::string HLTOniaSource::hltProcessName_ [private] |
Definition at line 56 of file HLTOniaSource.h.
Referenced by beginRun(), and HLTOniaSource().
std::map<std::string, MonitorElement *> HLTOniaSource::massME_ [private] |
Definition at line 72 of file HLTOniaSource.h.
Referenced by beginRun(), bookOniaTriggerInvariantMassMEs(), and fillInvariantMass().
std::map<std::string, MonitorElement *> HLTOniaSource::muonME_ [private] |
Definition at line 70 of file HLTOniaSource.h.
Referenced by analyze(), and beginRun().
std::vector<edm::InputTag> HLTOniaSource::oniaMuonTag_ [private] |
Definition at line 59 of file HLTOniaSource.h.
Referenced by analyze(), beginRun(), and HLTOniaSource().
std::map<std::string, MonitorElement *> HLTOniaSource::pixelAfterFilterME_ [private] |
Definition at line 67 of file HLTOniaSource.h.
Referenced by analyze(), and beginRun().
std::map<std::string, MonitorElement *> HLTOniaSource::pixelME_ [private] |
Definition at line 69 of file HLTOniaSource.h.
Referenced by analyze(), and beginRun().
edm::InputTag HLTOniaSource::pixelTag_ [private] |
Definition at line 63 of file HLTOniaSource.h.
Referenced by analyze(), beginRun(), and HLTOniaSource().
std::vector<edm::InputTag> HLTOniaSource::pixelTagsAfterFilter_ [private] |
Definition at line 60 of file HLTOniaSource.h.
Referenced by analyze(), beginRun(), and HLTOniaSource().
std::string HLTOniaSource::subsystemFolder_ [private] |
Definition at line 55 of file HLTOniaSource.h.
Referenced by beginRun(), and HLTOniaSource().
std::map<std::string, MonitorElement *> HLTOniaSource::trackAfterFilterME_ [private] |
Definition at line 68 of file HLTOniaSource.h.
Referenced by analyze(), and beginRun().
std::map<std::string, MonitorElement *> HLTOniaSource::trackME_ [private] |
Definition at line 71 of file HLTOniaSource.h.
Referenced by analyze(), and beginRun().
edm::InputTag HLTOniaSource::trackTag_ [private] |
Definition at line 65 of file HLTOniaSource.h.
Referenced by analyze(), beginRun(), and HLTOniaSource().
std::vector<edm::InputTag> HLTOniaSource::trackTagsAfterFilter_ [private] |
Definition at line 61 of file HLTOniaSource.h.
Referenced by analyze(), beginRun(), and HLTOniaSource().
std::vector< std::string> HLTOniaSource::triggerPath_ [private] |
Definition at line 58 of file HLTOniaSource.h.
Referenced by checkHLTConfiguration(), and HLTOniaSource().
Definition at line 62 of file HLTOniaSource.h.
Referenced by analyze(), and HLTOniaSource().
1.7.3