Doublet tuning and new vertex algo

Calibration/HcalAlCaRecoProducers/plugins/AlCaIsoTracksFilter.cc

@@ -60,8 +60,8 @@
 
 namespace AlCaIsoTracks {
   struct Counters {
-    Counters() : nAll_(0), nGood_(0), nRange_(0) {}
-    mutable std::atomic<unsigned int> nAll_, nGood_, nRange_;
+    Counters() : nAll_(0), nGood_(0), nRange_(0), nHigh_(0) {}
+    mutable std::atomic<unsigned int> nAll_, nGood_, nRange_, nHigh_;
   };
 }
 
@@ -95,12 +95,12 @@ class AlCaIsoTracksFilter : public edm::stream::EDFilter<edm::GlobalCache<AlCaIs
   const double                   eIsolate_;
   const double                   hitEthrEB_, hitEthrEE0_, hitEthrEE1_;
   const double                   hitEthrEE2_, hitEthrEE3_;
-  const double                   pTrackLow_, pTrackHigh_;
-  const int                      preScale_;
+  const double                   pTrackLow_, pTrackHigh_, pTrackH_;
+  const int                      preScale_, preScaleH_;
   const std::string              theTrackQuality_;
   spr::trackSelectionParameters  selectionParameter_;
   double                         a_charIsoR_;
-  unsigned int                   nRun_, nAll_, nGood_, nRange_;
+  unsigned int                   nRun_, nAll_, nGood_, nRange_, nHigh_;
   edm::EDGetTokenT<trigger::TriggerEvent>  tok_trigEvt_;
   edm::EDGetTokenT<edm::TriggerResults>    tok_trigRes_;
   edm::EDGetTokenT<reco::TrackCollection>  tok_genTrack_;
