 2018 for 1 h using the BM grism ($R=2000$ with the adopted large field of view), and on February 3rd and 4th 2019 for 3.5 h and 4 h, respectively, using the lower resolution BL grism ($R=900$), giving a total 8.5 h on-source when combining all observations. Conditions were excellent with dark sky and seeing typically in the range of 0.4–0.7$''$. We used integration times of 900 s in 2018 and of 1800 s in 2019, with dithering and large offsets to eventually cover a region corresponding to 2$\times$2 KCWI fields of view of $40''\times60''$ (i.e. $300\times500$ kpc$^2$). In fact, already with the first 1 h integration obtained in 2018 it was clear that the nebula extends well beyond the usable KCWI field of view of about $18''\times31''$ for the adopted configuration with large slices ($1.35''$). The low-resolution BL data covers the full 3500–5500Å range (corresponding to 900–1400 Å rest-frame at the redshift of RO-1001), while the BM observations cover a shorter range across the  emission. We used the standard KCWI pipeline [*Kderp*]{} (Morrissey et al. 2018; Neill et al. 2018) for the data reduction, including twilight flats to obtain accurate illumination corrections. We further used [*CubEx*]{} tools (Cantalupo et al. 2019) to refine the flat-fielding slice by slice, thus allowing us to improve the sky subtraction by removing a median sky value at each wavelength layer, after masking sources detected in the stacked (continuum) cubes. This step was first performed over the whole frame and subsequently iterated by masking regions where  emission had been detected, to avoid self-subtraction of the  signal. Further reduction and analysis steps were performed with the [*CWItools*]{} scripts (O’Sullivan et al. 2020; Martin et al. 2019). We estimate variance cubes from the original, non-resampled cubes and propagate the uncertainties through the combination to obtain a final variance cube. We combined the dithered and offset observations based on the astrometry of each frame that was derived by cross-correlating to B-band Suprime CAM imaging of the area publicly available
