Origin/add pt cut for pfb tagging

plugins/BTagHLTAnalyzer.cc

@@ -460,9 +460,14 @@ void BTagHLTAnalyzerT<IPTI,VTX>::analyze(const edm::Event& iEvent, const edm::Ev
   if ( !foundNames ) edm::LogError("TriggerNamesNotFound") << "Could not get trigger names!";
   if ( trigRes->size() != triggerList.size() ) edm::LogError("TriggerPathLengthMismatch") << "Length of names and paths not the same: "
 											  << triggerList.size() << "," << trigRes->size() ;
-
-  bool passTrig = processTrig(trigRes, triggerList);
-  if(!passTrig) return;
+
+  //debug: check if passes trigger
+  //bool passTrig = processTrig(trigRes, triggerList);
+  if(!passTrig){
+    //cout << "Error, failed TRIG"<<endl;
+
+    return;
+    }
 
   //------------------------------------------------------
   // Jet info

python/bTagAnalyzerHLT_cfi.py

@@ -2,21 +2,25 @@
 from RecoBTag.PerformanceMeasurements.bTagAnalyzerCommon_cff import *
 from RecoBTag.PerformanceMeasurements.variables_cfi import *
 from RecoBTag.PerformanceMeasurements.varGroups_cfi import *
-bTagHLTAnalyzer = cms.EDAnalyzer("BTagHLTAnalyzer",
-                                 bTagAnalyzerCommon,
-                                 variableSet,
-                                 groupSet,
-                                 HLTprimaryVertexColl = cms.InputTag('hltVerticesPFFilter'),
-                                 CaloJets             = cms.InputTag('hltAK4CaloJetsCorrectedIDPassed'),
-                                 CaloJetTags          = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsInfosCalo'),
-                                 CaloSVs              = cms.InputTag('hltInclusiveSecondaryVertexFinderTagInfos'),
-                                 CaloJetCSVTags       = cms.InputTag('hltCombinedSecondaryVertexBJetTagsCalo'),
-                                 CaloJetDeepCSVTags   = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsCalo:probb'),
-                                 PFJets               = cms.InputTag('hltAK4PFJetsLooseIDCorrected'),
-                                 PFJetTags            = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsInfos'),
-                                 PFSVs                = cms.InputTag('hltDeepSecondaryVertexTagInfosPF'),
-                                 PFJetCSVTags         = cms.InputTag('hltCombinedSecondaryVertexBJetTagsPF'),
-                                 PFJetDeepCSVTags     = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsPF:probb'),
-                                 HLTTriggerPathNames  = cms.vstring("HLT_Mu12_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_*",
-                                                                    "HLT_Mu23_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_*"),
-)
+
+def bTagAnalyzer_func(trigPaths=[]):
+    bTagHLTAnalyzer = cms.EDAnalyzer("BTagHLTAnalyzer",
+                                     bTagAnalyzerCommon,
+                                     variableSet,
+                                     groupSet,
+                                     HLTprimaryVertexColl = cms.InputTag('hltVerticesPFFilter'),
+                                     CaloJets             = cms.InputTag('hltAK4CaloJetsCorrectedIDPassed'),
+                                     CaloJetTags          = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsInfosCalo'),
+                                     CaloSVs              = cms.InputTag('hltInclusiveSecondaryVertexFinderTagInfos'),
+                                     CaloJetCSVTags       = cms.InputTag('hltCombinedSecondaryVertexBJetTagsCalo'),
+                                     CaloJetDeepCSVTags   = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsCalo:probb'),
+                                     PFJets               = cms.InputTag('hltAK4PFJetsLooseIDCorrected'),
+                                     PFJetTags            = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsInfos'),
+                                     PFSVs                = cms.InputTag('hltDeepSecondaryVertexTagInfosPF'),
+                                     PFJetCSVTags         = cms.InputTag('hltCombinedSecondaryVertexBJetTagsPF'),
+                                     PFJetDeepCSVTags     = cms.InputTag('hltDeepCombinedSecondaryVertexBJetTagsPF:probb'),
+                                     HLTTriggerPathNames  = cms.vstring(trigPaths),
+    )
+    print(trigPaths)
+
+    return bTagHLTAnalyzer

python/customise.py

@@ -112,7 +112,7 @@ def customize_HLTDeepCSVPF(process,ptVal=0.9):
                            seedMin3DIPSignificance = cms.double(1.2),
                            seedMin3DIPValue = cms.double(0.005)
                          ),
-        fitterRatio = cms.double(0.25),
+                         fitterRatio = cms.double(0.25),
                          fitterSigmacut = cms.double(3.0),
                          fitterTini = cms.double(256.0),
                          maxNTracks = cms.uint32(30),
