Fork me on GitHub

Changeset 190 in svn for trunk


Ignore:
Timestamp:
Jan 26, 2009, 3:37:43 PM (16 years ago)
Author:
Xavier Rouby
Message:

howto for pile-up added

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/README

    r168 r190  
     1** travail en cours**
     2 -> faire un  bfield plus rapide
     3 -> Track : x,y,z,eta,phi : vertex
     4          : x,y,z,eta,phi : calo
     5 -> les constructeurs de BlockClasses sont trÚs mal écrits!
     6 -> [11:05:29] LoÃŒst Un Ange says:
     7        Warning in <TVector3::PseudoRapidity>: transvers momentum = 0! return +/- 10e10
     8        ca lui arrive 6 fois, je pense que c'est pas grave mais je ne sais pas si pour vous il faut y faire qq choe
     9 -> TrackPropagation::TrackPropagation(const string DetDatacard) {
     10      if(DetDatacard="") { }; ///// TODO
     11 -> faire un "default" dans SmearUtil, qui sera appelé par le constructeur par défaut, si pas de data card donnée. Changer le ReadDatacard en constructeur
     12
     13
    114To compile :
    215 *** ./genMakeFile.tcl > Makefile
     
    922 - mettre une carte des efficacités en eta/phi ?
    1023 - quid de la construction du détecteur comme un assemblage de petites classes (comme un lego ou comme le beamline d'Hector)
    11 
    12 /// !!!
    13 * pour les temps de vol, on suppose que les particules sont hautement relativistes = elles voyages à la vitesse de la lumiÚre
    1424
    1525Pour les legos:
     
    4858 ligne 60 de SmearUtil.cc
    4959
    50 * Delphes.cpp: ligne 196:
    51   for(entry = 0; entry < 200; ++entry)
    52   -> 200???
     60
     61code pile-up FAMOS
     62 // How many pile-up events?
     63174  int PUevts = (int) random->poissonShoot(averageNumber_);
     64175 
     65176  // Get N events from random files
     66177  for ( int ievt=0; ievt<PUevts; ++ievt ) { 
     67178     
     68179    // Draw a file in a ramdom manner 
     69180    unsigned file = (unsigned) (theNumberOfFiles * random->flatShoot());
     70181    /*
     71182    if ( debug ) 
     72183      std::cout << "The file chosen for event " << ievt 
     73184                << " is the file number " << file << std::endl; 
     74185    */
     75186 
     76187    // Smear the primary vertex and express it in mm (stupid GenEvent convention...)
     77188    theVertexGenerator->generate();
     78189    HepMC::FourVector smearedVertex = 
     79190      HepMC::FourVector(theVertexGenerator->X()*10.,
     80191                        theVertexGenerator->Y()*10.,
     81192                        theVertexGenerator->Z()*10.,
     82193                        0.);
     83194    HepMC::GenVertex* aVertex = new HepMC::GenVertex(smearedVertex);
     84195    evt->add_vertex(aVertex);
     85196 
     86197    // Some rotation around the z axis, for more randomness
     87198    double theAngle = random->flatShoot() * 2. * 3.14159265358979323;
     88199    double cAngle = std::cos(theAngle);
     89200    double sAngle = std::sin(theAngle);
     90201 
Note: See TracChangeset for help on using the changeset viewer.