• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A ROOT-based detector test system

    2021-11-13 01:30:42HuangKaiWuChenLi
    Nuclear Science and Techniques 2021年10期

    Huang-Kai Wu · Chen Li

    Abstract Most detectors for nuclear physics experiments are detector arrays composed of numerous units. Testing each detector unit is a major part of the research work on detector arrays. To save time and simplify the research process, a ROOT-based detector test system was designed for detector unit testing. The test system is a general purpose and expandable software system that can support most of the hardware devices in the market. Users can easily build a complete detector test system using the required hardware devices. The software is based on the ROOT framework and is operated on the Linux platform. The software of the test system consists of four parts: the controller, data acquisition(DAQ), high-voltage power supply, and online monitoring and analysis. In addition, a user-friendly graphical user interface (GUI) was designed for user convenience. Moreover, the online analysis function of the software can implement automatic peak searching and spectrum fitting for different radioactive sources, and the results under different conditions can be shown automatically. The completion of the test system could greatly simplify the development process of the detector.

    Keywords Detector test · Data analysis · ROOT

    1 Introduction

    Detectors in domestic and foreign laboratories are usually detector arrays,which are composed of many identical detector units. For example, the Advanced Gamma Tracking Array (AGATA) is the next generation γ-ray spectrometer of the European project, which is composed of 180 HPGe detector units, covering 4π solid angles[1].The ELI-NP GANT-Thermal Neutrons array (ELIGANT-TN) is a 4π neutron detector consisting of 303He counters embedded in polyethylene, which is used for photo-neutron cross section measurements [2]. The ELINP Silicon Strip Array(ELISSA)consists of three rings of twelve single-sided X3 Si strip detectors,that are designed for photo-dissociation reactions relevant to Big Bang nucleosynthesis (BBN), supernova explosions, and p-process studies[3].Other types of detector units include HPGe detector arrays [4-7], neutron detector arrays [8-11], silicon detector arrays[12-15],and the liquid scintillator dark matter detector [16-18]. Consequently, unit tests are an important process in the development of detector arrays.To obtain the best operating parameters, time-consuming and repeated unit tests are often required. During the development of the CsI detector array designed by the Shanghai Institute of Applied Physics encountered the same problem[19]. To overcome this problem, a ROOT-based detector test system was developed that can automatically adjust the working parameters of the detector and analyze the data.Moreover,owing to the modular design of this test system,it can be easily expanded into a versatile test system to meet different test requirements.

    ROOT, a common data analysis software package used in nuclear physics, has both a powerful data analysis function and a flexible graphical user interface(GUI)design function. ROOT has been widely used in many projects. Examples of such projects include the DAQ system for the silicon pixel detector, NDAQ system of the Heavy Ion Medical Machine(HIMM)of Lanzhou Institute of Modern Physics, and the General-Purpose Digital Data Acquisition System (GDDAQ)of Peking University[20-23]. Therefore, ROOT [24] is suitable for implementation of the detector test system. Along with the development of advanced hardware, many newer nuclear electronic modules support remote control.The maturity of hardware and software makes it possible to develop a test system that can simplify the test process of the detector unit.

    Compared with the common DAQ system, our detector test system has the following characteristics. First, the test system can support more types of hardware devices than the DAQ system. The test system not only includes the DAQ function, but also the functions of controlling the main amplifier and high-voltage power-supply module.Second, for detector testing, the test system runs quickly,sometimes requiring only 10 min.After sufficient statistics are acquired, the DAQ loop can be completed. Third, to obtain the best parameters, the voltage and gain of the detector should be treated as variables under each test.Since the radioactive source is known, the analysis algorithm is universal. Therefore, the test system can integrate the data analysis function into the software. In contrast,common DAQ systems do not integrate data analysis functions in the software since different experiments have different experimental purposes.

    The remainder of this paper is organized as follows:Sect. 2 presents the hardware architecture of the test system.Section 3 describes the software framework and GUI.Finally,the feasibility and stability of the system are tested in Sect. 4, and a summary is provided in Sect. 5.

    2 Hardware support of test system

    Detector test systems usually involve several functional modules, such as a high-voltage power supply module,preamplifier, main amplifier, and analog-to-digital conversion module(ADC, TDC, flash ADC, etc.). Depending on the type of interface,modules can be divided into three types: NIM, VME, and desktop modules. Therefore, the software of the test system needs to support these three types of hardware devices. Usually, the connection between a device and a PC is through a USB interface or an optical fiber link. As shown in Fig. 1, the hardware of the test system can be divided into four parts: the desktop module, NIM, VME, and a PC for control and data acquisition. The instruction communication between various hardware devices is implemented uniformly through the controller module of the software.Data can be obtained not only by modules with VME interfaces, but also by desktop modules, such as DT5720. The trigger signal can be provided externally via the front panel trigger input, as well as via the software. Alternatively, each channel can generate a self-trigger signal when the digitized input pulse exceeds a configurable threshold set through the register address. DT5720 houses USB 2.0, and the optical link interfaces. USB 2.0, allows data transfers up to 30 MB/s.The optical link supports a transfer rate of 80 MB/s and offers a daisy chain capability.This was successfully added to the test system.Normally,NIM devices can also be used for data acquisition. However, because of fewer functions and the existence of alternative products,NIM devices with DAQ functions were not considered for the test system design.

    A CsI detector test system was built to verify the feasibility of the test system. The principal components of hardware for this system are the VME and NIM modules,as well as the desktop high voltage module with independent USB interfaces. The VME modules used were ADC(V785) and MADC-32 [25, 26]. In addition, the data acquisition module also includes DT5720 [27], which is linked to the PC via the USB interface. The NIM module includes the shaping amplifier MSCF-16 [27] and the gate and delay generator GG8020 manufactured by ORTEC.The VME controller selects the V2718 manufactured by CAEN and is connected to the PCI card using optical fibers. The maximum data transfer rate is 80MB byte/s.DT5533 is a desktop module housing four independent high voltage power supply channels with an independent USB interface, which can be controlled and monitored by software [28]. The specific test process is presented in Sec. 4.

    3 Description of test system

    The following software was developed to conveniently obtain and handle the raw data: the Controller module,HV_module, DAQ_module, and Online_module. Instruction communication and data transmission between different components are implemented through shared memory. The Controller_module is the main control part,which is used to send data acquisition instructions and feedback status parameters of hardware devices. The data acquisition software, DAQ_module, is responsible for implementing data acquisition, transmission, and storage functions according to the received instructions. The HV_module is responsible for adding voltage to the detector according to the working parameters set by the Controller_module. The Online_module’s monitoring and analysis components can implement the online data monitoring and analysis functions.

    3.1 Framework of software

    The relationships between the four parts are shown in Fig. 2. The Controller_module consists of two parts:DAQ_Control and HV_Control. DAQ_Control is responsible for sending data acquisition instructions and setting data acquisition parameters. The corresponding instructions are transmitted to the shared memory,received, and executed by the DAQ. HV_Control is responsible for setting the working parameters of highvoltage devices. Similarly, the parameters are written into the shared memory and read by HV_Control. After HV_Control receives the working parameters from the shared memory, the high-voltage module DT5533 is initialized with the setting value, and the corresponding voltage is applied to the detector. Each time HV_Control completes a voltage addition operation; it feeds back the status parameters to the Controller_module. The control function of the Controller_module is realized by calling the shared memory class and the control GUI class defined in the software. The parameter passed to the shared memory class is the defined structure, which includes the working parameters of the high-voltage devices,the data acquisition instructions, and the corresponding feedback parameters.The implementation of HV_Control is based on the shared memory class and HV_class defined in the software. The HV_class calls the CAENHV Wrapper library, which is a software library for power supply control,issued by CAEN to realize the control of DT5533.

    Fig. 2 Frame diagram of software

    The data acquisition module (DAQ_module) receives commands from the Controller_module via shared memory.When the DAQ_module receives the start command,it controls the VME bus to acquire the data. After each data acquisition loop ends, an increasing voltage signal is fed back to the HV_module which waits for the next instruction from the Controller_module. Once the trigger logic triggers a physical event for data acquisition,the events are recorded in the register of the front-end electronics(V785,V775, MADC-32,etc.), and then read by the external firstin first-out (FIFO) memory. Finally, the data in FIFO are transmitted to the PC via optical fiber and stored in the hard disk in binary format for future offline analysis and query.Simultaneously, the raw data will also be transmitted to shared memory for online monitoring through the write buffer. Alternatively, the data can also be acquired by a digitizer[29].The input pulses,which come from detectors coupled with a charge-sensitive preamplifier (such as HPGe,silicon,etc.)or detectors coupled with PMT,SiPM,etc.,. are digitized into many samples by the flash ADC.Based on the sampling frequency of the digitizer,users can set the number of sampling points according to the width of the waveform by the configuration file. For example, the sampling frequency of DT5720 is 250 MS/s which means sampling a point every 4 ns if a waveform with a width of 1 μs, 250 points were sampled. Each sampled data set can be analyzed to extract the information of time,energy,and so on using the software. Because of the differences in the interface, the base class of the data acquisition modules in the software is divided into two types:the base class of the module with the VME interface and USB interface. The modules with VME_interface, such as ADC (V785,MADC-32, etc.), TDC (V775, MTDC-32, etc.), and QDC(V792,MQDC-32,etc.)are inherited from the base class of the VME interface. The flash ADC contains both modules with a VME interface and a USB interface, which are inherited from two different base classes. Although the flash ADC is divided into two categories, the basic functions are the same. The obtained data are transmitted through the shared memory class, and the parameters passed to the shmget function in the shared memory class are transferred as an array of float type or int type.

    The Online_module is mainly composed of two parts:monitoring and analysis. Periodically (every 10 s in the present test system), raw data stored in the front-end electronic registers (e.g., the event number, time stamp,and each channel’s energy) can be monitored by online processes through the shared memory mechanism (Online,Fig. 2).The received raw data has a defined data structure,such as the existence of a header(including channel, ID,event count, etc.), datum, and end of event mark. Therefore, to achieve an online monitoring function, it is necessary to call the RawData class(defined in the software)to extract the necessary information from the raw data. In addition to the real-time display of the spectrum of the obtained data, the Online_module can automatically select the analysis algorithm to search for the peak position and set the fitting parameters to fit different spectra according to the radioactive source used in the test. The analysis function of the Online_module first converts the binary data into a ROOT format file after the data acquisition ends, and then analyzes the ROOT file using the corresponding algorithm. During detector unit tests, all data acquisition parameters and analysis processes are saved as a log file for future queries. Furthermore, the analysis results of energy resolution, peak position, and other measurements will also be stored in the disk in txt format.

    3.2 Graphical user interface (GUI)

    To facilitate user operation, a graphical user interface(GUI) was designed. The GUI is divided into two parts:the main control interface and the online monitor and analysis interface.The GUI was compiled based on ROOT under the Linux system. As shown in Fig. 3, the Controller_module serves as the hub for instruction transfer and parameter adjustment. It is composed of two parts:DAQ_Control and HV_Control, which are associated with the DAQ_module and HV_module, respectively.

    The DAQ_Control part of the GUI shown on the left side of Fig. 3a, contains the function of setting the data acquisition parameters. Users can run the start/stop command via the button in the top ’’control’’ region. In the middle ’’setup’’ region, users can set the file name header and running time of every data acquisition. The high voltage value added to the detector is automatically used as a supplement to complete the file name. All raw data were stored in binary format. The information region shows a simple log of the test system. The HV control part shown on the right side of Fig. 3b, allows users to control the high-voltage channel on/off in the left option box, as well as select the current monitoring range and accuracy. The right part of the HV control shows two high-voltage channel control interfaces. The control interface of each channel displays the default working parameters (modifiable) of DT5533, the monitoring voltage, and current values of the detector. The bottom of each control interface was used to set the voltage range and step interval of the detector test.

    Fig. 3 Interface of Control’s GUI:The figure shows the main control interface, through which all the acquisition commands and working parameters can be completed

    The main interface of the Online_module contains an online monitor and analysis function. The monitor part reads the raw data from the shared memory in real time.The obtained raw data are decoded first, and then transmitted to the online monitor for spectrum display. The monitored spectrum was updated every ten seconds.Moreover, the monitored channels of each module can be selected in the online GUI. There are also some buttons with specific functions on the online GUI interface. For example,the Online_module can realize an automatic data analysis function through the button ’’Analysis’’ on the online GUI. The details of the automatic data analysis function are presented in the next section.

    4 System test by CsI detector

    The test system was designed for the detector unit test.The feasibility and stability of the system were tested using the CsI detector unit.The CsI detector is one unit of the CsI detector array developed by the Shanghai Institute of Applied Physics,Chinese Academy of Sciences.This array is composed of 42 CsI detector units, which are divided into external and internal layers. The outer layer is composed of 18 units, and the inner layer is divided into two parts;each part consists of 12 units.This array is originally designed to measure the neutron capture (n,γ) crosssection.

    4.1 Structure of CsI detector and test

    The CsI(TI) crystal of each detector unit of the detector array is a hexagonal prism with a length of 200 mm, and the distance between the opposite sides of the hexagon is 55 mm.The crystal was produced by the Shanghai Institute of Ceramics, Chinese Academy of Sciences, and covered with 1 mm polytetrafluoroethylene(PTFE) and 2 mm thick aluminum as shielding materials. The photomultiplier tube connected to one side of the crystal was R1828-01 (Hamamatsu). The detailed construction of the detector unit can be found in Ref. [19].

    To verify the feasibility of the test system, two types of radioactive sources,137Cs and60Co, were used to test the performance of the CsI(TI) detector unit. During the test process,the measurement voltage range was set from 1600 to 2100 V,and the measurement interval was set to 100 V.The acquisition time for each condition was set to 10 min.When the Controller_module sends the start instruction to the DAQ_module, it enters the acquisition loop and waits for the feedback signal from the high voltage module(HV_Control). The HV_Control module receives the working parameters from the shared memory. After receiving the corresponding parameters,a high voltage was added to the detector at an initial value of 1600 V. During the voltage raising process, HV_Control feeds back the control status parameters of DT5533 control in real time.Until the voltage-raising operation is completed, DAQ_-module starts to control the VME bus to obtain a set of data under 1600 V.When the running time reaches the set time of 10 min, DAQ_module ends the acquisition loop and returns a raising-voltage signal to HV_Control again.After receiving this signal, HV_Control will raise the voltage to the next target value at a measurement interval of 100 V and feedback the status parameters to control in real time.Then, the Controller_module sends the start command again, and the DAQ_module controls the VME bus to obtain another set of data under 1700 V. The above steps are repeated until the test system obtains the data under 2100 Vat.

    Figure 4 shows the energy spectrum of137Cs and60Co obtained by CsI(TI) detector under many different voltages,which are acquired by ADC module V785 or MADC-32. As mentioned in Sect. 3.1, the digitizer, such as DT5720, has been added into the test system. Therefore,the data can also be acquired by digitizer.In this study,the trigger signal of the digitizer is generated by a self-trigger mode. The raw waveform is recorded as 300 samples according to the actual width of waveform output by the CsI detector. Although a single raw data file obtained by DT5720 will take up a large amount of storage space, the advantage is that it retains all the information of the output signal of the detector. Several algorithms have been integrated into the test system to extract energy and time information from sampling data. For example, users can extract the time information via the method of constant fraction discrimination (CFD) or leading-edge discrimination (LED). Since the spectrum obtained by DT5720 is consistent with that obtained by V785 and MADC-32,it is no longer presented in this paper.

    4.2 Data online analysis

    During the test process of the detector arrays, a great amount of data is acquired resulting in large files.This can result in a time-consuming and tedious offline data analysis. Because the spectrum of radioactive sources used for detector tests is usually known [30], it is possible to program automatic peak searching and spectrum fitting according to different radioactive sources. Users then only need to declare the appropriate radioactive source in the file name, and the Online_module can select the corresponding algorithm to process different raw data.

    Typically,the raw data obtained by the analog-to-digital conversion module(ADC,TDC,etc.)are a series of integer values, called channel numbers, which are determined by the resolution of the module. This needs to be converted into the corresponding energy by an energy calibration.Generally, the calibration function of the scintillator detector is nonlinear, and the calibration functions of the silicon detectors and HPGe detectors are linear. Therefore,at least three calibration points are required for the energy calibration.In this study,137Cs and60Co sources were used to test the performance of the CsI detector. There are two peaks in the γ-spectrum of60Co, and the corresponding energies are 1.173 and 1.332 MeV, respectively. There is one peak in the γ-spectrum of137Cs,and the corresponding energy was 0.662 MeV.It meets the minimum requirement for energy calibration.

    For online monitoring and analysis, the raw data must first be decoded,and then the binary file converted into the ROOT file.After this,the Online_module fits and searches for the peak of the spectrum to obtain the channel numbers corresponding to the peaks of the two types of radioactive sources under various voltages, as shown in Fig. 4. Generally, the fitting methods for the137Cs and60Co spectraare different.Because the energy spectrum of137Cs has an exponential background and a single Gaussian peak, the fitting algorithm is carried out using a single Gaussian plus an exponential function. However, for the spectrum of60Co,the fitting algorithm adopts double Gaussian function plus an exponential function because of the existence of double Gaussian peaks and an exponential shape. The corresponding fitting algorithm is integrated into the software.Users only need to declare a radioactive source in the file name.

    Fig. 4 Energy spectrum of CsI(TI) detector under 137Cs and 60Co radioactive sources. a 137Cs 1600 V-2100 V; b 60Co 1600 V-2100 V

    After fitting the spectrum, the calibration curve of the two sources under different voltages can be obtained, as shown in Fig.5.The horizontal axis represents the channel number, and the vertical axis represents the energy. When the γ-ray energy was lower than 1.332 MeV, there was a good linear relationship between the energy and channel numbers. As presented above, the calibration curve of scintillator detectors is usually nonlinear. However,because the tested CsI(TI) detector is a large-size detector,high-energy γ photons can deposit all their energy in a CsI(TI) crystal. This makes the calibration curve of the CsI(TI)detectors linear.A detailed discussion can be found in found in [19]. Figure 6 shows one of the results of the analysis. Based on this figure, it can be inferred that a maximum voltage limit can be added to the detectors to prevent the DAQ from exceeding the effective range to excessive working voltage.

    Fig. 6 Curve of peak position versus voltage

    To compare the performance of the detector under different voltages, it is necessary to calculate the absolute energy resolution under different conditions. As shown in Figs. 4 and 5, the spectrum fitting and energy calibration under different voltages were completed. Therefore, the FWHM of each peak can be calculated according to the Gaussian fitting parameter using the formula FWHM =2.355σ,and then it can be converted into energy according to the calibration curve. Figure 7 shows the absolute energy resolution of the CsI(TI) detector under different voltages.It can be inferred from the figure that,when the γray energy is lower than 1.332 MeV,the energy resolution improves with an increase in voltage. In addition to the graphs shown above, other analysis results, such as the energy calibration function, Gaussian fitting parameter σ,absolute energy resolution, and relative energy resolution,will be recorded in the txt file, which is convenient for users to archive and query.

    Fig. 5 Energy calibration on different voltage: 1600 V-2100 V

    Fig. 7 Absolute energy resolution with a gain of 3.3

    5 Summary

    This paper reports a ROOT-based detector test system that can automatically adjust the working parameters of the detector and online data analysis. The software of the test system is compiled based on ROOT and runs under the Linux operating system.A user friendly GUI helps users to monitor and control the detector test system easily and efficiently.The test system was successfully applied to the test of a CsI(TI) detector array.

    Owing to the modular design of the test system, the system can be easily expanded to match different types of detectors. At present, this study only implements the automatic adjustment of the voltage parameters. However,the test modes supported by the test system can be expanded by adding different function modules. For example,control of the shaping amplifier MSCF-16 can be added to implement the automatic adjustment of gain through the software. A stepping motor slide rail can also be added to the test system to implement a change in the incident position of the radioactive source by controlling the motor slide rail.

    In summary, the feasibility and stability of the test system were verified through a CsI(TI) unit test. This shows that the design idea of this test system is feasible and can greatly simplify the test process of the detector unit.Because of the modular architecture of the software, the test system has versatility, scalability, and many useful application prospects.

    Author ContributionsAll authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Huang-Kai Wu and Chen Li.The first draft of the manuscript was written by Huang-Kai Wu and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.

    在线观看www视频免费| 91大片在线观看| 1024香蕉在线观看| 757午夜福利合集在线观看| 亚洲欧美一区二区三区黑人| 亚洲精品在线美女| 亚洲国产中文字幕在线视频| 一个人免费在线观看的高清视频| av免费在线观看网站| 欧美丝袜亚洲另类 | 日本五十路高清| 成熟少妇高潮喷水视频| 亚洲片人在线观看| 亚洲国产日韩欧美精品在线观看 | 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 嫩草影院精品99| 免费看a级黄色片| 成人高潮视频无遮挡免费网站| 国产69精品久久久久777片 | 一本大道久久a久久精品| 老司机午夜十八禁免费视频| 国产精品九九99| 老熟妇仑乱视频hdxx| www.精华液| 极品教师在线免费播放| 天天躁夜夜躁狠狠躁躁| 久久精品国产99精品国产亚洲性色| 色综合站精品国产| 久久精品亚洲精品国产色婷小说| 亚洲成人久久爱视频| 大型黄色视频在线免费观看| 99久久综合精品五月天人人| 亚洲精品粉嫩美女一区| 日韩国内少妇激情av| 国产v大片淫在线免费观看| 国产高清视频在线观看网站| 一边摸一边抽搐一进一小说| 正在播放国产对白刺激| 69av精品久久久久久| 男人舔女人下体高潮全视频| 午夜老司机福利片| 亚洲精品色激情综合| 亚洲人成网站在线播放欧美日韩| 国产单亲对白刺激| 国产欧美日韩精品亚洲av| 精品少妇一区二区三区视频日本电影| 美女免费视频网站| 日韩av在线大香蕉| 日韩免费av在线播放| 人妻丰满熟妇av一区二区三区| 真人一进一出gif抽搐免费| 亚洲国产欧洲综合997久久,| 妹子高潮喷水视频| 国产探花在线观看一区二区| 欧美黑人巨大hd| 哪里可以看免费的av片| 不卡av一区二区三区| 啦啦啦免费观看视频1| 身体一侧抽搐| 波多野结衣巨乳人妻| 日本在线视频免费播放| 久久中文看片网| 看黄色毛片网站| 精品一区二区三区四区五区乱码| 男女之事视频高清在线观看| 久久中文看片网| 久久精品国产亚洲av香蕉五月| 欧美日韩国产亚洲二区| 久久久久国产一级毛片高清牌| 午夜亚洲福利在线播放| 亚洲精品在线观看二区| 日日爽夜夜爽网站| 国产高清videossex| 国产成人精品久久二区二区91| 久久亚洲精品不卡| 在线播放国产精品三级| 国产精品日韩av在线免费观看| 精品欧美一区二区三区在线| 国语自产精品视频在线第100页| 国产精品一区二区免费欧美| 99久久无色码亚洲精品果冻| or卡值多少钱| 日本精品一区二区三区蜜桃| 亚洲专区中文字幕在线| 中文字幕久久专区| 国产成人av教育| 99热只有精品国产| 亚洲熟女毛片儿| 久9热在线精品视频| 村上凉子中文字幕在线| 97碰自拍视频| 久久久久亚洲av毛片大全| 好男人电影高清在线观看| 欧美不卡视频在线免费观看 | 亚洲一区二区三区色噜噜| 国产一级毛片七仙女欲春2| 又紧又爽又黄一区二区| 欧美三级亚洲精品| ponron亚洲| 国内精品一区二区在线观看| 国产精品 欧美亚洲| 看免费av毛片| 亚洲精品在线美女| 国产精品亚洲一级av第二区| 美女大奶头视频| www.999成人在线观看| 亚洲,欧美精品.| 免费看美女性在线毛片视频| 免费在线观看亚洲国产| 日本 欧美在线| 久久久久性生活片| 夜夜躁狠狠躁天天躁| 99国产精品一区二区三区| 日韩免费av在线播放| 久久中文看片网| 亚洲精品久久国产高清桃花| 欧美中文综合在线视频| 欧美一级a爱片免费观看看 | 国产三级在线视频| 九色成人免费人妻av| 青草久久国产| 欧美最黄视频在线播放免费| 日韩有码中文字幕| 国产免费av片在线观看野外av| 亚洲专区中文字幕在线| 国产精品99久久99久久久不卡| 精品一区二区三区四区五区乱码| 久久性视频一级片| 久久久久国产精品人妻aⅴ院| 老司机深夜福利视频在线观看| 欧美成人一区二区免费高清观看 | 在线观看免费视频日本深夜| 久久久久久免费高清国产稀缺| 色av中文字幕| 国产片内射在线| 亚洲精品色激情综合| 国产精品一区二区三区四区免费观看 | а√天堂www在线а√下载| 麻豆成人午夜福利视频| 精品久久久久久久久久久久久| 欧美国产日韩亚洲一区| 又黄又粗又硬又大视频| 黄色 视频免费看| 久久久国产欧美日韩av| 久久人人精品亚洲av| 色播亚洲综合网| 亚洲一区高清亚洲精品| 成年人黄色毛片网站| aaaaa片日本免费| 国内精品久久久久精免费| 久久精品亚洲精品国产色婷小说| 在线观看免费视频日本深夜| 欧美午夜高清在线| 宅男免费午夜| 欧美色视频一区免费| 国产成人啪精品午夜网站| 特级一级黄色大片| 欧美+亚洲+日韩+国产| 久久天堂一区二区三区四区| 亚洲精品久久国产高清桃花| 亚洲欧洲精品一区二区精品久久久| 极品教师在线免费播放| 亚洲人成伊人成综合网2020| 啦啦啦韩国在线观看视频| 日本黄大片高清| 欧美+亚洲+日韩+国产| av中文乱码字幕在线| 天堂√8在线中文| 麻豆av在线久日| 亚洲一卡2卡3卡4卡5卡精品中文| 色av中文字幕| 在线观看66精品国产| 男人舔奶头视频| 亚洲av中文字字幕乱码综合| 久久国产精品影院| 巨乳人妻的诱惑在线观看| 一夜夜www| 日本一二三区视频观看| 欧美最黄视频在线播放免费| 精品电影一区二区在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产欧美网| 精品午夜福利视频在线观看一区| 50天的宝宝边吃奶边哭怎么回事| 亚洲av成人一区二区三| www日本黄色视频网| 国产精品98久久久久久宅男小说| 国产精品久久久av美女十八| 亚洲av第一区精品v没综合| 亚洲精品色激情综合| 给我免费播放毛片高清在线观看| 一本久久中文字幕| 校园春色视频在线观看| 91在线观看av| 在线视频色国产色| 在线观看美女被高潮喷水网站 | 两性夫妻黄色片| 51午夜福利影视在线观看| 日韩精品中文字幕看吧| 又紧又爽又黄一区二区| 在线观看日韩欧美| 国产视频内射| 午夜精品久久久久久毛片777| 欧美日韩一级在线毛片| 午夜老司机福利片| 国产精品日韩av在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | 我的老师免费观看完整版| 欧美日韩瑟瑟在线播放| 亚洲中文字幕日韩| 村上凉子中文字幕在线| 9191精品国产免费久久| 老司机在亚洲福利影院| 国产成人影院久久av| 亚洲成人精品中文字幕电影| 久久精品aⅴ一区二区三区四区| 久久久水蜜桃国产精品网| 听说在线观看完整版免费高清| 一本大道久久a久久精品| 欧美激情久久久久久爽电影| 国产精品亚洲一级av第二区| 男女那种视频在线观看| 精品无人区乱码1区二区| 狂野欧美激情性xxxx| 美女 人体艺术 gogo| 亚洲国产欧美人成| 亚洲av五月六月丁香网| 天堂影院成人在线观看| 日本撒尿小便嘘嘘汇集6| 国产人伦9x9x在线观看| 美女 人体艺术 gogo| 法律面前人人平等表现在哪些方面| 精品无人区乱码1区二区| 18美女黄网站色大片免费观看| 欧美成人一区二区免费高清观看 | 婷婷精品国产亚洲av| av欧美777| 19禁男女啪啪无遮挡网站| cao死你这个sao货| 国产亚洲av嫩草精品影院| 久久久精品大字幕| 亚洲一区二区三区不卡视频| 成人三级黄色视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲熟女毛片儿| 老汉色av国产亚洲站长工具| 国产真实乱freesex| 99热6这里只有精品| 日本一本二区三区精品| 好男人电影高清在线观看| 日本免费a在线| 国产激情久久老熟女| 精品少妇一区二区三区视频日本电影| 女人高潮潮喷娇喘18禁视频| 最近最新中文字幕大全免费视频| 国产伦人伦偷精品视频| www.999成人在线观看| 女同久久另类99精品国产91| 脱女人内裤的视频| 成人一区二区视频在线观看| 精品少妇一区二区三区视频日本电影| 巨乳人妻的诱惑在线观看| 精品久久久久久久毛片微露脸| 免费看美女性在线毛片视频| 欧美黑人巨大hd| 精品一区二区三区四区五区乱码| 精品一区二区三区四区五区乱码| 国产成人欧美在线观看| www国产在线视频色| 午夜老司机福利片| 国产在线精品亚洲第一网站| 国产精品久久视频播放| 99久久精品国产亚洲精品| 久久久精品国产亚洲av高清涩受| 男女床上黄色一级片免费看| 午夜福利18| 国产精品久久久人人做人人爽| 国产蜜桃级精品一区二区三区| 亚洲,欧美精品.| 亚洲中文av在线| 亚洲成人国产一区在线观看| 免费av毛片视频| 国产精品99久久99久久久不卡| 日本五十路高清| av有码第一页| av超薄肉色丝袜交足视频| 深夜精品福利| 波多野结衣高清作品| 岛国在线免费视频观看| 欧美色视频一区免费| 久久久国产成人免费| 欧美国产日韩亚洲一区| 国产精品乱码一区二三区的特点| 99国产精品99久久久久| tocl精华| 最近在线观看免费完整版| 欧美成狂野欧美在线观看| 国产成人aa在线观看| 制服诱惑二区| 日韩大尺度精品在线看网址| 久久热在线av| 成人三级做爰电影| 夜夜夜夜夜久久久久| 欧美在线一区亚洲| 男插女下体视频免费在线播放| 婷婷精品国产亚洲av在线| 国产精品1区2区在线观看.| 日韩中文字幕欧美一区二区| 国产aⅴ精品一区二区三区波| 黑人巨大精品欧美一区二区mp4| 午夜免费成人在线视频| 三级男女做爰猛烈吃奶摸视频| 色综合婷婷激情| 国产精品98久久久久久宅男小说| 国产一区二区在线观看日韩 | 亚洲精品在线观看二区| 亚洲av成人av| 成人国语在线视频| 日本撒尿小便嘘嘘汇集6| 国产精品99久久99久久久不卡| 一区二区三区国产精品乱码| 久久国产精品人妻蜜桃| av国产免费在线观看| 国产伦在线观看视频一区| 成人一区二区视频在线观看| 亚洲 欧美 日韩 在线 免费| 国产一区二区在线观看日韩 | 一进一出抽搐动态| 亚洲专区国产一区二区| 岛国视频午夜一区免费看| 老鸭窝网址在线观看| 国产欧美日韩一区二区精品| 欧洲精品卡2卡3卡4卡5卡区| 黄色毛片三级朝国网站| 一a级毛片在线观看| 精品国产美女av久久久久小说| av福利片在线| 亚洲免费av在线视频| 操出白浆在线播放| 白带黄色成豆腐渣| 国产一区二区激情短视频| 亚洲av电影在线进入| 国产91精品成人一区二区三区| 美女黄网站色视频| 最近最新中文字幕大全免费视频| 亚洲av熟女| 精品一区二区三区av网在线观看| 国产亚洲av高清不卡| 男女做爰动态图高潮gif福利片| 国产91精品成人一区二区三区| 亚洲国产精品成人综合色| 50天的宝宝边吃奶边哭怎么回事| 成人国产综合亚洲| 欧美黄色片欧美黄色片| 精品久久蜜臀av无| 老司机福利观看| 男人的好看免费观看在线视频 | 村上凉子中文字幕在线| а√天堂www在线а√下载| 亚洲一码二码三码区别大吗| 日日夜夜操网爽| 欧美国产日韩亚洲一区| 99久久久亚洲精品蜜臀av| 亚洲片人在线观看| 可以在线观看的亚洲视频| 亚洲欧美日韩东京热| 国产亚洲精品av在线| 亚洲国产欧美网| 波多野结衣高清作品| 老司机深夜福利视频在线观看| 天天一区二区日本电影三级| 精品高清国产在线一区| 一边摸一边做爽爽视频免费| 精品一区二区三区视频在线观看免费| 日韩成人在线观看一区二区三区| 久久久久久久午夜电影| 一个人观看的视频www高清免费观看 | 国产99久久九九免费精品| 怎么达到女性高潮| www.自偷自拍.com| 麻豆成人av在线观看| 小说图片视频综合网站| 亚洲av成人不卡在线观看播放网| 老汉色∧v一级毛片| 岛国视频午夜一区免费看| 黄色视频,在线免费观看| 国产精品久久视频播放| 久久精品91无色码中文字幕| 99精品久久久久人妻精品| 国产av一区在线观看免费| aaaaa片日本免费| 中文字幕av在线有码专区| 人人妻人人看人人澡| 女警被强在线播放| 黄频高清免费视频| 国产精品亚洲美女久久久| 在线国产一区二区在线| 欧美日韩瑟瑟在线播放| 一边摸一边做爽爽视频免费| 亚洲精品中文字幕在线视频| 欧美成人性av电影在线观看| 一本久久中文字幕| 中文字幕最新亚洲高清| 亚洲国产精品sss在线观看| 久久欧美精品欧美久久欧美| 日韩大码丰满熟妇| 男人的好看免费观看在线视频 | 熟女少妇亚洲综合色aaa.| 亚洲av成人一区二区三| 亚洲精品美女久久av网站| 亚洲av日韩精品久久久久久密| 国产黄色小视频在线观看| 国内精品久久久久精免费| 婷婷六月久久综合丁香| 18禁裸乳无遮挡免费网站照片| 不卡av一区二区三区| 久久久国产欧美日韩av| 精华霜和精华液先用哪个| 天堂动漫精品| 国产成人av教育| 国产片内射在线| 欧美成人午夜精品| 亚洲中文av在线| 一级片免费观看大全| 精品久久久久久久久久免费视频| 国内精品一区二区在线观看| 欧美日韩中文字幕国产精品一区二区三区| 丰满人妻一区二区三区视频av | 亚洲电影在线观看av| 亚洲国产欧美人成| 欧美一级a爱片免费观看看 | 国产精品一区二区三区四区免费观看 | 激情在线观看视频在线高清| 久久精品国产清高在天天线| 亚洲熟女毛片儿| 久久久久久久午夜电影| 精品国产乱子伦一区二区三区| 999精品在线视频| 一边摸一边抽搐一进一小说| 亚洲av成人一区二区三| 亚洲av中文字字幕乱码综合| 国产野战对白在线观看| 最近在线观看免费完整版| 麻豆国产av国片精品| 久久国产精品影院| 久久精品国产亚洲av香蕉五月| 欧美一级毛片孕妇| 此物有八面人人有两片| 午夜福利在线观看吧| 国产av又大| 国产精品亚洲一级av第二区| 国产一区在线观看成人免费| 一边摸一边做爽爽视频免费| 国产精品电影一区二区三区| av免费在线观看网站| 性色av乱码一区二区三区2| 两人在一起打扑克的视频| 成人永久免费在线观看视频| 精品久久久久久,| 日韩大码丰满熟妇| 国产成人啪精品午夜网站| 亚洲精品一区av在线观看| 麻豆国产av国片精品| 精品欧美一区二区三区在线| 不卡一级毛片| 美女 人体艺术 gogo| 伦理电影免费视频| 国产三级在线视频| 黄色片一级片一级黄色片| 成人特级黄色片久久久久久久| 亚洲精品中文字幕一二三四区| 欧美一级毛片孕妇| 青草久久国产| 亚洲国产精品成人综合色| a级毛片在线看网站| 亚洲欧美激情综合另类| 国产日本99.免费观看| 老司机福利观看| 精品久久久久久久久久久久久| 亚洲中文字幕一区二区三区有码在线看 | 国产精品久久电影中文字幕| 成人特级黄色片久久久久久久| 国产欧美日韩精品亚洲av| x7x7x7水蜜桃| 欧美精品亚洲一区二区| 亚洲国产欧美人成| 最近最新免费中文字幕在线| 午夜激情福利司机影院| 国产乱人伦免费视频| www.www免费av| √禁漫天堂资源中文www| 日韩免费av在线播放| 亚洲成人久久爱视频| 国产aⅴ精品一区二区三区波| 日本 av在线| 国产麻豆成人av免费视频| 欧美一区二区精品小视频在线| 一边摸一边抽搐一进一小说| 亚洲色图av天堂| 国产精品久久久久久亚洲av鲁大| 啦啦啦免费观看视频1| 成年免费大片在线观看| 日韩欧美在线二视频| 美女免费视频网站| 亚洲欧美精品综合久久99| 精品一区二区三区av网在线观看| 2021天堂中文幕一二区在线观| 男人舔女人下体高潮全视频| 丰满的人妻完整版| 丰满人妻熟妇乱又伦精品不卡| 欧美激情久久久久久爽电影| 精品久久久久久久末码| 国产乱人伦免费视频| avwww免费| 欧美3d第一页| 日韩欧美国产一区二区入口| 久久久久久大精品| 香蕉av资源在线| 国产aⅴ精品一区二区三区波| 亚洲成人国产一区在线观看| 亚洲片人在线观看| 美女扒开内裤让男人捅视频| 亚洲精品久久成人aⅴ小说| 少妇的丰满在线观看| 亚洲精品中文字幕在线视频| 日日夜夜操网爽| 99精品在免费线老司机午夜| 亚洲午夜精品一区,二区,三区| 亚洲七黄色美女视频| 日韩大码丰满熟妇| 国产精品永久免费网站| 欧美日韩瑟瑟在线播放| 大型黄色视频在线免费观看| 日韩成人在线观看一区二区三区| 三级毛片av免费| 日韩精品青青久久久久久| 五月玫瑰六月丁香| 可以免费在线观看a视频的电影网站| 桃色一区二区三区在线观看| 少妇粗大呻吟视频| 亚洲欧美日韩高清专用| 不卡一级毛片| 哪里可以看免费的av片| 色尼玛亚洲综合影院| 欧美乱码精品一区二区三区| 国产精品日韩av在线免费观看| 五月玫瑰六月丁香| 天天添夜夜摸| 免费一级毛片在线播放高清视频| 国产真人三级小视频在线观看| 啪啪无遮挡十八禁网站| 夜夜夜夜夜久久久久| 久久人妻av系列| 亚洲激情在线av| 狠狠狠狠99中文字幕| 一级作爱视频免费观看| 国产精品 欧美亚洲| 久久久久久免费高清国产稀缺| 久久人人精品亚洲av| av在线播放免费不卡| 国产成人精品久久二区二区91| 看片在线看免费视频| 亚洲熟妇中文字幕五十中出| 怎么达到女性高潮| 女人爽到高潮嗷嗷叫在线视频| 国产伦在线观看视频一区| 99久久久亚洲精品蜜臀av| 国产精品 国内视频| 久9热在线精品视频| 久久久精品大字幕| 亚洲美女黄片视频| 亚洲 欧美一区二区三区| 国产v大片淫在线免费观看| 老司机在亚洲福利影院| 夜夜爽天天搞| 国产精品国产高清国产av| 2021天堂中文幕一二区在线观| 免费在线观看亚洲国产| 国产久久久一区二区三区| 又紧又爽又黄一区二区| 丰满人妻一区二区三区视频av | av国产免费在线观看| 男女午夜视频在线观看| 淫妇啪啪啪对白视频| 国产成人啪精品午夜网站| 亚洲国产高清在线一区二区三| 精品少妇一区二区三区视频日本电影| 哪里可以看免费的av片| 一进一出抽搐动态| 亚洲 欧美一区二区三区| 国产精品久久视频播放| 国产一区二区三区在线臀色熟女| 亚洲av成人av| 一区福利在线观看| 婷婷六月久久综合丁香| 他把我摸到了高潮在线观看| 99精品久久久久人妻精品| 手机成人av网站| 村上凉子中文字幕在线| 亚洲五月天丁香| 国产成人一区二区三区免费视频网站| 国产v大片淫在线免费观看| www.熟女人妻精品国产| 一进一出好大好爽视频| 成在线人永久免费视频| 亚洲五月婷婷丁香| 日日干狠狠操夜夜爽| 亚洲国产中文字幕在线视频| 超碰成人久久| 免费看美女性在线毛片视频| 国产精品乱码一区二三区的特点| 伊人久久大香线蕉亚洲五| 99久久精品热视频| 国产欧美日韩一区二区精品|