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, 1 month ago) by pablom
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Error occurred while calculating annotation data.
Log Message:	SimPa implemented
#	Content
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	}