t.c

Go to the documentation of this file.

int t()
{
  
  for(Int_t i=0;i<339;i)
  {
    TFile* f1 = new TFile("ppdc.root");
    TH1F* h1 = (TH1F*)f1->Get("htemp");
    //set parameters
    Double_t pMag = h1->GetRandom();
    cout<<"pMag: "<<pMag<<endl;
    f1->Close();
 
    TFile* f2 = new TFile("theta.root");
    TH1F* h2 = (TH1F*)f2->Get("htemp");
    //randomize cos(theta)
    Double_t theta =(Double_t)h2->GetRandom();
    cout<<"theta: "<<theta<<endl;
    G4double costheta = cos(theta);
    f2->Close();
 
    TFile* f3 = new TFile("phi.root");
    TH1F* h3 = (TH1F*)f3->Get("htemp");
    //randomize phi
    Double_t phi = (Double_t)h3->GetRandom();
    cout<<"phi: "<<phi<<endl;
    f3->Close();
  }
}