Ejemplo n.º 1
0
int main (void)
{
    union {
       double result;
       unsigned long int_result;
    } u;

    double result;
    long t0, t1;
    int i;
    setStats(1);

    printf("Running %s...\n", __FILE__);
    for( int i = 0; i < 4; i++ ) {
      asm volatile("rdcycle %0":"=r"(t0));
      asm volatile("fence");
      edp_clear_cr();
      edp_prepare_dp(&input_data1, &input_data2);
      edp_run_dp(DATA_SIZE);
      asm volatile("fence");
      asm volatile("rdcycle %0":"=r"(t1));

      u.int_result = edp_read_double();
      printf("%s Run %d: Cycle Duration = %ld\n", __FILE__, i, (t1 - t0));

    }

    setStats(0);
    if( u.int_result != verify_data ) {
      printf("Error on Test %s\n", __FILE__);
      printf("Expected result = 0x%lx\n", verify_data);
      printf("Results received = 0x%lx\n", u.int_result);
    }
    return 0;
}
Ejemplo n.º 2
0
int main( int argc, char* argv[] )
{
  int results_data[DATA_SIZE];

  // Output the input array

#if HOST_DEBUG
  printArray( "input",  DATA_SIZE, input_data  );
  printArray( "verify", DATA_SIZE, verify_data );
#endif

  // If needed we preallocate everything in the caches

#if ( !HOST_DEBUG && PREALLOCATE )
  //median_asm( DATA_SIZE, input_data, results_data );
  median( DATA_SIZE, input_data, results_data );
#endif

  // Do the filter

  setStats(1);
  median( DATA_SIZE, input_data, results_data );
  setStats(0);

  // Print out the results

#if HOST_DEBUG
  printArray( "results", DATA_SIZE, results_data );
#endif

  // Check the results

  finishTest(verify( DATA_SIZE, results_data, verify_data ));

}
Ejemplo n.º 3
0
int main( int argc, char* argv[] )
{
  int results_data[DATA_SIZE];

  // Output the input array
  printArray( "input1", DATA_SIZE, input1_data );
  printArray( "input2", DATA_SIZE, input2_data );
  printArray( "verify", DATA_SIZE, verify_data );

#if PREALLOCATE
  // If needed we preallocate everything in the caches
  vvadd( DATA_SIZE, input1_data, input2_data, results_data );
#endif

  // Do the vvadd
  setStats(1);
  vvadd( DATA_SIZE, input1_data, input2_data, results_data );
  setStats(0);

  // Print out the results
  printArray( "results", DATA_SIZE, results_data );

  // Check the results
  return verify( DATA_SIZE, results_data, verify_data );
}
Ejemplo n.º 4
0
int main (void)
{

    union {
       double result;
       unsigned long int_result;
    } u;
    u.int_result = verify_data;
    double c = 0.0;
    double sum = 0.0;

    setStats(1);
    int i;
    for (i = 0; i < DATA_SIZE; i++) {
      double prod = input_data1[i]*input_data2[i];
      double y = prod - c;
      double t = sum + y;
      c = (t - sum) - y;
      sum = t;
    }
    setStats(0);

    double error = abs(sum - u.result);
    if( error < 1e-15 )
      return 0;
    else
      return 1;
}
int main (int argc, char * argv[])
{
    union {
       double result;
       unsigned long int_result;
    } u;

    u.int_result = verify_data; 
    setStats(1);
    double result = 0;
    int i;
    for (i = 0; i < DATA_SIZE; i++) {
      input_data1[i] = input_data1[i]*input_data2[i];
    }
    //Sort the products
    qsort(input_data1, DATA_SIZE, sizeof(input_data1[0]), abscomp);
    for (i = 0; i < DATA_SIZE; i++) {
        result += input_data1[i];
    }
    setStats(0);

    double error = abs(result - u.result);
    if( error < 1e-15 )
      return 0;
    else
      return 1;
}
Ejemplo n.º 6
0
int main( int argc, char* argv[] )
{
  int i;
  int results_data[DATA_SIZE];

  // Output the input arrays
  printArray( "input1",  DATA_SIZE, input_data1  );
  printArray( "input2",  DATA_SIZE, input_data2  );
  printArray( "verify", DATA_SIZE, verify_data );

#if PREALLOCATE
  for (i = 0; i < DATA_SIZE; i++)
  {
    results_data[i] = multiply( input_data1[i], input_data2[i] );
  }
#endif

  setStats(1);
  for (i = 0; i < DATA_SIZE; i++)
  {
    results_data[i] = multiply( input_data1[i], input_data2[i] );
  }
  setStats(0);
  
  // Print out the results
  printArray( "results", DATA_SIZE, results_data );

  // Check the results
  return verify( DATA_SIZE, results_data, verify_data );
}
Ejemplo n.º 7
0
int main( int argc, char* argv[] )
{
  // BENCHMARK #1:  Try to maximize store-data forwarding.

  setStats(1);
  benchmark();
  setStats(0);

  // Check the results
 

  return 0; // You'll probably want to call your verify function here.
            // Return 0 for success, anything else for failure.
}
Ejemplo n.º 8
0
    void Reset()
    {
        GC_Timer = 0;
        Heal_Timer = 0;
        Lesser_Healing_Timer = 0;
        Self_Lesser_Healing_Timer = 0;
        Flame_Shock_Timer = 20;
        Lightning_Bolt_Timer = 60;
        Earth_Shock_Timer = 150;
        Others_Heal_Timer = 0;
        Oom_timer = 0;

        Earth_Totem_Timer = 0;
        Fire_Totem_Timer = 0;
        Water_Totem_Timer = 0;
        Wind_Totem_Timer = 0;

        Potion_Timer = 0;
        Rez_Timer = 0;

        opponent = NULL;

        if (master) {
            setStats(CLASS_SHAMAN, m_creature->getRace(), master->getLevel());
       }
    }
Ejemplo n.º 9
0
void Player::shapeShift(Forms newForm)
{
	if (!(collisions[LEFT] || collisions[RIGHT]) && level > formLevels[newForm]) {
		form = newForm;
		setStats(newForm);
	}
}
Ejemplo n.º 10
0
        void Reset()
        {
            Pyroblast_cd = 0;
            FireBlast_cd = 0;
            DragonBreath_cd = 0;
            Living_Bomb_cd = 0;
            //Combustion_cd = 30000;//30 sec for reset
            Ice_Barrier_cd = 0;
            Iceblock_cd = 0;
            ConeofCold_cd = 0;
            Blizzard_cd = 10000;
            CounterSpell_cd = 0;
            Evocation_cd = 0;
            Blink_cd = 0;
            Bolt_cd = 0;
            Nova_cd = 0;
            polyCheckTimer = 0;
            fmCheckTimer = 0;
            Polymorph = false;
            clearcast = false;
            BOLT = 0;
            NOVA = 0;

            if (master)
            {
                setStats(CLASS_MAGE, me->getRace(), master->getLevel(), true);
                ApplyClassPassives();
                ApplyPassives(CLASS_MAGE);
           }
        }
Ejemplo n.º 11
0
    void Reset()
    {
        FirstAid_cd = 0;
        GC_Timer = 0;

        conflagarate_cd = 0;
        chaos_bolt_cd = 0;

        oom_spam = false;

        uint8 state = 1;
        next_state = 0;
        next_state_timer = 0;

        opponent = NULL;

        if (master) {
            setStats(CLASS_WARLOCK, m_creature->getRace(), master->getLevel());

            if (!m_creature->HasAura(56235)) m_creature->AddAura(56235,m_creature);  // Glyph of Conflagrate
            if (!m_creature->HasAura(63302)) m_creature->AddAura(63302,m_creature);  // Glyph of Haunt
            if (!m_creature->HasAura(17834)) m_creature->AddAura(17834,m_creature);  // Improved Immolation
            if (!m_creature->HasAura(17814)) m_creature->AddAura(17814,m_creature);  // Improved Corruption

        }
    }
