Changeset 190 in svn for trunk

Timestamp:

Jan 26, 2009, 3:37:43 PM (16 years ago)

Author:

Xavier Rouby

Message:

howto for pile-up added

File:

• trunk/README (modified) (3 diffs)

trunk/README

@@ +1 @@
+** travail en cours**
+ -> faire un  bfield plus rapide
+ -> Track : x,y,z,eta,phi : vertex
+          : x,y,z,eta,phi : calo
+ -> les constructeurs de BlockClasses sont trÚs mal écrits!
+ -> [11:05:29] LoÃŒst Un Ange says: 
+        Warning in <TVector3::PseudoRapidity>: transvers momentum = 0! return +/- 10e10
+        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
+ -> TrackPropagation::TrackPropagation(const string DetDatacard) { 
+      if(DetDatacard="") { }; ///// TODO
+ -> 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 
+
+
 To compile :
  *** ./genMakeFile.tcl > Makefile
@@ -9 +22 @@
  - mettre une carte des efficacités en eta/phi ?
  - quid de la construction du détecteur comme un assemblage de petites classes (comme un lego ou comme le beamline d'Hector)
-
-/// !!! 
-* pour les temps de vol, on suppose que les particules sont hautement relativistes = elles voyages à la vitesse de la lumiÚre
 
 Pour les legos:
@@ -48 +58 @@
  ligne 60 de SmearUtil.cc
 
-* Delphes.cpp: ligne 196:
-  for(entry = 0; entry < 200; ++entry)
-  -> 200???
+
+code pile-up FAMOS
+ // How many pile-up events? 
+174  int PUevts = (int) random->poissonShoot(averageNumber_); 
+175  
+176  // Get N events from random files 
+177  for ( int ievt=0; ievt<PUevts; ++ievt ) {  
+178      
+179    // Draw a file in a ramdom manner  
+180    unsigned file = (unsigned) (theNumberOfFiles * random->flatShoot()); 
+181    /* 
+182    if ( debug )  
+183      std::cout << "The file chosen for event " << ievt  
+184                << " is the file number " << file << std::endl;  
+185    */ 
+186  
+187    // Smear the primary vertex and express it in mm (stupid GenEvent convention...) 
+188    theVertexGenerator->generate(); 
+189    HepMC::FourVector smearedVertex =  
+190      HepMC::FourVector(theVertexGenerator->X()*10., 
+191                        theVertexGenerator->Y()*10., 
+192                        theVertexGenerator->Z()*10., 
+193                        0.); 
+194    HepMC::GenVertex* aVertex = new HepMC::GenVertex(smearedVertex); 
+195    evt->add_vertex(aVertex); 
+196  
+197    // Some rotation around the z axis, for more randomness 
+198    double theAngle = random->flatShoot() * 2. * 3.14159265358979323; 
+199    double cAngle = std::cos(theAngle); 
+200    double sAngle = std::sin(theAngle); 
+201