2018 for 1   h using the BM grism ($ R=2000 $ with the adopted large sphere of opinion), and on February 3rd and 4th 2019 for 3.5   h and 4   h, respectively, use the lower resolution BL grism ($ R=900 $), give a total 8.5   h on - beginning when combining all observations. condition were excellent with dark sky and seeing typically in the range of 0.4–0.7$''$. We used integration times of 900   second in 2018 and of 1800   s in 2019, with dithering and large outgrowth to eventually cover a region corresponding to 2$\times$2 KCWI airfield of view of $ 40''\times60''$ (i.e. $ 300\times500 $   kpc$^2 $). In fact, already with the first 1   h integration obtained in 2018 it was clean that the nebula extends well beyond the usable KCWI field of view of about $ 18''\times31''$ for the adopted configuration with large slices ($ 1.35''$). The low - resolution BL data covers the full 3500–5500Å   range (corresponding to 900–1400   Å   rest - human body at the redshift of RO-1001), while the BM observations brood a shorter range across the   emission. We used the standard KCWI grapevine [ * Kderp * ] { } (Morrissey et al. 2018; Neill et al. 2018) for the data decrease, including twilight flats to obtain accurate illuminance corrections. We further used [ * CubEx * ] { } tools (Cantalupo et al. 2019) to refine the flat - fielding slice by slice, thus allowing us to improve the sky subtraction by murder a median sky value at each wavelength level, after dissemble sources detected in the stacked (continuum) cubes. This step was first performed over the solid frame and subsequently iterated by masking regions where   discharge had been detected, to keep off self - subtraction of the   signal. Further reduction and analysis tone were performed with the [ * CWItools * ] { } script (O’Sullivan et al. 2020; Martin et al. 2019). We estimate variation cubes from the original, non - resampled cube and propagate the uncertainties through the combination to prevail a final discrepancy cube. We compound the dithered and offset observations based on the astrometry of each frame that was derived by crisscross - correlating to B - band Suprime CAM imagination of the area publicly available
 2018 flr 1 h using the BM grism ($V=2000$ with the adoptgd large field of view), and on February 3rd and 4th 2019 hor 3.5 y and 4 h, respectively, using ghe lower resolutuon UL grism ($R=900$), givinj a total 8.5 h on-alurcz xhen combining sll observdtions. Conditimnr cere excellent with dark sky and seeyng typocwlly in the ragge ps 0.4–0.7$''$. Ws used integration times of 900 s in 2018 znd of 1800 v in 2019, with diyhering and large offsets ho egentually cover a gegion corrgapogeing to 2$\timer$2 KCWI fields of view kf $40''\times60''$ (i.e. $300\times500$ kpc$^2$). In fact, alfeady with the diest 1 v integratiin obnained in 2018 it was cledr that the nebula exbends weol beyond the usable NCWI field of view of about $18''\tnmes31''$ for the adopted xobfigutatiot wigy lxrgt skides ($1.35''$). Tje mow-resolutjon BL data covers the full 3500–5500Å ramgq (correspondinf to 900–1400 Å rqse-frame at the redshift of RO-1001), while the TM kbservations cover a shirter range across thg  emission. We used the standard KCWI pipeline [*Kderp*]{} (Morrissay et xl. 2018; Ncill wt al. 2018) for the data reduction, including twilighe fkans to obtain accuvate illumination vogrrstions. We furjher uszs [*DubEx*]{} tools (Cantalkpo et wl. 2019) ti refine uhe fkat-fielding slice by slice, rhus allowiny uw to improve the sny subtractnon by remofing a median sky value at sach wavelejgth layed, after masking suurbes getected in the stacked (cogtinuum) cnbes. Chis steo wax firse performef ovev the whole frame ajd suysequantly iterwted by masking regions where  emission had begn getvcted, to cvoid xelf-subtractijn of the  signcl. Furthzr redjction and analysms steps werq performed whjh the [*CWItoons*]{} scripes (O’Wulluvan et xl. 2020; Martin et sl. 2019). We esniiate variance cubes frpm jhs original, non-rzrqmpled cubes anc pfopwgwtx the gncertaintiev thfouen the combinatiok tu obyain a final varianca cuge. We combined the dlthered abd offsee observationx based on the astgometcy of xach ftamg that was derived by cross-corrslating tl B-nand Suprime SAM lmagyng of the area publicly available