@@ -146,9 +146,11 @@ AlCaIsoTracksFilter::AlCaIsoTracksFilter(const edm::ParameterSet& iConfig, const
   hitEthrEE3_(iConfig.getParameter<double>("EEHitEnergyThreshold3") ),
   pTrackLow_(iConfig.getParameter<double>("momentumRangeLow")),
   pTrackHigh_(iConfig.getParameter<double>("momentumRangeHigh")),
+  pTrackH_(iConfig.getParameter<double>("momentumHigh")),
   preScale_(iConfig.getParameter<int>("preScaleFactor")),
+  preScaleH_(iConfig.getParameter<int>("preScaleHigh")),
   theTrackQuality_(iConfig.getParameter<std::string>("trackQuality")),
-  nRun_(0), nAll_(0), nGood_(0), nRange_(0) {
+  nRun_(0), nAll_(0), nGood_(0), nRange_(0), nHigh_(0) {
   //now do what ever initialization is needed
   const double isolationRadius(28.9);
   // Different isolation cuts are described in DN-2016/029
@@ -200,7 +202,8 @@ AlCaIsoTracksFilter::AlCaIsoTracksFilter(const edm::ParameterSet& iConfig, const
 			       <<"\t eIsolate_ "       << eIsolate_ << "\n"
 			       <<"\t Precale factor "  << preScale_
 			       <<"\t in momentum range " << pTrackLow_
-			       <<":" << pTrackHigh_;
+			       <<":" << pTrackHigh_ << " and prescale factor "
+			       << preScaleH_ << " for p > " << pTrackH_;
   for (unsigned int k=0; k<trigNames_.size(); ++k)
     edm::LogInfo("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
 } // AlCaIsoTracksFilter::AlCaIsoTracksFilter  constructor
@@ -323,7 +326,7 @@ bool AlCaIsoTracksFilter::filter(edm::Event& iEvent, edm::EventSetup const& iSet
 			 trkCaloDirections, false);
 
       std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
-      unsigned int nTracks(0), nselTracks(0), ntrin(0), ntrout(0);
+      unsigned int nTracks(0), nselTracks(0), ntrin(0), ntrout(0), ntrH(0);
       for (trkDetItr = trkCaloDirections.begin(),nTracks=0; 
 	   trkDetItr != trkCaloDirections.end(); trkDetItr++,nTracks++) {
 	const reco::Track* pTrack = &(*(trkDetItr->trkItr));
@@ -381,6 +384,7 @@ bool AlCaIsoTracksFilter::filter(edm::Event& iEvent, edm::EventSetup const& iSet
 					   << ":" << eIsolation;
 	  if (t_p>pTrackMin_ && eMipDR<eEcalMax_ && hmaxNearP<eIsolation) {
 	    if (t_p > pTrackLow_ && t_p < pTrackHigh_) ntrin++;
+	    else if (t_p > pTrackH_)                   ntrH++;
 	    else                                       ntrout++;
 	  }
 	}
@@ -391,6 +395,11 @@ bool AlCaIsoTracksFilter::filter(edm::Event& iEvent, edm::EventSetup const& iSet
 	if      (preScale_ <= 1)         accept = true;
 	else if (nRange_%preScale_ == 1) accept = true;
       }
+      if (!accept && ntrH > 0) {
+	++nHigh_;
+	if      (preScaleH_ <= 1)        accept = true;
+	else if (nHigh_%preScaleH_ == 1) accept = true;
+      }
     }
   }
   // Step 4:  Return the acceptance flag
@@ -404,12 +413,14 @@ void AlCaIsoTracksFilter::endStream() {
   globalCache()->nAll_   += nAll_;
   globalCache()->nGood_  += nGood_;
   globalCache()->nRange_ += nRange_;
+  globalCache()->nHigh_  += nHigh_;
 }
 
 void AlCaIsoTracksFilter::globalEndJob(const AlCaIsoTracks::Counters* count) {
   edm::LogVerbatim("HcalIsoTrack") << "Selects " << count->nGood_ << " in " 
 				   << count->nAll_ << " events and with "
-				   << count->nRange_ << " events in the p-range";
+				   << count->nRange_ << " events in the p-range"
+				   << count->nHigh_ << " events with high p";
 }
 
 
@@ -456,7 +467,7 @@ void AlCaIsoTracksFilter::fillDescriptions(edm::ConfigurationDescriptions& descr
   desc.add<double>("minimumTrackP",20.0);
   // signal zone in ECAL and MIP energy cutoff
   desc.add<double>("coneRadiusMIP",14.0);
-  desc.add<double>("maximumEcalEnergy",2.0);
+  desc.add<double>("maximumEcalEnergy",100.0);
   // following 3 parameters are for isolation cuts and described in the code
   desc.add<double>("maxTrackP",8.0);
   desc.add<double>("slopeTrackP",0.05090504066);
@@ -470,7 +481,9 @@ void AlCaIsoTracksFilter::fillDescriptions(edm::ConfigurationDescriptions& descr
   // Prescale events only containing isolated tracks in the range
   desc.add<double>("momentumRangeLow",20.0);
   desc.add<double>("momentumRangeHigh",40.0);
-  desc.add<int>("preScaleFactor",2);
+  desc.add<int>("preScaleFactor",10);
+  desc.add<double>("momentumHigh",60.0);
+  desc.add<int>("preScaleHigh",2);
   descriptions.add("alcaIsoTracksFilter",desc);
 }
 

Calibration/HcalAlCaRecoProducers/python/alcaisotrk_cfi.py

@@ -38,6 +38,5 @@
 # Prescale events only containing isolated tracks in the range
                          MomentumRangeLow  = cms.double(20.0),
                          MomentumRangeHigh = cms.double(40.0),
-                         PreScaleFactor    = cms.int32(2)
+                         PreScaleFactor    = cms.int32(10)
 )
-

CondTools/RunInfo/python/LHCInfoPopConAnalyzerStartFill.py

@@ -70,13 +70,14 @@
                                record = cms.string('LHCInfoRcd'),
                                name = cms.untracked.string('LHCInfo'),
                                Source = cms.PSet(fill = cms.untracked.uint32(6417),
-                                   startTime = cms.untracked.string('2018-05-31 00:00:00.000'),
-                                   #endTime = cms.untracked.string('2018-05-31 20:13:59.000'),
+                                   startTime = cms.untracked.string('2018-04-01 00:00:00.000'),
+                                   #endTime = cms.untracked.string('2018-03-25 05:00:00.000'),
                                    samplingInterval = cms.untracked.uint32( 600 ),
                                    endFill = cms.untracked.bool(False),
                                    connectionString = cms.untracked.string("oracle://cms_orcon_adg/CMS_RUNTIME_LOGGER"),
                                    ecalConnectionString = cms.untracked.string("oracle://cms_orcon_adg/CMS_DCS_ENV_PVSS_COND"),
                                    DIPSchema = cms.untracked.string("CMS_BEAM_COND"),
+                                   #authenticationPath =  cms.untracked.string("."),
                                    debug=cms.untracked.bool(False)
                                                  ),
                                loggingOn = cms.untracked.bool(True),

CondTools/RunInfo/src/LHCInfoPopConSourceHandler.cc

@@ -830,29 +830,25 @@ void LHCInfoPopConSourceHandler::getNewObjects() {
 			   << " ); from " << m_name << "::getNewObjects";
   }
 
-  cond::Time_t lastIov = tagInfo().lastInterval.first; 
-  if( lastIov == 0 ){
+  cond::Time_t lastSince = tagInfo().lastInterval.first; 
+  if( lastSince == 0 ){
     // for a new or empty tag, an empty payload should be added on top with since=1
     addEmptyPayload( 1 );
   } else {
     edm::LogInfo( m_name ) << "The last Iov in tag " << tagInfo().name 
-			   << " valid since " << lastIov
+			   << " valid since " << lastSince
 			   << "from " << m_name << "::getNewObjects";
   }
 
-  cond::Time_t targetIov = cond::time::from_boost( boost::posix_time::second_clock::local_time() );  
-  cond::Time_t endIov = cond::time::MAX_VAL;
-  //if( m_endFill ){
-    if( !m_startTime.is_not_a_date_time() ){
-      targetIov = cond::time::from_boost(m_startTime);
-    }
-    if( lastIov > targetIov ) targetIov = lastIov;
+  boost::posix_time::ptime executionTime = boost::posix_time::second_clock::local_time();
+  cond::Time_t targetSince = 0;  
+  cond::Time_t endIov = cond::time::from_boost( executionTime );  
+  if( !m_startTime.is_not_a_date_time() ){
+    targetSince = cond::time::from_boost(m_startTime);
+  }
+  if( lastSince > targetSince ) targetSince = lastSince;
 
-    if( !m_endTime.is_not_a_date_time() ){
-      endIov = cond::time::from_boost(m_endTime);
-    }
-    //}
-  edm::LogInfo(m_name) <<"Starting sampling at "<<boost::posix_time::to_simple_string(cond::time::to_boost(targetIov));
+  edm::LogInfo(m_name) <<"Starting sampling at "<<boost::posix_time::to_simple_string(cond::time::to_boost(targetSince));
 
   //retrieve the data from the relational database source
   cond::persistency::ConnectionPool connection;
@@ -875,36 +871,36 @@ void LHCInfoPopConSourceHandler::getNewObjects() {
     session3.transaction().commit();
   }
 
+  bool iovAdded = false;
   while( true ){
-    if( targetIov >= endIov ){
+    if( targetSince >= endIov ){
       edm::LogInfo( m_name ) <<"Sampling ended at the time "<<boost::posix_time::to_simple_string(cond::time::to_boost( endIov ));
       break;
     }
     bool updateEcal=false;
-    boost::posix_time::ptime targetTime = cond::time::to_boost( targetIov );
+    boost::posix_time::ptime targetTime = cond::time::to_boost( targetSince );
     boost::posix_time::ptime startSampleTime; 
     boost::posix_time::ptime endSampleTime;
     if( !m_endFill and m_prevPayload->fillNumber() and m_prevPayload->endTime()==0ULL){
       // execute the query for the current fill 
       session.transaction().start(true);
+      edm::LogInfo( m_name ) <<"Searching started fill #"<<m_prevPayload->fillNumber();
       bool foundFill = getFillData( session, m_prevPayload->fillNumber() );
       session.transaction().commit();
       if(!foundFill ){
 	edm::LogError( m_name )<<"Could not find fill #"<<m_prevPayload->fillNumber();
 	break;
       }
       updateEcal = true;
-      startSampleTime = cond::time::to_boost(lastIov);
-      // if fill is going, no further query to do...
-      if(m_fillPayload->endTime()==0ULL) endIov = targetIov;
+      startSampleTime = cond::time::to_boost(lastSince);
     } else {
       session.transaction().start(true);
       edm::LogInfo( m_name ) <<"Searching new fill after "<<boost::posix_time::to_simple_string(targetTime);
       bool foundFill = getNextFillData( session, targetTime, m_endFill );
       session.transaction().commit();
       if ( !foundFill ){
 	edm::LogInfo( m_name )<<"No fill found - END of job.";
-	addEmptyPayload( targetIov );
+	if( iovAdded ) addEmptyPayload( targetSince );
 	break;
       }
       startSampleTime = cond::time::to_boost(m_fillPayload->createTime());
@@ -914,11 +910,12 @@ void LHCInfoPopConSourceHandler::getNewObjects() {
     unsigned short lhcFill = m_fillPayload->fillNumber();
     if( endFillTime == 0ULL ){
       edm::LogInfo( m_name ) <<"Found ongoing fill "<<lhcFill<<" created at "<<cond::time::to_boost(startFillTime);
-      endSampleTime = cond::time::to_boost(targetIov);
+      endSampleTime = executionTime;
+      targetSince = endIov;
     } else {
       edm::LogInfo( m_name ) <<"Found fill "<<lhcFill<<" created at "<<cond::time::to_boost(startFillTime)<<" ending at "<<cond::time::to_boost(endFillTime);
       endSampleTime = cond::time::to_boost(endFillTime);
-      targetIov = endFillTime;
+      targetSince = endFillTime;
     }
 
     session.transaction().start(true);
@@ -937,7 +934,9 @@ void LHCInfoPopConSourceHandler::getNewObjects() {
     size_t niovs = LHCInfoImpl::transferPayloads( m_tmpBuffer, m_payloadBuffer, m_to_transfer, m_prevPayload );
     edm::LogInfo( m_name ) <<"Added "<<niovs<<" iovs within the Fill time";
     m_tmpBuffer.clear();
-    if(m_prevPayload->fillNumber() and m_prevPayload->endTime()!=0ULL) addEmptyPayload( m_fillPayload->endTime() );
+    iovAdded = true;
+    //if(m_prevPayload->fillNumber() and m_prevPayload->endTime()!=0ULL) addEmptyPayload( m_fillPayload->endTime() );
+    if(m_prevPayload->fillNumber() and m_fillPayload->endTime()!=0ULL) addEmptyPayload( m_fillPayload->endTime() );
   }
 }
 

Configuration/AlCa/python/autoCond.py

@@ -24,13 +24,13 @@
     # GlobalTag for MC production (p-Pb collisions) with realistic alignment and calibrations for Run2
     'run2_mc_pa'        :   '102X_mcRun2_pA_v2',
     # GlobalTag for Run1 data reprocessing
-    'run1_data'         :   '102X_dataRun2_v3',
+    'run1_data'         :   '102X_dataRun2_v4',
     # GlobalTag for Run2 data reprocessing
-    'run2_data'         :   '102X_dataRun2_v3',
+    'run2_data'         :   '102X_dataRun2_v4',
     # GlobalTag for Run2 data relvals: allows customization to run with fixed L1 menu
-    'run2_data_relval'  :   '102X_dataRun2_relval_v4',
+    'run2_data_relval'  :   '102X_dataRun2_relval_v5',
     # GlobalTag for Run2 data 2016H relvals only: Prompt Conditions + fixed L1 menu (to be removed)
-    'run2_data_promptlike' : '102X_dataRun2_PromptLike_v4',
+    'run2_data_promptlike' : '102X_dataRun2_PromptLike_v7',
     # GlobalTag for Run1 HLT: it points to the online GT
     'run1_hlt'          :   '101X_dataRun2_HLT_frozen_v6',
     # GlobalTag for Run2 HLT: it points to the online GT
@@ -40,19 +40,19 @@
     # GlobalTag for Run2 HLT for HI: it points to the online GT
     'run2_hlt_hi'       :   '101X_dataRun2_HLTHI_frozen_v7',
     # GlobalTag for MC production with perfectly aligned and calibrated detector for Phase1 2017 (and 0,0,~0-centred beamspot)
-    'phase1_2017_design'       :  '102X_mc2017_design_IdealBS_v3',
+    'phase1_2017_design'       :  '102X_mc2017_design_IdealBS_v4',
     # GlobalTag for MC production with realistic conditions for Phase1 2017 detector
-    'phase1_2017_realistic'    : '102X_mc2017_realistic_v3',
+    'phase1_2017_realistic'    : '102X_mc2017_realistic_v4',
     # GlobalTag for MC production (cosmics) with realistic alignment and calibrations for Phase1 2017 detector, Strip tracker in DECO mode
-    'phase1_2017_cosmics'      : '102X_mc2017cosmics_realistic_deco_v3',
+    'phase1_2017_cosmics'      : '102X_mc2017cosmics_realistic_deco_v4',
     # GlobalTag for MC production (cosmics) with realistic alignment and calibrations for Phase1 2017 detector, Strip tracker in PEAK mode
-    'phase1_2017_cosmics_peak' : '102X_mc2017cosmics_realistic_peak_v3',
+    'phase1_2017_cosmics_peak' : '102X_mc2017cosmics_realistic_peak_v4',
     # GlobalTag for MC production with perfectly aligned and calibrated detector for full Phase1 2018 (and 0,0,0-centred beamspot)
-    'phase1_2018_design'       : '102X_upgrade2018_design_v6',
+    'phase1_2018_design'       : '102X_upgrade2018_design_v7',
     # GlobalTag for MC production with realistic conditions for full Phase1 2018 detector
-    'phase1_2018_realistic'    : '102X_upgrade2018_realistic_v11',
+    'phase1_2018_realistic'    : '102X_upgrade2018_realistic_v12',
     # GlobalTag for MC production (cosmics) with realistic conditions for full Phase1 2018 detector,  Strip tracker in DECO mode
-    'phase1_2018_cosmics'      :   '102X_upgrade2018cosmics_realistic_deco_v10',
+    'phase1_2018_cosmics'      :   '102X_upgrade2018cosmics_realistic_deco_v11',
     # GlobalTag for MC production with perfectly aligned and calibrated detector for Phase1 2019
     'phase1_2019_design'       : '102X_postLS2_design_v5', # GT containing design conditions for postLS2
     # GlobalTag for MC production with perfectly aligned and calibrated detector for Phase1 2019

Configuration/Eras/python/Modifier_pf_badHcalMitigation_cff.py

@@ -0,0 +1,3 @@
+import FWCore.ParameterSet.Config as cms
+
+pf_badHcalMitigation =  cms.Modifier()

Configuration/Eras/python/Modifier_run2_nanoAOD_94X2016_cff.py

@@ -0,0 +1,3 @@
+import FWCore.ParameterSet.Config as cms
+
+run2_nanoAOD_94X2016 = cms.Modifier()

Configuration/Generator/python/PSweightsPythia/PythiaPSweightsSettings_cfi.py

@@ -5,11 +5,45 @@
         'UncertaintyBands:doVariations = on',
 # 3 sets of variations for ISR&FSR up/down
 # Reduced sqrt(2)/(1/sqrt(2)), Default 2/0.5 and Conservative 4/0.25 variations
+# 32 decorrelated variations of muR and non-singular terms (cNS) for each branching type
         'UncertaintyBands:List = {\
 isrRedHi isr:muRfac=0.707,fsrRedHi fsr:muRfac=0.707,isrRedLo isr:muRfac=1.414,fsrRedLo fsr:muRfac=1.414,\
-isrDefHi isr:muRfac=0.5, fsrDefHi fsr:muRfac=0.5,isrDefLo isr:muRfac=2.0,fsrDefLo fsr:muRfac=2.0,\
-isrConHi isr:muRfac=0.25, fsrConHi fsr:muRfac=0.25,isrConLo isr:muRfac=4.0,fsrConLo fsr:muRfac=4.0}',
-
+isrDefHi isr:muRfac=0.5,fsrDefHi fsr:muRfac=0.5,isrDefLo isr:muRfac=2.0,fsrDefLo fsr:muRfac=2.0,\
+isrConHi isr:muRfac=0.25,fsrConHi fsr:muRfac=0.25,isrConLo isr:muRfac=4.0,fsrConLo fsr:muRfac=4.0,\
+fsr_G2GG_muR_dn fsr:G2GG:muRfac=0.5,\
+fsr_G2GG_muR_up fsr:G2GG:muRfac=2.0,\
+fsr_G2QQ_muR_dn fsr:G2QQ:muRfac=0.5,\
+fsr_G2QQ_muR_up fsr:G2QQ:muRfac=2.0,\
+fsr_Q2QG_muR_dn fsr:Q2QG:muRfac=0.5,\
+fsr_Q2QG_muR_up fsr:Q2QG:muRfac=2.0,\
+fsr_X2XG_muR_dn fsr:X2XG:muRfac=0.5,\
+fsr_X2XG_muR_up fsr:X2XG:muRfac=2.0,\
+fsr_G2GG_cNS_dn fsr:G2GG:cNS=-2.0,\
+fsr_G2GG_cNS_up fsr:G2GG:cNS=2.0,\
+fsr_G2QQ_cNS_dn fsr:G2QQ:cNS=-2.0,\
+fsr_G2QQ_cNS_up fsr:G2QQ:cNS=2.0,\
+fsr_Q2QG_cNS_dn fsr:Q2QG:cNS=-2.0,\
+fsr_Q2QG_cNS_up fsr:Q2QG:cNS=2.0,\
+fsr_X2XG_cNS_dn fsr:X2XG:cNS=-2.0,\
+fsr_X2XG_cNS_up fsr:X2XG:cNS=2.0,\
+isr_G2GG_muR_dn isr:G2GG:muRfac=0.5,\
+isr_G2GG_muR_up isr:G2GG:muRfac=2.0,\
+isr_G2QQ_muR_dn isr:G2QQ:muRfac=0.5,\
+isr_G2QQ_muR_up isr:G2QQ:muRfac=2.0,\
+isr_Q2QG_muR_dn isr:Q2QG:muRfac=0.5,\
+isr_Q2QG_muR_up isr:Q2QG:muRfac=2.0,\
+isr_X2XG_muR_dn isr:X2XG:muRfac=0.5,\
+isr_X2XG_muR_up isr:X2XG:muRfac=2.0,\
+isr_G2GG_cNS_dn isr:G2GG:cNS=-2.0,\
+isr_G2GG_cNS_up isr:G2GG:cNS=2.0,\
+isr_G2QQ_cNS_dn isr:G2QQ:cNS=-2.0,\
+isr_G2QQ_cNS_up isr:G2QQ:cNS=2.0,\
+isr_Q2QG_cNS_dn isr:Q2QG:cNS=-2.0,\
+isr_Q2QG_cNS_up isr:Q2QG:cNS=2.0,\
+isr_X2XG_cNS_dn isr:X2XG:cNS=-2.0,\
+isr_X2XG_cNS_up isr:X2XG:cNS=2.0}',
+
+        'UncertaintyBands:nFlavQ = 4', # define X=bottom/top in X2XG variations
         'UncertaintyBands:MPIshowers = on',
         'UncertaintyBands:overSampleFSR = 10.0',
         'UncertaintyBands:overSampleISR = 10.0',

Configuration/Generator/python/VBFToH_Pow_JHU4l_LHE_13TeV_cff.py

@@ -5,7 +5,7 @@
 # https://github.com/cms-sw/genproductions/blob/060e6d3363a78ecfae90f7f5c46c968992820a56/bin/Powheg/production/V2/13TeV/Higgs/VBF_H_JHUGen_HZZ4L_NNPDF30_13TeV/JHUGen_VBF_H_ZZ4L.input
 
 externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
-    args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/slc6_amd64_gcc630/13TeV/Powheg/V2/RelValidation/VBFH/VBF_H_slc6_amd64_gcc630_CMSSW_9_3_9_patch1_VBF_new.tgz'),
+    args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/slc6_amd64_gcc630/13TeV/Powheg/V2/RelValidation/VBFH/VBF_H_slc6_amd64_gcc630_CMSSW_9_3_9_patch1_VBF_new_reducedPdf.tgz'),
     nEvents = cms.untracked.uint32(1000),
     numberOfParameters = cms.uint32(1),
     outputFile = cms.string('cmsgrid_final.lhe'),

Configuration/Generator/python/VHToH_Pow_LHE_13TeV_cff.py

@@ -4,7 +4,7 @@
 # https://github.com/cms-sw/genproductions/blob/329fda9f8d07c2d4d4e75c9a00279dcd6e78cda7/bin/Powheg/production/VH_from_Hbb/HWplusJ_HanythingJ_NNPDF30_13TeV_M125.input
 
 externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
-    args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/slc6_amd64_gcc630/13TeV/Powheg/V2/RelValidation/VH/HWJ_slc6_amd64_gcc630_CMSSW_9_3_9_patch1_VH_new.tgz'),
+    args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/slc6_amd64_gcc630/13TeV/Powheg/V2/RelValidation/VH/HWJ_slc6_amd64_gcc630_CMSSW_9_3_9_patch1_VH_new_reducedPdf.tgz'),
     nEvents = cms.untracked.uint32(5000),
     numberOfParameters = cms.uint32(1),
     outputFile = cms.string('cmsgrid_final.lhe'),