Skip to content

Instantly share code, notes, and snippets.

View pung96's full-sized avatar

Beomsu Chang pung96

  • University of Jyvaskyla
View GitHub Profile
username: '[email protected]'
password: ''
album_id: ''
const double kJPi = TMath::Pi();
const double kJTwoPi = kJPi*2;
void dijetcor(){
TH1D *hCor = new TH1D ( "hCor", "hCor", 300, -1/3.0, ( 2-1/3.0 ) );
for( int ievt=0;ievt<1000;ievt++ ){
double phi1 = gRandom->Uniform( kJTwoPi );
//double phi2 = phi1 + kJPi ;
double phi2 = gRandom->Gaus( phi1 + kJPi, 0.3 ) ;
double dphi1 = 0.1*kJPi;