Fork me on GitHub

source: svn/trunk/README@ 227

Last change on this file since 227 was 206, checked in by Xavier Rouby, 16 years ago

howto publish a new version of DELPHES (tar.gz, website, CVS)

File size: 3.9 KB
RevLine 
[190]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
[206]13*** pour sortir une nouvelle version de DELPHES
14 - 1) incrémenter le numéro de version dans VERSION
15 - 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!
16 - 3) mettre sur la page web + change logs
17 - 4) cvs tag Delphes_V_x-x
[190]18
[2]19
20* check all resolution number and put some reference on them
21* 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?
22
[168]23delires :
[2]24 - tenir compte des longueurs de radiation ?
25 - mettre une carte des efficacités en eta/phi ?
26 - quid de la construction du détecteur comme un assemblage de petites classes (comme un lego ou comme le beamline d'Hector)
27
[7]28Pour les legos:
29
30chaque sous détecteur
31 - quel type
32 - couverture en eta
33 - résolution
34-> ok pour calos et trackers
35-> quid des RP ?
36
37
38CASTOR arXiv:0706.2641
39 p0 = 1059
40
41ZDC = linearity < 2%
42 \sigma/E = 0.7/\sqrt(E) \oplus 0.08 <- positrons
43 \sigma/E = 1.38/\sqrt(E) + 0.13 <- pions
44
45stdhep.ps -> page9 = HEPCHG(ID) : get charge for particle ID
46
[20]47ne pas oublier de vérifier que la fonction SmearHadron doit avoir un sqrt à la fin
[28]48
49si on a une liste de fichier input avec différents types de fichiers (root et hep par exemple), ça va crasher !
[36]50
51attention à CalTower.hh qui est hyper spécifiqueà CMS, aprÚs les changements. comment introduire ça dans la carte?
[71]52
53
[168]54* attention : Trigger::
55 * void TriggerBit::PrintTrigCondition(string LogName,int i)
56 --> il y a a déjà un i défini ailleurs dans la fonction
57 * 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
[71]58
[168]59* dans RESOLution
60 TRACK_eff = 100; ???
61 ligne 60 de SmearUtil.cc
[71]62
[190]63
64code pile-up FAMOS
65 // How many pile-up events?
66174 int PUevts = (int) random->poissonShoot(averageNumber_);
67175
68176 // Get N events from random files
69177 for ( int ievt=0; ievt<PUevts; ++ievt ) {
70178
71179 // Draw a file in a ramdom manner
72180 unsigned file = (unsigned) (theNumberOfFiles * random->flatShoot());
73181 /*
74182 if ( debug )
75183 std::cout << "The file chosen for event " << ievt
76184 << " is the file number " << file << std::endl;
77185 */
78186
79187 // Smear the primary vertex and express it in mm (stupid GenEvent convention...)
80188 theVertexGenerator->generate();
81189 HepMC::FourVector smearedVertex =
82190 HepMC::FourVector(theVertexGenerator->X()*10.,
83191 theVertexGenerator->Y()*10.,
84192 theVertexGenerator->Z()*10.,
85193 0.);
86194 HepMC::GenVertex* aVertex = new HepMC::GenVertex(smearedVertex);
87195 evt->add_vertex(aVertex);
88196
89197 // Some rotation around the z axis, for more randomness
90198 double theAngle = random->flatShoot() * 2. * 3.14159265358979323;
91199 double cAngle = std::cos(theAngle);
92200 double sAngle = std::sin(theAngle);
93201
[197]94
95
96In src/BField
97-> la pseudorapidité est calculée à la main via
98 Eta_t = log( (p+Pz_t)/(p-Pz_t) )/2.;
99et on en déduit Theta_t
100 if(p>0) Theta_t = acos(Pz_t/p); else Theta_t=0;
101
102
103Todo:
104 - dans ParticleUtil
105 eta,phi et eta_calo,phi_calo
106
Note: See TracBrowser for help on using the repository browser.