source: svn/trunk/README@ 218

** 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 

*** pour sortir une nouvelle version de DELPHES
 - 1) incrémenter le numéro de version dans VERSION
 - 2) make tar <--> attention! si de nouveaux fichiers dans le répertoire principal doivent être inclus dans le tar.bz, il faut updater le fichier genMakefile!
 - 3) mettre sur la page web + change logs
 - 4) cvs tag Delphes_V_x-x


* check all resolution number and put some reference on them
* in smearing functions, should not be some upper value on the eta ? e.g. : if eta < 3 then it is the central calo, else it is the forward calo => put eta<5?

delires :
 - tenir compte des longueurs de radiation ?
 - 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 legos:

chaque sous détecteur
 - quel type
 - couverture en eta
 - résolution
-> ok pour calos et trackers
-> quid des RP ?


CASTOR arXiv:0706.2641
 p0 = 1059

ZDC = linearity < 2%
        \sigma/E = 0.7/\sqrt(E) \oplus 0.08    <- positrons
        \sigma/E = 1.38/\sqrt(E) + 0.13 <- pions

stdhep.ps -> page9 = HEPCHG(ID) : get charge for particle ID

ne pas oublier de vérifier que la fonction SmearHadron doit avoir un sqrt à la fin

si on a une liste de fichier input avec différents types de fichiers (root et hep par exemple), ça va crasher !

attention à CalTower.hh qui est hyper spécifiqueà CMS, aprÚs les changements. comment introduire ça dans la carte?


* attention : Trigger::
   * void TriggerBit::PrintTrigCondition(string LogName,int i)
       --> il y a a déjà un i défini ailleurs dans la fonction
   * Dans la lecture de la carte, on ne vérifie pas qu'il y a NTOWER+1 valeurs pour les TOWER_ETA_EDGE et NTOWER pour TOWER_DPHI

* dans RESOLution
 TRACK_eff         = 100;  ???
 ligne 60 de SmearUtil.cc


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  


In src/BField
-> la pseudorapidité est calculée à la main via
       Eta_t = log( (p+Pz_t)/(p-Pz_t) )/2.;
et on en déduit Theta_t 
        if(p>0) Theta_t = acos(Pz_t/p); else Theta_t=0;


Todo:
 - dans ParticleUtil
        eta,phi et eta_calo,phi_calo