Ejemplo n.º 12
0
//Used for copy constructor
void Projectile::operator=(Projectile& in)
{
	spriteData = in.spriteData;
	textureManager = in.textureManager;
	setStats(in.damage, in.minRange, in.maxRange, in.muzzelVelocity);

}
Ejemplo n.º 13
0
Player::Player(int x, int y) : Character(x, y, "data/spritesheets/player_human.png")
{
	setStats(HUMAN);
	health = MaxHealth;
	maxStamina = stamina;

	printf("Player created at pos (%i,%i) with stats: \n\tHP: %i\n\tStamina: %i\n", x, y, MaxHealth, stamina);
}
Ejemplo n.º 14
0
int main( int argc, char* argv[] )
{
  int i;
  int results_data[DATA_SIZE];

  // Output the input arrays

#if HOST_DEBUG
  printArray( "input1",  DATA_SIZE, input_data1  );
  printArray( "input2",  DATA_SIZE, input_data2  );
  printArray( "verify", DATA_SIZE, verify_data );
#endif

#if ( !HOST_DEBUG && PREALLOCATE )
  for (i = 0; i < DATA_SIZE; i++)
  {
    results_data[i] = multiply( input_data1[i], input_data2[i] );
  }
#endif

#if HOST_DEBUG
  for (i = 0; i < DATA_SIZE; i++)
  {
    results_data[i] = multiply( input_data1[i], input_data2[i] );
  }
#else
  setStats(1);
  for (i = 0; i < DATA_SIZE; i++)
  {
    results_data[i] = multiply( input_data1[i], input_data2[i] );
  }
  setStats(0);
#endif
  
  // Print out the results

#if HOST_DEBUG
  printArray( "results", DATA_SIZE, results_data );
#endif

  // Check the results

  finishTest(verify( DATA_SIZE, results_data, verify_data ));

}
Ejemplo n.º 15
0
int main(int argc, char *argv[]){
	// Instantiate tallied variables.
	int i, j, width, correct, hintCount, guesses, hints, trys;
	//Instantiate stats array
	int stats[5]  = {0};
	// Determine width of Matrix	
	width = getWidth(argc, argv);	
	// Generate matrix based on width
	complex** M = generateMatrix(width);	
	// Create randmom number generator.
	srand((unsigned)time(NULL));	
	// Create string array for user input.
	char choice[20];	
	//	Create array of type complex to store generated row column pairs
	complex uniques[width * width];	
	// Continue game until user quits.
	while(strncmp(choice, "Q", 1) != 0){
		// Generate random row column pairs for ser to guess.
		i = rand()  % width;
		j = rand()  % width;
		// Reset counting variable for each generated row column pair.
		correct = hints = guesses = 0;
		while(correct == 0 && strncmp(choice, "Q", 1) != 0){
			// Count number of hints used for this specific pair.
			hintCount = 0;
			//	Display message to user displaying options.
			printf("M[0][0]=%p. M[i][j]=%p What's i and j?\n(Q to Quit or H or
			 HH or HHH for hints.): ", &M[0][0], &M[i][j]);
			// Recieve and store user input.
			fgets(choice, 20, stdin);
			// Check if user requested hints
			hintCount += checkHints(choice, M, width);
			// If user did not want hints check is valid answer submitted
			if (hintCount == 0){
				//	Check answer
				correct = checkAnswer(choice, i, j);
			}
			//	Accumulate hints used.
			hints += hintCount;
			// Accumulate guesses used.
			guesses++;
		}
		//	If user has not quit calculate stats for pair.
		if (strncmp(choice, "Q", 1) != 0){
			setUnique(uniques, i, j, stats);
			setStats(stats,hints, guesses);
		}
	}
	// Print stats of the game after user has quit.
	printStats(stats);
	// terminate matrix.
	killMatrix(M, width);
}
Ejemplo n.º 16
0
        void Reset()
        {
            Torment_cd = 0;

            if (master && m_creatureOwner)
            {
                setStats(master->getLevel(), PET_TYPE_VOIDWALKER, true);
                ApplyPassives(PET_TYPE_VOIDWALKER);
                ApplyClassPassives();
                SetBaseArmor(162 * master->getLevel());
            }
        }
Ejemplo n.º 17
0
Mob::Mob(){
	setSTR(1);
	setDEF(1);
	setSPD(9);
	setINT(1);
	setStats(1,10,10,10, 1, 1);
	Name = "The Game";
	TTF_Font* font;
	SDL_Color fgColor = {255,255,0};
	SDL_Color fgColor1 = {255,255,255};
	font = TTF_OpenFont("../Fonts/coolvetica.ttf",25);
	enemyText[0] = TTF_RenderText_Blended(font,getName().c_str(),fgColor);
	enemyText[1] = TTF_RenderText_Blended(font,getName().c_str(),fgColor1);
}
Ejemplo n.º 18
0
        void Reset()
        {
            ArcaneShot_cd = 0;
            ChimeraShot_Timer = 0;
            SilencingShot_Timer = 0;
            AimedShot_Timer = 0;
            Feign_Death_Timer = 0;

            if (master)
            {
                setStats(CLASS_HUNTER, me->getRace(), master->getLevel(), true);
                ApplyClassPassives();
                ApplyPassives(CLASS_HUNTER);
             }
        }
Ejemplo n.º 19
0
int main( int argc, char* argv[] )
{
  // Output the input array
  printArray( "input", DATA_SIZE, input_data );
  printArray( "verify", DATA_SIZE, verify_data );

#if PREALLOCATE
  // If needed we preallocate everything in the caches
  sort(DATA_SIZE, verify_data);
  if (verify(DATA_SIZE, input_data, input_data))
    return 1;
#endif

  // Do the sort
  setStats(1);
  sort( DATA_SIZE, input_data );
  setStats(0);

  // Print out the results
  printArray( "test", DATA_SIZE, input_data );

  // Check the results
  return verify( DATA_SIZE, input_data, verify_data );
}
Ejemplo n.º 20
0
        void Reset()
        {
            Rain_of_fire_cd = 0;
            Haunt_cd = 0;
            conflagarate_cd = 0;
            chaos_bolt_cd = 0;
            fear_cd = 0;

            if (master)
            {
                setStats(CLASS_WARLOCK, me->getRace(), master->getLevel(), true);
                //TODO: passives
                ApplyClassPassives();
                ApplyPassives(CLASS_WARLOCK);
            }
        }
Ejemplo n.º 21
0
		void Reset()
		{
			oom_spam = false;

			opponent = NULL;

			GC_Timer = 0;
			ArcaneShot_cd = 0;
			ChimeraShot_Timer = 0;
			SilencingShot_Timer = 0;
			AimedShot_Timer = 0;

			if (master) {
				setStats(CLASS_HUNTER, m_creature->getRace(), master->getLevel());
		   }
		}