@@ -146,5 +146,83 @@ def customize_HLTDeepCSVPF(process,ptVal=0.9):
 
   return process
 
+#Hello world!
+def customize_CaloJet(process,ptVal=0.9):
+
+  ptStr = str(ptVal).replace(".","p")
+
+  setattr(process, "hltImpactParameterTagInfosPtCut" + ptStr,
+          cms.EDProducer("TrackIPProducer", 
+                computeGhostTrack = cms.bool(True),
+                computeProbabilities = cms.bool(True),
+                ghostTrackPriorDeltaR = cms.double(0.03),
+                jetDirectionUsingGhostTrack = cms.bool(False),
+                jetDirectionUsingTracks = cms.bool(False),
+                jetTracks = cms.InputTag("hltFastPixelBLifetimeL3Associator"), 
+                maximumChiSquared = cms.double(5.0),
+                maximumLongitudinalImpactParameter = cms.double(17.0),
+                maximumTransverseImpactParameter = cms.double(0.2),
+                minimumNumberOfHits = cms.int32(3), 
+                minimumNumberOfPixelHits = cms.int32(2),
+                minimumTransverseMomentum = cms.double(ptVal),    
+                primaryVertex = cms.InputTag("hltVerticesL3","WithBS"), 
+                useTrackQuality = cms.bool(False)
+          )
+  )
+
+  process.hltInclusiveSecondaryVertexFinderTagInfos.trackIPTagInfos = "hltImpactParameterTagInfosPtCut" + ptStr 
+
+  process.HLTBtagDeepCSVSequenceL3.replace(
+    process.hltImpactParameterTagInfos,
+    getattr(process, "hltImpactParameterTagInfosPtCut" + ptStr)
+    )
+
+  setattr(process, "hltInclusiveVertexFinderPtCut" + ptStr,
+          cms.EDProducer("InclusiveVertexFinder",
+                         beamSpot = cms.InputTag("hltOnlineBeamSpot"),
+                         clusterizer = cms.PSet(
+                                        clusterMaxDistance = cms.double(0.05),
+                                        clusterMaxSignificance = cms.double(4.5),
+                                        clusterMinAngleCosine = cms.double(0.5),
+                                        distanceRatio = cms.double(20.0),
+                                        seedMax3DIPSignificance = cms.double(9999.0),
+                                        seedMax3DIPValue = cms.double(9999.0),
+                                        seedMin3DIPSignificance = cms.double(1.2),
+                                        seedMin3DIPValue = cms.double(0.005)
+                                        ),
+                         fitterRatio = cms.double(0.25),
+                         fitterSigmacut = cms.double(3.0),
+                         fitterTini = cms.double(256.0),
+                         maxNTracks = cms.uint32(30),
+                         maximumLongitudinalImpactParameter = cms.double(0.3),
+                         maximumTimeSignificance = cms.double(3.0),
+                         minHits = cms.uint32(8),
+                         minPt = cms.double(ptVal), #original 0.8
+                         primaryVertices = cms.InputTag("hltVerticesL3"),
+                         tracks = cms.InputTag("hltMergedTracksForBTag"),
+                         useDirectVertexFitter = cms.bool(True),
+                         useVertexReco = cms.bool(True),
+                         vertexMinAngleCosine = cms.double(0.95),
+                         vertexMinDLen2DSig = cms.double(2.5),
+                         vertexMinDLenSig = cms.double(0.5),
+                         vertexReco = cms.PSet(
+                             finder = cms.string('avr'),
+                             primcut = cms.double(1.0),
+                             seccut = cms.double(3.0),
+                             smoothing = cms.bool(True)
+                         )
+          )
+  )
+
+  process.hltInclusiveSecondaryVertices.secondaryVertices = "hltInclusiveVertexFinderPtCut" + ptStr
+
+  process.HLTBtagDeepCSVSequenceL3.replace(
+    process.hltInclusiveVertexFinder,
+    getattr(process,"hltInclusiveVertexFinderPtCut"+ptStr)
+  )
+
+
+  return process
 
 
+

test/runHLTBTagAnalyzer_cfg.py

@@ -52,6 +52,15 @@
     "Global Pt Seed  (default is 0.4)"
 )
 
