OS notes

Status making tables

• Jake is cleaning up code and will commit to CVS
• will wait for him to finish, then make tables with him and put on twiki

FR included in tables

• FR prediction is included in looper
• need to add results to "table-maker"

DY prediction

• learn how to do it from Dave
• what are the restrictions in our special corner of phase space

Status Significance estimate

Look into Z_Bi

• looked at talk - too little info, need to read paper - arXiv: physics/0702156v4 , pub in Nuc. Ins. and Meth in Phys. Research A 595 (2008) 480-501 arXiv link, Jeff R. at CMS week: talk
• example code:

  {
    double n_on     = 140.; // total events in signal region (S+B)
    double mu_b_hat = 83.33;  // mean num of BG events expected in sig. region
    double sigma_b  = 8.333; // uncertainty of mu_b_hat
    double tau      = mu_b_hat / (sigma_b*sigma_b); // scale factor to corresp. Noff/Non
    double n_off    = tau*mu_b_hat;
    double P_Bi     = TMath::BetaIncomplete(1./(1.+tau), n_on, n_off+1);
    double Z_Bi     = sqrt(2.)*TMath::ErfInverse(1 - 2.*P_Bi);
  
    cout  <<"  n_on     " <<n_on      <<endl
          <<"  mu_b_hat " <<mu_b_hat  <<endl
          <<"  sigma_b  " <<sigma_b   <<endl
          <<"  tau      " <<tau       <<endl
          <<"  n_off    " <<n_off     <<endl
          <<"  P_Bi     " <<P_Bi      <<endl
          <<"  Z_Bi     " <<Z_Bi      <<endl;
  
  
  }
  
      

Generation of CMSSW / Tas tuples for the LHA files

SLHA++ part

1. SLHA files have spectroscopy. couplings, mixing, branching fractions link to files: srm://bsrm-1.t2.ucsd.edu:8443/srm/v2/server?SFN=/hadoop/cms/store/user/spadhi/SLHA/msugra_10_0_0_shla.tar.gz
2. "add pdfs" for pp xsecs, done with PYTHIA for all points => each SLHA file gets one number (xsec) link to xecs: png plot, Sanjay will send mail with txt values.
3. see below

CMSSW / TAS part

1. FastSim / TaSntuple config:
   • one config for entire chain from reading the SLHA file, reading the xsecs and producing the tas tuples - see below.

script magic to send jobs to crab

1. at which point Sanjay remarks that he is preparing to prepare RECO files for each of the SLHE files

Branching point: -> Either:

1. wait for Sanjay to have produced RECO files -> he says deadline is Saturday March 27th.
2. run std TaS ntuple jobs on the RECO files in local DBS

-> OR:

1. set up scripting magic
2. run FastSim + TaS
3. postprocess and store

-- IngoBloch - 19-Mar-2010