In this paper, we propose a mathematical model to assess the impacts of using face masks, hospitalization of symptomatic individuals and quarantine of asymptomatic individuals in combating the COVID-19 pandemic in India. We calibrate the proposed model to fit the four data sets, viz. data for the states of Maharashtra, Delhi, Tamil Nadu and overall India, and estimate the rate of infection of susceptible with symptomatic population and recovery rate of quarantined individuals. We also estimate basic reproduction number to illustrate the epidemiological status of the regions under study. Our simulations infer that the infective population will be on increasing curve for Maharashtra and India, and settling for Tamil Nadu and Delhi. Sophisticated techniques of sensitivity analysis are employed to determine the impacts of model parameters on basic reproduction number and symptomatic infected individuals. Our results reveal that to curtail the disease burden in India, specific control strategies should be implemented effectively so that the basic reproduction number is decreased below unity. The three control strategies are shown to be important preventive measures to lower disease transmission rate. The model is further extended to its stochastic counterpart to encapsulate the variation or uncertainty observed in the disease transmissibility. We observe the variability in the infective population and found their distribution at certain fixed time, which shows that for small populations, the stochasticity will play an important role. c 2021 the Author(s).