Calculating predicted probabilities using a Cox model
There is a way of obtaining prediction out of a Cox model, as survival probability at time $t$ ($S(t)$) depends on your cox model like so:
$S(t) = e^{-H_0(t) * exp(LP)}$
in this formula $H_0(t)$ is called the baseline hazard at time $t$; and $LP$ is the linear predictor.
If $X_i$ are the predictor variables $1,2,...,i$ in the Cox model, and $β_i$ are the corresponding coefficients from the Cox model, then the linear predictor is calculated like so:
$LP = X_1*β_1 + X_2*β_2 + ... + X_i*β_i$
You might already be familiar with what a linear predictor is, but I added this for clarity sake.
The baseline hazard is a little harder to obtain. Basically it is a function of time, which shows the hazard for an event of an individual who has a $LP$ of 0. Due to the way cox regression works this value is not estimated (look at this CV question or others, which sheds more light on baseline hazard). In R however (which I can see the OP uses), there is a function available called 'basehaz()' in the survival package, which will let you extract the baseline hazard at a specific timepoint based on your model's fit.
If you extract this baseline hazard you can complete the formula above for any individual in your data, and for unseen data, as long as you are able to calculate the linear predictor for these individuals.
Remember this formula results in the probability of survival or not having had the event at time $t$. If you want to know the probability of an event at time $t$, simply subtract the probability from 1:
$P(event|LP) = 1-S(t) = 1 - e^{-H_0(t) * exp(LP)}$
Final remark: Note that the basehaz function gives you a baseline hazard at time $t$ based on your specific data. So, as is the case for the coefficients, extrapolating this to new cases might not result in a good fit/prediction, due to overfitting, bias, etc.
Calibration & Discrimination
With the predicted probabilities and the event status at one year you can calculate various calibration statistics.
For the discrimination statistic (the c-index) there is a specific kind of rank correlation for a censored response variable. In R, the rcorr.cens function from the Hmisc package can provide you with it.
Whether these statistics are the most suitable for your research I can not say, as that depends mostly on the specifics (e.g. are you building a new and first model for this setting, or are you comparing to previously developed ones?)