Ejemplo n.º 22
0
NPC::NPC(){
	TTF_Init();
	setStats(1,1,1,1,1,1); // need to be set because it is still a Entity instance, but these will never be used.
	npcName = "Default NPC";
	//speach.push_back("I don't know what kinda of NPC I am. :(");
	// default image - should not be used. But just in case
	SDL_Surface* tempSurface = IMG_Load("../Images/npc/default.png");
	npcImage = SDL_DisplayFormatAlpha(tempSurface);
	SDL_FreeSurface(tempSurface);
	setFrames(); // set the frames
	sentence = 0;
	font = TTF_OpenFont("../Fonts/coolvetica.ttf",25);
	fgColor.r = 255; 
	fgColor.g = 255;
	fgColor.b = 255;
	surfaceTextRect.x = 40;
	surfaceTextRect.y = 500;
}
Ejemplo n.º 23
0
int main()
{
   Creature test(10, 10, 100, 5, 0, "Joe");

   cout << "Displaying original stats..." << endl;
   displayStats(test);

   cin.get();

   cout << "Changing the stats..." << endl;
   changeStats(test);

   cin.get();

   cout << "Displaying stats after the changes..." << endl;
   displayStats(test);

   cin.get();

   cout << "Setting the stats..." << endl;
   setStats(test);

   cin.get();

   cout << "Displaying stats after setting..." << endl;
   displayStats(test);

   cin.get();

   cout << "Using an item..." << endl;
   useAnItem(test);

   cin.get();

   cout << "Displaying stats after using the item..." << endl;
   displayStats(test);

   cin.get();


   return 0;
}
Ejemplo n.º 24
0
// restore from serialization
void
ReadGroupStatsSet::
load(
    const char* filename)
{
    clear();

    assert(NULL != filename);
    std::ifstream ifs(filename);
    boost::archive::xml_iarchive ia(ifs);

    int numGroups;
    ia >> boost::serialization::make_nvp("numGroups", numGroups);
    ReadGroupStatsExporter se;
    for (int i=0; i<numGroups; i++)
    {
        ia >> boost::serialization::make_nvp("bogus", se);

        setStats(KeyType(se.bamFile.c_str(), se.readGroup.c_str()), se.groupStats);
    }
}
Ejemplo n.º 25
0
void plotBuilding(TCanvas *canvas, TH1F **s, TH1F **r, int n,TText* te,
		  char * option, double startingY, double startingX = .1,bool fit = false, unsigned int logx=0){
  canvas->Divide(2,(n+1)/2); //this should work also for odd n
  for(int i=0; i<n;i++){
    s[i]->SetMarkerStyle(20);
    r[i]->SetMarkerStyle(21);
    s[i]->SetMarkerColor(2);
    r[i]->SetMarkerColor(4);
    s[i]->SetMarkerSize(0.7);
    r[i]->SetMarkerSize(0.7);
    s[i]->SetLineColor(1);
    r[i]->SetLineColor(1);
    s[i]->SetLineWidth(1);
    r[i]->SetLineWidth(1);

    TPad *pad=canvas->cd(i+1);
    setStats(s[i],r[i], -1, 0, fit);
    if((logx>>i)&1)pad->SetLogx();
    r[i]->Draw("lpf");
    s[i]->Draw("sameslpf");
  }

}
runLooper_Iso_MinBias_skimdata_mc_EE( TString iso_opt = "ran_ecalIso03_mu",  bool jura = false){
 

  gStyle->SetHistLineWidth(2);
  gROOT->ForceStyle();
  gStyle->SetOptStat(1101);

  //gStyle->SetOptStat(0);
  TCanvas *c1 = new TCanvas("can", "can");
 
  TChain *MinBias_skimdata = new TChain("Events");
 
  TChain *MinBias_mc = new TChain("Events");
 
 
  MinBias_mc->Add("MC_skimmednTuple.root"); 
  MinBias_skimdata->Add("Data_skimmednTuple.root"); 
  
  
  TString draw_opt_MinBias_skimdata = Form("%s%s", iso_opt.Data(), ">> _MinBias_skimdata(50,-0.5,2)");
  MinBias_skimdata->Draw(draw_opt_MinBias_skimdata.Data(), baseline_minbias_EE);
  TH1F* Iso03_MinBias_skimdata = (TH1F*)gDirectory->Get("_MinBias_skimdata");
 
  Iso03_MinBias_skimdata->Sumw2();
  Iso03_MinBias_skimdata->SetTitle(iso_opt.Data());
  Iso03_MinBias_skimdata->SetLineColor(kRed);
  //Iso03_MinBias_skimdata->SetAxisRange(0,1,"Y");
  Iso03_MinBias_skimdata->GetYaxis()->SetRangeUser(-1000, 1);
//   Iso03_MinBias_skimdata->SetLineWidth(3);
 
  int   nbins_MinBias_skimdata = Iso03_MinBias_skimdata->GetNbinsX();
  float n_MinBias_skimdata = Iso03_MinBias_skimdata->Integral(0, nbins_MinBias_skimdata+1);
  Iso03_MinBias_skimdata->Scale(1/n_MinBias_skimdata);
  Iso03_MinBias_skimdata->SetBinContent(nbins_MinBias_skimdata, Iso03_MinBias_skimdata->Integral(nbins_MinBias_skimdata, nbins_MinBias_skimdata+1));
  Iso03_MinBias_skimdata->SetBinContent(nbins_MinBias_skimdata+1, 0);

  TString draw_opt_MinBias_mc = Form("%s%s", iso_opt.Data(), ">> _MinBias_mc(50,-0.5,2)");
  MinBias_mc->Draw(draw_opt_MinBias_mc.Data(), baseline_minbias_EE, "sames");
  TH1F* Iso03_MinBias_mc = (TH1F*)gDirectory->Get("_MinBias_mc");
  Iso03_MinBias_mc->SetTitle(iso_opt.Data());
  Iso03_MinBias_mc->GetYaxis()->SetTitle("Events");
  Iso03_MinBias_mc->GetXaxis()->SetTitle("GeV");
  
  //Iso03_MinBias_mc->SetAxisRange(0,1,"Y");
  // Iso03_MinBias_mc->GetYaxis()->SetRangeUser(-0.5, 1.2);
  //  Iso03_MinBias_mc->Sumw2();
  Iso03_MinBias_mc->SetLineColor(kBlue);
//   Iso03_MinBias_mc->SetLineWidth(3);

  int   nbins_MinBias_mc = Iso03_MinBias_mc->GetNbinsX();
  float n_MinBias_mc = Iso03_MinBias_mc->Integral(0, nbins_MinBias_mc+1);
  Iso03_MinBias_mc->Scale(1/n_MinBias_mc);
  Iso03_MinBias_mc->SetBinContent(nbins_MinBias_mc, Iso03_MinBias_mc->Integral(nbins_MinBias_mc, nbins_MinBias_mc+1));
  Iso03_MinBias_mc->SetBinContent(nbins_MinBias_mc+1,0);
  c1->SetLogy(1);
  Double_t xl1=.2, yl1=0.75, xl2=xl1+.2, yl2=yl1+.125;
  TLegend *leg = new TLegend(xl1,yl1,xl2,yl2);
  leg->SetFillStyle(0);
  leg->SetFillColor(0);
  //leg->SetLineStyle(0);
  // leg->SetLineColor(0);
 
  leg->AddEntry(Iso03_MinBias_skimdata,"MinBias_skimdata","L");   // h1 and h2 are histogram pointers
  leg->AddEntry(Iso03_MinBias_mc,"MinBias_mc","L");   // h1 and h2 are histogram pointers

  

  setStats(Iso03_MinBias_skimdata,Iso03_MinBias_mc, 0.6, 0.6, false);
 
 
  Iso03_MinBias_mc->Draw("");
  Iso03_MinBias_skimdata->Draw("esames");

 
 
 //  // leg->Draw();
     if(!jura){
     TString save_opt1= Form("%s%s%s","../../../plots/900GeV/EE_skimdata_mc_",iso_opt.Data(),"_V03-00-23.gif" );
     TString save_opt2= Form("%s%s%s","../../../plots/900GeV/EE_skimdata_mc_",iso_opt.Data(),"_V03-00-23.eps" );
   }
   else{
     TString save_opt1= Form("%s%s%s","../../../plots/900GeV/EE_skimdata_mc_",iso_opt.Data(),"_V03-00-23_jura.gif" );
     TString save_opt2= Form("%s%s%s","../../../plots/900GeV/EE_skimdata_mc_",iso_opt.Data(),"__V03-00-23_jura.eps" ); 

   }
   c1->SaveAs(save_opt1.Data());
   c1->SaveAs(save_opt2.Data());
   //   can->Range(-0.5142857,-3.438716,6.0,0.4504303);
 
 
}
Ejemplo n.º 27
0
void decovspeak(bool printgif = false){
  
  //load macros and set style-----------------------------------
  //gROOT->SetStyle("Plain");
  gStyle->SetPalette(1);
  //gStyle->SetOptStat("mr");
  //gStyle->SetOptFit(0);
  gStyle->SetOptStat(0);
  int colors[]       = {2,4,6,1};
  
  vector<char*> filenames;
  vector<char*> types;
  
  //filenames.push_back("crabjobs/lpc/Commissioning10-GOODCOLL-v8_mintrkmom1_ALLPEAK/merged.root");
  filenames.push_back("Commissioning10-GOODCOLL-v8_ALLPEAK/res/merged.root");
  
  //filenames.push_back("crabjobs/lpc/Commissioning10-GOODCOLL-v8_mintrkmom1_copy/merged_2.root");
  filenames.push_back("Commissioning10-GOODCOLL-v8/res/merged.root");

  filenames.push_back("Commissioning10-GOODCOLL-v8_toblatebp06/res/merged.root");

  //filenames.push_back("crabjobs/lpc/Spring10-START3X_V26A_356ReReco-v1_standard_geom_mintrkmom1/merged.root");
  //filenames.push_back("crabjobs/lpc/Spring10-START3X_V26A_356ReReco-v1_standard_geom/merged.root");
  
  
  //filenames.push_back("crabjobs/lpc/Commissioning10-GOODCOLL-v8_mintrkmom1_toblatebp063/merged.root");
  
  types.push_back("peak");
  types.push_back("deco");
  types.push_back("deco (BP)");  
  //types.push_back("MC");

  vector<char*> cantitles;
  //cantitles.push_back("du_dw");         //delta u, delta w TH1s
  //cantitles.push_back("duvsdth");         //v+/- split du vs. delta tan(theta) TH2s
  cantitles.push_back("duvsdth_tgraph");  //v+/- split du vs. delta tan(theta) TGraphs
  
  const unsigned int nfiles = filenames.size();
  const unsigned int ncan   = cantitles.size();

  TCanvas *can[ncan];
  int idx = 0;

  TFile* file[nfiles];
  TH1F*  hdu[nfiles];
  TH1F*  hdw[nfiles];
  TH2F*  hduthp[nfiles];
  TH2F*  hduthm[nfiles];

  //dummy legend
  TLegend *leg1=new TLegend(0.15,0.65,0.35,0.85);

  for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++){
    TH1F *hdummy=new TH1F(Form("hdummy_%i",ifile),Form("hdummy_%i",ifile),1,0,1);
    hdummy->SetLineColor(colors[ifile]);
    hdummy->SetMarkerColor(colors[ifile]);
    leg1->AddEntry(hdummy,types.at(ifile));
  }

  leg1->SetBorderSize(1);
  leg1->SetFillColor(0);

  cout << "Getting histos... " << endl;

  for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++){
    cout << "Opening " << types.at(ifile) <<" file   " << filenames.at(ifile) << endl;
    file[ifile] = TFile::Open(filenames.at(ifile));
   
    hdu[ifile] = (TH1F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/du_TOB");
    hdw[ifile] = (TH1F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/dw_TOB");

    //    if( strcmp(types.at(ifile),"deco (BP)") == 0){
      hduthp[ifile] = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vp_wp_all");
      TH2F* htp     = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vp_wm_all");
      hduthp[ifile] -> Add(htp);
      
      hduthm[ifile] = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vm_wp_all");
      TH2F* htm     = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vm_wm_all");
      hduthm[ifile] -> Add(htm);
    // }else{
//       hduthp[ifile] = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vp_wp");
//       TH2F* htp     = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vp_wm");
//       hduthp[ifile] -> Add(htp);
      
//       hduthm[ifile] = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vm_wp");
//       TH2F* htm     = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vm_wm");
//       hduthm[ifile] -> Add(htm);
//     }
     
  }
  
  
  cout << "Plotting du, dw histos... " << endl;

  if(draw("du_dw",cantitles)){
    can[idx]=new TCanvas(Form("can_%i",idx),"du_dw",1200,450);
    can[idx]->Divide(2,1);
    
    can[idx]->cd(1);
    setStats(hdu,nfiles,0.8,0.65,0.15);  
    plotHists(hdu,nfiles,"TOB","#Delta u [#mum]",-1000,1000,10,3);
    leg1->Draw();

    can[idx]->cd(2);
    setStats(hdw,nfiles,0.8,0.65,0.15);  
    plotHists(hdw,nfiles,"TOB","#Delta w [#mum]",-2000,2000,10,3);
    leg1->Draw();
        
    idx++;
  }

  cout << "Plotting du vs. delta tan theta split v+/v- histos... " << endl;
  //du vs. delta tan(theta) split v+/v- histos-----------------------------------
  
  if(draw("duvsdth",cantitles)){

    can[idx]=new TCanvas(cantitles.at(idx),cantitles.at(idx),1200,450*nfiles);
    can[idx]->Divide(nfiles,2);
    
    for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++){
      
      plotHistsTH2(hduthp[ifile],
                   Form("TOB %s (v+)",types.at(ifile)),"tan(#theta_{trk})-tan(#theta_{L})",
                   "#Delta u [#mum]",-0.3,0.5,-200,200,0,can[idx],ifile+1);
      
      plotHistsTH2(hduthm[ifile],
                   Form("TOB %s (v-)",types.at(ifile)),"tan(#theta_{trk})-tan(#theta_{L})",
                        "#Delta u [#mum]",-0.5,0.3,-200,200,0,can[idx],ifile+2);
      
    }
  
    idx++;

  }
  
  cout << "Plotting du vs. delta tan theta split v+/v- TGraphs... " << endl;

  //TGraph du vs. delta tan theta----------------------------------------

  if(draw("duvsdth_tgraph",cantitles)){

    can[idx]=new TCanvas(cantitles.at(idx),cantitles.at(idx),1200,450);
    can[idx]->Divide(2,1);
  
    vector<float> thetabinsp;
    vector<float> thetabinsm;

    for(int ibin=14;ibin<=30;ibin++)    thetabinsp.push_back(ibin*0.05-1);
    for(int ibin=10;ibin<=26;ibin++)    thetabinsm.push_back(ibin*0.05-1);
    
 
    TLatex *t=new TLatex();
    t->SetNDC();

    TGraphErrors *gduthp[nfiles];
    TGraphErrors *gduthm[nfiles];
        
    TF1 *fduthp[nfiles];
    TF1 *fduthm[nfiles];

    stringstream sduthp1[nfiles];
    stringstream sduthp2[nfiles];
    
    stringstream sduthm1[nfiles];
    stringstream sduthm2[nfiles];

    bool fit = false;

    for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++ ){
      
      can[idx]->cd(1);
      
      hduthp[ifile]->SetName(Form("hduthp_%i",ifile));
      gduthp[ifile] = getTGraphFromTH2(hduthp[ifile],thetabinsp, colors[ifile] );
      fduthp[ifile]=new TF1(Form("fduthp_%i",ifile),"pol1");
      fduthp[ifile]->SetLineColor( colors[ifile] );
      if(fit) gduthp[ifile]->Fit(fduthp[ifile]);
         
      float dwp     = fduthp[ifile]->GetParameter(1);
      float dwperr  = fduthp[ifile]->GetParError(1);
      float bp      = fduthp[ifile]->GetParameter(0);
      float bperr   = fduthp[ifile]->GetParError(0);
      
      sduthp1[ifile] << "#DeltaW = " << fround(dwp,3) << " #pm " << fround(dwperr,3) << " #mum" << endl;
      sduthp2[ifile] << "#Deltatan(LA) = " << fround(bp/(235.-dwp),3) << " #pm " << fround(bperr/(235.-dwp),4) << endl;
      
      gduthp[ifile]->SetTitle("TOB (V+)");
      gduthp[ifile]->Draw("AP");
      gduthp[ifile]->GetXaxis()->SetTitle("<tan(#theta_{trk})-tan(#theta_{LA})>");
      gduthp[ifile]->GetYaxis()->SetTitle("<#Deltau> [#mum]");
      //gduthp[ifile]->GetXaxis()->SetLimits(-1,1);
      gduthp[ifile]->SetMinimum(-25);
      gduthp[ifile]->SetMaximum(20);
      
      
      can[idx]->cd(2);

      hduthm[ifile]->SetName(Form("hduthm_%i",ifile));
      gduthm[ifile] = getTGraphFromTH2(hduthm[ifile],thetabinsm, colors[ifile] );
      fduthm[ifile]=new TF1(Form("fduthm_%i",ifile),"pol1");
      fduthm[ifile]->SetLineColor( colors[ifile] );
      if(fit) gduthm[ifile]->Fit(fduthm[ifile]);
         
      float dwm     = fduthm[ifile]->GetParameter(1);
      float dwmerr  = fduthm[ifile]->GetParError(1);
      float bm      = fduthm[ifile]->GetParameter(0);
      float bmerr   = fduthm[ifile]->GetParError(0);
      
      sduthm1[ifile] << "#DeltaW = " << fround(dwm,3) << " #pm " << fround(dwmerr,3) << " #mum" << endl;
      sduthm2[ifile] << "#Deltatan(LA) = " << fround(bm/(235.-dwm),3) << " #pm " << fround(bmerr/(235.-dwm),4) << endl;
      
      gduthm[ifile]->SetTitle("TOB (V-)");
      gduthm[ifile]->Draw("AP");
      gduthm[ifile]->GetXaxis()->SetTitle("<tan(#theta_{trk})-tan(#theta_{LA})>");
      gduthm[ifile]->GetYaxis()->SetTitle("<#Deltau> [#mum]");
      //gduthm[ifile]->GetXaxis()->SetLimits(-1,1);
      gduthm[ifile]->SetMinimum(-35);
      gduthm[ifile]->SetMaximum(15);
      //leg1->Draw();
    }

    for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++ ){
      
      can[idx]->cd(1);
      if( ifile == 0 )   gduthp[ifile]->Draw("AP");
      else               gduthp[ifile]->Draw("sameP");
      leg1->Draw();

      if(fit){
        t->SetTextColor( colors[ifile] );
        t->DrawLatex(0.15,0.85-ifile*0.2,sduthp1[ifile].str().c_str());
        t->DrawLatex(0.15,0.75-ifile*0.2,sduthp2[ifile].str().c_str());
      }

      can[idx]->cd(2);
      if( ifile == 0 )   gduthm[ifile]->Draw("AP");
      else               gduthm[ifile]->Draw("sameP");
      //leg1->Draw();

      if(fit){
        t->SetTextColor( colors[ifile] );
        t->DrawLatex(0.15,0.85-ifile*0.2,sduthm1[ifile].str().c_str());
        t->DrawLatex(0.15,0.75-ifile*0.2,sduthm2[ifile].str().c_str());
      }
    
    }

    TCanvas *cdiff = new TCanvas("cdiff","duvsdtanthetadiff",1200,450);
    cdiff->Divide(2,1);

    cdiff->cd(1);


    TGraphErrors *gduthpdiff = diffTGraph(gduthp[0],gduthp[2],"TOB DECO - PEAK (v+)","<tan(#theta_{trk})-tan(#theta_{LA})>","<#Deltau> (#mum)");
    TF1* fduthpdiff=new TF1("fduthpdiff","pol1",-0.5,0.5);
    gduthpdiff->Fit(fduthpdiff,"R");
    gduthpdiff->Draw("AP");
    
    float dwp     = fduthpdiff->GetParameter(1);
    float dwperr  = fduthpdiff->GetParError(1);
    float bp      = fduthpdiff->GetParameter(0);
    float bperr   = fduthpdiff->GetParError(0);
    
    stringstream sduthp1diff;
    stringstream sduthp2diff;

    sduthp1diff << "#DeltaW = " << fround(dwp,3) << " #pm " << fround(dwperr,3) << " #mum" << endl;
    sduthp2diff << "#Deltatan(LA) = " << fround(bp/(235.-dwp),3) << " #pm " << fround(bperr/(235.-dwp),4) << endl;
  
    t->SetTextColor(1);
    t->DrawLatex(0.15,0.85,sduthp1diff.str().c_str());
    t->DrawLatex(0.15,0.75,sduthp2diff.str().c_str());

    cdiff->cd(2);

    TGraphErrors *gduthmdiff = diffTGraph(gduthm[0],gduthm[2],"TOB DECO - PEAK (v-)","<tan(#theta_{trk})-tan(#theta_{LA})>","<#Deltau> (#mum)");
    TF1* fduthmdiff=new TF1("fduthmdiff","pol1",-0.5,0.5);
    gduthmdiff->Fit(fduthmdiff,"R");
    gduthmdiff->Draw("AP");

    float dwm     = fduthmdiff->GetParameter(1);
    float dwmerr  = fduthmdiff->GetParError(1);
    float bm      = fduthmdiff->GetParameter(0);
    float bmerr   = fduthmdiff->GetParError(0);
    
    stringstream sduthm1diff;
    stringstream sduthm2diff;

    sduthm1diff << "#DeltaW = " << fround(dwm,3) << " #pm " << fround(dwmerr,3) << " #mum" << endl;
    sduthm2diff << "#Deltatan(LA) = " << fround(bm/(235.-dwm),3) << " #pm " << fround(bmerr/(235.-dwm),4) << endl;
  
    t->DrawLatex(0.15,0.85,sduthm1diff.str().c_str());
    t->DrawLatex(0.15,0.75,sduthm2diff.str().c_str());

    idx++;
  }

  
  for(unsigned int ican=0;ican<ncan;ican++){
    can[ican]->Modified();
    can[ican]->Update();
    if(printgif) can[ican]->Print(Form("plots/%s_TOB.gif",cantitles.at(ican)));
  }
  
}
Ejemplo n.º 28
0
void plot6histos(TCanvas *canvas, 
		TH1F *s1,TH1F *r1, TH1F *s2,TH1F *r2, 
		TH1F *s3,TH1F *r3, TH1F *s4,TH1F *r4,
		TH1F *s5,TH1F *r5, TH1F *s6,TH1F *r6,
		TText* te,
	       char * option, double startingY, double startingX = .1,bool fit = false){
  canvas->Divide(2,3);

  s1->SetMarkerStyle(20);
  r1->SetMarkerStyle(21);
  s1->SetMarkerColor(2);
  r1->SetMarkerColor(4);
  s1->SetMarkerSize(0.7);
  r1->SetMarkerSize(0.7);
  s1->SetLineColor(2);
  r1->SetLineColor(4);
  s1->SetLineWidth(2);
  r1->SetLineWidth(2);

  s2->SetMarkerStyle(20);
  r2->SetMarkerStyle(21);
  s2->SetMarkerColor(2);
  r2->SetMarkerColor(4);
  s2->SetMarkerSize(0.1);
  r2->SetMarkerSize(0.1);
  s2->SetLineColor(2);
  r2->SetLineColor(4);
  s2->SetLineWidth(2);
  r2->SetLineWidth(2);

  s3->SetMarkerStyle(20);
  r3->SetMarkerStyle(21);
  s3->SetMarkerColor(2);
  r3->SetMarkerColor(4);
  s3->SetMarkerSize(0.7);
  r3->SetMarkerSize(0.7);
  s3->SetLineColor(2);
  r3->SetLineColor(4);
  r3->SetLineWidth(2);
  s3->SetLineWidth(2);

  s4->SetMarkerStyle(20);
  r4->SetMarkerStyle(21);
  s4->SetMarkerColor(2);
  r4->SetMarkerColor(4);
  s4->SetMarkerSize(0.7);
  r4->SetMarkerSize(0.7);
  s4->SetLineColor(2);
  r4->SetLineColor(4);
  r4->SetLineWidth(2);
  s4->SetLineWidth(2);

  s5->SetMarkerStyle(20);
  r5->SetMarkerStyle(21);
  s5->SetMarkerColor(2);
  r5->SetMarkerColor(4);
  s5->SetMarkerSize(0.7);
  r5->SetMarkerSize(0.7);
  s5->SetLineColor(2);
  r5->SetLineColor(4);
  r5->SetLineWidth(2);
  s5->SetLineWidth(2);

  s6->SetMarkerStyle(20);
  r6->SetMarkerStyle(21);
  s6->SetMarkerColor(2);
  r6->SetMarkerColor(4);
  s6->SetMarkerSize(0.7);
  r6->SetMarkerSize(0.7);
  s6->SetLineColor(2);
  r6->SetLineColor(4);
  r6->SetLineWidth(2);
  s6->SetLineWidth(2);


  //setStats(r1,s1, startingY, startingX, fit);
  canvas->cd(1);
  setStats(s1,r1, -1, 0, false);
  r1->Draw();
  s1->Draw("sames");
  double kstest = mykolmo(s1,r1);
  TPad *c1_1 = canvas->GetPad(1);
  if(kstest<0.7){
     c1_1->SetFillColor(kYellow);}

  canvas->cd(2);
  gPad->SetLogy(); 
  setStats(s2,r2, -1, 0, false);
  s2->Draw();
  r2->Draw("sames");
  double kstest = mykolmo(s2,r2);
  TPad *c2_2 = canvas->GetPad(2);
  if(kstest<0.7){
     c2_2->SetFillColor(kYellow);}

  canvas->cd(3);
  setStats(s3,r3, -1, 0, false);
  r3->Draw();
  s3->Draw("sames");
  double kstest = mykolmo(s3,r3);
  TPad *c3_3 = canvas->GetPad(3);
  if(kstest<0.7){
     c3_3->SetFillColor(kYellow);}

  canvas->cd(4);
  gPad->SetLogy(); 
  setStats(s4,r4, -1, 0, false);
  s4->Draw();
  r4->Draw("sames");
  double kstest = mykolmo(s4,r4);
  TPad *c4_4 = canvas->GetPad(4);
  if(kstest<0.7){
     c4_4->SetFillColor(kYellow);}

  canvas->cd(5);
  setStats(s5,r5, -1, 0, false);
  r5->Draw();
  s5->Draw("sames");
  double kstest = mykolmo(s5,r5);
  TPad *c5_5 = canvas->GetPad(5);
  if(kstest<0.7){
     c5_5->SetFillColor(kYellow);}

  canvas->cd(6);
  gPad->SetLogy(); 
  setStats(s6,r6, -1, 0, false);
  r6->Draw();
  s6->Draw("sames");
  double kstest = mykolmo(s6,r6);
  TPad *c6_6 = canvas->GetPad(6);
  if(kstest<0.7){
     c6_6->SetFillColor(kYellow);}

}
Ejemplo n.º 29
0
void createPlot(TCanvas *canvas, TFile *file,  TString type, TString name, int nbins, double xMin, double xMax, 
		bool runSecTrackColl, TText* te, char * option, 
		double startingY, double startingX = .1, bool fit = false, 
		bool logx=false, bool logy=false, int normScale){
  TTree *tree;
  tree = file->GetObject("tree",tree);
  
  TH1F *h1; // for ctf trk
  h1 = new TH1F("h1", "h1", nbins, xMin, xMax);
  
  TH1F *h2; // for sectrk
  h2 = new TH1F("h2","h2",nbins,xMin,xMax);
  

  TString x_title = "";
  TString y_title = "Number of Tracks";
  
  // If the histogram is filled per track, cut on the event and track.
  // Otherwise, cut only on the event level

  TString basecut="isTechBit40&&!isBeamHalo&&hasGoodPvtx==1";
  TString evtcutstring=basecut;
  evtcutstring.Append("&&EVTSELECTION");
  TString trkcutstring=basecut;
  trkcutstring.Append("&&EVTSELECTION&&TRKSELECTION");
  
  TString ctf_name =  "ctf_";

  if( name.Contains("cluster_charge",TString::kExact) ) {
    ctf_name = "ctf"; 
    ctf_name.Append(name);   
    tree->Project("h1",ctf_name, evtcutstring);
  }
  else if (name == "recz_pxlpvtx") {
    tree->Project("h1",name, evtcutstring);
    y_title="Number of pixelVertices";
  }
  else if( name.Contains("pvtx")) {
    if(name.Contains("rec")) {  
      y_title="Number of Real PrimaryVertices";
      TString vtxcutstring=evtcutstring;
      tree->Project("h1",name, vtxcutstring.Append("&&isFake_pvtx==0")); 
    }
    else { 
      y_title="Number of Valid PrimaryVertices";
      tree->Project("h1",name, evtcutstring);
    }
  }
  else if ( name.Contains("nVertices",TString::kExact) 
	    || name.Contains("hasGoodPvtx",TString::kExact) 
	    || name.Contains("nPixelVertices",TString::kExact) 
	    || name.Contains("isTechBit40",TString::kExact)
            || name.Contains("isBSC",TString::kExact)  
	    || name.Contains("bsX0",TString::kExact)
	    || name.Contains("bsY0",TString::kExact)
	    || name.Contains("bsZ0",TString::kExact)
	    || name.Contains("glob_ls",TString::kExact) 
	    || name.Contains("glob_bx",TString::kExact)   
	    ) {
    tree->Project("h1",name, evtcutstring);
    y_title = "Number of Events";
  }
  else { 
    ctf_name.Append(name); 
    tree->Project("h1",ctf_name, trkcutstring);
  }

  if(runSecTrackColl && ctf_name.Contains("ctf")) {
    TString sectrk_name = ctf_name.ReplaceAll("ctf","sectrk");
    tree->Project("h2",sectrk_name,evtcutstring);
  }
 
  // Set the x_tile 
  if( name == "n" ) {
    x_title = "Tracks per Event" ;
    y_title = "Number of Events";
  }
  
  // All Clusters
  if( name.Contains("clusterCharge_all",TString::kExact) )  x_title = "All Cluster Charge (ADC Counts)";
  if( name.Contains("clusterCharge_TIB",TString::kExact) )  x_title = "All TIB Cluster Charge (ADC Counts)";
  if( name.Contains("clusterCharge_TID",TString::kExact) )  x_title = "All TID Cluster Charge (ADC Counts)";
  if( name.Contains("clusterCharge_TOB",TString::kExact) )  x_title = "All TOB Cluster Charge (ADC Counts)";
  if( name.Contains("clusterCharge_TEC",TString::kExact) )  x_title = "All TEC Cluster Charge (ADC Counts)";
  
  if( name.Contains("cluster_charge",TString::kExact) ) {
    x_title = " Strip Cluster Charge (ADC Counts)";
    y_title = "Number of Clusters";
  }
    
  // Histograms filled per Track
  if( name.Contains("nHit",TString::kExact) )  x_title = "Number of Valid Hits per Track" ;
  if( name.Contains("nLostHit",TString::kExact) )  x_title = "Number of Lost Hits per Track" ; 
  if( name.Contains("nPXBhit",TString::kExact) )  x_title = "Number of PixelBarrel Hits per Track" ;
  if( name.Contains("nPXFhit",TString::kExact) )  x_title = "Number of PixEndcap Hits per Track" ;
  if( name.Contains("nTIBhit",TString::kExact) )  x_title = "Number of TIB Hits per Track" ;
  if( name.Contains("nTOBhit",TString::kExact) )  x_title = "Number of TOB Hits per Track" ;
  if( name.Contains("nTIDhit",TString::kExact) )  x_title = "Number of TID Hits per Track" ;
  if( name.Contains("nTEChit",TString::kExact) )  x_title = "Number of TEC Hits per Track" ;
  if( name.Contains("nPixelHit",TString::kExact) )  x_title = "Number of Pixel Hits per Track" ;

  if( name.Contains("nLayers",TString::kExact) )  x_title = "Number of Layers per Track" ;
  if( name.Contains("nLayers3D",TString::kExact) )  x_title = "Number of 3D Layers per Track" ;
  if( name.Contains("nPXBLayers",TString::kExact) )  x_title = "Number of PixelBarrel Layers per Track" ;
  if( name.Contains("nPXFLayers",TString::kExact) )  x_title = "Number of PixEndcap Layers per Track" ;
  if( name.Contains("eta",TString::kExact) )  x_title = "Track Pseudorapidity" ;
  if( name.Contains("pt",TString::kExact) )  x_title = "p_{T} (GeV)" ;
  if( name.Contains("phi",TString::kExact) )  x_title = "Track Azimuthal Angle" ;
  if( name.Contains("dxy",TString::kExact) )  x_title = "Track dxy wrt (0,0,0) (cm)" ;
  if( name.Contains("dz",TString::kExact) )  x_title = "Track dz wrt (0,0,0) (cm)" ;
  if( name.Contains("dxyCorr",TString::kExact) )  x_title = "Track dxy wrt BS (cm)" ;  
  if( name.Contains("dzCorr",TString::kExact) )  x_title = "Track dz wrt BS (cm)" ;  
  if( name.Contains("xPCA",TString::kExact) )  x_title = "x PCA (cm)" ;
  if( name.Contains("yPCA",TString::kExact) )  x_title = "y PCA (cm)" ;
  if( name.Contains("zPCA",TString::kExact) )  x_title = "z PCA (cm)" ;
  if( name.Contains("chi2ndof",TString::kExact) )  x_title = "#chi^{2}/ndf" ;
  if( name.Contains("algo",TString::kExact) )  x_title = "Track Algorithm" ;  
  if( name.Contains("isHighPurity",TString::kExact) )  x_title = "isHighPurity" ;  

  // pvtx stuff
  if( name.Contains("recx_pvtx",TString::kExact) )  x_title = "PrimaryVertex Position X (cm)" ;
  if( name.Contains("recy_pvtx",TString::kExact) )  x_title = "PrimaryVertex Position Y (cm)" ;
  if( name.Contains("recz_pvtx",TString::kExact) )  x_title = "PrimaryVertex Position Z (cm)" ;
  if( name.Contains("isFake_pvtx",TString::kExact) )  x_title = "PrimaryVertex isFake()";  
  if( name.Contains("recz_pxlpvtx",TString::kExact) )  x_title = "pixelVertices Position Z (cm)" ;

  //beamspot
  if( name.Contains("bsX0",TString::kExact) )  x_title = "BeamSpot X (cm)" ;
  if( name.Contains("bsY0",TString::kExact) )  x_title = "BeamSpot Y (cm)" ;
  if( name.Contains("bsZ0",TString::kExact) )  x_title = "BeamSpot Z (cm)" ;
   
  

  if( name.Contains("nVertices",TString::kExact) )  x_title = "Number of PrimaryVertices";
  if( name.Contains("nPixelVertices",TString::kExact) )  x_title = "Number of pixelVertices";
  if( name.Contains("hasGoodPvtx",TString::kExact) )  x_title = "hasRealPrimaryVertex";
  if( name.Contains("isTechBit40",TString::kExact) )  x_title = "isTechBit40";
  if( name.Contains("isBSC",TString::kExact) )  x_title = "isBSC";
  if( name.Contains("glob_bx",TString::kExact) )  x_title = "bunchCrossing";

  if(logx) gPad->SetLogx();
  else gPad->SetLogx(0);
  if(logy) gPad->SetLogy();
  else gPad->SetLogy(0);
  
  h1->SetLineColor(4);
  h1->SetMarkerColor(4);
  h1->SetMarkerStyle(21);
  h1->SetMarkerSize(0.7);
  //h1->StatOverflows(kTRUE); 
  h1->SetLineWidth(3);
  h1->GetXaxis()->SetTitle(x_title); 
  h1->GetYaxis()->SetTitle(y_title);
  
  if(runSecTrackColl) {
    h2->SetLineColor(2);
    h2->SetMarkerColor(2);
    h2->SetMarkerStyle(21);
    h2->SetMarkerSize(0.7);
    //h2->StatOverflows(kTRUE); 
    h2->SetLineWidth(3);
    h2->GetXaxis()->SetTitle(x_title);  
    h2->GetYaxis()->SetTitle(y_title);
  }

  if(runSecTrackColl) {
    setStats(h1, h2, startingY,  startingX, fit);
    NormalizeHistograms(h1, h2, normScale);
    if(!logy && !logx)
      fixRangeY(h1,h2); //I don't know why, with logy or logx, this line creates trouble in overlaying two plots if there is an Y=0 bin
  }
  else {
    setStats(h1, startingY,  startingX, fit);
  }
    
  gPad->Update();
  h1->Draw();
  if(runSecTrackColl) h2->Draw("sames");
  canvas->Update();
  
  TLegend *leg = new TLegend(0.16,0.93,0.60,0.995);
  leg->SetTextSize(0.04);
  leg->SetTextFont(42);
  leg->SetFillColor(0); 
  leg->SetBorderSize(1);
  leg->AddEntry(h1, "TRACKCOLLECTION" , "L" );
  if(runSecTrackColl)
    leg->AddEntry(h2, "SECTRKCOLLECTION" , "L" );
  if(type == "perTrk" || name == "n" || name == "cluster_charge") 
    leg->Draw("SAME");
  
  TString pngdirname ="PNGDIR/SAMPLE/CUTSTRING/";
  TString pngfilename = pngdirname;
  pngfilename.Append(name);
  pngfilename.Append(".png");
  canvas->Print(pngfilename);

  TString epsdirname ="EPSDIR/SAMPLE/CUTSTRING/";
  TString epsfilename = epsdirname;
  epsfilename.Append(name);
  epsfilename.Append(".eps");
  canvas->Print(epsfilename);

  delete h1;
  delete h2; 
}
Ejemplo n.º 30
0
void decovspeak(bool printgif = false){
  
  //load macros and set style-----------------------------------
  //gROOT->SetStyle("Plain");
  gStyle->SetPalette(1);
  //gStyle->SetOptStat("mr");
  //gStyle->SetOptFit(0);
  gStyle->SetOptStat(0);
  int colors[]       = {2,4,6,1};
  
  vector<string> filenames;
  vector<string> types;
  
  //-----------files---------------------
  filenames.push_back("MinimumBias_Commissioning10-GOODCOLL-Jun9thSkim_v1_PEAK/res/merged.root");
  filenames.push_back("MinimumBias_Commissioning10-GOODCOLL-Jun9thSkim_v1/res/merged.root");

  types.push_back("peak");
  types.push_back("deco");

  //----------variables to plot----------
  vector<string> cantitles;
  //cantitles.push_back("du_dw");           //delta u, delta w TH1s
  //cantitles.push_back("duvsdth");         //v+/- split du vs. delta tan(theta) TH2s
  //cantitles.push_back("duvsdth_tgraph");  //v+/- split du vs. delta tan(theta) TGraphs
  cantitles.push_back("dwTEC");             //dw histos for TEC ring/zside/petalside
  
  const unsigned int nfiles = filenames.size();
  const unsigned int ncan   = cantitles.size();

  TCanvas *can[ncan];
  int idx = 0;

  TFile* file[nfiles];
  TH1F*  hdu[nfiles];
  TH1F*  hdw[nfiles];
  TH2F*  hduthp[nfiles];
  TH2F*  hduthm[nfiles];

  //dummy legend
  TLegend *leg1=new TLegend(0.15,0.65,0.35,0.85);

  for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++){
    TH1F *hdummy=new TH1F(Form("hdummy_%i",ifile),Form("hdummy_%i",ifile),1,0,1);
    hdummy->SetLineColor(colors[ifile]);
    hdummy->SetMarkerColor(colors[ifile]);
    leg1->AddEntry(hdummy,types.at(ifile).c_str());
  }

  leg1->SetBorderSize(1);
  leg1->SetFillColor(0);

  cout << "Getting histos... " << endl;

  for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++){
    cout << "Opening " << types.at(ifile) <<" file   " << filenames.at(ifile) << endl;
    file[ifile] = TFile::Open(filenames.at(ifile).c_str());
   
    hdu[ifile] = (TH1F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/du_TOB");
    hdw[ifile] = (TH1F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/dw_TOB");


      hduthp[ifile] = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vp_wp_all");
      TH2F* htp     = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vp_wm_all");
      hduthp[ifile] -> Add(htp);
      
      hduthm[ifile] = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vm_wp_all");
      TH2F* htm     = (TH2F*) file[ifile]->Get("PeakDecoResiduals/PeakDecoResiduals/duvsdtantheta_TOB_vm_wm_all");
      hduthm[ifile] -> Add(htm);
  }
  
  
  if(draw("dwTEC",cantitles)){

    gStyle->SetOptFit(11);

    TCanvas *dwTEC_can[2][2];
    TH1F* hdw_TEC[2][2][7];
    TF1* hgaus[2][2][7];
    float mean[2][2][7];
    float meanerr[2][2][7];
    TGraphErrors *g00[nfiles];
    TGraphErrors *g01[nfiles];
    TGraphErrors *g10[nfiles];
    TGraphErrors *g11[nfiles];

    can[idx]=new TCanvas(Form("can_%i",idx),"du_dw",1200,600);
    can[idx]->Divide(2,1);

    for( unsigned int ifile = 0 ; ifile < nfiles ; ++ifile ){

      for(int i=0;i<2;i++){
        for(int j=0;j<2;j++){

          dwTEC_can[i][j] = new TCanvas(Form("dwTEC_can_%i_%i",i,j),"",1200,600);
          dwTEC_can[i][j]->Divide(4,2);
        
          for(int k=0;k<7;k++){

            dwTEC_can[i][j]->cd(k+1);
  
            hdw_TEC[i][j][k] = (TH1F*) file[ifile]->Get(Form("PeakDecoResiduals/PeakDecoResiduals/hdw_TEC_z%i_petal%i_ring%i_stereo0",i,j,k+1));
        
            hgaus[i][j][k] = new TF1(Form("gaus_%i_%i_%i",i,j,k+1),"gaus",-500,500);
            hgaus[i][j][k] -> SetLineColor(2);
            hdw_TEC[i][j][k]->Rebin(10);
            hdw_TEC[i][j][k]->Draw();
            hdw_TEC[i][j][k]->Fit( hgaus[i][j][k] ,"R");
            mean[i][j][k]    = hgaus[i][j][k]->GetParameter(1);
            meanerr[i][j][k] = hgaus[i][j][k]->GetParError(1);
          }
        }
      }

      can[idx]->cd(ifile+1);

      float ring[7];
      float ringerr[7];
      float mean00[7];
      float mean01[7];
      float mean10[7];
      float mean11[7];
      float meanerr00[7];
      float meanerr01[7];
      float meanerr10[7];
      float meanerr11[7];

      for(int k=0;k<7;k++){
        ring[k] = k + 1;
        mean00[k] = mean[0][0][k];
        mean01[k] = mean[0][1][k];
        mean10[k] = mean[1][0][k];
        mean11[k] = mean[1][1][k];
        ringerr[k] = 0;
        meanerr00[k] = meanerr[0][0][k];
        meanerr01[k] = meanerr[0][1][k];
        meanerr10[k] = meanerr[1][0][k];
        meanerr11[k] = meanerr[1][1][k];
      }

      g00[ifile] = new TGraphErrors(7,ring,mean00,ringerr,meanerr00);
      g00[ifile]->SetMarkerStyle(20);
      g00[ifile]->GetXaxis()->SetTitle("Ring");
      g00[ifile]->GetYaxis()->SetTitle("<#Deltaw> (#mum)");
      g00[ifile]->SetTitle(Form("MinBias data (%s)",types.at(ifile).c_str()));
      g00[ifile]->SetMinimum(-50);
      g00[ifile]->SetMaximum(50);
      g00[ifile]->Draw("AP");

      g01[ifile] = new TGraphErrors(7,ring,mean01,ringerr,meanerr01);
      g01[ifile]->SetMarkerStyle(20);
      g01[ifile]->SetMarkerColor(2);
      g01[ifile]->SetLineColor(2);
      g01[ifile]->Draw("sameP");

      g10[ifile] = new TGraphErrors(7,ring,mean10,ringerr,meanerr10);
      g10[ifile]->SetMarkerStyle(20);
      g10[ifile]->SetMarkerColor(4);
      g10[ifile]->SetLineColor(4);
      g10[ifile]->Draw("sameP");

      g11[ifile] = new TGraphErrors(7,ring,mean11,ringerr,meanerr11);
      g11[ifile]->SetMarkerStyle(20);
      g11[ifile]->SetMarkerColor(6);
      g11[ifile]->SetLineColor(6);
      g11[ifile]->Draw("sameP");

      TLegend *leg=new TLegend(0.15,0.7,0.35,0.9);
      leg->AddEntry(g00[ifile],"Z- back-petal","p");
      leg->AddEntry(g01[ifile],"Z- front-petal","p");
      leg->AddEntry(g10[ifile],"Z+ back-petal","p");
      leg->AddEntry(g11[ifile],"Z+ front-petal","p");
      leg->SetBorderSize(1);
      leg->SetFillColor(0);
      leg->Draw();
    }

    if( nfiles > 1 ){
      TCanvas *dwdiff = new TCanvas("dwdiff","dwdiff",600,600);
      dwdiff->cd();

      TGraphErrors *gdiff00 = diffTGraph(g00[0],g00[1],"MinBias data deco-peak", "Ring", "<#Deltaw> (#mum)");
      gdiff00->SetMarkerStyle(20);
      gdiff00->SetMinimum(-25);
      gdiff00->SetMaximum(25);
      gdiff00->Draw("AP");

      TGraphErrors *gdiff01 = diffTGraph(g01[0],g01[1],"MinBias data deco-peak", "Ring", "<#Deltaw> (#mum)");
      gdiff01->SetMarkerStyle(20);
      gdiff01->SetMarkerColor(2);
      gdiff01->SetLineColor(2);
      gdiff01->Draw("sameP");

      TGraphErrors *gdiff10 = diffTGraph(g10[0],g10[1],"MinBias data deco-peak", "Ring", "<#Deltaw> (#mum)");
      gdiff10->SetMarkerStyle(20);
      gdiff10->SetMarkerColor(4);
      gdiff10->SetLineColor(4);
      gdiff10->Draw("sameP");

      TGraphErrors *gdiff11 = diffTGraph(g11[0],g11[1],"MinBias data deco-peak", "Ring", "<#Deltaw> (#mum)");
      gdiff11->SetMarkerStyle(20);
      gdiff11->SetMarkerColor(6);
      gdiff11->SetLineColor(6);
      gdiff11->Draw("sameP");

    }

  }
  
  if(draw("du_dw",cantitles)){
    cout << "Plotting du, dw histos... " << endl;

    can[idx]=new TCanvas(Form("can_%i",idx),"du_dw",1200,450);
    can[idx]->Divide(2,1);
    
    can[idx]->cd(1);
    setStats(hdu,nfiles,0.8,0.65,0.15);  
    plotHists(hdu,nfiles,"TOB","#Delta u [#mum]",-1000,1000,10,3);
    leg1->Draw();

    can[idx]->cd(2);
    setStats(hdw,nfiles,0.8,0.65,0.15);  
    plotHists(hdw,nfiles,"TOB","#Delta w [#mum]",-2000,2000,10,3);
    leg1->Draw();
        
    idx++;
  }


  //du vs. delta tan(theta) split v+/v- histos-----------------------------------
  
  if(draw("duvsdth",cantitles)){
    cout << "Plotting du vs. delta tan theta split v+/v- histos... " << endl;

    can[idx]=new TCanvas(cantitles.at(idx).c_str(),cantitles.at(idx).c_str(),1200,450*nfiles);
    can[idx]->Divide(nfiles,2);
    
    for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++){
      
      plotHistsTH2(hduthp[ifile],
                   Form("TOB %s (v+)",types.at(ifile).c_str()),"tan(#theta_{trk})-tan(#theta_{L})",
                   "#Delta u [#mum]",-0.3,0.5,-200,200,0,can[idx],ifile+1);
      
      plotHistsTH2(hduthm[ifile],
                   Form("TOB %s (v-)",types.at(ifile).c_str()),"tan(#theta_{trk})-tan(#theta_{L})",
                        "#Delta u [#mum]",-0.5,0.3,-200,200,0,can[idx],ifile+2);
      
    }
  
    idx++;

  }
  


  //TGraph du vs. delta tan theta----------------------------------------

  if(draw("duvsdth_tgraph",cantitles)){

    cout << "Plotting du vs. delta tan theta split v+/v- TGraphs... " << endl;

    can[idx]=new TCanvas(cantitles.at(idx).c_str(),cantitles.at(idx).c_str(),1200,450);
    can[idx]->Divide(2,1);
  
    vector<float> thetabinsp;
    vector<float> thetabinsm;

    for(int ibin=14;ibin<=30;ibin++)    thetabinsp.push_back(ibin*0.05-1);
    for(int ibin=10;ibin<=26;ibin++)    thetabinsm.push_back(ibin*0.05-1);
    
 
    TLatex *t=new TLatex();
    t->SetNDC();

    TGraphErrors *gduthp[nfiles];
    TGraphErrors *gduthm[nfiles];
        
    TF1 *fduthp[nfiles];
    TF1 *fduthm[nfiles];

    stringstream sduthp1[nfiles];
    stringstream sduthp2[nfiles];
    
    stringstream sduthm1[nfiles];
    stringstream sduthm2[nfiles];

    bool fit = false;

    for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++ ){
      
      can[idx]->cd(1);
      
      hduthp[ifile]->SetName(Form("hduthp_%i",ifile));
      gduthp[ifile] = getTGraphFromTH2(hduthp[ifile],thetabinsp, colors[ifile] );
      fduthp[ifile]=new TF1(Form("fduthp_%i",ifile),"pol1");
      fduthp[ifile]->SetLineColor( colors[ifile] );
      if(fit) gduthp[ifile]->Fit(fduthp[ifile]);
         
      float dwp     = fduthp[ifile]->GetParameter(1);
      float dwperr  = fduthp[ifile]->GetParError(1);
      float bp      = fduthp[ifile]->GetParameter(0);
      float bperr   = fduthp[ifile]->GetParError(0);
      
      sduthp1[ifile] << "#DeltaW = " << fround(dwp,3) << " #pm " << fround(dwperr,3) << " #mum" << endl;
      sduthp2[ifile] << "#Deltatan(LA) = " << fround(bp/(235.-dwp),3) << " #pm " << fround(bperr/(235.-dwp),4) << endl;
      
      gduthp[ifile]->SetTitle("TOB (V+)");
      gduthp[ifile]->Draw("AP");
      gduthp[ifile]->GetXaxis()->SetTitle("<tan(#theta_{trk})-tan(#theta_{LA})>");
      gduthp[ifile]->GetYaxis()->SetTitle("<#Deltau> [#mum]");
      //gduthp[ifile]->GetXaxis()->SetLimits(-1,1);
      gduthp[ifile]->SetMinimum(-25);
      gduthp[ifile]->SetMaximum(20);
      
      
      can[idx]->cd(2);

      hduthm[ifile]->SetName(Form("hduthm_%i",ifile));
      gduthm[ifile] = getTGraphFromTH2(hduthm[ifile],thetabinsm, colors[ifile] );
      fduthm[ifile]=new TF1(Form("fduthm_%i",ifile),"pol1");
      fduthm[ifile]->SetLineColor( colors[ifile] );
      if(fit) gduthm[ifile]->Fit(fduthm[ifile]);
         
      float dwm     = fduthm[ifile]->GetParameter(1);
      float dwmerr  = fduthm[ifile]->GetParError(1);
      float bm      = fduthm[ifile]->GetParameter(0);
      float bmerr   = fduthm[ifile]->GetParError(0);
      
      sduthm1[ifile] << "#DeltaW = " << fround(dwm,3) << " #pm " << fround(dwmerr,3) << " #mum" << endl;
      sduthm2[ifile] << "#Deltatan(LA) = " << fround(bm/(235.-dwm),3) << " #pm " << fround(bmerr/(235.-dwm),4) << endl;
      
      gduthm[ifile]->SetTitle("TOB (V-)");
      gduthm[ifile]->Draw("AP");
      gduthm[ifile]->GetXaxis()->SetTitle("<tan(#theta_{trk})-tan(#theta_{LA})>");
      gduthm[ifile]->GetYaxis()->SetTitle("<#Deltau> [#mum]");
      //gduthm[ifile]->GetXaxis()->SetLimits(-1,1);
      gduthm[ifile]->SetMinimum(-35);
      gduthm[ifile]->SetMaximum(15);
      //leg1->Draw();
    }

    for(unsigned int ifile = 0 ; ifile < nfiles ; ifile++ ){
      
      can[idx]->cd(1);
      if( ifile == 0 )   gduthp[ifile]->Draw("AP");
      else               gduthp[ifile]->Draw("sameP");
      leg1->Draw();

      if(fit){
        t->SetTextColor( colors[ifile] );
        t->DrawLatex(0.15,0.85-ifile*0.2,sduthp1[ifile].str().c_str());
        t->DrawLatex(0.15,0.75-ifile*0.2,sduthp2[ifile].str().c_str());
      }

      can[idx]->cd(2);
      if( ifile == 0 )   gduthm[ifile]->Draw("AP");
      else               gduthm[ifile]->Draw("sameP");
      //leg1->Draw();

      if(fit){
        t->SetTextColor( colors[ifile] );
        t->DrawLatex(0.15,0.85-ifile*0.2,sduthm1[ifile].str().c_str());
        t->DrawLatex(0.15,0.75-ifile*0.2,sduthm2[ifile].str().c_str());
      }
    
    }

    if( filenames.size() > 1 ){
    TCanvas *cdiff = new TCanvas("cdiff","duvsdtanthetadiff",1200,450);
    cdiff->Divide(2,1);

    cdiff->cd(1);


    TGraphErrors *gduthpdiff = diffTGraph(gduthp[0],gduthp[1],"TOB DECO - PEAK (v+)","<tan(#theta_{trk})-tan(#theta_{LA})>","<#Deltau> (#mum)");
    TF1* fduthpdiff=new TF1("fduthpdiff","pol1",-0.5,0.5);
    gduthpdiff->Fit(fduthpdiff,"R");
    gduthpdiff->Draw("AP");
    
    float dwp     = fduthpdiff->GetParameter(1);
    float dwperr  = fduthpdiff->GetParError(1);
    float bp      = fduthpdiff->GetParameter(0);
    float bperr   = fduthpdiff->GetParError(0);
    
    stringstream sduthp1diff;
    stringstream sduthp2diff;

    sduthp1diff << "#DeltaW = " << fround(dwp,3) << " #pm " << fround(dwperr,3) << " #mum" << endl;
    sduthp2diff << "#Deltatan(LA) = " << fround(bp/(235.-dwp),3) << " #pm " << fround(bperr/(235.-dwp),4) << endl;
  
    t->SetTextColor(1);
    t->DrawLatex(0.15,0.85,sduthp1diff.str().c_str());
    t->DrawLatex(0.15,0.75,sduthp2diff.str().c_str());

    cdiff->cd(2);

    TGraphErrors *gduthmdiff = diffTGraph(gduthm[0],gduthm[1],"TOB DECO - PEAK (v-)","<tan(#theta_{trk})-tan(#theta_{LA})>","<#Deltau> (#mum)");
    TF1* fduthmdiff=new TF1("fduthmdiff","pol1",-0.5,0.5);
    gduthmdiff->Fit(fduthmdiff,"R");
    gduthmdiff->Draw("AP");

    float dwm     = fduthmdiff->GetParameter(1);
    float dwmerr  = fduthmdiff->GetParError(1);
    float bm      = fduthmdiff->GetParameter(0);
    float bmerr   = fduthmdiff->GetParError(0);
    
    stringstream sduthm1diff;
    stringstream sduthm2diff;

    sduthm1diff << "#DeltaW = " << fround(dwm,3) << " #pm " << fround(dwmerr,3) << " #mum" << endl;
    sduthm2diff << "#Deltatan(LA) = " << fround(bm/(235.-dwm),3) << " #pm " << fround(bmerr/(235.-dwm),4) << endl;
  
    t->DrawLatex(0.15,0.85,sduthm1diff.str().c_str());
    t->DrawLatex(0.15,0.75,sduthm2diff.str().c_str());

    idx++;
    }
  }

  
  for(unsigned int ican=0;ican<ncan;ican++){
    can[ican]->Modified();
    can[ican]->Update();
    if(printgif) can[ican]->Print(Form("plots/%s_TOB.gif",cantitles.at(ican).c_str()));
  }
  
}