** travail en cours** -> les constructeurs de BlockClasses sont très mal écrits! -> 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 -1b) changer le numéro et la date dans src/Smearutil, dans les cout - 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; ievtflatShoot()); 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);