various_assignments/SimPa/CallSimulatorT5zz.C

#include "Modules/SimulatorT5zz.C"
#include <string>

int main(int argc, char *argv[]) {
  
  float massGluHi=1200;
  float massGluStep=25;
  
  float massLSP = 200;
  
  float xlow=0.1;
  float xstep=0.05;
  float xhi=0.99;
  
  size_t NeventsPerPoint=20000;
  
  for(float massGlu = massGluHi; massGlu > massLSP + 100; massGlu -= massGluStep) {
    for(float x = xlow; x < xhi; x += xstep) {
      float massChi = massLSP+x*(massGlu-massLSP);
      if(massChi < massLSP + 100.0 || massGlu < massChi + 2) continue;
      char outputname[50];
      int massChi_ = (int)massChi;
      int massGlu_ = (int)massGlu;
      int massLSP_ = (int)massLSP;
      sprintf(outputname, "output/output_%d_%d_%d.root", massGlu_, massChi_, massLSP_);
      cout << outputname << endl;
      SimulatorT5zz *a = new SimulatorT5zz(NeventsPerPoint, outputname, massGlu, massLSP, x);
      delete a;
    }
  }

  return 0;
}
Revision:	1.1
Committed:	Tue Apr 10 14:38:42 2012 UTC (13 years ago) by pablom
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Log Message:	SimPa implemented
#	User	Rev	Content
1	pablom	1.1	#include "Modules/SimulatorT5zz.C"
2			#include <string>
3
4			int main(int argc, char *argv[]) {
5
6			float massGluHi=1200;
7			float massGluStep=25;
8
9			float massLSP = 200;
10
11			float xlow=0.1;
12			float xstep=0.05;
13			float xhi=0.99;
14
15			size_t NeventsPerPoint=20000;
16
17			for(float massGlu = massGluHi; massGlu > massLSP + 100; massGlu -= massGluStep) {
18			for(float x = xlow; x < xhi; x += xstep) {
19			float massChi = massLSP+x*(massGlu-massLSP);
20			if(massChi < massLSP + 100.0 \|\| massGlu < massChi + 2) continue;
21			char outputname[50];
22			int massChi_ = (int)massChi;
23			int massGlu_ = (int)massGlu;
24			int massLSP_ = (int)massLSP;
25			sprintf(outputname, "output/output_%d_%d_%d.root", massGlu_, massChi_, massLSP_);
26			cout << outputname << endl;
27			SimulatorT5zz *a = new SimulatorT5zz(NeventsPerPoint, outputname, massGlu, massLSP, x);
28			delete a;
29			}
30			}
31
32			return 0;
33			}