// main function with the same name of this file
// it will be executed when you type ".x interpolation.C" in ROOT
void mySpline()
{
// Section 1. Draw the points on a canvas
TCanvas *c1 = new TCanvas("c1","interpolation",0,0,1000,800);
TGraph *g1 = new TGraph(n, xarray, yarray);// just plots points
g1->SetMarkerStyle(20);
g1->SetMarkerSize(2);
g1->GetXaxis()->SetLimits(-1, 3); // set real range
g1->GetXaxis()->SetRangeUser(-0.5, 2.5); // set visible range
g1->GetXaxis()->SetTitle("X");
g1->GetXaxis()->CenterTitle();
g1->GetYaxis()->SetLimits(-1, 2.0);
g1->GetYaxis()->SetRangeUser(-0.4, 1.6);
g1->GetYaxis()->SetTitle("Y");
g1->GetYaxis()->CenterTitle();
g1->Draw("ap"); // options to draw a graph are described on
// http://root.cern.ch/root/html/TGraph.html#TGraph:PaintGraph
// Section 3. Draw the Cubic Spline to the same canvas
TSpline3 *sp = new TSpline3("Cubic Spline", xarray, yarray, n, "b2e2", 0, 0);
// refer to http://root.cern.ch/root/html/Tspline3.html for the usage of TSpline3
// "b2e2" together with the last two "0" means that the second derivatives
// of the begin and end points equal to zero
sp->SetLineColor(kRed);
sp->Draw("lsame");
//____________________________________________________________________________
double LogChi2(double *x, double *par)
{
double ChiSqr = 0;
for (int i=0; i<n; i++)
ChiSqr = ChiSqr + (yarray[i] - x[0]*exp(-x[1]*xarray[i]))**2;
return log(ChiSqr);
}}
//Main of the program
void randSpline()
{
// making vectors
vector <double> xVector, yVector, xErrorVector, yErrorVector;
const int n = 6;
// filling vectors with random numbers
FillRandVectors(xVector, yVector, n);
// making graph
TGraph *gr = LoadGraphFromVectors(xVector, yVector);
//Plot
TCanvas *c1 = new TCanvas("c1","My Awesome Test Graph!!",200,10,700,500);
gr->Draw("apz");
// Draw the Cubic Spline to the same canvas
TSpline3 *sp = new TSpline3("Cubic Spline", gr,"b2e2", 0, 0);
sp->SetLineColor(kRed);
sp->Draw("lsame");
c1->Update();
}
void buildFakeAngTree(const Char_t* outtag,
const Float_t timereso, // ns
const UInt_t simevts=1,
const Float_t thetaOpt=400, // deg
const Float_t phiOpt=400, // deg
const Float_t coneOpt=400, // deg
const UInt_t rseed=23192,
const Float_t norm=100.0, // mV
const Float_t noise=20.0, // mV
const Char_t* outdir="/data/users/cjreed/work/simEvts",
const Char_t* infn="/w2/arianna/jtatar/nt.sigtemps.root",
const Char_t* geofn="/data/users/cjreed/work/"
"BounceStudy/Stn10/"
"CampSiteGeometry.root") {
// if any of the angles (thetaOpt, phiOpt, coneOpt) > 360, a random
// value will be used instead
//
// expect angles in the Templates tree to be in degrees
//
// expect the waveforms in the Templates tree to have amplitude 1
TRandom3 rnd(rseed);
geof = TFile::Open(geofn);
gg = dynamic_cast<TGeoManager*>(geof->Get("CampSite2013"));
site = dynamic_cast<const TSnGeoStnSite*>(gg->GetTopVolume());
TVector3 pos[NSnConstants::kNchans], nvec[NSnConstants::kNchans];
for (UChar_t ch=0; ch<NSnConstants::kNchans; ++ch) {
site->SetLPDAPosition(ch, pos[ch]);
site->SetLPDANormalVec(ch, nvec[ch]);
Printf("pos ch%d:",ch);
pos[ch].Print();
Printf("normal ch%d:",ch);
nvec[ch].Print();
}
TArrayD zeros(6);
inf = TFile::Open(infn);
nnt = dynamic_cast<TTree*>(inf->Get("Templates"));
TString infns(infn);
TString indir;
Int_t fl(0);
if (infns.Contains('/')) {
fl = infns.Last('/') + 1;
indir = infns(0, fl-1);
}
TString plaininfn = infns(fl, infns.Length()-fl);
TString outfn = Form("%s/FakeEvts.%s.%s", outdir, outtag,
plaininfn.Data());
outf = TFile::Open(outfn.Data(),"recreate");
outf->cd();
TParameter<Float_t> trp("TimeResolution", timereso);
trp.Write();
TParameter<Float_t> nmp("Normalization", norm);
nmp.Write();
TParameter<Float_t> nop("NoiseRMS", noise);
nop.Write();
TParameter<UInt_t> rsp("RandomSeed", rseed);
rsp.Write();
TSnCalWvData* wave = new TSnCalWvData;
Float_t eang(0), hang(0), hpf(0), limiter(0), coneang(0);
Bool_t bice(kFALSE);
nnt->SetBranchAddress("wave.",&wave);
nnt->SetBranchAddress("EAng",&eang);
nnt->SetBranchAddress("HAng",&hang);
nnt->SetBranchAddress("hpf",&hpf);
nnt->SetBranchAddress("limiter",&limiter);
nnt->SetBranchAddress("coneAng",&coneang);
nnt->SetBranchAddress("bIce",&bice);
// to look up waveform for EAng, HAng
nnt->BuildIndex("EAng + (1000*HAng)","coneAng");
// find the max angles
Printf("finding allowed angles...");
std::set<Float_t> Eangs, Hangs, Cangs;
const Long64_t nnents = nnt->GetEntries();
for (Long64_t i=0; i<nnents; ++i) {
nnt->GetEntry(i);
Eangs.insert(eang);
Hangs.insert(hang);
Cangs.insert(coneang);
}
#ifdef DEBUG
std::set<Float_t>::const_iterator ang, end = Eangs.end();
Printf("EAngs:");
for (ang=Eangs.begin(); ang!=end; ++ang) {
Printf("%g",*ang);
}
Printf("HAngs:");
for (ang=Hangs.begin(), end=Hangs.end(); ang!=end; ++ang) {
Printf("%g",*ang);
}
Printf("ConeAngs:");
for (ang=Cangs.begin(), end=Cangs.end(); ang!=end; ++ang) {
Printf("%g",*ang);
}
#endif
Float_t theta(0), phi(0), cone(0);
Float_t EAng[NSnConstants::kNchans],
HAng[NSnConstants::kNchans];
Float_t CAng(0);
TSnCalWvData* evdat = new TSnCalWvData;
TSnEventMetadata* meta = new TSnEventMetadata;
TSnEventHeader* hdr = new TSnEventHeader;
//ot = nnt->CloneTree(0);
//ot->SetName("SimTemplEvts");
ot = new TTree("SimTemplEvts","simulated events from templates",1);
ot->SetDirectory(outf);
ot->Branch("EventMetadata.",&meta);
ot->Branch("EventHeader.",&hdr);
ot->Branch("EAng",&(EAng[0]),Form("EAng[%hhu]/F",NSnConstants::kNchans));
ot->Branch("HAng",&(HAng[0]),Form("HAng[%hhu]/F",NSnConstants::kNchans));
ot->Branch("CAng",&CAng,"CAng/F");
ot->Branch("theta",&theta,"theta/F");
ot->Branch("phi",&phi,"phi/F");
ot->Branch("NuData.",&evdat);
// some useful aliases
TString an;
for (UChar_t ch=0; ch<NSnConstants::kNchans; ++ch) {
// to use as a cut for a particular channel:
an = Form("Ch%d",ch);
ot->SetAlias(an.Data(),
Form("(Iteration$>=(%hhu*%hhu)) && (Iteration$<(%hhu*%hhu))",
NSnConstants::kNsamps, ch,
NSnConstants::kNsamps,
static_cast<UChar_t>(ch+1)));
// to use as a variable showing the sample number [0,127] for any chan
an = Form("SmpCh%d",ch);
ot->SetAlias(an.Data(),
Form("Iteration$-%u", static_cast<UInt_t>(ch)
*static_cast<UInt_t>(NSnConstants::kNsamps)));
// e.g. Draw("RawData.fData:SmpCh2","EventHeader.fNum==21 && Ch2","l")
}
Printf("generating events...");
TStopwatch timer;
timer.Start();
for (UInt_t i=0; i<simevts; ++i) {
if ( (i%1000)==0 ) {
fprintf(stderr,"Processing %u/%u ... \r",i,simevts);
}
// choose angles
theta = (thetaOpt>360.) ? TMath::ACos( rnd.Uniform(-1.0, 0.0) )
: thetaOpt * TMath::DegToRad();
phi = (phiOpt>360.) ? rnd.Uniform(0.0, TMath::TwoPi())
: phiOpt * TMath::DegToRad();
cone = (coneOpt>360.)
? rnd.Uniform(*(Cangs.begin()), *(Cangs.rbegin()))
: coneOpt; // leave this one in degrees (as in the tree)
CAng = findNearestAllowedAngle(Cangs, cone);
#ifdef DEBUG
Printf("--- theta=%g, phi=%g, cone=%g",
theta*TMath::RadToDeg(), phi*TMath::RadToDeg(), cone);
#endif
// calculate channel shifts
TArrayD pwdt = NSnChanCorl::GetPlaneWaveOffsets(theta,
phi,
zeros,
pos,
kNgTopFirn);
TVector3 dir;
dir.SetMagThetaPhi(1.0, theta, phi);
#ifdef DEBUG
TObjArray graphs;
graphs.SetOwner(kTRUE);
TCanvas* c1 = new TCanvas("c1","c1",800,700);
c1->Divide(2,2);
#endif
for (UChar_t ch=0; ch<NSnConstants::kNchans; ++ch) {
// look up the EAng, fhang for this antenna
Float_t feang(0), fhang(0);
findEangHang(nvec[ch], dir, feang, fhang);
feang = TMath::Abs(TVector2::Phi_mpi_pi(feang));
fhang = TMath::Abs(TVector2::Phi_mpi_pi(fhang));
feang *= TMath::RadToDeg();
fhang *= TMath::RadToDeg();
// find closest allowed angle
EAng[ch] = findNearestAllowedAngle(Eangs, feang);
HAng[ch] = findNearestAllowedAngle(Hangs, fhang);
const Long64_t ni =
nnt->GetEntryNumberWithIndex(EAng[ch] + (1000*HAng[ch]), CAng);
#ifdef DEBUG
Printf("EAng=%g (%g), HAng=%g (%g), CAng=%g, ni=%lld",
EAng[ch],feang,HAng[ch],fhang,CAng,ni);
#endif
if (ni>-1) {
nnt->GetEntry(ni);
#ifdef DEBUG
c1->cd(ch+1);
TGraph* och = wave->NewGraphForChan(0, kTRUE);
const Int_t ochnp = och->GetN();
Double_t* ochy = och->GetY();
for (Int_t k=0; k<ochnp; ++k, ++ochy) {
*ochy *= norm;
}
graphs.Add(och);
och->SetLineColor(kBlack);
och->SetMarkerColor(kBlack);
och->SetMarkerStyle(7);
och->Draw("apl");
#endif
// first calculate the shift between chans due to the angle
// ch0 is always unshifted; other chans shifted w.r.t. ch0
// jitter the shift by the specified timing resolution
const Double_t shift =
rnd.Gaus( (ch==0) ? 0.0
: -pwdt.At( TSnRecoChanOffsets::IndexFor(ch, 0) ),
timereso);
// get a graph of the waveform
// data only in channel 0 of the template
TGraph* gch = wave->NewGraphForChan(0, kTRUE);
// "fit" the graph with an spline interpolation
TSpline3* gsp = new TSpline3("stmp", gch);
// evaluate the spline at the new sample positions
// (shifted, but NOT wrapped)
// and save that into the event data waveform
Float_t* d = evdat->GetData(ch);
const Float_t tstep = 1.0 / NSnConstants::kSampRate;
const Float_t tlast = static_cast<Float_t>(NSnConstants::kNsamps-1)
/ NSnConstants::kSampRate;
Float_t xloc = shift;
for (UChar_t s=0; s<NSnConstants::kNsamps; ++s, ++d, xloc+=tstep) {
if ( (xloc<0.0) || (xloc>=tlast) ) {
*d = 0.0;
} else {
*d = gsp->Eval( xloc );
}
}
#ifdef DEBUG
Printf("ch%hhu: shift=%g, dt=%g", ch, shift,
(ch==0) ? 0.0
: pwdt.At( TSnRecoChanOffsets::IndexFor(ch, 0) ));
TGraph* fch = evdat->NewGraphForChan(ch, kTRUE);
Double_t* y = gch->GetY();
Double_t* fy = fch->GetY();
for (UChar_t s=0; s<NSnConstants::kNsamps; ++s, ++y, ++fy) {
*y *= norm;
*fy *= norm;
}
gch->SetLineColor(kRed+1);
gch->SetMarkerColor(kRed+1);
gch->SetMarkerStyle(7);
gch->Draw("pl");
delete gsp;
gsp = new TSpline3("stmp",gch);
gsp->SetLineColor(kAzure-6);
gsp->SetMarkerColor(kAzure-6);
gsp->SetMarkerStyle(7);
gsp->Draw("pl same");
graphs.Add(fch);
fch->SetLineColor(kOrange+7);
fch->SetMarkerColor(kOrange+7);
fch->SetMarkerStyle(7);
fch->Draw("pl");
#endif
d = evdat->GetData(ch);
// finally add noise to the waveform
for (UChar_t s=0; s<NSnConstants::kNsamps; ++s, ++d) {
*d = rnd.Gaus( (*d) * norm, noise );
}
#ifdef DEBUG
TGraph* nch = evdat->NewGraphForChan(ch, kTRUE);
graphs.Add(nch);
nch->SetLineColor(kGreen+2);
nch->SetMarkerColor(kGreen+2);
nch->SetMarkerStyle(7);
nch->Draw("pl");
#endif
// cleanup
#ifdef DEBUG
graphs.Add(gch);
graphs.Add(gsp);
#else
delete gch;
delete gsp;
#endif
}
} // end channel loop
#ifdef DEBUG
TObject* o(0);
while ( (o=c1->WaitPrimitive())!=0 ) {
gSystem->ProcessEvents();
}
delete c1;
#endif
// save this event
ot->Fill();
} // end event loop
fprintf(stderr,"\n");
timer.Stop();
Printf("Finished generating events in:");
timer.Print();
outf->Write();
Printf("Wrote [%s]",outf->GetName());
delete outf; outf=0; // close file
}
void CommandMSUGRA(TString plotName_,Int_t tanBeta_, Bool_t plotLO_){
gStyle->SetOptTitle(0);
gStyle->SetOptStat(0);
gStyle->SetPalette(1);
gStyle->SetTextFont(42);
//convert tanb value to string
std::stringstream tmp;
tmp << tanBeta_;
TString tanb( tmp.str() );
// Output file
cout << " create " << plotName_ << endl;
TFile* output = new TFile( plotName_, "RECREATE" );
if ( !output || output->IsZombie() ) { std::cout << " zombie alarm output is a zombie " << std::endl; }
//-----------------------------------
//set old exclusion Limits
//-----------------------------------
TGraph* LEP_ch = set_lep_ch(tanBeta_);
TGraph* LEP_sl = set_lep_sl(tanBeta_); //slepton curve
TGraph* TEV_sg_cdf = set_tev_sg_cdf(tanBeta_); //squark gluino cdf
TGraph* TEV_sg_d0 = set_tev_sg_d0(tanBeta_); //squark gluino d0
//TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(tanBeta_); //trilepton cdf
//TGraph* TEV_tlp_d0 = set_tev_tlp_d0(tanBeta_); //trilepton d0
TGraph* stau = set_tev_stau(tanBeta_); //stau
//TGraph* NoEWSB = set_NoEWSB(tanBeta_);
TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(tanBeta_);
TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(tanBeta_);
//-----------------------------------
// constant sqquark and gluino lines
//-----------------------------------
const unsigned int nlines = 4;
TF1* lnsq[nlines];
TF1* lngl[nlines];
TLatex* sq_text[nlines];
TLatex* gl_text[nlines];
for(unsigned int i = 0; i < nlines; i++){
lnsq[i] = constant_squark(tanBeta_,i);
sq_text[i] = constant_squark_text(i,*lnsq[i],tanBeta_);
lngl[i] = constant_gluino(tanBeta_,i);
gl_text[i] = constant_gluino_text(i,*lngl[i]);
}
//-----------------------------------
// Legends
//-----------------------------------
TLegend* legst = makeStauLegend(0.05,tanBeta_);
TLegend* legexp = makeExpLegend( *TEV_sg_cdf,*TEV_sg_d0,*LEP_ch,*LEP_sl,*TEV_sn_d0_1,0.035,tanBeta_);
//TLegend* legNoEWSB = makeNoEWSBLegend(0.05,tanBeta_);
//-----------------------------------
// make Canvas
//-----------------------------------
TCanvas* cvsSys = new TCanvas("cvsnm","cvsnm",0,0,800,600);
gStyle->SetOptTitle(0);
cvsSys->SetFillColor(0);
cvsSys->GetPad(0)->SetRightMargin(0.07);
cvsSys->Range(-120.5298,26.16437,736.0927,500);
//cvsSys->Range(-50.5298,26.16437,736.0927,500);
cvsSys->SetFillColor(0);
cvsSys->SetBorderMode(0);
cvsSys->GetPad(0)->SetBorderSize(2);
cvsSys->GetPad(0)->SetLeftMargin(0.1407035);
cvsSys->GetPad(0)->SetTopMargin(0.08);
cvsSys->GetPad(0)->SetBottomMargin(0.13);
cvsSys->SetTitle("tan#beta="+tanb);
output->cd();
//and now
//the exclusion limits
TGraphErrors* First ;
//TGraphErrors* FirstDummy ;
TGraphErrors* Second;
TGraphErrors* Third;
TGraphErrors* Second_up;
TGraphErrors* Second_low;
TGraphErrors* expband;
TGraphErrors* obs2010;
if (tanBeta_ == 3) {
//First = getObserved_NLOunc();
//FirstDummy = getObserved_NLOunc();
} else {
//First = getNLOobsTanbeta10();
//First = getNLOobsTanbeta10_smooth();
//First = getNLOobsTanbeta10_funky();
First = getNLOobsTanbeta10();
//FirstDummy = getObserved_NLOunc();
//Second = getNLOexpTanbeta10();
Second = getNLOexpTanbeta10();
Second_up = getNLOexpUpTanbeta10();
Second_low = getNLOexpDownTanbeta10();
Third = getNLOexpTanbeta10();
expband = getNLOexpTanbeta10_band();
obs2010 = getNLOobsTanbeta10_2010();
//Second_up = getExpected_NLO_tanBeta3_up();
//Second_low = getExpected_NLO_tanBeta3_low();
}
//Third = getExpected_NLOunc();//getLO_jetMultis();
//Second = getLO_signalCont();
TGraph *grObserved3p5_shape = getObserved3p5_shape();
TGraph *grExpected3p5_shape = getExpected3p5_shape();
// First->SetMarkerColor(kWhite);
// First->GetXaxis()->SetRangeUser(2.,500.);
// First->GetYaxis()->SetRangeUser(80,500);
// if(tanBeta_ == 50) First->GetXaxis()->SetRangeUser(200,500);
// First->GetXaxis()->SetTitle("m_{0} (GeV)");
// First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
// First->GetYaxis()->SetTitleOffset(0.8);
double m0min = 0;
if (tanBeta_ == 50) m0min=200;
TH2D* hist = new TH2D("h","h",100,m0min,m0max,100,120,m12max);
hist->Draw();
hist->GetXaxis()->SetTitle("m_{0} (GeV/c^{2})");
hist->GetYaxis()->SetTitle("m_{1/2} (GeV/c^{2})");
hist->GetYaxis()->SetTitleOffset(1.);
hist->GetXaxis()->SetNdivisions(506);
// if (tanBeta_ == 50) hist->GetXaxis()->SetNdivisions(504);
hist->GetYaxis()->SetNdivisions(506);
//int col[]={2,3,4};
//TFile *f = TFile::Open("exclusion_Spring11_CLs.root");
//TFile *f = TFile::Open("exclusion_Fall10_tcmet_JPT.root");
//TFile *f = TFile::Open("exclusion_Fall10_pfmet_pfjets.root");
//TFile *f = new TFile("exclusion_Fall10_pfmet_pfjets_CLs.root");
//TH2F* h = (TH2F*) f->Get("hexcl_NLO_obs");
//TH2F* h = (TH2F*) f->Get("hexcl_NLO_exp");
//TH2F* h = (TH2F*) f->Get("hexcl_NLO_expp1");
//TH2F* h = (TH2F*) f->Get("hexcl_NLO_expm1");
//h->SetMaximum(3);
//h->Draw("samecolz");
TSpline3 *sFirst = new TSpline3("sFirst",First);
sFirst->SetLineColor(kRed);
sFirst->SetLineWidth(3);
First->SetLineColor(kRed);
First->SetLineWidth(3);
TSpline3 *sSecond = new TSpline3("sSecond",Second);
sSecond->SetLineColor(kBlue);
sSecond->SetLineStyle(2);
sSecond->SetLineWidth(3);
Second->SetLineColor(kBlue);
Second->SetLineStyle(2);
Second->SetLineWidth(3);
TSpline3 *sSecond_up = new TSpline3("sSecond_up",Second_up);
sSecond_up->SetLineColor(kCyan);
sSecond_up->SetLineStyle(1);
sSecond_up->SetLineWidth(3);
Second_up->SetLineColor(kBlue);
//Second_up->SetLineColor(1);
Second_up->SetLineWidth(2);
TSpline3 *sSecond_low = new TSpline3("sSecond_low",Second_low);
sSecond_low->SetLineColor(kCyan);
sSecond_low->SetLineStyle(1);
sSecond_low->SetLineWidth(3);
Second_low->SetLineColor(kBlue);
//Second_low->SetLineColor(1);
Second_low->SetLineWidth(2);
Third->SetLineColor(kCyan);
Third->SetLineWidth(30);
// TSpline3 *sThird = new TSpline3("sThird",Third);
// sThird->SetLineColor(kGreen+2);
// sThird->SetLineStyle(4);
// sThird->SetLineWidth(3);
// Third->SetLineColor(kGreen+2);
// Third->SetLineStyle(4);
// Third->SetLineWidth(3);
// First->Draw("AP");
/*
for(vector<TH1F*>::iterator at = exclusionPlots.begin();at != exclusionPlots.end();++at){
(*at)->SetContour(2);
if(n == 0){
(*at)->DrawCopy();
(*at)->SetTitle("tan#beta="+tanBeta_);
}
cout << " n " << n << endl;
(*at)->DrawCopy("same");
// (*it)->Write();
cout << " here " << endl;
n++;
}*/
TLegend* myleg;
if( plotLO_ ) myleg = new TLegend(0.3,0.75,0.54,0.9,NULL,"brNDC");
else myleg = new TLegend(0.25,0.75,0.54,0.9,NULL,"brNDC");
myleg->SetFillColor(0);
myleg->SetShadowColor(0);
myleg->SetTextSize(0.03);
myleg->SetBorderSize(0);
TH1F* hdummy = new TH1F();
hdummy->SetLineColor(4);
hdummy->SetFillColor(4);
hdummy->SetFillStyle(3002);
hdummy->SetLineWidth(2);
hdummy->SetLineStyle(2);
// myleg->AddEntry(sSecond,"NLO Expected Limit","L");
if (tanBeta_ == 3 && plotLO_) {
myleg->AddEntry(sSecond,"LO Observed Limit","L");
myleg->AddEntry(sFirst,"NLO Observed Limit","L");
} else {
//myleg->AddEntry(sFirst,"CMS OS Dilepton Limit","L");
myleg->AddEntry(sFirst,"NLO observed limit","L");
//myleg->AddEntry(hdummy,"NLO expected limit","LF");
myleg->AddEntry(hdummy,"NLO expected limit","L");
//myleg->AddEntry(sSecond,"NLO expected limit","L");
//myleg->AddEntry(sSecond_up,"NLO expected limit (+/-1#sigma)","L");
myleg->AddEntry(obs2010,"2010 NLO observed limit","L");
}
//sSecond_up->Draw("h same");
//sSecond_low->Draw("h same");
//constant squark and gluino mass contours
for (unsigned int it=1;it<nlines;it++) {
lngl[it]->Draw("same");
lnsq[it]->Draw("same");
sq_text[it]->Draw();
gl_text[it]->Draw();
}
sSecond_up->SetFillStyle(4010);
sSecond_up->SetFillColor(kCyan-10);
sSecond_low->SetFillStyle(1001);
sSecond_low->SetFillColor(10);
//expected and observed (LO & NLO) contours
//sFirst->Draw("same");
//sSecond->Draw("same");
//sThird->Draw("same");
//Third->Draw("samec");
//expband->Draw("samecf"); // summer11 expected band
//First->Draw("samec"); // summer11 observed exclusion
//First->SetMarkerColor(1);
//First->Draw("samep");
//Second->Draw("samec"); // summer11 expected limit
obs2010->Draw("samec");
//Second_up->Draw("samec");
//Second_low->Draw("samec");
grObserved3p5_shape->SetLineColor(2);
grExpected3p5_shape->SetLineColor(4);
grExpected3p5_shape->SetLineStyle(2);
grObserved3p5_shape->Draw("same");
grExpected3p5_shape->Draw("same");
// if (tanBeta_ == 3) Third->Draw("samec");
//if (tanBeta_ == 3 && plotLO_) Second->Draw("samec");
//exclusion limits previous experiments
if(tanBeta_ == 3){
TEV_sn_d0_1->Draw("fsame");
TEV_sn_d0_2->Draw("fsame");
}
LEP_ch->Draw("fsame");
if (tanBeta_ != 50) LEP_sl->Draw("fsame");
//remove CDF/D0 excluded regions
TEV_sg_cdf->Draw("fsame");
TEV_sg_d0->Draw("same");
TEV_sg_d0->Draw("fsame");
//other labels
Double_t xpos = 0;
Double_t xposi = 0;
Double_t ypos = 0;
if(tanBeta_ == 50) xposi = 100;
if(tanBeta_ == 50) xpos = 200;
if(tanBeta_ == 50) ypos = -10;
//TLatex* lumilabel = new TLatex(135.+xposi,510.,"L_{int} = 34 pb^{-1}, #sqrt{s} = 7 TeV");
//TLatex* lumilabel = new TLatex(305.+xposi + 100,510.,"L_{int} = 976 pb^{-1}, #sqrt{s} = 7 TeV");
TLatex* lumilabel = new TLatex(490,m12max+15,"#sqrt{s} = 7 TeV, #scale[0.6]{#int} L dt = 3.5 fb^{-1}");
lumilabel->SetTextSize(0.05);
lumilabel->Draw("same");
TLatex* cmslabel = new TLatex(10.,m12max+15,"CMS Preliminary");
cmslabel->SetTextSize(0.05);
cmslabel->Draw("same");
TString text_tanBeta;
//text_tanBeta = "tan#beta = "+tanb+", A_{0} = 0, sign(#mu) > 0";
text_tanBeta = "tan#beta = "+tanb+", A_{0} = 0, #mu > 0";
//TLatex* cmssmpars = new TLatex(70.+xpos,340.+ypos,text_tanBeta);
TLatex* cmssmpars = new TLatex(120,540,text_tanBeta);
//TLatex* cmssmpars = new TLatex(200,370,text_tanBeta);
cmssmpars->SetTextSize(0.045);
cmssmpars->Draw("same");
//LM points
TMarker* LM0 = new TMarker(200.,160.,20);
TMarker* LM1 = new TMarker(60.,250.,20);
TMarker* LM3 = new TMarker(330.,240.,20);
TMarker* LM6 = new TMarker(80.,400.,20);
LM0->SetMarkerSize(1.2);
LM1->SetMarkerSize(1.2);
TLatex* tLM0 = new TLatex(205.,160.," LM0");
tLM0->SetTextSize(0.035);
TLatex* tLM1 = new TLatex(80.,245.,"LM1");
tLM1->SetTextSize(0.035);
//TLatex* tLM3 = new TLatex(350.,235.,"LM3 (tan#beta=20)");
TLatex* tLM3 = new TLatex(350.,235.,"LM3");
tLM3->SetTextSize(0.035);
TLatex* tLM6 = new TLatex(100.,395.,"LM6");
tLM6->SetTextSize(0.035);
// if (tanBeta_ != 50){
// LM0->Draw("same");
// tLM0->Draw("same");
// LM1->Draw("same");
// tLM1->Draw("same");
// }
if (tanBeta_ == 10){
LM1->Draw("same");
tLM1->Draw("same");
LM3->Draw("same");
tLM3->Draw("same");
LM6->Draw("same");
tLM6->Draw("same");
}
/*
Int_t n = 0;
for(vector<TH1F*>::iterator at = exclusionPlots.begin();at != exclusionPlots.end();++at){
(*at)->SetContour(2);
if(n == 0){
(*at)->DrawCopy("same");
(*at)->SetTitle("tan#beta=3");
}
cout << " n " << n << endl;
(*at)->DrawCopy("same");
// (*it)->Write();
cout << " here " << endl;
n++;
}
*/
//stau=LSP contour
stau->Draw("fsame");
//NoEWSB->Draw("fsame");
//legends
legexp->Draw();
legst->Draw();
myleg->Draw();
//legNoEWSB->Draw();
//First->Draw("samec");
// if (tanBeta_ == 3) Third->Draw("samec");
//if (tanBeta_ == 3 && plotLO_) Second->Draw("samec");
hist->Draw("sameaxis");
cvsSys->RedrawAxis();
cvsSys->Update();
cvsSys->Write();
if( plotLO_ ){
cvsSys->SaveAs("RA6_ExclusionLimit_tanb"+tanb+"_LO.pdf");
cvsSys->SaveAs("RA6_ExclusionLimit_tanb"+tanb+"_LO.png");
}else{
cvsSys->SaveAs("RA6_ExclusionLimit_tanb"+tanb+".eps");
cvsSys->SaveAs("RA6_ExclusionLimit_tanb"+tanb+".pdf");
cvsSys->SaveAs("RA6_ExclusionLimit_tanb"+tanb+".png");
}
output->Write();
//output->Close();
//delete output;
}