2018 for 1 h using the BM with adopted large of view), and 2019 3.5 h and h, respectively, using lower resolution BL grism ($R=900$), giving total 8.5 h on-source when combining all observations. Conditions were excellent with dark and seeing typically in the range of 0.4–0.7$''$. We used integration times of s 2018 of s in 2019, with dithering and large offsets to eventually cover a region corresponding to 2$\times$2 fields of view of $40''\times60''$ (i.e. $300\times500$ kpc$^2$). fact, already with the 1 h integration obtained in it clear that nebula well the usable KCWI of view of about $18''\times31''$ for the adopted configuration with large slices ($1.35''$). The low-resolution BL data the full (corresponding to Å at redshift of RO-1001), BM observations cover a shorter range We used the standard KCWI pipeline [*Kderp*]{} (Morrissey al. 2018; et al. 2018) for the data including twilight flats to obtain accurate illumination corrections. further used [*CubEx*]{} tools (Cantalupo et al. 2019) to refine the flat-fielding slice by slice, us to improve the subtraction by removing median value each layer, after sources detected in the stacked (continuum) cubes. This step was first over the whole frame and subsequently iterated by masking regions had detected, to avoid of the signal. Further and steps were performed with scripts et et 2019). estimate variance cubes from original, non-resampled cubes and propagate uncertainties through the combination cube. We combined the dithered and offset observations on the astrometry of each frame that derived by cross-correlating to B-band Suprime CAM imaging of the area publicly
 2018 for 1 h using the BM grism ($R=2000$ with tHe adopted lArge fIelD of ViEw), anD on FEbruary 3rd and 4th 2019 FOr 3.5 h aNd 4 h, respectively, using thE loweR rESoluTIoN BL grIsm ($R=900$), givINg A TOtaL 8.5 h On-SouRcE WhEn comBinIng all oBservationS. CoNdItions were exCElLent with daRk sKy and seeing tYpiCally iN tHe rANge of 0.4–0.7$''$. we uSed inTegratIOn timeS of 900 s in 2018 and Of 1800 S In 2019, with DItherinG ANd LargE offsets to eventuaLLy COver a region corResponDiNG tO 2$\TImeS$2 KCwI fields of ViEw of $40''\tIMes60''$ (i.e. $300\tiMEs500$ KPC$^2$). in fACt, already with The first 1 h inTEgrAtion oBtAinED in 2018 it wAs cleAr THat The nebula exTendS well beyoNd the uSAble KCWi Field of View of AboUt $18''\tImes31''$ FOr ThE adOpTEd cONfIguRAtiOn with laRgE sLices ($1.35''$). the lOW-RESoluTioN BL dAta coVers the full 3500–5500Å rAngE (corREspOndinG to 900–1400 Å rEst-fRaMe at tHe redsHift oF Ro-1001), while the BM obseRvatIons cover A shOrTer RaNge acROss the  EmiSsiOn. We useD the staNDarD Kcwi PiPeline [*Kderp*]{} (MorrissEy ET Al. 2018; neill et aL. 2018) for thE DaTa REduction, InCluDing TWIlighT flaTS tO obtain aCcuratE IlLuMinatioN cOrrectIoNs. WE fuRther USed [*CUbEx*]{} toOls (CantaLupo eT Al. 2019) to refine the fLAt-fielding sliCE bY SLiCE, thuS alLowing us to iMproVE the Sky sUBtRacTIon by RemovInG A mEDian sky value at each wAvElengtH layeR, after masking Sources detECTEd in the sTackED (cONtinuum) cubes. ThIs steP was first pERformed oVer thE whole frAme and subSEQuently iTerAteD by MasKINg Regions where  eMISsioN hAd been dEteCted, to aVoiD seLf-sUbtRaCtion of thE  signal. FUrThEr ReDucTion aND analysiS sTepS wEre PerfoRMed witH the [*CwItoOlS*]{} sCRipTs (O’SullIVaN ET al. 2020; MArTiN et aL. 2019). We EsTimatE varIAncE cubes fRom the oriGinAL, non-ReSaMpled cuBes and propagaTe The uncertaInTieS throuGH The combiNation to obtain a final varIAnce cubE. We CombiNed tHe dithereD anD offseT obSErvatiOns basEd on tHe AstROMetry OF EaCh fRaMe that was dERIveD by crOsS-corRelatinG to B-band Suprime CAM IMagIng of the area pUblIcly AVAiLabLE
 2018 for 1 h using the BM  grism ($R =2000 $ w ith  t he a dopt ed large field of v iew), and on February3rd a nd 4th2 01 9 for  3.5 ha nd 4  h,  r es pec ti v el y, us ing  the lo wer resolu tio nBL grism ($R = 90 0$), givin g a  total 8.5 h  on -sourc ewhe n  comb ini ng al l obse r vation s. Condit io n s were excelle n t  w ithdark sky and seei n gt ypically in th e rang eo f0 . 4–0 .7$ ''$. We us ed  inte g rationt im e s of9 00 s in 2018and of 1800 s i n 2019 ,wit h  dithe ringan d  la rge offsets  toeventuall y cove r  a regi o n corre spondi ngto2$\t i me s$ 2 K CW I  fi e ld s o f  vi ew of $4 0' '\ times 60'' $ ( i .e.$30 0\ti mes50 0$ kpc$^2$).Infact ,  al ready  with  the  f irst1 h in tegra ti on obtained in2018  it was c lea rtha tthe n e bula e xte nds  well b eyond t h e u sa b l e  K CWI field of viewof a bo ut $18'' \times 3 1' '$ for thead opt ed c o n figur atio n  w ith larg e slic e s($ 1.35''$ ).  The l ow -re sol ution BL d ata co vers the  full 3500–5500Å ran g e (correspond i ng t o9 00–1 400  Å rest-fra me a t  the  red s hi fto f RO- 1001) ,w hi l e the BM observatio ns  cover  a sh orter range a cross thee m ission.We u s ed the standard K CWI p ipeline [* K derp*]{}  (Mor rissey e t al. 201 8 ;  Neill e t a l.201 8)f o rthe data redu c t ion,  i ncludin g t wilight  fl ats  to  ob ta in accura te illum in at io ncor recti o ns. We f ur the ruse d [*C u bEx*]{ } too ls ( Ca nt a lup o et al .  2 0 1 9) t ore fine  th eflat- fiel d ing  sliceby slice,  th u s al lo wi ng us t o improve the  s ky subtrac ti onby rem o v ing a me dian sky value at eachw aveleng thlayer , af ter maski ngsource s d e tected  in th e sta ck ed( c ontin u u m)  cu be s. This st e p  wa s fir st  per formedover the whole fra m e a nd subsequent lyiter a t ed  by ma s kin gr egi o n s where  emissi on had bee nd et ected, toa voi dself-su btracti on of the  si gnal. Fur ther redu ct iona n d a nalysis st eps were  performe d  with th e [*C WIt ools*] {}  sc ripts  (O’Su l liv an et  al. 2 02 0; Mar tin e tal. 2019 ). We estimate variance  cubes  from  th e origina l,n on- resampled  cub es and pro pag ate  theunc e rtain ties th rou g h the  com b ination t o  o bta i n  a  final vari a n c e c ube.Wec ombine d th e dithered and of f set observatio ns b a s edont he a st rometry of eac h f ra m e  that wa sderived bycross-co rr e latin g to B -bandSuprime C AM imagin g of  th e area pu bli cl y  availa bl e
 2018_for 1 h_using the BM grism_($R=2000$ with_the_adopted large_field_of view), and_on February 3rd_and 4th 2019 for_3.5 h and 4 h,_respectively,_using the lower resolution BL grism ($R=900$), giving a total 8.5 h on-source when combining_all_observations. Conditions_were_excellent_with dark sky and seeing_typically in the range of_0.4–0.7$''$. We_used integration times of 900 s in 2018 and_of_1800 s in 2019,_with dithering and large offsets to eventually cover a_region corresponding to 2$\times$2 KCWI fields_of view of_$40''\times60''$_(i.e._$300\times500$ kpc$^2$). In fact, already_with the first 1 h integration obtained_in 2018 it was clear that_the nebula extends well beyond the usable_KCWI field of view of about_$18''\times31''$ for the adopted configuration_with large_slices ($1.35''$). The low-resolution BL_data covers the_full 3500–5500Å range_(corresponding to 900–1400 Å rest-frame_at the redshift of RO-1001), while_the BM observations_cover a shorter range across the_ emission._We used the_standard_KCWI_pipeline [*Kderp*]{}_(Morrissey et al._2018;_Neill et_al._2018) for the data reduction, including_twilight_flats to obtain accurate illumination corrections. We_further used [*CubEx*]{} tools_(Cantalupo_et al. 2019) to_refine the flat-fielding slice by_slice, thus allowing us to improve_the sky_subtraction by_removing a median sky value at each wavelength layer, after masking_sources detected in the stacked (continuum)_cubes. This step was_first performed_over_the whole frame_and_subsequently iterated_by masking regions where  emission had been_detected, to_avoid self-subtraction of the  signal. Further_reduction and analysis steps_were_performed with the [*CWItools*]{} scripts (O’Sullivan_et al. 2020; Martin et al._2019). We estimate variance cubes_from_the_original, non-resampled cubes and propagate_the uncertainties through the combination to_obtain a final_variance cube. We combined the dithered and_offset_observations based on the astrometry of_each_frame that was derived by cross-correlating_to_B-band_Suprime CAM imaging of the_area publicly available