1 | samples higgs tautau:
|
---|
2 | scp higgs*root h*tau*root server.fynu.ucl.ac.be:/nfs/cms/mass8/CMSDATA/CMSSW/root_files_delphes/
|
---|
3 |
|
---|
4 |
|
---|
5 | BUG:
|
---|
6 | * erreur de compilation si dans le répertoire courant un fichier porte le nom d'un des #include <...>
|
---|
7 | par exemple, si les fichiers "list" ou "vector" existent (quel que soit le contenu)
|
---|
8 |
|
---|
9 | ** travail en cours**
|
---|
10 | -> les constructeurs de BlockClasses sont trÚs mal écrits!
|
---|
11 | -> faire un "default" dans SmearUtil, qui sera appelé par le constructeur par défaut, si pas de data card donnée.
|
---|
12 | Changer le ReadDatacard en constructeur
|
---|
13 |
|
---|
14 | *** pour sortir une nouvelle version de DELPHES
|
---|
15 | - 1) incrémenter le numéro de version dans VERSION
|
---|
16 | -1b) changer le numéro et la date dans src/Smearutil, dans les cout
|
---|
17 | - 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!
|
---|
18 | - 3) mettre sur la page web + change logs
|
---|
19 | - 4) cvs tag Delphes_V_x-x
|
---|
20 |
|
---|
21 |
|
---|
22 | * check all resolution number and put some reference on them
|
---|
23 | * 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?
|
---|
24 |
|
---|
25 | delires :
|
---|
26 | - tenir compte des longueurs de radiation ?
|
---|
27 | - mettre une carte des efficacités en eta/phi ?
|
---|
28 | - quid de la construction du détecteur comme un assemblage de petites classes (comme un lego ou comme le beamline d'Hector)
|
---|
29 |
|
---|
30 | Pour les legos:
|
---|
31 |
|
---|
32 | chaque sous détecteur
|
---|
33 | - quel type
|
---|
34 | - couverture en eta
|
---|
35 | - résolution
|
---|
36 | -> ok pour calos et trackers
|
---|
37 | -> quid des RP ?
|
---|
38 |
|
---|
39 |
|
---|
40 | CASTOR arXiv:0706.2641
|
---|
41 | p0 = 1059
|
---|
42 |
|
---|
43 | ZDC = linearity < 2%
|
---|
44 | \sigma/E = 0.7/\sqrt(E) \oplus 0.08 <- positrons
|
---|
45 | \sigma/E = 1.38/\sqrt(E) + 0.13 <- pions
|
---|
46 |
|
---|
47 | stdhep.ps -> page9 = HEPCHG(ID) : get charge for particle ID
|
---|
48 |
|
---|
49 | ne pas oublier de vérifier que la fonction SmearHadron doit avoir un sqrt à la fin
|
---|
50 |
|
---|
51 | si on a une liste de fichier input avec différents types de fichiers (root et hep par exemple), ça va crasher !
|
---|
52 |
|
---|
53 | attention à CalTower.hh qui est hyper spécifiqueà CMS, aprÚs les changements. comment introduire ça dans la carte?
|
---|
54 |
|
---|
55 |
|
---|
56 | * attention : Trigger::
|
---|
57 | * void TriggerBit::PrintTrigCondition(string LogName,int i)
|
---|
58 | --> il y a a déjà un i défini ailleurs dans la fonction
|
---|
59 | * 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
|
---|
60 |
|
---|
61 | * dans RESOLution
|
---|
62 | TRACK_eff = 100; ???
|
---|
63 | ligne 60 de SmearUtil.cc
|
---|
64 |
|
---|
65 |
|
---|
66 | code pile-up FAMOS
|
---|
67 | // How many pile-up events?
|
---|
68 | 174 int PUevts = (int) random->poissonShoot(averageNumber_);
|
---|
69 | 175
|
---|
70 | 176 // Get N events from random files
|
---|
71 | 177 for ( int ievt=0; ievt<PUevts; ++ievt ) {
|
---|
72 | 178
|
---|
73 | 179 // Draw a file in a ramdom manner
|
---|
74 | 180 unsigned file = (unsigned) (theNumberOfFiles * random->flatShoot());
|
---|
75 | 181 /*
|
---|
76 | 182 if ( debug )
|
---|
77 | 183 std::cout << "The file chosen for event " << ievt
|
---|
78 | 184 << " is the file number " << file << std::endl;
|
---|
79 | 185 */
|
---|
80 | 186
|
---|
81 | 187 // Smear the primary vertex and express it in mm (stupid GenEvent convention...)
|
---|
82 | 188 theVertexGenerator->generate();
|
---|
83 | 189 HepMC::FourVector smearedVertex =
|
---|
84 | 190 HepMC::FourVector(theVertexGenerator->X()*10.,
|
---|
85 | 191 theVertexGenerator->Y()*10.,
|
---|
86 | 192 theVertexGenerator->Z()*10.,
|
---|
87 | 193 0.);
|
---|
88 | 194 HepMC::GenVertex* aVertex = new HepMC::GenVertex(smearedVertex);
|
---|
89 | 195 evt->add_vertex(aVertex);
|
---|
90 | 196
|
---|
91 | 197 // Some rotation around the z axis, for more randomness
|
---|
92 | 198 double theAngle = random->flatShoot() * 2. * 3.14159265358979323;
|
---|
93 | 199 double cAngle = std::cos(theAngle);
|
---|
94 | 200 double sAngle = std::sin(theAngle);
|
---|