diff --git a/docs/part3/commonstatsmethods.md b/docs/part3/commonstatsmethods.md
index cb5ccdc13dc..4633d3933d2 100644
--- a/docs/part3/commonstatsmethods.md
+++ b/docs/part3/commonstatsmethods.md
@@ -920,8 +920,11 @@ This will produce a plot like the one below,
 
 ![](images/2D_LHScan.png)
 
-Similarly, 1D scans can be drawn directly from the tree, however for 1D likelihood scans, there is a python script from the [`CombineHarvester/CombineTools`](http://cms-analysis.github.io/HiggsAnalysis-CombinedLimit/#combine-tool) package [plot1DScan.py](https://github.com/cms-analysis/CombineHarvester/blob/113x/CombineTools/scripts/plot1DScan.py) that can be used to make plots and extract the crossings of the `2*deltaNLL` - e.g the 1σ/2σ boundaries.
+However for 1D likelihood scans, a python script to make plots and extract crossing of the `2*deltaNLL`, e.g the 1σ/2σ boundaries, is already available `plot1DScan.py`
 
+```
+plot1DScan.py higgsCombineTest.MultiDimFit.mH125.root -o output
+```
 
 ### Useful options for likelihood scans