/////////////////////////////////////////////////////////
// This class has been automatically generated on
// Thu Oct  1 20:13:36 2015 by ROOT version 5.34/10
// from TTree Clusters/TbUT nTuple
//////////////////////////////////////////////////////////

#ifndef ClusterWithTrackAna_h
#define ClusterWithTrackAna_h

#include <TROOT.h>
#include <TChain.h>
#include <TFile.h>
#include <TProfile.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TLegend.h>
#include <TBox.h>
#include <TLine.h>

#include "AnalysisBase.h"
#include "CMS.h"

//#include "ClusterWithTrackAna_Inputs.h"

// Header file for the classes stored in the TTree if any.
#include <vector>
#include <iostream>
#include <algorithm>

// Fixed size dimensions of array or collections stored in the TTree if any.

class ClusterWithTrackAna  : public AnalysisBase{
public :
   
  CMS            *fCMS;

  ClusterWithTrackAna(int biasVal);
  //ClusterWithTrackAna(TTree *tree=0);
   virtual ~ClusterWithTrackAna();
   virtual void     Loop();
  virtual TString getFileBase(TString scan="Bias", TString board="", TString bias="", TString angle="0", TString sector="1");  

  TFile *fout;
};

#endif

#ifdef ClusterWithTrackAna_cxx
ClusterWithTrackAna::ClusterWithTrackAna(int biasVal){
  //ClusterWithTrackAna::ClusterWithTrackAna(TTree *tree){
// if parameter tree is not specified (or zero), connect the file
// used to generate this class and read the Tree.
  TString filename;
  if(biasVal != 0) m_bias = Form("%d",biasVal);
  TTree* tree = 0;
  
   if (tree == 0) {
     //
     // lhcb-dev
     //
     TString scanType = "Run_Bias";
     filename = m_fileIndir+"Board"+m_board+"/"+scanType+"_Scan-"+getFileBase("Bias",m_board, m_bias, m_angle,m_sector)+"_Tuple_Tracks.root"; 
     if(m_scanType.Contains("Angle")){
       scanType = "Run_Angle";
       filename = m_fileIndir+"Board"+m_board+"_Angle/"+scanType+"_Scan-"+getFileBase("Bias",m_board, m_bias, m_angle,m_sector)+"_Tuple_Tracks.root"; 
     }
     
     cout << "===> Opening file: " << filename << endl;
     

     // lxplus
     //
     //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA8/Run_Bias_Scan-B8-A-296-13332_Tuple_Tracks.root"; // A8 - s1
     //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA8/Run_Bias_Scan-B8-A-324-13355_Tuple_Tracks.root"; // A8 - s2
     //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA8/Run_Bias_Scan-B8-A-359-13386_Tuple_Tracks.root";  // A8 - s3
     //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA6/Run_Bias_Scan-B6-A-212-8358_Tuple_Tracks.root";   // A6
     //

     TFile *f = new TFile(filename);
     tree = (TTree*)f->Get("Clusters");
     TH1F *hn = (TH1F*)f->Get("hWidthNoise");
     for(int i=0;i<nChan;i++){
       noise[i] = hn->GetBinContent(i+1);
     }

     // Get mean noise and width
     hMeanNoise 	= (TH1F*)f->Get("hMeanNoise"); 
     hWidthNoise 	= (TH1F*)f->Get("hWidthNoise"); 
     cout << "===> Opening file: " << filename << endl;

     TString filename2 = filename.ReplaceAll("_Tracks","");
     TFile * f2 = new TFile(filename2);
     if(f2) {
       TTree * tree2 = (TTree*) f2->Get("TbUT/CMS");
       if(tree2) {
         tree->AddFriend(tree2);
         fCMS = new CMS(tree2);
       }else{
         fCMS=0;  
       }
     } else {
       cout << "WARNING: Could not find CMS data file: " << filename2 << endl;
       fCMS =0 ;
     }
   }
   Init(tree);
   


}

