Coulomb-corelated pairs

{
  gROOT->ProcessLine(".x loadlibs.cxx");

  FILE *f = fopen("job","r");
  if (!f || ferror(f) ) {
    fprintf( stderr, "Failed to find job numbern" );
    exit();
  };

  char bufstr[255];
  fgets( bufstr, sizeof(bufstr), f );
  int job = strtol( bufstr, 0, 10 );
  fclose(f);

  fprintf( stderr, "Job %i started...n", job );

  gRandom = new TRandom3( job );

  TGenVertex *vertex = new TGenVertex( 97.6e-4 );

  TF1 *f1 = new TF1( "PlabforCC", 
                     "x^2*exp(-1.1*x)*(1-exp(-0.4066-0.2864*x))",
		     2.5, 10.0 );
  TGenPlabTF1 *Plab = new TGenPlabTF1(f1);

  TGenQcmsCoulomb *Qcms = new TGenQcmsCoulomb( 0.13957018, 0.13957018 );
  Qcms->fselect = kselQregion::kQtcms | kselQregion::kQlcms;
  Qcms->SetQRange( 0.0, 0.030 );
  Qcms->SetQtRange( 0.0, 0.015 );
  Qcms->SetQlRange( -0.025, 0.025 );

  TGenPair *pair = new TGenPair( vertex, Plab, Qcms );
  pair->fcheckQcms = kTRUE;
  pair->fcheckacceptance = kTRUE;

  TKfactorCounter *kcount = new TKfactorCounter( Qcms->GetMass(1), Qcms->GetMass(2), kSPHERE );
  kcount->fselect = kselQregion::kQcms;
  kcount->SetQRange( 0.0, 0.002 );
  kcount->fifPsumtest = kTRUE;
  kcount->SetPsumRange( 3.0, 8.4 );
 
  pair->SetKfactorCounter(kcount);

  TTree *t = pair->GenXPRandom( 50000 );
  WriteGeantInputTTree( t, "kine.dat" );
  delete t;

  pair->Print();

  delete pair;
  delete kcount;
  delete vertex;
  delete Plab;
  delete Qcms;
}