Fork me on GitHub

Changeset 404 in svn for trunk/src/SmearUtil.cc


Ignore:
Timestamp:
May 20, 2009, 10:38:03 AM (16 years ago)
Author:
Xavier Rouby
Message:

new parameters for RomanPots in datacard ; print_header added

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/SmearUtil.cc

    r399 r404  
    194194  //for very forward detectors
    195195  RP_offsetEl_s     = 120;    // distance of beam separation point, from IP
    196   RP_offsetEl_x     = 0.097;
     196  RP_offsetEl_x     = -0.097; // half distance of separation of beams in horizontal plan, in m
     197  RP_offsetEl_y     = 0;      // half distance of separation of beams in vertical plan, in m
    197198  RP_cross_x        = -500;   // IP offset in horizontal plane, in micrometers
    198199  RP_cross_y        = 0.0;    // IP offset in vertical plane, in micrometers
     
    312313  RP_offsetEl_s     = DET.RP_offsetEl_s;
    313314  RP_offsetEl_x     = DET.RP_offsetEl_x;
     315  RP_offsetEl_y     = DET.RP_offsetEl_y;
    314316  RP_cross_x        = DET.RP_cross_x;
    315317  RP_cross_y        = DET.RP_cross_y;
     
    449451  RP_offsetEl_s     = DET.RP_offsetEl_s;
    450452  RP_offsetEl_x     = DET.RP_offsetEl_x;
     453  RP_offsetEl_y     = DET.RP_offsetEl_y;
    451454  RP_beam1Card      = DET.RP_beam1Card;
    452455  RP_beam2Card      = DET.RP_beam2Card;
     
    529532    else if(strstr(temp_string.c_str(),"RP_offsetEl_s"))    {curstring >> varname >> value; RP_offsetEl_s     = value;}
    530533    else if(strstr(temp_string.c_str(),"RP_offsetEl_x"))    {curstring >> varname >> value; RP_offsetEl_x     = value;}
     534    else if(strstr(temp_string.c_str(),"RP_offsetEl_y"))    {curstring >> varname >> value; RP_offsetEl_y     = value;}
    531535    else if(strstr(temp_string.c_str(),"RP_cross_x"))       {curstring >> varname >> value; RP_cross_x        = value;}
    532536    else if(strstr(temp_string.c_str(),"RP_cross_y"))       {curstring >> varname >> value; RP_cross_y        = value;}
     
    724728    f_out << left << setw(35) <<"* Datacard for beam 2: "<<""
    725729          << left << setw(25) <<RP_beam2Card            <<""<< right << setw(10)<<"*"<<"\n";
    726     f_out << left << setw(44) <<"* Beam separation, in meters:     "<<""
     730    f_out << left << setw(44) <<"* Beam separation, in meters(hor):"<<""
    727731          << left << setw(6) << RP_offsetEl_x          <<""<< right << setw(10)<<"*"<<"\n";
     732    f_out << left << setw(44) <<"* Beam separation, in meters(ver):"<<""
     733          << left << setw(6) << RP_offsetEl_y          <<""<< right << setw(10)<<"*"<<"\n";
    728734    f_out << left << setw(44) <<"* Distance from IP for Beam separation (m):"<<""
    729735          << left << setw(6) <<RP_offsetEl_s           <<""<< right << setw(10)<<"*"<<"\n";
     
    14041410    // ID name   chg       mass    total width   lifetime
    14051411    //  1 d      -1      0.33000     0.00000   0.00000E+00
     1412    //  in the table, the charge is in units of e+/3
     1413    //  the total width is in GeV
     1414    //  the lifetime is ctau in mm
    14061415    curstring >> ID >> name >> charge >> mass >> width >> lifetime;
    1407     PdgParticle particle(ID,name,mass,charge/3.,width,lifetime);
     1416    PdgParticle particle(ID,name,mass,charge/3.,width,lifetime/1000.);
    14081417    PDGtable.insert(ID,particle);
    14091418    //PdgTable.insert(pair<int,PdgParticle>(ID,particle));
     
    15221531
    15231532}
     1533
     1534void print_header() {
     1535  cout << endl << endl;
     1536 
     1537  cout <<"*********************************************************************"<< endl;
     1538  cout <<"*********************************************************************"<< endl;
     1539  cout <<"**                                                                 **"<< endl;
     1540  cout <<"**                            Welcome to                           **"<< endl;
     1541  cout <<"**                                                                 **"<< endl;
     1542  cout <<"**                                                                 **"<< endl;
     1543  cout <<"**     .ddddddd-            lL             hH                      **"<< endl;
     1544  cout <<"**     -Dd` `dD:           Ll             hH`                      **"<< endl;
     1545  cout <<"**     dDd   dDd    eeee.  lL  .pp+pp    Hh+hhh`   -eeee- `sssss   **"<< endl;
     1546  cout <<"**    -Dd   `DD   ee. ee  Ll   .Pp. PP   Hh. HH.  ee. ee  sSs      **"<< endl;
     1547  cout <<"**    dD`   dDd  eEeee:   lL.  pP.  pP  hH   hH` eEeee:` -sSSSs.   **"<< endl;
     1548  cout <<"**   .Dd  :dd    eE.     LlL   PpppPP   Hh  Hh   eE         sSS    **"<< endl;
     1549  cout <<"**   dddddd:.     eee+:  lL.  pp.      hh.  hh    eee+  sssssS     **"<< endl;
     1550  cout <<"**                            Pp                                   **"<< endl;
     1551  cout <<"**                                                                 **"<< endl;
     1552  cout <<"**           Delphes, a framework for the fast simulation          **"<< endl;
     1553  cout <<"**                 of a  generic collider experiment               **"<< endl;
     1554  cout <<"**                    arXiv:0903.2225v1 [hep-ph]                   **"<< endl;
     1555  cout <<"**                                                                 **"<< endl;
     1556  cout <<"**                  --- Version 1.6 of Delphes ---                 **"<< endl;
     1557  cout <<"**                Last date of change:    7 May 2009               **"<< endl;
     1558  cout <<"**                                                                 **"<< endl;
     1559  cout <<"**                                                                 **"<< endl;
     1560  cout <<"**     This package uses:                                          **"<< endl;
     1561  cout <<"**     ------------------                                          **"<< endl;
     1562  cout <<"**     FastJet algorithm: Phys. Lett. B641 (2006) [hep-ph/0512210] **"<< endl;
     1563  cout <<"**     Hector: JINST 2:P09005 (2007) [physics.acc-ph:0707.1198v2]  **"<< endl;
     1564  cout <<"**     FROG:                               [hep-ex/0901.2718v1]    **"<< endl;
     1565  cout <<"**                                                                 **"<< endl;
     1566  cout <<"**-----------------------------------------------------------------**"<< endl;
     1567  cout <<"**                                                                 **"<< endl;
     1568  cout <<"**   Main authors:                                                 **"<< endl;
     1569  cout <<"**   -------------                                                 **"<< endl;
     1570  cout <<"**                                                                 **"<< endl;
     1571  cout <<"**              Séverine Ovyn               Xavier Rouby           **"<< endl;
     1572  cout <<"**       severine.ovyn@uclouvain.be      xavier.rouby@cern         **"<< endl;
     1573  cout <<"**       Center for Particle Physics and Phenomenology (CP3)       **"<< endl;
     1574  cout <<"**       Universite Catholique de Louvain (UCL)                    **"<< endl;
     1575  cout <<"**       Louvain-la-Neuve, Belgium                                 **"<< endl;
     1576  cout <<"**                                                                 **"<< endl;
     1577  cout <<"**-----------------------------------------------------------------**"<< endl;
     1578  cout <<"**                                                                 **"<< endl;
     1579  cout <<"**   Former Delphes versions and documentation can be found on :   **"<< endl;
     1580  cout <<"**               http://www.fynu.ucl.ac.be/delphes.html            **"<< endl;
     1581  cout <<"**                                                                 **"<< endl;
     1582  cout <<"**                                                                 **"<< endl;
     1583  cout <<"**   Disclaimer: this program is a beta version of Delphes and     **"<< endl;
     1584  cout <<"** therefore comes without guarantees. Beware of errors and please **"<< endl;
     1585  cout <<"**        give us your feedbacks about potential bugs              **"<< endl;
     1586  cout <<"**                                                                 **"<< endl;
     1587  cout <<"*********************************************************************"<< endl;
     1588  cout <<"*********************************************************************"<< endl;
     1589
     1590}
Note: See TracChangeset for help on using the changeset viewer.