ClusterWithTrackAna::~ClusterWithTrackAna()
{
   if (!fChain) return;
   delete fChain->GetCurrentFile();
}
TString ClusterWithTrackAna::getFileBase(TString scan, TString board, TString bias, TString angle, TString sector) {

  TString tag = "";

  if(board.Contains("A6")){
    if(sector == "1") {
      if ( bias == "300" ) tag = "B6-A-212-8358";
      if ( bias == "250" ) tag = "B6-A-213-8359";
      if ( bias == "200" ) tag = "B6-A-214-8361";
      if ( bias == "150" ) tag = "B6-A-215-8362";
      if ( bias == "100" ) tag = "B6-A-216-8363";
      if ( bias == "75" )  tag = "B6-A-219-8365";
      if ( bias == "50" )  tag = "B6-A-220-8366";
    } else if(sector == "2"){
      if ( bias == "300" ) tag = "B6-A-242-8389";
      if ( bias == "250" ) tag = "B6-A-245-8392";
      if ( bias == "200" ) tag = "B6-A-246-8393";
      if ( bias == "150" ) tag = "B6-A-247-8394";
      if ( bias == "100" ) tag = "B6-A-248-8395";
      if ( bias == "75" )  tag = "B6-A-249-8396";
      if ( bias == "50" )  tag = "B6-A-250-8397";
    } else if(sector == "3"){
      if ( bias == "300" ) tag = "B6-A-293-8425";
      if ( bias == "250" ) tag = "B6-A-294-8426";
      if ( bias == "200" ) tag = "B6-A-295-8427";
      if ( bias == "150" ) tag = "B6-A-296-8428";
      if ( bias == "100" ) tag = "B6-A-297-8429";
      if ( bias == "75" )  tag = "B6-A-298-8430";
      if ( bias == "50" )  tag = "B6-A-299-8431";
    } else if(sector == "4"){
      if ( bias == "300" ) tag = "B6-A-326-8452";
      if ( bias == "250" ) tag = "B6-A-327-8453";
      if ( bias == "200" ) tag = "B6-A-328-8454";
      if ( bias == "150" ) tag = "B6-A-329-8455";
      if ( bias == "100" ) tag = "B6-A-330-8456";
      if ( bias == "75" )  tag = "B6-A-331-8457";
      if ( bias == "50" )  tag = "B6-A-332-8458";
    } else if(sector == "5"){
      if ( bias == "300" ) tag = "B6-A-378-8494";
      if ( bias == "250" ) tag = "B6-A-381-8497";
      if ( bias == "200" ) tag = "B6-A-382-8498";
      if ( bias == "150" ) tag = "B6-A-383-8499";
      if ( bias == "100" ) tag = "B6-A-384-8500";
      if ( bias == "75" )  tag = "B6-A-385-8501";
      if ( bias == "50" )  tag = "B6-A-386-8502";
    } else if(sector == "6"){
      if ( bias == "300" ) tag = "B6-A-408-8523";
      if ( bias == "250" ) tag = "B6-A-409-8524";
      if ( bias == "200" ) tag = "B6-A-410-8525";
      if ( bias == "150" ) tag = "B6-A-411-8526";
      if ( bias == "100" ) tag = "B6-A-412-8527";
      if ( bias == "75" )  tag = "B6-A-413-8528";
      if ( bias == "50" )  tag = "B6-A-414-8529";
    }
  } else if (board.Contains("A4")) {
    if(sector == "1") {
      if(bias == "400" ) tag = "B4-A-210-8552";
      if(bias == "350" ) tag = "B4-A-211-8553";
      if(bias == "300" ) tag = "B4-A-212-8554";
      if(bias == "250" ) tag = "B4-A-213-8555";
      if(bias == "200" ) tag = "B4-A-214-8556";
      if(bias == "150" ) tag = "B4-A-215-8557";
      if(bias == "100" ) tag = "B4-A-216-8558";
      if(bias == "75" )  tag = "B4-A-217-8559";
      if(bias == "50" )  tag = "B4-A-218-8560";
    }else if(sector == "2") {
      if(bias == "400" ) tag = "B4-A-242-8583";
      if(bias == "350" ) tag = "B4-A-243-8584";
      if(bias == "300" ) tag = "B4-A-246-8586";
      if(bias == "250" ) tag = "B4-A-247-8587";
      if(bias == "200" ) tag = "B4-A-248-8588";
      if(bias == "150" ) tag = "B4-A-249-8589";
      if(bias == "100" ) tag = "B4-A-250-8590";
      if(bias == "75" )  tag = "B4-A-251-8591";
      if(bias == "50" )  tag = "B4-A-252-8592";
    } else if(sector == "3"){
      if(bias == "400" ) tag = "B4-A-275-8615";
      if(bias == "350" ) tag = "B4-A-276-8616";
      if(bias == "300" ) tag = "B4-A-277-8617";
      if(bias == "250" ) tag = "B4-A-278-8618";
      if(bias == "200" ) tag = "B4-A-279-8619";
      if(bias == "150" ) tag = "B4-A-280-8620";
      if(bias == "100" ) tag = "B4-A-281-8621";
      if(bias == "75" )  tag = "B4-A-282-8622";
      if(bias == "50" )  tag = "B4-A-283-8623";
    } else if(sector == "4"){
      if(bias == "400" ) tag = "B4-A-217-8651";
      if(bias == "350" ) tag = "B4-A-218-8652";
      if(bias == "300" ) tag = "B4-A-219-8653";
      if(bias == "250" ) tag = "B4-A-220-8654";
      if(bias == "200" ) tag = "B4-A-221-8655";
      if(bias == "150" ) tag = "B4-A-222-8656";
      if(bias == "100" ) tag = "B4-A-223-8657";
      if(bias == "75" )  tag = "B4-A-224-8658";
      if(bias == "50" )  tag = "B4-A-225-8659";
    } else if(sector == "5"){
      if(bias == "400" ) tag = "B4-A-246-8680";
      if(bias == "350" ) tag = "B4-A-247-8681";
      if(bias == "300" ) tag = "B4-A-248-8682";
      if(bias == "250" ) tag = "B4-A-249-8683";
      if(bias == "200" ) tag = "B4-A-250-8684";
      if(bias == "150" ) tag = "B4-A-251-8685";
      if(bias == "100" ) tag = "B4-A-252-8686";
      if(bias == "75" )  tag = "B4-A-253-8687";
      if(bias == "50" )  tag = "B4-A-254-8688";
    } else if(sector == "6"){
      if(bias == "400" ) tag = "B4-A-227-8711";
      if(bias == "350" ) tag = "B4-A-228-8712";
      if(bias == "300" ) tag = "B4-A-229-8713";
      if(bias == "250" ) tag = "B4-A-230-8714";
      if(bias == "200" ) tag = "B4-A-231-8715";
      if(bias == "150" ) tag = "B4-A-232-8716";
      if(bias == "100" ) tag = "B4-A-233-8717";
      if(bias == "75" )  tag = "B4-A-234-8718";
      if(bias == "50" )  tag = "B4-A-235-8719";
      tag = tag.ReplaceAll("B4","B1");
    }
  } else if (board.Contains("A8")) {
    if(sector == "1"){
      if(bias == "500") tag = "B8-A-296-13332";
      if(bias == "400") tag = "B8-A-297-13333";
      if(bias == "350") tag = "B8-A-298-13334";
      if(bias == "300") tag = "B8-A-299-13335";
      if(bias == "250") tag = "B8-A-300-13336";
      if(bias == "200") tag = "B8-A-301-13337";
      if(bias == "150") tag = "B8-A-302-13338";
      if(bias == "100") tag = "B8-A-303-13339";
      if(bias == "75")  tag = "B8-A-304-13340";
      if(bias == "50")  tag = "B8-A-305-13341";
      if(m_scanType=="Angle" && angle=="10" && bias == "300") tag = "B8-A-378-13399";
      if(m_scanType=="Angle" && angle=="20" && bias == "300") tag = "B8-A-379-13400";
    }else if (sector == "2"){
      //if(bias == "500") tag = "B8-A-321-13354";
      if(bias == "400") tag = "B8-A-324-13355";
      if(bias == "350") tag = "B8-A-325-13356";
      if(bias == "300") tag = "B8-A-327-13357";
      if(bias == "250") tag = "B8-A-328-13358";
      if(bias == "200") tag = "B8-A-329-13359";
      if(bias == "150") tag = "B8-A-330-13360";
      if(bias == "100") tag = "B8-A-331-13361";
      if(bias == "75")  tag = "B8-A-332-13362";
      if(bias == "50") tag = "B8-A-333-13363";
    }else if (sector == "3"){
      //if(bias == "500") tag = "B8-A-358-13385";
      if(bias == "400") tag = "B8-A-359-13386";
      if(bias == "350") tag = "B8-A-360-13387";
      if(bias == "300") tag = "B8-A-361-13388";
      if(bias == "250") tag = "B8-A-365-13389";
      if(bias == "200") tag = "B8-A-366-13390";
      if(bias == "150") tag = "B8-A-367-13391";
      if(bias == "100") tag = "B8-A-368-13392";
      if(bias == "75")  tag = "B8-A-369-13393";
      if(bias == "50")  tag = "B8-A-370-13394";
    }else if (sector == "4"){
      //if(bias == "500") tag = "B8-A-309-13344";
      if(bias == "400") tag = "B8-A-311-13345";
      if(bias == "350") tag = "B8-A-312-13346";
      if(bias == "300") tag = "B8-A-314-13347";
      if(bias == "250") tag = "B8-A-315-13348";
      if(bias == "200") tag = "B8-A-316-13349";
      if(bias == "150") tag = "B8-A-317-13350";
      if(bias == "100") tag = "B8-A-318-13351";
      if(bias == "75")  tag = "B8-A-319-13352";
      if(bias == "50")  tag = "B8-A-320-13353";
    }else if (sector == "5"){
      //if(bias == "500") tag = "B8-A-334-13364";
      if(bias == "400") tag = "B8-A-337-13365";
      if(bias == "350") tag = "B8-A-339-13367";
      if(bias == "300") tag = "B8-A-340-13368";
      if(bias == "250") tag = "B8-A-341-13369";
      if(bias == "200") tag = "B8-A-342-13370";
      if(bias == "150") tag = "B8-A-343-13371";
      if(bias == "100") tag = "B8-A-344-13372";
      if(bias == "75")  tag = "B8-A-345-13373";
      if(bias == "50")  tag = "B8-A-346-13374";
      if(bias != "400") tag = tag.ReplaceAll("B8","B1");
    }else if (sector == "6"){
      //if(bias == "500") tag = "B8-A-347-13375";
      if(bias == "400") tag = "B8-A-348-13376";
      if(bias == "350") tag = "B8-A-349-13377";
      if(bias == "300") tag = "B8-A-350-13378";
      if(bias == "250") tag = "B8-A-352-13379";
      if(bias == "200") tag = "B8-A-353-13380";
      if(bias == "150") tag = "B8-A-354-13381";
      if(bias == "100") tag = "B8-A-355-13382";
      if(bias == "75")  tag = "B8-A-356-13383";
      if(bias == "50")  tag = "B8-A-357-13384";
      tag = tag.ReplaceAll("B8","B1");
    }
  } else if (board.Contains("D5")) {
    if(sector == "1"){
      if(bias == "50") tag = "B5-D-14-13066";
      if(bias == "75") tag = "B5-D-13-13065";
      if(bias == "100") tag = "B5-D-12-13064";
      if(bias == "150") tag = "B5-D-11-13063";
      if(bias == "200") tag = "B5-D-10-13062";
      if(bias == "250") tag = "B5-D-6-13058";
      if(bias == "300") tag = "B5-D-19-13071";
      if(bias == "350") tag = "B5-D-4-13056";
      if(bias == "400") tag = "B5-D-3-13055";
      if(bias == "500") tag = "B5-D-16-13068";
    }else if(sector == "2"){
      if(bias == "50")  tag = "B5-D-119-13165";
      if(bias == "75")  tag = "B5-D-118-13164";
      if(bias == "100") tag = "B5-D-117-13163";
      if(bias == "150") tag = "B5-D-116-13162";
      if(bias == "200") tag = "B5-D-115-13161";
      if(bias == "250") tag = "B5-D-114-13160";
      if(bias == "300") tag = "B5-D-113-13159";
      if(bias == "350") tag = "B5-D-111-13157";
      if(bias == "400") tag = "B5-D-110-13156";
      if(bias == "500") tag = "B5-D-109-13155";
    }else if(sector == "3"){
      if(bias == "50")  tag = "B5-D-167-13209";
      if(bias == "75")  tag = "B5-D-166-13208";
      if(bias == "100") tag = "B5-D-165-13207";
      if(bias == "150") tag = "B5-D-164-13206";
      if(bias == "200") tag = "B5-D-163-13205";
      if(bias == "250") tag = "B5-D-162-13204";
      if(bias == "300") tag = "B5-D-159-13202";
      if(bias == "350") tag = "B5-D-158-13201";
      if(bias == "400") tag = "B5-D-157-13200";
      if(bias == "500") tag = "B5-D-156-13199";
    }else if(sector == "4"){
      if(bias == "50")  tag = "B5-D-65-13114";
      if(bias == "75")  tag = "B5-D-64-13113";
      if(bias == "100") tag = "B5-D-63-13112";
      if(bias == "150") tag = "B5-D-62-13111";
      if(bias == "200") tag = "B5-D-61-13110";
      if(bias == "250") tag = "B5-D-59-13108";
      if(bias == "300") tag = "B5-D-58-13107";
      if(bias == "350") tag = "B5-D-57-13106";
      if(bias == "400") tag = "B5-D-56-13105";
      if(bias == "500") tag = "B5-D-55-13104";
    }else if(sector == "5"){
      if(bias == "50")  tag = "B5-D-218-13253";
      if(bias == "75")  tag = "B5-D-217-13252";
      if(bias == "100") tag = "B5-D-216-13251";
      if(bias == "150") tag = "B5-D-215-13250";
      if(bias == "200") tag = "B5-D-214-13249";
      if(bias == "250") tag = "B5-D-213-13248";
      if(bias == "300") tag = "B5-D-212-13247";
      if(bias == "350") tag = "B5-D-211-13246";
      if(bias == "400") tag = "B5-D-210-13245";
      if(bias == "500") tag = "B5-D-209-13244";
    }else if(sector == "6"){
      if(bias == "50")  tag = "B5-D-260-13294";
      if(bias == "75")  tag = "B5-D-259-13293";
      if(bias == "100") tag = "B5-D-258-13292";
      if(bias == "150") tag = "B5-D-257-13291";
      if(bias == "200") tag = "B5-D-256-13290";
      if(bias == "250") tag = "B5-D-255-13289";
      if(bias == "300") tag = "B5-D-254-13288";
      if(bias == "350") tag = "B5-D-253-13287";
      if(bias == "400") tag = "B5-D-252-13286";
      if(bias == "500") tag = "B5-D-251-13285";
    }
  } else if (board.Contains("D7")) {
    if(sector == "1"){
      if(bias == "50") tag = "B1-D-347-9360";
      if(bias == "75") tag = "B1-D-346-9359";
      if(bias == "100") tag = "B1-D-344-9356";
      if(bias == "150") tag = "B1-D-345-9357";
      if(bias == "200") tag = "B1-D-343-9355";
      if(bias == "300") tag = "B1-D-343-9354";
      if(bias == "400") tag = "B1-D-341-9353";
      if(bias == "500") tag = "B1-D-340-9272";
    }else if(sector == "2"){
      if(bias == "50") tag = "B1-D-471-9659";
      if(bias == "75") tag = "B1-D-470-9658";
      if(bias == "100") tag = "B1-D-469-9657";
      if(bias == "150") tag = "B1-D-468-9656";
      if(bias == "200") tag = "B1-D-467-9655";
      if(bias == "300") tag = "B1-D-466-9654";
      if(bias == "400") tag = "B1-D-465-9653";
      if(bias == "500") tag = "B1-D-464-9652";
    }else if(sector == "3"){
      if(bias == "50") tag = "B1-D-472-9660";
      if(bias == "75") tag = "B1-D-473-9661";
      if(bias == "100") tag = "B1-D-474-9662";
      if(bias == "150") tag = "B1-D-475-9663";
      if(bias == "200") tag = "B1-D-476-9664";
      if(bias == "300") tag = "B1-D-477-9665";
      if(bias == "400") tag = "B1-D-478-9666";
      if(bias == "500") tag = "B1-D-479-9667";
    }else if(sector == "4"){
      if(bias == "50") tag =  "B1-D-379-9424";
      if(bias == "75") tag =  "B1-D-378-9423";
      if(bias == "100") tag = "B1-D-377-9422";
      if(bias == "150") tag = "B1-D-376-9421";
      if(bias == "200") tag = "B1-D-375-9419";
      if(bias == "300") tag = "B1-D-374-9417";
      if(bias == "400") tag = "B1-D-373-9415";
      if(bias == "500") tag = "B1-D-372-9413";
    }else if(sector == "5"){
      if(bias == "50")  tag = "B1-D-489-9763";
      if(bias == "75")  tag = "B1-D-488-9762";
      if(bias == "100") tag = "B1-D-487-9761";
      if(bias == "150") tag = "B1-D-484-9674";
      if(bias == "200") tag = "B1-D-483-9673";
      if(bias == "300") tag = "B1-D-482-9672";
      if(bias == "400") tag = "B1-D-481-9671";
      if(bias == "500") tag = "B1-D-480-9670";
    }else if(sector == "6"){
      if(bias == "50") tag =  "B1-D-321-9252";
      if(bias == "75") tag =  "B1-D-320-9251";
      if(bias == "100") tag = "B1-D-319-9250";
      if(bias == "150") tag = "B1-D-318-9249";
      if(bias == "200") tag = "B1-D-317-9248";
      if(bias == "300") tag = "B1-D-316-9247";
      if(bias == "400") tag = "B1-D-315-9246";
      if(bias == "500") tag = "B1-D-314-9245";
    }
  } else if (board.Contains("A1")) {
    if(sector == "1"){    
      if(bias == "350") tag = "B1-A-1022-13947";
      if(bias == "300") tag = "B1-A-1021-13946";
      if(bias == "250") tag = "B1-A-1053-13971";
      if(bias == "200") tag = "B1-A-1052-13970";
      if(bias == "150") tag = "B1-A-1015-13939";
      if(bias == "100") tag = "B1-A-1014-13938";
      if(bias == "75")  tag = "B1-A-1013-13937";
      if(bias == "50")  tag = "B1-A-1012-13936";
    }else if(sector == "2"){
      if(bias == "350") tag = "B1-A-1062-13980";
      if(bias == "300") tag = "B1-A-1061-13979";
      if(bias == "250") tag = "B1-A-1060-13978";
      if(bias == "200") tag = "B1-A-1059-13977";
      if(bias == "150") tag = "B1-A-1058-13976";
      if(bias == "100") tag = "B1-A-1057-13975";
      if(bias == "75")  tag = "B1-A-1056-13974";
      if(bias == "50")  tag = "B1-A-1055-13973";
    }    
  } else if (board.Contains("A2")) {
    if(sector == "1"){    
      if(bias == "340") tag = "B2-A-1095-14030";
      if(bias == "325") tag = "B2-A-1094-14029";
      if(bias == "300") tag = "B2-A-1084-14017";
      if(bias == "250") tag = "B2-A-1083-14016";
      if(bias == "200") tag = "B2-A-1082-14015";
      if(bias == "150") tag = "B2-A-1081-14014";
      if(bias == "100") tag = "B2-A-1078-14013";
      if(bias == "75")  tag = "B2-A-1077-14012";
      if(bias == "50")  tag = "B2-A-1076-14011";
      if(bias == "300" && angle=="-1") tag = "B2-A-1085-14018";
      if(bias == "300" && angle=="1") tag = "B2-A-1088-14021";
    }else if(sector == "2"){
      if(bias == "340") tag = "B2-A-1099-14034";
      if(bias == "325") tag = "B2-A-1098-14033";
      if(bias == "300") tag = "B2-A-1097-14032";
      if(bias == "250") tag = "B2-A-1100-14035";
      if(bias == "200") tag = "B2-A-1101-14036";
      if(bias == "150") tag = "B2-A-1102-14037";
      if(bias == "100") tag = "B2-A-1103-14038";
      if(bias == "75")  tag = "B2-A-1104-14039";
      if(bias == "50")  tag = "B2-A-1105-14040";
    }    
  }

  return tag;
}


#endif // #ifdef ClusterWithTrackAna_cxx