+options.register('trigNames',  'HLT_Mu12_trkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_*,HLT_Mu23_trkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_*',
+    VarParsing.multiplicity.list,
+    VarParsing.varType.string,
+    "Trigger Names (defaults are HLT_Mu12_trkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_* and HLT_Mu23_trkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_*)"
+)
+
+    #HLT_ZeroBias_Beamspot_v*
+
+
 options.register('globalTag', '110X_mcRun3_2021_realistic_v6',
     VarParsing.multiplicity.singleton,
     VarParsing.varType.string,
@@ -8543,7 +8552,7 @@
 )
 
 
-process.hltDeepCombinedSecondaryVertexBJetTagsInfosCalo = cms.EDProducer("TrackDeepNNTagInfoProducer",
+process.hltDeepCombinedSecondaryVertexBJetTagsInfosCalo = cms.EDProducer("TrackDeepNNTagInfoProducer", #check this
     computer = cms.PSet(
         SoftLeptonFlip = cms.bool(False),
         charmCut = cms.double(1.5),
@@ -8689,7 +8698,7 @@
     extSVCollection = cms.InputTag("hltDeepInclusiveMergedVerticesPF"),
     extSVDeltaRToJet = cms.double(0.3),
     minimumTrackWeight = cms.double(0.5),
-    trackIPTagInfos = cms.InputTag("hltDeepBLifetimeTagInfosPF"),
+    trackIPTagInfos = cms.InputTag("hltDeepBLifetimeTagInfosPF"), #FIXME
     trackSelection = cms.PSet(
         a_dR = cms.double(-0.001053),
         a_pT = cms.double(0.005263),
@@ -22741,7 +22750,7 @@
 
 
 process.noFilter_CaloDeepCSV = cms.Path(process.HLTBeginSequence+process.hltPrenoFilterCaloDeepCSV+process.HLTAK4CaloJetsSequence+process.HLTBtagDeepCSVSequenceL3+process.HLTEndSequence)
-
+#HELLO CaloDeepCSV
 
 
 
@@ -23038,7 +23047,9 @@
 
 #-------------------------------------
 from RecoBTag.PerformanceMeasurements.BTagHLTAnalyzer_cff import *
-process.btagana = bTagHLTAnalyzer.clone()
+bta = bTagAnalyzer_func(trigPaths=options.trigNames)
+process.btagana = bta.clone()
+
 
 #------------------
 #Handle groups
@@ -23090,6 +23101,8 @@
 from RecoBTag.PerformanceMeasurements.customise import customize_HLTDeepCSVPF
 process = customize_HLTDeepCSVPF(process,options.trackPtSeed)
 
+from RecoBTag.PerformanceMeasurements.customise import customize_CaloJet 
+process = customize_CaloJet(process,options.trackPtSeed)
 
 ## Define analyzer sequence
 process.analyzerSeq = cms.Sequence( )

test/submit_allHLTMC.sh

@@ -1,10 +1,87 @@
 #!/bin/bash
 
+cmsenv
+#voms-proxy-init -rfc -voms cms -valid 192:00
+source /cvmfs/cms.cern.ch/crab3/crab.sh
+
 echo "!!!! WARNING: Submitting for MC!!!"
 python submit_allHLT.py \
   runHLTBTagAnalyzer_cfg.py \
   -f CRAB/tosubmit_MC.txt \
   -s T3_US_FNALLPC \
-  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=2.0 \
-  -o /store/user/johnda/BTagNTuples/Run3/ \
-  -v crab_projects_run3_trackPt2p0
+  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=0.4 trigNames="HLT_ZeroBias_Beamspot_v*"\
+  -o /store/user/wterrill/BTagNTuples/Run3/ \
+  -v crab_projects_run3_trackPt0p4_trigChange
+
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=0.9 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt0p9_trigChange
+#
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=1.5 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt1p5_trigChange
+#
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=2.0 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt2p0_trigChange
+#
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=2.5 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt2p5_trigChange
+#
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=5.0 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt5p0_trigChange
+#
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=7.5 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt7p5_trigChange
+#
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=10.0 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt10p0_trigChange
+#
+#echo "!!!! WARNING: Submitting for MC!!!"
+#python submit_allHLT.py \
+#  runHLTBTagAnalyzer_cfg.py \
+#  -f CRAB/tosubmit_MC.txt \
+#  -s T3_US_FNALLPC \
+#  -p groups="HLTEventInfo,HLTJetInfo,HLTTagVar,HLTJetTrack,HLTJetSV,HLTCSVTagVar" runOnData=False globalTag="110X_mcRun3_2021_realistic_v6" trackPtSeed=20.0 trigNames="HLT_ZeroBias_Beamspot_v*"\
+#  -o /store/user/wterrill/BTagNTuples/Run3/ \
+#  -v crab_projects_run3_trackPt20p0_trigChange
+#